Рис. 4. Геометрия исследуемых профилей
Система (4) интегрируется методом Рунге–Кутты четвертого по-
рядка с постоянным шагом
Δ
τ
. Система (5) интегрируется методом
Эйлера первого порядка с шагом
Δ
t
с проверкой условия “проника-
ния”: вихревые элементы, попавшие внутрь профиля между контроль-
ными точками, возвращаются в поток. Цикл вычислений продолжается
до достижения конечного времени вычислений
t
K
.
На основе описанного алгоритма автором настоящей статьи создан
программный комплекс MDVDD, при помощи которого были прове-
дены расчеты переходных режимов.
Результаты моделирования.
Исследование проводилось на мо-
дельных задачах с безразмерными параметрами. В качестве выдвигае-
мого элемента конструкции был взят прямоугольный профиль, закруг-
ленный с одного торца; его размеры показаны на рис. 4, на котором
также показаны профили и размеры оснований.
Использовались следующие безразмерные исходные данные:
|
~v
∞
|
= 1
;
p
∞
= 1
,
ρ
∞
= 1
;
v
0
Y
= 0
,
1
;
l
(0) = 1
,
5
;
Y
0max
= 1
,
3
;
EJ
= 5
;
ρ
= 3
,
125
;
C
1
=
C
2
= 100
;
ε
i
= 0
,
004
;
Δ
t
= 0
,
04
;
Δ
τ
= Δ
t/
100
;
t
K
= 13
;
N
K
= 2000
. . .
2200
(по мере выдвижения профиля из осно-
вания число панелей увеличивается).
Были проведены последовательно расчеты нескольких схем, обо-
значенных далее как:
схема 1 — профиль выдвигается из бесконечного экрана, на котором
не происходит генерация завихренности [7];
54 ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2008. № 4
