Сходимость в моделях с shell181

[dimkasat 2]

Имею следующую проблему.

Модель состоит из листового металла, поэтому смоделирована с спользованием surface bodies (shell181 elements type) . В результате нагружения в стыках наблюдаются сингулярные элементы: в нодах, принадлежащих двум листам одновременно, как правило в углах, напряжения нереальные и не сходятся при улучшении сетки. Погуглив вопрос, узнал, что это неизбежный недостаток shell элементов. Включение пластики далеко не всегда помогает. Выход: исключать элементы, содержащие эти ноды из постпроцечссинга.

 

Вопрос: есть ли у кого пример АPDL сниппета, исключающего элементы, отвечающие какому либо геометрическими признаку (принадлежность двум плоскостями например)?

 

Буду благодарен! В APDL новичок я.

[Борман 2 375]

Все что можно формализовать, можно запрограммировать на APDL.

[hr4d 87]

Для начала просто почитайте в хелпе про селект. И самое просто это тупо анселектнуть ненужные элементы.

upd. Гребаный онлайн хелп, так бы точно сказал. Что-то типа selection logic в mechanical apdl.

А вообще конечно это старое, но актуальности не потеряло. Я по ней учился апдльке

spravochnik-programmista-ansys-apdl-439421702.pdf

Изменено 14 декабря 2018 пользователем hr4d

[dimkasat 2]

В 14.12.2018 в 17:48, hr4d сказал:

Для начала просто почитайте в хелпе про селект. И самое просто это тупо анселектнуть ненужные элементы.

upd. Гребаный онлайн хелп, так бы точно сказал. Что-то типа selection logic в mechanical apdl.

А вообще конечно это старое, но актуальности не потеряло. Я по ней учился апдльке

spravochnik-programmista-ansys-apdl-439421702.pdf

Видимо надо отсеять элементы, прилегающие к edges (lines в apdl если я правильно понимаю).

[soklakov 1 475]

В 14.12.2018 в 18:48, hr4d сказал:

Гребаный онлайн хелп

можно ж поставить локальный

[hr4d 87]

7 часов назад, soklakov сказал:

можно ж поставить локальный

Нет админ прав(( Во вторник поставлю, на другом компе был 15-й, а тут обновили до 19.1

8 часов назад, dimkasat сказал:

Видимо надо отсеять элементы, прилегающие к edges (lines в apdl если я правильно понимаю).

Просто выкиньте элементы. Можно скриптик простой написать, но тут больше Борман или Федо поможет.

Изменено 16 декабря 2018 пользователем hr4d

[Борман 2 375]

Конкретно то что надо?

[kolo666 253]

3 часа назад, hr4d сказал:

Просто выкиньте элементы.

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

Я отступаю всегда на толщину+катет сварки, при этом этот отступ состоит минимум из двух элементов. Есть еще идеи?

Изменено 16 декабря 2018 пользователем kolo666

[Борман 2 375]

5 минут назад, kolo666 сказал:

Есть еще идеи?

Отказываться от критериев прочности в напряжениях, и переходить к силам-моментам, как это делается в СП Стальные конструкции (номер не помню, 16.*** кажись).

[Fedor 1 597]

Это возврат от Коши к Галилею  :)

[dimkasat 2]

8 часов назад, hr4d сказал:

Нет админ прав(( Во вторник поставлю, на другом компе был 15-й, а тут обновили до 19.1

Просто выкиньте элементы. Можно скриптик простой написать, но тут больше Борман или Федо поможет.

Не вижу пока как определить принадлежность ноды  к элементам разных тел (part/group). Разве что пересмотреть все ноды по каждому телу и отсеять ноды с повторяющимися id? У них же есть id! ?

Изменено 16 декабря 2018 пользователем dimkasat

[dimkasat 2]

20 часов назад, kolo666 сказал:

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

Я отступаю всегда на толщину+катет сварки, при этом этот отступ состоит минимум из двух элементов. Есть еще идеи?

Ну это и есть т.н. hot spot stress анализ: выбирают точки на расстоянии 0.4 и 1 t (t-толщина листа) и по подтвержденным экспериментально завиимостям аппроксимируют значение напряжения на краю, где его невозможно определить в принципе не зная точной геометрии сварного шва и его мех свойств. Однако этот способ невозможно автоматизировать: каждый раз надо отыскивать опасное место, определять направление нормали к шву, вдоль нормали отбирать указанные точки, определять значения в них. Если сборка большая и концентраторов несколько, или, например, речь идет об определении несущей способности, то это все становится слишком трудоемко

21 час назад, Борман сказал:

Конкретно то что надо?

Модель состоит из shell181. Необходимо в workbench написать скрипт в command объекте, который при постпроцессинге исключал из рассмотрения элементы, содержащие ноды, принадлежащие двум телам одновременно (ибо в них из-за усреднения значения напряжений неверные).

[soklakov 1 475]

4 часа назад, dimkasat сказал:

Необходимо в workbench написать скрипт в command объекте, который при постпроцессинге исключал из рассмотрения элементы, содержащие ноды, принадлежащие двум телам одновременно (ибо в них из-за усреднения значения напряжений неверные).

мне кажется, это можно сделать с помощью Named Selection и worksheet

ну просто в качестве альтернативы) да и удобнее

ну может не удобнее))

[dimkasat 2]

В 17.12.2018 в 13:49, soklakov сказал:

мне кажется, это можно сделать с помощью Named Selection и worksheet

ну просто в качестве альтернативы) да и удобнее

ну может не удобнее))

Да, но не могу подобрать критерий принадлежности элемента определенному телу сборки

[dimkasat 2]

В 17.12.2018 в 13:49, soklakov сказал:

мне кажется, это можно сделать с помощью Named Selection и worksheet

ну просто в качестве альтернативы) да и удобнее

ну может не удобнее))

Да, но не могу подобрать критерий принадлежности элемента определенному телу сборки

А конкретно, не могу понять как получить I'd элементов, находящихся на одной line.

[soklakov 1 475]

1 час назад, dimkasat сказал:

А конкретно, не могу понять как получить I'd элементов, находящихся на одной line.

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

[soklakov 1 475]

В 14.12.2018 в 17:16, dimkasat сказал:

исключающего элементы, отвечающие какому либо геометрическими признаку (принадлежность двум плоскостями например определенным линиям из выборки)?

 

[Борман 2 375]

3 часа назад, dimkasat сказал:

Да, но не могу подобрать критерий принадлежности элемента определенному телу сборки

vol_ID=123 ! номер определенного тела сборки

vol_Name='SOMEBODY' ! или имя определенного тела сборки

elem_ID=321 ! номер элемента

 

allsel,all

vsel,s,,,vol_ID (или) cmsel,s,vol_Name

 

eslv,s

 

*if,esel(elem_ID),eq,1,then

ЭЛЕМЕНТ ВЫБРАН

*endif

 

*if,esel(elem_ID),eq,-1,then

ЭЛЕМЕНТ НЕ ВЫБРАН

*endif

 

3 часа назад, dimkasat сказал:

не могу понять как получить I'd элементов, находящихся на одной line.

Это пока тяжело понять.

[hr4d 87]

В 16.12.2018 в 13:16, kolo666 сказал:

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

Я отступаю всегда на толщину+катет сварки, при этом этот отступ состоит минимум из двух элементов. Есть еще идеи?

Я обычно смотрю на nodal difference и принимаю решения исходя из этой инфы, не знаю как в фемапе эта функция называется.

[kolo666 253]

10 часов назад, hr4d сказал:

Я обычно смотрю на nodal difference и принимаю решения исходя из этой инфы, не знаю как в фемапе эта функция называется.

У нас это Error Estimate, но какое отношение имеет ошибка осреднения к сингулярности?

Изменено 20 декабря 2018 пользователем kolo666