基于大变形梁模型的风力机叶片气动弹性分析

以水平轴风力机为研究对象分析其叶片大幅值气动弹性动态响应问题.基于非线性欧拉梁模型和风力机叶片动态入流、动态失速气动力模型进行气动弹性建模.气弹方程的动态平衡求解采用Newton-Raphson迭代,分别通过牛顿下山法和同伦沿拓法这两种数值手段改善其迭代收敛性,并对这两种算法的有效性进行比较,结果表明同伦沿拓法的收敛性相对较好.以两种典型的水平轴风机为算例,在额定工作状态以及某些极端风况下计算叶片动态响应.根据叶片根部弯矩载荷幅值定义系统是否达到临界稳定状态,并针对关停状态的风机研究了叶片临界风速与来流方向角的关系.结果表明,某些极端风况和工作状态下,风机叶片根部弯曲载荷较大,相应临界稳定的风速较低.     

基于解析梯度的经典Lambert问题迭代求解方法

针对现有求解模型复杂、收敛速度慢等问题,在将经典Lambert转移问题转化为超越方程的基础上,提出一种基于解析梯度的Lambert问题迭代求解算法。选择转移轨道的真近点角为迭代变量,导出转移时间关于真近点角的解析梯度,构造一种基于解析梯度的牛顿迭代算法,降低了算法计算复杂度。理论分析表明该算法具有二阶以上的收敛速度。依据偏心率向量与转移轨道形状的关系,通过几何方法分析得到转移轨道在初始位置处的速度约束条件,推导转移轨道真近点角的最大值和最小值的解析表达式,并采用线性插值方法确定迭代初值,进一步提高了迭代算法的收敛速度。数学仿真结果表明在各种转移条件下算法均能快速收敛,采用所给出的初值选取方法初值确定精度高,进而能够加快收敛速度,而与较割线法相比较收敛速度快、计算量小,验证了所提出算法的有效性。     

天基对地手打击动能武器再入解析预测制导技术

天基对地打击动能武器用于从太空对地面高价值战略目标进行快速、准确的打击.为了在各种干扰因素的作用下仍能保证足够的命中精度,动能弹头必须实施再入制导.针对实时性要求,探讨了一种解析预测制导方法.首先详细推导了零攻角再入弹道参数的三维解析解,在此基础上借鉴牛顿迭代法的思想设计了速度倾角与航向角的迭代修正算法,并最终将其用于制导指令的生成.仿真结果表明,解析预测制导方法能有效提高再入弹头的落点精度,且实时性强.此外,通过对不同制导参数下的制导性能进行分析还发现,制导步长取450～600米、精度参数取20～50米最为合适.     

飞船舱-伞系统稳降阶段动力学模型的建立

建立一种思路清晰、求解方便易行的降落伞返回舱系统的动力学模型是研究降落伞-返回舱系统动力学特性的基础。文章从一般刚体动力学方程入手,结合刚体的运动学方程,建立了稳降及着陆阶段舱-伞系统的动力学模型,并与空投试验的一些测试曲线进行对比,证明了本模型的正确性和有效性。     

飞机三维结冰模型及其数值求解方法

 对飞机结冰外形进行纯三维数值模拟,是目前飞机结冰预测中的一大难点。为了建立三维结冰数值模拟方法,基于Messinger的二维结冰模型,提出了一种考虑液态水溢流效应的三维结冰计算模型,并针对该模型建立了表面单元内溢流水流动的分配方案、发展了相应的迭代求解方法。采用本文方法对MS-317后掠翼结冰进行了计算,并与实验和Lewice3D的计算结果进行了对比。研究结果显示：本文方法具有较好的收敛性,计算的霜冰和明冰外形均与Lewice3D计算的冰形一致;对于霜冰和结冰时间较短的明冰,本文计算的冰形与实验吻合较好;对于结冰时间较长的明冰,本文计算的冰形与实验对比还有一定差异,但冰生长的总体趋势和大致体积与实验一致。同时,对比了MS-317无后掠翼与有后掠翼结明冰的外形,发现机翼后掠导致的三维溢流效应对结冰外形有明显影响,因此,对于三维结冰分析,如果用二维截面的结果来代替三维结果,其合理性还需验证。     

飞机结冰中水滴撞击特性的欧拉法准确性分析

过冷水滴撞击特性计算是飞机结冰冰形预测与防除冰系统性能分析的基础，常用方法为欧拉法与拉格朗日法，2种方法的结果通常是一致的，但在某些部件上会出现差异。通过对欧拉法与拉格朗日法进行比较，分析2种方法结果的异同，进而讨论欧拉法对于飞机结冰中水滴撞击特性计算的准确性。以NACA 0012翼型、冰风洞风道、S形进气道与发动机气膜帽罩为对象，采用欧拉法与拉格朗日法计算获得水滴运动及局部水收集系数。结果表明：当水滴运动未受到上游效应影响时，欧拉法与拉格朗日法的结果一致；当水滴发生偏转后，欧拉法速度的单一性使水滴流线不能相交，而拉格朗日法能捕捉水滴轨迹的交叉，导致2种方法的预测产生差异，且欧拉法结果与水滴不碰撞及聚并的假设存在冲突；水滴在上游部件空气绕流或气流吹袭作用下都会发生偏转，使得欧拉法与拉格朗日法得到的下游表面水收集系数不相符，欧拉法对于S形进气道与发动机气膜帽罩的水滴运动及撞击特性计算存在误差。研究成果对飞机结冰冰形的准确预测及防除冰系统的设计有重要参考价值。     

基于数据库的助推滑翔导弹助推段弹道诸元快速解算方法

针对助推滑翔导弹快速发射的任务需求，提出了一种基于诸元数据库的助推段弹道诸元快速解算方法。首先，利用滑翔段弹道解析估算，得到影响射程的助推段终端参数。然后，对助推段终端空域进行网格化剖分，针对每个高度速度特征点，利用牛顿迭代法对弹道诸元进行解算并形成诸元数据库。最后，基于诸元数据库利用临近点双线性插值方法计算弹道诸元并开展了仿真验证。结果表明，所提方法在保证良好精度的同时可以大大缩短弹道诸元计算时间，从而提高弹道设计的效率。     

升力式飞行器集结轨迹实时规划方法

针对升力式飞行器在强干扰环境下集结的问题，提出一种可达时间快速计算方法以及集结轨迹规划方法。在满足平衡滑翔条件的前提下，提出能量航程剖面上最快到达与最慢到达对应的两种轨迹模式，推导了两种轨迹满足的条件，基于双层牛顿割线法，构建一种可达时间实时求解方法，形成了升力式飞行器集结轨迹实时规划方法。仿真结果表明，提出的集结轨迹规划方法能够完成未知强干扰作用下升力式飞行器的实时集结轨迹规划。     

冰晶粒子运动过程中的相变特性

针对冰晶在温暖环境下运动时固液汽耦合的相变现象，基于欧拉法建立了冰晶粒子的运动相变模型和计算方法。计算了冰晶粒子在强迫对流环境下的融化相变过程，与实验结果对比验证了运动相变模型和计算方法的准确性。针对NACA0012翼型计算了冰晶绕流运动时的热力学特性，得到了冰晶粒子到达撞击表面时的融化状态与收集系数。研究了冰晶粒径大小、初始球形度、气流相对湿度和温度对运动相变的影响。结果表明：冰晶粒子运动相变模型可以有效地评估冰晶结冰风险，冰晶粒子的融化速率主要取决于粒子直径、球形度、气流温度、湿度等因素，环境温度为288 K时冰晶粒子的融化时间为27.5 s，而相同条件下环境温度为302 K时的融化时间仅有5.2 s。