跨声速母弹空腔流的数值模拟

将ＳＬＩＰ格式和ＬＵ－ＳＳＯＲ隐式迭代方法结合，数值模拟了尖头母弹在开舱前后的跨声速定常流场。采用修正的Ｂ－Ｌ湍流模型，用轴对称Ｎ－Ｓ方程求解零攻角时母弹在Ｍａｃｈ数０．９４～０．９８ 的临界气动特性和空腔流特性。开舱前母弹计算结果与实验以及Ｓａｈｕ 的结果能够很好吻合， 开舱后空腔流的计算结果也给出了合理的理论分析。     

DARPA2潜艇模型定常流动 黏性流场和水动力计算

为了进一步验证SA湍流模式模拟分离流动问题的能力,采用SA(C_υ=30)湍流模型,数值求解定常不可压缩的RANS方程,计算DARPA2潜艇模型定常流动黏性流场和水动力;速度-压力耦合采用SIMPLE方法处理,对流项采用二阶迎风格式,隐式迭代求解方程组。得到的数值计算结果与实验值符合较好。潜艇不同部分对水动力的贡献及受攻角的影响是不同的,升力系数随攻角的增加,艉附体的作用明显增强,而指挥台的作用不是很明显;指挥台和艉附体分别存在时,水动力系数随攻角变化的趋势与全附体时是一致的,曲线不存在奇异性变化。     

导弹大攻角高机动飞行特性分析与仿真

以新一代无翼式气动布局导弹为背景,在大攻角气动特性理论基础上,分析了在导弹主动段、攻角在180°范围内变化的导弹飞行特性,得出了90°内大攻角最大转弯能力和180°攻角附近最大降速效果的结论。进行了俯仰平面S型轨迹机动飞行仿真,仿真结果说明了大攻角飞行的高机动特性。     

发动机内流场湍流数值模拟

以可压缩的全Ｎ－Ｓ方程为控制方程,用显式有限差分格式对燃烧室－－喷管轴对称湍流内流场进行了数值模拟。计算中采用ＭａｃＣｏｒｍａｃｋ两步显格式和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,得到了管状装药发动机内流场的数值解,取得了令人满意的结果。此方法对发动机的设计有实际的工程参考价值。     

测试投稿，大家不要管哈~

采用非定常双时间步长方法、刚性旋转动网格技术,结合空间隐式格式数值分析了高超声速钝锥体定常与非定常流场特征,并运用最小二乘法辨识出钝锥体的俯仰阻尼导数。通过与相应的实验数据、工程估算值进行比较,分析并验证了壁面网格雷诺数、钝头比、攻角与质心位置对钝锥体俯仰阻尼导数的影响。同时采用五种差分格式,其中包含五阶精度WENO与高精度WNND格式,对钝锥体俯仰阻尼导数进行了计算,结果表明不同无粘差分格式对俯仰静导数影响较小,但不同格式之间俯仰动导数结果存在差异,分析认为差分格式之间不同的粘性分辨率应是差异的主要来源。     

高超声速钝锥体俯仰阻尼导数影响因素分析

采用非定常双时间步长方法、刚性旋转动网格技术,结合空间隐式格式数值分析了高超声速钝锥体定常与非定常流场特征,并运用最小二乘法辨识出钝锥体的俯仰阻尼导数。通过与相应的实验数据、工程估算值进行比较,分析并验证了壁面网格雷诺数、钝头比、攻角与质心位置对钝锥体俯仰阻尼导数的影响。同时采用五种差分格式,其中包含五阶精度WENO与高精度WNND格式,对钝锥体俯仰阻尼导数进行了计算,结果表明不同无粘差分格式对俯仰静导数影响较小,但不同格式之间俯仰动导数结果存在差异,分析认为差分格式之间不同的粘性分辨率应是差异的主要来源。     

舱门开启状态下UUV小攻角水动力特性数值模拟

为了掌握舱门开启状态下UUV小攻角水动力特性,首先基于风洞模型实验验证了所用数值计算方法的准确性与有效性,进而采用该方法对舱门开启状态UUV阻力、升力及力矩特性进行了数值模拟研究,分析了舱门开启后载荷舱及舱门对UUV表面压力分布的影响.计算结果表明,载荷舱前后端面形成的较大压力差,导致舱门开启状态UUV阻力显著提升,且阻力关于零攻角表现为不对称性,主要原因在于负攻角条件下载荷舱粘压阻力随攻角绝对值的增大而增大;由于载荷舱的流体阻滞作用,零攻角工况下舱门下表面平均压力高于上表面,UUV零攻角升力为正,零升力攻角约为-3.,且零攻角力矩为正力矩,零力矩攻角约为-1°.     

钝倒锥体动导数数值工程模拟

分别采用求解非定常薄层近似的N－S方程数值模拟有粘扰动流场及非定常内伏牛顿流理论计算了钝倒锥外形的高超声速俯仰阻尼导数Cma＋Cmq。在非定常计算中通过引入“亚迭代”过程提高非定常流场的计算精度。两种方法计算结果一致，数值计算可应用于3＜M＜20，0°＜α＜20°时的一般流态俯仰阻尼导数计算。     

再入飞行器高超声速三维粘性流场的数值模拟

用数值方法研究大钝头倒锥体再入飞行器高超车速完全气体有攻角绕流的流场特性。控制方程为完全气体三维Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,数值方法采用激波捕捉的ＮＮＤ隐式算法,计算区域为含底部的全流场。文中给出了流场特性和气功力系数,并与实验结果进行了对比。     

飞船再入舱三维化学非平衡流数值模拟

对飞船再入舱高超声速化学非平衡三维流场进行数值模拟。控制方程为含化学反应源项的全Ｎａｖｉｅｒ－ Ｓｔｏｋｅｓ 方程组, 化学模型为高温空气的七组元模型, 其组元成分为Ｎ２ 、Ｏ２、ＮＯ、Ｎ、Ｏ、ＮＯ＋ 、ｅ－ 。差分格式采用张涵信院士提出的ＮＮＤ激波捕捉格式并用全耦合方法处理流动方程及化学反应方程, 对化学反应源项采用全隐式处理, 运用时间预处理技术以加快收敛。同时, 也对完全气体模型进行计算, 以分析真实气体效应。与实验数据对比, 说明本文结果是可信的。     

三维高超音速干扰流场数值模拟与实验验证

采用Osher Chakraverthy的TVD格式、Baldwin Lomax湍流模型和LU SSOR隐式方法求解了完全NS方程 ,数值模拟了三维高超音速绕流与横向喷流干扰流场。并在高超音速炮风洞中开展了喷流实验研究 ,对该喷流流场的数值计算结果进行实验验证     

Investigation of normal,lateral, and oblique impact of microscale projectiles into unidirectional glass/epoxy composites

Spacesuits and spacecraft must endure high velocity impacts from micrometeoroids. This work considers the impact of 100 μm diameter projectiles into composite targets at velocities from 0.5 km/s to 2 km/s. This work begins by presenting an energy-based theoretical model relating depth of penetration (DoP) and impact force to impact velocity, characteristic time, and threshold velocity and force. Next, this work compares numerical simulations of normal impact on composites to the theoretical model. Numerical simulations are conducted with LS-DYNA and the well-known composite model, MAT-162. The numerical models consider unidirectional S2-glass fiber reinforced SC-15 epoxy composite laminates. The numerical model shows good correlation with the theoretical model. The numerical model also investigates lateral impact, parallel to the fiber direction, and oblique impact at angles from 30° to 82.5°. This work decomposes oblique impact into normal and lateral components, and compares them with normal and lateral impact results. The results show good correlation of the normal component of oblique results with the theoretical model. This numerical and theoretical study focuses on DoP, velocity, and penetration resistance force as functions of time. The theoretical model and numerical simulations are used to determine new DoP parameters: characteristic time of depth of penetration and threshold impact velocity. These models are a first step in developing the capability to predict DoP for oblique, microscale, high-speed impact on composite materials.     

逃逸飞行器喷流干扰流场的分区数值模拟

基于拼接网格系统 ,应用三维区域边界格式 ;通过区域边界格式对区域边界上的通量进行守恒性处理 ,发展了一个超音速复杂流场分区算法 ;对不同攻角下的逃逸飞行器喷流干扰流场进行了分区数值模拟 ,取得了较好的结果     

The non-isothermal gravimetric method for study the thermal decomposition kinetic of HNBB and HNS explosives

The hexanitrostilben (HNS) is a thermally stable explosive that can be prepared from hexanitrobibenzyl (HNBB). Therefore, the investigation of thermal stability of HNBB can be important in the yield of preparation of HNS. The decomposition kinetic of HNBB and HNS are studied by non-isothermal gravimetric method. The TG/DTG curves in non-isothermal method are obtained in range of 25°C–400 °C at heating rates of 3 °C/min, 5 °C/min, 8 °C/min, 10 °C/min and 12 °C/min. The data of weight-temperature are used for calculation of activation energy (E_a) of thermal decomposition reactions by methods of Ozawa, Kissinger, Ozawa-Flynn-Wall (OFW) and Kissinger-Akahira-Sunose (KAS) as model-free methods and Strink's equation as model-fitting method. The compensation effect is used for prediction of mechanism and determination of pre-exponential factor (lnA) of the decomposition reaction. A reduction 60 kj/mol for the average of activation energy of thermal decomposition reaction of HNBB is obtained versus HNS. This result shows the lower thermal stability of HNBB in comparison to HNS. The Avrami equation (A_3/2) with function f(α) = 3/2(1-α)[-ln(1-α)]^1/3 indicates the predicted mechanism for thermal decomposition reaction both explosives.     

高超声速飞行器低可探测性滑翔弹道优化方法

为进一步提高高超声速飞行器的突防性能,提出高超声速飞行器低可探测性滑翔弹道优化方法。考虑飞行器180°×360°方向的雷达散射截面,针对原数据尖峰多、收敛难的难题,运用高斯滤波法对其进行预处理,既不改变原数据趋势又加以平滑,提高优化问题收敛性能。为使计算所用雷达散射截面数据具备较强的保真性,采用三次样条插值方法调用离散数据计算实时雷达散射截面。完成了高超声速飞行器低可探测性滑翔弹道优化问题的建模,以探测概率为目标函数,运用hp自适应Radau伪谱法优化求解,采用逐步计算策略进一步提高优化效率和收敛性能。与传统最短飞行时间弹道对比表明,该方法有效降低了飞行器被雷达发现的概率。     

圆柱形装药驱动轴向预制破片飞散特性

为了增强柱形战斗部轴向威力,在无壳柱形战斗部底面布置单层离散的预制破片。开展圆柱形TNT装药驱动轴向预制破片飞散试验,获得预制破片的最大初速、飞散角等特征参数;运用LS-DYNA软件对装药驱动预制破片过程进行数值模拟,阐述预制破片群飞散过程;对装药驱动整体平板理论计算公式进行改进,获得预制破片的最大初速。结果表明:破片初速理论计算结果、数值计算结果和试验结果吻合良好;随着与装药底部中心距离的加大,破片初速、径向飞散角分别近似呈“抛物线”减小、增大;试验实测、理论计算得到的破片最大初速值超过2500 m/s,试验实测的径向飞散角最大约为22°,而周向飞散角则普遍较小,均值在5°以内。     

基于梯度直方图的交叉点检测方法

提出基于梯度直方图的交叉点检测方法:统计像素点邻域梯度信息生成特征描述子,通过判定准则筛选候选交叉点,基于欧氏距离对候选交叉点进行聚类后,利用灰度加权法定位图像交叉点坐标。仿真实验验证了算法的正确性与可行性,结果表明:该方法适用于交叉角为65°～90°范围内的交叉点检测,定位精度优于0.6像素,检测率大于90%;该方法可以很好地抵抗图像旋转,定位精度和检测率基本不受图像旋转影响。该方法具有良好的抗噪性能,可满足一般透视变换下的图像交叉点检测需求,适用于结构光三维重建中的系统标定和形貌测量。     

Dynamic modeling and parameter identification of a gun saddle ring

In this study, a theoretical nonlinear dynamic model was established for a saddle ring based on a dynamic force analysis of the launching process and the structure according to contact-impact theory. The ADAMS software was used to build a parameterized dynamic model of the saddle ring. A parameter identification method for the ring was proposed based on the particle swarm optimization algorithm. A loading test was designed and performed several times at different elevation angles. The response histories of the saddle ring with different loads were then obtained. The parameters of the saddle ring dynamic model were identified from statistics generated at a 50° elevation angle to verify the feasibility and accuracy of the proposed method. The actual loading history of the ring at a 70° elevation angle was taken as the model input. The response histories of the ring under these working conditions were obtained through a simulation. The simulation results agreed with the actual response. Thus, the effectiveness and applicability of the proposed dynamic model were verified, and it provides an effective method for modeling saddle rings.     

前缘和转轴影响翼摇滚特性的数值模拟

采用耦合求解N-S方程和Euler刚体动力学方程的数值模拟方法研究80°后掠三角翼摇滚特性,讨论机翼前缘构型和滚转轴位置对三角翼摇滚特性的影响。N-S方程的离散采用Roe格式和含双时间步的LU-SGS方法,刚体动力学方程的离散采用二阶精度的差分方法,通过交替求解流动和运动控制方程组的耦合策略,模拟80°后掠三角翼自激摇滚的非定常过程。针对转轴安装于上表面的细长三角翼,研究前缘上削尖、下削尖和双面削尖三种构型的翼摇滚特性,分析前缘构型对摇滚振幅及分岔攻角的影响。考虑到上削尖和下削尖构型在翻转后相当于前缘构型互换但转轴位置不同,进一步考察了滚转轴位置对三角翼摇滚特性的影响。结果表明:在一定的攻角下,三种不同前缘构型的三角翼均能形成大幅自维持的摇滚现象,其中上削尖前缘的三角翼振幅最大,双面削尖前缘次之,下削尖前缘三角翼振幅最小;攻角增大到30°以后,双面削尖和上削尖前缘的三角翼会发生翻转现象,翻转180°以后维持等幅振荡。     

Effect of bending temperatures on the microstructure and springback of a TRIP steel sheet

Transformation-induced plasticity (TRIP) steel possesses high strength and formability, enabling the use of a thinner gauge material and allowing for the fabrication of complex shapes. In this research, we measured the effect of bending temperatures on the microstructure and air-bending springback angle of TRIP steel at temperatures from 25 to 600 °C. Real-time in situ X-ray diffraction and scanning electron microscopy were used for pre- and postbending analysis. As the prebending temperature increased from 25 °C to 600 °C, the retained austenite (RA) volume fraction decreased, and the RA transformed to bainite at temperatures above 400 °C. The springback angle was positively correlated with the prebending RA volume fraction, with the smallest springback angle achieved at 400 °C. Additionally, the springback angle was positively correlated with the bending angle, because the RA transformation ratio contributed to increased strain hardening. Further microstructure analysis revealed that the RA became elongated in the tension direction as the bending temperatures increased.