Ошибки при решении

[Wise 0]

У вас просто предупреждения о том, что форма некоторых элементов не очень нравится.

Я на такую ерунду даже внимания не обращаю. Сложная задача не червонец, чтобы все нравиться.

Ну будет возможно число обусловленности чуть побольше, мелочь это.

Мне тоже кажется, что ничего криминального не произошло. Так а почему все таки тогда возникает эта ошибка.

[Борман 2 375]

Wise

Скиньте лист всех применяемых меатериалов. Прямо из ансиса.

[Fedor 1 597]

А вот это и надо найти. Тупо попробуйте переразбить покрупнее и долбитесь частями.

Отладка она и есть отладка. Без нее не бывает

Меньше думайте, больше проверяйте

Ошибка обычно там, где ее меньше всего ожидаешь, в очевидном и под рукой

Проверьте хоть десяток элементов то решается, если выкинуть остальные ?

[Fedor 1 597]

"Основное время уходит на перемножение разреженной матрицы на вектор" - ее никогда не будет,

коль нет указателей. Такие дела, обычно ручками приходится.

Наверное что-то вроде dll или соm на выходе. Не переписываете же всю программу

Главное скалярные произведения. Наверное надо выбирать согласованные ненулевые и отправлять, выборка много дешевле умножений, должно быть ускорение.

Ну а раскидать по процессорам наверное можно через блочные разбиения как в литературе описано. по разреженным матрицам.

Это при разложении на множители, ну а там где есть итерации все проще, молоти да молоти

"Никаких ближних и дальних указателей нет" - да эта архаика при нынешних разрядностях операционок не нужна. На 16 было актуально Просто вспомнилось.

[Fedor 1 597]

Да мне не до этого, просто из любопытства

[sergeyd 3]

перед solve

наберите

/outp,test,outp

появится файл test.outp

(либо в меню найдите switch output to file...)

советую так делать всегда. иначе не разобраться.

(откуда, скажем вы судите, что ДИАгностики не было?)

его и нужно вывесить, если нужна помощь.

в принципе можно отключить проверку пивотов (есть недокумументированная опция),

но лучше разобраться. иначе будет забивание шурупов молотком.

версии еще есть, но гадать не хочу. вывешивайте файл.

если москвич, пишите в личку. привозите всю задачу.

[Wise 0]

Wise

Скиньте лист всех применяемых меатериалов. Прямо из ансиса.

EVALUATE MATERIAL PROPERTIES FOR MATERIALS 1 TO 1 IN INCREMENTS OF 1

MATERIAL NUMBER = 1 EVALUATED AT TEMPERATURE OF 0.0000

EX = 0.18500E+12

NUXY = 0.30000

PRXY = 0.30000

[Wise 0]

вывешивайте файл.

AT ALL PICKED OR COMPONENT AREAS,

SET DESIRED ELEMENT EDGE LENGTH TO 0.25000E-02

GENERATE NODES AND ELEMENTS IN ALL PICKED VOLUMES

Initiating meshing of volume 1 exterior.

Meshing of volume 1 exterior is complete. Facet count = 9036.

Initiating meshing of volume 1 interior.

Estimated number of elements to fill this volume = 66300.

Estimated number of nodes = 82900.

Begin meshing of volume 1

Completed pre-check of triangle facets

*** WARNING *** CP = 929.844 TIME= 09:43:50

Mid-nodes of some elements have been modified to lie on straight edges

because of distortion with the original mid-nodes. This condition is

sometimes eliminated by tetrahedron element improvement when enabled.

Issue ESEL,,STRAIGHTENED to select such elements for listing or

plotting.

Meshing of volume 1 is complete.

Now storing 66012 nodes and 56291 elements.

Initiating meshing of volume 2 exterior.

Meshing of volume 2 exterior is complete. Facet count = 3170.

Initiating meshing of volume 2 interior.

Estimated number of elements to fill this volume = 13700.

Estimated number of nodes = 17100.

Begin meshing of volume 2

Completed pre-check of triangle facets

Meshing of volume 2 is complete.

Now storing 13731 nodes and 12631 elements.

Initiating meshing of volume 3 exterior.

Meshing of volume 3 exterior is complete. Facet count = 4662.

Initiating meshing of volume 3 interior.

Estimated number of elements to fill this volume = 20100.

Estimated number of nodes = 25100.

Begin meshing of volume 3

Completed pre-check of triangle facets

*** WARNING *** CP = 931.578 TIME= 09:43:51

Mid-nodes of some elements have been modified to lie on straight edges

because of distortion with the original mid-nodes. This condition is

sometimes eliminated by tetrahedron element improvement when enabled.

Issue ESEL,,STRAIGHTENED to select such elements for listing or

plotting.

Meshing of volume 3 is complete.

Now storing 21826 nodes and 19799 elements.

Initiating meshing of volume 5 exterior.

Meshing of volume 5 exterior is complete. Facet count = 4072.

Initiating meshing of volume 5 interior.

Estimated number of elements to fill this volume = 32700.

Estimated number of nodes = 40800.

Begin meshing of volume 5

Completed pre-check of triangle facets

Meshing of volume 5 is complete.

Now storing 37731 nodes and 31326 elements.

NUMBER OF VOLUMES MESHED = 4

MAXIMUM NODE NUMBER = 178967

MAXIMUM ELEMENT NUMBER = 120047

------------------------------------------------------------------------------

<<<<<< SHAPE TESTING SUMMARY >>>>>>

<<<<<< FOR NEW OR MODIFIED ELEMENTS >>>>>>

------------------------------------------------------------------------------

--------------------------------------

| Element count 120047 SOLID92 |

--------------------------------------

Test Number tested Warning count Error count Warn+Err %

---- ------------- ------------- ----------- ----------

Aspect Ratio 120047 0 0 0.00 %

Maximum Angle 120047 0 0 0.00 %

Jacobian Ratio 120047 1 0 0.00 %

Any 120047 1 0 0.00 %

------------------------------------------------------------------------------

*** WARNING *** CP = 932.547 TIME= 09:43:52

Shape testing revealed that 1 of the 120047 new or modified elements

violate shape warning limits. To review test results, please see the

output file or issue the CHECK command.

PRODUCE ELEMENT PLOT IN DSYS = 0

вывешивайте файл.

***** ANSYS SOLVE COMMAND *****

TRANSFER SOLID MODEL BOUNDARY CONDITIONS TO FINITE ELEMENT MODEL

*** NOTE *** CP = 472.078 TIME= 09:30:30

Nodes on symmetry/antisymmetry surfaces will be rotated as required.

CONSTRAINTS TRANSFERRED FROM AREAS = 1134

SYMM/ASYM CONSTRAINTS FROM AREAS = 4242

SURFACE LOADS TRANSFERRED FROM AREAS = 168

*** WARNING *** CP = 472.672 TIME= 09:30:31

Previous testing revealed that 1 of the 120047 selected elements

violate shape warning limits. To review warning messages, please see

the output or error file, or issue the CHECK command.

*** NOTE *** CP = 472.672 TIME= 09:30:31

The model data was checked and warning messages were found.

Please review output or errors file ( D:\Ansys\Stoika mel\Stoika

sigma R.err ) for these warning messages.

S O L U T I O N O P T I O N S

PROBLEM DIMENSIONALITY. . . . . . . . . . . . .3-D

DEGREES OF FREEDOM. . . . . . UX UY UZ

ANALYSIS TYPE . . . . . . . . . . . . . . . . .STATIC (STEADY-STATE)

GLOBALLY ASSEMBLED MATRIX . . . . . . . . . . .SYMMETRIC

*** NOTE *** CP = 472.766 TIME= 09:30:31

Present time 0 is less than or equal to the previous time.

Time will default to 1.

*** NOTE *** CP = 472.875 TIME= 09:30:31

The conditions for direct assembly have been met. No .emat or .erot

files will be produced.

L O A D S T E P O P T I O N S

LOAD STEP NUMBER. . . . . . . . . . . . . . . . 1

TIME AT END OF THE LOAD STEP. . . . . . . . . . 1.0000

NUMBER OF SUBSTEPS. . . . . . . . . . . . . . . 1

STEP CHANGE BOUNDARY CONDITIONS . . . . . . . . NO

PRINT OUTPUT CONTROLS . . . . . . . . . . . . .NO PRINTOUT

DATABASE OUTPUT CONTROLS. . . . . . . . . . . .ALL DATA WRITTEN

FOR THE LAST SUBSTEP

SOLUTION MONITORING INFO IS WRITTEN TO FILE= Stoika sigma R.mntr

Element Formation Element= 1000 Cum. Iter.= 1 CP= 475.406

Time= 1.0000 Load Step= 1 Substep= 1 Equilibrium Iteration= 1.

Range of element maximum matrix coefficients in global coordinates

Maximum= 3.109870721E+09 at element 22626.

Minimum= 106356870 at element 104024.

*** ELEMENT MATRIX FORMULATION TIMES

TYPE NUMBER ENAME TOTAL CP AVE CP

1 120047 SOLID92 96.734 0.000806

Time at end of element matrix formulation CP= 487.484375.

SPARSE MATRIX DIRECT SOLVER.

Number of equations = 531605, Maximum wavefront = 309

Memory available for solver = 356.212 MB

Memory required for in-core = 4206.110 MB

Optimal memory required for out-of-core = 222.728 MB

Minimum memory required for out-of-core = 57.037 MB

*** ERROR *** CP = 918.344 TIME= 09:32:08

Large negative pivot value ( -232646489 ) in Eqn.system. May be

because of a bad temperature-dependent material property used in the

model.

Sparse solver maximum pivot= 2.989864602E+09 at node 14201 UZ.

Sparse solver minimum pivot= -232646489 at node 113220 UY.

Sparse solver minimum pivot in absolute value= 23184775 at node 110489

UZ.

*** WARNING *** CP = 918.375 TIME= 09:32:10

Small negative equation solver pivot term encountered at UZ DOF of node

110489. Check for an insufficiently constrained model.

Изменено 30 июня 2009 пользователем Wise

[Fedor 1 597]

" 14201 113220 110489 " - вот и посмотрите, что-то около этих узлов барахлит.

Еще попробуйте 186 элемент, он получше.

Обусловленность матрицы хорошая, грубо оценить число потерянных знаков в точности

можно как десятичный логарифм от отношения максимального к минимальному элементу на диагонали

Коль просто трехмерка, без оболочек и балок, то просто переразбейте грубой сеткой и проверьте

[Wise 0]

" 14201 113220 110489 " - вот и посмотрите, что-то около этих узлов барахлит.

Элементы содержащие эти узлы unsellect ом убираю. Ставлю на счет, опять эта же ошибка. Выводит еще не очень хорошие узлы. Убираю их то же и снова на счет, опять эта ошибка. Такое ощущение, что это будет продолжаться пока у меня вообще элементов не останется.

Разбивал грубо, ошибка та же самая.

[Fedor 1 597]

Тогда очевидно, что где-то на поверхностях разреза числятся одинаковые по виду но по сути не совпадающие поверхности.

Включите их нумерацию и покрутите, в таких областях обычно цвета рябят.

Посмотрите аналогично линии, где-то неустойчивый цвет меняющийся будет, обычно видно.

Так бывает при неаккуратном моделировании. В итоге узлы элементов не все стыкуются и будут или дыры или вот как у Вас

Обычно булевскими операциями удается поправить.

[Wise 0]

Тогда очевидно, что где-то на поверхностях разреза числятся одинаковые по виду но по сути не совпадающие поверхности.

Включите их нумерацию и покрутите, в таких областях обычно цвета рябят.

Посмотрите аналогично линии, где-то неустойчивый цвет меняющийся будет, обычно видно.

Так бывает при неаккуратном моделировании. В итоге узлы элементов не все стыкуются и будут или дыры или вот как у Вас

Обычно булевскими операциями удается поправить.

Перед разбиение использовал Merg а затем GLU.Разве может остся не совпадающие поверхности?

Узлы на которые ругается находятся не на одной поверхности,а разбросаны по всей модели.

Изменено 30 июня 2009 пользователем Wise

[Fedor 1 597]

Легко. Например одна поверхность может иметь n линий, а другая m они будуть внешне совпадать, но будут по сути разными.

"Оверлапить" надо, или иногда склеивать при соприкасающихся линиях поверхностей.

В общем с высокой вероятностью проблема в геометрической модели.

Ищите частями с грубым разбиением для скорости

"Такое ощущение, что это будет продолжаться пока у меня вообще элементов не останется" - останутся те, которые в связном графе с краевыми условиями находятся

[Wise 0]

Что-то странное, построил обыкновенный кубик, задал перемещение четырем точкам лежащие в одной плоскости перемещение по одной оси. Появляется такая же ошибка. Когда креплю по одной поверхности по всем направлениям и нагружаю давлением то все нормально считает. Это разве нормально?

[Fedor 1 597]

"задал перемещение четырем точкам лежащие в одной плоскости перемещение по одной оси" а по другим пусть летят куда хотят?

[Wise 0]

"задал перемещение четырем точкам лежащие в одной плоскости перемещение по одной оси" а по другим пусть летят куда хотят?

Закрепил в остальных направлениях. Та же самая ошибка.

[Fedor 1 597]

Точки то на одной прямой лежат?

И уменьшите давление на несколько порядков

[Wise 0]

Точки то на одной прямой лежат?

И уменьшите давление на несколько порядков

Подождите.

Когда задается перемещение четырем точкам,то давления не прикладываю.Я хочу что бы у меня кубик переместился.

[Wise 0]

Точки то на одной прямой лежат?

И уменьшите давление на несколько порядков

Закрепление 4-х точек по всем направлениям и приложения давления дает такой же не хороший результат. И уменьшение давления то же плохой результат.

[Fedor 1 597]

Закрепление 4-х точек по всем направлениям и приложения давления дает такой же не хороший результат. И уменьшение давления то же плохой результат.

Задайте трем точкам не лежащих на одной прямой какие-нибудь одинаковые вектора перемещений, сколько-нибудь осмысленные, например 1,2,3