Избранные задачи аэроупругости. Гришанина Т.В - 14 стр.

UptoLike

Рубрика: 

27
α
ρ
=
1M
4
2
2
2
U
p
; (3.3)
2) при
M > 4 по линеаризованной поршневой теории с учетом
влияния толщины симметричных профилей
h(x,z) и начального угла
атаки
0
α :
α
α+
+
+γ
+
ρ
=
2
0
2
2
M
4
M
M
4
1
1
M
4
2 x
h
x
hU
p
. (3.4)
Здесь
= aU /M число Маха;
a скорость звука в невозмущен-
ном потоке воздуха;
4.1
=
γ
показатель адиабаты.
В этих случаях приращение давления пропорционально при-
ращению местного угла атаки в той же точке и оператор P имеет вид
коэффициента:
1M
4
2
=
P
при 4Mcos/1 <<χ ;
α+
+
+γ
+=
2
0
2
M
4
M
M
4
1
1
M
4
x
h
x
h
P
при M > 4. (3.5)
Гипотезу стационарности можно использовать, если частота
колебаний
ω достаточна мала, так что выполняется условие [5]
2
π
<<
ω
U
b
. (3.6)
Уравнения колебаний стабилизатора в потоке могут быть со-
ставлены на основании уравнений Лагранжа в обобщенных коорди-
натах
)(
1
tq
,
)(
2
tq
:
28
,2,1,0])([
2
1
==+++
=
iqbkqdqm
j
j
ijijjijjij
&&&
(3.7)
где
dSmm
S
jiij
∫∫
ϕϕ= ;
dSP
U
d
S
jiij
∫∫
ϕϕ
ρ
= ][
2
;
dS
x
P
Ua
bbMb
S
j
iijijij
∫∫
ϕ
ϕ
ρ
==
2
;
. (3.8)
3.2. Вычисление коэффициентов
Рассмотрим сплошной однородный стабилизатор с ромбовид-
ной формой профиля постоянной относительной толщины
const/
=
= bcc , hm
0
ρ
=
;
0
ρ
плотность материала (рис. 3.1).
Коэффициенты инерции
ij
m :
;
48
1
)tg)((
2
1
4
0
22
00
2
11
dzbcdzbezdcJdSmxm
le
e
le
e
z
S
∫∫
++
ρ+χ+ρ===
zdzbezdcJdSmxzmm
le
e
xz
S
2
002112
)tg)((
2
1
∫∫
+
χ+ρ==== ;
(3.9)
dzzbcJdSmzm
le
e
x
S
22
0
2
22
2
1
∫∫
+
===
ρ
.
При вычислении аэродинамических коэффициентов
ij
d и
ij
b
учтем, что оператор
P в рассматриваемом случае является коэффици-
ентом, который определяется по формулам (3.5).
Коэффициенты аэродинамического демпфирования
ij
d :
                                      ρU 2      4                                  2
                               ∆p =
                                       2
                                                     α;        (3.3)              ∑ [m q&&     ij   j   + d ij q& j + (k ij + bij )q j ] = 0,           i = 1, 2,          (3.7)
                                              M −1
                                                2
                                                                                   j =1


         2) при M > 4 по линеаризованной поршневой теории с учетом                                                                  ρU
                                                                       где      mij = ∫∫ mϕi ϕ j dS ;                      d ij =        ∫∫ ϕ P[ϕ ]dS ;
                                                                                                                                              i     j
влияния толщины симметричных профилей h(x,z) и начального угла                             S
                                                                                                                                     2   S

атаки α 0 :                                                                                                          ρUa∞               ⎡ ∂ϕ j ⎤
                 ρU 2 4 ⎧⎪ γ + 1 ⎡ ∂h M ⎛ ∂h ⎞    ⎤ ⎫⎪    2                      bij = Mbij ;                bij =
                                                                                                                       2       ∫∫S ⎢⎣ ∂x ⎥⎦dS .
                                                                                                                                  ϕ i P                                    (3.8)
            ∆p =         ⎨1 +   M ⎢ + ⎜ ⎟ + Mα 02 ⎥ ⎬α . (3.4)
                  2 M ⎪⎩      4   ⎣⎢ ∂x 4 ⎝ ∂x ⎠  ⎦⎥ ⎪⎭
                                                                                                    3.2. Вычисление коэффициентов
Здесь M = U / a∞ – число Маха; a ∞ – скорость звука в невозмущен-
                                                                               Рассмотрим сплошной однородный стабилизатор с ромбовид-
ном потоке воздуха; γ = 1.4 – показатель адиабаты.
                                                                       ной    формой                профиля          постоянной              относительной          толщины
         В этих случаях приращение давления пропорционально при-
                                                                       c = c / b = const , m = ρ0 h ; ρ 0 – плотность материала (рис. 3.1).
ращению местного угла атаки в той же точке и оператор P имеет вид
                                                                               Коэффициенты инерции mij :
коэффициента:
                                                                                                                     e+l                                             e+l
                  4                                                                            1                                      1
       P=                    при 1 / cos χ < M < 4 ;                   m11 = ∫∫ mx 2 dS = J z = ρ0c ∫ (d 0 + ( z − e) tg χ) 2 b 2 dz + ρ0c ∫ b 4 dz;
               M2 − 1                                                         S
                                                                                               2    e
                                                                                                                                      48   e

                                                                                                                                              e+l
          4 ⎧⎪ γ + 1 ⎡ ∂h M ⎛ ∂h ⎞
                                     2
                                       ⎤ ⎫⎪                                                                     1
       P = ⎨1 +     M ⎢ + ⎜ ⎟ + Mα 02 ⎥ ⎬ при M > 4. (3.5)                   m12 = m21 = − ∫∫ mxzdS = − J xz = − ρ0c ∫ (d 0 + ( z − e) tg χ)b 2 zdz ;
                                                                                                                2
          M ⎪⎩   4    ⎣⎢ ∂x 4 ⎝ ∂x ⎠   ⎦⎥ ⎪⎭
                                                                                            S                        e

                                                                                                                                                                           (3.9)
         Гипотезу стационарности можно использовать, если частота                                                              e+l
                                                                                                                     1
колебаний ω достаточна мала, так что выполняется условие [5]                 m22 = ∫∫ mz 2 dS = J x =                  ρ 0 c ∫ b 2 z 2 dz .
                                                                                       S
                                                                                                                     2       e
                                             ωb   π
                                                << .           (3.6)
                                             U    2                            При вычислении аэродинамических коэффициентов d ij и bij
                                                                       учтем, что оператор P в рассматриваемом случае является коэффици-
         Уравнения колебаний стабилизатора в потоке могут быть со-
                                                                       ентом, который определяется по формулам (3.5).
ставлены на основании уравнений Лагранжа в обобщенных коорди-
                                                                              Коэффициенты аэродинамического демпфирования d ij :
натах q1 (t ) , q 2 (t ) :

                                                                 27    28