Исследуем приземление самолета с шасси на воздушной подушке

нет производных в правой части (нормальный вид).

Последовательность

В нормальном виде система дифференциальных уравнений 8 порядка:

ВВ’ – задается (для “аккорда-2001” ВВ’=1,9м)

где

Уравнение (3) в виде, не готовом для интегрирования

a11 = 1,

a12 = 0,

a13 = ,

a14 = ,

a15 = 0,

a16 = 0,

b1 = .

4)

a21 = - r1,

a22 = r1 – r2,

a23 = Рб – Рвп,

a24 = 0,

a25 = - Рб,

a26 = 0,

b2 = 0.

5)

a31 = 0,

a32 = 0,

a33 = sinφ1,

a34 = r1 cosφ1,

a35 = sinφ2,

a36 = r2 cosφ2,

b3 = 0.

6)

a41 = 0,

a42 = 0,

a43 = 1 – cosφ1,

a44 = r1sinφ1,

a45 = - (1 – cosφ2),

a46 = - r2sinφ2,

b4 = 0.

7)

a51 = 0,

a52 = - ,

a53 = φ1,

a54 = r1,

a55 = ,

a56 = r2,

b5 = 0.

8)

a61 = 0,

a62 = ,

a63 = ,

a64 = ,

a65 = ,

a66 = ,

b6 = 0.

Алгоритм решения (до контакта).

Исходные данные:

H0, , Рвп0, Рб0, r10, r20, φ10, φ20 – начальные условия.

H0 = 1м,

= 1м,

Н = -3м/с,

Рб0 = 1000Н/м2,

Ризб = Рвп = 0,

r1 = r2 = 0.21м,

φ1 = φ2 = 2.3рад,

ТЕ = 100000Н/м2.

Геометрические данные и массовые данные: S0, BB’, m.

m = 270 кг/м,

ВВ’ = 1.9м,

S0 = 0.123375236м2,

L0нат = 0,96629434м.

Параметры состояния:

Ратм = 100000Н/м2,

ρатм = 1,25кг/м3.

Интегрируется система дифференциальных уравнений методом Эйлера:

Hi+1 = Hi +∆t

Hi+1 = Hi +∆t