环肋圆柱壳模型在水下爆炸响应试验中的边界条件

为验证试验研究中环肋圆柱壳模型的边界条件,运用大型商业有限元程序MSC.DYTRAN分别设置3种边界条件,对试验过程进行了数值仿真研究。壳体塑性变形的仿真过程显示,迎爆面中部首先变形,随后侧爆面和背爆面相继变形,壳体最终形成数个凸凹相间的轴向凹凸带,最大塑性变形发生在迎爆面中心肋位。壳体塑性变形的仿真计算结果与试验结果的对比结果表明:固支边界条件符合试验实际情况。该结论适用于试验研究中模型边界条件的确立。     

二元翼型跨音速非定常压力计算

本文描述了用Ｅ—Ｏ型的ＡＤＩ隐式差分格式计算二维、无粘、定常与非定常低频小扰动跨音速流的方法。应用附加时间导数项的物面、尾涡面的边界条件及压力系数表达式和远场无反射边界条件,计算翼型绕其１／４弦长处轴作正弦俯仰振动时,翼面上的非定常压力分布。所得结果与文献［１］比较,具有较好的一致性。     

不同边界条件下水下爆炸气泡的数值模拟

采用有限元程序MSC.Dytran对不同边界条件下水下爆炸气泡的脉动过程进行了数值模拟。对数值模拟中的一些重要问题,如炸药药包的建模、网格密度、状态方程和静压力场的设定等问题进行了探讨。对无限水域中的水下爆炸气泡以及自由面边界、水平刚性面边界和垂直刚性面边界三种边界条件下的水下爆炸气泡进行了数值模拟,分析了水下爆炸气泡的半径、外形变化、气泡脉动压力、流场中的速度分布,结果发现边界的存在对气泡外形的影响较大,刚性边界面会使气泡脉动压力增大较多,而自由面边界对气泡脉动压力的影响较小。将数值模拟的结果与理论模型、经验公式和试验结果进行了比较分析,结果吻合较好。     

滑移边界条件的收敛性分析及应用

用速度滑移与温度跳跃边界条件代替通常假定的无滑移边界条件,可有效地提高计算流体力学模型对高空滑移流区域流动的预测精度。应用Maxwell滑移边界条件时,通过直接计算速度梯度及温度梯度而得到速度滑移和温度跳跃量的处理方法在网格较密的时候会出现迭代计算发散的问题。理论分析表明,直接计算梯度的方法使边界条件的时间推进过程等价于雅克比迭代过程,因此必须满足相应的收敛性条件。为了消除收敛性条件的限制,给出了一种在任意网格密度下均收敛的边界条件处理方法并通过数值算例验证了该方法的正确性。针对高空高超声速流动,以空天飞机为例,对比了滑移/无滑移边界条件所得结果的差异,分析了滑移效应对飞行器气动特性及热环境的影响。     

合成孔径声纳中RD和ω-k成像算法的比较

从聚焦能力和处理相关的空变特性的能力两方面 ,对合成孔径声纳中几种成像算法进行比较 .给出了RD(RangeDoppler)和ω -k算法的对水池实验数据成像的结果 .     

等离子体粒子模拟中的电磁场边界条件

针对等离子体粒子模拟中电磁场边界条件的选取问题 ,对简单吸收边界条件、Lindman边界条件、超吸收边界条件等进行了详细的讨论。结果表明 ,简单吸收边界条件计算量最小、误差最大 ,超吸收边界条件计算量最大、误差最小 ,Lindman边界条件误差较小且其反射系数对入射角的依赖性不大     

三体船旋转水动力的势流计算与分析

以三体船操纵性能预报为目的,采用势流理论三维面元法对相对较难计算的三体船旋转水动力进行了计算。根据涡元法在船体表面和尾涡面布置均匀涡格,使涡格强度满足边界条件和尾缘库塔条件,进而求解船体表面附着流场分布,并依此计算出船舶所受水动力。通过算例验证了该算法的正确性和适用性,并得出了侧体布局对三体船操纵性旋转水动力影响极大的结论。     

用面元法预报吊舱推进器水动力性能

面元法的各个环节有不同的数值处理方法,不同的方法计算结果有所不同,通过对各种方法的计算结果进行比较,得到适合螺旋桨计算的数值处理方法。用面元法对吊舱螺旋桨的水动力性能进行了计算,得到了吊舱推进器的敞水曲线和轴承力。     

面内载荷下复合材料帽型加筋板结构稳定性分析

为探究复合材料帽型加筋板结构稳定性问题,首先根据层合板理论和等效刚度法,通过分析边界条件得出复合材料帽型加筋板在面内压载作用下局部失稳的特征方程;然后,引入算例进行了数值计算,定量分析了结构参数变化对帽型加筋板临界载荷和屈曲模态的影响规律;最后,结合算例进行了相应的试验分析,试验结果和计算结果吻合较好。分析表明：帽型加筋板在面内压载作用下优先发生局部失稳,通过合理设计帽型筋条可有效避免结构整体失效。该研究结果对工程结构设计应用具有一定的理论参考价值。     

锥壳稳定方程的一般解

在本文中,用应力函数和中面法向位移表示的锥壳稳定方程被变化成为一对无奇异性的简化方程,然后得出了方程的幂级数形式的一般解。计算结果与已有的实验结果符合良好。计算分析了８种边界情况下边界条件对静水外压临界载荷的影响及解的收敛性。     

燃烧室-喷管流场数值解

本文用时间相关法计算了固体火箭发动机燃烧室─喷管亚跨声速流场数值解,控制方程用ＭａｃＣｏｒｍａｃｋ二步显格式：边界点参数用物理边界条件和参考平面上的特征方程计算。计算表明,达到收敛的数值积分步数比纯喷管的跨声速计算要多得多。虽然喷管壁上和轴线上的马赫数分布与纯喷管计算类似,但喷管中的等马赫线分布与纯喷管计算的结果［３］相差较远。     

Natural convection effects on TNT solidification inside a shaped charge mold

High Explosive Anti-Tank (HEAT) warheads and ammunitions are frequently produced by explosive casting inside an axis-symmetric mold with an inverted conical geometry in the basis. In order to prevent manufacturing defects, the solidification process must be controlled. In this study, a dimensionless solidification model has been proposed to investigate the heat transfer considering the natural convection inside the liquid explosive and the numerical simulations were performed by using COMSOL Multiphysics and Modeling Software, employing trinitrotoluene (TNT) thermophysical properties. The effect of three different boundary conditions on the top of the mold have been evaluated: convection, adiabatic and isothermal. It has been observed that solidification process was faster for convection case and slower for isothermal case, while an intermediary total solidification time value was found for adiabatic case. Moreover, liquid explosive was completely surrounded by solid explosive during the solidification process for convection case and also for adiabatic case through the end of the process. Otherwise, it was not observed for isothermal case. The natural convection effects promoted a vortex inside the liquid explosive, accelerating the heat transfer process. It has been concluded that isothermal mold top boundary condition should be preferred to prevent manufacturing defects, avoiding high thermal stress.     

预处理方法求解低马赫数流场的边界条件

通过求解预处理后的二维可压缩Euler方程 ,数值试验了常见的三种进口边界条件的提法及其对稳定性和收敛性的影响。数值试验结果表明 ,预处理后 ,在低马赫数时基于黎曼不变量的边界条件为反射边界条件 ,而另两种边界条件为无反射边界条件     

钝体高超声速气动加热与结构热传递耦合的数值计算

气动加热与结构热传递耦合问题在航天和工程应用领域非常重要。分别采用松耦合与紧耦合方法,数值模拟了高超声速二维圆管绕流的流场与结构传热耦合的非定常过程。在紧耦合方法中,流场部分采用基于Navier Stokes方程的有限体积法,将AUSM+格式与时间方向的显式多步Runge Kutta法结合;结构传热部分采用基于二维热传导方程的Galerkin有限元法。流场与结构区通过交界面的热流和温度边界条件实现耦合。计算结果分别与实验、文献做了对比,结构内部温度变化关系以及壁面的热流分布均较好地吻合。两种耦合方式的计算结果对比表明,对于流场特征时间远小于结构传热特征时间的问题,松耦合方法计算效率高,精度与紧耦合方法接近。     

潜艇螺旋桨区域腐蚀电场模型特征及奇异峰成因分析

针对边界元法推算潜艇腐蚀电场分布时边界条件影响电场信号特征,导致螺旋桨区域产生奇异峰现象的问题,以实测船用921合金钢和镍青铜螺旋桨材料极化曲线作为建模边界条件,重点对比分析了恒电位和非线性边界条件下的电场特征,并通过建立阻抗谱参数下的船壳-螺旋桨电化学阻抗等效电路,分析了潜艇腐蚀电场螺旋桨区域产生奇异峰的原因。仿真结果表明:潜艇电场特征分布及螺旋桨区域奇异峰现象与船体材料电化学极化状态有关,合理设定非线性极化边界参数可达到削弱奇异峰现象、平滑腐蚀电场模型的效果。     

具有时滞和变消耗率的比率型Chemostat模型稳定性分析

研究了一类具有时滞和变消耗率的比率型Chemostat模型,讨论了边界平衡点和正平衡点的局部稳定性以及Hopf分支的存在性;利用比较原理和构造Lyapunov泛函,给出了边界平衡点和正平衡点的全局稳定性的充分性条件;并通过数值模拟验证了理论结果。