一种改进的结构拓扑优化水平集方法

针对有限的结构设计区域,界面力的大小和作用的位置以及部分给定位移的结构边界始终保持不变的要求,为了解决基于水平集方法的结构边界演化停滞等问题,导出了基于水平集演化的结构拓扑优化所需的法向速度场.研究了由该法向速度场导出的序列解的收敛特性,并建立了一套新的结构柔顺度为目标函数,体积为约束条件的水平集演化算法,给出的算例验证了该方法的正确性和有效性.     

基于结构拓扑随机变异的水平集优化方法

基于结构拓扑导数的水平集优化方法,利用结构拓扑导数信息和通过不断减少结构体积的方式来确定需开孔的结构位置,需选用最大设计区域结构作为初始优化结构.该方法不适合求解结构体积等式约束的优化问题.为了解决上述问题和克服水平集方法不能在结构内产生新孔洞的困难,提出了一种基于结构拓扑随机变异的水平集优化方法.引进了以小概率随机方式进行结构拓扑变异的结构优化方案.设计了变异算子,讨论了提出方法的收敛性情况.最后,结合考虑结构最大设计区域限制的结构拓扑优化的水平集方法,建立了一套新的涉及结构柔顺度作为目标函数,体积作约束条件的水平集演化算法.给出的算例验证了该方法的正确性和有效性.     

N-S方程在非结构网格下的求解

在Ｒｏｅ的矢通量差分分裂的基础上,吸收了ＮＮＤ格式的优点,提出了一种非结构网格下求解Ｅｕｌｅｒ方程和Ｎ－Ｓ方程的高分辨率高精度迎风格式．这种格式具有捕捉强激波和滑移线的良好性能．在时间方向上采用了显式和隐式两种解法．文中还给出了自适应技术．最后,成功地完成了ＧＡＭＭ超音速前台阶绕流、二维平板无粘激波反射、三维Ｈｏｂｓｏｎ叶栅流动、ＶＫＩ叶栅流动、Ｃ３Ｘ叶栅流动的数值模拟,得到了满意的结果     

壳装炸药殉爆实验和数值模拟

进行了壳装固黑铝炸药殉爆实验,通过观察残留炸药、壳体和见证板变形,判断被发炸药的爆炸情况,得到了炸药临界殉爆距离。建立了壳装炸药殉爆实验计算模型,采用非线性有限元计算方法,对壳装固黑铝炸药殉爆实验进行了数值模拟。计算中采用预设壳体单元破片方法描述主发炸药壳体破片的形成和破片对被发炸药的撞击起爆,炸药临界殉爆距离的计算结果与实验结果基本一致。主要是主发炸药中部的壳体破片撞击到被发炸药,被发炸药起爆位置也在装药中部。炸药壳体厚度主要影响破片速度和质量、被发炸药的防护性能,从而影响炸药临界殉爆距离。 更多还原     

PBX-01炸药水中爆轰产物状态方程研究

提出通过水中实验确定炸药的水中爆轰产物JWL状态方程参数的方法;选择PBX-01高能炸药进行水中实验,利用ANSYS/LS-DYNA程序建立炸药的水中实验模型,将实验结果与数值计算结果进行对比,确定PBX-01炸药水中爆轰产物的JWL状态方程参数。研究结果显示,圆筒实验确定的JWL参数在反映炸药水中爆轰产物的膨胀状态时有所不足,水中实验确定的JWL状态方程参数能够更准确地描述PBX-01炸药水中爆轰产物的膨胀过程,因此对水中爆炸的研究需要通过水中爆炸实验建一套状态方程参数。

     

对碰波作用下金属圆管材料参数对其运动特征影响的数值模拟

基于两点起爆实验,采用动力学有限元程序,在不考虑圆管破裂的前提下,详细分析了同一炸药爆轰下金属圆管各材料力学参数对其运动的影响规律。计算结果表明:弹性模量越大,泊松比越大,金属圆管的膨胀速度越大,但它们对圆管膨胀的位移和速度影响极小,在工程中可以忽略;屈服应力对金属圆管的运动有一定影响,屈服应力越大,对碰部位的鼓包范围越宽,鼓包峰值越小,但圆管各处的膨胀量均落后于屈服应力较小的情况;密度对圆管运动有较大影响,密度越大,受鼓包影响范围越宽,鼓包峰高越小。     

JB-9014炸药超压爆轰产物的状态方程

根据P.K.Tang等提出的对炸药爆轰产物超压状态方程建模时只对JWL状态方程CJ等熵线中高压指数项做修正的研究思路,首先给定超压状态下内能等熵线的修正项,再根据热力学定律对内能等熵线求微分而得到沿压力等熵线的修正项。对JB-9014炸药超压爆轰冲击Hugoniot实验数据和声速实验数据同时进行拟合,得到了3个JWL状态方程在超压爆轰状态下的修正项,并进行了分析与比较。得到超压修正项的方法简单,3组超压修正项与P.K.Tang的修正项一样,都能很好地拟合超压Hugoniot数据。在实验数据范围外,对超压状态下的声速-压力实验数据的拟合精度有所差别。     

不同升温速率下复合药柱烤燃实验与数值模拟研究

为了研究不同升温速率下复合炸药JO-9159/JB-9014烤燃实验的热反应规律,建立了复合炸药的烤燃模型,利用有限元程序LS-DYNA3D对不同结构的复合药柱在烤燃过程中的热响应情况进行了数值模拟,并利用实验进行了验证,结果显示模拟结果可信。利用已建立的模型对5 K/h、3 K/min和10 K/min等3种不同升温速率下复合药柱烤燃过程进行了数值模拟,结果表明：升温速率和装药结构的不同对复合药柱的点火时间和位置有较大影响,随着升温速率的增大,点火时间变短,点火位置由药柱的中心处逐渐移至药柱的两端边缘,升温速率较小时,复合药柱的热安定性取决于内部高能炸药的特性,升温速率较大时,复合药柱的热安定性与单一钝感药柱性能近似。因此,只有在较大的升温速率下,钝感炸药内部嵌入高能炸药才能既提高整体药柱的威力,又保证其具有较好的热安定性。     

强爆轰驱动飞片的数值模拟研究

在一维流体动力学程序中,使用了JWLT状态方程,对强爆轰驱动飞片的实验模型进行了数值模拟研究,得到了炸药爆轰产物的压力和密度分布,验证了强爆轰的存在,同时得到了二级飞片自由面的速度曲线,与实验测量结果符合较好,为解释实验结果和设计新的强爆轰实验提供了计算依据,同时验证了程序的正确性。     

PBX-9501强爆轰产物状态方程的数值模拟

利用P.K.Tang的JWL强爆轰产物状态方程改造了DEFEL二维弹塑性流有限元源程序,并用改造后的DEFEL程序模拟了一系列飞片撞击PBX-9501炸药的一维强爆轰过程,给出了炸药驱动二级飞片(铝)的速度。计算得到的产物状态参数与实验测量结果符合较好,验证了增加修正项的新状态方程可较好地描述高能炸药的超压爆轰的高压流体动力学过程。     

以DSD理论和LS方法为基础的程序燃烧法

给出了二维弯曲爆轰波后产物流场计算方法。爆轰波阵面传播规律满足Detonation Shock Dynamics (DSD)理论并用level set (LS)方法计算,波阵面传播规律与波后流场的耦合通过程序燃烧法实现,反应进程变量可作为LS函数的函数给出。爆轰波从刚性细管道向粗管道传播产生绕射的二维计算结果表明,化学反应速率不影响波后流场分布,只影响反应区结构。此方法可用于钝感炸药的驱动计算问题。     

高浓度氩气稀释气体爆轰波临界管径和临界间距关系

建立圆管及环形管道系统研究临近极限下爆轰波在管道内传播失效机理。选用C₂H₂+2.5O₂+70%Ar气体,采用光纤探针测量爆轰波在管道内传播速度,用烟迹法记录管道内爆轰波胞格结构。结果表明：初始压力远大于爆轰极限压力时,爆轰波在管道内以稳定速度传播;随着初始压力的减小,爆轰波速度逐渐降低;当初始压力一定时,爆轰波速度随着管道尺寸的减小而逐渐减小;当初始压力达到临界压力时,爆轰波在进入到管道内后其速度会逐渐衰减直至爆轰波完全失效。对于不同几何尺寸的圆管与环管,通过引入无量纲参数d/λ及w /λ(d为圆管管径,w为环管间距,λ为爆轰胞格尺寸)得出,爆轰波在管道内传播的临界圆管直径为环形间距的2倍,与理论模型结果相吻合,验证了稳态气体基于爆轰波波面曲率的失效机理。