摘要
随着航天科技的迅速发展以及航天发射任务的急剧增加,大量的空间碎片散落在近地轨道上,构成了复杂的空间碎片环境,对航天器的安全构成了极大的威胁。如何确保航天器的安全,成为当前航天器设计领域中一个十分重要的问题。由于超高速碰撞问题的复杂性,碰撞碎片云的特性还不是很清楚,防护结构的设计还有待进一步研究。
光滑粒子流体动力学(简称SPH)方法作为最早提出的无网格方法,现在已逐渐被应用到许多流体力学相关领域。它具有逻辑简单,可以避免有限元方法中的网格缠绕和扭曲问题等优点,因而特别适合模拟超高速碰撞问题。
本文系统介绍了SPH方法,并对该方法数值计算相关问题作了较深入的探讨。在此基础上,采用SPH方法对空间碎片超高速碰撞碎片云特性以及Whipple防护结构设计问题作了数值研究。
本文的主要工作及其结论如下:
1、利用Fortran语言自行编写了基于直角坐标系的二维和三维SPH计算程序,并成功引入Johnson-Cook拉伸累积损伤本构模型、Johnson-Cook粘塑性本构模型以及SCG本构模型,由此得到的数值结果与实验较为吻合,说明自编程序能较好的模拟各类碰撞动力学问题,特别是超高速碰撞问题。
2、从计算精度、拉伸不稳定性、虚粒子方法、初始光滑长度和非均匀粒子间距界面等五个方面讨论了SPH方法数值计算相关问题。对三种拉伸不稳定性处理方法作了比较,结果表明,与守恒光滑法相比较,Monaghan型拉伸不稳定处理方法以及基于Lagrangian核函数的稳定化粒子方法能更好地消除拉伸不稳定性。虚粒子方法是处理固壁边界的一种有效方法。然而,如果在对称结构的对称面处(为减少计算量,对称面通常被当作固壁边界处理)采用虚粒子方法,则在对称面附近会存在一个粒子数明显偏少的低密度区域,这显然是不合理的。研究表明,采用XSPH速度纠错公式在一定程度上可以改善计算结果。分析了初始光滑长度对数值模拟结果的影响。结果表明,初始光滑长度越大,弹丸材料表现出的刚性越小,h0的合理取值范围应为d0 3、从碎片云形状、质量、速度和动量等方面研究了超高速正碰撞碎片云特性。结果表明,球形弹丸的碎片云前端呈半球形分布,柱状弹丸的碎片云前端有一锥状突出,当长径比足够小时圆盘形弹丸形成的碎片云内部会产生柱状结构。提出了识别碎片的自然边界判定方法,结合近邻粒子搜索技术,得到了碎片云的质量分布情况。给出了累积碎片百分数与碰撞速度和碎片质量的近似函数关系。通过在SPH程序中加入速度矢量子程序,获得了碎片云速度矢量分布。碎片云速度特性分析表明前端速度与碰撞速度成线性关系。分析了弹丸初始半径范围内的碎片云总动量特性,结果表明球形弹丸的动量衰减幅度大于其它两种弹丸。
4、在超高速正碰撞的基础上,进一步对超高速斜碰撞碎片云特性作了数值研究。球形和柱状弹丸的碎片云形状基本为倾斜的椭球,柱状弹丸形成的椭球形碎片云更加倾斜,圆盘形弹丸形成的碎片云类似于斜锥体。得到了超高速斜碰撞碎片云速度矢量分布。给出了单位面积累积动量MA的获取方法,应用该方法得到了三种弹丸MA分布图。累积动量分布表明,在超高速斜碰撞条件下,柱状弹丸和圆盘形弹丸的破坏能力明显强于球形弹丸,且柱状弹丸的破坏能力最强。
5、讨论了物态方程和本构模型对超高速碰撞数值模拟结果的影响。碎片云形状特征的比较表明,Mie-Grüneisen物态方程得到的碎片云图像与实验图像吻合得较好,而Tillotson物态方程得到的碎片分布则较为分散。从弹坑深度和靶孔直径的角度看,Tillotson物态方程的数值结果更接近实验数据,Mie-Grüneisen物态方程得到的结果则明显大于实验结果。本构模型对碎片的整体空间分布影响较小,在局部区域,Johnson-Cook本构模型得到的碎片分布更加符合物理事实。
6、基于弹道极限方程,分别从缓冲屏和后墙间距、缓冲屏密度、缓冲屏厚度、后墙厚度以及后墙材料等五个方面讨论了Whipple结构的设计问题,给出了相应的结构设计方案,数值结果表明结构设计后的防护结构的防护性能得到了较大提高。通过比较不同长径比的非球形弹丸和等效直径球形弹丸超高速正碰撞和斜碰撞Whipple结构的数值模拟结果,讨论了非球形弹丸对Whipple结构防护性能的影响。结果表明,长径比大于1的柱状弹丸对后墙的破坏能力都明显强于圆盘形弹丸和球形弹丸。
7、通过引入弹道极限质量比的概念,定量地比较了非球形弹丸和等效直径球形弹丸破坏能力的大小。超高速斜碰撞分析表明,后墙受到不同长径比的弹丸超高速碰撞后最有可能发生破坏的位置是不同的,而同一长径比、不同等效直径、不同碰撞速度的非球形弹丸对后墙造成破坏的位置基本相同。通过数据拟合,得到了非球形弹丸对后墙造成破坏的位置LBC的经验表达式。
Along with the quick development of aerospace sci-tech and sharply increment of aerospace mission, a great deal of space debris scatter on low earth orbit. This constitutes the complicated space debris environment which has been threatening the safety of the spacecraft seriously. How insure the safety of the spacecraft, which becomes a very important problem of the current spacecraft design. Because of the complexity of hypervelocity impacts, characteristics of debris clouds are not still clear. Therefore the design of spacecraft shielding structure still needs the further research.
Smoothed particle hydrodynamics(SPH) method is the pioneer of meshless particle methods. Nowadays it has been applied to many realms with related fluid mechanics. It has simple logic, can avoid the mesh enlacement and distort of finite element method. So it is suitable for simulating hypervelocity impacts especially.
SPH method is introduced systemically in this paper. And some kinds of problems related with numerical simulations are discussed deeply. The characteristics of debris clouds produced by hypervelocity impacts are studied using SPH method. Hypervelocity impacts of Whipple shields are simulated and analyzed.
The main achievements and conclusions in this paper are presented in the following:
1. Based on Cartesian coordinate system, the 2D and 3D SPH program is composed on one's own. Johnson-Cook tensile cumulative damage constitutive model, Johnson-Cook visco-plastic constitutive model and Steinberg-Cochran-Guinan (SCG) constitutive model are imported in the program successfully. The simulated results obtained by SPH program are in accord with the experiments. This shows that the SPH program is capable of simulating different impact dynamics problems, especially hypervelocity impacts.
2. Five problems related with numerical simulations, such as the computational accuracy, tensile instability, ghost particle method, initial smoothing length and material interface made by different particle spacing, are studied. The three tensile instability treatment methods are compared. The calculated results show that Monaghan method and stable particle methods based on Lagrangian kernels are better than Conservative Smoothing Approach (CSA). The ghost particle method is an effective method to deal with the rigid boundary condition. But a low density region will be produced in the neighborhood of symmetric plane if the ghost particle method is used to deal with the symmetric boundary. This result is unreasonable. The research shows that the simulated results obtain some extent improvement by using the XSPH correction to the velocity of a particle. The effects of initial smoothing lengths on SPH simulated results are analyzed. The results show that the bigger initial smoothing length, the less rigid projectile. And the rational range of initial smoothing length is d0 3. The characteristics of debris clouds obtained by normal hypervelocity impacts, which include sharp, mass, velocity and momentum, are studied. For the spherical projectile, the front of debris clouds shows a hemisphere. And the front of debris clouds has a cone-shaped jag for the rod. The interior of debris clouds produced by the disk will come into being the column structure when its ratio of length to diameter is enough small. A natural boundary judgement approach of the identification of fragments is presented. Combining the neighbour searching method, the mass distribution of debris clouds is achieved. Based on the numerical fit, the cumulative percentage of fragments as a function of the fragment mass and impact velocity is presented. The velocity vector figures of debris clouds are presented through adding the subroutine of velocity vector in SPH program. The analysis of the velocity characteristics shows that the front velocity has a linear change with the impact velocity. The characteristics of the total momentum of debris clouds within initial projectile radius are analyzed. The results show that the attenuation extent of the sphere is larger than that of two other projectiles.
4. Following the analysis of the normal hypervelocity impact, the characteristics of debris clouds obtained by oblique hypervelocity impacts are simulated. The shapes of debris clouds produced by the sphere and the rod are an inclined ellipsoid. However the ellipsoid formed by the rod is more inclined than the sphere. And debris clouds formed by the disk are similar to an inclined centrum. The velocity vector figure of debris clouds of oblique hypervelocity impacts is presented. The acquisition approach of cumulative momentum unit area MA is presented. The MA distributions of three projectiles are obtained by the approach. The results show that the executions of the disk and the rod are more powerful than that of the sphere. And the execution of the rod is strongest.
5. The effects of equations of state and constitutive models on SPH simulated results are discussed. By comparing the shapes of debris clouds, the results obtained by Mie-Grüneisen EOS agree with the experiments. However debris clouds obtained by Tillotson EOS tend to expanding outwards. In the perspective of the depth of the crater and the diameter of the hole of the target, the simulated results from Tillotson EOS are closer to the experimental data than Mie-Grüneisen EOS. The latter results are larger than the experiment clearly. The constitutive models don’t affect the whole distribution of fragments. However debris clouds obtained by Johnson-Cook constitutive model are more reasonable than SCG constitutive model.
6. Based on ballistic limit equations (BLE), the design of Whipple shields is discussed from five aspects which include shield standoff distance, bumper density, bumper thickness, rear wall thickness and rear wall material. The corresponding design project is presented. The results show that the performance of the shielding structure obtains considerable improvements. The effects of non-spherical projectiles with different ratios of length to diameter on shielding performance are studied. Normal and oblique hypervelocity impacts are simulated. The calculated results of non-spherical projectiles are compared with ballistic limit nerves of the sphere. The results show the executions of the rod are stronger than that of the disk and the sphere.
7. By introducing the conception of ballistic limit mass ratio, the executions of non-spherical projectiles are made a qualitative comparison with that of spherical projectiles. The analysis of oblique hypervelocity impacts show that most possible locations on which the rear wall is destroyed relate to the ratio of length to diameter nearly. The projectiles with same ratio of length to diameter have almost same damage position even if they have different equivalent diameter or impact velocity. An empirical expression of the damage position LBC is presented by fitting simulated results.
引文
[1] Lucy L B. A numerical approach to the testing of the fission hypothesis. The Astron J, 1977, 8(12):1013~1024.
[2] Gingold R A, Monaghan J J. Smoothed particle hydrodynamics: theory and applications to non-spherical stars.Mon Not Roy Astrou Soc, 1977, 18:375~389.
[3]黄建国等.等离子驱动超高速微小碎片加速研究.第四届全国空间碎片专题研讨会,南京,2007.
[4] Crews J L, Christiansen E L. The NASA JSC hypervelocity impact test facility(HIT-F).AIAA paper, 92-1640, 1992.
[5] Satoshi Nonaka, Kazuyoshi Takayama, Seishirou Kibe. Hypervelocity impact tests with bumper-walled structure against space debris. AIAA paper, 98-0800, 1998.
[6]柳森,李毅,黄洁,罗锦阳,谢爱民,石安华.用于验证数值仿真的Whipple屏超高速撞击试验结果.宇航学报,2005,26(4):505~508.
[7] Yoshinori Isomoto Oka, Kazuo Nagahashi, Yukimasa Ishii, Yasuo Kobayashi, Toshinori Tsumura. Damge behaviour of metallic materials caused by subsonic to hypervelocity particle impact. Wear 258, 2005:100~106.
[8] Zhou J S, Zhen L, Yang D Z, Li H T. Macro- and microdamage behaviors of the 30CrMnSiA steel impacted by hypervelocity projectiles. Mat Sci and Eng , A282, 2000: 177~182.
[9] Fortov V E, Kim V V, Lomonosov I V, Matveichev A V, Ostrik A V. Numerical modeling of hypervelocity impacts. Int J Impact Engng, 2006, 33: 244~253.
[10] Fortov V E, Lebedev E F, Luzganov S N, Kozlov V V, Medin S A, Parshikov A N, Polistchook V P, Shurupov A V. Railgun experiment and computer simulaion of hypervelocity impact of lexan projectile on aluminum target. Int J Impact Engng, 2006, 33: 253~262.
[11] Bashurov V V, Bebenin G V, Belov G V, Bukharev Y N, Zhukov V I, Ioilev A G, Lapichev N V, Mikhailov A L, Smirnov G S, Fateev Y A, Schlyapnikov G P. Experimental modeling and numerical simulation of high- and hypervelocity space debris impact to spacecraft shield protection. Int J Impact Engng, 1997, 20:69~78.
[12] Orphal D L. Highly oblique impact and penetration of thin targets by steel spheres. Int J Impact Engng, 1999, 23:687~698.
[13] Piekutowski A J. Debris Clouds Produced by the Hypervelocity Impact of Nonspherical Projectiles. Int J of Impact Engng, 2001, 26: 613~624.
[14] Piekutowski A J. Holes Produced in Thin Aluminum Sheets by the Hypervelocity Impact of Aluminum Spheres. Int J of Impact Engng, 1999, 23: 711~722.
[15] Christiansen E L, Kerr J H. Projectile shape effects on shielding performance at 7km/s and 11km/s. Int J Impact Engng, 1997, 20:165~172.
[16] Cour-palais B G. The shape effect of non-spherical projectiles in hypervelocity impacts. Int J Impact Engng, 2001, 26:129~143.
[17] Palmieri D, Faraud M, Destefanis R, Marchetti M. Whipple shield ballistic limit at impact velocities higher than 7 km/s. Int J Impact Engng, 2001, 26:579~590.
[18] Destefanis R, Sch?fer F, Lambert M, Faraud M. Selecting enhanced space debris shields for manned spacecraft. Int J Impact Engng, 2006, 33:219~230.
[19] Schonberg W P. Characterizing secondary debris impact ejecta. Int J Impact Engng, 2001, 26:713~724.
[20] Schonberg W P. Modelling oblique hypervelocity impact phenomena using elementary shock physics. Int J Impact Engng, 1999, 23:823~834.
[21] Schonberg W P. Characterizing material states in orbital debris impacts.SPIE, 2483/31.
[22] Schonberg W P. A first-principles based model characterizing the debris cloud created in a hypervelocity impact. AIAA 94-4487, 1994.
[23] Corvvonato E, Destefanis R, Faraud M. Integral model for the description of the debris cloud structure and impact. Int J Impact Engng, 2001, 21:115~128.
[24] Cohen L J. A debris cloud cratering model. Int J Impact Engng, 1995, 17:229~240.
[25] Fahrenthold E P, Horban B A. An improved hybrid particle-element method for hypervelocity impact simulation. Int J Impact Engng, 2001, 26:169~178.
[26]李长生等.金属塑性加工过程无网格数值模拟方法.沈阳:东北大学出版社,2004.
[27] Groenenboom P H L. Numerical simulation of 2D and 3D hypervelocity impact using the SPH option in PAM-SHOCK. Int J Impact Engng, 1997, 20:309~323.
[28] Hiermaier S, K?nke D, Stilp A J, Thoma K. Computaional simulation of the hypervelocity impact of Al-spheres on thin plates of different materials. Int J Impact Engng, 1997, 20:363~374.
[29] Bashurov V V, Bebenin G V, Ioilev A G. Numerical simulation of rod particles hypervelocity impact effectiveness at various attack angles. Int J Impact Engng, 1997, 20:79~88.
[30] Sekine H, Ito R, Shintate K. A single energy-based parameter to assess protection capability of debris shields. Int J Impact Engng, 2007, 34:958~972.
[31] Carrasco C, Melchor-Lucero O, Osequeda R, Espino L, Fernandez A. Damage-potential comparison of spherical and cylindrical projectiles impacting on a system of bumper plates. Int J Impact Engng, 2006, 33:143~157.
[32] Faraud M, Destefanis R, Palmieri D, Marchetti M. SPH simulations of debris impacts using two different computer codes. Int J Impact Engng, 1999, 23:249~260.
[33] Medina D F, Chen J K. Three-dimensional simulations of impact induced damage in composite structures using the parallelized SPH method. Composites: Part A, 2000, 31:853~860.
[34] Chen J K, Medina D F. The effects of projectile shape on laminated composite perforation. Composites Science and Technology, 1998, 58:1629~1639.
[35] Hertel E S, Mcintosh R L, Patterson B C. A comparison of phase change phenomena in CTH experimental data. Int J Impact Engng, 1995, 17:399~408.
[36] Johnson G R, Beissel S R, Gerlach C A, Stryk R A, Johnson A A, Holmquist T J. Advanced computations for ballistic impact problems. Network Computing services, Inc. Minneapolis, Minnesota, 2004.
[37]闫晓军,张玉珠,聂景旭.空间碎片超高速碰撞数值模拟的SPH方法.北京航空航天大学学报,2005,31(3): 351~354.
[38]贾光辉,黄海,胡震东.超高速撞击数值仿真结果分析.爆炸与冲击,2005,25(1): 47~53.
[39]张伟,庞宝君,贾斌,曲焱喆.弹丸超高速撞击防护屏碎片云数值模拟.高压物理学报,2004,18(1): 47~52.
[40]张伟,马文来,管公顺,庞宝君.非球弹丸超高速撞击航天器防护结构数值模拟.爆炸与冲击,2007,27(3): 240~245.
[41]韩旭,杨刚,强洪夫.光滑粒子流体动力学—一种无网格粒子法.长沙:湖南大学出版社,2005.
[42]张锁春.光滑质点流体动力学(SPH)方法(综述).计算物理,1996,13(4): 385~397.
[43]贝新源,岳宗五.三维SPH程序及其在斜高速碰撞问题的应用.计算物理,1997,14(2): 155~166.
[44] Libersky L D, Petschek A G, Carney T C, Hipp J R, Allahdadi F A. High strain Lagrangian hydrodynamics: a three–dimensional SPH code for dynamic material response. J Comput Phys, 1993, 109:67~75.
[45] Johnson G R, Stryk R A, Beissel S R. SPH for high velocity impact computations. Comput Methods Appl Mech Engrg, 1996, 139:347~373.
[46] Parshikov A N, Medin S A, Loukashenko I I, Milekhin V A. Improvements in SPH method by means of interparticle contact algorithm and analysis of perforation tests at moderate projectile velocities. Int J Impact Engng, 2000, 24:779~796.
[47] Swegle J W. An analysis of smoothed particle hydrodynamics. SAND93-2513, Sandia National Laboratories, Albuquerque, NM 1994.
[48] Randles P W, Libersky L D. Smoothed particle hydrodynamics: some recently improvements and applications. Comput Methods Appl Mech Engrg, 1996, 139:375~408.
[49] Dyka C T, Ingel R P. Addressing tension instability in SPH methods. NRL/MR/6384-94 -7641, Naval Research Laboratory, Washington, 1994.
[50] Seo S, Min O. Axisymmetric SPH simulation of elasto-plastic contact in the low velocity impact. Computer Physics Communications, 2006, 175: 583~603.
[51] Seo S, Min O, Lee J. Application of an improved contact algorithm for penetration analysis in SPH. Int J Impact Engng, 2008, 35: 578~588.
[52] Rabczuk T, Belytschko T, Xiao S P. Stable particle methods based on Lagrangian kernels. Comput Methods Appl Mech Engrg, 2004, 193:1035~1063.
[53] Randles P W, Libersky L D. Recent improvements in SPH modeling of hypervelocity impact. Int J Impact Engng, 1997, 20:525~532.
[54] Krongauz Y, Belytschko T. Consistent pseudo derivatives in meshless methods. Comput Methods Appl Mech Engrg, 1997, 146:371~386.
[55] Bonet J, Kulasegaram S. A simplified approach to enhance the performance of smooth particle hydrodynamics methods. Applied Mathematics and Computation, 2002, 126: 133~155.
[56] Bonet J, Lok T S L. Variational and momentum preservation aspects of smooth particle hydrodynamics formulations. Comput Methods Appl Mech Engrg, 1999, 180:97~115.
[57] Oger G, Doring M, Alessandrini B, Ferrant P. An improved SPH method:towards higher order convergence. J Comput Phys, 2007, 225: 1472~1492.
[58] Monaghan J J. SPH without a tensile instability. J Comput Phys, 2000, 159: 290~311.
[59]丁桦,龙丽平,伍彦峰. SPH方法在模拟线弹性波传播中的应用.计算力学学报,2005,22(3): 320~325.
[60]徐志宏,汤文辉,张若棋.改进的接触算法及其在光滑粒子流体动力学中的应用.国防科技大学学报,2006,28(4): 32~36.
[61]徐志宏.光滑粒子流体动力学方法的改进及其应用:博士论文.长沙:国防科学技术大学,2006.
[62] Benz W, Asphaug E. Simulations of brittle solids using smoothed particle hydrodynamics. Computer Physics Communications, 87:253~265.
[63] Bonet J, Kulasegaram S. Correction and stabilization of smoothed particle hydrodynamics method with applications in metal forming simulations. International Journal for Numerical Methods in Engineering, 2000, 47:1149~1214.
[64] Cleary P W, Prakash M, Ha J. Novel applications of smoothed particle hydrodynamics(SPH) in metal forming. Journal of Materials Processing Technology, 2006, 177:41~48.
[65] Antoci C, Gallati M, Sibilla S. Numerical simulation of fluid-structure interaction by SPH. Computers and Structures, 2007, 85: 879~890.
[66] Beissel S R, Gerlach C A, Johnson G R. Three-dimensional impact simulations by conversion of finite elements to meshfree particles. AIP, 2004.
[67] Hu X Y, Adams N A. An incompressible multi-phase SPH method. J Comput Phys, 2007, 227:264~278.
[68] Jiang F, Oliveira M S A, Sousa A C M. Mesoscale SPH modeling of fluid flow in isotropic porous media. Computer Physics Communications, 2007 , 176:471~480.
[69]宋顺成,才鸿年.弹丸侵彻混凝土的SPH算法.爆炸与冲击,2003,23(1):56~60.
[70]王裴,秦承森,张树道,刘超. SPH方法对金属表面微射流的数值模拟.高压物理学报,2004,18(2):149~155.
[71] Liu M B, Liu G R, Zong Z, Lam K Y. Numerical simulation of underwater explosion by SPH. In Atluri SN & Brust FW(Eds): Advances in computational Engineering & Science, 2000, pp.1475~1480.
[72] Monaghan J J. Why particle methods work(hydrodynamics). SIAM Journal on Scientific and Statistical Computing, 1982, 3:422~433.
[73] Morris J P. A study of the stability properties of SPH. Applied Mathematics Reports and preprints, Monash University, 1994.
[74] Swegle J W, Hicks D L, Attaway S W. Smoothed particle hydrodynamicsstability analysis. J Comput Phys, 1995, 116: 123~134.
[75] Lisio R D, Astronomico I, Grenier E, Pulvirenti M. The convergence of the SPH method. Computers Math Applic, 1998, 35(2): 95~102.
[1]张雄,刘岩.无网格法.北京:清华大学出版社,2004.
[2]韩旭等.光滑粒子流体动力学—一种无网格粒子法.长沙:湖南大学出版社,2005.
[3]张刚明.光滑粒子法及其在冲击动力学中的应用:博士论文.合肥:中国科技大学,2002.
[4] Monaghan J J. Particle methods for hydrodynamics. Computer Physics Report, 1985, 3:71~124.
[5] Morris J P. A study of the stability properties of SPH. Applied Mathematics Reports and Preprints, Monash University, 1994.
[6] Liu W, Jun S, Adee J, Belytschko T. Reproducing kernel particle methods for structural dynamics. Int J Numer Meth Eng, 1995, 38: 1655~1679.
[7] Bonet J, Kulasegaram S. A simplified approach to enhance the performance ofsmooth particle hydrodynamics methods. Applied Mathematics and Computation, 2002, 126: 133~155.
[8] Bonet J, Lok T S L. Variational and momentum preservation aspects of smooth particle hydrodynamics formulations. Comput Methods Appl Mech Engrg, 1999, 180:97~115.
[9] Wang Z Q, Lu Y, Hao H, Chong K. A full coupled numerical analysis approach for buried structures subjected to subsurface blast. Computers and Structures, 2005, 83:339~356.
[10] Monaghan J J. An introduction to SPH. Computer Physics Communications, 1988, 48:89~96.
[11] Liu M B, Liu G R, Lam K Y. Constructing smoothing functions in smoothed particle hydrodynamics with applications. Journal of Computational and Applied Mathematics, 2003, 155: 263~284.
[12] Bonet J, Rodríguez-Paz M X. Hamiltonian formulation of the variable-h SPH equations. J Comput Phys, 2005, 209:541~558.
[13] Benz W. Smoothed particle hydrodynamics: a review. Les: NATO Workshop, 1989.
[14] Hernquist L, Katz N. TREESPH: a unification of SPH with the hierarchical tree method. Astrophys J Suppl Ser, 1989, 70:419~446.
[15] B?rve S, Omang M, Trulsen J. Regularized smoothed particle hydrodynamics with improved multi-resolution handling. J Comput Phys, 2005, 208:345~367.
[16] von Neumann J, Richtmyer R D. A method for the numerical calculation of hydrodynamics shocks. J ournal of Applied Physics, 1950, 21:232~247.
[17] Monaghan J J, Gingold R A. Shock simulation by the particle method of SPH. J Comput Phys, 1983, 52: 374~381.
[18] Balsara D S. von Neumann stability analysis of smoothed particle hydrodynamics—suggestions for optimal algorithms. J Comput Phys, 1995, 121:357~372.
[19] Morris J P, Monaghan J J. A switch to reduce SPH viscosity. J Comput Phys, 1997, 136:41~50.
[20] Owen J M. A tensor artificial viscosity for SPH. J Comput Phys, 2004, 201:601~629.
[21] Noh W F. Errors for calculations of strong shocks using an artificial viscosity and an artificial heat flux. J Comput Phys, 1987, 72:78~120.
[22] Heat conduction with discontinuous conductivity. Applied Mathematics Reports and Preprints, Monash university, 1995.
[23] Johnson G R, Stryk R A, Beissel S R. SPH for high velocity impact computations. Comput Methods Appl Mech Engrg, 1996, 139: 347~373.
[24] Campell J, Vignjevic R, Libersky L. A contact algorithm for smoothed particle hydrodynamics. Comput Methods Appl Mech Engrg, 2000, 184: 49~65.
[25] Libersky L D, Petschek A G, Carney T C, Hipp J R, Allahdadi F A. High strain Lagrangian hydrodynamics: a three–dimensional SPH code for dynamic material response. J Comput Phys, 1993, 109: 67~75.
[1] Liu M B, Liu G R, Lam K Y. Constructing smoothing functions in smoothed particle hydrodynamics with applications. Journal of Computational and Applied Mathematics, 2003(155): 263~284.
[2] Liu W, Jun S, Adee J, Belytschko T. Reproducing kernel particle methods for structural dynamics. Int J Numer Meth Engrg, 1995, 38: 1655~1679.
[3] Bonet J, Lok T S L. Variational and momentum preservation aspects of smooth particle hydrodynamic formulations. Comput Methods Appl Mech Engrg, 1999, 87: 97~115
[4] Randles P W, Libersky L D. Recent improvements in SPH modeling of hypervelocity impact. Int J Impact Engng, 1997, 20:525~532.
[5] Monaghan J J. SPH without a tensile instability. J Comput Phys, 2000, 159: 290~311.
[6] Randles P W, Libersky L D. Smoothed Particle Hydrodynamics: Some recent improvement and applications[J]. Comput Methods Appl Mech Engrg 1996, 139: 375~408.
[7]汤文辉,张若棋.物态方程理论及计算概论(第二版).北京:高等教育出版社, 2008.
[8]宋顺成等.战斗部对混凝土先侵彻后爆轰的数值模拟.兵工学报,2006,27(2): 230~234.
[9] B?rvik T, Hopperstad O S, Berstad T, Langseth M. Perforation of 12mm thick steel plates by 20mm diameter projectiles with flat, hemispherical and conical noses part II: numerical simulations. Int J of Impact Engng, 2002, 27: 37~64.
[10] Chen D, Sarumi M, AL-Hassani S T S, Gan S, Yin Z H. A model for erosion at normal impact.Wear, 1997, 205: 32~39.
[11] AL-Hassani S T S, Chen D and Sarumi M. A simple non-local spallation failure model. Int J of Impact Engng, 1997, 19: 493~501.
[12] Rabczuk T, Belytschko T, Xiao S P. Stable particle methods based on Lagrangian kernels. Comput Methods Appl Mech Engrg, 2004, 193:1035~1063.
[13]龙丽平.无网格光滑粒子法及其在爆炸振动传播中的应用:博士论文.北京:北京航空航天大学,2007.
[14] Johnson G R, Stryk R A, Beissel S R, Holmquist T J. An algorithm to automatically convert distorted finite elements into meshless particles during dynamic deformations. Int J Impact Engng, 2002, 27:997~1013.
[15] Anderson Jr C E, Holer V, Walker J D, Stilp A J. Time-resolved penetration of long rods into steel targets. Int J Impact Engng, 1995, 16: 1~18.
[16] Swegle J W, Attaway S W, Heinstein M W, Mello F J, Hicks D L. An analysis ofsmoothed particle hydrodynamics. Sandia National Lab, Report No.SAND93-2513UC-705, Albuquerque, 1994.
[17] Swegle J W, Hicks D L ,Attaway S W. Smoothed particle hydrodynamics stability analysis. Journal of Computational Physics, 1995, 116:123~134.
[1]张庆明,黄风雷.超高速碰撞动力学引论.北京:科学出版社,2000.
[2] Schonberg W P. Characterizing material states in orbital debris impacts. SPIE, 1995.
[3]郭效忠,贾光辉,黄海.超高速撞击碎片云中大碎片特性SPH模拟分析.第四届全国空间碎片专题研讨会,南京, 2007.
[4] Piekutowski A J. Properties of largest fragment produced by hypervelocity impact of aluminum spheres with thin aluminum sheets, AIAA 92-1588.
[5] Benz W, Asphaug E. Simulations of brittle solids using smoothed particle hydrodynamics. Computer Physics Communications, 1995, 87:253~265.
[6] Petschek A G,Libersky L D. Cylindrical Smoothed Particle Hydrodynamics. J Comput Phys, 1993, 109: 76~80.
[7] Piekutowski A J. Debris Clouds Produced by the Hypervelocity Impact of Nonspherical Projectiles. Int J of Impact Engng, 2001, 26: 613-624.
[8] Petschek A G,Libersky L D. Cylindrical Smoothed Particle Hydrodynamics. J Comput Phys, 1993, 109: 76~80.
[9] Livingstone H G, Verolme K, Hayhurst C J. Predicting the Fragmentation Onset Velocity for Different Metallic Projectiles Using Numerical Simulations. Int J of Impact Engng, 2001, 26: 453~464.
[10] Schonberg W P. Modelling oblique hypervelocity impact phenomena using elementary shock physics. Int J Impact Engng, 1999, 23: 823~834.
[11] Petschek A G, Libersky L D. Cylindrical Smoothed Particle Hydrodynamics. JComput Phys, 1993, 109: 76~80.
[12]汤文辉,张若棋.物态方程理论及计算概论(第二版).北京:高等教育出版社,2008.
[13] Randles P W, Libersky L D. Smoothed Particle Hydrodynamics: Some recent improvement and applications. Comput Methods Appl Mech Engrg, 1996, 139: 375~408.
[14]王裴等. SPH方法对金属表面微射流的数值模拟.高压物理学报,2004,18(2):149~156.
[15] Hiermaier S, Konke D, Stilp A J, Thoma K. Computational Simulation of the Hypervelocity Impact of Al-shperes on Thin Plates of different Materials. Int J of Impact Engng, 1997, 20: 363~374.
[16]彭建祥等.冲击压缩下铝、铜、钨的剪切模量和屈服强度与压力和温度的相关性.物理学报,2005, 54(5): 2194~2197.
[17] Campbell J, Vignjevic R. Development of Lagrangian Hydrocode Modelling for Debris Impact Damage Prediction. Int J of Impact Engng, 1997, 20: 143~152.
[18] Piekutowski A J. Debris Clouds Produced by the Hypervelocity Impact of Nonspherical Projectiles. Int J of Impact Engng, 2001, 26: 613-624.
[1]霍玉华,姜军,满广龙,赵海涛,刘志栋.空间碎片对航天器影响初步分析.第四届全国空间碎片专题研讨会,南京,2007.
[2]龚自正,杨继运,张文兵,童靖宇,向树红,庞贺伟.航天器空间碎片超高速撞击防护的若干问题.航天器环境工程,2007,24(3):125~130.
[3] Kalinski M E. Hypervelocity impact analysis of international space station Whipple and enhanced stuffed Whipple shields. Naval Postgraduate School, Monterey, Master’s thesis, 2004.
[4] Christiansen E L, Crews J L, Williamsen J E, Robinsen J H, Nolen A M. Enhanced meteoroid and orbital debris shielding. Int J Impact Engng, 1995, 17:217~218.
[5] Palmieri D, Faraud M, Destefanis R, Marchetti M. Whipple shield ballistic limit at impact velocities higher than 7km/s. Int J Impact Engng, 2001, 26:579~590.
[6] Chhabildas L C, Reinhart W D, Thornhill T F, Brown J L. Shock-induced vaporization in metals. Int J Impact Engng, 2006, 33:158~168.
[7]韩增尧.空间碎片防护技术研究的若干问题.第四届全国空间碎片专题研讨会,南京,2007.
[8] Christiansen E L, Kerr J H. Projectile shape effects on shielding performance at 7km/s and 11km/s. Int J Impact Engng, 1997, 20:165~172.