基于sph方法的聚能射流侵彻混凝土靶板数值模拟-强洪夫.pdf

《基于sph方法的聚能射流侵彻混凝土靶板数值模拟-强洪夫.pdf》由会员分享，可在线阅读，更多相关《基于sph方法的聚能射流侵彻混凝土靶板数值模拟-强洪夫.pdf（9页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第36卷第4期2016年7月爆炸与冲击EXPLOSl0N AND SHOCK WAVESV0136，No4Jul。2016基于SPH方法的聚能射流侵彻混凝土靶板数值模拟。强洪夫，范树佳，陈福振，刘 虎(火箭军工程大学动力工程系，陕西西安710025)摘要：在完全变光滑长度SPH(smoothed particle hydrodynamics)方法的基础上，利用FOtt等提出的修正SPH方法处理在求解多介质大密度问题时的数值不稳定性问题，运用HolmquistJohnson-Cook本构模型处理混凝土在冲击载荷下的变形和损伤问题，对聚能装药射流侵彻混凝土靶板的过程进行了数值模拟，同时利用LS-D

2、YNA非线性有限元程序进行对比，分析了2种方法得到的混凝土von Mises应力变化、射流头部特定节点处的速度变化及裂纹演变，验证了SPH方法的准确性。分析了另外2种不同尺寸的靶板在射流侵彻作用下的破坏形式，结果符合射流侵彻物理规律，表明该方法适合模拟聚爆炸与冲击等大变形破坏等问题。关键词：爆炸力学；侵彻；光滑粒子动力学；混凝土；破坏中图分类号：0389 国标学科代码：1303530 文献标志码：A二级串联随进战斗部n1作为攻击跑道目标的有效手段，在现在战争中发挥着重要作用。其第1级战斗部即采用聚能装药爆破后形成高速射流，对混凝土目标侵彻形成孔道，以便于第2级战斗部跟进，继续侵彻一定深度后爆破

3、，从而对目标造成不可恢复的破坏。聚能射流侵彻混凝土的主要特点是：炸药能量高，聚能效应明显，形成高速射流，混凝土在高速冲击作用下由于其拉压不等的特性其破坏的模式复杂等。因此，对于该过程的数值模拟对于串联随进战斗部的研究与设计具有重要意义。传统用于计算该过程的方法主要是基于现有商业软件中的成熟算法，如多物质Euler法、Lagrange法以及ALE(arbitrary LagrangianEulerian)法等，这些方法大都采用有限差分法(finite differencemethod，FDM)、有限元法(finite element method，FEM)或有限体积法(finite volume

4、 method，FVM)等基于网格的数值方法进行离散求解，计算中不可避免地出现网格的扭曲和缠绕(针对Lagrange网格离散)或界面追踪复杂、精度较低(针对Euler网格离散)的问题。光滑粒子流体动力学(smoothed particlehydrodynamics，SPH)方法作为一种无网格粒子方法，在对射流进行弹塑性流体动力学计算，对混凝土进行大应变、高应变率的变形计算时，可避免网格重分及算法耦合，因此非常适合此类问题的求解。最早采用SPH方法模拟爆炸的可追溯到JWSwegle等口，而LDLibersky等3q1最先将SPH方法运用于高速冲击领域，随后，MBLiu等61采用SPH方法模拟了聚

5、能装药的爆轰过程，Qiang Hongfu等71运用FOtt等83提出的修正SPH方法对聚能装药射流过程进行模拟，分析了不同的起爆方式对射流的影响。本文中，拟在强洪夫等910提出的完全变光滑长度SPH方法的基础上，结合修正SPH方法8，处理爆炸和冲击过程中密度和光滑长度变化剧烈的问题以及多介质界面问题。采用SPH方法计算混凝土损伤大多引入JohnsonCook11|、SCG1幻等本构模型，与混凝土的特性不符，未能捕捉混凝土在高速冲击下出现的裂纹扩展等细节，且采用SPH方法计算包括炸药的爆轰、药罩的挤压、射流的形成及发展、射流高速侵彻混凝土、混凝土损伤破碎等过程在内的复杂全过程问题未见报道。因此

6、，本文中将HolmquistJohnsonCook本构模型1 3引入SPH方法，处理混凝土在高速冲击作用下的变形损伤问*收稿151期：20141203；修回日期：20150324基金项目：国家自然科学基金项目(51276192)；国家重点基础研究发展计划(973计划)基金项目(61338)；火箭军工程大学创新型基金项目(EPXY0806)第一作者：强洪夫(1963一)，男，博士，教授，博士生导师；通信作者：范树佳，fan_shu_jia163com。万方数据第4期 强洪夫，等：基于SPH方法的聚能射流侵彻混凝土靶板数值模拟 517题，对聚能炸药爆炸挤压形成射流，进而侵彻混凝土靶板的整个过程进行

7、数值模拟，并对侵彻过程中射流头部特定点处的速度变化进行分析，采用传统网格方法进行对比验证。同时为更充分地描述混凝土靶板在高速射流侵彻作用下的破坏损伤效应，对另外2种不同厚度的靶板进行设计，并开展数值实验，探讨SPH方法用于模拟诸如聚能射流侵彻混凝土等涉及爆炸与冲击大变形多介质问题的适用性。1 SPH基本方程在SPH方法中，通过对问题域的粒子离散来求解系统的状态，系统中物理量，(r)及其导数V厂(r)的SPH形式为：一re,f(r，)w。 (1)卜)一E mf(r，)vW。 (2)一 。 ，)V 。 (2)卜】式中：Dt、D、r分别为粒子的质量、密度和位置矢量；W。，一w(r。一r，h)为核函数

8、，h为光滑长度，通常选用三次样条核函数1“。为更好地模拟爆炸与冲击、大变形大扭曲等密度和光滑长度变化剧烈的问题，本文中采用完全变光滑长度SPH方法，同时为了很好地解决药型罩和爆轰气体间密度梯度较大带来的间断面不稳定的问题，采用修正SPH方法101对聚能装药射流过程进行模拟，结合后的方程组如下：ddlo。imi妻i=1V。w。+丢(警+鲁)筹面dvi一一善N m，(磐地)V。w。仁(1+ihi军ow，o，一 (3)警一帆)警警邓式中：W。一W(x。一x，)为核函数，它的选取直接影响计算的误差和稳定性，通常选用三次样条核函数；V。W。为核函数对粒子i坐标的空间导数，v一v。一q；仃为总应力张量；i

9、为人工黏度；h为插值核宽度的一种度量，称为光滑长度，表示w不显著为零的取值范围，控制着SPH粒子的影响域，通常设定l X。一X，lkh，点一2时，W一0，这里h为空间和时间的函数，在连续方程中将其对时间求导，即：面dhi一一上鲁监 (4)q-ddt一一 ，d _0。 一式中：d为空间维数。对于光滑长度变化率dh，dt与密度变化率dp，dt之间相互耦合而带来的难以显式求解的问题，本文中采用迭代法93求解密度方程和光滑长度。在计算过程中，为使粒子分布更有序，消除由于分布不均匀引来的粒子非物理聚集的问题，采用JJMonaghan143提出的XSPH对速度场进行光滑：华一，一y当，。w。 (5)T=，

10、一夕，。，w lJ Lo，Q _lDv式中：(0e1)为常数，通过施加邻近粒子的影响使自身的运动速度与邻近粒子的平均速度相近，一般e取03。2炸药及药型罩材料模型本文模拟中，炸药选用TNT，其爆轰速度由实验153测得为6 930 ms。对于爆轰气体，使用标准万方数据爆 炸 与 冲 击 第36卷JWL状态方程叫得到气体压力：舭，u)一A(卜卷)exp(-跏)+B(卜惫)exp(-R2v)+半 (6)式中：po=1 630 kgm3为TNT炸药的初始密度，v=p。p，P为爆轰气体的密度，e为爆轰气体的比内能；A、B、R。、R：和硼为由实验结果拟合得到的系数17，且A一3712 GPa，B一3231

11、 GPa，R1415，R2095，叫一03。对于药型罩材料，选用铜，状态方程选用MieGriineisen状态方程16：P(p，e)一(1一专r7)pH(10)+昂P 7一点一1 (7)声Hfn。r+bo叩2+f。叩3 】7o (8)【no 7 叩：0a。一伽C：， b。一a。E1+2(S。一1)， f。一a。1+2(S。一1)+3(S。一1)2式中：po一8 530 kgm3为铜的初始密度，r一199，C。一3 940 ms，S。一i489。3 混凝土本构模型及状态方程HJC本构模型1 3分为强度模型、累计损伤模型和状态方程3部分，模型曲线分别如图13所示。31屈服强度模型归一化的屈服强度函

12、数为：d一EA(1一D)+Bp。N(1+Clnl。) (9)式中：盯=af：，仃为真实等效强度；P。=pf：为归一化压力，P为真实的压力，f为准静态单轴抗压强度；D为损伤变量(0D10)；+一；。为等效应变率，；为应变率，；。一10 S-1为参考应变率；A为给定静水压力下完全断裂强度与无损伤强度的比值，B为压力强化因数，N为压力硬化因数，C为应变率敏感因数。屈服强度模型曲线如图1所示，其中S为混凝土所能够达到的归一化最大强度，T+一Tf：为归一化最大抗拉静水压力，丁为材料所能承受的最大抗拉静水压力。32累计损伤模型损伤值为由等效塑性应变和塑性体积应变引起的损伤累积起来的形式：Dy A6，p+A

13、，zp (10)一e；十卢；式中：。和卢。分别为等效塑性应变和塑性体积应变在一个计算周期内的增量；：+卢：为恒定压力下损伤塑性应变，其表达式为：+肛：一D1(P+丁+)Dz (11)矿一D=O(Undamaged)-弋p图1混凝土屈服强度模型曲线Fig1 Curves for yield strength model of concrete群州?f，圣!j唑，7式p+图2混凝土累计损伤模型曲线Fig2 Curve for cumulative damage model of concrete式中：D，和Dz为混凝土的损伤常数。累计损伤模型曲线如图2所示。33 状态方程混凝土靶板在冲击载荷下的压

14、缩变形分为弹性阶段、过渡阶段和压实阶段3个阶段，状态方程曲线万方数据第4期 强洪夫，等：基于SPH方法的聚能射流侵彻混凝土靶板数值模拟 519如图3所示，图中的P-。t和P。n分别为下式中的P和P。，210ck、1crush和川。k分别为下式中的pI、2c和1pIo(1)弹性阶段(O。)在这一阶段，混凝土发生可以恢复的弹性变形，加载与卸载的方程相同，为：P=K。口 (12)式中：弹性体积模量K。一P。2。，P。和2。分别为单轴强度抗压实验中得到的压碎压力和压碎体积应变，肚=pp。一1为单元体积应变，lD和JD。分别为单元的密度和初始密度。(2)过渡阶段(1。卢岸。1)在这一阶段中混凝土中的气体

15、逐渐被挤压出去，混凝土产生破碎性断裂，加载时：死lD)图3混凝土状态方程曲线Fig3 Curves for equation of state of concreteP=P。+K。(2一。) (13)式中：K，一(PlP。)(2。-一P。)为塑性体积模量，pl为压实压力；户。lPgpo一1为压实体积应变，对应于P处的体积应变，pg为颗粒密度。该区的卸载是通过相邻区域间插值的一条路径进行的，卸载的状态方程为：P一(1一F)K。+FK。 (14)式中：插值因子F一(2。一2。)(刖一2。)，2。为卸载前达到的最大体积应变。(3)压实阶段(卢岸，1)该阶段混凝土完全被压碎，加载的状态方程为：PK12

16、+K2口2+K313 (15)式中：2一(口一21)(1+121)为修正的体积应变；K-、Kz和K。为混凝土材料常数。卸载的状态方程为：PKl2 (16)混凝土所能承受的最大拉伸强度为T(1一D)，状态方程曲线如图3所示。4聚能射流侵彻混凝土靶板的数值模拟计算中，采用的XSPH参数一03；核函数选用三次样条核函数1 4】，混凝土靶板材料参数1胡为：剪切模量G，148 GPa；初始密度IDo，2 450 kgm3；A，790 MPa；B，1 600；N，061；o，10；C，0007；，48 MPa；S，70；Dl，004；D2，10；混凝土破碎的最小塑性应变：。，001；户。，160 MPa；

17、1。，0001；户I，800 MPa；21，01；T，25 MPa；Kl，85 GPa；K2，一171 GPa；K3，208 GPa；。数值模拟中采用的聚能装药模型及混凝土模型如图4(a)所示，药柱宽度为4000 mm，装药头长度为2075 mm，药孔长度为4200 mm，张角为60。，药型罩厚度为300 mm，炸高为4000 mm，混凝土靶板尺寸为100 mm30 mm。炸药采用点起爆的方式，起爆点为(0，007)，具体粒子配置如图4(b)所示，粒子间距均为05 mm，其中TNT炸药粒子数为26 766，药型罩粒子数为4 480，混凝土粒子数为12 000。光滑长度取15倍粒子间距，时间积分

18、采用蛙跳格式，时间步长为01”s。本文中利用LSDYNA非线性有限元程序进行对比验证。计算模型由炸药、药型罩、空气和混凝土靶板4部分组成。其中，对炸药、药型罩和空气采用Euler网格建模，对混凝土靶板采用Lagrange网格建模，在混凝土靶板与空气和药型罩材料间采用耦合算法。聚能装药是线性的，考虑到模型的对称性，为了减少计算时间，将模型简化为二维平面对称问题并建立其12模型。在对称面上施加对称边界约束，在混凝土靶板边界施加固定约束，把靶板侧面设定为无反射边界。对炸药、药型罩、空气和混凝土靶板均选用与SPH数值模拟中参数相同的本构模型和状态方程。万方数据520 爆 炸 与 冲 击 第36卷图5给

19、出了采用SPH方法和LSDYNA程序计算过程中炸药爆轰及药罩被挤压、射流形成过程的对比。图6给出了采用SPH方法和LS-DYNA程序计算过程中射流侵彻混凝土形成开孔、混凝土损伤及混凝土完全破坏的过程的对比。可以看出，炸药爆轰波逐渐将药型罩向中轴线方向挤压形成射流和杵体，随后混凝土在射流的侵彻下发生压缩和剪切变形，产生剥落，靶板弹坑的上部呈漏斗状。随着射流的不断侵彻，靶板背面首先由于反射波到达壁面后的反射拉伸作用，导致背面出现拉伸损伤，部分脱落，开孔周围则逐渐出现不同程度的层裂，并且裂纹逐渐向外扩展，最终形成一个漏斗状损伤区域，中间部分出现粗细均匀的孔道，裂纹形态与图7中的实验结果基本相符。通过

20、对比验证，发现利用2种数值模拟方法得到的结果均符合射流侵彻物理规律，但LS-DYNA程序通过定义失效删除单元产生的裂纹非常规则，数量少，基本没有在裂纹上产生小裂纹的状况，且不能模拟侵彻后碎片的物理特性，裂纹模拟效果受到网格划分大小的影响，SPH方法则恰恰克服了以上困难。由图6可知，在侵彻过程中，前2段射流发挥了主要作用。因此，分别选取SPH方法及LSDYNA程序计算结果中前2段射流的头部顶点粒子A、B进行分析。图8为2种方法的计算结果中A点和B点粒子速度随时间的变化曲线，由此可知，药型罩在002 ms时受到炸药爆轰波的挤压开始逐渐加速，在约010 ms时，在SPH方法的计算结果中A点粒子速度达

21、到最大值3 2479 ms，在LS-DYNA程序的计算结果中A点粒子速度达到最大值3 4574 ms。射流头部在014 ms时开始与靶板接触，速度大幅降低，造成混凝土的损伤，随着侵彻的深入，速度逐渐减小为零。2种方法的不同之处在于，SPH方法的计算结果中第1段射流在020 ms时速度基本降为零的同时第2段射流开始侵彻，而在LS-DYNA程序的计算结果中，第1段射流在014 ms时与混凝土靶板接触，速度逐步降低，在017 ms左右时，第2段追上第1段射流，2段射流合为一段继续侵彻，直到速度基本降至零。l】黟藿姜亘呈iConcrete图4算例模型结构Fig4 Construction of the

22、 modelISPH图5聚能射流形成过程对比Fig5 Comparison of the formation of jet flows万方数据第4期 强洪夫，等：基于SPH方法的聚能射流侵彻混凝土靶板数值模拟V。n Mises stress脚a【二二二卫20 49 78 110 140 160 190 220 250 280(a)t=O20ms 一。 (b)t=O25 ms(c)t-O30 msSPHLSDYNA SPH图6聚能射流侵彻过程中不同时刻，混凝土靶板中Von Mises应力的分布Fig6 Von Mises stress distribution in concrete targe

23、t penetrated by shaped charge jet at different times图7聚能射流侵彻的混凝土靶Fig7 Concrete target penetrated by shaped charge jetO-05一l0，b一15豆吾一2o一2S一3O-35-o1 O O1 O2 O3tms图8射流头部特定点处速度一时间曲线Fig8 Velocity-time curves at special points of jet head相对于LSDYNA中通过删除网格单元形成裂纹的处理方法，SPH不删除粒子单元而是通过显示混凝土靶板的损伤值D来表述材料失效形成的裂纹，损

24、伤后的SPH粒子按无黏性的流体粒子进行计算，可以完整地再现裂纹的损伤演化过程，损伤状态及裂纹扩展分布情况如图9所示，计算得到的靶板最终损伤状态如图10所示。由图10可知：采用LS-DYNA得到的计算结果中弹坑直径与采用SPH方法得到的计算结果中弹坑直径差别较大，其主要原因在于LS-DYNA中理想空气与混凝土靶板的重叠区域小；采用SPH方法得到的计算结果中靶板左右壁面出现裂纹，而采用LS-DYNA得到的计算结果中没有出现，其主要原因在于SPH中靶板的左右边界是自由边界，冲击波到达左右壁面时由于反射波的作用开始出现拉伸损伤，而LS-DYNA靶板定义的是无反射边界，反射波不能到达左右壁面。万方数据5

25、22 爆 炸 与 冲 击 第36卷100(a)t=O20 ms-50 0 50 100xlmmHl()一50(b)t-O30 ms-50 0 50 100图9在聚能射流的侵彻下不同时刻混凝土靶板的损伤情况Fig9 Damage in concrete target penetrated by shaped charge jet at different timesLSDYNA SPH图10不同方法模拟得到的裂纹长度及分布情况Fig10 Length and distribution of crack simulated by different methods为了更充分地描述混凝土靶板在高速射

26、流侵彻作用下的破坏损伤效应，采用上述SPH方法对另2种尺寸的靶板进行了数值实验，粒子数分别为4 000和40 000，计算结果见图11，图12为选取的射流头部特定点进行追踪得到的速度一时间曲线。从图1112可以看出，当靶板尺寸较小时，在较短的时间内只需射流头部一段即可将靶板完全击穿，射流孔处损伤均匀；靶板较厚时，需要3段或更多段射流才可将靶板击穿，而且在射流孔处不会形成大面积的损伤，而是形成不同程度的损伤裂纹延伸至靶板内。80604020至 o一20一406()-80fal 100mm10mill at 022ms一604()-20 0 20 40 60xmmfbl 200mm50mm at0

27、401TlglI)0 -50 0 50 100xnun图11另外2种不同尺寸的混凝土靶板在聚能射流侵彻作用下的损伤Fig1 1 Damage in other two concrete targets with different sizes penetrated by shaped charge jetsDI磊两弱钻驺筋I三踟o加踟一一一一lD踮巧的L导站：Gb覆万方数据第4期 强洪夫，等：基于SPH方法的聚能射流侵彻混凝土靶板数值模拟 523图12另外2种不同尺寸的靶板聚能射流侵彻作用下特定节点的速度一时间曲线Fig1 2 Velocitytime curves at special po

28、ints of other two concrete targets with different sizespenetrated by shaped charge jets5 结 论运用无网格SPH方法计算得到了聚能装药射流形成、发展及高速冲击混凝土整个过程，并与传统网格进行了对比验证，得知：(1)采用SPH方法能有效解决在高速冲击、大变形问题的计算中不可避免地出现Euler网格界面追踪复杂或Lagrange网格扭曲和缠绕的问题。基于完全变光滑长度SPH方法与修正SPH方法相结合的方法，有效解决了爆炸与冲击中光滑长度变化剧烈的问题以及多介质界面上大密度间断带来的计算不稳定问题。(2)运用HJ

29、C本构模型和状态方程描述混凝土在高速冲击下的变形损伤特性，不仅得到了破孔的形状，而且对由于壁面反射波造成的靶板背面拉伸损伤及破孑L周围裂纹的扩展等变形损伤细节描述精确，与实际物理现象符合。采用2种方法得到的计算结果基本一致，验证了SPH方法的准确性。同时，采用SPH方法计算得到的裂纹形态与实验结果吻合更好，更合理。(3)随着靶板尺寸的增大，射流孔径逐渐减小，裂纹数目逐渐增加，裂纹的细度及延伸的长度逐渐增加，破坏的能量逐渐向周围扩散。参考文献：Eli段建，杨黔龙，周刚，等串联随进战斗部侵彻混凝土靶实验研究J爆炸与冲击，2007，27(4)：364369Duan Jian，Yang Qianlon

30、g，Zhou Gang，et a1Experimental studies of a tandem follow-through warhead penetrating concrete targetJExplosion and Shock Waves，2007，27(4)：364369E2 Swegle J W，Attaway S WOn the feasibility of using smoothed particle hydrodynamics for underwater explosioncalculationsJComputational Mechanics，1995，17(3)

31、：15卜1683Libersky L D，Petscheck A GSmoothed particle hydrodynamics with strength of materialsCTrease H，FrittsJ，Crowley WProceedings of the Next Free Lagrange ConferenceNY：SpringerVerlag，1991，395：2482574Libersky L D，Petscheck A G，Carney T C，et a1High strain Lagrangian hydrodynamics：A threedimensionalS

32、PH code for dynamic material responseJJournal of Computational Physics，1993，109(1)：67755Liu M B，Liu G R，Zong Z，et a1Computer simulation of high explosive explosion using smoothed particle hydrodynamics methodologyJComputers&Fluids，2003，32(3)：3053226Liu M B，Liu G R，Lam K Y，et a1Meshfree particle simu

33、lation of the detonation process for high explosives inshaped charge unlined cavity configurationsJShock Waves，2003，12(6)：509520万方数据524 爆 炸 与 冲 击 第36卷7 Qiang Hongfu，Wang Kunpeng，Gao Weiran Numerical simulation of shaped charge j et using multiphase SPHmethodJTransactions of Tianjin University，2008，1

34、4(1)：495499r8ott F，Schnetter EA modified SPH approach for fluids with large density differencesJArxiv Physics Eprints，2003：31129强洪夫，高巍然完全变光滑长度SPH法及其实现J计算物理，2008，25(5)：569575Qiang Hongfu，Gao WeiranSPH method with fully variable smoothing lengths and implementationJChineseJournal of Computational Phys

35、ics，2008，25(5)：56957510强洪夫，王坤鹏，高巍然基于完全变光滑长度SPH方法的HE爆轰过程的数值试验J含能材料，2009，17(1)：2731Qiang Hongfu，Wang Kunpeng，Gao Weiran Numerical study of high explosive detonation process using SPHmethod with fully variable smoothing lengthsJChinese Journal of Energetic Materials，2009，17(1)：2731r11Johnson G R，Cook W

36、 HA constitutive model and data for metals subjected to large strains，high strain ratesand high temperaturesCProceedings of the Seventh International Symposium on BallisticsHagueNetherlands，1983：571-57412 Steinberg D J，Cochran S G，Guinan M wA constitutive model for metals applicable at high strain r

37、ateJJournal of Applied Physics。1980，51(3)：1498150413HoImquist T J，Johnson G R，Cook W HA computational constitutive model for concrete subjected to largestrains，high strain rates，and high pressurescProceedings of the 14th International Symposium on BallisticsQuebec，Canada，1993：59160014Monaghan J JSmo

38、othed particle hydrodynamics-JReports on Progress in Physics，2005，68(8)：17031759153 Liu G R，Liu M B光滑粒子流体动力学：一种无网格粒子法M韩旭，杨刚，强洪夫，译长沙：湖南大学出版社，2005：19519716Livermore Software Technology CorporationLS-DYNA keyword users manualMLivermore：LivermoreSoftware Technology Corporation，2012：17451 7Liu M B，Liu G

39、R，Zong Z，et a1Computer simulation of high explosive explosion using smoothed particle hydrodynamics methodologyJ-Computers and Fluids，2003，32(3)：305322Numerical simulation on penetration of concrete target by shaped charge j etwith SPH methodQiang Hongfu，Fan Shujia，Chen Fuzhen，Liu Hu(Power Engineeri

40、ng Department，Xian HiTech Institute，Xian 710025，Shaanxi，China)Abstract：On the basis of the smoothed particle hydrodynamics(SPH)method with fully variablesmoothing lengths，the modified OttSchnetter SPH method was used to cope with the computationalinstability in the large density gradient problemsThe

41、 HolmquistJohnsonCook constitutive modelwas used to cope with the deformation and damage of concrete under impactBy comparing with thecorresponding ones simulated by the nonlinear finite element program LS-DYN A，the changing courses of the von-Mises stress and cracks in the concrete targets，and the

42、velocities at the some specialpoints of the jet head were analyzed，which proved the feasibility and accuracy of the SPH methodOther two concrete plates with different sizes penetrated by shaped charge jets were also simulatedThe results are in good agreement with the physical principleSo this method can be used to deal withthe multimaterial and large deformation problems such as detonation and impactKey words：mechanics of explosion；penetration；smoothed particle hydrodynamics；concrete；damage(责任编辑张凌云)万方数据