基于可控翼伞回收技术的多体动力学仿真研究

基于可控翼伞回收技术的多体动力学仿真研究

近年来，随着中国空间技术的快速发展，太空发射活动日益增多。虽然地面推进器的传输范围避免了人口稠密地区，但地面推进器的残余和残留物仍然对一些村庄和城市居民的生活和财产损失有重大影响。这不仅增加了人员分配的成本，而且降低了宇宙运输能力，并造成了不利的政治影响。因此近年来,火箭及其组件的可重复使用技术已成为领域内的研究热点。翼伞不仅具有普通降落伞折叠包装方便、体积小和质量轻的优点,而且具有良好的滑翔性、操纵性和稳定性,能够使火箭助推器落点散布范围减小,在航空航天飞行器和分离体的精确定点回收方面有着很高的应用价值。火箭助推器翼伞回收系统,如图1所示,主要由翼伞(包括伞衣、伞绳和吊带等)、操纵控制系统和系统集成结构(控制平台)组成,火箭助推器作为其有效载荷,控制系统设备布置在控制平台或火箭助推器上,翼伞系统在火箭发射时包装在火箭助推器的头锥部位。目前国内外尚无采用翼伞系统对火箭助推器回收的先例。为实现火箭助推器的定点回收任务,要对翼伞系统的飞行轨迹进行规划和控制,前提是要对翼伞系统的动力学响应特性进行深入研究。然而,翼伞系统与传统航空飞行器不同,翼伞与有效载荷或控制平台之间存在明显相对运动,尤其是开伞、转弯机动或雀降机动过程中更加明显,因此要从多体系统动力学的角度对翼伞-载荷系统进行动力学建模和分析;另外,翼伞是一种超轻结构,在分析翼伞系统的运动特性时,需要考虑翼伞表观质量的影响目前用于翼伞系统多体动力学特性研究的仿真模型均将伞体和有效载荷分别看作两个独立的刚体,不同之处在于连接约束模型,归纳起来主要分为3类:第一类是“两体+弹簧”,两体(控制平台简化到火箭助推器固定结构上)间通过弹簧相连,连接点的空间位置显式求解,弹簧参数设置依赖于设计和经验,增加了求解自由度且难于收敛;以Vishnyak第二类是“两体+约束方程”,两体间通过约束方程建立,保证了连接点空间位置的一致性,建模过程程式化,便于模型的扩展,但增加了求解自由度;以wolf第三类采用“共铰点两体+扭簧”,两体间共用一个坐标系原点,相对姿态关系通过扭簧建立关系,其参数的设置依赖于实际系统和仿真经验,为最小解集系统。以Barrows本研究采用扩展性好的拉格朗日乘子法(第二类模型)建立系统动力学方程,对两种吊挂方式(如图1所示)和火箭助推器外形气动力对系统滑翔和转弯特性的影响进行研究。1翼伞系统动力学建模本文主要研究对象是火箭助推器翼伞回收系统在开伞和充气张满到着陆过程的系统动力学行为,不研究流场作用机理,因此对一般可控翼伞-控制平台-载荷系统建立了两体(控制平台简化到火箭助推器固定结构上)8DOF(自由度)和三体10DOF动力学仿真模型,采用Barrows表观质量估算方法1.1dof气动力在建立系统动力学仿真模型时,对翼伞系统做了如下基本假设:1)飞行水平距离短,大地是水平的,忽略地球自转,且重力加速度恒定;2)翼伞完成充气张满后几何形状不变,左右对称,后缘操纵只影响气动力;3)翼伞、控制平台和火箭助推器均看作6DOF刚体,两体间通过刚性吊带连接约束;4)控制平台和火箭助推器的质量特性恒定不变,但翼伞含有表观质量,且随大气密度变化;5)火箭助推器的气动力模型简化为圆柱体,忽略控制平台的气动力影响;6)忽略翼伞操纵过程的系统响应延迟。1.2坐标系定义1x-qp-4控制平台坐标系惯性坐标系:固定于地球表面的东-北-天坐标系,坐标原点在着陆目标点地面位置,x轴指向当地东向,y轴指向当地北向,z轴由右手定则确定;翼伞坐标系:固定翼伞上的本体坐标系,坐标原点位于翼伞的等效平均气动弦1/4处(c/4),x轴在翼伞纵向对称面内与平均气动弦线平行并指向开口方向,z轴在翼伞纵向对称面上与x轴正交并指向上表面,y轴由右手定则确定,如图1所示;控制平台坐标系:固定在有效载荷平台上的本体坐标系,坐标原点位于有效载荷平台的质心位置,x轴与稳定飞行前方一致,z轴与x轴垂直并指向上方,y轴由右手定则确定;有效载荷坐标系:固定在火箭助推器上的本体坐标系,坐标原点位于火箭助推器的质心位置,x轴指向尾喷口或头锥,z轴与x轴垂直并处于迎风区或吊点一侧,y轴由右手定则确定,如图1所示;翼伞气流坐标系:固定在翼伞上的本体坐标系,坐标原点位于翼伞等效c/4处,x轴指向气流的来流方向,z轴在翼伞的纵向对称面内且与x轴垂直,y轴由右手定则确定。2序、攻角和侧滑角从惯性坐标系到翼伞/控制平台/有效载荷坐标系采用z-y-x转序的欧拉姿态角(偏航、俯仰和滚转)描述,记为(ψ,θ,φ),通过下角标区分,转换矩阵为S从气流坐标系到翼伞坐标系采用z-y-x转序的攻角α和侧滑角β描述,转换矩阵为C1.3系统的约束模型翼伞系统中,翼伞、控制平台和火箭助推器均看作6DOF刚体,对于垂直吊挂的两体(控制平台简化到火箭助推器上)系统模型,系统完整约束的多体系统动力学方程形式为其中:T为系统动能;q为广义坐标;λ为拉格朗日乘子;Q为广义力。为求得方程(1)左侧的前两项,取惯性系下广义坐标和广义速度其中:q为12×1列阵,不能直接积分;下标P表示翼伞;C表示控制平台或有效载荷。从本体系到惯性系的坐标转换矩阵为S,有系统的动能T可表示为其中,I为惯量矩阵,H为惯量交叉矩阵。对于一般物体,质量矩阵M=mE=S翼伞系统多体间的约束模型是一个复合约束副,如图2所示。1)翼伞-控制平台的连接特征线A其中,ε=(-sinθ2)翼伞-有效载荷(控制平台)连接特征线A3)翼伞-有效载荷间沿A其中,r对于三体系统,控制平台与有效载荷间约束与式(8)和式(9)同,可直接写出结果:为得到方程(1)第3项,并直接求解出拉格朗日乘子,对约束方程求导其中,ue2e1为拉格朗日方程广义坐标,其中,对于三体系统,相应的3个余项与R两体系统的动力学方程可以表示为其中,S从而最终得到的系统动力学方程为以上方程可通过数值积分方法进行求解。1.4表内质量评价本研究所述翼伞沿展向展开后为矩形的圆弧形翼伞,采用Barrows给出的翼伞的表观惯量矩阵的估算方法其中,h1.5气动载荷变化翼伞气动力计算采用经典飞行力学方法,计算式及相关气动系数如表1所示。升阻力气动系数参考了lingard在考虑有效载荷气动力时,将火箭助推器简化为圆柱体扰流模型,不考虑两端圆形端面的气动力,圆柱面受到气动阻力和气动升力作用,如图2所示,采用经典飞行力学计算方法,主要与来流总攻角有关,阻力系数和升力系数如图3所示2动力学模拟研究2.1模拟条件表明翼型Clark-Y,弦长c=5.55m,翼展b=14.4m,厚度e=0.15,系统长度R=11.5m,安装角θ2.2载荷变化对翼伞液压系统的影响本节采用了两体8DOF仿真模型。图4给出了载荷气动力对滑翔运动影响的仿真分析结果。图5给出了载荷气动力对转弯运动影响的仿真分析结果,翼伞在40s时左侧施加60%的单边转弯操纵。分析得到以下基本结论:1)从图4(a)、图4(b)和图5(a)和图5(c)看出,载荷外形的气动力使系统水平速度微弱减小,而下降速度增大,从而滑翔性能降低。图4中滑翔比降低约10%左右;2)从图4(d)、图4(c)和图5(d)看出,载荷气动力会导致滑翔和转弯时的翼伞配平攻角微弱减小,对转弯时的侧滑角影响很小;3)从图4(e)、图4(f)和图5(f)、图5(g)看出,载荷气动力会导致翼伞和有效载荷俯仰角增大,而翼伞相对于有效载荷的俯仰角变化不大,相当于增大了翼伞的安装角;4)从图5(f)看出,在转弯运动过程中,载荷气动力会引起有效载荷俯仰角和偏航角振荡,而对滚转角影响很小。2.3有效载荷气动力本节采用了三体10DOF仿真模型。图6给出了有效载荷气动力对滑翔运动影响的仿真分析结果。图7给出了有效载荷气动力对转弯运动影响的仿真分析结果,在40s时左侧施加60%控制量。分析得到以下结论:1)从图6(a)、图6(b)、图7(a)~图7(c)可以看出,有效载荷外形的气动力对系统的水平速度、下降速度、转弯半径、滑翔比和转弯半径等影响很小;2)从图6(c)、图7(d)看出,有效载荷气动力使对翼伞配平攻角微弱减小,对侧滑角影响很小;3)从图6(d)~图6(f)和图7(e)~图7(g)看出,有效载荷气动力会引起其与控制平台的俯仰角微弱增大,发生小角度振荡,同时引起翼伞俯仰角小幅增大,而翼伞相对平台俯仰角变化很小;4)从图6(e)看出,在转弯运动过程中,载荷气动力会引起有效载荷、控制平台和翼伞发生相对运动,偏航振荡显著,对滚动角影响很小,对导航测量与控制会产生不利影响。3气动载荷的影响分析本文以火箭助推器翼伞回收完成充气张满后的多物体飞行动力学系统为研究对象,将翼伞简化为刚体模型,连接吊带简化为刚性连接杆,采用拉格朗日乘子法建立了两体8DOF和三体10DOF动力学仿真模型,并考虑了翼伞的表观质量特性。对垂直吊挂方式和水平吊挂方式中火箭助推器外形气动力对飞行性能进行了仿真分析。火箭助推器外形气动力对滑翔和转弯性能的影响分析表明,有效载荷火箭助推器外形会使翼伞和控制平台的滑翔攻角小幅增大,对水平