============= /INTER/TYPE24 ============= Ключевое слово формата блока TYPE24 — это общий контакт между узлами и поверхностью. интерфейс с использованием метода штрафа. Штрафная жесткость постоянна и, следовательно, шаг по времени не изменяется (для стандартная штрафная жесткость). Твердым элементам задается нулевой зазор. Три типа входные контакты могут быть определены: одна поверхность, поверхность к поверхности или узлы к поверхность. Этот контактный интерфейс может заменить интерфейс TYPE3, TYPE5 или TYPE7. Для неявное решение, этот интерфейс TYPE24 доступен только с SMP. Формат ------ .. csv-table:: :header: "(1)", "(2)", "(3)", "(4)", "(5)", "(6)", "(7)", "(8)", "(9)", "(10)" :widths: 10, 10, 10, 10, 10, 10, 10, 10, 10, 10 "/INTER/TYPE24/inter_ID/unit_ID", "/INTER/TYPE24/inter_ID/unit_ID", "/INTER/TYPE24/inter_ID/unit_ID", "/INTER/TYPE24/inter_ID/unit_ID", "/INTER/TYPE24/inter_ID/unit_ID", "/INTER/TYPE24/inter_ID/unit_ID", "/INTER/TYPE24/inter_ID/unit_ID", "/INTER/TYPE24/inter_ID/unit_ID", "/INTER/TYPE24/inter_ID/unit_ID", "/INTER/TYPE24/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", "", "", "Иедж", "Edge_angle", "Edge_angle", "Gap_max_s", "Gap_max_s", "Gap_max_m", "Gap_max_m" "Стмин", "Стмин", "Стмакс", "Стмакс", "Igap0", "Ипен0", "Ipen_max", "Ipen_max", "Stfacm", "Stfacm" Обязательные поля .. 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", "", "", "Бездействие", "ВИС", "ВИС", "", "", "Тпрессфит", "Тпрессфит" "Ифрик", "Ифильтр", "частота", "частота", "", "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", "", "", "", "", "", "", "", "" Определение ----------- .. csv-table:: :header: "Поле", "Содержание", "Пример единицы СИ" :widths: 33, 33, 33 "inter_ID", "Интерфейс идентификатор.(Целое число, максимум 10 цифр)", "" "unit_ID", "Идентификатор объекта.(Целое число, максимум 10 цифр)", "" "inter_title", "Интерфейс заголовок.(Символ, максимум 100 символов)", "" "surf_ID1", "Первый идентификатор поверхности. 1 (Целое число)", "" "surf_ID2", "Вторая поверхность идентификатор.(Целое число)", "" "истф", "Жесткость интерфейса флаг определения (используется только тогда, когда Istf > 1 и Istf < 7). 2 3 = 0 Установите значение, определенное в /DEFAULT/INTER/TYPE24. = 2 Жесткость интерфейса представляет собой среднее значение основных и вторичная жесткость. = 3 Жесткость интерфейса – максимальная из основных и вторичная жесткость. = 4 Жесткость интерфейса – минимальная из основных и вторичная жесткость. = 5 Жесткость интерфейса — это основная и вторичная жесткость. в серии. = 6 Жесткость интерфейса – минимальная из основных и вторичная жесткость со специальной регулировкой для улучшения сходимость для неявных решений. =7 Жесткость интерфейса основана только на стабильности состояние. 13 = 12 Метод Ниче используется со средним значением основных и вторичная жесткость. = 13 Метод Ниче применяется с максимумом основных и вторичная жесткость. = 14 Метод Ниче используется с минимумом основных и вторичная жесткость. = 1000 По умолчанию, если /DEFAULT/INTER/TYPE24 нет определен. Жесткость интерфейса рассчитывается только по основной стороне. жесткость. (Целое число)", "" "Irem_i2", "Деактивация флага для вторичный узел, если была использована та же пара контактов (узел/сегмент). определено в интерфейсе TYPE2. =0 Установите значение, определенное в /DEFAULT/INTER/TYPE24. =1 По умолчанию, если /DEFAULT/INTER/TYPE24 нет определенный Вторичные узлы в /INTER/TYPE2 привязаны контакты удаляются из этого контакта. =3 Никаких изменений во вторичных узлах.", "" "Идель", "Удаление узла и сегмента флаг. = 0 Установите значение, определенное в /DEFAULT/INTER/TYPE24. = 1 Когда все элементы (4-узловые оболочки, 3-узловые оболочки, твердые вещества), связанные с одним сегментом, удаляются, сегмент удален с основной стороны интерфейса. Он также удаляется в случае явного удаления с помощью Ключевое слово Radioss Engine /DEL в файле Engine. Кроме того, из списка удаляются несвязанные узлы. вторичная сторона интерфейса. = 2 Когда 4-узловая оболочка, 3-узловая оболочка или сплошной элемент удаляется, соответствующий сегмент удаляется из основная часть интерфейса. Он также удаляется в случай явного удаления с использованием ключевого слова Radioss Engine /DEL в файле Engine. Кроме того, из списка удаляются несвязанные узлы. вторичная сторона интерфейса. = 1000 По умолчанию, если /DEFAULT/INTER/TYPE24 нет определен. Никакого удаления. Примечание. Использование Idel приводит к более высокая стоимость процессора.", "" "Ипстиф", "Добавить флаг контактной жесткости исходя из условия устойчивости: 13 = 0 (по умолчанию) Никакой дополнительной жесткости, основанной на стабильности. состояние. = 1 (по умолчанию при использовании Истф=7) Дополнительная жесткость, основанная на условии устойчивости, равна используется с массами, поступающими из сетки (плотность, только размер элемента). = 2 Дополнительная жесткость, основанная на условии устойчивости, равна используется с исходными массами (плотностью, размером элемента, добавленная масса, массовое масштабирование). (Целое число)", "" "grnd_IDs", "Группа узлов идентификатор. 1 Если определено, группа узлов будет добавлена как вторичная. узлы.(Целое число)", "" "Иедж", "Флаг контакта от края до края. = 0 (по умолчанию) Установите значение, определенное в /DEFAULT/INTER/TYPE24. = 1 Контакт от края до края активируется с помощью внешнего края границ из surf_ID1 и surf_ID2 и острые края между контактными сегментами. 9 = 1000 По умолчанию, если /DEFAULT/INTER/TYPE24 нет определенный Нет контакта края к краю. (Целое число)", "" "Edge_angle", "Угол кромки. Используется только если Иэдж =1. Острые края включены в контакт с кромкой если угол между двумя сегментами, имеющими одно и то же ребро, равен меньше значения Edge_angle. 9По умолчанию = 135° (реальный)", ":math:`[deg]`" "Gap_max_s", "Вторичные максимальные разрывы. 4По умолчанию = 1030 (Реал)", ":math:`[m]`" "Gap_max_m", "Основные максимальные разрывы. 4 По умолчанию = 1030 (реальное)", ":math:`[m]`" "Стмин", "Минимальная жесткость (используется только когда Istf > 1 и Iстф < 7). 2 (Реал)", ":math:`[\frac{N}{m}]`" "Стмакс", "Максимальная жесткость (используется только когда Istf > 1 и Istf < 7). 2 По умолчанию = 1030 (реальное)", ":math:`[\frac{N}{m}]`" "Igap0", "Флаг модификации пробела для вторичные узлы оболочки на свободных краях. = 0 (по умолчанию) Установите значение, определенное в /DEFAULT/INTER/TYPE24. = 1 Установите зазор на ноль для узлов вторичной оболочки. = 1000 По умолчанию, если /DEFAULT/INTER/TYPE24 нет определенный Никаких изменений. (Целое число)", "" "Ипен0", "Первоначальное проникновение флаг обнаружения. 8 = 0 (по умолчанию) Установите значение, определенное в /DEFAULT/INTER/TYPE24. = 1 Включая самовоздействие в каждой части. = 1000 По умолчанию, если /DEFAULT/INTER/TYPE24 нет определенный Исключая самовоздействие в каждой части. (Целое число)", "" "Ipen_max", "Максимальный начальный проникновение. 7 По умолчанию = 0 (реальный)", ":math:`[m]`" "Stfacm", "Масштабный коэффициент жесткости на основе условия устойчивости (Ipstif > 0 или Istf =7). 13 = 0,1 (по умолчанию) Стандартная жесткость интерфейса только с использованием Istf. (Настоящий)", "" "Стфак", "Шкала жесткости интерфейса фактор. 2 По умолчанию = 1,0 (реальное)", "" "Фрич", "Кулоновское трение. 5 (Реал)", "" "Тстарт", "Время начала. 9 (Реал)", ":math:`[s]`" "Тстоп", "Временная деактивация время. 9По умолчанию = 1030 (Реал)", ":math:`[s]`" "IBC", "Флаг деактивации граничные условия при ударе. (логическое значение)", "" "Бездействие", "Флаг начального проникновения. 7 = -2 Учитываются все первоначальные проникновения и жесткость для удаления пересечений основана на Istf = 2 и Стфак. = -1 Учитываются все первоначальные проникновения и жесткость для удаления пересечений основана на Istf = 4 и Стфак. = 0 (по умолчанию) Установите значение, определенное в /DEFAULT/INTER/TYPE24. = 5 Основной сегмент смещается при первоначальном проникновении ценность :math:`P_{0}` . Если :math:`P\geP_{0}` , тогда :math:`P'=P−P_{0}` , где :math:`P_{0}` является начальным проникновение. = 1000 По умолчанию, если /DEFAULT/INTER/TYPE24 нет определенный Во внимание будут приниматься только незначительные первоначальные проникновения. счет. (Целое число)", "" "ВИС", "Критическое демпфирование коэффициент жесткости интерфейса. По умолчанию = 0,05. (Реал)", "" "Тпрессфит", "Срок подачи заявки контактные силы из-за начального проникновения (запрессовка). Только б/у только с явным решением Inacti=-1. 7По умолчанию = время, соответствующее 10000 циклов (реальный)", ":math:`[s]`" "Ифрик", "Флаг формулировки трения. 5Используется только если fric_ID не определен. = 0 (по умолчанию) Статический закон трения Кулона. = 1 Обобщенный закон вязкого трения. = 2 (Модифицированный) Дармстадский закон трения. = 3 Закон трения Ренара. = 4 Экспоненциальный закон трения затухания. (Целое число)", "" "Ифильтр", "Флаг фильтрации трения. 6 = 0 (по умолчанию) Фильтр не используется. = 1 Простой числовой фильтр. = 2 Стандартный фильтр -3 дБ с периодом фильтрации. = 3 Стандартный фильтр -3 дБ с частотой среза. (Целое число)", "" "частота", "Коэффициент фильтрации. 6 По умолчанию = 1,0 (реальное)", "" "Дстиф", "Шаг времени, используемый для контакта жесткость, основанная на условии устойчивости. (Ipstif > 0 или Istf = 7). 13По умолчанию = 0 (Реал)", ":math:`[s]`" "sens_ID", "Идентификатор датчика для активировать/деактивировать интерфейс. 9 (целое число)", "" "fric_ID", "Идентификатор трения для определение трения для выбранных пар деталей. = 0 (по умолчанию) Используйте параметры трения, определенные в этом интерфейсе. ≠ 0 Используйте /FRICTION/fric_ID. (Целое число)", "" "C1 - C6", "Коэффициент закона трения. 5(Реал)", "См. Таблицу 1." Флаги для деактивации граничных условий: IBC -------------------------------------------- .. csv-table:: :header: "(1)-1", "(1)-2", "(1)-3", "(1)-4", "(1)-5", "(1)-6", "(1)-7", "(1)-8" :widths: 12, 12, 12, 12, 12, 12, 12, 12 "", "", "", "", "", "IBCX", "IBCY", "IBCZ" Определение ----------- .. csv-table:: :header: "Поле", "Содержание", "Пример единицы СИ" :widths: 33, 33, 33 "IBCX", "Флаг деактивации X граничное условие при ударе. =0 Бесплатная глубина резкости =1 Фиксированная глубина резкости (логическое значение)", "" "IBCY", "Флаг деактивации Y граничное условие при ударе. =0 Бесплатная глубина резкости =1 Фиксированная глубина резкости (логическое значение)", "" "IBCZ", "Флаг деактивации Z граничное условие при ударе. =0 Бесплатная глубина резкости =1 Фиксированная глубина резкости (логическое значение)", "" Комментарии ----------- 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)) с противоположными направлениями нормалей будет сгенерировано: .. image:: images/inter_type24_starter_r_inter_type24.png :alt: inter_type24 *(Рисунок 1.)* В случае SPMD каждый основной сегмент определяется surf_IDi (i=1, 2) должно быть связанный с элементом (возможно, с пустым элементом). Поверхность определение /SURF/PART/ALL недоступно с ТИП24. 2. Контактная жесткость, K определяется как: :math:`K=\frac{max[St_{min},min(St_{max},K_{n})]}{2}` Где, :math:`K_{n}` зависит от I стф - If I стф = 1000 , :math:`K_{n}=\frac{K_{m}}{2}` (Жесткость по умолчанию) - I стф = 2 , :math:`K_{n}=\frac{K_{m}+K_{s}}{2}` - I стф = 3 , :math:`K_{n}=max⁡(K_{m},K_{s})` - I стф = 4 , :math:`K_{n}=min⁡(K_{m},K_{s})` - I стф = 5 , :math:`K_{n}=\frac{K_{m}⋅K_{s}}{K_{m}+K_{s}}` - I стф = 6 , :math:`K_{n}=min⁡(K_{m},K_{s})` , Мягкая жесткость. Этот вариант только доступен с неявным решением. Для каждого контакта сделать нелинейная итерационная сходимость проще, меньшая начальная жесткость использованный; по функции реакции контактирующих частей (при увеличение проникновения или отскока), жесткость будет отрегулирована, но всегда меньше входной жесткости. - I стф = 7 , :math:`K_{n}=K_{msdt}` . Видеть Комментарий 13 . - I стф = 12 , Ниче метод с :math:`K_{n}=\frac{K_{m}+K_{s}}{2}` - I стф = 13 , Ниче метод с :math:`K_{n}=max(K_{m},K_{s})` - I стф = 14 , Ниче метод с :math:`K_{n}=min(K_{m},K_{s})` :math:`K_{m}` : жесткость основного сегмента - :math:`K_{m}=Stfac⋅E⋅t` , когда основной отрезок лежит на оболочка - :math:`K_{m}=max⁡(Stfac⋅E⋅t,Stfac⋅B⋅\frac{S^{2}}{V})` , когда основной сегмент используется оболочкой совместно и твердый - :math:`K_{m}=Stfac⋅B⋅\frac{S^{2}}{V}` , когда основной сегмент лежит на твердом теле. :math:`K_{s}` : Жесткость вторичного узла эквивалентна узловая жесткость рассматривается для интерфейса TYPE24 и рассчитывается как: - :math:`K_{s}=Stfac⋅E⋅t` , когда узел подключен к оболочке элемент, - :math:`K_{s}=Stfac⋅B⋅\sqrt[3]{V}` , когда узел подключен к твердому телу элемент. Где, :math:`S` Область сегмента :math:`V` Объем твердого тела :math:`B` Объемный модуль :math:`t` Толщина корпуса Значение Stfac может быть больше, чем 1.0. Коэффициент жесткости не имеет ограничений (значение, большее чем 1,0, может уменьшить начальный шаг по времени). При использовании /PROP/VOID и /MAT/VOID, материал необходимо ввести свойства и толщину материала VOID; в противном случае, контактная жесткость пустотных элементов будет равна нулю. Это особенно важно, если элементы оболочки VOID имеют общие элементы с твердотельными элементами в качестве жесткость элементов оболочки используется в контактном расчете. Если Необходимо учитывать неявный анализ (определяется /IMPLICIT), значение по умолчанию. istf = 4; Iстф = 6 есть рекомендуется для гибких или изгибающихся конструкций. 3. Метод Ниче 1 представляет собой контактный алгоритм, используемый для решения контактные уравнения. Этот метод можно использовать вместо традиционного наказания. контактный метод. Ниче-метод лучше предотвращает контактные проникновения особенно в случаях контакта компонентов из разных материалов жесткость, такая как контакт между резиной и металлом. Для этих случаев это сложно выбрать штрафную жесткость, поскольку использование минимальной жесткости может приводят к большим проникновениям, а использование максимальной контактной жесткости может привести к сильно деформированные элементы и уменьшающийся шаг по времени. Метод Ниче может улучшить моделирование таких проблем, поскольку контактные силы вычисляются с использованием контактная жесткость, проникновения и напряжения элемента. Метод Ниче доступен только для контакта одной поверхности или поверхности к поверхности между твердые и не доступны для контакта от края до края. 4. Зазор рассчитывается автоматически (аналогично переменному зазору, I зазор = 1 ТИПА 7) для каждого воздействия как: :math:`gap=g_{m}+g_{s}` Пока, - :math:`g_{m}` : разрыв основного элемента: :math:`g_{m}=\frac{t}{2}` , с :math:`t` это толщина основного элемента для элементов оболочки :math:`g_{m}=0` , для кирпичных элементов - :math:`g_{s}` : разрыв вторичного узла: :math:`g_{s}=0` , если вторичный узел не соединяется с любым элементом или соединяется только с кирпичом или пружиной элементы. :math:`g_{s}=\frac{t}{2}` , если вторичный узел соединенный с элементом оболочки, с :math:`t` является наибольшей толщиной элементы оболочки, соединенные со вторичным узлом. :math:`g_{s}=\frac{\sqrt{S}}{2}` , если вторичный узел соединенные с элементами фермы или балки, с :math:`S` являющийся поперечным сечением 1D элемент. Если вторичный узел подключен к нескольким оболочкам и/или балок или ферм, используется наибольший рассчитанный вторичный зазор. :math:`g_{m}` и :math:`g_{s}` ограничиваются отдельно Gap_max_m и Gap_max_s перед зазор рассчитывается. Ограничение относительно Igap0=1: Флаг модификации пробела для вторичных узлов оболочки на свободные края не влияет, если вторичный узел определен через необязательная группа узлов (grnd_IDs). 5. If fric_ID определено, контактное трение определяется в /FRICTION и входы трения ( I Фриц , C 1 и т. д.) в этом входе карты не используются. Силы трения являются: :math:`F_{t}^{new}=min⁡(\muF_{n},F_{adh})` Сила сцепления рассчитывается как: :math:`F_{adh}=F_{t}^{old}+\Delta F_{t}` с :math:`\Delta F_{t}=K⋅V_{t}⋅dt` Где, :math:`\mu` представляет собой коэффициент кулоновского трения и определяется как: - Для флага I Фриц по умолчанию: :math:`\mu=Fric` с :math:`F_{T}\le\mu⋅F_{N}` (Кулоновское трение) - Для флага I Фриц > 1, новое трение представлены модели. При этом коэффициент трения задается по функции: :math:`\mu=\mu(p,V)` Где, :math:`p` Давление нормальной силы на главный сегмент :math:`V` Тангенциальная скорость вторичного узла относительно в основной сегмент В настоящее время коэффициенты C1–C6 используются для определения переменный коэффициент трения :math:`\mu` для новых формул трения. Следующие доступны составы: - I Фриц = 1 (Обобщенный закон вязкого трения): :math:`\mu=Fric+C_{1}⋅p+C_{2}⋅V+C_{3}⋅p⋅V+C_{4}⋅p^{2}+C_{5}⋅V^{2}` - I Фриц = 2 (Модифицированный закон Дармстада): :math:`\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 (закон Ренарда): :math:`\mu=C_{1}+(C_{3}−C_{1})⋅\frac{V}{C_{5}}⋅(2−\frac{V}{C_{5}})` если :math:`V\in[0,C_{5}]` :math:`\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}}))` если :math:`V\in[C_{5},C_{6}]` :math:`\mu=C_{2}−\frac{1}{\frac{1}{C_{2}−C_{4}}+(V−C_{6})^{2}}` если :math:`V\geC_{6}` Где, :math:`C_{1}=\mu_{s}` :math:`C_{2}=\mu_{d}` :math:`C_{3}=\mu_{max}` :math:`C_{4}=\mu_{min}` :math:`C_{5}=V_{cr1}` :math:`C_{6}=V_{cr2}` - Первая критическая скорость :math:`V_{cr1}=C_{5}` должно отличаться от 0 ( :math:`C_{5}\ne0` ). - Первая критическая скорость :math:`V_{cr1}=C_{5}` должно быть меньше, чем вторая критическая скорость :math:`V_{cr2}=C_{6}(C_{5} 0, условие устойчивости жесткости :math:`K_{msdt}` используется, если оно больше :math:`K_{n}` , жесткость определяется флагом Istf (см. комментарий 2). :math:`K=\frac{maxSt_{min},minSt_{max},max(K_{n},K_{msdt})}{2}` По умолчанию, в случае контактной жесткости, основанной на условие устойчивости (если Istf = 7 или Ipstif > 0), узловые массы, используемые для :math:`K_{msdt}` первоначально рассчитываются с использованием материала плотность и размер элементов. Он не учитывает добавленную массу и может привести низкой жесткости и контактной неустойчивости. В этом случае Ipstif= 2 можно использовать для рассмотрения изначально добавлен масс эффект. 1 Шули, Франц и Ив Ренар. «Явное интегрирование по времени Верле для аппроксимации упругодинамики, основанной на Ниче. проблемы с контактами». Расширенное моделирование и симуляция в технике Науки 5, no. 1 (2018): 31