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

[Wise 0]

Здравствуйте! Помогите пожалуйста. При моделировании з-х мерной задачи и с использованием элемента solid 187.В процессе решения системы уравнений всплывает следующие сообщение:

large negative pivot value in Eqn.system.May be because of a bad temperature-dependent material property used in the model.

Далее решение останавливается.Подскажите пожалуйста что делать в таком случае.

[Форумный боец 1]

а в модели есть

temperature-dependent material property

?

но вообще такая ошибка указывает на расходимость решения. это может быть вызвано некорректным закреплением, некорректными параметрами модели, некорректными параметрами решателя, некорректным выбором решателя и т.д.

я не уверен, что написал не чушь. и я не уверен, что в данном случае что-то такое :)

[sergeyd 3]

задача линейная?

просто проверь закрепления.

или (вероятнее) не соединены части модели.

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

нелинейная?

мб что угодно. нужны подробности.

например просто разрушилось в пластике все вокруг какого-то узла

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

[Wise 0]

Задача линейная. Задача была уже решена, только на 32-х битной системе, хотелось сделать сетку по корректнее и решить на 64-х битной, может в этом дело? Хотя на сколько я знаю это не должно повлиять. После перестроения сетки при решении возникает такая ошибка. После чего попробовал старую базу данных решить заново на 64-х битной системе, опять тоже самое. При создании модели использовал операцию склеивание, так что разлететься не может. Крепление были старые, но может быть я что-то забыл. Попробую изменить крепления.

Спасибо большое.

[Crippling 0]

Крепление были старые, но может быть я что-то забыл. Попробую изменить крепления.

В смысле, перемешивал модель и не перезадал закрепления и силы?

[Wise 0]

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

[Crippling 0]

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

Ну так сетку то меняли..., номера узлов и элементов у вас изменились. Был к примеру, закреплен узел № 1547489, после перемешивания, положения, номера и количество узлов изменились. Узел этот 1547489, скорее всего остался, и закреплен, как и было, но к модели он уже не имеет ни какого отношения, т.е. он просто висит в воздухе, а модель с уже другими узлами и элементами, не закреплена. Проверьте закрепления.

[Wise 0]

Крепления были переназначены, на новой модели.

К тому же закрепления были приложены к твердотельной модели(ТТМ). И при запуске на счет они переносятся на конечно-элементную модель (КЭМ), если ТТМ не изменилась, то и граничные условия остались без изменения и перенесутся соответственно на новую КЭМ. Даже если я пытаюсь запустить на решение старую базу данных,( без изменения сетки) то результат тот же самый.

Импортировал геометрию, создавал сетку, прикладывал нагрузки, еще раз, что бы проверить корректность поставленной задачи, результат тот же самый.

Очень странная ситуация, даже не могу предположить с чем это может быть связано. Попробовал перерешать еще одну модель,которая так же была посчитана ранее, результат аналогичный. Изменяю граничные условия такой ошибки нет, но они далеки от реальных.

Самое странное, что это уже было решено, но на другой машине и на другой платформе (32 бит), и проблем подобного рода не возникали.

[Crippling 0]

А с материалом там все нормально?

У меня выдавалась точно такая ошибка, но в двух случаях. Либо была модель не закреплена (или узлы не слиты). Либо когда материал задавался с очень маленьким модулем упругости (получались большие перемещения)

[Crippling 0]

Задача линейная. Задача была уже решена, только на 32-х битной системе, хотелось сделать сетку по корректнее и решить на 64-х битной, может в этом дело? Хотя на сколько я знаю это не должно повлиять. После перестроения сетки при решении возникает такая ошибка. После чего попробовал старую базу данных решить заново на 64-х битной системе, опять тоже самое.

Т.е. если я правильно понял, даже если просто берешь модель, которая на 32 битной машине решалась, - на 64 битной она не решается? Попробуй выбрать другой решатель, поиграть точностью.

Еще сейчас вспомнил, такая фигня у меня была, когда геометрия модели очень говенная. Если просто смотришь на модель, то кажеться все нормально, но в одном месте, ни с того ни с сего сетка загущена неправдоподобно, а присмотришся, увеличишь в 50 раз, там или микроотверстие или очень острый край поверхности порядка 1-5 градусов. И как раз такая ошибка у меня тогда вылетала, на таких моделях, причем через раз!!! Нагрузку уменьшишь чуть-чуть (или точность изменишь), - все нормально, увеличишь - опять эта ошибка вылетает. Причем в тех случаях когда решалось, на диаграмме деформаций в этих местах были нереально большие значения. Лечилось исправлением геометрии.

P/s. У меня 64-битная система тоже.

[Wise 0]

Материал стандартная сталь Е=2е11, коэффициент Пуассона 0,29.

Модель состоит из одной симметричной детали, которая закреплена по поверхности и по всем направлениям, а так же условия симметрии. Вот к стати, когда убираю условия симметрии и плоскость симметрии так же закрепляю по всем направлениям, то все считает.

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

[sergeyd 3]

так помочь невозможно...

наверняка были диагностики при формировании матриц.

что в *.Err?

направь вывод в текст и вывесь файл.

непонятно: когда всю задачу считаешь, как крепишь?

думаю, напутано с ГУ.

частоты посчитай!!!! формы посмотри!

+

метод sparse или pcg?

ну или вывесь всю задачу, только архивируй (cdwr).

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

[Crippling 0]

частоты посчитай!!!! формы посмотри!

А что это даст? Просто мне тоже интересно.

[Борман 2 385]

А что это даст? Просто мне тоже интересно.

Если детать свободна, то будут формы с нулевой частотой без всяких деформаций.

[Fedor 1 605]

Самое первое стоит отключить все условия и нагрузки и задать простой поворот. То есть шарнирно закрепить пару точек и третьей не лежащей на прямой задать перемещение побольше в ортогональном направлении. Тело должно повернуться и увеличить свои размеры, но напряжения и деформации должны быть нулевыми. Если этого не произошло, значит напахано в элементах, иначе в условиях и нагрузках

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

[Wise 0]

так помочь невозможно...

наверняка были диагностики при формировании матриц.

что в *.Err?

направь вывод в текст и вывесь файл.

непонятно: когда всю задачу считаешь, как крепишь?

думаю, напутано с ГУ.

частоты посчитай!!!! формы посмотри!

+

метод sparse или pcg?

ну или вывесь всю задачу, только архивируй (cdwr).

1. Стояла опция Program chosen solver.Изменил на PCG ,такой ошибки не возникло. Поменял на Sparse возникает эта ошибка.

2. Частоты посчитал, нулевых нет.

3.При создании сетки не ругался, хотя сетка грубовата.

4. Креплю по поверхности и по всем направлениям.

[Wise 0]

То есть шарнирно закрепить пару точек

Я же использую solid92,а у него есть только линейные перемещения.Т.е. в любом случае закрепленная точка не будет ни перемещаться,ни поворачиваться.

Как тело может увеличить свои размеры,но при этом деформации останутся нулевыми?

[sergeyd 3]

я думал, ошибка при PCG..

первый раз наблюдаю обратное. Чудеса!

кстати, время счета сравнил? (Pcg при средней точности (1e-6) раз в 10 быстрее на больших задачах!)

сколько итераций?

повторю:выходной файл вывешивать-то будешь?

[Fedor 1 605]

Это не зависит от элементов. Кососимметричные преобразования, а поворот как раз к таковым

относится, дают нулевые деформации при ненулевых перемещениях. Это свойство теории упругости.

Легко проверить аналитически.

Если конечный элемент таким свойством не обладает, значит в нем нарушается "полнота" базизных функций и стоит его остерегаться

Как простейший вариант можете задать постоянные одинаковые перемещения в трех не лежащих на прямой точках.

Тоже позволяет легко проверить корректность матрицы

Кососимметричные со сдвигом тоже

"Изменил на PCG ,такой ошибки не возникло. Поменял на Sparse возникает эта ошибка" - при разложении есть деление

на диагональный или извлечение квадратного корня иногда. Скорее всего что-то не связано и тогда это обычно возникает.

При сопряженных градиентах ошибки округления превращают его в итерационный метод да еще с предобуславливанием.

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

[sergeyd 3]

Уважаемый Федор!

я вроде бы должен Вас понимать, а не пойму...

слова и термины по-отдельности все-все понимаю, а в целом- ну никак!

молодежь совсем запутается...

<

сопряженные градиенты не нужно превращать, они и так ВСЕГДА итерационный метод...

и предобуславливание (оно тоже бывает разное) обычно проводится

с использованием прямого разложения матрицы, близкой (в той или иной степени) исходной матрице

и делается ДО сопряженных градиентов.

ну и т.д.

элементы (ЧТО 92 ЧТО 187) все проверены на бенчмарках.

причем здесь кососимметрия...

автор!

еще версия- пром.узлы где-то не посередине ребер.

кстати если есть 12 версия-

удалите промежуточные узлы emid,remove и примените 285 тип тетраэдры 4-узловые

интересно, что будет!

диагностики не было или не заметили?

теперь основная версия- стыкуются элем с оч. большими разницами жесткостей (из-за размеров эл)

где-то есть крохотульки...

если доли мм, то в СИ коэфф будут типа 1E20, а это не есть хорошо.

давайте файл ! там пишут о подобном!

еще версия- разница в 32 и 64 показывает на обнаруженную нехорошесть разбивки

+

температура нигде не задавалась?

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