基于多体动力学的旋翼模型与气弹稳定性

针对无轴承等先进旋翼结构特点,建立了基于柔性多体系统动力学方法的旋翼气弹分析模型.旋翼各构件建模相互独立,便于组成不同构型旋翼,适合于新构型旋翼气弹分析.分析模型集成了适合于旋翼气弹分析的隐式大变形桨叶模型,构件动力学方程增加了铰链动力学方程,统一约束方程形式,改进了角运动约束方程,避免连续旋转引起的奇点问题.应用分析模型计算无铰式旋翼和无轴承旋翼的气弹稳定性,分析结构参数对旋翼气弹稳定性的影响.分析结果表明模型能准确计算结构弹性变形的耦合及非线性,提高旋翼气弹稳定性分析精度.      

基于等效迎角的气动融合轨道直接积分计算无控航天器轨道衰降研究

针对低地球轨道航天器服役期满离轨陨落、特别是功能失效处于无控飞行的航天器轨道衰降难以模拟的问题,在航天器轨道动力学方程直接积分计算框架中,把稀薄气体动力学与轨道动力学结合起来,发展了基于跨流域空气动力精细数值模拟驱动的低轨航天器气动特性一体化快速算法,提出了无控航天器轨道衰降失稳自旋飞行姿态等效迎角模拟法,初步建立了基于等效迎角的气动融合轨道直接积分高精度计算模型。对大型航天器A受控陨落开展计算分析,验证了本文发展提出的轨道计算、气动力一体化计算方法正确性与高精度。在大型航天器B无控失稳自旋飞行轨道衰降计算分析中,得到了与外测轨道星历数据吻合一致的长弧段和短弧段预报结果,在停止外测轨道数据供给5 h内,地心惯性系位置预报偏差低于1.5 km,验证了基于等效迎角的气动融合轨道直接积分计算模型对无控航天器轨道衰降预报的合理性。     

大迎角非线性飞行动力学（一）

本文较为系统地介绍了大迎角非线性飞行动力学研究的基础知识——非线性系统特性,大迎角时空气动力特性,大迎角飞行时特殊现象,非线性系统的近似研究方法及其在非线性飞行动力学中的应用.     

电控旋翼气弹动力学建模研究

研究了电控旋翼的气弹动力学建模方法.以有限元旋翼气弹分析程序LORA01为基础,首先给出了一种能计及襟翼移轴补偿影响的时域综合非定常气动力模型,并将伺服襟翼作为额外引入的质量体计入结构动力学模型;之后结合有限元法,根据Hamilton原理导出由广义力表示的带襟翼桨叶非线性运动方程,采用基于Newmark方法的隐式数值积分法对桨叶运动方程进行求解.最后,利用主动控制襟翼旋翼的试验数据对模型进行了比较验证,证明了所建模型的正确性.     

电控旋翼直升机飞行动力学建模与配平特性研究

首先基于Peters He广义动态尾迹理论,建立了电控旋翼动态尾迹入流模型,进一步结合电控旋翼襟翼操纵与桨叶变距之间的关系、桨叶挥舞运动方程和带襟翼翼型非定常气动力模型建立了适用于飞行力学分析的电控旋翼气动力模型.在此基础上,结合机身、尾桨、尾面的气动力模型,建立了完整的电控旋翼直升机飞行动力学分析模型.以Z-11直升机为基准改造为电控旋翼直升机作为算例,计算了前飞状态下电控旋翼直升机的诱导速度分布和桨盘迎角分布,对比了电控旋翼与常规旋翼的气动特性差异;在此基础上,进一步分析了电控旋翼直升机的配平特性随前飞速度的变化规律以及与常规直升机的差异.     

基于嵌套网格的内埋物单侧投放影响研究

内埋物投放分离技术是新一代飞行器研制的关键技术之一。首先，基于动态嵌套非结构网格耦合求解流场非定常N-S方程和刚体动力学方程，实现了内埋物投放的精确模拟；然后，研究载机不同飞行姿态对内埋舱内左侧无挂载、仅右侧挂弹投放过程的影响。研究结果表明：正迎角工况对弹体投放分离影响较小，负迎角工况将引起弹体有较大的反方向的横向运动和偏航；负侧滑工况对弹体投放影响较小，正的小侧滑角工况下，在投放分离过程中弹体气动力和气动力矩存在较强的脉动变化；随着侧滑角进一步增大，弹体气动特性的脉动变化减小。     

大展弦比飞机运动稳定性建模及分析

随着飞机设计技术的提高和新型材料的应用，现代飞机气动弹性效应越来越显著，刚体运动和弹性体振动互相解耦的关系不再明显，在分析弹性飞机相关问题时，应该综合地、统一地考虑飞行动力学与气动弹性力学问题。本文在一般体轴系下综合考虑了飞机刚体运动自由度和弹性振动自由度，建立弹性飞机刚弹耦合运动的方程，通过在配平状态下引入小扰动运动假设，以及利用有理函数拟合技术将偶极子格网法计算所得到的频域非定常气动力转化为时域形式，建立了大柔性飞机刚弹耦合小扰动状态空间方程，并对纵向运动稳定性进行分析计算。针对文中构型飞机，运用本文所采用的方法计算所得结果与传统的线性及刚体计算所得结果进行对比分析，非线性计算方法所得颤振速度更小，而且失稳模态亦发生变化，对于飞行力学模态而言，长周期模态稳定性变差，在计算速度范围内出现失稳的现象，由此说明机翼大变形对飞机稳定性所带来的影响。     

周期平均控制下固定翼双旋弹角运动特性分析

针对应用在固定翼双旋弹上的一种周期平均控制方法，通过对弹体修正组件滚转角进行傅里叶变换，得到了双旋弹在周期平均控制下角运动方程的模型，进而分析了周期平均控制下固定翼双旋弹的角运动响应，并根据双旋弹角运动特性和执行机构的性能限制给出了周期平均控制角频率的设计方法。引进2种具体的周期平均控制方案作为分析对象，并通过六自由度弹道仿真进行了验证。结果表明，对固定翼双旋弹角运动响应的分析能够反映周期平均控制的特性，周期平均控制角频率的设计方法可以设计出合理的角频率。     

倾转旋翼机的多体运动建模与仿真(英文)

现有的倾转旋翼机建模着重于研究其复杂的气动特性,在飞行动力学运动方程方面则研究不足,均采用了普通固定翼飞机的六自由度刚体动力学方程。然而,倾转旋翼机在过渡模式,旋翼和短舱的倾转使飞机构型产生动态变化,除了造成气动力干扰以外,旋翼和短舱的倾转还会引起重心变化,而且由于高速旋转的旋翼倾转会产生陀螺力矩效应,造成对机体的姿态干扰,采用六自由度刚体动力学方程不能反映由此造成的多体惯性耦合问题。为此,本文提出将旋翼、短舱、机体视为不同运动实体,建立倾转旋翼机的多体动力学方程,以完善其飞行动力学建模。建模过程中首先对各个微元运用牛顿定律和角动量定理,然后对各实体积分得到了惯性系下多体动力学方程,最后通过转换得到机体坐标系下的多体动力学方程组。该方程组在形式上为隐式非线性微分方程组,需首先求解状态微分量的相容初值。然后采用Runge-Kutta-Felhberg积分算法,对隐式非线性微分方程组进行了数值仿真。以XV-15数据为例对模型和仿真算法进行检验,结果表明模型和仿真算法正确。建立的模型和采用的仿真算法对研究倾转旋翼机飞行动力学特性、飞行控制技术、飞行安全问题有重要的理论和实际意义。     

重型直升机飞行动力学刚弹耦合建模及空中共振稳定性分析

针对重型直升机（HLH）大重量、低转速的固有特性,提出了一种适用于重型直升机的飞行动力学刚弹耦合建模方法。该方法结合传统直升机飞行动力学与旋翼机体耦合动力学,将传统飞行力学的分析频段拓展到了5 Hz,额外考虑了桨叶和机体的弹性变形,基于阻抗匹配法推导出了显式的旋翼/机体耦合动力学方程,模拟了真实飞行状态下的直升机气弹耦合特性,利用该模型计算并分析了算例重型直升机的悬停飞行特性和空中共振稳定性。结果表明：旋翼机体耦合导致摆振前进型和机体弹性模态的阻尼-转速曲线先相互靠近至同一点再分离,可能引起直升机的高频瞬态振动;在摆振等效阻尼不足时,旋翼摆振后退型是不稳定的,但随着等效阻尼增加,摆振二阶周期型模态和机体弹性模态会出现耦合;桨叶弹性变形与机体弯曲模态及挥舞集合型耦合,但不会引起明显的不稳定现象。     

高速飞行弹箭目标表面动态热辐射

针对红外系统实时跟踪捕获高速飞行弹箭目标的关键问题,提出了一种用于求解其全弹道动态红外辐射（IR）特性的方法.以典型155mm口径无控弹箭为研究对象,基于Simulink建立了模块化的6自由度(DOF)刚体弹道仿真模型,数值计算并分析了弹道诸元的变化规律.应用热网络法建立了弹体表面耦合换热动态热辐射场的物理模型,推导了弹箭高速旋转飞行的气动加热计算模型,并利用蒙特卡洛(M-C)法考虑了环境热辐射的影响.运用Runge-Kutta法耦合求解节点热平衡方程组,得到了动态气动表面传热系数、温度场以及红外辐射场的分布规律,对比分析了目标在整个飞行过程、不同部位、不同波段内的红外辐射特征.结果表明:目标发射后,其表面温度迅速升高,越靠近弹头部,温度升高速率越快,峰值温度越高；在飞行前20s内,其红外辐射特征明显；随着飞行速度衰减,热量散失较快,辐射强度较弱,且主要集中于8～14μm波段.      

基于拉格朗日力学的伞-弹系统动力学模型

为了回避牛顿力学体系在多体动力学建模过程中存在的复杂约束力求解问题,以拉格朗日力学体系为基础,从分析力学的角度出发,建立了铅垂平面内降落伞-弹系统五自由度动力学模型。详细介绍了拉格朗日力学体系动力学建模过程,推导了铅垂平面内伞-弹系统动能公式,解决了伞-弹系统广义力求解问题,给出了伞-弹所受空气动力的广义力表示。算例将建立的伞-弹系统模型和牛顿力学体系下简化的伞-弹模型进行对比,两种模型计算结果吻合,验证了本文模型的正确性。降落伞-弹系统动力学模型验证了拉格朗日力学在常规飞行器动力学建模中的可行性,为使用ADAMS等基于拉格朗日力学思想建模的商业软件进行伞-弹系统动力学求解提供理论依据,可以用于指导伞-弹系统的分析和设计。     

飞机结冰后非线性气动力建模及动态响应特性研究

当前结冰后的动力学仿真主要基于线性气动力模型及结冰影响模型进行,对于结冰后非线性气动力模型的研究较少,然而这对于飞机结冰后在大迎角等非线性区域的动力学分析极为重要.研究了飞机结冰后的非线性气动力模型,将结冰后线性气动力模型的估算方法扩展到非线性气动力模型上,并在结冰后大迎角区域风洞试验数据的基础之上,对失速及过失速区的气动参数变化进行修正.研究了飞机在不同升降舵脉冲信号作用下的动态响应.仿真结果表明,计算结果能够较好地反映飞机在结冰状态下的动力学特性,提出的非线性结冰影响模型能够为结冰条件下研究飞机失速及过失速区域的动力学特性提供理论支撑.     

高超声速飞行器非线性纵向运动稳定性研究

针对研究高超声速飞行器飞行动力学稳定性时忽略方程非线性影响的问题,建立了保留非线性特征的六自由度动力学方程组和非线性活塞理论耦合模型。研究了初始迎角和活塞理论非线性项对高超声速飞行器纵向稳定性的影响,并在稳定巡航中施加瞬时扰动,分析其姿态响应。仿真结果表明:高超声速飞行器姿态收敛速度比常规低速飞行器慢很多;对迎角施加扰动后,飞行器姿态变化主要影响短周期模态,且偏离平衡位置越远,收敛越快;改变气动力非线性项时,其姿态变化主要影响长周期模态,且收敛较慢;两种扰动均对振荡周期无明显影响。     

固定翼双旋弹动力学分岔特性分析

针对一种滚转稳定的固定翼双旋弹,对其非线性动力学进行了分岔特性分析,并在此基础上研究了各系统参数对其动力学分岔特性的影响。根据固定翼双旋弹非对称的特点,通过数值计算方法研究其飞行过程中平衡点随同向鸭翼安装角的变化规律,通过系统的分岔图得知系统具有三组稳定平衡点,其中只有一组平衡点为理想可行的稳定平衡点,因此需限定同向鸭翼安装角的范围以使固定翼双旋弹保持稳定飞行。在此基础上针对固定翼双旋弹弹道修正组件周期旋转和转角固定两种工作模式,通过各系统参数下的系统分岔图总结了固定翼双旋弹结构及气动力参数对其动力学系统分岔特性的影响。仿真结果表明,固定翼双旋弹的各气动力参数及飞行速度均对系统的分岔特性具有较大影响,应合理选定这些系统参数以使其具有良好的气动特性。     

基于闭环控制的带翼导弹虚拟飞行数值模拟

通过耦合求解非定常Navier-Stokes方程、刚体六自由度（6DOF）运动方程及舵偏控制律,发展了一套虚拟飞行数值模拟平台,对基本带翼导弹外形进行了虚拟飞行数值模拟研究。流场求解采用基于结构网格的有限体积法,时间推进采用双时间步法,使用刚性重叠动网格技术模拟舵与弹身的相对运动,采用四阶Runger-Kutta方法求解运动方程组,舵偏控制律采用基于比例-积分-微分（PID）的反馈控制。模拟结果表明：发展的虚拟飞行数值模拟平台能够处理包含复杂运动边界的非定常运动问题,可以为非定常气动力建模提供数据、预示飞行器稳定性及检验飞行控制律,该模拟平台具备一定的工程应用价值。     

鸭式翼弹道修正火箭弹修正能力分析

采用十字型鸭式翼二维弹道修正火箭弹弹体布局,系统地研究了修正火箭弹的修正弹道理论,建立了修正质点弹道的数学模型。在Matlab中编制了四自由度弹道模型,仿真分析可用于分析无控和有控情况下的弹道修正能力问题。大射角下舵机45s起控纵向修正效果最好,起控时间越早,横偏修正越好。舵偏角越大,修正能力越强。     

旋转弹体马格努斯效应数值模拟方法研究

旋转的弹体通常会存在马格努斯效应,严重影响飞行稳定性和弹道轨迹。提出基于旋转壁面法的多对称面弹体旋转效应数值模拟方法,无需借助动网格和多参考系模型,仅需要将旋转速度附加壁面,即可使用定常方法模拟马格努斯效应;同时,使用SOCBT标准弹体模型对该方法进行验证分析。结果表明:计算值与试验值吻合较好;由于旋转引起的附面层堆积以及大攻角分离效应、激波的干扰,计算时附面层的网格节点要足够密,并且应该针对不同的攻角选择精度与效率兼顾的计算模型。     

风雨对飞机飞行安全性的影响

风雨严重影响飞机的飞行安全。基于动量定理建立了一种根据降雨条件、风场特性、飞机特征和飞行状态计算雨滴对飞机产生的撞击力和力矩的方法,研究了风雨对飞机气动特性的影响,进而建立了飞机在风雨中飞行的运动方程。通过引入驾驶员的操纵模型,对飞机在风雨中飞行的运动特性进行了数值仿真研究。仿真结果表明:降雨会使飞机的升力减小,阻力增大,平衡迎角增加,不利于飞机的飞行安全;当降雨过程中伴随有风时,风会改变飞机的飞行迎角,使得飞机更容易出现因失速而造成的飞行事故;降雨时飞机的抗风能力降低。     

Impact point prediction guidance based on iterative process for dual-spin projectile with fixed canards

The fixed canards configuration of a dual-spin projectile makes it difficult to apply the traditional guidance law. In this study, a modified impact point prediction guidance strategy based on an iterative process was developed for a class of dual-spin projectiles with fixed canards, to reduce the impact point dispersion. The guidance strategy is dependent on the modified projectile linear theory to rapidly predict the flight states and the impact point. For projectiles with control applied to the trajectory, the modified projectile linear theory method is known to achieve poor impact point prediction. To improve the prediction accuracy, improvements were made to the modified projectile linear theory by considering the products of the yaw rate and other small quantities. The guidance strategy is based on the iterative process for the continuous adjustment of the expected output of the roll angle of the course correction fuze, to minimize the direction error between the predicted impact point and target location. Studies were conducted on a model dual-spin projectile configuration to demonstrate the guidance details. The numerical simulations indicate that the proposed guidance strategy can effectively reduce the projectile impact point dispersion.