contact, remesh, implicit

[aleksei_babichev 0]

Коллеги, помогите настроить перестройку в контактной задаче в неявной постановке. В другом пакете считал это на раз два MSC.Marc, в Ls-Dayna новичок, а хотелось бы считать задачи с большим шагом в статике с большими нелинейностями.

[aleksei_babichev 0]

2 часа назад, aleksei_babichev сказал:

Коллеги, помогите настроить перестройку в контактной задаче в неявной постановке. В другом пакете считал это на раз два MSC.Marc, в Ls-Dayna новичок, а хотелось бы считать задачи с большим шагом в статике с большими нелинейностями.

Расчет заканчивается ошибкой Error 41436 (SOL+1436) Too many random files opened

[AlexKaz 1 823]

@aleksei_babichev , если прикрепите пример проекта чтобы наглядно посмотреть проблему, ответ найдётся быстрее.

[Victoria 104]

В статике?

[aleksei_babichev 0]

29.12.2023 в 22:10, AlexKaz сказал:

@aleksei_babichev , если прикрепите пример проекта чтобы наглядно посмотреть проблему, ответ найдётся быстрее.

https://disk.yandex.ru/d/LmpjYw6GvRKEww

 

Спасибо за отклик и всех с праздником прошедшим!

30.12.2023 в 03:30, Victoria сказал:

В статике?

Да, к сожалению, в статике, времена расчетов просто 1е14 сек.

[AlexKaz 1 823]

@aleksei_babichev ,

 

12_10_2023_TOPA_3_mod.7z

*control_timestep -> dt2ms ставьте больше 1, например, 1e5 и выше.

*damping -> global можно значительно увеличить.

[AlexKaz 1 823]

37 минут назад, AlexKaz сказал:

*control_timestep -> dt2ms ставьте больше 1, например, 1e5 и выше.

Ммм. Не раскурил. До time=1e13-1e16 дойти не удаётся...

IMHO,  расчёт фактически запирается на 1e7-1e8 по времени т.к. в *control_adaptive время обновления сетки задано константой 1e7. Чтобы скорректировать время перестроения, надо бы явно задать времена перестроения по нарастающей через таблицу, там же параметр LCADP.

 

Как мне кажется, лучше всё же поменять систему единиц, перейдя от секунд к дням-годам. 

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

Чтобы явно задать статику/динамику в неявной схеме на всякий случай всё же нужна карта *CONTROL_DYNAMICS.

[aleksei_babichev 0]

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

@aleksei_babichev ,

 

12_10_2023_TOPA_3_mod.7z 55 \u041a\u0431 · 1 скачивание

*control_timestep -> dt2ms ставьте больше 1, например, 1e5 и выше.

*damping -> global можно значительно увеличить.

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

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

Ммм. Не раскурил. До time=1e13-1e16 дойти не удаётся...

IMHO,  расчёт фактически запирается на 1e7-1e8 по времени т.к. в *control_adaptive время обновления сетки задано константой 1e7. Чтобы скорректировать время перестроения, надо бы явно задать времена перестроения по нарастающей через таблицу, там же параметр LCADP.

 

Как мне кажется, лучше всё же поменять систему единиц, перейдя от секунд к дням-годам. 

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

Чтобы явно задать статику/динамику в неявной схеме на всякий случай всё же нужна карта *CONTROL_DYNAMICS.

 

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

Ммм. Не раскурил. До time=1e13-1e16 дойти не удаётся...

У меня этот файл доходит до 1.5609e9, число шагов 2156 шагов. Дело в решателе скорее всего, НО !!!! замечено мной и не только один и тот же к-файл один решатель на разных компах и даже дисках решает по разному, или я что-то не так делал.

 

Насчет перестройки: В других нелинейных решателях есть куча настроек, где число, период перестройки меняется по разным признакам-искажение угла, накопленная деформация и т.д. В данном случае в LS Dyna этого нет. По этой причине тут особо нет смысла по кривой заставлять его перестраивать сетку. Пробовал по разному, с этими настройками он дальше всего проходит, ну или я делаю что-то не то.

 

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

 

Может быть стоит уйти от постановки с конечными элементами уже и перейти на SPH? Или что-то другое более "пробивное", ясно что в явной постановке все будет решаться, но не с моими временами. 

[AlexKaz 1 823]

2 часа назад, aleksei_babichev сказал:

Переход к годам или часам .... можно, но эта модель только старт

Материалы в СИ, а размерность геометрии какая? Они согласованы?

[aleksei_babichev 0]

16 часов назад, AlexKaz сказал:

Материалы в СИ, а размерность геометрии какая? Они согласованы?

Да, конечно!

[AlexKaz 1 823]

Добирался до 1.69е9.

10.01.2024 в 14:12, aleksei_babichev сказал:

Может быть стоит уйти от постановки с конечными элементами уже и перейти на SPH?

Можно посмотреть в сторону EFG. Потом, если не подойдёт, SPH.

[aleksei_babichev 0]

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

Добирался до 1.69е9.

Можно посмотреть в сторону EFG. Потом, если не подойдёт, SPH.

Спасибо! Вроде же на EFG нет перестройки или я не прав? Смотрел демонстрацию где штамп, большие деформации, но у меня еще больше могут быть и перестройка лучше бы была для моих задач. SPH конечно круто, но незнакома она мне, как и многим... Ну что делать, придется возиться с ней, надеюсь на помощь сообщества. 

[AlexKaz 1 823]

2 часа назад, aleksei_babichev сказал:

на EFG нет перестройки или я не прав

На сайте dynaexamples.com как раз есть пример с EFG c авторемэшингом.

Изменено 12 января пользователем AlexKaz

[aleksei_babichev 0]

15 часов назад, AlexKaz сказал:

На сайте dynaexamples.com как раз есть пример с EFG c авторемэшингом.

Вспомнил, EFG с remesh идет только в 3D постановке, начинал именно с него, задача большой становится , времени много на отладку, сперва бы 2D освоить.

[aleksei_babichev 0]

12.01.2024 в 20:03, AlexKaz сказал:

На сайте dynaexamples.com как раз есть пример с EFG c авторемэшингом.

А не могли бы Вы вкратце рассказать про контакт, философию в LS DYNA? Никак не могу понять, там в контактных картах надо указывать контактные пары (когда у контактного тела может быть несколько пар или минимум одна) или все со всеми или множества с множествами, как лучше?

[AlexKaz 1 823]

2 часа назад, aleksei_babichev сказал:

про контакт, философию в LS DYNA?

Нет никакой философии.

В 3д наиболее корректный тип для нетребовательных пользователей automatic node to surface.

Для склейки - tied node to surface offset https://cccp3d.ru/topic/112014-tied-contact-для-склеивания-2d-и-3d/

 

 

2 часа назад, aleksei_babichev сказал:

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

Для automatic single surface  можно указать только один раз поверхности или их набор.

Для всех остальных обычно достаточно выбрать подтип (SSTYP, MSTYP) part-part (№3) и указать соответствующие part'ы. Другие подтипы (наборы узлов, наборы part'ов, наборы граней и т.п.) обычно работают. Если какой-то подтип не работает - в логе возникает ошибка, решатель останавливается и выдаёт в логе сообщение что и как нужно исправить.

Все остальные настройки встречаются редко на практике, и для того, чтобы их понять, достаточно справки или мануалов JDay https://ftp.lstc.com/anonymous/outgoing/jday/

Если мануалы Jday не открываются - постучите мне в личку.

[aleksei_babichev 0]

Спасибо за совет! По ссылке файлы не качаются пишет-"Not a file". И в личку не удается постучать, может быть по тому что я новичок, пишет "Вам разрешено отправлять 0 сообщений в день...". Может в телеге @Babichev_Alex

[aleksei_babichev 0]

23 часа назад, AlexKaz сказал:

Нет никакой философии.

В 3д наиболее корректный тип для нетребовательных пользователей automatic node to surface.

Для склейки - tied node to surface offset https://cccp3d.ru/topic/112014-tied-contact-для-склеивания-2d-и-3d/

 

 

Для automatic single surface  можно указать только один раз поверхности или их набор.

Для всех остальных обычно достаточно выбрать подтип (SSTYP, MSTYP) part-part (№3) и указать соответствующие part'ы. Другие подтипы (наборы узлов, наборы part'ов, наборы граней и т.п.) обычно работают. Если какой-то подтип не работает - в логе возникает ошибка, решатель останавливается и выдаёт в логе сообщение что и как нужно исправить.

Все остальные настройки встречаются редко на практике, и для того, чтобы их понять, достаточно справки или мануалов JDay https://ftp.lstc.com/anonymous/outgoing/jday/

Если мануалы Jday не открываются - постучите мне в личку.

Я имел ввиду, если к примеру в задаче 4 тела и контакт к примеру осуществляется между 1-2, 1-3, 2-3, 3-4 надо заводить контакты в карте CONTACT_AUTOMATIC_NODES_TO_SURFACE для каждой контактной пары тел или можно просто указать два множества А и В? 

[AlexKaz 1 823]

Можно. SSTYP/MSTYP eq 2.

[aleksei_babichev 0]

16.01.2024 в 13:28, AlexKaz сказал:

достаточно справки или мануалов JDay https://ftp.lstc.com/anonymous/outgoing/jday/

Если мануалы Jday не открываются - постучите мне в личку.

Не получается в личку стучать. Как быть?