光滑粒子流体动力学法(smoothed particle hydrodynamics, SPH)是一种Lagrange无网格粒子算法，超高速撞击中的强冲击波使材料表现出如流体一样的性质，伴随着材料大变形甚至破碎，传统网格方法难以求解，而基于流体动力学控制方程的SPH方法非常适用[1]。在SPH数值模拟前的处理中，必须设定材料模型，包括状态方程、强度模型和失效模型等。材料状态方程不可或缺，描述材料密度、压力和内能之间的关系；材料强度一般不可忽略，引入强度模型可计算材料偏应力张量和材料当前应变、应变率、温度等对屈服强度的影响。薄板超高速撞击中材料破碎形成碎片云的过程，被认为是稀疏波导致的拉应力超过材料断裂应力后材料不断层裂破坏的过程，而状态方程和强度模型均不能表征材料断裂方面的特征，因此需要考虑在计算中引入失效模型。

利用SPH方法对碎片云进行数值模拟，计算结果与实验吻合很好，但同时发现材料模型的选择对计算结果有影响。Groenenboom[2]发现SESAME状态方程在压力和径向速度方面稍小于多项式状态方程；张伟等[3]发现Tillotson和Puff状态方程相比于Shock状态方程，碎片云径向速度稍大，碎片云颗粒更细小且分散，而强度模型的选择对结果影响不大；Higashide等[4]发现不包含材料相变的状态方程计算得到的碎片最大，考虑材料液化和气化的热力学完全状态方程对相态分布计算较准确，其中考虑了亚稳态相态情况的状态方程计算结果最准确；李宝宝等[5]通过比较GRAY三相物态方程与Tillotson物态方程，也发现材料相变对碎片云外形尺寸影响较大。然而，SPH方法失效模型的选择尚无具体结论，有研究者认为SPH方法中不需要材料失效模型，另一些研究者则采用最大拉应力或Grady失效模型。另外，现有对材料模型方案的研究大多针对碎片云外形或运动状态，未能很好地讨论对碎片质量和数量分布的影响，以及影响程度与撞击工况的相关性，难以实际参考应用。

本文中，利用AUTODYN软件中的SPH模块，针对无失效模型、Grady失效模型和最大拉应力失效模型3种方案，考察失效模型对碎片云外形尺寸及速度、粒子点失效表现和碎片数量分布的影响，讨论失效模型的影响程度随撞击工况的变化，并比较Grady失效模型下不同失效阈值导致的碎片云差异。本文旨在补充SPH方法中材料模型的选择，对碎片云数值模拟特别是碎片分布进行深入分析。

中国画是由内而外生发的艺术，内心的清浊、深浅、静躁、广狭必然影响和决定作品的雅俗和境界高低。所以，历来中国画家都视学识修养的积累为创作的生命之本 ，在“陶冶性灵，变化气质”的过程中提升自己的精神品格，充实自己的内心世界。有大度则华贵，有朴素则真淳，有纯粹则清新，有宁静则淡雅，有潇散则冷逸。写意画如此，工笔画亦然。

1 碎片云几何表现的差异

1.1 计算模型设置

本文中所用算例共2部分，第1部分参照实验情况[6]，取直径为9.53 mm的球形弹丸，材料为2017-T4铝合金，采用Shock状态方程和Johnson-Cook强度模型；薄板材料为6061-T6铝合金，采用Shock状态方程和Steinberg-Guinan强度模型；弹丸正撞击薄板，表1中列出了共3种工况7个算例的薄板厚度、撞击速度及所采用的失效模型方案。

最大拉应力模型即材料拉应力超过设定的失效阈值后材料失效，本文中通过主应力失效实现，设定6061-T6铝合金的失效阈值为2.6 GPa[3]，2017-T4铝合金的失效阈值为2.5 GPa[6]，不考虑最大剪切力失效情况。

表1 各算例的设置 Table 1 Settings of various cases

工况算例薄板厚度/mm撞击速度/(km·s-1)失效模型1010.4676.62无1020.4676.62Grady模型1030.4676.62最大拉应力模型2040.4674.71Grady模型2050.4674.71最大拉应力模型3060.9726.72Grady模型3070.9726.72最大拉应力模型

首先对比有、无失效模型的情况，在算例01和算例02中设置测量点A和B，图3绘制了测量点压力-时间历程。在冲击波作用段内同一测量点曲线完全重合，表明材料被冲击过程与失效模型无关。稀疏波到达一段时间后开始出现小的分歧，进入负压后曲线出现明显分歧。观察右上的放大图可见，算例01中测量点负压均超过6 GPa，并在达到峰值后圆滑回落；算例02中A点负压在到达2.86 GPa后突然跳到0 GPa并保持，B点负压在到达峰值2.12 GPa后圆滑回落，并在后期继续变化。上述差异正是由失效模型导致的，不引入失效模型时粒子点在强拉力下不会失效，从而完整地传递稀疏波，但这和工程材料拉伸断裂的实际表现明显不符；引入失效模型后，拉伸主应力超过失效阈值即会判为失效，粒子不再能承受拉应力或剪切力，如算例02中A点表现所示，负压跳变为0 GPa。

中年级的小学生在朗读的过程中，不但应该学会掌握文本的字面意义，还需要对朗读基调进行把握。要达到这一朗读教学目的，还需要教师带领小学生初步走进文中去体会作者向读者传达的情感，进而引导其清晰地在脑海中勾画文章层次结构与情感脉络。而这也是小学生把握朗读语气及其节奏的有效前提，可使其在朗读时把握好生动灵活的朗读技巧和方法。教师应该在学生朗读前为其进行朗读示范，将学生应掌握的朗读技巧融入范读之中，使其有所感悟地去通过不同朗读基调来对文中的情感进行表达，提高其朗读的水平。

1.2 不同失效模型方案下碎片云几何表现的比较

从材料失效考虑，当材料破碎不充分时，稀疏波强度与失效阈值相当，因此失效阈值上细微的差别就能导致很大的碎片分布差异；材料破碎充分时，稀疏波强度远大于失效阈值，失效模型影响便较弱。另外，Grady模型中失效阈值主要随材料屈服强度改变，撞击较弱时应变、应变率等较小，导致失效阈值计算值偏小。Grady模型下材料本就更易失效，更加偏小的失效阈值就会和最大拉应力模型下恒定的阈值差距更大，加剧两种模型下碎片分布差异。

如图1对照所示，有、无失效模型对碎片云外形影响严重，特别是对弹丸碎片云部分；而Grady失效模型和最大拉应力失效模型的选择并无明显差别。测量碎片云整体尺寸，发现：算例01径向尺寸为算例02的81%，轴向尺寸为算例02的96%；算例03和算例02碎片云整体尺寸相差很小。

图1 算例01～03数值模拟结果 Fig.1 Simulation results of cases 01-03

针对弹丸碎片云部分，图2为算例01～07弹丸碎片主体碎片识别结果放大图，其中碎片识别为只显示包含有未失效粒子的碎片。对比各图可见，算例01中粒子聚集紧密形成一个大碎片，包含了弹丸中94.8%的粒子点；算例02和算例03结果相近，碎片孤立并高度分散，但算例02中碎片更多，且最前端平整，更接近实验图像。算例01中碎片云尺寸明显偏小。分析算例02～07碎片云特征速度[6]发现，最大拉应力模型下碎片扩散程度不及Grady模型，即前端和径向速度偏小但后端速度偏大，不过两者相差不大，都和实验很接近。

综上，不引入失效准则的计算结果明显与实验不符，有失效模型的计算结果均与实验符合，但最大拉应力模型下碎片云的扩散程度稍弱于Grady模型下。

Lesage从1858年起家，存有近60吨的刺绣辅料，有包括从1858年开始的四万本样本存档，每年有几件新纹样的设计。

图2 算例01～07弹丸碎片云碎片识别结果与实验对照[6] Fig.2 Main parts of projectile fragments of cases 01-07 compared with experiments[6]

2 材料失效细观表现的差异

失效模型方案不同导致碎片云计算结果有差异，为分析其内在机理，在计算模型中设置测量点，提取粒子点压力-时间历程作对比分析。

2.1 有、无失效模型计算结果的差异

图3 测量点A、B压力时间历程 Fig.3 Pressure-time curves of gauge points A and B

Grady失效模型同样基于主应力判断失效，即主应力超过失效阈值ps后材料失效，不同的是失效阈值ps非固定值，而由下式计算得到：

酒店管理专业是应用型专业，学生须掌握必需的“专业能力”和“职业能力”。现在的问题是高校教学内容与行业存在脱节，学生知识、能力和素养与酒店的要求有差距[4]。本专业已建设起酒店与校园融合、国内与国外平台融合的“教学、实践一体化”实践教学模式，这一努力和尝试无疑有助于为行业培养更合格更优秀的专业人才，具有一定推广意义。

综上，不引入失效模型也能得到相对孤立的碎片，但只是形式上表现为碎片云，材料表现与实际不符；引入失效模型后粒子点失效表现才符合物理实际，得到的碎片分布符合实验结果。

2.2 最大拉应力及Grady失效模型结果比较

在算例02和算例03中取不同位置处C、D、E等3个测量点，如图5所示为3个点处压力时间历程，比较2种失效模型下粒子点失效表现的差异。C点粒子在Grady和最大拉应力模型下同时失效，但最大拉应力模型下负压更大；D点粒子在最大拉应力模型下负压也更大，而且比Grady模型下更晚失效；E点粒子在Grady模型下失效，但在最大拉应力模型下未失效。

图4 材料失效阈值变化及对应失效规则触发的示例 Fig.4 An instance of material failure regulation and its failure threshold-time curve

进一步考虑，最大拉应力模型中材料失效阈值恒定，但Grady模型中随着材料当前状态的改变，阈值ps的计算值也会变化。图4示例了2种失效模型下某粒子点处失效阈值变化及主应力失效规则的触发，可见Grady模型下失效阈值在冲击波到达后开始剧烈变化，事实上主要随屈服强度变化，密度和体积声速随时间变化不大，而粒子点在最大主应力超过失效阈值后立即被判为失效，应力值归零。

为了给游客带来更好的旅游体验，同时在主营农业的基础上拓展副业，在政府的支持下，村民们利用闲置田地，在景区附近建立农家乐、休闲农庄等场所，给游客提供住宿、餐饮等服务。此外，还包括手工特色农产品零售，如淌(米)面、糍粑等。这些基于乡村旅游的业务提升了村民们的收入，有利于农村的经济增长。

本工程位于陕西省汉中市石马路与傥骆路十字东北侧地块，连通西成高铁汉中站，交通十分便利。本工程总建筑面积155 536m2，地上建筑面积 97 728m2，含东、西翼塔楼（22 层，高100m）和裙房（3～5 层，最高 38m），包括百货、超市、电影院、兼具餐饮功能。地下建筑面积57 808m2，为2层地下室。各单体之间在地上均设缝相互分开独立，2层地下室联成一体，结构模型如图1所示。

由图3可见算例02中B点未失效，但其负压峰值远低于无失效模型的算例01下，可见其他粒子的失效导致的压力变化会影响B点压力，失效准则的引入降低了粒子点负压；无失效模型情况可视为粒子最难失效情况(不失效)，结合图5中最大拉应力模型下粒子点负压大于Grady模型下，由此推断：材料更难失效的失效模型下粒子负压更大。从粒子失效时间可得到相近结论，更难失效的模型下粒子更晚失效甚至不会失效。

综上，从粒子点压力表现来看，在本文所采用失效参数下，最大拉应力模型下材料更难失效。

图5 测量点C、D、E压力时间历程 Fig.5 Pressure-time curves of gauge points C, D and E

3 碎片数量分布差异

3.1 碎片数量分布的整体比较

粒子失效后成为断裂点，孤立碎片的边界必然全是已失效粒子，AUTODYN软件基于此进行碎片识别，获得SPH粒子和碎片的对应关系及碎片特征。本文中利用AUTODYN碎片识别功能，统计数值模拟结果中正向运动的弹丸碎片数量、质量。由于不引入失效模型的方案明显与实验不符，此处只对算例02～07进行比较，分析最大拉应力模型和Grady模型下结果的差异。

式中：ρ为材料密度，由材料连续性方程或核函数更新；Y为屈服强度，由材料强度模型根据材料当前应变、应变率和温度计算更新，变化幅度较大，对ps计算值影响严重；c0为材料体积声速，由材料体积模量和当前密度决定；εc为临界失效应变[7]，材料常数，铝合金材料一般取εc=0.15。可见，在不同区域、不同时刻，ps的取值均可能不同。

统计模拟结果中所有碎片，发现3个工况中Grady模型下碎片总数分别比最大拉应力模型下多1.78%、4.48%、0.42%，碎片越多意味着失效粒子越多，这说明Grady模型下材料更易失效。因此，从碎片总数来看，在本文所采用失效参数下，最大拉应力模型下材料更难失效。

在登月探索中和在爱情中一样，时机是关键。理论上来说，如今对于iSpace公司来说是个证明技术能让“月球谷”成真的好时机。最近，该公司的登月航天器设计图通过了初步设计审查，接受26位专家组成的小组的详细检查。该公司计划在2019年夏天开始组装硬件。他们的月球车设计图正在接受升级和更新，而且是基于iSpace公司为了X奖早已研发和建造起的架构。

统计质量大于0.01 mg的弹丸碎片，得图6所示对数坐标下碎片无量纲质量m/M与累计数量Nc的关系，其中Nc为质量不小于m的碎片累计数量，M为0.01 mg以上碎片总质量，3种工况下依次取683、936、340 mg。需要注意，对数坐标下坐标分布非线性，小质量碎片段差异被压缩。不同模型下曲线变化速度有差异，导致不同质量处碎片累计数量相对大小不一，lg(m/M)=-2.9标记线起不同模型下曲线基本一致，直到lgNc=1.1所标记的大碎片段位置。

从曲线整体来看，2种模型下曲线变化趋势一致，各工况中Grady模型下曲线均接近直线，而最大拉应力模型下曲线波动大，偏离直线。不同工况下材料强度及冲击波强度等不同，Grady模型能够很好地与撞击工况自适应，对应的调整失效阈值而在对数坐标下呈现出更好地线性分布，表现更佳。

3.2 大质量碎片分布差异比较

图6中2种模型下曲线lgNc=1.1的大碎片段差异明显，3种工况中均呈现Grady模型下碎片质量先更大后更小的规律。分析碎片位置发现Grady模型下将最大的若干个大碎片的部分材料分给了稍小的大碎片，使大碎片段质量分布更平缓。最大拉应力模型下粒子更聚集，导致最大的个别碎片偏大，其他大碎片偏小。另外，各工况中均为Grady模型下碎片总数更多，但在分界线处碎片数量基本一致，说明最大拉应力模型下是小碎片数量偏少。这同样是最大拉应力模型下材料难失效导致，小碎片附着在大碎片上未剥离，小碎片数量大幅减小，同时增加大碎片质量。

图6 0.01 mg以上碎片累计数量分布 Fig.6 Distribution of cumulative number of debris above 0.01 mg

表2 最大碎片球等效直径计算值与实验值对比(单位：mm) Table 2 Comparison of simulation and experiment results of equivalent diameters of the largest fragments (unit in mm)

工况Grady模型最大拉应力模型实验值12.822.512.9526.235.926.3631.301.281.45

最大碎片的侵彻性能可代表碎片云的侵彻极限，比较发现：各工况中Grady模型下最大碎片的质量分别比最大拉应力模型下小9.84%、18.5%、19.2%，而轴向速度基本一致。综上可推断：Grady模型下得到的碎片云侵彻极限更小，对后续靶板的损伤更均匀；最大拉应力模型下后续靶板损伤更严重。

测量计算结果中最大碎片在各坐标方向尺寸并计算其球等效直径，如表2所示为最大碎片球等效直径计算值与实验值的对比。虽然最大拉应力模型下最大碎片的质量较大，但由于其粒子聚集度高，碎片等效直径却稍小。不难看到，所有工况中Grady模型下最大碎片等效直径计算值更接近实验值，Grady模型计算结果更符合实验。

4 失效模型影响随撞击工况变化

根据各工况板厚及撞击速度可知，工况1下材料初步破碎，工况3下材料充分破碎，工况2介于前两者之间。观察图6中各工况下两曲线间的差异，材料破碎越充分，曲线间差异越小，失效模型方案对计算结果影响越弱，从碎片总数也能得到此结论。

为分析失效准则对碎片云计算结果的影响，设计了算例01～03作对比，算例01不引入失效模型，算例02采用Grady模型，算例03采用最大拉应力模型。另外，文中碎片云图像及分布统计均基于撞击后20 μs时的计算结果。

撞击速度对稀疏波强度影响显著，工况2撞击速度远低于工况1，而工况3的稍大于工况1。因此，工况2相比于工况1和3，不同失效模型下碎片数量曲线之间差异明显，且最大碎片显著减小。板厚主要改变强稀疏波作用时长和作用区域大小，严重影响小碎片数量，但弹丸材料总量不变，导致工况3下碎片总数显著增大，曲线斜率明显变化，但对曲线间差异影响不明显。

从细观上分析失效模型方案对粒子点失效的影响，在弹丸和薄板模型中轴线遍布测量点，图7为只在一种失效模型下失效而在另一种模型下不失效的测量点的轴向坐标Z，0点到弹丸尾端和薄板背面等距。图7中不同模型下失效表现相反的最远粒子的坐标，表征了失效模型对材料失效表现产生明显影响的起始位置。工况2中表现相反的最远粒子坐标最大，工况3中则坐标最小，正是上述的稀疏波强度和失效阈值的比较。工况2下初始稀疏波强度即和失效阈值相当，边缘位置处失效模型的影响已展现出；工况3下需要稀疏波不断衰弱直至中心附近才能表现出失效模型的影响；工况1下最远粒子坐标与工况3接近而与工况2差距较大，主要是由撞击速度的差异导致的。

值得关注的是，图6中大碎片段碎片质量出现明显分歧，但开始分歧位置均在lgNc=1.1标记线附近，说明撞击工况对大碎片段差异大小影响弱。大碎片多靠近弹丸中心，材料不断层裂破碎同时降低稀疏波强度，导致大碎片位置时稀疏波强度均较低，因此撞击工况对大碎片段差异影响不显著。

考虑失效模型对碎片云宏观表现的影响，当粒子点失效后不再承受拉应力，效果即为断裂点，会阻碍稀疏波的传播，从而导致断裂点两侧粒子运动状态差异较大，碎片高度分散。大量失效粒子会构成断裂面，一个孤立碎片的表面即为失效粒子组成的断裂面，碎片内部粒子表现如同算例02中B点，应力波在碎片内部振荡导致其压力值不为零。反之，不引入失效模型，则不会出现断裂点，通过SPH算法核函数近似，近邻粒子之间无障碍地相互影响，导致粒子点速度矢量比较接近从而集中分布，出现图2中算例01所示聚集现象。但核函数近似的影响域有限大，若2个粒子速度差异较大而距离较远，则难以相互影响，也会产生相对孤立的碎片，形式上表现为碎片云。

云南省11月下旬公布《云南省深化国有企业改革三年行动方案（2018—2020年）》，提出将组建云南省国有股权运营管理有限公司。根据整合重组和混改等情况，逐步将主业处于充分竞争领域的省属企业国有股权注入该公司，打造具有云南特色的“类淡马锡”模式。

图7 模型中轴线上不同失效模型下失效表现不同的粒子点 Fig.7 Particles with different failure performance under Grady and max-tension models

5 失效阈值对碎片云的影响

Grady模型中临界应变常数εc直接影响失效阈值计算值，其对所有材料被建议取值0.15[7]，对于铝合金材料已被验证该取值较为适用。考虑到铜材料实验失效应力为1.0～2.5 GPa，而εc=0.15时Grady模型计算值仅为1.0 GPa[7]，因此本文以OFHC为例，尝试εc=0.15,0.30,0.45和0.60，考察失效阈值对碎片云的影响，同时确认铜材料合适的εc取值。

图4表示不同的断层错动量对管道最大轴向拉伸应变和最大轴向压缩应变的影响。从图4可以看出，当断层错动量较小时，管子的最大轴向拉伸和最大轴向压缩应变曲线变化都较平缓，最大轴向应变值较小，这是因为当断层错动较小时，管子抵抗断层变形和自身塑性形变的能力稍强。随着断层位移量的增加，最大轴向拉应变和最大轴向压应变均在增加，但轴向压缩应变的变化值始终小于轴向拉伸应变的变化值[7]。由此可知在此工况下，管子的破坏形式以受拉破坏为主。由于管子有较小的压应变故应当对管子进行适当的构造措施保护，以防管子受压发生破坏。

第2部分算例参照实验[8]取直径7.72 mm、厚2.36 mm的OFHC圆盘正撞击1 mm厚的6061-T6铝薄板，弹速为6.39 km/s，圆盘轴线与速度方向呈4.6°夹角，采用Shock状态方程和Steinberg-Guinan强度模型。

图8为圆盘中心粒子失效时刻压力的变化曲线，随着εc增大，粒子负压增大且更晚失效。临界应变常数εc越大，失效阈值越大，材料越难失效。比较碎片云特征速度发现：随着εc增大，碎片云最前端速度减小1.3%，最外侧径向速度减小3.5%，失效阈值的增大导致碎片云扩散程度下降。但由于边缘位置稀疏波强度很高，失效阈值对碎片云速度影响非常小。事实上，各εc取值下模拟的碎片云扩散程度均不及实验，因此εc=0.15时速度结果最接近实验，但不同取值时碎片云特征速度变化不超过3%，合适的εc取值需要根据具体算例中其他材料参数和人工黏性等确定。

图8 圆盘中心粒子失效时刻压力的变化 Fig.8 Pressure-time curves of the center particles of the disc projectile

表3 不同临界应变常数下碎片统计信息 Table 3 Detail information related to debris statistics under different critical strains

εc碎片总数最大碎片/mg0.15177 83214.430.30161 61622.230.45152 95631.210.60145 48034.14

图9 碎片云碎片识别结果侧视图 Fig.9 Side view of debris recognition result of debris clouds

图10 0.01 mg以上碎片累计数量分布 Fig.10 Distribution of cumulative number of debris above 0.01 mg

临界应变常数εc不同直接导致材料破碎结果差异，图9为各εc取值下碎片识别结果，明显εc增大后碎片数量增多，大碎片变大。然而从表3可见εc增大碎片总数减少，图10中0.01 mg以上碎片累计数量随εc增大增多，可见失效阈值变大导致0.01 mg以下微小碎片聚集为小碎片。图10中不同失效阈值下曲线趋势一致，但由于微小碎片的聚集，曲线斜率随阈值增大而增大。表3中最大碎片质量也随失效阈值增大而显著增大，但由于Grady模型下失效阈值与相关，随着εc线性增大,碎片质量的增加幅度越来越小。

综上，增大Grady模型下临界应变常数εc直接导致材料更难失效，宏观表现为碎片云扩散程度小幅减弱，碎片总数减少但0.01 mg以上碎片增加，最大碎片质量明显增大，和前文中采用最大拉应力模型导致的材料更难失效表现基本一致。

做好鸭圈建设。鸭圈在鱼凼之上建设，主要铺设材料为木板与竹条，顶部铺设采用的是彩钢瓦，并设置与田水相连通道，按照每只鸭0.1平方米进行鸭圈建设。

6 结 论

比较发现，不引入失效模型时，材料在强负压下不断裂的表现不符合物理实际，弹丸碎片中粒子高度聚集，与实验结果严重不符。因此，碎片云数值模拟中必须引入材料失效模型。有失效模型时，从碎片识别图像、碎片分布曲线及最大碎片尺寸来看，Grady模型计算结果更符合实验。不过失效模型对碎片云的影响程度与撞击工况相关，材料破碎越充分，不同失效模型下碎片分布曲线之间差异越小。

在本文所用失效参数下，最大拉应力模型下材料更难失效，导致碎片云扩散程度小幅减弱，小碎片聚集为大碎片，大碎片质量变大，碎片云侵彻极限提高。增大失效阈值同样导致材料更难失效，亦有上述表现。

参考文献：

[1] LIU G R, LIU M B. 光滑粒子流体动力学:一种无网格粒子法[M].韩旭,杨刚,强洪夫,译.长沙:湖南大学出版社,2005:12-32.

[2] GROENENBOOM P H L. Numerical simulation of 2D and 3D hypervelocity impact using the SPH option in PAM-SHOCKTM[J]. International Journal of Impact Engineering, 1997,20(1/2/3/4/5):309-323. DOI:10.1016/S0734-743X(97)87503-3.

[3] 张伟,庞宝君,贾斌,等.弹丸超高速撞击防护屏碎片云数值模拟[J].高压物理学报,2004,18(1):47-52. DOI:10.3969/j.issn.1000-5773.2004.01.009.

ZHANG Wei, PANG Baojun, JIA Bin, et al. Numerical simulation of debris cloud produced by hypervelocity impact of projectile on bumper[J]. Chinese Journal of High Pressure Physics, 2004,18(1):47-52. DOI:10.3969/j.issn.1000-5773.2004.01.009.

[4] HIGASHIDE M E, KHISHCHENKO K V, LEVASHOV P R. Hypervelocity impact modeling with different equations of state[J]. International Journal of Impact Engineering, 2006,33(1/2/3/4/5/6/7/8/9/10/11/12):625-633. DOI:10.1016/j.ijimpeng.2006.09.078.

[5] 李宝宝,汤文辉,冉宪文.两种物态方程在超高速碰撞数值模拟中的比较[J].航天器环境工程,2010,27(5):576-580. DOI:10.3969/j.issn.1673-1379.2010.05.006.

LI Baobao, TANG Wenhui, RAN Xianwen. Comparison of two kinds of equation of state in numerical simulation of hypervelocity impact[J]. Spacecraft Environment Engineering, 2010,27(5):576-580. DOI:10.3969/j.issn.1673-1379.2010.05.006.

[6] PIEKUTOWSKI A J. Formation and description of debris cloud produced by hypervelocity impact: NASA CR-4707[R]. USA: NASA, 1996.

[7] GRADY D E. The spall strength of condensed matter[J]. Journal of the Mechanics and Physics of Solids, 1988,36(3):353-384. DOI:10.1016/0022-5096(88)90015-4.

[8] PIEKUTOWSKI A J. A simple dynamic model for the formation of debris clouds[J]. International Journal of Impact Engineering, 1990,10(1/2/3/4):453-471. DOI:10.1016/0734-743X(90)90079-B.