Ls-Dyna встроенная в WB и Explict

[rasta89 113]

Всем привет.

В WB 2021 есть встроенная Ls-Dyna, это очень круто, что можно пользоваться этим решателем, используя привычный интерфейс механикла, хоть и решатель отстает по скорости от трушной Ls-Dyna.

Решил сделать тест, проверить сверление в Explict и в Dyna (встроенная), проверить скорость расчета, по идее скорость расчетов следующая:

1. Быстрее всех решит высокочастотная трушная Ls-Dyna

2. На втором месте встроенная Ls-Dyna

3. Долгий Explict.

 

Если у кого-то есть туториалы по встроенной Ls-Dyna, поделитесь пожалуйста?

 

Возникла масса вопросов.

1. Ls-Dyna

Задал параметры по-минимуму стандартными инструментами:

- сверло сталь, плашка алюминий

- у сверла и плашки своя коор система (для перемещения и вращения)

- body interaction, задавал как frictional 0.4 (и dynamics 0.2) так и frictionless

- сетка 1 мм сверло, и 0,7 плашка

- вращаю сверло 250 об/мин

- перемещаю плашку на 30 мм вдоль сверла

не считает, спустя 40-60 мин расчета, все стоит на месте, нужно использовать rigid body?

делал два расчета, как описывал эту задачу выше и с некоторыми настройками в anslysis settings и использованием rigid body, выложил оба проекта

export_ls1.wbpzexport_ls.wbpz

 

2. Explict

Задача такая же, отличие в том, что добавил контакт:

export_exp.wbpz

 

Но возникает ошибка: Time step too small, почитал форумы, пишут это ошибка минимального тайм степа, протестил 1000 вариаций, не помогло((

 

Если по счастливой случайности будет время у кого-нибудь посмотреть, буду очень рад. Спасибо.

[Victoria 104]

СИ - так более привычно, хотя очень много примеров для Дайны из=за бугра, где тонны и т.д.

[rasta89 113]

6 минут назад, Victoria сказал:

СИ - так более привычно, хотя очень много примеров для Дайны из=за бугра, где тонны и т.д.

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

 

По теме

[Victoria 104]

А масса и не задается - есть плотность, а программа на основе конечно-разностной модели посчитает объем

[rasta89 113]

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

А масса и не задается - есть плотность, а программа на основе конечно-разностной модели посчитает объем

В системе СИ

Время расчета 2 сек

Перемещение

Вращение

Материал

 

 

[Victoria 104]

?

[rasta89 113]

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

?

Перевел в СИ и теперь сверло просто вращается на месте

 

попробую с перемещением поиграться конечно

сейчас перемещение такое, в течении 2 сек на 6 мм

время расчета такое же

[rasta89 113]

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

попробую с перемещением поиграться конечно

установил значение 6, вместо 0,006, ради интереса

[Victoria 104]

Нужен файл - что-то тут не так...

[rasta89 113]

25 минут назад, Victoria сказал:

Нужен файл - что-то тут не так...

приложил

wb_prep_model_solid_solve.key

[Bonus 70]

Длину у Вас задает сетка, расстояние между узлами. Переход из одной системы в другую требует масштабирование модели. Вы это делали?

 

Цитата

Размерность каждой единицы согласовывается по правилу: сила=масса*ускорение, ускорение=длина/время2.

 

Изменено 20 июля 2022 пользователем Bonus

[rasta89 113]

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

Длину у Вас задает сетка, расстояние между узлами. Переход из одной системы в другую требует масштабирование модели. Вы это делали?

 

 

не делал, но проверил размеры, все соотв

[Bonus 70]

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

не делал, но проверил размеры, все соотв

выше

Цитата

Перевел в СИ и теперь сверло просто вращается на месте

 

попробую с перемещением поиграться конечно

сейчас перемещение такое, в течении 2 сек на 6 мм

 

Если у вас модель построена в мм, то это не система СИ, так как СИ это м. Так что если у вас ранее был диаметр 10мм, как расстояние между крайними узлами. То переведя в СИ, т.е. в м, нужно масштабировать модель, что бы расстояние между теми же узлами было 0,01мм. Отсюда если у вас 6мм, то либо вес в тоннах (плотность 7,8Е-9), либо время в мс. Но это в любом случае не СИ. СИ это перемещение 0,006мм.

[rasta89 113]

16 минут назад, Bonus сказал:

выше

 

Если у вас модель построена в мм, то это не система СИ, так как СИ это м. Так что если у вас ранее был диаметр 10мм, как расстояние между крайними узлами. То переведя в СИ, т.е. в м, нужно масштабировать модель, что бы расстояние между теми же узлами было 0,01мм. Отсюда если у вас 6мм, то либо вес в тоннах (плотность 7,8Е-9), либо время в мс. Но это в любом случае не СИ. СИ это перемещение 0,006мм.

Да, не в си, а как вы говорите, был тест выше по форуму, тоже не работает(

[Bonus 70]

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

Да, не в си, а как вы говорите, был тест выше по форуму, тоже не работает(

 

Цитата

*MAT_RIGID_TITLE
drill
$#     mid        ro         e        pr         n    couple         m     alias
         3    7850.02.00000E11       0.3       0.0       0.0       0.0       

Плотность у вас в кг/м3 модуль в Па. Это из вашего файла. Значит расстояние должно быть в метрах. Это СИ. Растояние между узлами считается в метрах. Смотрим узлы:

Цитата

*NODE
$#   nid               x               y               z      tc      rc  
       1        -2.62263        1.456643           -3.12       0       0
       2        2.749506         1.20009           -3.12       0       0
       3       -2.410133        1.786409           -3.12       0       0
       4       -2.156423        2.085627           -3.12       0       0
       5       -1.865837        2.349181           -3.12       0       0
       6       -1.543346        2.572564           -3.12       0       0
       7       -1.194463        2.751955           -3.12       0       0
       8      -0.8251551        2.884288           -3.12       0       0
       9      -0.4417368          2.9673           -3.12       0       0
      10      -0.0507647         2.99957           -3.12       0       0
      11       0.3410754        2.980548           -3.12       0       0
      12       0.7270831        2.910558           -3.12       0       0
      13        1.100658        2.790798           -3.12       0       0
      14        1.455411        2.623315           -3.12       0       0
      15        1.785277        2.410972           -3.12       0       0
      16        2.084614        2.157403           -3.12       0       0
      17        2.348304        1.866941           -3.12       0       0
      18        2.571838        1.544554           -3.12       0       0
      19        2.438794       0.4525096           -3.12       0       0
      20        2.584131       0.8305165           -3.12       0       0
      21        2.467266         -1.0688           -3.12       0       0
      22        2.357823      -0.6959516           -3.12       0       0
      23        2.316013      -0.3096288           -3.12       0       0
      24        2.343255       0.0779938           -3.12       0       0
      25         2.62263       -1.456643           -3.12       0       0
      26       -2.749506        -1.20009           -3.12       0       0
      27        2.410133       -1.786409           -3.12       0       0
      28        2.156423       -2.085627           -3.12       0       0
      29        1.865837       -2.349181           -3.12       0       0
      30        1.543346       -2.572564           -3.12       0       0
      31        1.194463       -2.751955           -3.12       0       0
      32       0.8251551       -2.884288           -3.12       0       0
      33       0.4417368         -2.9673           -3.12       0       0

У вас габариты модели по 3 метра?

[rasta89 113]

6 минут назад, Bonus сказал:

 

Плотность у вас в кг/м3 модуль в Па. Это из вашего файла. Значит расстояние должно быть в метрах. Это СИ. Растояние между узлами считается в метрах. Смотрим узлы:

У вас габариты модели по 3 метра?

а, когда в си переводил значения не подумал об этом((

значит первый вариант был верный, не могли бы вы его поглядеть?

приложил его еще раз, как раз он у меня на ЯД есть

https://disk.yandex.ru/d/_84B-w4rKNKQdw

[Bonus 70]

10 минут назад, rasta89 сказал:

а, когда в си переводил значения не подумал об этом((

значит первый вариант был верный, не могли бы вы его поглядеть?

приложил его еще раз, как раз он у меня на ЯД есть

https://disk.yandex.ru/d/_84B-w4rKNKQdw

Здесь у Вас 4 схема, по той таблице что Вы приводили выше. Сетка походу в мм.

 

Проверьте время в мсек, усилия в кН, я это не проверял.

 

Цитата

сейчас перемещение такое, в течении 2 сек на 6 мм

Т.е. возможно это 2 мсек. А не сек. Но проверяйте сами.

[rasta89 113]

21 минуту назад, Bonus сказал:

Здесь у Вас 4 схема, по той таблице что Вы приводили выше. Сетка походу в мм.

 

Проверьте время в мсек, усилия в кН, я это не проверял.

 

Т.е. возможно это 2 мсек. А не сек. Но проверяйте сами.

Для варианта в системе си, это секунды, для варианта не в системе си это мс. Проверял сто раз(

Там в варианте не си, все работает но сверло отскакивает от плашки. Я конечно ещё контакт и add_erosion не гонял на разных вариантах, но может быть дело и не в них. В общем не получается довести до состояния сверления.

[Bonus 70]

23 минуты назад, rasta89 сказал:

Для варианта в системе си, это секунды, для варианта не в системе си это мс. Проверял сто раз(

Там в варианте не си, все работает но сверло отскакивает от плашки. Я конечно ещё контакт и add_erosion не гонял на разных вариантах, но может быть дело и не в них. В общем не получается довести до состояния сверления.

Вообще, такую задачку может есть смысл считать в Fluid-Structure Interaction тогда обрабатываемый материал будет сетка Эйлера (сходная с CFD) там надо сетку сделать на материал и немного воздуха вокруг, а сверло Лагранжа (привычная нам структурная сетка). Вообще, рассмотрите вариант когда сверло абсолютно твердое. Тогда его вообще можно задать только поверхностью. Fluid-Structure Interaction сложнее в моделировании. Но тогда отпадет вопрос с деформированием/искажением до недопустимости сетки, да и с контактами будет как то по другому.

 

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

 

З.Ы. в Дайне сто лет не работал, зашел сюда просто так.

[rasta89 113]

13 минут назад, Bonus сказал:

Вообще, такую задачку может есть смысл считать в Fluid-Structure Interaction тогда обрабатываемый материал будет сетка Эйлера (сходная с CFD) там надо сетку сделать на материал и немного воздуха вокруг, а сверло Лагранжа (привычная нам структурная сетка). Вообще, рассмотрите вариант когда сверло абсолютно твердое. Тогда его вообще можно задать только поверхностью. Fluid-Structure Interaction сложнее в моделировании. Но тогда отпадет вопрос с деформированием/искажением до недопустимости сетки, да и с контактами будет как то по другому.

 

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

 

З.Ы. в Дайне сто лет не работал, зашел сюда просто так.

сверло итак уже Rigid тело)

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

[Victoria 104]

Я повторюсь об использовании автором темы  в своей работе метод ТЫКА. Это для начала

Размеры сверла и плашки заданы в метрах. В препосте  Дайнв можно просто масштабировать с к-том 1000 - дело 20 секунд.

Дальше лучше. В материале плашки не задан критерий разрушения - в модели материала к-ты Di - (их 5 штук) везде стоят 0. Хотя достаточно и одного, но нужно выбрать какой будет критерий разрушения.

Можно было-бы оставить и все 0, но в другой альтернативной карте  *MAT_ADD_EROSION тоже все нули.

Из руководства (том 2) "Многие определяющие модели программы LS-DYNA не содержат разрушения и эрозии материала. Опция ADD_EROSION предоставляет способ включения разрушения в эти модели, а также может использоваться с определяющими моделями, которые содержат другие критерии разрушения/эрозии."

 

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

лежащие на поверхности. Можно задать как взаимодействие двух частей

$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 
*CONTACT_ERODING_SURFACE_TO_SURFACE_ID
$#     cid                                                                 title
        12
$#    ssid      msid     sstyp     mstyp    sboxid    mboxid       spr       mpr
         1         2         3         3         0         0         0         0
$#      fs        fd        dc        vc       vdc    penchk        bt        dt
       0.4       0.1       0.0       0.0      10.0         0       0.0   1.00E20
$#     sfs       sfm       sst       mst      sfst      sfmt       fsf       vsf
       1.0       1.0       0.0       0.0       1.0       1.0       1.0       1.0
$#    isym    erosop      iadj    
         0         1         1
$#    soft    sofscl    lcidab    maxpar     sbopt     depth     bsort    frcfrq
         2       0.1         0     1.025       2.0         2         0         1
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8

Желательно ввести карту *CONTROL_CONTACT и еще кучу мелких исправлений

Критерии не выбирал - модель материала - с потолка