球锥体高超声速绕流的俯仰阻尼导数的数值计算

采用Navier Stokes方程描述球锥体外形非定常振动流场 ,在Etkin理论下给出俯仰阻尼导数的计算公式。对定常流场的计算采用ADI形式的NND格式 ,对非定常流场的计算采用四步Runge Kutta方法 ,并引入变系数残值光顺技术加速收敛。将俯仰阻尼导数计算结果与实验及工程计算进行了比较 ,并数值研究了其随振动中心的变化规律。     

履带车辆静动液辅助制动系统仿真与实验研究

针对履带车辆静动液辅助制动系统液力辅助制动和液力液压联合辅助制动工况,建立对应工况下的数学模型进行仿真研究,并根据相似原理进行静动液辅助制动系统实验．仿真和实验结果表明,所建立的数学模型能够准确反映系统的实际制动性能,且该辅助制动系统能够为车辆提供可靠的制动．     

不同攻角下AUV回收过程中的水动力规律

通过对AUV(Autonomous Underwater Vehicle)回收过程中运动规律的研究和对纳维—斯托克斯方程(Navier-Stokes equations)及其求解方法的分析,基于CFD(Computational Fluid Dynamics)理论建立了AUV回收运动的三维计算模型,并使用CFD软件完成其物理模型的建立、网格的划分及求解设置。在此基础上对AUV回收过程中的阻力、升力干扰进行分析计算,得到不同攻角下AUV所受到的水动力干扰与AUV和回收平台之间距离的相互联系。     

掠飞末敏弹耦合运动动导数计算

以新型掠飞末敏弹为研究对象,提出了基于欧拉转动定理和滑移网格技术的复杂角运动模拟方法,利用著名的罗德里格斯转换矩阵插值求得弹箭在每个时间步的角速度修正值,并指定给球形滑移网格区。通过对非定常气动参数进行求解辨识,分析了不同马赫数下掠飞末敏弹滚转运动对其俯仰组合动导数和升力系数动导数的影响规律。结果表明:所提角运动模拟方法可有效消除姿态角计算的累积误差,实现对弹箭任意给定角运动的准确模拟;弹箭滚转运动对俯仰组合导数和升力系数动导数的辨识结果均存在显著影响,在进行弹箭动导数计算和稳定性分析时需充分考虑俯仰耦合效应的影响。     

液压水锤效应引起液体喷溅特性及其影响因素试验研究

为研究液压水锤效应引起的液体喷溅特性及其影响因素,进行了高速破片撞击充液容器的试验,测试了液体内的空腔振荡特性、压力分布特性和容器外的液体喷溅特性。试验结果表明:液压水锤效应引起的液体喷溅有两个不同的阶段。第一个阶段发生在空腔生长到最大体积后的400~700 μs内,喷溅液体的头部呈箭状;在第二个阶段出现多次形状相似的脉动喷溅,且单次脉动均发生在每次空腔溃灭之后,喷溅液体具有伞状头部与线状尾部。液体内的压力和侵彻孔的形态共同影响液体的喷溅速度,相对喷溅速度与破片的撞击速度成反比。喷溅液体在运动过程中存在速度波动,且随着喷溅次数的增加,液体喷溅速度逐渐下降,速度拐点出现的位置逐渐接近侵彻孔。     

改进的接触算法及其在光滑粒子流体动力学中的应用

采用以黎曼解描述粒子间相互作用的接触算法,模拟一维情况下包含间断的流场。在弱间断的条件下,采用弱波近似的黎曼解来描述粒子间的相互作用;在强间断的条件下,则引入非迭代黎曼解来描述粒子间的相互作用。并引入Taylor展开思想,提高了接触算法在自由边界的计算精度。该算法无需引入人为粘性,程序结构简洁。数值计算表明,改进后的接触算法能很好地描述强间断,并有效改善了原始接触算法在自由面的计算缺陷。     

旋转弹体入水弹道的计算方法

借助已有的实验资料与经验 ,给出了弹体入水空泡的数学模型 ;分析了旋转弹体与空泡壁相互作用引起的弹体上的流体动力特性 ,提出了一种预测旋转弹体入水弹道的计算方法     

What you should know about simulation and derivatives

Derivatives (or gradients) are important for both sensitivity analysis and optimization, and in simulation models, these can often be estimated efficiently using various methods other than brute‐force finite differences. This article briefly summarizes the main approaches and discusses areas in which the approaches can most fruitfully be applied: queueing, inventory, and finance. In finance, the focus is on derivatives of another sort. © 2008 Wiley Periodicals, Inc. Naval Research Logistics, 2008     

用动网格法计算理想平板的颤振导数

提出了一种数值模拟振动的理想平板绕流场 ,并提取其气动导数的方法。时间相关的不可压N S方程采用Projetion 2格式解耦 ,关于中间速度的动量方程的时间和空间离散采用二阶半隐格式 ,压力Possion方程迭代用多层网格法加速收敛。分别计算平板作竖向强迫振动或扭转强迫振动的气动力 ,用动网格法考虑平板和气流的耦合作用。由计算得到的气动力用最小二乘法确定 8个气动导数。计算结果和理想平板的Theodorsen理论值有较好的一致性     

机体/推进一体化飞行器动导数频域计算方法

运用Etkin非定常气动力模型,将谐波平衡法应用于复杂外形动导数辨识。计算结果表明:谐波平衡法和时域方法的动导数辨识结果一致,与实验值吻合,验证了程序的正确性和谐波平衡法的可靠性。采用谐波平衡法对类X-51高超声速机体/推进一体化飞行器WR-A进行动导数辨识,并与时域方法的计算结果比较,考察了该方法在复杂外形下数值模拟非定常流动和动导数辨识的能力。实验结果表明:对于WR-A这种复杂外形飞行器,谐波平衡法也能在保证精度的同时具有较高的计算效率。