[Cfytrr]

Siemens PLM Femap with NX Nastran "Real FEA made easy"




Тема посвящается вопросам, возникающим при работе в среде программного комплекса Femap

Примеры выполнения некоторых расчетов можно посмотреть:
на домашней страничке

на сайте www.iberisa.com http://iberisa.wordpress.com,
Форум (англ.)
Линейный анализ; Нелинейный анализ
Форум на домашней страничке

[Студентик]

Сообщение от Cfytrr Студентик, 1. У Вас очень мелкая сетка,думаю, что 650000 элементов многовато для этой задачи. Для пластины достаточно 5 элементов по толщине, а в плане элементы могут быть раз в 10 побольше. 2. В gap элементах задан начальный зазор, пока он не выберется, пластина висит в воздухе, следовательно имеем механизм и соответствующую ошибку. 3. В gap элементах не указана сдвиговая жесткость и коэфф трения покоя, и нет закреплений у пластины которые помешают ей уехать в бок, следовательно опять механизм. 4. Жесткость Gap элементов не должна сильно отличаться от жесткости того что они соединяют. 5. Gap элементы штука капризная (имхо) и чем их больше, тем хуже, для Вашей задачи надо использовать контакт.

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

[Cfytrr]

Это значит назначить два контактных региона, и коннектор между ними. В этой теме уже не раз обсуждалось. На домашней страничке есть видео с примером.

[Студентик]

Сообщение от Cfytrr Это значит назначить два контактных региона, и коннектор между ними. В этой теме уже не раз обсуждалось. На домашней страничке есть видео с примером.

Спасибо большое за совет,буду пробовать

[konotop]

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

[Cfytrr]

konotop,
1. У Вас в схеме есть несшитые узлы, они не сшиваются так как имеют разные Output Csys

это происходит, например, когда копируем (зеркалим) узлы вместе с нагрузками и закреплениями.Лечится так:
Modify - Update Other - Output Csys и сшить все повторно.
2. Опять таки Gap имеют начальный зазор, это приводит к тому, что плита находится в фазе полета. Gap при линейном расчете работают как простые связи заданной жесткости. Жесткость лучше сделать повыше.
3. При наложении закреплений на солиды, указывать закрепление от поворота напрасное дело, так как они его игнорируют.

Я не понял работу схемы в целом, Вы что, просто плющите торец плиты?
Файл чуть подправленный:

[konotop]

Cfytrr, спасибо
на этом этапе - да. Просто плющу :-)
далее в пустоты буду вставлять заглушки, а также арматурные каркасы, но пока что надо было разобраться с ошибкой...

От чего зависит размер файла? у меня сейчас файл, который весил 1.5 Мб весит 180 Мб. Удалил результаты расчета, лишние группы, лишние элементы и узлы... но существенно это на ситуацию не повлияло... Как можно привести его к человеческому размеру ?

[golnr]

Здравствуйте) уже долгое время пытаюсь решить простую нестационарную задачу на плоскости XY (аналитическое решение: T=x+t ):
d^2(T)/(dx)^2 + d^2(T)/(dy)^2 + Q = dT/dt. где Q=1, Kxx=1, Kyy=1, dT/dt=1

Для простоты задал на границе условие 1го рода.


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

link: http://dl.dropbox.com/u/45772592/A1.rar

[selega]

Сообщение от golnr уже долгое время пытаюсь решить простую нестационарную задачу на плоскости XY (аналитическое решение: T=x+t ): d^2(T)/(dx)^2 + d^2(T)/(dy)^2 + Q = dT/dt. где Q=1, Kxx=1, Kyy=1, dT/dt=1 Для простоты задал на границе условие 1го рода. В результатах температура то увеличивается то уменьшается с каждым шагом времени, скажите кто знает почему?

golnr, вероятно, я что-то не понял. Поясните, пожалуйста - в чем заключается Ваша нелинейность? И зачем Вы задавали нагрузки в узлах? Судя по дифференциальному уравнению у Вас имеется функциональная зависимость температуры от времени. Значит, Вам нужно было описать эту функцию уравнением или набором значений в характерных точках. У Вас же функция задана единичным значением, при этом Вы задали набор отдельных нагружений. И решатель выбран стационарный.
Изобразите от руки суть Вашей задачи, тогда можно будет сказать что-то более внятное.

[golnr]

http://dl.dropbox.com/u/45772592/07052012536.JPG

При решении задачи я взял отрезок времени (0; 8)

менеджере анализа задал: Analysis Set -> Analisys Type = 21..Transient Heat Transfer

Nonlinear control option-> Number of steps = 8 Initial Time Increment = 1; Output Step Interval = 0;

[selega]

Сообщение от golnr http://forum.dwg.ru/attachment.php?a...1&d=1336398402

По Вашей ссылке

[golnr]

я исправил ссылку

[selega]

Сообщение от golnr я исправил ссылку

Спасибо. А можно обычными словами, чтобы не морочиться над нюансами теоретического тепломассобмена? В смысле, что изначально задано и что происходит

[golnr]

Сообщение от selega golnr, вероятно, я что-то не понял. Поясните, пожалуйста - в чем заключается Ваша нелинейность? И зачем Вы задавали нагрузки в узлах? Судя по дифференциальному уравнению у Вас имеется функциональная зависимость температуры от времени. Значит, Вам нужно было описать эту функцию уравнением или набором значений в характерных точках. У Вас же функция задана единичным значением, при этом Вы задали набор отдельных нагружений. И решатель выбран стационарный. Изобразите от руки суть Вашей задачи, тогда можно будет сказать что-то более внятное.

так мне надо эту задачу через "20..Steady-State Heat Transfer" надо было делать ?


Вообще мне надо решить задачу так: из известного аналитического уравнения получаем производные и значения уравнения на границе (граничные условия), и источник тепла..
..и если правильно задать начальные и граничные условия, то график результата должен соответствовать аналитическому решению

[selega]

Сообщение от golnr так мне надо эту задачу через "20..Steady-State Heat Transfer" надо было делать ?

Прошу прощения, по поводу решателю все правильно. Но вот в остальном...

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

В учебнике Рудакова этот вопрос рассмотрен, стр. 122-131.

[golnr]

нету там ничего полезного ( в Рудакове )

[Voyager2007]

konotop, после стирания результатов нужно перестроить базу данных командой File - Rebuild. После чего запишите файл, при этом его размер сильно уменьшится. Жесткость GAP элементов должна быть на 2-3 порядка выше жесткости соединяемых деталей в направлении контакта.

[konotop]

Voyager2007, спасибо

[selega]

Сообщение от golnr нету там ничего полезного ( в Рудакове )

Шимкович, изд. 2008 г. На стр. 650 есть пример расчета нестационарного температурного поля.

[golnr]

Здравствуйте) хочу сравнить решение на своей программе с расчётами на Nastran

задача такая: квадратная анизотропная (Kxx=1, Kyy=20) пластина с координатами x=(0,3) и y=(0,3); на границе известна температура T(верхн)=x*x+73; T(справа)=109-y*y*y; T(снизу)=x*x+100; T(слева)=100-y*y*y; Ещё есть источник: Q=120y-2


создал сетку 23х23 точки(разбил на треугольные элементы 22x22)

задаю Гранич условия: Model->Load->on curve->выбрал верхнюю границу->выбрал Temperature->задал Value=1, выбрал Method=variable->нажал advanced->в окне ввел !x*!x+73;

задаю Гранич условия: Model->Load->on curve->выбрал правую границу->выбрал Temperature->задал Value=1, выбрал Method=variable->нажал advanced->в окне ввел 109-!y*!y*!y;

задаю Гранич условия: Model->Load->on curve->выбрал нижнюю границу->выбрал Temperature->задал Value=1, выбрал Method=variable->нажал advanced->в окне ввел !x*!x+100;

задаю Гранич условия: Model->Load->on curve->выбрал левую границу->выбрал Temperature->задал Value=1, выбрал Method=variable->нажал advanced->в окне ввел 100-!y*!y*!y;



задаю источник: Model->Load->on surface->выбрал пластину->выбрал Heat generation->задал Value=-1/529 т. к. точек 529, выбрал Method=variable->нажал advanced->в окне ввел 120*!y-2;

посчиталось.., только график не тот что должен быть, правильно ли я задал условия??


может источник надо домножать, потому что он получается маленьким.. но тогда я не знаю на какой коэффициент..??

[selega]

Сообщение от golnr Здравствуйте) хочу сравнить решение на своей программе с расчётами на Nastran

Очень благородная цель! Только файл забыли приложить.