赵军民, 聂聪, 常冠男, 吕梅柏, 李新国
(1.西北工业大学 航天学院, 陕西 西安 710072; 2.西安现代控制技术研究所, 陕西 西安 710065)
全捷联式导引头直接固连于弹体,取消了机械常平架及其控制机构,降低了系统结构复杂度,减小了尺寸及质量,提高了可靠性;消除了俯仰、方位通道之间的交叉耦合;跟踪视线角速率也不再受限制。其诸多优点成为未来空地导弹小型化、智能化、低成本化的发展趋势[1-2]。
全捷联导引头虽然带来了诸多的优点,但也引入了一些技术难题[3]:无法直接测量视线角速度;瞬时大视场引入了大的测量噪声,降低了探测信噪比,进而影响制导精度;在中/末制导交接班时,受瞬时视场的限制,导弹攻角不能大于导引头的半瞬时视场角,对交接班截获目标提出了严峻的考验;过载受限,否则导引头容易脱锁甚至丢失目标。此外,对于侵彻型空地导弹,还需兼顾制导精度与落角。上述技术难点制约了全捷联式导引头在空地导弹上的工程化应用。
针对全捷联导引头的技术特点,一些学者展开了研究。LI等[4]提出了一种带有落角和视场角约束的自适应最优制导律,通过引入IAW因子调节制导律以满足相关约束;黄诘等[5]提出了一种偏置比例导引的间接撞击角度控制方法,能够适应视场角和可用过载等约束条件;Yan等[6]研究了一种考虑自动驾驶仪延迟的滑模制导律,对于目标机动具有一定的鲁棒性,能够补偿自动驾驶仪的延迟。王佩等[7]提出了将跟踪-微分器应用于全捷联制导角速度信息的提取的方法,张韬等[8]提出了一种基于USKF的全捷联导引头信息提取算法,李富贵等[9]提出了一种空地导弹打击移动目标的方法,但该方法并未充分考虑导引头视场角、机动过载约束等问题。
我们注意到,上述研究成果大多是针对全捷联式导引头工程应用中的某一个问题进行研究,采用单一的制导律完成全弹道约束控制,并未考虑空地导弹攻击过程中不同阶段所需解决的主要问题。本文针对全捷联式空地导弹视场小、过载能力有限、落角控制对时敏机动目标打击的弹道制导律问题进行研究,提出了一种分阶段制导方法。该方法通过对弹道进行合理规划,确保中/末制导交接班时能截获目标,末制导段采取一种兼顾稳定跟踪与命中精度的制导方法,保证导引头始终锁定目标并以期望落角命中。该方法可解决全捷联式空地导弹攻击地面移动目标的工程化应用问题,对采用全捷联导引头的精确制导武器弹道规划及制导律设计具有较强的工程化指导意义。
空地导弹离机后,在一定时间内保持舵面处于零位,以保证机弹分离安全;导弹稳定姿态,延时一定时间后,进入最大升阻比飞行阶段,目的是为了提高导弹可攻击距离;随后进入落角约束段,给截获目标创造有利条件,该段将导弹以期望落角导引至虚拟目标,当飞行至目标点上方时,导弹已具备一定的侵彻角度,视线偏差角也较小,初步具备截获目标的条件;在满足一定的弹目相对条件后,进入截获调整段,进一步修正视线偏差角,以确保目标进入导引头视场并趋于中心;当导弹飞行至导引头截获距离时,进入末制导段,导引头开机,采用模板匹配或特征识别方法截获目标,输出导引信号,引导空地导弹攻击目标,末制导段导弹应在有限过载内稳定跟踪目标;当飞行至导引头盲区时,保持停控前指令,直至命中目标。
综上所述,全捷联导引头主要受视场、可用过载、末端落角等条件约束,为解决其工程化应用问题,应将弹道进行合理规划。基于此考虑,将攻击弹道分为滑翔段、落角约束段、截获调整段、末制导段和盲区段。图1所示为分段制导示意图。
图1 分阶段制导示意图
图2所示为虚拟目标位置示意图。
图2 虚拟目标位置示意图
虚拟目标位置坐标表达式如(1)式所示。
(1)
式中:(xT,yT)T为真实目标坐标;(xVT,yVT)T为虚拟目标坐标;Rc为导引头截获距离;θd为期望落角。
设置上述虚拟目标点,当其运动至允许截获点时,虚拟目标点刚好与允许截获点重合,此时可同时满足落角、导引头视场角约束,顺利截获目标。
围绕带落角约束的制导问题,已提出了多种形式的制导律[10-13]。本文设计了一种快速非奇异终端滑模变结构制导律以满足落角约束。弹目相对运动关系示意图如图3所示。
图3 弹目相对运动关系
由图3可知,弹目运动学关系式为
(2)
选取状态变量
(3)
(4)
根据平行接近原理,设计制导律使得x1,x2在有限时间趋近于零,即能保证空地导弹以期望的落角精确命中目标。
对于制导系统(4),选取非奇异终端滑模面
s=x2+k1x1+k2|x1|γsgn(x1)
(5)
式中:k1=const>0;k2=const>0;1<γ=p/q<2,p,q均为大于零的正奇数。
设计自适应趋近律为
(6)
式中:α=const>0;η=const>0;ε=Δ+η=
const>0。
从(6)式可以看出,当R较大时,趋近速度较慢,当R→0时,趋近速度迅速增加,这样可以保证制导系统进入滑动模态。
对(5)式求导,并结合(6)式,设计制导律为
(7)
定理1对于制导系统(4)式、滑模变量(5)式、自适应趋近律(6)式、设计制导律(7)式,则系统状态将在有限时间内收敛至零,并且整个制导过程中不存在奇异的问题。
证明选取Lyapnov函数
(8)
对(8)式求导,有
(9)
因此,满足滑模到达条件。根据文献[14]的引理,系统是有限时间收敛的,有限收敛时间tr1满足
tr1 (10) 式中, s(0)=x2(0)+k1x1(0)+k2|x1(0)|γsgn(x1(0))。系统到达滑模面后,将沿滑模面运动直到收敛至原点,此时滑模面满足 s=x2+k1x1+k2|x1|γsgnx1=0 (11) (12) 选取Lyapnov函数 (13) 对(13)式求导,结合(12)式,有 (14) 同上,根据引理,系统是有限时间收敛的,有限收敛时间tr2满足 tr2≤ξ|x1(tr1)|(1-γ)/2 (15) 式中,ξ=2-(1+γ)/2(1-γ)k2。 由(10)式和(15)式可知,制导系统将在有限时间内收敛至滑动模态,然后沿着滑模面运动,当x1收敛至零时,x2也收敛至零。证毕。 1) 滑翔段 该段的作用是增加射程,采用最大升阻比方案。 α(t)=α*(Ma,h) (16) 式中,α*为当前飞行环境的马赫数和高度下最大升阻比所对应的攻角。 2) 落角约束段 y (17) 式中,hys为转入落角约束段的高度,其取值根据大量仿真结果确定。 3) 截获调整段 RMVT (18) 式中:RMVT为空地导弹与虚拟目标的距离;Rtz为调整段转换条件参数。 4) 末制导阶段 该阶段采用比例导引与姿态追踪的复合制导策略。一方面,比例导引能够提高打击精度,但其较大的前置角和攻角容易造成导引头丢失目标,而姿态追踪能够减小体视线偏差。因此,从兼顾精度和视场约束的角度考虑,采用比例导引与姿态追踪的复合制导律。其表达式如(19)式所示。 (19) 图4 过载指令选择示意图 5) 盲区段 当R≤Rcmin时,制导系统停控,该阶段制导指令为进入盲区时刻的制导指令,直至命中目标。 中末制导由于2种制导律给出的过载指令差异,可能引起姿态抖动,对导引头截获目标造成影响。因此,采用线性过渡方法实现平稳过渡。中末制导交接班过载指令表达式如下 aM=Kjatz+(1-Kj)amzd (20) 式中:atz截获调整段过载指令;amzd为末制导段过载指令,切换系数Kj的表达式如下: (21) 式中:Tgd=0.5 s为过渡时间;t0为末制导开始时刻。 全捷联导引头无法直接测量视线角,需对视线角进行重构,再经适当的微分网络,估计出视线角速度。本文采用跟踪微分器作为视线角的微分网络,图5为视线角速度估计方案结构图。 图5 全捷联导引头视线角速度估计方案结构图 视线角重构是利用视线坐标系、体视线坐标系、弹体坐标系、惯性坐标系之间的旋转关系,建立体视线角qB至惯性视线角qI的关系式。文献[7,15]对重构方法有详细描述,本文不再敷述,直接给出视线角重构表达式 qI=arcsinξ2 (22) 若对(22)式所重构的视线角直接进行微分,会极大地放大噪声。为解决非线性滤波问题,韩京清等[16]创造性地提出了非线性跟踪-微分器,该方法广泛应用于非线性滤波、非线性控制等领域中。根据全捷联导引头视线角重构模型及测量噪声的特性,本文采用一种非线性跟踪微分器对视线角速度进行提取。其形式如下: (23) 式中 r,α决定了跟踪的快速性;β影响噪声抑制效果。通过调整r,α,β3个参数,结合(22)式就可估计视线角速度。本文选取r=50,α=0.7,δ=5h,β=1.7,h为积分步长。 以某型全捷联空地导弹攻击地面目标为例,采用上述分段制导方案、视线角速度估计方法进行仿真以验证方案的有效性,仿真条件如下: 初始条件: 空地导弹初始位置(0,10 000)Tm,可用过载为2g,速度VM0=240 m/s,倾角θM0=0°,导引头截获距离Rc=3 000 m,导引头最小作用距离Rcmin=200 m,导引头视场ψFOV=±3°,期望落角θd=-79°,目标初始位置为(10 000,0)Tm,目标速度为VT=20 m/s。 测量误差: GPS测量误差/m:Δx~(0,100),Δy~(0,152);导引头测角误差/(°):ΔqB~(0,(0.286 5)2);姿态角测角误差/(°):Δϑ~(0,0.52)。 制导律转换条件:hys=9 000 m,Rtz=400 m; 末制导过程中,2种制导策略切换可能导致落角约束不能满足要求,实际应用中在落角约束段应选择适当的落角指令,留有余量。本文仿真选取指令落角为-85°。 仿真图6、图7和图9中,曲线中的标志点“·”,按顺序依次为落角约束段、截获调整段、末制导段和盲区段起始点。仿真图中的局部放大图为从截获调整段起始至进入导引头盲区截止。 图6 弹目运动轨迹 图7 弹道倾角 图8 指令过载与实际过载 图9 视线角偏差qE 图估计值与真实值 从图6空地导弹运动轨迹可以看出,空地导弹在末制导段的弹道较为平直。图7所示为空地导弹弹道倾角变化曲线,从图中可以看出,经过落角约束段之后,实现了落角约束要求,并留有余量,至末制导结束倾角为在-82°,随后进入盲区段,至仿真结束倾角为-79.65°,终端满足落角约束要求。进入末制导阶段,由于导引头探测的视线角qB、弹体姿态角ϑ均存在测量误差,采用(22)式进行惯性视线重构,并经过(23)式的滤波,有效降低了噪声带来的影响,从图8的过载指令与实际响应可以看出,过载波动响应全程在2g内,满足过载约束。图9所示,视线角从仿真开始时的45°,至姿态调整段开始时刻的2.54°,已进入导引头视场范围,经过调整段的视线偏差修正之后,qE从2.54°减小至1.50°,确保了导弹到达允许截获点时,目标处于导引头半视场并接近视场中心;从图9可以看出,在末制导段,当目标即将脱离导引头视场线性区时,姿态追踪制导律迅速调整姿态,使qE往趋于零的方向运动,以确保整个末制导段目标处于导引头视场内。图10可以看出,经视线角重构,利用非线性跟踪微分器估计的视线角速度能够满足末制导的需要。 本文针对全捷联空地导弹在导引头视场、可用过载、落角等约束条件下对地攻击移动目标的问题,提出了一种适合于工程化应用的弹道规划方案,并根据各阶段主要任务有针对性的设计了相对应的制导律。文中设计了带约束的非奇异快速终端滑模制导律实现了导弹以期望落角快速接近目标,通过引入虚拟目标和在截获调整段的姿态调整,实现了中末交接班时导引头准确截获目标,并保证目标始终处于视场内,最终以期望落角精确命中目标。仿真结果表明,该方法具有较强鲁棒性,对于采用全捷联制导体制导弹的工程化应用问题,该方法具有一定的指导意义。1.3 分段制导律设计
1.4 中末制导交接班
2 全捷联导引头视线角速度估计
2.1 视线角重构模型
2.2 视线角速度提取
3 弹道仿真验证
4 结 论