/INTER/TYPE7

Интерфейс ключевых слов блочного формата TYPE7 — это многофункциональный интерфейс,

моделирование контакта между основной поверхностью и группой второстепенных узлов. Также возможно Рассмотрим теплообмен и тепловое трение.

Описание

При этом решены все ограничения, возникшие с интерфейсами TYPE3, TYPE4 и TYPE5.

интерфейс:

  • Узел может одновременно быть второстепенным и основным узлом.

  • Каждый вторичный узел может влиять на каждый основной сегмент; кроме случаев, когда это

подключен к этому сегменту.

  • Узел может влиять на более чем один сегмент.

  • Узел может воздействовать с двух сторон, на края и на

углы каждого сегмента.

  • Это быстрый алгоритм поиска без ограничений.

Основные ограничения этого интерфейса: - Шаг по времени уменьшается в случае высокой скорости удара или контактов с малым зазором. - Он не работает должным образом при использовании твердого тела с высокой скоростью удара или твердого тела с небольшим зазором. - Это не решает проблему контакта края к краю (чтобы решить эту проблему, /INTER/TYPE11

следует использовать вместе с TYPE7).

Формат

/INTER/TYPE7/inter_ID/unit_ID

/INTER/TYPE7/inter_ID/unit_ID

/INTER/TYPE7/inter_ID/unit_ID

/INTER/TYPE7/inter_ID/unit_ID

/INTER/TYPE7/inter_ID/unit_ID

/INTER/TYPE7/inter_ID/unit_ID

/INTER/TYPE7/inter_ID/unit_ID

/INTER/TYPE7/inter_ID/unit_ID

/INTER/TYPE7/inter_ID/unit_ID

/INTER/TYPE7/inter_ID/unit_ID

inter_title

inter_title

inter_title

inter_title

inter_title

inter_title

inter_title

inter_title

inter_title

inter_title

grnd_IDs

surf_IDm

истф

Это

Игап

Ибаг

Идель

Icurv

Ядм

Fscalegap

Fscalegap

Гапмакс

Гапмакс

Фпенмакс

Фпенмакс

ITIED

Стмин

Стмин

Стмакс

Стмакс

%mesh_size

%mesh_size

dtmin

dtmin

Irem_gap

Irem_i2

Вставьте, если I кривая = 1 or 2 .. csv-table:

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

"node_ID1", "node_ID2", "", "", "", "", "", "", "", ""

Обязательные поля .. csv-table:

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

"Стфак", "Стфак", "Фрич", "Фрич", "Гапмин", "Гапмин", "Тстарт", "Тстарт", "Тстоп", "Тстоп"
"IBC", "", "", "Бездействие", "ВИС", "ВИС", "VISF", "VISF", "Бумульт", "Бумульт"
"Ифрик", "Ифильтр", "частота", "частота", "Яформа", "sens_ID", "fct_IDF", "АскальФ", "АскальФ", "fric_ID"

If I Фриц > 0 .. csv-table:

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

"C1", "C1", "C2", "C2", "C3", "C3", "C4", "C4", "C5", "C5"

If I Фриц > 1 .. csv-table:

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

"C6", "C6", "", "", "", "", "", "", "", ""

If I адм = 2 .. csv-table:

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

"Нрадм", "Падм", "Падм", "Англадм", "Англадм", "", "", "", "", ""

If I тот = 1 .. csv-table:

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

"Кте", "Кте", "fct_IDK", "", "оттенок", "оттенок", "Ithe_form", "АскалеК", "АскалеК", ""
"Фрад", "Фрад", "Драд", "Драд", "Подвиги", "Подвиги", "Фитм", "Фитм", "", ""

Определение

Поле

Содержание

Пример единицы СИ

inter_ID

Идентификатор интерфейса.(Целое число, максимум 10 цифр)

unit_ID

(Необязательно) Единица измерения идентификатор.(Целое число, максимум 10 цифр)

inter_title

Название интерфейса.(Символ, максимум 100 символов)

grnd_IDs

Группа вторичных узлов идентификатор.(Целое число)

surf_IDm

Основная поверхность идентификатор.(Целое число)

истф

Флаг определения жесткости интерфейса. 3 Для SPH только Istf=0, Доступны 1 и 1000. = 0 Установите значение, определенное в /DEFAULT/INTER/TYPE7. = 1 Жесткость интерфейса вводится как Stfac. = 2 Жесткость интерфейса представляет собой среднее значение основной и вторичной жесткости. жесткость. = 3 Жесткость интерфейса — максимальная из основных и второстепенных. жесткость. = 4 Жесткость интерфейса – минимальная из основных и второстепенных. жесткость. = 5 Жесткость интерфейса — это основная и вторичная жесткость последовательно. = 1000 По умолчанию, если /DEFAULT/INTER/TYPE7 не определен Жесткость интерфейса основана только на жесткости главной стороны. (Целое число)

Это

Флаг теплового контакта. = 0 Отсутствие теплопередачи и теплового трения. = 1 Активируется теплопередача или тепловое трение. (Целое число)

Игап

Флаг опции пробела/элемента. 11 = 0 Установите значение, определенное в /DEFAULT/INTER/TYPE7. = 1 Переменный зазор варьируется в зависимости от характеристик затронутой магистрали. поверхность и воздействующий вторичный узел. = 2 Переменный зазор + коррекция масштаба вычисленного зазора. = 3 Переменный зазор + коррекция масштаба вычисленного зазора + размер сетки принимается во внимание, чтобы избежать первоначального проникновения во время самоконтакта. = 1000 По умолчанию, если /DEFAULT/INTER/TYPE7 не определено Постоянный разрыв; равен минимальному зазору Gapmin. (Целое число)

Ибаг

Флажок закрытия вентиляционных отверстий подушки безопасности в чехле контакта. Этот флаг также используется для активации частиц Sol2SPH, если соответствующий твердый элемент находится в контакте. = 0 Установите значение, определенное в /DEFAULT/INTER/TYPE7. = 1 Закрытие. = 2 По умолчанию, если /DEFAULT/INTER/TYPE7 не определено Никакого закрытия. (Целое число)

Идель

Флаг удаления узла и сегмента. = 0 Установите значение, определенное в /DEFAULT/INTER/TYPE7. = 1 Когда все элементы (4-узловые оболочки, 3-узловые оболочки, тела), связанные с удаляется один сегмент, сегмент удаляется с основной стороны интерфейс. Он также удаляется в случае явного удаления с использованием ключевого слова Radioss Engine /DEL в файле Engine. Кроме того, со вторичной стороны удаляются несвязанные узлы. интерфейс. = 2 При удалении 4-узловой оболочки, 3-узловой оболочки или твердотельного элемента соответствующий сегмент удален с основной стороны интерфейса. Это также удаляется в случае явного удаления с использованием ключевого слова Radioss Engine /DEL в Файл двигателя. Кроме того, со вторичной стороны удаляются несвязанные узлы. интерфейс. = 1000 По умолчанию, если /DEFAULT/INTER/TYPE7 не определено Никакого удаления. = -1 То же, что = 1, за исключением того, что несвязанные узлы не удаляются. со вторичной стороны интерфейса. = -2 То же, что =2, за исключением того, что несвязанные узлы не удаляются из вторичная сторона интерфейса. (Целое число)Примечание: Idel = 1 и -1 имеет более высокую стоимость процессора по сравнению с Idel = 2 и -2.

Icurv

Конверт с зазором и кривизной. 5 = 0 (по умолчанию) Никакой кривизны. = 1 Сферическая кривизна. = 2 Цилиндрическая кривизна. = 3 Автоматическая бикубическая поверхность. (Целое число)

Ядм

Вычисление флага локальной кривизны для адаптивная сетка. 6 7 8 = 0 (по умолчанию) Не активирован. = 1 Обновление интерфейса в соответствии с размером сетки. = 2 Обновление интерфейса в соответствии с размером сетки, проникновением и углом. (Целое число)

Fscalegap

Масштабный коэффициент зазора (используется только тогда, когда Igap = 2 и 3).По умолчанию = 1,0 (Реальное)

Гапмакс

Максимальный зазор (используется только тогда, когда Igap = 2 и 3).(Реальный)

\([m]\)

Фпенмакс

Максимальная доля начального проникновения. 13(Реал)

ITIED

Флаг узлов контакта связи. = 0 (по умолчанию) Узлы не связаны из-за шага по времени. = 1 Вторичный узел связан с возможным отскоком при шаге времени контакта меньше /DT/INTER/DEL или определено dtmin в контакте. = 2 Вторичный узел завязывается без возможного отскока при шаге времени контакта меньше /DT/INTER/DEL или dtmin определился в контакте. (Целое число)

Стмин

Минимальная жесткость (используется только при Istf = 2, 3, 4 или 5).(Реальный)

\([\frac{N}{m}]\)

Стмакс

Максимальная жесткость (используется только при Istf = 2, 3, 4 или 5). По умолчанию = 1030 (Реал)

\([\frac{N}{m}]\)

%mesh_size

Процент размера сетки (используется только тогда, когда Igap = 3).По умолчанию = 0,4 (Реал)

dtmin

Минимальный временной шаг интерфейса. 23 (Реал)

\([s]\)

Irem_gap

Флаг деактивации вторичных узлов, если размер элемента < величина зазора, в случае самоударного контакта. 14 = 0 Установите значение, определенное в /DEFAULT/INTER/TYPE7. = 1 По умолчанию, если /DEFAULT/INTER/TYPE7 не определено Никакой деактивации вторичных узлов. = 2 Деактивация вторичных узлов. (Целое число)

Irem_i2

Флаг для деактивации вторичного узел, если в интерфейсе TYPE2 определена одна и та же контактная пара (узлы). = 0 Установите значение, определенное в /DEFAULT/INTER/TYPE7. = 1 Вторичные узлы в связанных контактах /INTER/TYPE2 удалены. из этого контакта. = 3 По умолчанию, если /DEFAULT/INTER/TYPE7 не определено Никаких изменений во вторичных узлах.

node_ID1

Первый узел идентификатор.(Целое число)

node_ID2

Идентификатор второго узла (игнорируется, если Icurv = 1).(Целое число)

Стфак

Масштабный коэффициент жесткости, применяемый к основному сторона интерфейса (если Istf \(\ne\) 1).По умолчанию = 1,0 (Реальное)

Жесткость интерфейса (если Istf = 1). По умолчанию = 1,0. (Реал)

\([\frac{N}{m}]\)

Фрич

Кулоновское трение (если fct_IDF = 0). По умолчанию = 0,0. (Реал)

Масштабный коэффициент кулоновского трения (если fct_IDF \(\ne\) 0). 20По умолчанию = 1,0 (Реальное)

Гапмин

Минимальный зазор для ударной активации.Для значения по умолчанию, см. 10 (реальный)

\([m]\)

Тстарт

Время начала.(Реальное)

\([s]\)

Тстоп

Время для временного деактивация.(Реальная)

\([s]\)

IBC

Флаг деактивации граничных условий при ударе. (логическое значение)

Бездействие

Флаг деактивации жесткости в случае начальных проникновений. 13 = 0 Установите значение, определенное в /DEFAULT/INTER/TYPE7. = 1 Деактивация жесткости на узлах. = 2 Деактивация жесткости элементов. = 3 Измените координаты узла, чтобы избежать начальных проникновений. = 5 Разрыв является переменным со временем, и начальный разрыв корректируется следующим образом: \(gap_{0}=Gap−P_{0}\) , где \(P_{0}\) это начальное проникновение = 6 Зазор меняется со временем, но начальный зазор корректируется как (узел слегка углублённый): \(gap_{0}=Gap–P_{0}–5%⋅(Gap–P_{0})\) = 1000 По умолчанию, если /DEFAULT/INTER/TYPE7 не определено Никаких действий. (Целое число)

ВИС

Критический коэффициент демпфирования жесткость интерфейса. По умолчанию установлено значение 0,05 (реальное).

VISF

Критический коэффициент демпфирования трение интерфейса. 22 По умолчанию установлено значение 1,0 (реальное).

Бумульт

Фактор сортировки используется для ускорения алгоритм сортировки и зависит от машины. По умолчанию установлено значение 0,20. (Реал)

Ифрик

Флаг формулировки трения. 1617 Используется только если fric_ID не определен. = 0 (по умолчанию) Статический закон трения Кулона. = 1 Обобщенный закон вязкого трения. = 2 (Модифицированный) Дармстадский закон трения. = 3 Закон трения Ренара. = 4 Экспоненциальный закон трения затухания. (Целое число)

Ифильтр

Флаг фильтрации трения. 18 = 0 (по умолчанию) Фильтр не используется. = 1 Простой числовой фильтр. = 2 Стандартный фильтр -3 дБ с периодом фильтрации. = 3 Стандартный фильтр -3 дБ с частотой среза. (Целое число)

частота

Коэффициент фильтрации. Значение по умолчанию зависит от Ifiltr (реального)

Яформа

Тип формулировки штрафа за трение. 19 = 0 Установите значение, определенное в /DEFAULT/INTER/TYPE7. = 1 По умолчанию, если /DEFAULT/INTER/TYPE7 не определено Вязкая (общая) рецептура. = 2 Формулировка жесткости (приращения). (Целое число)

sens_ID

Идентификатор датчика для активации/деактивации интерфейс. 24Если определен датчик-идентификатор, активация/деактивация интерфейса основана на датчике, а не на Tstart или Tstop. (Целое число)

fct_IDF

Коэффициент трения в зависимости от температуры идентификатор функции. По умолчанию = 0 (Целое число)

АскальФ

Масштабный коэффициент абсцисс для fct_IDF.По умолчанию = 1,0 (реальный)

fric_ID

Идентификатор трения для трения определение для выбранных пар деталей. = 0 (по умолчанию) Используйте параметры трения, определенные в этом интерфейсе. ≠ 0 Используйте /FRICTION/fric_ID. (Целое число)

C1 - C6

Закон трения коэффициент.(Реальный)

См. Таблицу 1.

Кте

Коэффициент теплообмена (если fct_IDK = 0). По умолчанию = 0,0

\([\frac{W}{m^{2}K}]\)

Коэффициент теплообмена (если fct_IDK \(\ne\) 0). 20 По умолчанию = 1,0 (реальное)

fct_IDK

Идентификатор функции теплового нагрева определение обмена с контактным давлением. По умолчанию = 0 (Целое число)

АскалеК

Масштабный коэффициент абсцисс для fct_IDK.По умолчанию = 1,0 (реальный)

\([Pa]\)

оттенок

Температура интерфейса. 20(Реал)

\([K]\)

Ithe_form

Флаг формулировки теплового контакта. = 0 (по умолчанию) Обмен только между интерфейсом (постоянная температура) и оболочками (вторичный сторона). = 1 Теплообмен между всеми контактирующими деталями. (Целое число)

Нрадм

Количество элементов через 90 градусов радиус.(Целое число)

Падм

Критерии процентного содержания проникновение.По умолчанию = 1,0 (Реальное)

Англадм

Угловые критерии.(Реальные)

\([deg]\)

Фрад

Коэффициент радиации.(Реальный)

\([\frac{W}{m^{2}K^{4}}]\)

Драд

Максимальное расстояние для излучения расчет.(Реальный)

\([m]\)

Подвиги

Фактор трения вторичного нагрева. 21 (Реал)

Фитм

Коэффициент фрикционного нагрева магистрали. 21 (Реал)

Флаги для деактивации граничных условий: IBC

(1)-1

(1)-2

(1)-3

(1)-4

(1)-5

(1)-6

(1)-7

(1)-8

IBCX

IBCY

IBCZ

Определение

Поле

Содержание

Единица СИ Пример

IBCX

=1 Флаг деактивации граничного условия X при ударе. (логическое значение)

IBCY

=1 Флаг деактивации граничного условия Y при ударе. (логическое значение)

IBCZ

=1 Флаг деактивации граничного условия Z при ударе. (логическое значение)

Комментарии

  1. В случае SPMD каждый основной сегмент

определенный surf_IDm, должен быть связан с элементом (возможно, с пустой элемент).

  1. Чтобы узнать о флаге Ibag, обратитесь к опции контролируемого объема (Контролируемые объемы (подушки безопасности)).

  2. Контактная жесткость,

K рассчитывается как: Если Iстф =1: \(K=Stfac\) Если Iстф = 2, 3,

4 или 5:

\(K=\frac{max[St_{min},min(St_{max},K_{n})]}{2}\) Если Iстф =1000: \(K=\frac{K_{m}}{2}\) Где,

\(K_{n}\) рассчитывается на основе жесткости обоих основных сегментов \(K_{m}\) и жесткость вторичного узла \(K_{s}\)

:

Истф = 2,

\(K_{n}=\frac{K_{m}+K_{s}}{2}\)

Истф = 3,

\(K_{n}=max(K_{m},K_{s})\)

Истф = 4,

\(K_{n}=min(K_{m},K_{s})\)

Истф = 5,

\(K_{n}=\frac{K_{m}⋅K_{s}}{K_{m}+K_{s}}\)

\(K_{m}\)

– жесткость основного сегмента, рассчитываемая как:

Когда основной сегмент лежит на оболочке или является общим для оболочки и твердого тела: \(K_{m}=Stfac⋅E⋅t\) Когда основной сегмент лежит на твердом теле: \(K_{m}=Stfac⋅B⋅\frac{S^{2}}{V}\) Где, \(S\) Область сегмента \(V\) Объем твердого тела \(B\) Объемный модуль \(K_{s}\)

— эквивалентная узловая жесткость, рассматриваемая для интерфейса TYPE7 и рассчитывается как:

Когда узел подключен к элементу оболочки: \(K_{s}=E⋅t\) Когда узел соединен с твердотельным элементом: \(K_{s}=B⋅\sqrt[3]{V}\) Нет ограничений на значение масштабного коэффициента жесткости.

Stfac (но значение больше 1,0 может уменьшить начальный шаг по времени).

При использовании

/PROP/VOID и /MAT/VOID, свойства материала и толщину материала VOID необходимо ввести; в противном случае контактная жесткость пустые элементы будут равны нулю. Это особенно важно, если элементы оболочки VOID разделять элементы с твердотельными элементами, поскольку жесткость элементов оболочки используется в контактный расчет.

  1. Истф = 2, 3,

4 или 5 несовместимы с составом SPH.

  1. If

I кривая = 1 , сферическая кривизна

определяется для разрыва с

node_ID 1 (центр сферы). Если Icurv = 2, цилиндрический

кривизна определяется для зазора с node_ID1 и node_ID2 (по оси цилиндра).

If I кривая = 3 , основная поверхность

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

и нормальное от одного сегмента к другому. В случае быстрого и значительного изменения

кривизны, эта формулировка может стать нестабильной (будет улучшена в будущем).

версия).

  1. В случае адаптивной сетки и

I адм = 1 : Если контакт

происходит в зоне (основной стороне), радиус кривизны которой меньше размера элемента

(вторичная сторона), элемент на вторичной стороне будет разделен (если еще не на

максимальный уровень).

  1. В случае адаптивной сетки и

I адм = 2 : Если контакт произошел

в зоне (главной стороне), радиус кривизны которой меньше NRadm, умноженный на размер элемента (вторичная сторона), элемент на вторичной стороне будет разделен (если еще не на максимальный уровень).

Если контакт происходит в зоне (основной стороне), где углы

между нормалями больше, чем

Английский адм и процент

проникновение больше, чем

P адм ,

элемент на вторичной стороне будет разделен (если еще не на максимальном уровне).

images/inter_type7_starter_r_inter_type7_angladm.svg*(Рисунок3.)*
  1. Коэффициенты

Нрадм, Падм и Англадмы используются только в том случае, если они адаптивны. сетка и Iadm=2.

  1. Если Gapmax=0, максимального значения нет.

за пробел.

  1. If

Зазор мин = 0 или пусто, значение по умолчанию

рассчитывается как:

Если основными сегментами являются оболочечные и твердотельные элементы, Gapmin = min (

\(tm\)

,

\(\frac{l_{min}}{2}\)

).

Где, \(tm\) Средняя толщина основных элементов оболочки, для I зазор = 0 \(tm\) Минимальная толщина основных элементов оболочки, для I зазор = 1 , 2 , or 3 \(l_{min}\) Наименьшая длина стороны всех основных сегментов (ракушки или кирпича) Если все основные сегменты представляют собой сплошные элементы, Gapmin =

\(\frac{l_{min}}{10}\)

Где,

\(l_{min}\) самая маленькая сторона всего основного кирпича сегменты.

  1. Переменный зазор:

Если Igap =1, переменный разрыв равен

вычисляется как:

\(max[Gap_{min},(g_{s}+g_{m})]\) Если Igap =2, переменный разрыв равен

вычисляется как:

\(max{Gap_{min},min[Fscale_{gap}⋅(g_{s}+g_{m}),Gap_{max}]}\) Если Igap =3, переменный разрыв равен

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

\(max{Gap_{min},min[Fscale_{gap}⋅(g_{s}+g_{m}),%mesh_size⋅(g_{s_l}+g_{m_l}),Gap_{max}]}\) Где, - \(g_{m}\)

: разрыв основного элемента \(g_{m}=\frac{t}{2}\)

: с

\(t\) толщина основного элемента оболочки элементы

\(g_{m}=0\)

: для кирпичных элементов
  • \(g_{s}\) : разрыв вторичного узла \(g_{s}=0\)

: если вторичный узел не подключен ни к одному
элемент или соединяется только с кирпичными или пружинными элементами

\(g_{s}=\frac{t}{2}\)

: с

\(t\) наибольшая толщина элементов оболочки подключен к вторичному узлу

\(g_{s}=\frac{1}{2}\sqrt{S}\)

: для ферм и балочных элементов. С

\(S\) сечение элемента

  • \(g_{m_l}\) : длина меньшего края элемента

  • \(g_{s_l}\) : длина меньшего края элементов, соединенных с

    вторичный узел

Если вторичный узел соединен с несколькими оболочками и/или балками или

ферм, используется наибольший рассчитанный вторичный зазор.

Переменный зазор всегда находится на уровне

минимум равен Gapmin.

  1. Деактивация граничного условия

применяется к группе вторичных узлов (grnd_IDs).

  1. Бездействие

= 3 может создать первоначальный

энергии, если узел принадлежит пружинному элементу.

Бездействие = 6 вместо этого рекомендуется

из

Бездействие = 5 , чтобы избежать высокого

частотные эффекты в интерфейсе.

images/inter_type7_starter_r_inacti_flag_6.svg*(Рисунок4.)*

Если Fpenmax не равен нулю, жесткость узлов равна

деактивируется, если:

\(Penetration\geFpenmax⋅Gap\)

независимо от значения Inacti.

  1. С

Irem_gap = 2 , это позволяет иметь размер элемента меньше, чем значения зазора: .. image:: images/inter_type7_starter_r_inter_type7_master_seg.png

alt

inter_type7_master_seg

(Рисунок 5. Вторичные узлы, удаленные от узла к поверхности контакта)

Для самоударного контакта при криволинейном расстоянии (от узла главной

сегмент к вторичному узлу) меньше, чем \(\sqrt{2}⋅Gap\) (в исходной конфигурации) этот вторичный узел не будет учитываться этим основным сегментом и не будет удалена из контакта для остальных основных сегментов.

  1. Один узел может принадлежать двум

поверхности одновременно.

  1. If

fric_ID определено, контактное трение определяется в /FRICTION и

входы трения (

I Фриц , C 1 ,

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

Для формулировки трения: Если

флаг трения Ifric=0 (по умолчанию), старый используется формулировка статического трения:

\(F_{t}\le\mu⋅F_{n}\) с \(\mu\) коэффициент кулоновского трения - If

fct_ID F = 0 : Ифрик — это Кулон

трение.

\(\mu=Fric\)

  • If fct_ID F ≠ 0 : Fric становится масштабным коэффициентом Кулона.

коэффициент трения, который зависит от температуры.

\(\mu=Fric⋅f_{F}(Ascale_{F},T_{interface})\)

Пока

\(T_{interface}\) это температура, которую принимают за среднюю температуру второстепенных и основных:

\(T_{interface}=\frac{T_{secondary}+T_{main}}{2}\) Для флага Ifric > 0 новые модели трения

вводятся. В этом случае коэффициент трения задается функцией ( \(\mu=\mu (\rho,V)\)

).

Где, \(\rho\) Давление нормальной силы на основной сегмент \(V\) Тангенциальная скорость вторичного узла относительно основного сегмента

  1. В настоящее время коэффициенты

C 1 через C 6 используются для определения переменного трения

коэффициент

\(\mu\) для

новые формулы трения.

Доступны следующие составы: - I

Фриц = 1 (Обобщенный

Закон вязкого трения):

\(\mu=Fric+C_{1}.p+C_{2}⋅V+C_{3}.p⋅V+C_{4}⋅p^{2}+C_{5}⋅V^{2}\)

  • I Фриц = 2 (Модифицированный Дармстад

    закон):

    \(\mu=Fric+C_{1}.e^{(C_{2}V)}.p^{2}+C_{3}.e^{(C_{4}V)}.p+C_{5}.e^{(C_{6}V)}\)

  • I Фриц = 3 (Ренар

    закон):

    \(\mu=C_{1}+(C_{3}−C_{1})⋅\frac{V}{C_{5}}⋅(2−\frac{V}{C_{5}})ifV\in[0,C_{5}]\) \(\mu=C_{3}−((C_{3}−C_{4})⋅(\frac{V−C_{5}}{C_{6}−C_{5}})^{2}⋅(3−2⋅\frac{V−C_{5}}{C_{6}−C_{5}}))ifV\in[C_{5}C_{6}]\) \(\mu=C_{2}−\frac{1}{\frac{1}{C_{2}−C_{4}}+(V−C_{6})^{2}}ifV\geC_{6}\) Где для I Фриц = 3 : - \(C_{1}=\mu_{s}\)

    , статический коэффициент трения, должен быть \(\mu_{min}<\mu_{s}<\mu_{max}\)

    • \(C_{2}=\mu_{d}\) , динамический коэффициент трения, должен быть \(\mu_{min}<\mu_{d}<\mu_{max}\)

    • \(C_{3}=\mu_{max}\) , максимальный коэффициент трения

    • \(C_{4}=\mu_{min}\) , минимальный коэффициент трения

    • \(C_{5}=V_{cr1}\) , первая критическая скорость, должна быть > 0

    • \(C_{6}=V_{cr2}\) , вторая критическая скорость, должна быть \(>V_{cr1}\)

  • I Фриц = 4 (Экспоненциальный закон трения затухания) Фрикционный

предполагается, что коэффициент зависит от относительной скорости \(V\) поверхностей, соприкасающихся в соответствии с чтобы:

\(\mu=C_{1}+Fric−C_{1}⋅e^{−C_{2}V}\)

  1. Фрикционная фильтрация:

Если флаг Ifiltr

\(\ne\)

0, касательные силы сглаживаются с помощью фильтра:

\(F_{Tf}=\alphaF_{T}(t)+(1−\alpha)F_{Tf}(t−dt)\) Где, \(F_{Tf}\) Фильтрованная тангенциальная сила. \(F_{T}(t)\) Рассчитанная тангенциальная сила в момент времени t перед фильтрацией. \(F_{Tf}(t−dt)\) Отфильтрованная тангенциальная сила на предыдущем временном шаге \(t\) Текущее время моделирования \(dt\) Текущий временной шаг моделирования \(\alpha\) Коэффициент фильтрации Где \(\alpha\) коэффициент

рассчитано из:

  • If I фильтр = 1 : \(\alpha=X_{freq}\) , простой числовой фильтр со значением от 0 до

  • If I фильтр = 2 : \(\alpha=\frac{2⋅\pi}{X_{freq}}\) , стандартный фильтр -3дБ, с количеством временных шагов

    для фильтрации определяется как

    \(X_{freq}=\frac{dt}{T}\) , и \(T\) это период фильтрации.

  • If I фильтр = 3 : \(\alpha=2⋅\pi⋅X_{freq}⋅dt\) , стандартный фильтр -3 дБ, с X частота = частота резания.

  1. Формулировка штрафа за трение

I форма : - If

I форма = 1 (по умолчанию) вязкий

В формулировке силы трения равны:

\(F_{t}=min⁡(\muF_{n},F_{adh})\) Сила сцепления рассчитывается как: \(F_{adh}=C⋅V_{t}withC=VIS_{F}⋅\sqrt{2Km}\)

  • If I форма = 2 , жесткость

    формулировке), силы трения равны:

    \(F_{t}^{new}=min⁡(\muF_{n},F_{adh})\) В то время как адгезия рассчитывается как: \(F_{adh}=F_{t}^{old}+\Delta F_{t}\)

с \(\Delta F_{t}=K⋅V_{t}⋅dt\)

Где,

\(V_{t}\) - тангенциальная скорость вторичного узла относительно основного сегмента.

Iform = 2 рекомендуется для

неявный и явный анализ.

  1. Теплообмен:

By I тот = 1 (теплопередача

активирован) для учета теплообмена и теплового трения при контакте.

  • Если Ithe_form = 0, то теплообмен происходит между оболочкой и константой.

контакт температуры Оттенок.
  • Если Ithe_form = 1, то теплообмен происходит между всеми контактными деталями. Оттенок используется только тогда, когда Ithe_form= 0. В этом случае предполагается, что температура главной стороны равна

константа (равная Оттенку). Если Ithe_form=1, то оттенок не учитывается, т.к. будет учитываться узловая температура основной стороны.

Коэффициент теплообмена: - Если fct_IDK = 0, то Kthe – коэффициент теплообмена и теплоотдача.

обмен зависит только от поверхности теплообмена.
  • If fct_ID K ≠ 0 , затем K тот масштабный фактор и теплообмен

    будет зависеть от контактного давления:

    \(K=K_{the}⋅f_{K}(Ascale_{K},P)\)

  • Пока \(f_{K}\) это функция fct_ID K .

  1. Тепловое трение:

  • Энергия трения преобразуется в тепло, когда Ithe > 0 для интерфейса типа 7.

только.
  • Фит s и Фит m определяются как

    долю энергии трения и распределяется соответственно на вторичной стороне и

    основная сторона. Итак, в целом:

    \(Fheat_{s}+Fheat_{m}\le1.0\) Когда оба подвига

и Fheatm равны 0, преобразование энергии трения скольжения в тепло не происходит. активирован.

  • Тепло трения Q Фрич определяется: - If

    I форма = 2 (жесткость

    формулировка):

    Вторичная сторона: \(Q_{Fric}=Fheat_{s}⋅\frac{(F_{adh}−F_{t})}{K}⋅F_{t}\) Основная сторона: \(Q_{Fric}=Fheat_{m}⋅\frac{(F_{adh}−F_{t})}{K}⋅F_{t}\) ( I the_form = 1 )

    • If I форма = 1 (штраф

      формулировка):

      Вторичная сторона: \(Q_{Fric}=Fheat_{s}⋅C⋅V_{t}^{2}⋅dt\) Основная сторона: \(Q_{Fric}=Fheat_{m}⋅C⋅V_{t}^{2}⋅dt\) ( I the_form = 1 )

  1. Радиация:

Радиация считается

в контакте, если \(F_{rad}\ne0\) и расстояние, \(d_{}\)

, вторичного узла к основному сегменту
это:

\(Gap<d<D_{rad}\) Пока \(D_{rad}\) — максимальное расстояние для расчета излучения.

значение по умолчанию для

\(D_{rad}\) рассчитывается как максимум: - Верхнее значение разрыва (в момент времени 0) среди всех узлов - Наименьшая длина стороны вторичного элемента Рекомендуется не устанавливать слишком высокое значение для

\(D_{rad}\)

, что может снизить производительность Radioss Engine.

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

вычисляется как:

\(h_{rad}=F_{rad}(T_{m}^{2}+T_{s}^{2})⋅(T_{m}+T_{s})\) с \(F_{rad}=\frac{\sigma}{\frac{1}{\epsilon_{1}}+\frac{1}{\epsilon_{2}}−1}\) Где, \(\sigma=5.669\times10^{−8}[\frac{W}{m^{2}K^{4}}]\) Постоянная Стефана Больцмана \(\epsilon_{1}\) Коэффициент излучения вторичной поверхности \(\epsilon_{2}\) Коэффициент излучения основной поверхности

  1. Если временной шаг вторичного узла в

этот контакт становится меньше, чем dtmin, вторичный узел удаляется из контакта, и в выходной файл печатается предупреждающее сообщение. Это значение dtmin имеет приоритет над любым интерфейсом модели. минимальный шаг по времени введен в /DT/INTER/DEL.

  1. Когда sens_ID определен для активации/деактивации

интерфейса, Tstart и Tstop не учитываются.