模拟长杆弹侵彻混凝土靶的MCA方法

采用改进的移动元胞自动机法数值软件,对长杆弹侵彻混凝土靶进行了二维数值模拟,给出了弹靶的破坏变形过程,得到了侵彻深度与冲击速度的关系曲线。数值模拟结果与已有的实验现象吻合较好,说明该计算方法可以有效地计算和模拟高速侵彻问题。     

弹体高速侵彻钢筋混凝土的实验与数值模拟研究

为了得到钢筋混凝土目标在动能弹高速冲击作用下的破坏数据，基于大口径发射平台进行了100 mm口径卵形弹体高速侵彻钢筋混凝土靶体的实验，弹体质量为5.4 kg，靶体尺寸分为2 m × 2 m × 1.25 m 和 2 m × 2 m × 1.50 m两种，混凝土抗压强度为50 MPa，弹体侵彻速度为1 345～1 384 m/s，实验获得了弹体的侵彻深度及钢筋混凝土靶体的破坏数据。通过“钢筋混凝土全体单元分离式共节点建模方法”建立钢筋混凝土靶体模型，结合Riedel-Hiermaier-Thoma本构模型对实验工况进行计算。数值模拟给出了侵彻过程中钢筋的拉压力变化和分布规律，很好地模拟出贴近迎弹面钢筋在弹体高速冲击作用下伴随混凝土反向飞溅而产生的反向拉伸现象及靶体背面钢筋在混凝土崩落作用下发生的拉伸现象；数值模拟得到的弹体侵深数据、现象与实验结果吻合良好，实验验证了“钢筋混凝土全体单元分离式共节点建模方法”的可靠性。     

带模拟装药弹体高速冲击岩石靶时的断裂特性

针对高速侵彻过程中的弹体破碎断裂问题，本文中设计2种不同壁厚的试验弹，进行约1 000 m/s着速的高强度岩体侵彻试验，试验表明:在该高着速条件下，两种结构的试验弹体均发生完全破碎且未能有效侵入岩石靶，而岩石靶体仅在表层产生粉碎性破坏；另外，高速侵彻岩石靶的弹体头部破碎情况与侵彻金属薄靶有所区别。在试验基础上，利用Autodyn-3D建立了弹体侵彻岩石靶的物理模型，结合SPH算法与Mott失效模型对弹体破坏过程进行了数值模拟，可有效地揭示弹体破碎机理，并进一步讨论模拟装药和小范围内不同高速对弹体破坏的影响。试验结果和建立的数值模型可为研究高速侵彻中弹体结构安全提供参考。     

陶瓷靶抗射流侵彻的理论模型及数值模拟

以球型空腔膨胀理论为基础,提出了一个计算陶瓷靶板阻力的损伤模型,该模型考虑了损伤因子对陶瓷靶板弹道性能的影响.结合不可压缩流体力学理论,对射流侵彻陶瓷靶板的侵彻速度进行了理论值计算,并与未考虑损伤的侵彻速度进行了比较,该模型的计算结果更接近实验结果.建立了射流侵彻陶瓷靶板的数值计算模型,对铜射流侵彻陶瓷靶的动态破坏过程进行了研究,讨论了药型罩的锥角、壁厚对射流侵彻结果的影响,结果表明:相同锥角的药型罩,壁厚对陶瓷靶板孔径的影响较小;同壁厚的药型罩,随着锥角的增大,侵彻孔径增大.侵彻速度的数值模拟结果与理论结果进行了比较,得到了较好的一致性.     

药型罩锥角对射孔枪射流冲击过程的影响

采用ALE方法对射孔弹射流形成的过程及聚能射流对混凝土靶板的侵彻进行了数值模拟,对比了锥形药型罩的不同锥角对聚能射流形成和侵彻的影响.研究结果表明药型罩的锥角大小对聚能射流的速度和形状、射流头部和杵体的质量、侵彻的宽度和深度有着明显显著的影响.小锥角射流头部比重较小且杵体速度未达到侵彻临界值而无法起到很好的侵彻效果,大...     

靶板厚度对卵形弹丸垂直贯穿中等厚度混凝土靶的影响

为研究卵形弹丸贯穿中等厚度混凝土靶体的贯穿规律，开展直径60 mm尖卵形弹丸贯穿不同厚度混凝土靶体的侵彻实验，获得了不同撞击速度的弹丸贯穿不同厚度混凝土靶体的剩余速度规律。结合无网格SPH方法、RHT混凝土本构以及状态方程，对贯穿实验进行数值模拟，对不同工况下的弹丸过载规律以及靶体的损伤过程的分析发现:弹丸贯穿中等厚度混凝土靶体的贯穿过程分为开坑阶段、隧道稳定侵彻阶段以及靶背影响出靶阶段，在相同初始撞击速度下的靶背影响区的厚度随着靶体厚度的增加而增大。实验结果与数值模拟结果对比，表明模型能够有效模拟弹丸贯穿混凝介质问题，研究结果可为贯穿机理的研究提供参考。     

弹头部形状对侵彻影响的数值模拟研究

采用MOCL(Markoncellline)分界面跟踪算法,用二维多流体网格法的欧拉程序,对卵型头和平头动能弹侵彻混凝土平靶的过程进行了数值模拟研究,给出了侵彻过程中的物质变形情况及对应的破坏区域。基于计算结果,分析了弹头部形状对侵彻的影响。得出卵形头弹体的穿透过程是一种刺穿性模式,对靶孔周围破坏区域呈现出两头大、中间小的M型;而平头弹体则是一种挤凿性模式,对靶孔周围破坏区域较小。穿透同样的靶板,平头弹的侵彻耗能比卵形头弹侵彻耗能大,剩余速度小。     

钢纤维混凝土抗侵彻与贯穿特性的实验研究

为考察素混凝土及钢纤维混凝土（钢纤维为平直型与端钩型，体积分数0.01~0.05）抗侵彻贯穿特性，采用12.7 mm弹道炮-测速靶系统开展了初速297~848 m/s的弹道冲击实验，获得了弹丸着靶速度及对应的最大侵彻深度、弹坑直径、靶体（板）破坏形态等实验参数，并利用高速摄影系统记录了靶体（板）的动态破坏过程。实验结果的对比表明，靶板的抗侵彻贯穿性能和破坏模式与钢纤维类型及含量密切相关。当钢纤维体积分数为0.05时,端钩型钢纤维混凝土的侵彻深度与相同强度等级素混凝土相比降低约52%，且贯穿破坏后靶板碎片的数量及飞散角度大幅降低，显示了高含量异型钢纤维混凝土在抗侵彻贯穿方面的适用性。     

钢筋混凝土穿甲的数值模拟

比较了LS-DYNA中关于钢筋和混凝土相互作用的不同处理方式（即接触方法和耦合方法）对钢筋混凝土穿甲结果的影响。数值模拟S.J.Hanchak 等关于钢筋混凝土靶的穿甲实验，给出了弹靶的变形破坏过程。讨论了穿甲钢筋混凝土靶的撞击位置对终点弹道性能的影响，加深了钢筋对侵彻作用的认识。所确定的数值模拟方法可用于进一步讨论钢筋结构与尺寸、配筋率等参数对侵彻影响等问题的研究。     

间隔靶对射流侵彻影响的数值模拟和实验研究

对某聚能装药射流侵彻靶板的过程进行了数值模拟 ,得出其碰撞点附近应力分布与传统理论相符 ,侵彻深度与实验结果及工程计算结果基本相符 ;分别对该聚能装药侵彻连续靶和间隔靶的过程进行了数值模拟 ,数值模拟结果显示间隔靶的侵彻深度明显低于连续靶的侵彻深度 ,这说明侵彻开坑阶段的能耗侵深比远大于准定常阶段。为了验证间隔靶对射流侵彻的影响 ,用另一聚能装药分别对连续靶和间隔靶进行了侵彻实验 ,并排除了炸高的影响。实验结果也表明 ,间隔靶对射流侵彻的确存在着不利影响。还结合数值模拟及实验结果对传统的侵彻公式进行了修正。     

A coupled WC-TL SPH method for simulation of hydroelastic problems

A coupled weakly compressible (WC) and total Lagrangian (TL) smoothed particle hydrodynamics (SPH) method is developed for simulating hydroelastic problems. The fluid phase is simulated using WCSPH method, while the structural dynamics are solved using TLSPH method. Fluid and solid components of the method are validated separately. A sloshing water tank problem is solved to test the WCSPH method while oscillation of a thin plate and large deformation of a cantilever beam are simulated to test the TLSPH method. After validating each component, the coupled WC-TL SPH scheme is used to simulate two benchmark hydroelastic problems. The first test case shows the evolution of water column with an elastic boundary gate, and the second one investigates the breaking water column impact on elastic structures. The agreement between WC-TL SPH results and literature data shows the ability of the proposed method in simulating hydroelastic phenomena.     

SPH方法在聚能装药射流三维数值模拟中的应用

采用有限元软件LS-DYNA中的SPH算法实现聚能装药射流形成过程的三维数值模拟.将SPH方法和FEM方法与理论计算结果相比较.研究结果表明:SPH数值模拟计算过程十分稳定,避免了有限元数值模拟过程中的网格扭曲、缠绕和物质穿透等问题,而且计算精度比有限元方法更高；采用SPH方法计算的射流速度、射流长度与有限元计算结果和...     

一种考虑界面不连续的改进的有限粒子法

有限粒子法（finite particle method，FPM）作为SPH（smoothed particle hydrodynamics）方法的重要改进，有效提高了边界区域粒子的近似精度，但是当FPM处理多物理场时，在不连续界面附近的计算精度会大大降低，并且FPM必须满足的矩阵非奇异性也提高了对界面处理的要求。本文中基于DSPH（discontinuous SPH）方法，提出了一种考虑界面不连续的改进FPM—DSFPM（discontinuous special FPM）法，旨在改善FPM在界面不连续处的计算精度，从而进一步提高其计算效率和稳定性。首先，分析了DSFPM的核近似精度。其次，根据不同的工程问题，给出DSFPM处理小变形和大变形问题的算法流程。利用DSFPM、DSPH和FPM等3种方法对弹性铝块小变形碰撞冲击算例进行了模拟，通过对比分析铝块的速度和应力以及计算时间验证了DSFPM算法在非连续界面处计算精度和计算效率的优势。最后，通过结合DSFPM和DFPM（discontinuous FPM）实现了对于大变形问题的模拟。     

三维自适应FE-SPH耦合算法在多层间隔金属靶侵彻问题中的应用

鉴于有限元算法不能有效地模拟侵彻过程所产生的金属碎片, 本文中基于三维自适应FE-SPH耦合算法的基本理论, 自主开发了模拟多层间隔金属靶侵彻问题的三维FE-SPH耦合计算程序。该程序采用四面体单元对多层间隔金属靶侵彻模型进行初始离散, 计算过程中, 当四面体单元等效塑性应变超过某一设定值时, 单元自动转化为SPH粒子, 并引入有限单元-粒子接触算法和耦合算法, 实现大变形和破碎区域采用SPH方法计算, 克服有限元法单元畸变存在的问题。多层间隔靶侵彻算例分析表明, 三维FE-SPH耦合计算程序采用等效塑性应变作为转化判据计算结果较稳定, 并且能够有效地再现侵彻过程中所产生的碎片, 能够模拟侵彻碎片对后层靶的毁伤效应。     

RESEARCH DEVELOPMENTS OF SMOOTHED PARTICLE HYDRODYNAMICS METHOD AND ITS COUPLING WITH FINITE ELEMENT METHOD

拉格朗日型的有限元法和光滑粒子法在模拟材料大变形问题时各存优缺点, 而有限元与光滑粒子耦合算法实现了在小变形区域采用有限元法计算, 在局部的大变形区域采用光滑粒子法计算, 有效地综合了有限元法计算效率高和光滑粒子法能够自然地模拟材料大变形问题的特点.重点论述了有限元法、光滑粒子法以及有限元与光滑粒子耦合算法的研究现状及应用进展, 并讨论了各方法中需要进一步解决的问题.最后通过算例对3种方法的计算精度和计算效率进行了分析, 供研究人员参考.     

弹丸侵彻混凝土的SPH算法

给出了弹丸侵彻混凝土的数值计算,其中弹体作为刚体处理并划分成Lagrangian标准有限元网格,而混凝土划分成光滑粒子并经历大应变、高应变率和高压作用。为了描述混凝土的非线性变形及断裂特性,在计算中引入了Holmquist Johnson Cook本构模型及损伤模型。计算结果与实验结果的对比表明,利用光滑粒子流体动力学方法对混凝土材料进行大应变、高应变率的变形计算是有效的,并可避免网格重分或网格消蚀。此外,将光滑粒子流体动力学方法和有限元方法结合可以保持计算过程中材料界面的清晰。     

SPH-FEM接触算法在冲击动力学数值计算中的应用

为了充分发挥光滑粒子流体动力学方法(Smoothed Particle Hydrodynamics,SPH)易于处理大变形以及有限元(Finite Element Method,FEM)计算精度和效率高的优势,论文基于无网格粒子接触算法,在有限元节点处设置背景粒子,通过接触力的方式计算SPH粒子和有限单元之间的接触问题...     

固体介质中SPH方法的拉伸不稳定性问题研究进展

光滑粒子流体动力学法(smoothed particle hydrodynamics, SPH)是一种基于核估计的无网格Lagrange数值方法.它用粒子方程 离散流体动力学的连续方程, 既可以处理有限元难于处理的大变形和严重扭曲问题, 又可以 处理有限差分法不易处理的自由边界和材料界面的问题, 在固体力学中的冲击、爆炸和裂纹 模拟中具有广阔的发展前景.但是, 该算法的拉伸不稳定性(tensile instability)问题是它在固体力学领域中应用的最大障碍.对SPH稳定性分析表明, 算法不 稳定性的条件仅与应力状态和核函数的2阶导数有关.目前, 应力点法(stress points)、Lagrange核函数法、人工应力法(artificial stress)、修正光滑粒子法(corrective smoothed particle method, CSPM)和守恒光滑法(conservative smoothing)以及其他一些方法成功地改善了SPH的拉伸不稳定性, 但是每一种方法都不能彻 底解决SPH的拉伸不稳定性问题.本文介绍了SPH法的方程和Von Neumann稳定性分析的思想, 以及国内外在这几个方面的研究成果及其最新进展, 同时指出目前研究中存在的问题和研究的方向.     

流固耦合破坏分析的多分辨率PD-SPH方法

流固耦合破坏是一类涉及结构变形与破坏以及复杂自由表面现象的强非线性力学问题. 结合近场动力学(peridynamics, PD)与光滑粒子流体动力学(smoothed particle hydrodynamics, SPH)各自的优势并考虑其计算效率问题, 提出一种适用于分析流?固耦合破坏问题的多分辨率PD-SPH混合方法. 分别采用SPH和PD方法以不同的空间和时间分辨率对流体和结构进行离散与求解, 利用具有与流体粒子相同光滑长度的虚粒子处理流?固界面, 以高精度满足界面边界条件. 通过两个经典算例: 液柱静压力下弹性板的变形和溃坝流体冲击弹性闸门的变形问题, 表明提出的多分辨率PD-SPH方法兼具较高的计算精度和计算效率; 对含裂缝的Koyna重力坝水力劈裂问题进行模拟计算, 所得裂缝扩展路径与文献结果吻合, 说明该方法适用于涉及结构破坏的流固耦合问题仿真. 最后尝试采用该方法进行流体冲击作用下含裂纹混凝土板崩塌过程数值仿真, 准确描述混凝土板的断裂破坏和全过程中的流体运动. 多分辨率PD-SPH混合方法或可为流?固耦合作用下的结构损伤破坏仿真提供一种新选择.      

光滑粒子法与有限元的耦合算法及其在冲击动力学中的应用

扼要讨论了光滑粒子法的离散思想,充分利用光滑粒子法和有限元方法各自的优点,提出了一种初始时刻用有限元建模,计算过程中大变形单元自动转换为光滑粒子的耦合算法。高速碰撞的系列算例说明,耦合算法不但适宜于计算大变形冲击动力学问题,而且由于集两种方法的优点于一身,可以更高效地模拟一些高速碰撞问题,提高计算效率。