9、飞行力学第四章-运动方程

合集下载

9、飞行力学第四章-运动方程

9、飞行力学第四章-运动方程

共84页
飞行器飞行力学
18
四、在气流坐标系中的平动动力学方程
根据速度之间的关系
u = V cos α cos β ⎫ ⎪ v = V sin β ⎬ w = V sin α cos β ⎪ ⎭
可得
du dV dα dβ ⎫ V sin α cos β − V cos α sin β ⎪ cos α cos β − = dt dt dt dt ⎪ dv dV dβ ⎪ V cos β sin β + = ⎬ dt dt dt ⎪ dw dV dα dβ ⎪ sin α cos β + V cos α cos β − V sin α sin β ⎪ = dt dt dt dt ⎭
共84页
飞行器飞行力学
5
引言(续)
本章主要工作 (1) 应用牛顿第二定律、动量矩定律和运动学原理, 导出飞行器相对于动坐标轴系的一般运动方程组; (2) 在小扰动假设前提下,对非线性运动方程组进行 线性化; (3) 在特定条件下,将纵向运动和横航向向小扰动运 动分开,得到飞行器的纵向和横航运动方程组。
共84页 飞行器飞行力学 4
引言(续)
假设 1、飞行器为面对称刚体 2、静止地球与平面大地假设 3、标准大气假设 4、忽略旋转部件及液体晃动的影响。 运动方程 采用刚体六自由度数学模型,包括六个动力学方程: 其中三个用来描述质心的平动,三个描述飞行器绕质 心的转动。另外,还有六个运动学方程,分别用来描 述飞行器在空间的位置和姿态。
固连于飞行器的任意动系中质心的转动动力学方程
dh δh = +ω×h dt δt
得转动动力学方程为
⎡ 0 ⎢ 其中 ω × h = ⎢ ω z ⎢− ω y ⎣

飞行力学第四章

飞行力学第四章

解出
dV / dt, d / dt, d / dt
在、 不大的快速机动中,可近似认为速度不变,且
u V const v V w V
从而有
T cos Z b g d q p cos cos dt mV V Yb d g p r sin cos dt mV V
3、重力
2、气动力
Ax D Ay C L Az a
W x 0 W y m 0 g W z
g
航迹轴 系 ,
(无风时)
气流轴 系 ,
,, dVb Vb b Vb dt t Tb Lba Aa LbgW g
引言
研究内容
性能
飞行品质
质点动力学系统
质点系动力学系统
操纵性与稳定性:研究飞机在外力和外 力矩作用下运动参数的变化特性。
稳定性
指飞行器在受到外界瞬时扰动后,是否具有自动地恢复到原来平 衡状态的能力。
操纵性
指飞行器对驾驶员操纵或舵面指令输入的响应,即从一种飞行 状态过渡到另一种飞行状态的能力.
Nanjing University of Aeronautics and Astronautics
第四章
刚性飞行器运动方程
内容
引言 4.1 刚性飞行器动力学方程 4.2 刚性飞行器运动学方程 4.3 刚性飞行器运动方程讨论* 4.4 运动方程组线性化 4.5 纵向小扰动运动方程组 4.6 横航向小扰动运动方程组 小结
矩阵 形式
hx I x hy I xy h z I zx

飞行动力学飞机方程

 飞行动力学飞机方程

xydm Ixy
表示惯性积
依据假设 Ixy=Izy=0 ,H 的各分量
H
x
H y
pI x qI y
rI xz
代入
dH dt
1H
dH dt
H
H x dt
pI x rI xz
dH y dt
qI y
dH z dt
rI z pI xz
由于
i jk
H p q r i(qH z rH y ) j(rH x pH z ) k ( pH y qH x )
1.地轴系与机体轴系间的方向余弦表
o
xg
x
cos cos
y
cos sin sin- sincos
z
cos sin cos+sinsin-
yg sincos sin sin sin+cos cos sin sin cos-cos sin
zg -sin cos sin cos cos
表中,oxyz为机体轴系, oxgygzg为地轴系
—动坐标系对惯性系的总角速度向量
—表示叉积,向量积
1H —沿动量矩 H 的单位向量
dV , dH dt dt
—对动坐标系的相对导数
1.力方程
F
m
dV dt
dV dt
1V
dV dt
V
V 和 用机体坐标系上的分量(u,v,w;p,q,r)表示
V iu jv kw, ip jq kr
三个力方程 三个力矩方程 飞机六自由 度动力学
线性方程 增量方程
m
d u dt
( X u
)0 u
( X
)0
( X
)0
m

飞行动力学公式总结

飞行动力学公式总结

飞行动力学第二章公式总结空气动力:X=C x qS阻力公式Y=C y qS升力公式Z=C z qS侧向力公式动压公式q=ρV22升力:C y=f(Ma,α,δ)升力系数函数C y=C y0+C yαα+C yδzδz升力系数在攻角和舵偏角不大的情况下的表达式C y=C yαα+C yδzδz轴对称时Y=Y0+Yαα+Yδzδ升力在攻角和舵偏角不大的情况下的表达式α攻角不大情况下攻角变化引起的升力Yα=C yαρV22Yδ=C yδzρV2δz舵偏角不大的情况下舵偏角变化引起的升力2侧向力:C z=C zββ+C zδzδz侧向力因数在侧滑角和舵偏角不大的情况下的表达式-C zβ=C yα轴对称下成立(不大)-C yδz=C zδz轴对称下成立(不大)阻力:X= X0+X i阻力的组成由零升阻力和诱导阻力构成C x=C x0+C x i阻力因数由零升阻力因数和诱导阻力因数构成气动力矩:M x1=m x1qSL滚转力矩M y1=m y1qSL偏航力矩M z1=m z1qSL俯仰力矩M z =f(M a ,H,α,δz ,,ωz ,α̇, δz ) 俯仰力矩的函数M z = M z 0+M z αα+M z δz δz+ M z ωz ωz+ M z αα̇+M z δz δz参数不大的情况下升力表达式 m z = m z 0+m z αα+m z δz δz+ m z ωz ̅̅̅̅ωz ̅̅̅̅+ m z α̅α̇̅+m z δz ̅̅̅̅δz̅ 无量纲力矩因数表达式 δz ̅=δzL/V 舵偏角角速度对应的无量纲参数 α̇̅=α̇L/V 攻角角速度对应的无量纲参数 ωz ̅̅̅̅=ωzL/V 俯仰角角速度对应的无量纲参数M z α=C z αSqα(x g −x F )=m z αSqαL 升力力矩和里表达式之间的关系m z α=C z α(X g ̅̅̅−X F ̅̅̅̅) 攻角升力系数和攻角升力力矩系数之间的关系 m z δz =C z δz (X g ̅̅̅−X r ̅̅̅) 舵偏角升力系数和舵偏角升力力矩系数之间的关系 m z =m z αα+m z δz δz 轴对称定常直线飞行下的升力力矩系数表达式m z ααb +m z δz δz=0 "瞬时平衡假设"下的升力力矩平衡状态方程C b y =C b ααb +C b δz δzb =(C b α−C b δz m z αm z δz )αb “瞬时平衡”状态下平衡升力的表达式m z α|α=αb <0 纵向静稳定条件m z C y =ðm zðC y =(X g ̅̅̅−X F ̅̅̅̅) 稳定性的定量表示——静稳定度 ∆α=arctanrωz V 俯仰角角速度引起的下洗角度 M z ωz =M z ω̅z ω̅z qSL 俯仰阻尼力矩表达式t t t αεεαα•∆()=(()-)实际下洗角 偏航力矩:m y =m y ββ+m y δy δy +m y ω̅y ω̅y +m y ω̅x ω̅x +m y δ̅y δy +m y β̅β 偏航力矩系数表达式 ω̅y =ωy L/V偏航角速度对应的无纲量因数 δy=δy L/V 航向舵偏角速度对应的无纲量因数 β=βL/V 偏航角角速度对应的无量纲因数m x =m x0+m x ββ+m x δy δy +m x δx δx +m x ω̅x ω̅x +m x ω̅y ω̅y 滚转力矩因数的表达式 m x ββ<0 横向静稳定性的条件M ℎ=m ℎq t S t b t 铰链力矩模式表达式M ℎ=−Y t ℎcos(α+δz ) 铰链力矩实际表达式M ℎ≈M ℎαα+M ℎδz δz 铰链力矩的近似表达式 推力:P =m s μe +S a (P a −P ℎ) 推力的表达式 M p =R p ×P 推力力矩表达式重力:G=G 1+F e 重力表达式F e =mR e Ωe 2cosψe 离心惯性力的表达式 g =g 0R e 2(R e +H e )2 重力加速度随高度变化的表达式导弹建模基础:m dV dt =F质心移动的动力学公式 dH dt =M 绕质心转动的动力学公式导弹质心移动的动力学方程:m dV dt =m (ðV ðt +Ω×V)=F 用相对坐标系表示以绝对坐标系为基准的矢量变化率表示-力 ρ=V θ 曲率半径的计算公式a y2=Vθ 弹道法线加速度 导弹绕质心转动的动力学方程:dH dt =ðH ðt +ω×H =M用相对坐标系表示以绝对坐标系为基准的矢量变化率表示-力矩 H =J ∙ω动量矩M =J ∙α力矩 J ={J x1−J x1y1−J z1x1−J x1y1J y1−J y1z1−J z1x1−J y1z1J z1} 三维空间下转动惯量矩阵 dm dt =−m s (t)导弹质量流率方程 m =m 0−∫m s (t)dt tf t0 导弹质量方程角度几何关系:cosφ=cosα1cosα2+cosβ1cosβ2+cosγ1cosγ2 余弦定理α=ϑ−θ 无滚转无侧滑角度关系时β=ψ−ψv 无攻角无滚转时角度关系操纵关系方程:N =P +R 控制力为空气动力与推力的合力N =N n +N τ 控制力的切向与法向的分解N τ=P τ−X 切向控制力分解 N n =P n +Y +Z 法向控制力分解导弹飞行的运动方程组(轴对称型导弹,以地面为绝对坐标系): 质心移动的动力学方程(弹体->弹道坐标系):m dV dt =Pcosαcosβ−X −mgsinθ切向运动的动力学方程 mV dθdt =P (sinαcosγv +cosαsinβsinγv )+Ycosγv −Zsinγv −mgcosθ 竖直法向运动的动力学方程 −mVcosθdψv dt =P (sinαsinγv −cosαsinβcosγv )+Ysinγv +Zcosγv 水平法向运动的动力学方程 绕质心转动的动力学方程(弹体坐标系):J xdωx dt +(J z −J y )ωy ωz =M x 弹体x 轴力矩表达式 J ydωy dt +(J x −J z )ωz ωx =M y 弹体y 轴力矩表达式 J z dωz dt +(J y −J x )ωx ωy =M z 弹体z 轴力矩表达式质心移动的运动学方程(弹道->地面坐标系):dxdt=Vcosθcosψv地面坐标系x轴方向运动学方程dydt=Vsinθ地面坐标系y轴方向运动学方程dxdt=−Vcosθsinψv地面坐标系z轴方向运动学方程绕质心转动的运动学方程(弹体->地面坐标系):dϑdt=ωy sinγ+ωz cosγ俯仰角角速度表达式dψdt =1cosϑ(ωy cosγ+ωz sinγ)偏航角角速度表达式dγdt=ωx−tanϑ(ωy cosγ+ωz sinγ)滚转角角速度表达式质量方程:dmdt=−m s角度转换:sinβ=cosθ[cosγsin(ψ−ψv)+sinϑsinγcos(ψ−ψv)]−sinθcosϑsinγ侧滑角用其他角的表达关系cosα=[cosϑcosθcos(ψ−ψv)+sinϑsinθ]/cosβ俯仰角用其他角进行表示cosγv=[cosγcos(ψ−ψv)−sinϑsinγsin(ψ−ψv)]/cosβ速度滚转角的表示控制方程:ε1=0 俯仰方向的控制方程ε2=0 滚转方向的控制方程ε3=0 偏航方向的控制方程ε4=0 速度大小的控制方程描述导弹纵向运动的方程组(忽略z、β、ψ、ψv、ωy、γ、γv、ωx):质心移动的动力学方程:m dVdt=Pcosα−X−mgsinθ纵向平面内沿速度方向的动力学方程mV dθdt=Psinα+Y−mgcosθ纵向平面内速度纵法线方向的动力学方程绕质心转动的动力学方程:J z dωzdt=M z纵向平面内绕弹体z轴旋转的动力学方程质心移动的运动学方程:dxdt=Vcosθ纵向平面水平运动学方程dydt=Vsinθ纵向平面竖直运动学方程绕质心转动的运动学方程:dϑdt=ωz弹体绕z轴的转动质量方程:dmdt=−m s质量变化方程几何关系方程:α=ϑ−θ纵向平面俯仰角、弹道倾角、攻角之间的关系控制方程:ε1=0 俯仰方向的控制方程ε4=0 速度大小的控制方程侧向运动方程组(基于纵向运动方程组):质心移动的动力学方程:−mVcosθdψvdt=P(sinα+Y)sinγv−(Pcosαsinβ−Z)cosγv速度侧法向方向动力学方程绕质心转动的动力学方程:J x dωxdt=M x−(J z−J y)ωzωy绕弹体x轴转动的力矩守恒J y dωydt=M y−(J x−J z)ωxωz绕弹体y轴转动的力矩守恒质心移动的运动学方程:dzdt=−Vcosθsinψv地面坐标系下z轴方向的运动绕质心转动的运动学方程:dψdt =1cosϑ(ωy cosγ−ωz sinγ)偏航方向转动方程dγ=ωx−tanϑ(ωy cosγ−ωz sinγ)滚转方向转动方程dt几何关系方程:sinβ=cosθ[cosγsin(ψ−ψv)+sinϑsinγcos(ψ−ψv)]−sinθcosϑsinγ侧滑角用其他角的表达关系cosγv=[cosγcos(ψ−ψv)−sinϑsinγsin(ψ−ψv)]/cosβ速度滚转角的表示控制方程:ε2=0 侧滑角的控制方程ε3=0 滚转角的控制方程有侧滑无倾斜的水平运动方程组:条件:θ=0弹道倾角为零γ=γv=0滚转角为零ωx=0滚转角速度为零质心移动的动力学方程(弹体->弹道坐标系):=Pcosαcosβ−X切向运动的动力学方程m dVdtPsinα+Y=mg竖直法向运动的动力学方程−mVcosθdψv=−Pcosαsinβ+Z水平法向运动的动力学方程dt绕质心转动的动力学方程(弹体坐标系):=M y弹体y轴力矩表达式J y dωydt=M z弹体z轴力矩表达式J z dωzdt质心移动的运动学方程(弹道->地面坐标系):dx=Vcosψv地面坐标系x轴方向运动学方程dtdx=−Vsinψv地面坐标系z轴方向运动学方程dt绕质心转动的运动学方程(弹体->地面坐标系):dϑdt=ωz俯仰角角速度表达式dψdt =ωycosϑ偏航角角速度表达式质量方程:dmdt=−m s角度转换:α=ϑ俯仰方向角度关系β=ψ−ψv偏航方向角度关系控制方程:ε2=0 偏航方向的控制方程ε4=0 速度大小的控制方程导弹的质心运动:条件:m zααb+m zδzδzb=0攻角方向的力矩守恒m yββb+m yδyδyb=0侧滑角方向的力矩守恒ε1=0 ε2=0 ε3=0 ε4=0 俯仰、侧滑、滚转、速度方向上实现理想控制质心移动的动力学方程(弹体->弹道坐标系):m dVdt=Pcosαb cosβb−X b−mgsinθ切向运动的动力学方程mV dθdt=P(sinαb cosγv+cosαb sinβb sinγv)+Y b cosγv−Z b sinγv−mgcosθ竖直法向运动的动力学方程−mVcosθdψvdt=P(sinαb sinγv−cosαb sinβb cosγv)+Y b sinγv+Z b cosγv水平法向运动的动力学方程质心移动的运动学方程(弹道->地面坐标系):dxdt=Vcosθcosψv地面坐标系x轴方向运动学方程dydt=Vsinθ地面坐标系y轴方向运动学方程dxdt=−Vcosθsinψv地面坐标系z轴方向运动学方程质量方程:dmdt=−m s描述导弹质心铅锤平面内运动方程组:质心移动的动力学方程:m dVdt=Pcosα−X−mgsinθ纵向平面内沿速度方向的动力学方程mV dθdt=Psinα+Y−mgcosθ纵向平面内速度纵法线方向的动力学方程质心移动的运动学方程:dxdt=Vcosθ纵向平面水平运动学方程dydt=Vsinθ纵向平面竖直运动学方程质量方程:dmdt=−m s质量变化方程几何关系方程:δzb=−m zαm zδzαb控制方程:ε1=0 俯仰方向的控制方程ε4=0 速度大小的控制方程导弹质心在水平面内的运动方程组:条件:θ=0弹道倾角为零γ=γv=0滚转角为零ωx=0滚转角速度为零α->0攻角很小β->0侧滑角很小质心移动的动力学方程(弹体->弹道坐标系):=P−X b切向运动的动力学方程m dVdtPαb+Y=mg竖直法向运动的动力学方程−mVcosθdψv=−Pβb+Z b水平法向运动的动力学方程dt质心移动的运动学方程(弹道->地面坐标系):dx=Vcosψv地面坐标系x轴方向运动学方程dtdz=−Vsinψv地面坐标系z轴方向运动学方程dt质量方程:dm=−m sdt角度转换:ψ=ψv+βb偏航角、速度滚转角、侧滑角水平飞行时的几何关系ϑ=α水平飞行时俯仰角和攻角之间的几何关系m zααb+m zδzδzb=0攻角方向的力矩守恒m yββb+m yδyδyb=0侧滑角方向的力矩守恒控制方程:ε2=0 滚转方向的控制方程ε4=0 速度大小的控制方程过载:过载矢量的定义n=NGF i=nG i通过过载来求导弹任意部分的外力大小过载的投影:(Pcosαcosβ−X)速度坐标系x轴方向过载的投影n x3=1Gn y3=1(Psinα+Y)速度坐标系y轴方向过载的投影Gn z3=1G(Pcosαcosβ+Z)速度坐标系z轴方向过载的投影n x2=1G(Pcosαcosβ−X)弹道坐标系x轴方向过载的投影n y2=1G(cos(γv) (sin(α) P + Y) − sin(γv) (−sin(β) cos(α) P + Z))弹道坐标系y轴方向过载的投影n z2=1G(sin(γv) (sin(α) P + Y) + cos(γv) (−sin(β) cos(α) P + Z))弹道坐标系z轴方向过载的投影过载表示动力学方程:m dVdt=N x2+G x2沿速度方向的动力学方程mV dθdt=N y2+G y2沿速度法向纵向对称面内的动力学方程−mVcosθdψvdt=N z2+G z2沿速度法向横向动力学方程用V、θ、ψv来表示过载:n x2=1gdVdt+sinθn y2=Vgdθdt+cosθn z2=−Vgdψvdtcosθ根据过载判断飞行状态:n x2=sinθ等速飞行n y2=cosθ不做上下拐弯n z2=0不做左右拐弯曲率半径与过载之间的关系:ρy2=V2g(n y2−cosθ)竖直转弯曲率半径与过载之间的关系ρz2=V2cosθg(n z2)水平转弯曲率半径与过载之间的关系n L=1G(PsinαL+qSC ymax)极限过载表达式n L>n P>n R(LIMIT>P ASSABLE>REQUIRE)ε1=α−α∗=0 给定攻角下的理想控制关系式ε1=n y2−n y2∗=0 给定法向过载下的理想控制关系式α=n y2−(n y2b )α=0n y2αb 给定过载下小攻角的表达式式ε1=θ−θ∗=0 给定弹道倾角下的理想控制关系式ε1=ϑ−ϑ∗=0 给俯仰角下的理想控制关系式δz =K ϑ(ϑ−ϑ∗) 给定俯仰角下升降舵的偏转控制律θ=arcsin (1VdH ∗dt ) 给定弹道倾角的方案飞行可按给定高度飞行的方案弹道 α=mg P+Y α←[Psinα+Y =mg] 等高飞行下小攻角的表达式δz =−m z0+mgm zαP+Y αm z δz 等高飞行小攻角瞬时平衡假设下舵偏角表达式δz =δz0+K H (H −H 0)+K H ΔH等高飞行下升降舵的偏转控制律(微分项消除震荡) 侧滑转弯飞行情况下的飞行方案:3303()=y y b y b n n n ααα=- 平衡状态下的攻角的法向过载表达式303()1=y b y b n n ααα=- 平衡状态下无倾斜的攻角的法向过载表达式3031/cos ()=y v b y b n n αααγ=- 平衡状态下无侧滑的攻角的法向过载表达式水平面内给定弹道偏角下侧滑转弯飞行情况下的飞行方案: 2*0v v 给定弹道偏角的理想控制关系式dV dt =P−X m 切向方程303()1=y b y b n n ααα=- 竖直法向方程 −V gdψv dt n z3 b β=β 水平法向方程 dx dt=Vcosψv x 轴方向方程*()V V t 给定弹道倾角水平面内给定侧滑角或偏航角下侧滑转弯飞行情况下的飞行方案: φ:2*0v v 给定弹道偏角的理想控制关系式β:2*0v v 给定侧滑角的理想关系式dV dt =P−X m 切向方程303()1=y b y b n n ααα=- 竖直法向方程 dψv dt=1mV (Pβ−Z) 水平法向方程 dx dt=Vcosψv x 轴方向方程 dz dt =−Vsinψv z 轴方向方程φ:*()t 给定偏航角v =-水平飞行下侧滑、偏航、弹道偏角之间的几何关系 β:()*=t 给定侧滑角水平面内给定侧向过载下侧滑转弯飞行情况下的飞行方案:222*=n n ()0x x t 给定过载下的控制方程dV dt =P−X m 切向方程303()1=y b y b n n ααα=- 竖直法向方程 dψv dt=−g V n z2 水平法向方程dz dt =−Vsinψv z 轴方向方程 22b z z n n β角度和过载间关系 22*()z z n n t 给法向过载自动瞄准的相对运动方程组(极坐标系): cos cos T T drV V dt导弹与目标之间的矢径方向关系式 sin sin T T dq rV V dt 导弹与目标之间的角度方向关系式 q 导弹自身角度关系式q T T 目标角度关系式=0 导引关系式遥控导引的运动学方程组:d cos RV dt基站与导弹之间矢径方向关系式 sindR V dt 速度垂直于目标线方向上的关系式 航天器的开普勒轨道推导:3r r r 万有引力下的动力学方程 const h r r单位质量的角动量守恒 r rv h L 拉普拉斯常量-守恒 22v E const r 能量守恒 222=+2L Eh 三个守恒量之间的关系。

飞行器自动控制导论_第四章 纵向运动

飞行器自动控制导论_第四章 纵向运动

第四章 纵向运动 4.1 纵向运动线性化方程前面推导出来的线性化的纵向方程组重写如下:⎪⎪⎪⎩⎪⎪⎪⎨⎧∆+∆=∆-+∆+-∆-∆+∆=∆-+-∆--+∆-∆+∆=∆+∆-∆-T e q w w u T e qw w u T e w u T e T e T e M M dt d M dt d w M dt d M u M Z Z g dt d Z u w Z dt d Z u Z X X g w X u X dt dδδθδδθθδδθθδδδδδδ)()(]sin )[())1[()cos ()(22000 (4.1-1)其中e δ∆和T δ∆分别是空气动力控制项和推力控制项。

在工程实践中,力的导数q Z 和w Z 通常被忽略,因为它们对飞机响应的贡献非常小。

考虑到q ∆=∆θ,上面方程改写为状态空间的形式,得 ⎥⎦⎤⎢⎣⎡∆∆⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡+++⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡∆∆∆∆⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡+++-=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡∆∆∆∆T e w w w q ww w u w u wu w u T T ee Te Te Z M M Z M M Z Z X X q w u u M M Z M M Z M M u Z Z g X X q w uδδθθδδδδδδδδ0001000000(4.1-2)如果写成η B x A x+= 则有u w x q θ∆⎡⎤⎢⎥∆⎢⎥=⎢⎥∆⎢⎥∆⎣⎦,⎥⎦⎤⎢⎣⎡∆∆=T e δδη (4.1-3)⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡+++-=01000000u M M Z M M Z M M u Z Z g X X A w q ww w u w u wu w u (4.1-4) ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡++=00T TTeTe Z M M Z M M Z Z X X B w w δδδδδδδδ (4.1-5)矩阵中力和力矩的导数已经分别除以飞机的质量和惯性矩。

1 北航飞行力学_飞机性能计算的原始数据和质心运动方程

1 北航飞行力学_飞机性能计算的原始数据和质心运动方程
法向运动
xh
G d V Y Pky sin( P ) G cos g dt
北航 509
0

G
§1-3 飞机质心运动方程
几种特殊形式
•直线飞行(直线上升、下降等)
const , d / dt 0
•水平直线飞行(平飞加减速等)
G dV Pky Q G sin g dt Y G cos G dV Pky Q g dt Y G
喷气式发动机性能参数以及其高度特性、速度 特性、转速特性、特定油门状态 能画出铅垂平面内质心运动受力图,并推出各 种特殊运动状态下的质心运动方程
北航 509
2)最大状态:对应于最大许用转速(nmax)的发动机状态 。推力为非加力时的最 大值。只能连续工作5-10min,通常用于起飞、短时加速、爬升、空中机动等。 3)额定状态:对应于最大转速97% ,推力为最大状态的85-90%,可较长时间 工作(半小时~1小时),用于平飞、爬升、远航飞行等。
4)巡航状态:n巡90% n额,Pf巡 80%Pf额,耗油率最小,不限时,用于巡航。
最大可配平升力
Ymax
Y' LT ( ) max xA
Y2max Y1max
超音速时平尾平衡能力剧降形成飞行限制
C ymax
Cy
C ymax
最大允许升力系数
C ysx C yyx
C ydd C y max
C yyx min{ C ysx , C y max }
M
北航 509
f 0, 0 f 0, 0 一 般 f 0, 0 ( 0 f 0 0 0
0
Y 0
f 0, 0 0

飞机运动方程及小扰动方程推导

飞机运动方程及小扰动方程推导

飞机运动方程及小扰动方程推导飞机是一种重要的交通工具,它的运动可以通过一系列的方程来描述。

其中,飞机的运动方程是指描述飞机在空气中运动的基本方程,而小扰动方程是指描述飞机在受到微小扰动时的动力学方程。

飞机的运动方程是通过对牛顿运动定律的应用得到的。

根据牛顿第二定律,飞机的运动可以用以下方程来描述:F = ma在这个方程中,F代表飞机所受到的合力,m代表飞机的质量,a 代表飞机的加速度。

根据牛顿第二定律,飞机所受到的合力等于质量乘以加速度。

飞机所受到的合力可以分解为重力和气动力两部分。

重力是指地球对飞机的引力,可以用以下方程来描述:Fg = mg在这个方程中,Fg代表重力,m代表飞机的质量,g代表重力加速度。

重力的大小与飞机的质量成正比。

气动力是指飞机在空气中运动时所受到的阻力和升力。

阻力是指空气对飞机运动的阻碍力,可以用以下方程来描述:Fr = 0.5 * ρ * V^2 * Cd * A在这个方程中,Fr代表阻力,ρ代表空气密度,V代表飞机的速度,Cd代表阻力系数,A代表飞机的参考面积。

阻力的大小与空气密度、速度、阻力系数和参考面积有关。

升力是指空气对飞机垂直方向上的支持力,可以用以下方程来描述:Fl = 0.5 * ρ * V^2 * Cl * A在这个方程中,Fl代表升力,ρ代表空气密度,V代表飞机的速度,Cl代表升力系数,A代表飞机的参考面积。

升力的大小与空气密度、速度、升力系数和参考面积有关。

小扰动方程是指在飞机受到微小扰动时,飞机的动力学方程。

小扰动方程可以用以下方程来描述:m * δa = δF在这个方程中,m代表飞机的质量,δa代表飞机的微小加速度,δF代表飞机所受到的微小合力。

小扰动方程描述了飞机在受到微小扰动时的运动情况。

通过对飞机的运动方程和小扰动方程的推导,我们可以更好地理解飞机在空中的运动规律。

这些方程为飞机的设计、控制和性能评估提供了重要的理论基础。

同时,这些方程也为飞机的飞行安全和效率提供了指导,使得飞机能够更加稳定、安全地在空中飞行。

01_飞机的一般运动方程

01_飞机的一般运动方程

0 1 L qh 0 cos s 0 sin s
0 sin s cos s
coscos Ltq sin cos sin
sin cos 0
cos sin sin sin cos
2015/10/7 5
无人驾驶飞机:无人飞机和微型无人飞机
最大尺寸微型飞行器
英国的“Sender”无人机
微型飞行器和小尺寸无人机的尺寸对比
2015/10/7 6
“黑寡妇”微型飞机
“微星”微型飞机
2015/10/7
7
特殊航空器:微型扑翼和旋翼飞机
加州理工大学的“微型蝙蝠” 微型扑翼飞机
美国加州大学:扑翼机(翼展 200mm,总重11.5克,微型电 机驱动
10
三、飞机的主要组成部分及其功能
2015/10/7
11
机翼 :产生升力 ,机翼上一般有用于横向操 纵的副翼和扰流片;机翼前后缘部分还设有各 种形式的襟翼,增加升力 尾翼:水平尾翼和垂直尾翼;V型尾翼;水平尾 翼一般有水平安定面和升降舵组成;垂直尾翼 一般有垂直安定面和方向舵组成;超音速飞行 时通常采用全动水平尾翼(差动);鸭翼 机身:容纳人员、货物或其他载重和设备;要 求流线;飞翼式飞机取消机身。 起落架:起飞降落(机轮、滑撬、浮桶)
2015/10/7
17
半机体坐标系Oxbybzb :O在质心, Oxb沿飞 行速度矢量 V 在飞机对称平面投影方向, Oyb在对称平面内,垂直于Oxb向上(因而与 Oyq重合),Ozb垂直于飞机对称平面(与轴 Ozt重合)。
2015/10/7
18
2015/10/7
图2-2
19
航迹坐标系 Oxhyhzh : O 在质心, Oxh与 Oxq一 致,Oyh在包含飞行速度矢量V的铅垂面内, 指向上, Ozh 垂直于 Oxhyh(因而使水平的), 指向右。
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
共84页 飞行器飞行力学 21
4-2 飞机绕质心转动动力学方程
动量矩定理:
dh = ΣM dt
微元动量矩
Δh = r ×Vdm
微元绝对速度
V = VO + ω × r
全机动量矩 对质心
共84页
h = ∫ r ×Vdm =∫ rdm×VO + ∫ r × (ω × r )dm h = ∫ r × (ω × r )dm
可得标量方程组为:
共84页 飞行器飞行力学 9
质心动力学方程
dV x ⎫ m( + Vzω y − V yω z ) = Fx ⎪ dt ⎪ dV y ⎪ m( + V xω z − Vzω x ) = F y ⎬ dt ⎪ dVz ⎪ m( + V yω x − V xω y ) = Fz ⎪ dt ⎭
共84页
飞行器飞行力学
3
引言
研究内容 性能 飞行品质 质点动力学系统 质点系动力学系统
稳定性与操纵性:研究飞机受到扰动后及在 外力、外力矩作用下运动参数的变化特性。
稳定性 指飞行器在受到外界瞬时扰动后,是否具有自动地恢 复到原来平衡状态的能力。 操纵性 指飞行器对驾驶员操纵或舵面指令输入的响应,即从 一种飞行状态过渡到另一种飞行状态的能力.
共84页
⎡W x ⎤ ⎡0⎤ ⎢ ⎥ ⎢0⎥ Wy ⎥ = m⎢ ⎥ ⎢ ⎢W z ⎥ ⎢ g⎥ ⎣ ⎦ ⎣ ⎦g
飞行器飞行力学 16
航迹 轴系 χ, γ 矢量形式方程组 地面 轴系
μ (无风时)
气流 轴系 α,β 机体 轴系
ψ,θ,φ
⎛ δ Vb ⎞ dV b m = m⎜ + ω b × Vb ⎟ dt ⎝ δt ⎠ = Tb + Lba Aa + LbgW g
飞行器飞行力学
15
外力分量
1、推力
⎡Tx ⎤ ⎡ T cos ϕ ⎤ ⎢ ⎥ ⎢ ⎥ Ty ⎥ = ⎢ 0 ⎢ ⎥ ⎢Tz ⎥ ⎣ ⎦ ⎣ ⎦ b ⎢ − T sin ϕ ⎥
3、重力
2、气动力
⎡ Ax ⎤ ⎡− D⎤ ⎢ ⎥ ⎢ C ⎥ Ay ⎥ = ⎢ ⎢ ⎥ ⎢ Az ⎥ ⎢ − L⎥ ⎦ ⎣ ⎦a ⎣
转动惯 量矩阵
矩阵 形式
⎡ hx ⎤ ⎡ I x ⎢ ⎥ ⎢ ⎢ hy ⎥ = ⎢ − I xy ⎢ hz ⎥ ⎢ − I zx ⎣ ⎦ ⎣
− I xy Iy − I yz
− I zx ⎤ ⎡ω x ⎤ ⎥⎢ ⎥ − I yz ⎥ ⎢ω y ⎥ I z ⎥ ⎢ω z ⎥ ⎦⎣ ⎦
24
共84页
飞行器飞行力学
共84页
飞行器飞行力学
18
四、在气流坐标系中的平动动力学方程
根据速度之间的关系
u = V cos α cos β ⎫ ⎪ v = V sin β ⎬ w = V sin α cos β ⎪ ⎭
可得
du dV dα dβ ⎫ V sin α cos β − V cos α sin β ⎪ cos α cos β − = dt dt dt dt ⎪ dv dV dβ ⎪ V cos β sin β + = ⎬ dt dt dt ⎪ dw dV dα dβ ⎪ sin α cos β + V cos α cos β − V sin α sin β ⎪ = dt dt dt dt ⎭
− I zx ⎤ ⎡ω x ⎤ ⎡ M x ⎤ ⎥⎢ ⎥ ⎢ ⎥ − I yz ⎥ ⎢ω y ⎥ = ⎢ M y ⎥ I z ⎥ ⎢ω z ⎥ ⎢ M z ⎥ ⎦ ⎦⎣ ⎦ ⎣
25
飞行器飞行力学

dω x ⎫ 2 2 + ( I z − I y )ω y ω z + I yz (ω z − ω y ) + ⎪ Ix dt ⎪ dω y dω z ⎪ I xy (ω x ω z − ) − I zx (ω x ω y + ) = Mx⎪ dt dt ⎪ dω y ⎪ 2 2 + ( I x − I z )ω x ω z + I zx (ω x − ω z ) + ⎪ Iy ⎪ dt ⎬ dω z dω x I yz (ω x ω y − ) − I xy (ω y ω z + ) = M y⎪ ⎪ dt dt ⎪ dω z 2 2 + ( I y − I x )ω x ω y + I xy (ω y − ω x ) + ⎪ Iz ⎪ dt ⎪ dω y dω x I zx (ω y ω z − ) − I yz (ω z ω x + ) = Mz ⎪ ⎪ dt dt ⎭
共84页
飞行器飞行力学
6
4-1飞机质心动力学方程
基本原理: 牛顿第二运动定律
dV m =F dt
一、任意动坐标系中质心运动方程
速度和角速度在动坐标系的投影
V = V x i + V y j + Vz k
ω = ω xi + ω y j + ωzk
共84页 飞行器飞行力学 7
速度的微分
dV y dVz dV dV x j+ k i+ = dt dt dt dt di dj dk + Vx + V y + Vz dt dt dt
共84页 飞行器飞行力学 11
外力分量
1、推力
⎡Tx ⎤ ⎡ T cos ϕ ⎤ ⎢ ⎥ ⎢ ⎥ Ty ⎥ = ⎢ 0 ⎢ ⎥ ⎢Tz ⎥ ⎣ ⎦ ⎣ ⎦ b ⎢ − T sin ϕ ⎥
3、重力
2、气动力
⎡ Ax ⎤ ⎡− D⎤ ⎢ ⎥ ⎢ C ⎥ Ay ⎥ = ⎢ ⎢ ⎥ ⎢ Az ⎥ ⎢ − L⎥ ⎦ ⎣ ⎦a ⎣
共84页
飞行器飞行力学
14
三、在机体坐标系中的平动动力学方程
速度分量
⎡V x ⎤ ⎡u⎤ ⎢ ⎥ ⎢v ⎥ ⎢V y ⎥ = ⎢ ⎥ ⎢Vz ⎥ ⎣ ⎦ ⎣ ⎦ b ⎢w⎥
角速度分量
⎡ω x ⎤ ⎡ p⎤ ⎢ ⎥ ⎢q⎥ ⎢ω y ⎥ = ⎢ ⎥ ⎢ω z ⎥ ⎣ ⎦ ⎣ ⎦b ⎢ r ⎥
共84页
共84页
飞行器飞行力学
5
引言(续)
本章主要工作 (1) 应用牛顿第二定律、动量矩定律和运动学原理, 导出飞行器相对于动坐标轴系的一般运动方程组; (2) 在小扰动假设前提下,对非线性运动方程组进行 线性化; (3) 在特定条件下,将纵向运动和横航向向小扰动运 动分开,得到飞行器的纵向和横航运动方程组。
共84页
⎡W x ⎤ ⎡0⎤ ⎢ ⎥ ⎢0⎥ Wy ⎥ = m⎢ ⎥ ⎢ ⎢W z ⎥ ⎢ g⎥ ⎣ ⎦ ⎣ ⎦g
飞行器飞行力学 12
航迹 轴系 χ, γ 矢量形式方程组 地面 轴系
μ (无风时)
气流 轴系 α,β 机体 轴系
ψ,θ,φ
⎛ δ Vk ⎞ dV k m = m⎜ + ω k × Vk ⎟ dt ⎝ δt ⎠ = Lkb Tb + Lka Aa + LkgW g
飞行器飞行力学 22
因为
h = ∫ r × (ω × r )dm = ω ∫ r dm − ∫ r (ω ⋅ r )dm
2

ω = ω x i + ω y j + ωz k
r = xi + yj + zk

r = x + y +z ω ⋅ r = ω x x + ω y y + ωz z
2 2 2 2
共84页 飞行器飞行力学 19
解出
dV / dt , dα / dt , dβ / dt
共84页
飞行器飞行力学
20
在α、β 不大的快速机动中,可近似认为速度不 变,且
u ≈ V = const ⎫ ⎪ v ≈ Vβ ⎬ ⎪ w ≈ Vα ⎭
从而有
T sin ϕ − Z b g dα ⎫ = q − pβ − + cos φ cos θ ⎪ ⎪ dt mV V ⎬ Yb dβ g ⎪ = pα − r + + sin φ cos θ ⎪ dt mV V ⎭
共84页
飞行器飞行力学
10
二、在航迹坐标系中的平动动力学方程
速度分量
⎡V x ⎤ ⎡V ⎤ ⎢ ⎥ Vy ⎥ = ⎢ 0 ⎥ ⎢ ⎢ ⎥ ⎢Vz ⎥ ⎣ ⎦ ⎣ ⎦k ⎢ 0 ⎥
角速度分量
ω = ψ a + θa
⎡ω x ⎤ ⎡ 0 ⎤ ⎡ 0 ⎤ ⎡ − ψ a sin θ a ⎤ ⎡ − χ sin γ ⎤ ⎥ ⎢ ⎥ ⎢ 0 ⎥ + ⎢θ ⎥ = ⎢ ⎥= ⎢ γ θa ⎥ ⎢ω y ⎥ = Lkg ⎢ ⎥ ⎢ a ⎥ ⎢ ⎥ ⎢ ⎢ω z ⎥ ⎢ψ a ⎥ ⎢ 0 ⎥ ⎢ ψ a cos θ a ⎥ ⎢ χ cos γ ⎥ ⎦ ⎣ ⎦ ⎣ ⎦ ⎣ ⎦ ⎣ ⎣ ⎦k
共84页 飞行器飞行力学 17
标量形式方程组
du ⎫ m ( + qw − rv ) = T cos ϕ − D cos α cos β ⎪ dt ⎪ − C cos α sin β + L sin α − mg sin θ ⎪ ⎪ dv ⎪ m ( + ru − pw ) = − D sin β + C cos β + mg sin φ cos θ ⎬ dt ⎪ dw ⎪ m( + pv − qu) = −T sin ϕ − D sin α cos β ⎪ dt ⎪ − C sin α sin β − L cos α + mg co运动方程
内容
引言 动力学方程 平动动力学方程 转动动力学方程 运动学方程 平动运动学方程 转动运动学方程 运动方程线性化 线性化方程的分组 ——纵向/横航向
相关文档
最新文档