Расчет нестационарной гидродинамической нагрузки на крышку пускового устройства …
9
arccos
;
i
j
k
a a
sin 0,5
;
sin
j
k i
i
i
b
2arcsin
i
j k
b b
и индексы
i
,
j
и
k
образуют циклическую перестановку из индек-
сов 1, 2 и 3.
Составление СЛАУ, порядок вычислений.
Искомыми величи-
нами ГИУ (5) относительно
являются значения потенциала
M
в центрах плоских треугольных элементов, покрывающих
часть поверхности крышки в октанте
0,
0,
0
x y z
. Коэффициен-
ты при этих значениях образуют матрицу СЛАУ. Вектор правых ча-
стей СЛАУ вычисляется в результате интегрирования
и
n
вдоль
границ меридионального сечения поверхностей ЛА и ГП как объяс-
нено выше. Для численного решения СЛАУ использовался стандарт-
ный прямой метод [12].
Вычисленные значения
вдоль поверхности крышки численно
дифференцируются в касательном направлении для вычисления
квадрата скорости
2
v
частиц жидкости на поверхности крышки в
формуле (2) для вычисления избыточного давления. В быстропроте-
кающем нестационарном процессе эта составляющая давления имеет,
вообще говоря, второстепенное значение.
Основное значение имеет решение ГИУ (5) относительно
.
t
Искомыми величинами являются значения
M
t
в центрах тре-
угольных ГУ. Из ГИУ (5) следует, что матрица СЛАУ получается той
же, что и в случае
. Вектор правых частей второй СЛАУ вычис-
ляется в результате интегрирования
t
и
n t
вдоль границ
меридионального сечения поверхностей ЛА и ГП.
