Error 10125 key 125 ls dyna - Исправление ошибок и поиск оптимальных решений проблем

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

Положительное по знаку давление на поверхности элемента действует нормально поверхности «внутрь» элемента.

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

З а м е ч а н и е! Последовательность действий при выборе элементов 3D Tet-Solid 168 и определение опций полностью соответствует элементу

SOLID164.

Рисунок 76. Элемент SOLID168

Элементы Thin Shell 163

Элемент SHELL163 представляет собой плоский четырехугольник и используется для моделирования тонкостенных оболочечных конструкций. Данный элемент имеет четыре узла, но может вырождаться и принимать треугольную форму.

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

Рисунок 77. Формы элемента Shell163

Последовательность действий при определении опций элементов

Thin Shell 163

Выбираем пункт меню Main Menu> Preprocessor> Element Type> Add/ Edit/ Delete. При выборе данного пункта появится диалоговое окно Element Types (Рисунок 55).

Нажимаем кнопку Add… .

В окне Library of Element Types выбираем элемент Thin Shell 163.

Нажимаем OK.

Нажимаем кнопку Options… .
В окне Shell163 element type options (Рисунок 78)	устанавливаем
следующие опции:
1. Первый	раздел	окна	Shell163 element type options.
«По умолчанию»	в	раскрывающемся	списке
Element Formulation	выбрана опция Belychko(def). В
качестве	значения	Quadrature Rule «по умолчанию»
выбрана	Gauss (1-5 Pts) (с	наименьшим	числом	точек

интегрирования). Опция Trapezoidal выбирается когда точек интегрирования более 5, но менее 100. Это касается прежде всего многослойных композиционных материалов

2.Второй раздел Layered Composite Mode. «По умолчанию» выбрана Non-Composite. Это означает, что элемент не будет иметь слоев.

3. Третий раздел Integration rule ID. «По умолчанию» выбрана опция 0. Если в первом разделе в качестве значения Quadrature Rule выбрана опция Trapezoidal и во втором разделе выбрана опция Composite, то в третьем разделе Integration rule ID следует указать целое число кратное количеству слоев.

Нажимаем OK в окне Shell163 element type options. В окне Element Types нажимаем Close.

Рисунок 78. Окно SHELL163 element type options

Опции Element Formulation окна Shell163 element type options

Использование первой опции Hughes-Liu в раскрывающемся списке Element Formulation эффективно в случае больших деформаций конструкции. Однако выбор опции существенно увеличивает вычислительные затраты.

Использование опции Belytschko(def) позволяет наиболее быстро решать динамические задачи. Не следует использовать данную опцию при больших деформациях, т.к. она не точно отображает сильное изменение геометрии конечно-элементной сетки в отличие от первой опции Hughes-Liu.

Третья опция BCIZ Triangular основана на использование теории Kirchhoff. Использование опции несколько увеличивает время решения задачи.

Четвертая опция C0 Triangular основана на теории оболочек MindlinReissner.

Использование пятой опции Belytschko-Membr позволяет не рассматривать изгибную жесткость конструкции (безмоментная теория).

З а м е ч а н и е!	Кроме того, можно выбрать опции	S/R Hughes-Liu,
S/R co-rotational,	Belytschko-Levt’n, Full Integr Membr,	Belytschko-Wong,
Fast Hughes-Liu,	Full Intgr Shell, которые являются	модификациями

предыдущих опций.

Опции 1, 2, 6-12 могут быть использованы в случае переключения элементов (и соответственно решателя) Explicit-to-Implicit. Опции 10 и 12 рекомендуются для пластического формообразования (штамповки).

Если используется модель гиперэластичного материала для элементов SHELL163, то вне зависимости от опции, выбранной пользователем, в раскрывающемся меню Element Formulation будет использована модификация опции Belytschko-Tsay.

Последовательность действий при определении констант оболочечных элементов Thin Shell 163

Выбираем пункт меню Main Menu> Preprocessor> Real Constants. При выборе данного пункта появится диалоговое окно Real Constants (Рисунок 57). Далее:

Нажимаем кнопку Add… .

В появившемся окне Element Type for… выбираем (если элементов несколько) Shell163. Нажимаем OK.

В первом окне Real Constant Set Number 1, for Shell163

указываем номер («по умолчанию» будет всегда стоять номер по порядку). Нажимаем OK.

Во втором окне Real Constant Set Number 1, for Shell163 следует указать переменные SHRF, NIP, T1, T2, T3, T4:

1.SHRF — коэффициент сдвига. Рекомендуемое значение 5/6, «по умолчанию» (если пусто) значение равно 1.

2.NIP (число точек интегрирования, если окно пусто или введен 0, то «по умолчанию» значение 2), если ранее выбраны опции композиционного материала, то желательно чтобы было выполнено, чтобы NIP было кратно количеству слоев.

3.T1 — толщина оболочки в узле I (Рисунок 77).

4.T2 — толщина оболочки в узле J (Рисунок 77).

5.T3 — толщина оболочки в узле K (Рисунок 77).

6.T4 — толщина оболочки в узле L (Рисунок 77).

Нажимаем OK.

В окне Real Constants нажимаем Close.

Особенности определения свойств многослойных

композиционных материалов

В окне Shell163 element type options (Рисунок 78) устанавливаем следующие опции:

Первый раздел. В раскрывающемся списке Element Formulation необходимо выбрать одну из опций (см. предыдущий пункт). В качестве значения Quadrature Rule «по умолчанию» выбрана Gauss (1-5 Pts) (с наименьшим числом слоев). Опция Trapezoidal выбирается когда слоев более 5, но менее 100.

Следует	выбрать	Composite	во	втором	разделе
Layered Composite Mode.

Следует указать целое число кратное количеству слоев в композиционном материале третьем разделе Integration rule ID.

Нажимаем OK в окне Shell163 element type options. В окне Element Types нажимаем Close.

Во втором окне Real Constant Set Number 1, for Shell163 следует указать переменные SHRF, NIP, T1, T2, T3, T4:

SHRF — коэффициент сдвига. Рекомендуемое значение 5/6, «по умолчанию» (если пусто) значение равно 1.

NIP – опции композиционного материала, желательно чтобы было выполнено, что NIP кратно количеству слоев (параметр из окна Shell163 element type options, см. выше).

T1 – суммарная толщина композиционной слоистой оболочки в узле I (Рисунок 79).

T2 — суммарная толщина композиционной слоистой оболочки в узле J (Рисунок 79).

T3 — суммарная толщина композиционной слоистой оболочки в узле K (Рисунок 79).

T4 — суммарная толщина композиционной слоистой оболочки в узле L (Рисунок 79).

Нажимаем OK.

Рисунок 79. Схема элемента Shell163

В третьем окне Real Constant Set Number 1, for THIN SHELL163

(Рисунок 80) необходимо указать для каждой точки интегрирования следующие параметры:

1.Относительную координату S(i) точки интегрирования с номером i.

2.Весовой коэффициент WF(i), т.е. толщина связанная с точкой интегрирования, деленная на действительную толщину оболочки.

3.ID модели материала для каждого слоя в отдельности параметр MAT(i).

В окне Real Constants следует нажать кнопку Сlose.

Рисунок 80. Вид третьего окна Real Constant Set Number 1, for Shell163

Определение коэффициентов вязкости и демпфирования. Корректировка формы деформированных элементов

Определение объемной вязкости

Main Menu> Preprocessor> Material Props> Bulk Viscosity

Данный пункт меню предназначен для глобального переопределения значений коэффициентов объемной вязкости, заданных «по умолчанию». Это может быть рекомендовано для решения задач с учетом воздействия ударных волн (в частности, исследования фронта ударной волны).

При использовании указанного пункта появляется окно Bulk Viscosity

(Рисунок	81),	в	котором	необходимо	заполнить	поля
Quadratic Viscosity Coefficient (значение «по умолчанию» равно	1.5) и
Linear Viscosity Coefficient	(значение	«по умолчанию»	равно 0.06). Для
переопределения значений необходимо нажать OK.

Рисунок 81. Окно Bulk Viscosity
Определение демпфирования
Main Menu> Preprocessor> Material Props> Damping
При использовании данного пункта меню	появляется	окно
Damping Options for LS-DYNA Explicit (Рисунок 82),	в котором	можно

определить номер «группы» (part) модели с помощью раскрывающегося списка PART number, номер кривой определяющей демпфирование «группы» Curve ID. Данный параметр принимает значения:

0 — используется постоянный коэффициент демпфирования,

заданный переменной System Damping Constant.

N (N > 0) — номер созданной кривой (в этом случае значение

System Damping Constant не используется). Усилие демпфирования f определяется уравнением f = -d(t)m v, где d(t) – коэффициент демпфирования, как функция времени, m – масса, v — скорость.

После определения параметров необходимо нажать кнопку OK.

Рисунок 82. Окно Damping Options for LS-DYNA Explicit

Глобальный параметр корректировки искаженной формы элементов

Main Menu> Preprocessor> Material Props> Hourglass Ctrls> Global

При использовании данного пункта меню появляется окно Hourglass Controls (Рисунок 83), в котором можно переопределить значение параметра Set Hourglass Coefficient to установленного «по умолчанию» (0.1). Использование значения более 0.15 может привести к неустойчивости решения.

Рисунок 83. Окно Hourglass Controls

Локальные параметры корректировки искаженной формы элементов

Main Menu> Preprocessor> Material Props> Hourglass Ctrls> Local

При использовании данного	пункта	меню	появляется	окно
Define Hourglass Material Properties	(Рисунок	84),	в котором	можно

переопределить номер материала модели (Material Reference Number), задать тип подавления искаженной формы элементов (VAL1 Hourglass control type), который может принимать 6 значений:

0, 1 — Стандартная форма вязкости LS-DYNA («по умолчанию»).

2— форма вязкости Flanagan-Belytschko.

3— форма вязкости Flanagan-Belytschko с точным интегрированием по объему для объемных элементов.

4— форма жесткости по Flanagan-Belytschko.

5— форма жесткости по Flanagan-Belytschko с точным интегрированием по объему для объемных элементов.

За м е ч а н и е! Для объемных элементов типы подавления различаются, но для оболочечных элементов опции 1-3 и 4, 5 идентичны.

Кроме того, в окне Define Hourglass Material Properties необходимо определить следующие параметры:

Коэффициент подавления искажения формы VAL2 Hourglass coefficient. «По умолчанию» его значение равно 0.1, значения данного параметра больше чем 0.15 могут приводить к неустойчивости.

Квадратичный	коэффициент	объемной	вязкости
VAL3 Quadratic bulk viscosity coefficient равен «по умолчанию»
1.5.
Линейный	коэффициент	объемной	вязкости

VAL4 Linear bulk viscosity coefficient равен «по умолчанию» 0.06.

Коэффициент подавления при изгибе оболочечного элемента

VAL5 Hourglass coefficient for shell bending	равен
«по умолчанию» значению VAL2.

Коэффициент подавления при искривлении поперечного сечения оболочки VAL6 Hourglass coefficient for shell warping равен

«по умолчанию» значению VAL2.

После определения параметров необходимо нажать кнопку OK.

Рисунок 84. Окно Define Hourglass Material Properties

Локальные системы координат

Создание локальной системы координат

Main Menu> Preprocessor> Material Props> Local CS> Create Local CS

Main Menu> Preprocessor> LS-DYNA Options> Constraints> Apply> Local CS> Create Local CS

Main Menu> Solution> Constraints> Apply> Local CS> Create Local CS

При использовании одного из этих пунктов меню появляется окно

Create Local CS at Specified Location (Рисунок 85), в котором можно:
Определить	номер	координатной	системы	в	поле

CID Coordinate system ID Number.

В полях X1, Y1, Z1 задаются X, Y, и Z координаты точки на локальной X-оси.

В полях X2, Y2, Z2 задаются X, Y, и Z координаты точки на локальной XY плоскости.

В полях X3, Y3, Z3 задаются X, Y, и Z координаты точки начала локальной системы координат («по умолчанию» X3 = 0, Y3 = 0, Z3 = 0).

Для подтверждения необходимо нажать OK.

З а м е ч а н и е! Локальные системы координат могут применяться для определения свойств ортотропных материалов.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]

Источник

Здравствуйте!

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

Итак, вопросы:

1. В ANSYS/LS-DYNA работаю на двух-ядерном процессоре Core 2duo. При расчете оба процессора нагружаются лишь на 50%. Вопрос, как вы уже догадались в следующем — как увеличить нагрузку на процессоры, чтобы сократить время расчета?

2. Пожалуйста объясните смысл параметра «Shear factor» при задании «Real Constant» для элемента Shell163 и (там же) на что влияет параметр «Thiсkness» («Толщина») у оболочечных элементов? (Иногда при контакте, если задать слишком большую величину «Толщины», то элемент отскакивает не долетая до поверхности.. Т.е., напрашивается вывод, что это какая-то незримая глазу оболочка..?)

3. В чем будет заключаться разница в контакте (и расчете) при создании «компонентов» («Component») и при создании «частей» («PART-s)» в ЛС-Дайне? Что предпочтительнее?

4. В свойствах некоторых моделей материала есть параметр «Hardening Parameter» («Параметр Упрочнения»), в хелпе написано, что он может принимать значения от 0 до 1. Вопрос такой — в чем разница между значением «0» и «1»?

5. «Explicit» в дословном переводе означает «Явный», «Implicit» — «Не явный». Пожалуйста,приведите на примере разницу, между двумя этими методами, чтобы понять смысл.. И, правильно ли я понимаю, что в ЛС-Дайне рассчитывается ТОЛЬКО «Explicit dynamic».

6. Скажите, где и как после динамического расчета, посмотреть деформацию деформируемого элемента в размерностях принятой Системы Единиц? Например по оси «Х»?

7. Я не нашел в ЛС-Дайне «компонентных элементов», т.е. когда один слой состоит из одного элемента (напр. «Солид»), а второй слой из «Шелла» или наоборот — но в двоем они образуют один элемент, который деформируется послойно? Могу ли я вообще так сделать в ЛС-Дайне?

8. В чем будет заключаться разница, при задании скорости на «узлы» («Nods») элемента или же на «части» (PART-s) элемента?

Первая волна вопрос перед вами!

Источник

Форумы
Динамика и прочность
Явная динамика
Внезапное прекращение работы решателя LS-DYNA без ошибок и предупреждений

18 апреля, 2019 — 11:36

vikotis

Здравствуйте, в процессе расчета модели закрывается командное окно с LS-DYNA, причем внезапно и без предупреждения и ошибок. Такое происходит только на одной задаче, если уменьшить шаг через restart то считает дальше. Как можно этого избежать если вдруг придётся поставить большую задачу на выходные?

P.S. шаг рассчитывал по рекомендациям www.dynasupport.com

1102 просмотра
цитата
ответить

18 апреля, 2019 — 12:06

Rav

После прекращения счёта, просматривали ли Вы файлы ошибок (*.err) — там есть какие то комменатрии? или файл output?

Можете при ответе прикрепить k-файл в архиве.

цитата
ответить

18 апреля, 2019 — 13:05

vikotis

В том то и дело что там все выглядит так как будто просто завершили процесс, хотя на самом деле это не так.

Вот так выглядит конец файл messag.

цитата
ответить

18 апреля, 2019 — 14:46

Rav

Шаг с начала счета у Вас не меняется ( стабильно 9е-7)?

Если результаты расчета загрузить в LS PrePost, видны какие нибудь баги на последнем расчетном шаге?

цитата
ответить

18 апреля, 2019 — 14:47

Rav

Выходит что расчет ложится почти в самом начале, решатель просчитал всего 0,5 мс

цитата
ответить

18 апреля, 2019 — 14:48

Rav

Цитата:

Выходит что расчет ложится почти в самом начале, решатель просчитал всего 0,5 мс

даже меньше — 0,05 мс

цитата
ответить

18 апреля, 2019 — 14:49

Rav

как у Вас задана нагрузка? можете привести параметры *Define_curve ?

цитата
ответить

18 апреля, 2019 — 18:37

vikotis

В качестве нагружающего устройства стандартное Silverado 2017 с начальной скоростью 60 км/ч. Машина врезается в барьер.

Изначально карта контроля шага расчета имеет такой вид.

Далее, после «выбрасываний» через restart шаг уменьшал до -5e-8. Сразу это помогало и расчет потиху «двигался вперед», но теперь уже даже это не помогает.

цитата
ответить

18 апреля, 2019 — 18:42

vikotis

На последнем шаге расчета, как и на остальных, никаких «артефактов» в PrePost нет, да и чисто визуально все весьма правдоподобно)

цитата
ответить

19 апреля, 2019 — 12:56

#10

Rav

Добрый день!

Если правильно понял, то речь идет о задаче из соседней темы — про взаимодейсвтие твердого тела с тросом, k-файлы для анализа взяты оттуда.

Ваша задача считается, по времени перевалили за шаг, на котором у Вас слетало:

Пока назвать конкретную причину сложно, почему задача у Вас слетает. Самому привычнее работать с k-файлом напрямую, карта *CONTROL_TERMINATION немного смутила лишними параметрами.

Вот как она задана в Вашем к-файле

Предлагается в качестве первого шага отладки — почистить эту управляющую карту до состояния представленного на рисунке (оставить только первый параметр — время окончания счета)

В решении на первых шагах все хорошо

Но с увеличением времени счета, происходит обрыв тросов

Такое решение и требовалось получить или обрыв тросов нежелателен?

цитата
ответить

vikotis

Да, это происходило в схожей задаче с тросом, однако не в этой (https://cae-club.ru/forum/modelirovanie-stalnogo-trosa-v-ls-dyna). Там была проблема с потерей контакта, однако она не вылетала внезапно. К сожалению у меня так и не получилось настроить контакт *CONTACT_AUTOMATIC_BEAMS_TO_SURFACE таким образом чтобы он не слетал в процессе расчета.

Если вернуться проблеме с «вылетом» то меня больше интересует вопрос почему решатель не выводит в log файлы причину остановки расчета в общем, а не для конкретного случая. Пересмотрев свою модель несколько раз и сделав часть верификационных расчетов я пришел к выводу что расчет может слеть из-за проникновения тела в тело (BEAM в SHELL), добавил еще несколько контактов *CONTACT_AUTOMATIC_BEAMS_TO_SURFACE и изменил часть настроек, запустил расчет. Однако хотелось бы чтобы в log файлах решатель выводил причину остановки решения, чтобы можно было оперативно отладить модель в случае ошибки, а не разбирать ее по узлам.

цитата
ответить

1
2
вперёд

Добавить комментарий

Источник

Ansys / LS-DYNA Help Forum

5
Posts

3
Users

1
Likes

1,541
Views

February 26, 2021 11:19 am

(@felix_fedi)

Student

Joined: 2 years ago

Hi guys, i’m having problems with lsdyna. No matter what problem I input, lsdyna only runs for 1 second and then ends, like the codes i show below. I don’t think it’s a keyword problem, did you guys have this problem before and can help me out?

termination time reached
1 t 9.5256E-07 dt 9.53E-07 write d3dump01 file 02/26/21 12:07:03
1 t 9.5256E-07 dt 9.53E-07 write d3plot file 02/26/21 12:07:03

(@marianadu)

Unpaid Intern

Joined: 2 years ago

Hi Felix,

I know you said you don’t think it was a keyword problem, I am just wondering if you have tried running a working example to make sure?

For example, you could one of the examples here: https://www.dynaexamples.com/implicit/basic-examples . You can download, for example the keyword file «Basics I» and you should be able to see it print 10 or so d3plots on your screen.

If I saw something like that, my first steps would be to check that the keywords *Control_Termination, *Database_Binary, and the DT0 (initial time-step size) had the required inputs and were consistent.

Cheers,

Источник

Content:

Error message (Segmentation Violation)
Checknan option (Segmenation Violation)
Error message (MPI)

Error message (Segmentation Violation )

>>>>>  Process 18  <<<<<
 >>>>> Signal 11 : Segmentation Violation <<<<<

wikipedia: http://en.wikipedia.org/wiki/Segmentation_fault

Possible procedure for the debugging of the LS-DYNA Error:

1. Is it possible to reproduce the LS-DYNA Error?

Run the identical input with identical configurations

Compare the two runs: diff between both NODOUT and both GLSTAT files

For MPP Run: extract NODOUT from the BINOUT* via l2a

—> No comparability: Hardware problem? (cpu, memory, disk, network, …)

.

2. Identify the LS-DYNA Error

activate checknan option (*CONTROL_SOLUTION, ISNAN=1)
Check the energy progression (GLSTAT, MATSUM, if necessary SLEOUT)
- Can irregularities be detected immediately before the termination?
- Which structural parts are affected?

Do elements fail immediately before the termination?
- check-failed –l mes*
- check-failed –pid –l mes* (Which PIDs are affected?)

Which are the last messages in STDOUT / STDERR / mes00xx?
Are there any problems with viscoplasticity (VP=1, “…did not converge”)?
- check-hsp d3hsp –xcon mes*

Have unreasonable contact settings been applied? (e.g. large penalties- or thickness scaling etc.)
- check-hsp d3hsp –cfac

Can initial penetrations be detected? How many? How large?
- e.g. check-c mes* -buc –typ –typ

Are there any problems with spot welds or tied contacts? (e.g. untied nodes = cold spot welds)
- e.g. check-c mes* -select tied

Within which phase the error occurs?

Keyword processing, initialization, decomposition 1-3, solution

3. Debug-Run

Switch-off „Plotcprs“ (velocity, stress, strain … available)
Termination Cycle = <Error Cycle> — 1
investigate last plot state w.r.t. max. velocities, max. plastic strains, etc..

Checknan option (Segmentation Violation)

How can I activate the checknan option?

Command-Line: checknan=1
Keyword: *CONTROL_SOLUTION, ISNAN=1

Which output is produced by checknan option?

Output to screen (e.g. jobout (STDOUT)):

*** NaN detected.  Please check message file from processor for detail.
*** Error NaN detected on processor #         30  <<<

Output to messag-files (e.g. grep ‘out-of-range‘ mes0030):

*** termination due to out-of-range forces
    number of nodes has out-of-range forces       364
    Node list:
    37502068  37502162  37502192  37502206  37502357 ....

Visualization of output by checknan option (nodes)?

DYNA Tool: check-failed –xnan –a4nan
Animator: Read in the Session-File ´fail-ide.ses´

What does NaN mean?

NaN („Not a Number“): http://en.wikipedia.org/wiki/NaN

Node forces/ -moments/ -velocities with NaN occur due to different reasons – e.g. division by zero

http://www.dynasupport.com/howtos/general/not-a-number-nan-1

check-hsp, check-failed, check-c, plotcprs are DYNAmore tools (free for customers: http://www.dynamore.de/tools)

Error message (MPI)

MPI Application rank 0 exited before MPI_Finalize() with status 13
 forrtl: error (78): process killed (SIGTERM)

„…MPI_Finalize() with status 13“ is an MPI message due to a possible LS-DYNA Error
- Is not important for debugging of the LS-DYNA error.

gp 04/14

Источник

Bhupesh Katiyar

unread,

Aug 10, 2021, 6:52:00 PM8/10/21

to LS-DYNA2

Dear all,

I am trying to use *MAT_260A material card. Can anyone please help me find the reason for the following errors? All the Part IDs, sec IDs and MAT IDs are properly defined as other material card like *MAT_36, *MAT_122, *MAT_133 are running smoothly with the same model. Please help. Thanks in advance.

*** Error 10246 (KEY+246)

line contains improperly formatted data

reading *MAT_STOUGHTON_NON_ASSOCIATED_FLOW cards.

At line# 473 of file

C:UsersBHUPESHDesktop1_DE_NAF4_nafEDDAnaf.k

——————————————————————————

3.0

——————————————————————————

*** Error 10450 (KEY+450)

in keyword command:

0.0 0.0 0.0 0.0 0.0 0.0

At line# 475 of file

C:UsersBHUPESHDesktop1_DE_NAF4_nafEDDAnaf.k

*** Warning 10434 (KEY+434)

The following lines will be skipped:

1. 0.0 1.0 0.0 0.0 0.0 0.0

expanding memory to 640252 d 2007299

expanding memory to 794368 d 2007299

expanding memory to 819464 d 2007299

expanding memory to 845150 d 2044943

*** Error 10157 (KEY+157)

MAT 7 is not found

expanding memory to 868098 d 2044955

expanding memory to 1360480 d 2044955

expanding memory to 1589295 d 2044955

expanding memory to 1614686 d 2044955

expanding memory to 1640077 d 2044955

*** Error 10304 (KEY+304)

CHECKING MATERIAL INPUT Part ID= 5

PART ID 5 with

SECTION ID 1 and

MATERIAL ID 7 does not exist.

This is PART 4 in the order of input.

*** Error 10133 (KEY+133)

input data failed with: 7 errors

Warm Regards

Bhupesh

l…@schwer.net

unread,

Aug 10, 2021, 8:43:49 PM8/10/21

to Bhupesh Katiyar, LS-DYNA2

It would be helpful if you include the input lines for MAT_260A especially Lines 473 & 475 —len

Bhupesh Katiyar

unread,

Aug 10, 2021, 9:44:59 PM8/10/21

to LS-DYNA2

Hi Len

I am including the input lines for MAT_260A for your reference.

464. *MAT_STOUGHTON_NON_ASSOCIATED_FLOW_XUE_TITLE
465. NAF
466. $# mid ro e pr r00 r45 r90 sig00
467. 77.83000E-9 210000.0 0.3 1.58 1.42 1.7 161.42
468. $# sig45 sig90 sig_b lcids lcidv scale
469. 163.09 165.08 162.06 8 9 0.03
470. $# ef0 plim q gama m
471. 1.0 1.0 1.0 1.0 1.0
472. $# aopt
473. 3.0
474. $# xp yp zp a1 a2 a3
475. 0.0 0.0 0.0 0.0 0.0 0.0
476. $# v1 v2 v3 d1 d2 d3
477. 0.0 1.0 0.0 0.0 0.0 0.0

Best Regards

Bhupesh

l…@schwer.net

unread,

Aug 10, 2021, 10:39:31 PM8/10/21

to Bhupesh Katiyar, LS-DYNA2

The Card 3 parameter AOPT should be input as an integer,

i.e. 3 rather than 3.0 – see the User Manual Volume II

When you correct that error, the other error at Line 475 should go away,

Bhupesh Katiyar

unread,

Aug 10, 2021, 11:29:02 PM8/10/21

to L…@schwer.net, LS-DYNA2

Thank you Len. I have edited the .k File in word pad and it worked. Can you please tell me what could be the possible reasons for these values of AOPT ( 2.0 and 3.0) ? I am using LS-PrePost version 4.6. One thing, I also noticed that element material direction is absent for MAT_260A. The element material direction was there for other material cards. Thank you.

With best regards

l…@schwer.net

unread,

Aug 10, 2021, 11:59:37 PM8/10/21

to Bhupesh Katiyar, LS-DYNA2

I am not sure why LSPP does not output integer values for AOPT for MAT_260A. It could just be an oversight by the developers, or a change in the way AOPT was input for older material models requiring material directions, e.g. *MAT_022

I know nothing about MAT_260A – I just now common input errors.

—len

Bhupesh Katiyar

unread,

Aug 11, 2021, 10:27:29 AM8/11/21

to L…@schwer.net, LS-DYNA2

Thank you Len for your valuable insight. It was of much help to me.

Bhupesh Singh Katiyar

Research Scholar

Metal Forming Laboratory

Department of Mechanical Engineering

Indian Institute of Technology Kharagpur

Kharagpur, India — 721302

James M. Kennedy

unread,

Aug 11, 2021, 6:42:30 PM8/11/21

to L…@schwer.net, Bhupesh Katiyar, LS-DYNA2

Dear Len,

There is a mixture of data type I and F for the AOPT parameter in the LS-DYNA User’s Manual.

Most use F data type; however, there is a few I data types used.

A quick glance (could have missed some) of the different data types used for selected data materials.

—————————

AOPT – type I input

*MAT_ADD_PZELECTRIC

*MAT_122_3D

*MAT_CODAM2

*MAT_KINEMATIC_HARDENING_BARLAT2000

*MAT_STOUGHTON_NON_ASSOCIATED_FLOW

*MAT_TISSUE_DISPERSED

*MAT_SHAPE_MEMORY_ALLOY

*MAT_THERMAL_CHEMICAL_REACTION_ORTHOTROPIC

—————————

AOPT – type F input

*MAT_OPTIONTROPIC_ELASTIC

*MAT_EXTENDED_3-PARAMETER_BARLAT

*MAT_ANISOTROPIC_PLASTIC

*MAT_COMPOSITE_LAYUP

*MAT_COMPOSITE_MATRIX

*MAT_HILL_3R

*MAT_HILL_3R_TABULATED

*MAT_HEART_TISSUE

*MAT_BARLAT_YLD2000

*MAT_WTM_STM

*MAT_WTM_STM_PLC

*MAT_SPECIAL_ORTHOTROPIC

*MAT_ORTHOTROPIC_SMEARED_CRACK

*MAT_VEGTER

*MAT_VEGTER_STANDARD

*MAT_VECTER_2017

*MAT_TRANSVERSELY_ISOTROPIC_CRUSHABLE_FOAM

*MAT_WOOD

*MAT_RATE_SENSITIVE_COMPOSITE_FABRIC

*MAT_CONCRETE_EC2

*MAT_FLD_3-PARAMETER_BARLAT

MAT_BARLAT_YLD2004

*MAT_KINEMATIC_HARDENING_BARLAT89

*MAT_DRY_FABRIC

*MAT_CAZACU_BARLAT

*MAT_HILL_90

*MAT_RESULTANT_ANISOTROPIC

*MAT_ANISOTROPIC_THERMOELASTIC

*MAT_RC_SHEAR_WALL

*MAT_VISCOELASTIC_LOOSE_FABRIC

*MAT_MICROMECHANICS_DRY_FABRIC

*MAT_PML_OPTIONTROPIC_ELASTIC

*MAT_MOHR_NON_ASSOCIATED_FLOW

*MAT_LAMINATED_FRACTURE_DAIMLER_PINHO

*MAT_LAMINATED_FRACTURE_DAIMLER_CAMANHO

*MAT_TABULATED_JOHNSON_COOK_ORTHO_PLASTICITY

*MAT_EIGHT_CHAIN_RUBBER

*MAT_PAPER

*MAT_CF_MICROMECHANICS

*MAT_ANISOTROPIC_HYPERELASTIC

*MAT_THERMAL_ORTHOTROPIC

*MAT_THERMAL_ORTHOTROPIC_TD

*MAT_THERMAL_ORTHOTROPIC_TD_LC

*MAT_THERMAL_USER_DEFINED

—————————

Sincerely,

James M. Kennedy

KBS2 Inc.

August 11, 2021

Bhupesh Katiyar

unread,

Aug 12, 2021, 9:09:59 AM8/12/21

to James M. Kennedy, L…@schwer.net, LS-DYNA2

Dear James and Len,

Thank you for sharing the data types used with AOPT for the different material models. It will certainly help all of us. I want to ask one more question regarding the use of MAT_260 A. Sometimes I get a warning message plasticity algorithm did not converge for MAT_260A. I tried changing the time step still I am getting the same warning message. What could be the possible reason for the same? If you can please share your thoughts on this it will be great. Thank you.

Warm regards

Bhupesh

l…@schwer.net

unread,

Aug 12, 2021, 5:07:09 PM8/12/21

to Bhupesh Katiyar, James M. Kennedy, LS-DYNA2

Apparently this constitutive model uses an iterative technique, e.g. Newton-Raphson, to solve for some needed quantity, e.g. plastic strain increment; this happens with the Johnson-Cook model MAT015.

In the case of MAT015, if the strain increment is too large, the iteration scheme exceeds the internally fixed number of iteration before converging to the internally set tolerance. The only way to decrease the strain increment is to use a smaller time step.

As I said before, I know nothing about MAT_260A, only common error messages.

Bhupesh Katiyar

unread,

Aug 12, 2021, 7:27:12 PM8/12/21

to L…@schwer.net, James M. Kennedy, LS-DYNA2

Thank you Len. Keep up the good work. Cheers!

Источник

LS-DYNA является одним из моих самых любимых инструментов численного моделирования. Однако и у нее есть свои пределы.

Разработанная доктором Джоном О. Холквистом (John O. Hallquist) в Ливерморской национальной лаборатории имени Лоуренса в 1976 г.[wiki], она изначально предназначалась для моделирования воздействия ядерного взрыва. По прошествии множества лет развития ее эффективный и быстрый код приспособили для решения огромного числа задач механики твердых тел и жидкостей.

LS-DYNA очень сурова: у него вообще нет графического интерфейса, что уж там говорить о приятном пользовательском опыте и простоте работы. Пере-/постпроцессоры для нее редко покрывают полный функционал решателя или отличаются таким же страшным интерфейсом. Данное программное обеспечение не для слабых духом. Однако, для тех, кто смог перешагнуть через все трудности в освоении методик работы открываются поистине неограниченные возможности моделирования физических процессов. Одно только перечисление реализованных методов вызывает уважение:

FEM (Finite element method, метод конечных элементов)
BEM (Boundary element method, метод граничных элементов)
DEM (Discrete element method, метод дискретных элементов)
XFEM (Extended FEM, расширенный метод конечных элементов)
ALE (Arbitrary Lagrange Eulerian, обобщённый Лагранжево-Эйлеров решатель)
CFD (Computational fluid dynamics, вычислительная гидродинамика)
EFG (Element Free Galerkin, без сеточный метод Галеркина)
SPH (Smoothed Particle Hydrodynamics, гидродинамика сглаженных частиц)
EM (Electromagnetism Module, решение задач электромагнетизма)

Несмотря на всю имеющуюся мощь, во время выполнения одного из проектов я столкнулся с любопытной ошибкой: «forrtl: severe (170): Program Exception — stack overflow». Данная ошибка получилась при попытке выполнения расчета модели, содержащей чуть более 50 000 тел, большая часть из которых являлись балочными. Для справки, модели, содержащие более 10 000 тел, встречаются довольно редко. По всей видимости, разработчики не предполагали создания таких масштабных моделей. Использование самых свежих дистрибутивов с сайта разработчика не помогло решить данную проблему. Посему было решено пересобрать решатель с увеличенным размером стека.

Так как в качестве рабочей станции я использую компьютер под управлением операционных систем Windows, то и переборка будет осуществлять для этой же платформы. Первое, что понадобится для выполнения задачи, это подходящий компилятор. Сейчас LSTC говорит об официальной поддержке только Intel Fortran, так что выбор падает Intel Parallel Studio XE 2013 с Intel Compiler XE v13.1.

После установки компилятора пробуем вызвать команду конфигурирования рабочего окружения:

«%ProgramFiles(x86)%IntelComposer XE 2013bincompilervars_arch.bat» intel64

Получаем следующий ответ:

«ERROR: Visual Studio 2008, 2010, 2012 or 11 is not found in the system».

Система четко дает понять, что ей нужна MS Visual Studio, причем почти любая. В данном примере выбираем Visual Studio 2010 Pro. После установки VS вызываем командную строку конфигурирования окружения компилятора снова:

«%ProgramFiles(x86)%IntelComposer XE 2013bincompilervars_arch.bat» intel64

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

Intel® Parallel Studio XE 2013
Copyright © 1985-2012 Intel Corporation. All rights reserved.
Intel® Composer XE 2013 (package 089)
Setting environment for using Microsoft Visual Studio 2010 x64 tools.

Дополнительным подтверждением служит успешное выполнение команды вызова компилятора Fortran:

ifort

Команда распознается системой, система дает ответ вроде следующего:

Intel® Visual Fortran Intel® 64 Compiler XE for applications running on Intel® 64, Version 13.0.0.089 Build 20120731
Copyright © 1985-2012 Intel Corporation. All rights reserved.

ifort: command line error: no files specified; for help type «ifort /help»

Теперь имеется готовая среда для пересборки LS-DYNA. Конечно, для человека, разбирающегося в программировании, данные действия не новы и не удивительны, но специалист в LS-DYNA не всегда обязан иметь навыки программирования.

Для дальнейшей работы нам потребуется архив «LS-DYNA User-Defined Materials Library», из которого будет собираться новый исполняемый файл. Данный архив может быть скачан с официальных ресурсов разработчика при наличии действующей технической поддержки:

www.lstc.com/download/ls-dyna_(win)
ftp.lstc.com/user/pc-dyna

В наших предпочтение всегда отдается решателю с двойной точностью, поэтому будем использовать соответствующий набор: ls-dyna_smp_d_R711_winx64_ifort131_lib.zip

За размер стека LS-DYNA отвечает значение параметра «-F» для LFLAGS указанное в makefile. По умолчанию размер стека установлен равным 180000000. В рассматриваемом случаи стабильной работы решателя удается достичь при десятикратном увеличении размера стека. Редактируем makefile и запускаем сборку при помощи nmake. В стандартном выводе можно будет увидеть много сообщений, которые разработчик рекомендует проигнорировать. Об успехе выполнения операции сборки свидетельствует следующее завершающее сообщение:

-out:lsdyna.exe
-subsystem:console
-stack:1800000000
-force
-nodefaultlib:msvcrt.lib
-nodefaultlib:vcomp.lib
dyn21.obj
dyn21b.obj
couple2other_user.obj
libdyna.lib
libansys.lib
shell32.lib
user32.lib
comctl32.lib
comdlg32.lib
Iphlpapi.lib
gdi32.lib

Таким образом получаем новый исполняемый файл «lsdyna.exe» с увеличенным в 10 раз объемом стека. Теперь LS-DYNA готова решать еще более сложные задачи механики!

Источник