长直机翼及其坐标系如图2-1所示,其半展长为 l ,半弦长为 b 。取固支端与机翼刚轴的交点为原点建立坐标系, y 轴沿机翼轴线从翼根指向翼尖, x 轴沿机翼弦向由前缘指向后缘,与 y 轴正交, z 轴与 x 、 y 轴构成右手坐标系。在此坐标系下,机翼的弯扭振动微分方程可写为 [48]
其中,
h
为机翼弯曲振动位移;
α
为机翼扭转振动转角;
EI
为机翼抗弯刚度;
GJ
为机翼抗扭刚度;
m
为机翼单位长度质量;
I
α
为单位长度机翼绕弹性轴的转动惯量;
L
h
为机翼单位长度的升力;
T
α
为机翼单位长度的扭矩;
y
为机翼展向坐标值;
t
为时间。机翼剖面图如图2-2所示,
为机翼弹性轴到
z
轴的距离,
x
α
为机翼弹性轴到机翼横截面重心的距离。
图2-1 长直机翼及其坐标系
图2-2 机翼剖面图
在忽略机翼重力影响的条件下,机翼颤振时的外力只有气动力。本书采用片条理论进行非定常气动力计算。根据Theodorson理论,单位展长的非定常升力与相应的俯仰力矩按式(2-2)计算:
其中,
V
为空速;
ρ
为空气密度;
C
(
k
)为Theodroson函数;
k
为减缩频率,
k
=
ωb
/
V
;
ω
为圆频率;
为机翼弹性轴到机翼弦长中点的距离占半弦长的百分比;其他变量说明同式(2-1)。
将式(2-2)代入式(2-1),得到机翼颤振微分方程: