/FAIL/INIEVO

Ключевое слово формата блока. Этот критерий допускает двухэтапный отказ.: подход, разделенный на начальную фазу, в которой повреждение не влияет на стресс расчет и фаза развития повреждения, на которой может быть создано смягчение напряжений. Инициирование основано на пластической деформации как функции нескольких величин напряженного состояния.
При достижении критерия инициирования эволюция переменной повреждения при смягчении напряжений: срабатывает. Эта эволюция может быть основана на пластическом смещении в момент разрушения или от определенного значения энергии разрушения. Кроме того, форма может быть выбор между классическим линейным снижением напряжения или экспоненциальным падением напряжения. грузоподъемность. Можно использовать несколько пар критериев инициации/эволюции. объединены в одной входной карте. Этот критерий совместим как с твердыми веществами, так и с оболочках и может использоваться с нелокальной регуляризацией.

Формат

Карточка 1 – Номер критерия зарождения/развития, влияние компонентов сдвига, элемент

контроль удаления


/FAIL/INIEVO/mat_ID/unit_ID	/FAIL/INIEVO/mat_ID/unit_ID	/FAIL/INIEVO/mat_ID/unit_ID	/FAIL/INIEVO/mat_ID/unit_ID	/FAIL/INIEVO/mat_ID/unit_ID	/FAIL/INIEVO/mat_ID/unit_ID	/FAIL/INIEVO/mat_ID/unit_ID	/FAIL/INIEVO/mat_ID/unit_ID	/FAIL/INIEVO/mat_ID/unit_ID	/FAIL/INIEVO/mat_ID/unit_ID
NINIEVO	ISHEAR	ILEN					FAILIP	PTHICKFAIL	PTHICKFAIL

Карточка 2 — Прочитать NINIEVO (Число пар инициации/эволюции, в

минимум 1) карты


INITYPE	EVOTYPE	EVOSHAP	COMPTYP


TAB_ID	SR_REF	SR_REF	FSCALE	FSCALE	PARAM	PARAM


TAB_EL	EL_REF	EL_REF	ELSCAL	ELSCAL


DISP	DISP	ALPHA	ALPHA	ENER	ENER

Карточка 6 – Дополнительная строка .. csv-table:

:header: "(1)", "(2)", "(3)", "(4)", "(5)", "(6)", "(7)", "(8)", "(9)", "(10)"
:widths: 10, 10, 10, 10, 10, 10, 10, 10, 10, 10

"fail_ID", "", "", "", "", "", "", "", "", ""

Определение

Поле	Содержание	Пример единицы СИ
mat_ID	Идентификатор материала.(Целое число, максимум 10 цифр)
unit_ID	(Необязательно) Идентификатор устройства. (Целое число, максимум 10 цифр)
NINIEVO	Количество инициации/эволюции критерии.По умолчанию = 1 (Целое число)
ISHEAR	Флаг для принятия поперечного сдвига компоненты (только оболочки). = 0 (по умолчанию) Поперечные компоненты не учитываются. = 1 Учтены поперечные компоненты. (Целое число)
ILEN	Характерная длина элемента флаг формулировки. = 0 (по умолчанию) Начальная геометрическая формулировка. = 1 Начальная формулировка критического временного шага. = 2 Текущая геометрическая формулировка (для оболочек только). (Целое число)
FAILIP	Количество неудачных точек интеграции перед удалением твердого элемента. По умолчанию = 1 (Целое число)
PTHICKFAIL	Процент неудачных предыдущих слоев для удаления элемента оболочки.0.0 ≤ PTHICKFAIL ≤ 1,0По умолчанию = 0,0 (Реальное)
INITYPE	Инициирование повреждений сведено в таблицу тип критерия. 2 = 1 (по умолчанию) Пластическая деформация в зависимости от трехосности напряжений (по сравнению с скорость деформации). = 2 Пластическая деформация против влияния сдвига (в зависимости от деформации ставка). = 3 Модифицированная пластическая деформация по сравнению с основным пластиком коэффициент скорости деформации (от скорости деформации) MSFLD. = 4 Пластическая деформация в зависимости от основной скорости пластической деформации коэффициент (в зависимости от скорости деформации) FLD. = 5 Зависимость пластической деформации от параметра напряженного состояния (от скорость деформации). (Целое число)
EVOTYPE	Тип развития урона. = 1 (по умолчанию) Пластическое смещение при разрушении. = 2 На основе энергии разрушения. (Целое число)
EVOSHAP	Форма развития повреждений. = 1 (по умолчанию) Линейный = 2 Экспоненциальный (Целое число)
COMPTYP	Тип комбинации критериев (только если NINIEVO > 0). = 1 (по умолчанию) Максимум = 2 Умножение (Целое число)
TAB_ID	Таблица критериев возникновения отказа идентификатор.(Целое число)
SR_REF	Справочная скорость деформации для таблицы идентификатор.По умолчанию = 1,0 (реальный)	\([\frac{1}{s}]\)
FSCALE	Масштабный коэффициент возникновения отказа Таблица критериев. По умолчанию = 1,0 (Реальное)
PARAM	Параметр критерия возникновения отказа. Если INITYPE = 1 Не используется. Если INITYPE = 2 Параметр зависимости давления. Если INITYPE = 3 = 0,0 Прямая формулировка. = 1,0 Инкрементная формулировка Если INITYPE = 4 Не используется = 0,0 Прямая формулировка. = 1,0 Инкрементная формулировка Если INITYPE = 5 Параметр влияния трехосности. По умолчанию = 0,0 (реальное)
TAB_EL	Масштабирование размера элемента при сбое критерий инициирования.(Целое число)
EL_REF	Размер эталонного элемента для размера таблица масштабирования. По умолчанию = 1,0 (Реальное)	\([m]\)
ELSCAL	Масштабный коэффициент для размера элемента функция масштабирования. По умолчанию = 1,0 (Реальное)
DISP	Пластическое смещение при отказ.По умолчанию = 0,0 (реальный)	\([m]\)
ALPHA	Параметр экспоненциальной формы (не применимо для экспоненциальной энергетической эволюции). По умолчанию = 1.0
ENER	Энергия разрушения. По умолчанию = 0,0	\([\frac{J}{m^{2}}]\)
fail_ID	(Необязательно) Критерии отказа идентификатор. (Целое число, максимум 10 цифр)

Комментарии

The

INIEVO критерий отказа – двухэтапный отказ

критерий, основанный на пластической деформации и напряженном состоянии. Эти два шага состоят

при вычислении двух последовательных переменных ущерба:

1 st Шаг: фаза инициации повреждения, на которой

Обозначенная эволюция внутренней переменной вычисляется. Эта переменная представляет собой

чисто внутренняя ценность, не влияющая на стресс

расчет. Как только эта переменная инициации достигнет значения 1,0,

вторая стадия отказа инициируется вычислением

эволюция переменной повреждения, обозначенная

\(D\) .

2 nd шаг: переменная урона \(D\) эволюция вычисляется, генерируя

эффект смягчения напряжений вплоть до полного выхода элемента из строя

и его удаление.

Различные типы отказов

критерий инициации может быть использован в зависимости от

INITYPE значение параметра. - If

INITYPE = 1 : таблица

Карта пластической деформации при возникновении разрушения приведена по отношению к

трехосность напряжений и, опционально, скорость деформации.

\(\epsilon_{p}^{init}=f\eta,\dot{\epsilon}with\eta=−\frac{p}{\sigma_{VM}}\) Где, \(p\) Гидростатическое давление \(−Tr\sigma/3\) \(\sigma_{VM}\) стресс фон Мизеса

If INITYPE = 2 : таблица

Карта пластической деформации при возникновении разрушения приведена относительно

Переменная влияния сдвига, обозначаемая как и, необязательно, со скоростью деформации.

\(\epsilon_{p}^{init}=f\theta,\dot{\epsilon}with\theta=\frac{\sigma_{VM}+k_{s}p}{\tau}\) Где, \(k_{s}\) Параметр влияния давления \(\tau\) Максимальное напряжение сдвига \(\tau=\frac{\sigma_{major}−\sigma_{minor}}{2}\)

If INITYPE = 3 : таблица

Карта модифицированной пластической деформации при возникновении разрушения представлена с помощью

относительно основного коэффициента скоростей деформаций и, необязательно, с деформацией

ставка. Этот критерий инициирования также обозначается как

MSFLD.

\(\epsilon_{p}^{init}=f\alpha,\dot{\epsilon}with\alpha=\frac{\dot{\epsilon}_{minor}^{p}}{\dot{\epsilon}_{minor}^{p}}\approx\frac{s_{minor}}{s_{major}}\) Где,

\(s_{i}\) являются главными отклонениями стрессовая составляющая.
If INITYPE = 4 : таблица

Карта пластической деформации при возникновении разрушения приведена по отношению к

главное соотношение скоростей деформаций и, опционально, со скоростью деформации. Это

критерий инициирования также обозначается как

ФЛД.

\(\epsilon_{p}^{init}=f\alpha,\dot{\epsilon}with\alpha=\frac{\dot{\epsilon}_{minor}^{p}}{\dot{\epsilon}_{minor}^{p}}\approx\frac{s_{minor}}{s_{major}}\)
If INITYPE = 5 : таблица

Карта пластической деформации при возникновении разрушения приведена относительно

параметр напряженного состояния, обозначаемый как

\(\beta\) и, необязательно, со скоростью деформации. Этот критерий инициации

также обозначается как

ФЛД.

\(\epsilon_{p}^{init}=f\beta,\dot{\epsilon}with\beta=\frac{\sigma_{VM}+k_{d}p}{\sigma_{major}}\) .. note:
Зависимость скорости деформации, примененная к
 критерий разрушения можно использовать только с материальными законами, которые
 зависят от скорости деформации. Скорость деформации, используемая для
 конститутивный закон (общая скорость деформации, скорость девиаторной деформации
 или скорость пластической деформации), будет таким же, как и для
 критерий отказа.

Начало повреждения

переменная,

\(\omega_{D}\) вычисляется постепенно

как:

\(\omega_{D}=\sumt=0\infty\frac{\Delta \epsilon_{p}}{\epsilon_{p}^{init}}\) - If

INITYPE = 3 (MSFLD), а

используется модифицированная пластическая деформация, которая развивается только тогда, когда

\(\eta>0\) : \(\omega_{D}=\sumt=0\infty\frac{\Delta \epsilon_{p}}{\epsilon_{p}^{init}}\)

If INITYPE = 3 (MSFLD) и INITYPE = 4 (FLD) критерии,

в зависимости от

PARAM значение, прямая формулировка

можно использовать вместо инкрементного:

\(\omega_{D}=\frac{\epsilon_{p}}{\epsilon_{p}^{init}}\)

Чтобы принять во внимание

регуляризация размера элемента, таблица (

TAB_EL ) может

быть определен, описывая масштабный коэффициент для размера элемента, обозначаемого как

\(f_{size}\) карта относительно начального размера элемента \(L_{e}\) и, необязательно, переменная состояния стресса

используется в качестве критерия возникновения отказа в зависимости от

INITYPE ценить. Например, если INITYPE = 2 , \(f_{size}=f(L_{e},\theta)\) . Фактор размера элемента вводится в

переменная инициирования повреждения вычисляется как:

\(\omega_{D}=\sumt=0\infty\frac{\Delta \epsilon_{p}}{\epsilon_{p}^{init}⋅f_{size}}\)

При возникновении повреждения

переменная

\(\omega_{D}\) достигает значения 1 ,

Эволюция переменной урона

\(D\) срабатывает, вызывая смягчение напряжения.

Эта эволюция может быть основана на:

Пластическое смещение при разрушении ( DISP ), обозначенный \(u_{p}^{f}\) , if EVOTYPE = 1 .

Диссипированная энергия разрушения ( ENER ) обозначается \(G_{f}\) , if EVOTYPE = 2 .

EVOSHAP

Параметр позволяет определить форму развития повреждений. Для обоих

тип эволюции (пластическое перемещение или энергия разрушения), линейный и

можно выбрать экспоненциальную форму.

If EVOSHAP = 1 , ущерб

эволюция линейна:

If EVOTYPE = 1 : пластик

смещение при разрушении вводится непосредственно с помощью

DISP : \(\Delta D=\frac{L_{e}\Delta \epsilon_{p}}{u_{p}^{f}}\) Где, \(L_{e}\) это

начальная длина элемента.

пластическое смещение DISP при разрушении)*

If EVOTYPE = 2 : пластик

смещение при разрушении выводится из энергии разрушения

вход, который рассчитывается с помощью

ENER : \(\Delta D=\frac{L_{e}\Delta \epsilon_{p}}{u_{p}^{f}}withu_{p}^{f}=\frac{2G_{f}}{\sigma_{Y}^{0}}\) Где, \(\sigma_{Y}^{0}\) предел текучести при

срабатывание эволюции урона.

Note

Это выражение
считает, что пластическое поведение почти идеально подходит в начале повреждения, чтобы обеспечить рассеиваемая энергия равна \(G_{f}\)

.

входная энергия разрушения ENER)*

If EVOSHAP = 2 , ущерб

эволюция экспоненциальна.

If EVOTYPE = 1 : пластик

смещение при разрушении вводится непосредственно с помощью

DISP и форма экспоненты может

быть изменен с помощью

ALPHA . \(D=\frac{1−e^{−\alpha\frac{L_{e}(\epsilon_{p}−\epsilon_{p}^{0})}{u_{p}^{f}}}}{1−e^{−\alpha}}\) .. image:: images/fail_inievo_starter_r_fail_inievo_exponential_0.5mm.png

*(Рисунок 3. Экспоненциальная эволюция ущерба при

пластическое смещение DISP. при ошибке с другим параметром альфа ценность)*

If EVOTYPE = 2 : площадь

между кривой напряжения – пластического смещения

экспоненциальная функция соответствует энергии разрушения

ENER . \(D=1−e^{−\frac{E_{dis}}{G_{f}}}withE_{dis}=\intD=0D=1\sigma_{Y}⋅L_{e}\Delta \epsilon_{p}\) .. image:: images/fail_inievo_starter_r_fail_inievo_exponential_100mj.png

*(Рисунок 4. Экспоненциальная эволюция ущерба при

входная энергия разрушения ENER)*

Для того же материала

несколько пар критериев инициации/эволюции могут быть определены в одном и том же

карта ввода, установка значения

NINIEVO > 1 .

В этом случае другой параметр

COMPTYP можно использовать для

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

If COMPTYP = 1 : ущерб

переменная рассматриваемого критерия, обозначаемая как

\(D_{i}\) , сравнивается с

максимальное значение среди всех критериев с использованием

COMPTYP = 1 . \(D_{MAX}=maxD_{i},D_{MAX}withi\inN_{MAX}\) Где,

\(N_{MAX}\) это число карты инициации/эволюции с использованием COMPTYP = 1.

If COMPTYP = 2 : ущерб

переменная рассматриваемого критерия, обозначаемая как

\(D_{i}\) , накапливается мультипликативным

формула среди всех критериев с использованием

COMPTYP = 2 . \(D_{MULT}=1−\underset{i\inN_{mult}}{\prod}(1−D_{i})\) Где,

\(N_{MULT}\) это количество карта инициации/эволюции с использованием COMPTYP = 2.

Наконец, переменная глобального ущерба, используемая для расчета поврежденного объекта.

Тензор напряжений определяется как следующее максимальное значение:
\(D=maxD_{MAX},D_{MULT}\sigma=(1−D)\sigma_{eff}\) Где, \(\sigma\) Окончательный тензор поврежденных напряжений \(\sigma_{eff}\) Тензор эффективных напряжений (полученный из материального закона

после картирования возврата пластика)

FAILIP

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

Параметр

PTHICKFAIL — реальный параметр, используемый оболочкой. элементы. Если PTHICKFAIL пусто или установлено значение 0,0, значение PTHICKFAIL из свойства оболочки. Если PTHICKFAIL > 0,0, любой PTHICKFAIL значение, определенное в свойствах оболочки, игнорируется, а значение, введенное в используется эта модель отказа. Для значений PTHICKFAIL > 0,0, элементы оболочки выходят из строя и удаляются при соотношении сквозной толщины неудачные точки интеграции равны или превышают PTHICKFAIL.

ILEN

— это параметр-флаг, позволяющий вам выбирать между двумя возможными формулами для

вычисление характеристической длины элемента

\(L_{e}\) . ILEN = 0 Исходная геометрическая формулировка, где характерная длина

соответствует квадратному корню из начальной площади оболочек

\(L_{e}=\sqrt{A_{0}}\) и кубический

корень из начального объема для твердых тел

\(L_{e}=\sqrt[^{3}]{V_{0}}\) . ILEN = 1 Начальная формулировка критического временного шага, где характеристика

длина соответствует той, которая использовалась для вычисления начального

критический временной шаг элемента. Используемая формула может варьироваться в зависимости от

различные составы оболочек или твердых тел.

ILEN = 2 Текущая геометрическая формулировка, где характерная длина

соответствует квадратному корню из текущей площади для оболочек

\(L_{e}=\sqrt{A}\) .

Если нелокальный

используется регуляризация (

/NONLOCAL/MAT ), нелокальный

пластическая деформация используется для расчета переменной инициирования повреждения и

переменная повреждения (и переменная нестабильности, если она используется). Также в этом случае

начальная длина элемента

\(L_{e}\) во всех формулах заменяется нелокальным

параметр

\(L_{MAX}\) .