/ALE/GRID/STANDARD

Ключевое слово формата блока Описывает стандартную формулировку сетки ALE.

расчет скорости.

Это улучшенный вариант /ALE/GRID/SPRING рецептура на основе краевых пружин

и противосрезные пружины.

1

Формат

/ALE/GRID/STANDARD	/ALE/GRID/STANDARD	/ALE/GRID/STANDARD	/ALE/GRID/STANDARD	/ALE/GRID/STANDARD	/ALE/GRID/STANDARD	/ALE/GRID/STANDARD	/ALE/GRID/STANDARD	/ALE/GRID/STANDARD	/ALE/GRID/STANDARD
\(\alpha\)	\(\alpha\)	\(\gamma\)	\(\gamma\)	\(\eta\)	\(\eta\)	\(l_{c}\)	\(l_{c}\)
Пустой формат	Пустой формат	Пустой формат	Пустой формат	Пустой формат	Пустой формат	Пустой формат	Пустой формат	Пустой формат	Пустой формат

Определение

Поле	Содержание	Пример единицы СИ
\(\alpha\)	Масштабный коэффициент для максимальной жесткости. 2 По умолчанию = 0,9 (реальное)
\(\gamma\)	Коэффициент нелинейности для краевой пружины жесткость. 3По умолчанию = 1e-2 (Реальное)
\(\eta\)	Коэффициент демпфирования. 4По умолчанию = 1e-2 (Реальное)
\(l_{c}\)	Характеристика длина.(Реальная)	\([m]\)

Комментарии

1. Для управления сеткой на твердых элементах введены фиктивные пружины.

   скорости.

Эти пружины являются нелинейно-упругими вязкими. Чтобы обеспечить стабильность, их

жесткость рассчитывается по временному шагу. Два типа пружин: краевые и противосрезные.

пружины.

1. Краевые пружины

Силы, действующие на краевую пружину, являются функцией изменения ее длины.

во время.

\(\Delta F_{edge}=k(h)⋅(w_{2}−w_{1})dt\) Где, \(w_{1},w_{2}\)

— скорости сетки в узлах Н1 и Н2, соответственно.

\(h\)

расстояние N1 с противоположного лица

\(dt\)

это временной шаг

и

\(k(h)\) это жесткость пружины \(k(h)=k_{critical}\)

Если

\(h\) меньше характерной длины \(l_{c}\) и N1 движется тогда к противоположному лицу,

\(k(h)=\frac{1}{\lambda^{2}}[\gamma+(\gamma−1)(\frac{h−l_{c}}{l_{c}})^{3}]k_{critical}\) \(\frac{1}{\lambda^{2}}\) – коэффициент устойчивости с учетом коэффициента демпфирования \(\beta\) , масштабный коэффициент \(\alpha\) и шаг по времени \(dt\) 4 .. image:: images/ale_grid_standard_starter_r_ale_standard_towards.png

alt

ale_standard_towards

(Рисунок 2.)

в противном случае, \(k(h)=\frac{\gamma}{\lambda^{2}}k_{critical}\) .. image:: images/ale_grid_standard_starter_r_ale_standard_away.png

alt

ale_standard_eq

(Рисунок 3.)

2. Противосрезные пружины

Противосдвигающие силы

\(F_{shear}\) рассчитываются на основе проникновения узла. разрыв \(\frac{l_{c}}{s}\) с противоположного лица.

Стоимость

\(F_{shear}\) это:

\(F_{shear}=gap⋅k(h)\) и \(k(h)=\frac{\gamma}{\lambda^{2}}[\gamma+(\gamma−1)(\frac{h−l_{c}}{l_{c}})^{3}]k_{critical}\) .. image:: images/ale_grid_standard_starter_r_ale_standard_critical.png

alt

ale_standard_critical

(Рисунок 4.)

3. Вязкое демпфирование

Вязкие силы рассчитываются на основе критического демпфирования, соответствующего

до верхней границы жесткости: \(\frac{1}{\lambda^{2}}\)

\(F_{viscous}=\beta\alpha(\sqrt{1+\beta^{2}}−\beta)(w_{2}−w_{1})\)

4. Скорость сетки

Затем скорость сетки обновляется в соответствии с: \(w_{n+1}=w_{n}+\frac{(\Delta F_{edge}+F_{shear}+F_{viscous})dt}{m}\) Где,

\(m\) фиктивная масса на узле из пружин (автоматически вычисляется во время Стартера).

2. Увеличение

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

длина. Масштабный коэффициент гарантирует, что критическое значение стабильности не будет превышено (во избежание

уменьшение шага по времени).

3. Этот флажок действует на жесткость

   форма. Жесткость линейна с

\(\gamma\) = 0 . Более того, увеличение \(\gamma\) , нижняя граница жесткости краевой пружины будет увеличена. Пружины

иметь критическую жесткость при нулевой длине (это соответствует унитарному коэффициенту). Для

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

критическая жесткость, умноженная на

\(\gamma\) .

4. Рекомендуется использовать небольшие

   ценности для

\(\beta\) , в противном случае демпфирование может стать слишком критическим. Фактор стабильности

это:

\(\lambda=\frac{dt}{\alpha_{0}(\sqrt{1+\beta^{2}}−\beta)}\)

5. \(l_{c}\)

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

\(h<l_{c}\)

• Противосрезная пружина активируется: \(h<\frac{l_{c}}{5}\)

6. Все эти параметры могут быть

изменено во время перезапуска двигателя (/ALE/GRID/STANDARD).

7. Автокоррекция сетки. Это

   можно придать больший вес противосдвигающим силам одним из следующих способов:

• Параметр \(l_{c}\) параметр, близкий к размеру сетки

• Установка отрицательного значения для \(\alpha\) параметр (упругие силы на краях установлены равными 0 при

  первый цикл текущего запуска)

8. Этот метод предполагает однородность

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