/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 при ударе. (логическое значение) |
Комментарии
В случае SPMD каждый основной сегмент
определенный surf_IDm, должен быть связан с элементом (возможно, с пустой элемент).
Чтобы узнать о флаге Ibag, обратитесь к опции контролируемого объема (Контролируемые объемы (подушки безопасности)).
Контактная жесткость,
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 разделять элементы с твердотельными элементами, поскольку жесткость элементов оболочки используется в контактный расчет.
Истф = 2, 3,
4 или 5 несовместимы с составом SPH.
If
I кривая = 1 , сферическая кривизна
определяется для разрыва с
node_ID 1 (центр сферы). Если Icurv = 2, цилиндрический
- кривизна определяется для зазора с node_ID1 и node_ID2 (по оси цилиндра).
If I кривая = 3 , основная поверхность
форма получается с помощью бикубической интерполяции с соблюдением непрерывности координат
и нормальное от одного сегмента к другому. В случае быстрого и значительного изменения
кривизны, эта формулировка может стать нестабильной (будет улучшена в будущем).
версия).
В случае адаптивной сетки и
I адм = 1 : Если контакт
происходит в зоне (основной стороне), радиус кривизны которой меньше размера элемента
(вторичная сторона), элемент на вторичной стороне будет разделен (если еще не на
максимальный уровень).
В случае адаптивной сетки и
I адм = 2 : Если контакт произошел
в зоне (главной стороне), радиус кривизны которой меньше NRadm, умноженный на размер элемента (вторичная сторона), элемент на вторичной стороне будет разделен (если еще не на максимальный уровень).
Если контакт происходит в зоне (основной стороне), где углы
между нормалями больше, чем
Английский адм и процент
проникновение больше, чем
P адм ,
элемент на вторичной стороне будет разделен (если еще не на максимальном уровне).
![]()
Коэффициенты
Нрадм, Падм и Англадмы используются только в том случае, если они адаптивны. сетка и Iadm=2.
Если Gapmax=0, максимального значения нет.
за пробел.
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}\) самая маленькая сторона всего основного кирпича сегменты.
Переменный зазор:
Если 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.
Деактивация граничного условия
применяется к группе вторичных узлов (grnd_IDs).
Бездействие
= 3 может создать первоначальный
энергии, если узел принадлежит пружинному элементу.
Бездействие = 6 вместо этого рекомендуется
из
Бездействие = 5 , чтобы избежать высокого
частотные эффекты в интерфейсе.
![]()
Если Fpenmax не равен нулю, жесткость узлов равна
- деактивируется, если:
\(Penetration\geFpenmax⋅Gap\)
независимо от значения Inacti.
С
Irem_gap = 2 , это позволяет иметь размер элемента меньше, чем значения зазора: .. image:: images/inter_type7_starter_r_inter_type7_master_seg.png
- alt
inter_type7_master_seg
(Рисунок 5. Вторичные узлы, удаленные от узла к поверхности контакта)
Для самоударного контакта при криволинейном расстоянии (от узла главной
сегмент к вторичному узлу) меньше, чем \(\sqrt{2}⋅Gap\) (в исходной конфигурации) этот вторичный узел не будет учитываться этим основным сегментом и не будет удалена из контакта для остальных основных сегментов.
Один узел может принадлежать двум
поверхности одновременно.
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\) Тангенциальная скорость вторичного узла относительно основного сегмента
В настоящее время коэффициенты
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}\)
Фрикционная фильтрация:
Если флаг 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 частота = частота резания.
Формулировка штрафа за трение
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 рекомендуется для
неявный и явный анализ.
Теплообмен:
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 .
Тепловое трение:
Энергия трения преобразуется в тепло, когда 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 )
Радиация:
Радиация считается
в контакте, если \(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}\) Коэффициент излучения основной поверхности
Если временной шаг вторичного узла в
этот контакт становится меньше, чем dtmin, вторичный узел удаляется из контакта, и в выходной файл печатается предупреждающее сообщение. Это значение dtmin имеет приоритет над любым интерфейсом модели. минимальный шаг по времени введен в /DT/INTER/DEL.
Когда sens_ID определен для активации/деактивации
интерфейса, Tstart и Tstop не учитываются.