/INTER/TYPE25

Ключевое слово формата блока TYPE25 — это общие узлы для: поверхность контакта с использованием метода штрафов. Штрафная жесткость постоянна и поэтому шаг по времени не изменяется.
Твердые элементы имеют нулевую толщину контактного зазора. Контактные входы можно определить как: одна поверхность, поверхность к поверхности или узлы к поверхности.
Этот контактный интерфейс может заменить интерфейс TYPE3,: ТИП5, ТИП7, ТИП19 или ТИП24.

Этот интерфейс недоступен в неявном решении.

Формат


/INTER/TYPE25/inter_ID/unit_ID	/INTER/TYPE25/inter_ID/unit_ID	/INTER/TYPE25/inter_ID/unit_ID	/INTER/TYPE25/inter_ID/unit_ID	/INTER/TYPE25/inter_ID/unit_ID	/INTER/TYPE25/inter_ID/unit_ID	/INTER/TYPE25/inter_ID/unit_ID	/INTER/TYPE25/inter_ID/unit_ID	/INTER/TYPE25/inter_ID/unit_ID	/INTER/TYPE25/inter_ID/unit_ID
inter_title	inter_title	inter_title	inter_title	inter_title	inter_title	inter_title	inter_title	inter_title	inter_title
surf_ID1	surf_ID2	истф	Это	Игап	Irem_i2		Идель	Иедж	Ипстиф
grnd_IDs		Gap_scale	Gap_scale	%mesh_size	%mesh_size	Gap_max_s	Gap_max_s	Gap_max_m	Gap_max_m
Стмин	Стмин	Стмакс	Стмакс	Igap0	Ishape	Edge_angle	Edge_angle	Stfacm	Stfacm
Стфак	Стфак	Фрич	Фрич	Тпрессфит	Тпрессфит	Тстарт	Тстарт	Тстоп	Тстоп
IBC		IVIS2	Бездействие	ВИС	ВИС
Ифрик	Ифильтр	частота	частота		sens_ID	Дстиф	Дстиф		fric_ID

Прочтите этот ввод только в том случае, если 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"

Прочтите этот ввод только в том случае, если 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", "", "", "", "", "", "", "", ""

Прочтите этот ввод только в том случае, если ИВИС2 = -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

"ViscFluid", "ViscFluid", "СигМаксАдх", "СигМаксАдх", "ВискАдфакт", "ВискАдфакт", "", "", "", ""

Прочтите этот ввод только в том случае, если 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

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

Прочтите этот ввод только в том случае, если I зазор = 5 .. csv-table:

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

"Thickness_s", "Thickness_s", "Scale_s", "Scale_s", "Thickness_m", "Thickness_m", "Scale_m", "Scale_m", "", ""

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

Поле	Содержание	Пример единицы СИ
inter_ID	Интерфейс идентификатор.(Целое число, максимум 10 цифр)
unit_ID	Идентификатор объекта.(Целое число, максимум 10 цифр)
inter_title	Интерфейс заголовок.(Символ, максимум 100 символов)
surf_ID1	Первый идентификатор поверхности. 1 (Целое число)
surf_ID2	Вторая поверхность идентификатор.(Целое число)
истф	Жесткость интерфейса флаг определения. 2 = 0 Установите значение, определенное в /DEFAULT/INTER/TYPE25. = 2 Жесткость интерфейса представляет собой среднее значение основных и вторичная жесткость. = 3 Жесткость интерфейса – максимальная из основных и вторичная жесткость. = 4 Жесткость интерфейса – минимальная из основных и вторичная жесткость. = 5 Жесткость интерфейса — это основная и вторичная жесткость. в серии. = 7 Жесткость интерфейса основана только на стабильности состояние. 15 = 1000 По умолчанию, если /DEFAULT/INTER/TYPE25 нет определенный Жесткость интерфейса рассчитывается только по основной стороне. жесткость. (Целое число)
Это	Флаг теплопередачи. = 0 Нет теплопередачи = 1 Передача тепла активирована. Теплообмен недоступен для края край.(Целое число)
Игап	Флаг опции пробела/элемента. 3 = 0 Используйте значение, определенное в /DEFAULT/INTER/TYPE25. = 1 По умолчанию, если /DEFAULT/INTER/TYPE25 нет определенный Переменный зазор варьируется в зависимости от характеристик воздействующая основная поверхность и воздействующая вторичная узел. = 2 Переменный разрыв (аналогично Igap=1) и деактивация вторичного узлы, если размер элемента < значение зазора, в случае контакт с самовоздействием. = 3 Переменный зазор, при котором размер сетки (определенный в %mesh_size) считается избегайте первоначальных проникновений при контакте с самим собой. = 5 Постоянный разрыв. (Целое число)
Irem_i2	Деактивация флага для вторичный узел, если в интерфейс ТИП2. = 0 Используйте значение, определенное в /DEFAULT/INTER/TYPE25. = 1 По умолчанию, если /DEFAULT/INTER/TYPE25 нет определенный Вторичные узлы в /INTER/TYPE2 привязаны контакты удаляются из этого контакта. = 3 Никаких изменений во вторичных узлах.
Идель	Флаг удаления узла и сегмента. = 0 Используйте значение, определенное в /DEFAULT/INTER/TYPE25. = 1 Когда все элементы (4-узловые оболочки, 3-узловые оболочки, твердые вещества), связанные с одним сегментом, удаляются, сегмент удален с основной стороны интерфейса. Он также удаляется в случае явного удаления с помощью Ключевое слово Radioss Engine /DEL в двигателе файл. Кроме того, из списка удаляются несвязанные узлы. вторичная сторона интерфейса. = 2 Когда 4-узловая оболочка, 3-узловая оболочка или сплошной элемент удаляется, соответствующий сегмент удаляется из основная часть интерфейса. Он также удаляется в случай явного удаления с использованием ключевого слова Radioss Engine /DEL в файле Engine. Кроме того, из списка удаляются несвязанные узлы. вторичная сторона интерфейса. = 1000 По умолчанию, если /DEFAULT/INTER/TYPE25 нет определенный Никакого удаления.
Иедж	Варианты краевого контакта. Контакт происходит между главными и второстепенными краями, которые автоматически извлекается из surf_ID1 и surf_ID2. Острые кромки внешних твердых граней определяются с помощью угла Edge_angle. = 0 Установите значение, определенное в /DEFAULT/INTER/TYPE25 = 1 Второстепенный и главный края — это внешняя граница. края сегментов раковины. Нет краевого контакта для твердые элементы. = 11 Вторичные края — это острые края внешних сплошные сегменты и внешние края оболочки сегменты. Главными ребрами являются все ребра внешних сплошные сегменты и внешние края оболочки сегменты. = 13 Вторичные края — это острые края внешних сплошные сегменты и внешние края оболочки сегменты. Главными ребрами являются все ребра внешних сплошные сегменты и все сегменты оболочки. = 22 Вторичные и главные ребра — это все ребра из внешних сплошные сегменты и все ребра из сегментов оболочки. = 1000 По умолчанию, если /DEFAULT/INTER/TYPE25 нет определенный Нет контакта края к краю. (Целое число)
Ипстиф	Добавить флаг контактной жесткости исходя из условия устойчивости: 15 = 0 (по умолчанию) Никакой дополнительной жесткости, основанной на стабильности. состояние. = 1 (по умолчанию при использовании Istf=7) Дополнительная жесткость, основанная на условии устойчивости, равна используется с массами, поступающими из сетки (плотность, только размер элемента). = 2 Дополнительная жесткость, основанная на условии устойчивости, равна используется с исходными массами (плотностью, размером элемента, добавленная масса, массовое масштабирование). (Целое число)
grnd_IDs	Группа узлов идентификатор. 1Если определено, группа узлов будет добавлены как вторичные узлы.(Целое число)
Gap_scale	Масштабный коэффициент разрыва для всех Параметры Igap. По умолчанию = 1,0 (Реальное)
%mesh_size	Процент размера сетки (используется только при Igap = 3). По умолчанию = 0,4 (реальное)
Gap_max_s	Вторичные максимальные разрывы. 3По умолчанию = 1030 (Реал)	\([m]\)
Gap_max_m	Основные максимальные разрывы. 3По умолчанию = 1030 (Реал)	\([m]\)
Стмин	Минимальная жесткость (используется только когда Istf > 1 и Iстф < 7). 2(Реал)	\([\frac{N}{m}]\)
Стмакс	Максимальная жесткость (используется только когда Istf > 1 и Istf < 7). 2По умолчанию = 1030 (Реал)	\([\frac{N}{m}]\)
Igap0	Флаг модификации пробела для второстепенные узлы оболочки на свободных краях или элементах оболочки. 3 = 0 Установите значение, определенное в /DEFAULT/INTER/TYPE25. = 1 Установите зазор на ноль для узлов вторичной оболочки, которые по свободному краю. При контакте с кромкой оболочки свободные кромки сдвинуты так, чтобы край не выступал за пределы сегмент оболочки. = 1000 По умолчанию, если /DEFAULT/INTER/TYPE25 нет определенный Никаких изменений. (Целое число)
Ishape	Флаг, определяющий форму зазор вдоль внешней границы поверхности(й) в месте контакта узла с поверхностью. = 0 Установите значение, определенное в /DEFAULT/INTER/TYPE25. = 1 По умолчанию, если /DEFAULT/INTER/TYPE25 нет определенный Квадратный разрыв. = 2 Круглый зазор. (Целое число)
Edge_angle	Угол кромки.Только б/у с Iedge =11,13. Острые края есть включается в краевой контакт, если угол между двумя сегментами которые имеют одно и то же ребро, меньше, чем Значение Edge_angle. 6 По умолчанию = 135° (реальный)	\([deg]\)
Stfacm	Масштабный коэффициент жесткости на основе условия устойчивости (Ipstif > 0 или Istf =7). 15По умолчанию = 0,1 (Реал)
Стфак	Шкала жесткости интерфейса фактор. 2По умолчанию = 1,0 (Реал)
Фрич	Кулоновское трение (если fct_IDF= 0). По умолчанию = 0,0 (Реальное) Кулон. Масштабный коэффициент трения (если fct_IDF≠ 0). По умолчанию = 1,0 (реальное)
Тпрессфит	Срок подачи заявки контактные силы, возникающие из-за начального проникновения (запрессовки). Только б/у только с явным решением Inacti= -1. 9По умолчанию = время, соответствующее 10000 циклов (реальный)	\([s]\)
Тстарт	Время начала. 10 (Реал)	\([s]\)
Тстоп	Временная деактивация время. 10По умолчанию = 1030 (Реал)	\([s]\)
IBC	Флаг деактивации граничные условия при ударе. (логическое значение)
Бездействие	Флаг начального проникновения. = -1 Учитываются все первоначальные проникновения. = 0 Установите значение, определенное в /DEFAULT/INTER/TYPE25. = 5 Основной сегмент смещается при первоначальном проникновении ценность \(P_{0}\) . Если \(P\geP_{0}\) , тогда \(P'=P−P_{0}\) , где \(P_{0}\) является начальным проникновение. = 1000 По умолчанию, если /DEFAULT/INTER/TYPE25 нет определенный Во внимание будут приниматься только незначительные первоначальные проникновения. счет. (Целое число)
ВИС	Критическое демпфирование коэффициент жесткости интерфейса. По умолчанию = 0,05. (Реал)
Ифрик	Формулировка трения флаг. Используется только если fric_ID не определен. = 0 (по умолчанию) Статический закон трения Кулона = 1 Обобщенный закон вязкого трения = 2 (Модифицированный) Дармстадский закон трения = 3 Закон трения Ренара = 4 Экспоненциальный закон трения затухания. Для контакта от края до края только статический кулон доступен закон трения. (Целое число)
Ифильтр	Флаг фильтрации трения. 8 = 0 (по умолчанию) Фильтр не используется. = 1 Простой числовой фильтр. = 2 Стандартный фильтр -3 дБ с периодом фильтрации. = 3 Стандартный фильтр -3 дБ с частотой среза. (Целое число)
частота	Фильтрация коэффициент.По умолчанию = 1,0 (Реальный)
sens_ID	Идентификатор датчика для активировать/деактивировать интерфейс. (Целое число)
Дстиф	Шаг времени, используемый для контакта жесткость, основанная на условии устойчивости. (Ipstif > 0 или Istf = 7). 15По умолчанию = 0 (Реал)	\([s]\)
fric_ID	Идентификатор трения для определение трения для выбранных пар деталей. = 0 (по умолчанию) Используйте параметры трения, определенные в этом интерфейсе. ≠ 0 Используйте /FRICTION/fric_ID При контакте от края к краю действует только изотропное трение. считается. Если соответствующая модель ортотропна, то только учитывается коэффициент в направлении 1 контактирующей части учитывать контакт от края до края. (Целое число)
C1 - C6	Коэффициент закона трения. 5(Реал)	См. Таблицу 1.
IVIS2	Флаг адгезии интерфейса. 12 = -1 Включить поперечное сцепление и тангенциальную вязкость сила. = 0 (по умолчанию) Отсутствие сил адгезии на границе раздела. (Целое число)
ViscFluid	Вязкость жидкости при интерфейс. 12(Реал)	\([Pa⋅s]\)
СигМаксАдх	Максимальное поперечное адгезионное напряжение на границе раздела. 12(Реал)	\([Pa]\)
ВискАдфакт	Тангенциальная вязкая Масштабный коэффициент силы сопротивления. 12(Реал)
Кте	Кондуктивный теплообмен коэффициент (если fct_IDK = 0). 14По умолчанию = 0,0 (Реал)	\([\frac{W}{m^{2}K}]\)
Масштабный коэффициент теплообмена (если fct_IDK ≠ 0)По умолчанию = 0,0 (реальное)
fct_IDK	Идентификатор функции для определение теплообмена с контактным давлением. По умолчанию = 0 (Целое число)
АскалеК	Масштабный коэффициент по оси абсцисс fct_IDK.По умолчанию = 1,0 (реальное)	\([Pa]\)
оттенок	Температура интерфейса. 13(Реал)	\([K]\)
Ithe_form	Рецептура теплового контакта флаг. 13 = 0 Обмен только между интерфейсами (постоянная температура) и оболочки (вторичная сторона). = 1 Теплообмен между всеми контактирующими деталями. (Целое число)
Фрад	Радиационный фактор. 13(Реал)	\([\frac{W}{m^{2}K^{4}}]\)
Драд	Максимальное расстояние для расчет излучения.(Реальный)	\([m]\)
fct_IDc	Идентификатор функции для определение коэффициента кондуктивного теплообмена как функции расстояние. 14По умолчанию = 0 (Целое число)
Дконд	Максимальное расстояние для кондуктивный теплообмен. По умолчанию = 0,0 (Реал)	\([m]\)
Подвиги	Коэффициент нагрева трением вторичного. 14(Реал)
Фитм	Коэффициент нагрева трением основного. 14(Реал)
fct_IDF	Коэффициент трения при идентификатор температурной функции. По умолчанию = 0 (Целое число)
Thickness_s	Толщина контакта вторичной обмотки. Б/у только для Igap=5. (Настоящий)	\([m]\)
Scale_s	Разрыв шкалы средней. Используется только для Игап =5. По умолчанию = 1,0 (реальное)
Thickness_m	Толщина контакта основной. Только б/у для Igap =5. (Настоящий)	\([m]\)
Scale_m	Разрыв шкалы основного. Используется только для Igap =5. По умолчанию = 1,0 (реальное)

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

(1)-1	(1)-2	(1)-3	(1)-4	(1)-5	(1)-6	(1)-7	(1)-8	(1)-9	(1)-10
							IBCX	IBCY	IBCZ

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

Поле	Содержание	Пример единицы СИ
IBCX	Флаг деактивации X граничное условие при ударе. =0 Бесплатная глубина резкости =1 Фиксированная глубина резкости (логическое значение)
IBCY	Флаг деактивации Y граничное условие при ударе. =0 Бесплатная глубина резкости =1 Фиксированная глубина резкости (логическое значение)
IBCZ	Флаг деактивации Z граничное условие при ударе. =0 Бесплатная глубина резкости =1 Фиксированная глубина резкости (логическое значение)

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

Контактные пары главная/вторичная

можно определить тремя способами:

Только одна самоударная поверхность: surf_ID1 > 0 и surf_ID2 = 0.

Симметрично между поверхностями: surf_ID1 > 0 и surf_ID2 > 0.

Узлы на поверхности: grnd_IDs > 0,

surf_ID1 = 0 и surf_ID2 > 0
grnd_IDs > 0 используется

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

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

сегменты (сдвинутые на половину толщины (t)) с противоположными направлениями нормалей

будет сгенерировано:

В случае SPMD каждый основной сегмент определяется

surf_IDi (i=1, 2) должно быть связанный с элементом (возможно, с пустым элементом).

В случаях, когда

используются квадратичные элементы, поверхности рекомендуется определять по используя /SURF/PART/EXT, поскольку в этом случае средние узлы При контактной обработке используются квадратичные элементы.

Поверхность

определение /SURF/PART/ALL недоступно с ТИП25.

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

\(K\) рассчитывается как: \(K=\frac{max[St_{min},min(St_{max},K_{n})]}{2}\) Где, \(K_{n}\) зависит от I стф : - I

стф = 1000 , \(K_{n}=\frac{K_{m}}{2}\)

I стф = 2 , \(K_{n}=\frac{K_{m}+K_{s}}{2}\)

I стф = 3 , \(K_{n}=max⁡(K_{m},K_{s})\)

I стф = 4 , \(K_{n}=min⁡(K_{m},K_{s})\)

I стф = 5 , \(K_{n}=\frac{K_{m}⋅K_{s}}{K_{m}+K_{s}}\)

I стф = 7 , \(K_{n}=K_{msdt}\) . Видеть 15 .

\(K_{m}\) : жесткость основного сегмента и рассчитывается как: - \(K_{m}=Stfac⋅E⋅t\)

, когда основной отрезок лежит на

оболочка.

\(K_{m}=Stfac⋅B⋅\frac{S^{2}}{V}\) , когда основной сегмент лежит на

твердый.

\(K_{m}=max⁡(Stfac⋅E⋅t,Stfac⋅B⋅\frac{S^{2}}{V})\) , когда основной сегмент используется совместно

оболочка и твердая.

\(K_{s}\) : Жесткость вторичного узла эквивалентна

узловая жесткость рассматривается для интерфейса ТИПА 25 и рассчитывается как:

\(K_{s}=Stfac⋅E⋅t\) , когда узел подключен к оболочке

элемент,

\(K_{s}=Stfac⋅B⋅\sqrt[3]{V}\) , когда узел подключен к твердому телу

элемент.

Где, \(S\) Область сегмента \(V\) Объем твердого тела \(B\) Объемный модуль \(t\) Толщина корпуса Значение Stfac может быть больше, чем

1.0. Коэффициент жесткости не имеет ограничений (значение, большее чем 1,0, может уменьшить начальный шаг по времени).

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

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

Зазор рассчитывается

автоматически для каждого удара как:

If I зазор = 1 , переменный разрыв рассчитывается как: \(min(g_{s},Gap_max_s)+min(g_{m},Gap_max_m)\)

If I зазор = 2 , переменный разрыв рассчитывается как: \(min(g_{s},Gap_max_s)+min(g_{m},Gap_max_m)\) с отключением вторичных узлов при

размер элемента меньше значений зазора:

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

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

If I зазор = 3 , переменный разрыв рассчитывается как: \(min[min(g_{s},Gap_max_s)+min(g_{m},Gap_max_m),%mesh_size⋅(g_{s_l}+g_{m_l})]\) Где, - \(g_{m}\)

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

, с

\(t\) это толщина основной элемент для элементов оболочки

\(g_{m}=0\)

, для кирпича

элементы

\(g_{s}\) : разрыв вторичного узла: \(g_{s}=0\)

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

не подключен ни к какому элементу или подключен только к кирпичные или пружинные элементы.

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

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

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

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

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

соединяется с элементами фермы или балки, при этом \(S\) являющийся поперечным сечением элемента 1D.

Если флаг модификации разрыва

для узлов вторичной оболочки на свободных краях Igap0 устанавливается равным 1: \(g_{s}\) обнуляется, если вторичный узел лежит на свободных краях вторичного поверхность. Флаг модификации пробела для вторичной оболочки узлы на свободных краях не влияет, если вторичный узел определяется через необязательную группу узлов (grnod_IDs).

Если

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

\(gm_l\) : длина наименьшего края

основного сегмента.

\(gs_l\) : если вторичный узел

принадлежит основной поверхности,

\(gs_l\) это длина

наименьший край основных сегментов, соединенный со второстепенными

узел,

\(gs_l\) =1E+30, иначе. В любом

случай, \(g_{m}\) и \(g_{s}\) ограничены отдельно на Gap_max_m и Gap_max_s до того, как будет установлен пробел. вычислено.

Если вторичный узел не принадлежит

основная поверхность, зазор остается

\(min(g_{s},Gap_max_s)+min(g_{m},Gap_max_m)\)

If I зазор =5, переменный разрыв рассчитывается как: \(min(g_{s},Gap_max_s)+ming_{m},Gap_max_m\) Где, - \(g_{m}\)

: основной элемент

разрыв:

\(g_{m}=Scale_m·\frac{Thickness_m}{2}\) с \(t\) являющаяся толщиной

основной элемент для элементов оболочки.

\(g_{m}=0\)

для кирпича элементы.

\(g_{s}\) : вторичный узел

разрыв:

\(g_{m}=0\)

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

\(g_{s}=Scale_s·\frac{Thickness_s}{2}\) , если вторичный узел

соединенный с элементом оболочки, пружиной, фермой или балкой

элементы.

Для контакта узла с поверхностью,

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

внешняя граница.

I форма определяет, является ли форма этого зазора квадратной или круглой, и контактная сила

(нормальное) направление.

I форма имеет

не влияет на зазор и его форму при контакте кромок.

В зависимости от I форма промежуток, используемый для

контакт на внешней границе основной поверхности и результирующее направление силы.

Разрыв))*

Разрыв))*
Ishape

=1 недоступен при Igap =3 и затем будет сброшен на Ishape =2.

Для кромки элемента оболочки до

краевой контакт, зазор круглый. Основной боковой контактный зазор на свободном крае составляет

сдвинута так, чтобы край не выходил за пределы сегмента оболочки.

Поведение свободного края контактного зазора на вторичной стороне зависит от

ценность

I разрыв0 , как показано ниже. .. image:: images/inter_type25_starter_r_inter_type25_edge_contact_secondary.png

(Рисунок 6. Вторичная сторона краевого контакта)

Для твердых элементов, когда

I край = 11 и I край = 13 , вторичная сторона состоит только из острых кромок

с углом меньшим, чем

Edge_angle . Для I край = 22 , все ребра твердотельных элементов рассматриваются на

вторичная сторона. На главной стороне включены все ребра из твердотельных элементов.

все 3

I край случаи. .. image:: images/inter_type25_starter_r_inter_type25_secondary_side_effect.png

(Рисунок 7. Второстепенные боковые кромки для Iedge=11 и Iedge=13.)

If

fric_ID определено, контактное трение определяется в /FRICTION и входы трения ( I Фриц , C 1 и т. д.) в этом входе карты не используются. Силы трения

являются:
\(F_{t}^{new}=min⁡(\muF_{n},F_{adh})\) Сила сцепления рассчитывается как: \(F_{adh}=F_{t}^{old}+\Delta F_{t}\)

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

Где, \(\mu\) представляет собой коэффициент кулоновского трения и определяется как: - Для флага

I Фриц по умолчанию: \(\mu=Fric\)

с \(F_{T}\le\mu⋅F_{N}\) (Кулоновское трение)

Для флага I Фриц > 1, новое трение

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

по функции:

\(\mu=\mu(p,V)\) Где, \(p\) Давление нормальной силы на главный

сегмент

\(V\) Тангенциальная скорость вторичного узла относительно

в основной сегмент

В настоящее время коэффициенты C1–C6 используются для определения

переменный коэффициент трения \(\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}})\)

если \(V\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}}))\)

если \(V\in[C_{5},C_{6}]\)

\(\mu=C_{2}−\frac{1}{\frac{1}{C_{2}−C_{4}}+(V−C_{6})^{2}}\)

если \(V\geC_{6}\)

Где, - \(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}\ne0\) , первый критический

скорость, должна быть > 0

\(C_{6}=V_{cr2}\) , второй критический

скорость, должна быть

\(>V_{cr1}\)

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

вторая критическая скорость

\(V_{cr2}=C_{6}(C_{5}<C_{6})\) .

Коэффициент статического трения \(C_{1}\) и динамическое трение

коэффициент

\(C_{2}\) , должно быть меньше, чем

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

\(C_{3}\) ( \(C_{1}\leC_{3}\) и \(C_{2}\leC_{3}\) ).

Минимальный коэффициент трения \(C_{4}\) должно быть меньше, чем

статический коэффициент трения

\(C_{1}\) и динамическое трение

коэффициент

\(C_{2}\) ( \(C_{4}\leC_{1}\) и \(C_{4}\leC_{2}\) ).

I Фриц = 4 (Экспоненциальный закон трения затухания) Предполагается, что коэффициент трения зависит от

относительная скорость \(V\) поверхностей, соприкасающихся согласно:

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

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

Если Ифильтр = 1,

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

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

фильтр = 1 ➤ \(\alpha=X_{freq}\) , простой числовой фильтр

со значением между

0 и 1 .

I фильтр = 2 ➤ \(\alpha=\frac{2⋅\pi}{X_{freq}}\) , стандартный фильтр -3 дБ,

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

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

I фильтр = 3 ➤ \(\alpha=2⋅\pi⋅X_{freq}⋅dt\) стандартный фильтр -3дБ,

с

X частота =

частота резания

Бездействие

и Ipen_max ,

начальная обработка проникновения:

Бездействие = 1000 :

Начальные проникновения игнорируются: контактные силы не применяются, но

узлы не отключаются от контакта; если узел выходит из строя

контакт, а затем снова вступает в контакт, тогда контактные силы

применяется.

Неактивность = -1:

Начальные силы приложены ко всем проникающим узлам. Высокий начальный следует избегать проникновений, так как они могут привести к сильному контакту сил и привести к высокой энергетической ошибке в начале расчет. Контактные силы, вызванные первоначальным проникновением, равны увеличивается с нуля во время активации контакта, определенное Tstart или sens_ID, до максимального значения с течением времени Tpressfit после активация контакта. Увеличение контактных сил, вызванное начальным проникновение позволяет моделировать ситуации прессовой посадки в моделях. Чтобы избежать динамических эффектов, Tpressfit не должен быть слишком маленький. По умолчанию Tpressfit — это время, соответствующее до 10000 циклов. Жесткость определяется автоматически при Inacti= -1, Istf, Stmin и Stmax не оказывают никакого влияния, для масштабирования жесткости можно использовать только Stfac.

Бездействие = 5 :

Основной сегмент смещается на начальную величину проникновения (

\(P_{0}\) ); следовательно, в нулевой момент времени нет начального

применяются силы.

Положение основного сегмента восстанавливается только в случае отскока большего

чем \(P_{0}\)

.

В противоположном случае, когда вторичное

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

\(P'=P−P_{0}\) .. image:: images/inter_type25_starter_r_inter_type24_inacti=5.png

alt

inter_type24_inacti=5

(Рисунок 9.)

Пересечения и большое начальное проникновение (Inacti= -1 и

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

Когда sens_ID определен для

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

Для выходных сил: Когда

Тип контакта асимметричный, выходной нормальный контакт. силы в истории времени рассчитываются правильно, если две поверхности хорошо разделены.

ИВИС2

= -1 : используется для добавления

адгезия в нормальном направлении и силы вязкого сопротивления в тангенциальном направлении.

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

Когда

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

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

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

Gap_scale . .. image:: images/inter_type25_starter_r_inter_type25_ivis2.png

(Рисунок 10.)

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

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

Сила адгезии препятствует выходу узла из зоны адгезии.

зоне и применяется в нормальном направлении.

\(F_{N}=\frac{SigMaxAdh⋅Area}{\frac{1}{2}Gap}(\frac{1}{2}Gap−P_{adh})\) Где, \(Area\) Площадь вторичной поверхности \(P_{adh}\) Проникновение в зону адгезии \(Gap\) Контактный зазор, рассчитанный в комментариях. 3 Клейкая пружина разрывается при выходе узла из зоны адгезии.

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

вторичные узлы входят в зону адгезии. Вязкая тангенциальная

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

как:

\(F_{T}=−(ViscAdhFact)\frac{ViscFluid⋅Area}{\frac{1}{2}Gap}V_{rel}\) Где, \(Area\) Площадь вторичной поверхности \(V_{rel}\) Проникновение в зону адгезии \(Gap\) Контактный зазор .. image:: images/inter_type25_starter_r_inter_type25_contact_gap.png

(Рисунок 11.)

Теплообмен

By I тот = 1 (активируется передача тепла) для учета тепла

обмен и тепловое трение при контакте.

Если Ithe_form=0, то теплообмен происходит между вторичным

узлов и постоянной температуры контакта Тонировка.

Если Ithe_form=1, то теплообмен происходит между всеми контактами.

штук.
В этом случае Tint используется только тогда, когда Ithe_form=0. Предполагается, что температура основной стороны равна

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

Контактный теплообмен может включать: - Теплопроводность

Ithe = 1 нужен материал вторичной стороны

быть термическим материалом с использованием формулы конечных элементов для теплопередача (/HEAT/MAT).

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

попадает в пропасть:
\(gap=max[Gap_{min},min(Fscale_{gap}⋅g_{s},Gap_{max})]\) Коэффициент теплообмена - Если fct_IDK = 0, то Kthe – коэффициент теплообмена и теплообмен

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

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

контактное давление:

\(K=K_{the}⋅f_{K}(Ascale_{K},P)\) В то время как,

\(f_{K}\) является функцией fct_IDK.

Теплопроводность и излучение Когда

fct_IDc ≠ 0, коэффициент теплопередачи может меняться в зависимости от расстояние \(d\) когда \(Gap<d\leD_{cond}\)

. В этой зоне проводящие и

рассмотрены радиационные потоки теплопередачи.: Абсцизы

и ординаты этой функции fct_IDc должно быть между 0 и 1.

Коэффициент теплопередачи

вычисляется как:: \(K=K_{the}(P=0)⋅f_{c}(\frac{d−Gap}{Dcond−Gap})\) Максимальное значение

\(K\) равно значению Kthe, когда Kthe постоянно. В противном случае, если Kthe зависит от давления, максимум равен значению Kthe для контактного давления. \(P=0\)

.

\(K\) падает до нуля, когда расстояние равен Dконд.

Тепловое излучение Излучение считается контактным, если

\(F_{rad}\ne0\) и расстояние, \(d\)

, вторичного узла к

основной сегмент:

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

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

D рад рассчитывается как максимум: - верхнее значение разрыва (в момент времени 0) среди всех узлов - наименьшая длина стороны вторичного элемента .. note:

Рекомендуется не устанавливать слишком высокое значение для
 Drad, что может снизить производительность 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}\) Излучательная способность основной поверхности.

Тепловое трение

Фрикционный

энергия может быть преобразована в тепло, если

I тот > 0 для интерфейса. - Подвиги определены

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

Fheatm определен

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

Тепло трения

\(Q_{Fric}\) определяется для формулировки жесткости:

\(Q_{Fric}=Fheat⋅\frac{(F_{adh}−F_{t})}{K}⋅F_{t}\)

Контактная жесткость на основе

состояние устойчивости.

Этот метод интересно использовать, особенно в случаях

материалы имеют большую разницу в жесткости. Жесткость условия устойчивости вычисляется по массе и шагу по времени:

\(K_{msdt}=Stfacm⋅\frac{m_{m}⋅m_{s}}{m_{m}+m_{s}}⋅\frac{1}{\Delta t_{c}}^{2}\) Где, \(m_{m}\) и \(m_{s}\) Главные и вторичные узловые массы. \(\Delta t_{c}\) Временной шаг начального решения или введенное пользователем значение Дстиф . Если Istf = 7, контактная жесткость рассчитывается только на основе

условие устойчивости.
\(K_{n}=K_{msdt}\) В противном случае, если Ipstif > 0, условие устойчивости жесткости

\(K_{msdt}\) используется, если оно больше \(K_{n}\)

, жесткость определяется флагом Istf (см. комментарий 2).: \(K=\frac{maxSt_{min},minSt_{max},max(K_{n},K_{msdt})}{2}\) По умолчанию, в случае контактной жесткости, основанной на

условие устойчивости (если Istf = 7 или Ipstif > 0), узловые массы, используемые для \(K_{msdt}\) первоначально рассчитываются с использованием материала плотность и размер элементов. Он не учитывает добавленную массу и может привести низкой жесткости и контактной неустойчивости. В этом случае Ipstif= 2 можно использовать для рассмотрения изначально добавлен масс эффект.