牵顶伞在降落伞拉直过程中的作用

采用有限元方法,建立降落伞回收系统的多质点动力学模型,研究了牵顶伞在拉直过程中对主伞性能的影响.模型中将柔性大变形的伞绳、伞表、连接带离散成若干绳段,各绳段假设为质量集中在中心的质点,各质点间以阻尼弹簧连接,其运动由对应绳段的重力、气动力和张力确定.利用该模型对某批次的空投试验进行了仿真,给出了张力、速度、位移等参数的瞬态变化.经过分析表明,牵顶伞能够改善主伞张力分布,抑制伞衣项部横向摆动.     

基于流体速度势的非线性晃动问题的ALE有限元分析

基于流体速度势，本导出了非线性晃动问题的ALE有限元格式。通过引入ALE运动学描述，使网格的运动独立于流体的运动，通过采用与贮箱固连的非惯性参考系建立控制方程，使网格只须相对贮箱更新，因此运动边界和自由液面的跟踪算法不仅非常简单，而且不会引起网格畸变．本分别采用有限元法和有限差分法进行空间和时间离散，中的算例分别模拟了贮箱受横向激励和俯仰激励时流体的非线性晃动问题，所得结果证实了本方法的有效性与可靠性。中的方法还可进一步推广用于流体—结构耦合问题．     

火星探测器减速着陆过程中若干问题的研究

围绕火星探测器减速着陆过程中涉及的开伞控制技术、降落伞拉直过程的绳帆现象、降落伞充气性能以及物伞系统动力学特性等关键技术问题进行了研究。通过对探测器进入轨道的分析,提出了一种以动压为控制目标的自适应开伞控制方法,并对其原理和实现方法进行了研究。采用降落伞拉直过程多刚体、多质点、多自由度动力学模型,对绳帆现象进行了数值仿真分析,发现开伞攻角是绳帆现象的主要成因,并提出了设计控制措施。采用降落伞轴向-径向动量守恒充气模型研究了大气密度对降落伞充气性能的影响,得到了大气密度对降落伞充气时间、充气距离、充气过程中伞衣面积变化及开伞载荷的影响规律或特点。利用物伞系统动力学模型,研究了物伞系统减速下降过程的动力学特性。本文研究内容对于中国火星探测器减速着陆系统的工程设计、试验及性能评估等均具有较好的参考价值。     

二维网格局部动态细化的关键技术及其在切削仿真中的应用

为了解决仿真过程中网格畸变导致的计算终止以及计算精度与效率低等问题,本文利用Python对ABAQUS有限元仿真软件进行了二次开发,提出一种二维局部网格动态细化算法。该算法通过研究网格细化准则和细化方法以及物理场的传递过程,完成了整个算法程序的开发,并应用该算法实现了在二维切削仿真中局部网格动态细化。与全局加密网格的模型相比,在保证仿真精度误差为5%以内的情况下,采用局部网格动态细化时计算效率提高了210%。最后,通过试验验证了采用该算法所建立的二维切削仿真模型的准确性,仿真结果与试验结果基本吻合。     

高速碰撞中Lagrange有限元方法及其应用

根据连续介质力学守恒方程和有限元离散方法，针对三维高速碰撞现象的特性，应用瞬时最小闰能原理，完整地推导了高速碰撞有限元基本列式和离散方程，并给出了高速碰撞问题的中心差分求解算法，讨论了非线性Ｌａｇｒａｎｇｅ动力有限元程序数值计算中应重视的问题，本文给出的典型算例说明本方法对高速碰撞问题分析的有效性。     

板壳结构非线性动态随机响应的模态降阶分析

航空航天飞行器中的板壳结构在强噪声、高温和机械力载荷共同作用的环境中常会产生大振幅非线性随机振动响应,对非线性动态随机响应问题的模态降阶分析方法进行了综述。先介绍了模态降阶分析的基本原理,然后总结模态降阶技术的实现方法。分析每种方法的优缺点和适用范围,提出了目前研究存在的一些问题及今后的研究方向。     

Ahmed 模型的流固耦合数值计算方法探索与实验验证

流固耦合效应会影响汽车行驶安全性,对该效应的数值模拟方法的可信度需要验证。以 Ahmed 模型作为数值仿真和实验验证的对象,数值计算是流体计算与结构计算两部分耦合,配合风洞实验的位移、压力测量技术与流场显示技术(烟流法、丝线法)与之比较,数值模拟与实验所得位移结果之间的差异始终在20%以内。从流场结构角度分析得到流固耦合效应影响气动力的机理主要在于流场中涡的拓扑结构变化。为将来数值研究时变侧风下车辆的流固耦合问题提供实验及理论依据。     

Identification of Magnetic Bearing Stiffness and Damping Based on Hybrid Genetic Algorithm

Identifying the stiffness and damping of active magnetic bearings(AMBs)is necessary since those parameters can affect the stability and performance of the high-speed rotor AMBs system.A new identification method is proposed to identify the stiffness and damping coefficients of a rotor AMB system.This method combines the global optimization capability of the genetic algorithm(GA)and the local search ability of Nelder-Mead simplex method.The supporting parameters are obtained using the hybrid GA based on the experimental unbalance response calculated through the transfer matrix method.To verify the identified results,the experimental stiffness and damping coefficients are employed to simulate the unbalance responses for the rotor AMBs system using the finite element method.The close agreement between the simulation and experimental data indicates that the proposed identified algorithm can effectively identify the AMBs supporting parameters.     

CALCULATION OF THREE DIMENSIONAL BOUNDARY LAYER EQUATIONS USING DIFFERENTIAL QUADRATURE METHOD

ＣＡＬＣＵＬＡＴＩＯＮＯＦＴＨＲＥＥＤＩＭＥＮＳＩＯＮＡＬＢＯＵＮＤＡＲＹＬＡＹＥＲＥＱＵＡＴＩＯＮＳＵＳＩＮＧＤＩＦＦＥＲＥＮＴＩＡＬＱＵＡＤＲＡＴＵＲＥ　ＭＥＴＨＯＤＴａｎｇＤｅｎｇｂｉｎ；ＣｈｅｎｇＧｕｏｗｅｉ（ＤｅｐａｒｔｍｅｎｔｏｆＡｅｒｏ...     

密实织物伞形状和应力计算研究

建立了充满稳定状态下密实织物伞形状和应力计算的理论模型。将充满的伞看成柔软的壳体，用力学原理建立伞绳（伞衣上径向加强带）和伞衣幅中线的力平衡方程。伞绳和伞衣处理成非线性弹性构件。考虑伞衣中线上的子午向应力，得到了一个双轴应力模型。结合伞充满状态下几何大变形的几何关系式，得到一组含６个微分方程的非线性微分方程组。从伞顶孔开始，利用伞顶处的边界条件，离散积分至伞底，校对伞底处的另一边界条件，如不满足再从伞顶开始进行迭代，直至满足伞底边处的边界条件。由此编成进行非线性迭代数值求解程序ＣＳＬＡＰ。选用一具典型的平面圆形伞进行了验证，结果很好     

降落伞初始充气阶段数值模拟

根据降落伞的结构和其在充气过程中的受力特性,以某平面圆形伞为原型,建立了伞衣初始充气过程中的计算流体力学与结构动力学的耦合模型.首次考虑了充气过程中折叠伞表的张开问题,建立了更接近降落伞物理模型的初始充气阶段伞表质点结构和受力方程.对流场的变化采用了准定常假设,利用simple算法数值模拟求解RNG κ-ε湍流模型下的雷诺平均N-S方程以获得每一状态伞表张开部分与折叠部分交界处质点的压差系数.把数值计算结果和试验结果及经验值比较,得到如下结论:(1)初始充气阶段伞衣外形变化为:整个阶段,伞衣展开部分外形基本保持较光滑的直筒形状,而非喇叭形.与试验结果相比,计算结果较真实地反映了初始充气阶段伞衣外形的变化情况.(2)当无因次充气时间为0.27左右时.初始充气阶段结束,伞衣投影面积随充气时间呈线性变化,计算值与实验值接近.     

一种新型锥形降落伞的设计

介绍了两种锥形降落伞的设计。一种伞型使用了零透气量纺织材料，另一种伞型使用了微透气量纺织材料。两种伞型均使阻力系数值达到１．２，从而显著减少了伞衣面积，降低了伞衣重量。     

降落伞开伞过程的试验研究

探讨了降落伞充气过程中开伞动载的产生原理及其影响因素,提出了在风洞中模拟真实开伞动载的可行性。在此基础上,对小型平面圆形伞的开伞过程进行了试验,研究了开伞过程中伞衣形状和伞衣所受载荷之间的关系。试验表明:开伞载荷受降落伞阻力和附加质量变化率产生的力影响最大;在开伞过程中伞衣会发生“呼吸”现象,伞衣所受最大载荷出现在伞衣第一次完全充满时。同时,变速风洞的伞衣载荷研究表明,变速风洞可较好地模拟实际开伞情况。     

二维复值Ginzburg-Landau方程的一个高阶紧致ADI差分格式

对二维复值金兹堡朗道(Ginzburg-Landau,GL)方程提出一个基于时间分裂的高阶紧致交替方向隐式有限差分格式。本文通过时间分裂法将GL方程分裂成一个非线性子问题及两个线性子问题,对非线性子问题以及其中一个线性子问题均通过精确积分进行计算,并对另一线性子问题构造紧致交替方向隐式差分格式进行数值计算。实际计算中,在每一时间步,利用追赶法求解一族常系数三对角线性代数方程组,从而使得算法既具有较高精度又拥有较快的计算速度。数值实验表明该算法在时间和空间方向分别具有二阶和四阶精度,并模拟了方程的一些动力学行为。     

基于ALE／SUPG同步交替法求解流固耦合问题

该文基于预测一多步校正方法提出了流固耦合ALE／SUPG有限元同步交替求解法。流体动量方程采用SUPG有限元进行空间离散以消除对流项引起的数值振荡;时间域上的积分采用预测一多步校正（predictor—multicorrector）方法,对流体和结构同时进行预测;在多步校正的过程中对流体固体交替求解,从而在时间上达到同步推进。数值算例与Hughes,Chen等人的结果进行了比较,表明了该文方法的正确性。     

有限元结构分析并行计算的若干研究进展

有限元结构分析并行处理在大规模科学与工程计算中占有重要地位，本文简要介绍作者在这一领域内的研究工作及其成果。涉及单元刚度矩阵的并行计算、有限元方程组的并行直接解法、有限元方程组的并行迭代解法、结构动力分析并行直接积分法、广义特征值问题的并行算法以及EBE（Element-By-Element）技术在有限元结构分析并行处理中的应用等研究领域。特别提出了“伪单元向量”与“单元分组技术”等概念与新方法；首先考虑了矩阵向量积的多自由度问题；系统、深入而又全面地研究了EBE策略在结构分析并行处理中的应用，填补了国内空白并突破了国际现有成果。     

非线性双曲型守恒律的一类应用通量限制器的二步二阶格式

本文提出了计算非线性双曲型守恒律弱解的一类新的二步二阶精度有限差分格式,这类高精度非线性差分格式的构造应用了三点单调差分格式的数值通量和通量限制器。文章证明,在拟CFL条件限制下,格式为总变差衰减差分格式(TVD格式),因此不发生常用的诸如二步二阶Lax-Wendroff格式、Maccomark格式和其他古典二阶格式计算不连续解时在其周围出现的剧烈振荡现象。 为了减少Roe推广技术引起的计算复杂性,文中我们基于Davis,S.F.的思想,使用了Rusanov数值通量,把标量差分格式推广到方程组,构造了求解双曲型守恒律方程组的新的预测-修正格式,较之其他一般的预测-修正格式而言,仅增加了很少的计算复杂性,初步的数值试验表明,本文提出的格式无振荡和具有极好的激波分辨率。     

OPTIMAL DESIGN AND AERODYNAMIC CALCULATION OF WING CONFIGURATION OF CIVIL AIRCRAFT

ＯＰＴＩＭＡＬＤＥＳＩＧＮＡＮＤＡＥＲＯＤＹＮＡＭＩＣＣＡＬＣＵＬＡＴＩＯＮＯＦＷＩＮＧＣＯＮＦＩＧＵＲＡＴＩＯＮＯＦＣＩＶＩＬＡＩＲＣＲＡＦＴ￥ＷａｎｇＬｉａｎｇｙｉ（ＤｅｐａｒｔｍｅｎｔｏｆＡｅｒｄｄｙｎａｍｉｃｓ，ＮＵＡＡ２９ＹｕｄａｏＳｔｒｅ...     

基于CFD/CSD/CAA耦合的声气动弹性仿真技术在壁板颤振中的应用

提出了一种基于计算空气动力学/计算结构动力学/计算气动声学（CFD/CSD/CAA）耦合的声气动弹性时域仿真方法。空气动力学和气动声学分别采用大涡模拟和Lighthill声比拟理论进行计算。采用D’Alembert原理和有限元方法建立了二维简支壁板在超声速流作用下的非线性动力学方程,并在时域内进行求解。其中几何非线性采用Von Karman大变形理论描述。以二维弹性壁板作为研究对象,研究了超音速流场中噪声、气动力、和壁板之间的相互作用所产生的非线性颤振现象。     

Ejection Separation Characteristic Analysis of Parachute Container Cover from Return Capsule for Lunar Exploration

The parachute container cover ejection separation is the first and foremost motion for the return capsule recovery system, which is related to the success of a recovery system. Adopting the computational fluid dynamics （CFD） simulation and flight dynamics coupling method, the parachute container cover ejection separation is simu lated. The rationality of the ejection separation speed and dynamic characteristics of the separation process is ana- lyzed. Meanwhile, the influences of angle of attack, Mach number and ejection separation speed on the parachute container cover ejection are also investigated. Results show that the ejection separation speed design is reasonable. It has a certain design margin for parachute container cover to escape from the wake region, and to pull out the drag parachute completely. The results may provide a theoretical basis for recovery system engineering design of the lunar exploration project.