A new smoothed particle hydrodynamics method based on the pseudo-fluid model and its application in hypervelocity impact of a projectile on a thin plate
-
摘要: 引入颗粒动力学理论(拟流体模型)建立了适用于超高速碰撞的SPH新方法。将超高速碰撞中处于损伤状态的碎片等效为拟流体,在描述其运动过程中引入了碎片间相互作用和气体相对碎片的作用。采用该方法对球形弹丸超高速碰撞薄板形成碎片云的过程进行了数值模拟,得到了弹坑直径、外泡碎片云和内核碎片云的形状、分布,并与使用传统SPH方法、自适应光滑粒子流体动力学(ASPH)方法的模拟结果进行对比,结果显示:新方法在内核碎片云形状和分布上计算结果更加准确。同时对Whipple屏超高速碰撞问题进行了研究,分析了不同撞击速度下防护屏弹坑尺寸及舱壁损伤特性等特性,计算结果与实验吻合较好且符合Whipple防护结构的典型撞击极限曲线。Abstract: Based on the kinetic theory of granular flow (pseudo-fluid model), a new Smoothed Particle Hydrodynamics (SPH) algorithm suited for hyper velocity collision was presented in this paper. The damaged debris of the hyper velocity impact was equated with the pseudo-fluid and the effects of the debris' interaction and the effects of the gas on the formation process of the debris cloud were investigated. The new SPH algorithm was employed to simulate the 3D hyper velocity impact of an alloy projectile on thin target plates, and the numerical results of crater diameters, the structure and morphology characteristics of the debris cloud and the core debris cloud's shape and distribution were in good agreement with the experimental results. Compared with the simulations of the standard SPH and ASPH, the simulation of the new algorithm is more accurate in the core debris cloud's characteristics. Meanwhile, the hyper velocity impact of the Whipple shield problem was also simulated at different impact velocities. It was found that the crater diameters and the damage characteristics of the rear walls agree well with the experimental results, and that the simulation results are consistent with the typical ballistic limit curve of a Whipple shielding structure.
-
Key words:
- smoothed particle hydrodynamics /
- debris cloud /
- hypervelocity impact /
- Whipple shield
-
混凝土作为土木工程结构中最常用的建筑材料,其复杂的受力形式及力学性能对结构安全有着决定性作用。在城市环境内,火灾通常伴随偶然性爆炸同时发生,从而对建筑结构造成更严重的损伤和破坏[1]。经历高温后的混凝土承载能力缓慢下降,发生不同程度的损伤劣化;在爆炸冲击等高应变率荷载作用下混凝土材料内部发生一系列变化,甚至可能发生脆性破坏,与静态荷载作用下相比,其应力-应变关系更为复杂,强度等随应变率升高明显增大。实际工程中,结构破坏往往与混凝土抗拉性能有关,已有研究表明,混凝土单轴受压破坏过程在本质上是由垂直于压力作用方向的受拉损伤所控制[2]。因此,对火灾(温度效应)和冲击荷载(率效应)联合作用下混凝土拉伸力学性能的研究具有重要的工程和科学意义。
常温下关于混凝土材料动态破坏行为的应变率效应及机制,国内研究者开展了丰硕的试验及模拟工作。王政等[3]从极限面、状态方程以及损伤的定义及其演化等方面深入分析了用于数值模拟冲击问题的几个混凝土本构模型;Zhou等[4]采用细观尺度模型考虑不同的骨料粒径、骨料分布和材料性能分析了混凝土在高应变率下的动态拉伸行为;刘海峰等[5]运用细观力学方法推导了微裂纹对材料弹性模量的弱化作用以及微裂纹的损伤演化方程,建立了混凝土材料在冲击荷载作用下的一维动态本构模型;Lu等[6]研究了应变率为10−4 s−1到102 s−1范围内混凝土材料的抗拉强度的试验方法,对直接动态拉伸试验、动态劈裂试验和剥落试验进行了数值分析;秦川等[7]进行了混凝土SHPB劈拉试验并对其进行了细观数值模拟,观察了高应变率作用下混凝土的动态劈拉破坏模式;宋来忠等[8]进行了干燥与饱和混凝土的动态劈拉试验,研究了加载速率、湿度对混凝土劈拉强度的影响,结果表明:混凝土的劈拉强度随应变率的增大而提高,饱和混凝土的劈拉强度对应变率更敏感;王孝政等[9]对直径为150 mm的圆柱体混凝土试件在不同加载方式下进行劈拉试验,分析了混凝土力学性能随应变率的变化规律。
对于高温混凝土材料静态拉伸破坏行为,研究者亦做了大量卓有成效的工作。Ma等[10]统计了部分文献中无应力立方体试件高温后剩余抗压强度、抗弯强度、劈拉强度和弹性模量的试验结果;郭金纯等[11]通过对混凝土在不同温度及不同恒温时间下的劈拉试验,得到了温度、恒温时间、劈拉强度的耦合关系,建立了高温后混凝土劈拉强度的计算公式;金鑫等[12]对C40高性能混凝土立方体试件的劈拉强度进行试验研究,分析讨论了不同高温作用后混凝土材料劣化的机理以及对劈拉强度的影响,探讨了劈拉强度随温度变化的规律,结果表明:高温后对混凝土的劈拉强度影响最大的因素是最高作用温度。
目前,国内外的研究工作集中于火灾(高温)和爆炸、冲击荷载等独立作用,对于高温混凝土材料的动态力学行为研究极少,对二者联合作用下的“热-力”耦合作用内在机制的认识尚不清晰,缺乏深入探讨。鉴于此,Jin等[13]对高温下混凝土的动态压缩破坏行为及其细观损伤机制进行了分析探讨。本文中,针对混凝土的动态拉伸特性,考虑混凝土内部非均质性,从细观尺度出发,综合考虑高温下各细观组分力学性能的退化行为及动态加载的应变率效应,研究高温与应变率耦合作用下混凝土的劈拉性能及其细观损伤演化规律。
1. 混凝土细观尺度计算模型
为揭示高温下混凝土的动态劈裂拉伸破坏机理,从细观角度出发,考虑混凝土内部非均质性,将混凝土看作由骨料、砂浆基质及两者间界面过渡区组成的三相复合材料,建立细观尺度分析模型,对高温下混凝土动态劈裂拉伸行为进行数值模拟。模拟过程可分为两个步骤:(1)温度场模拟,模拟非均质混凝土的热传导行为,确定温度场空间分布;(2)力学反应模拟,将温度场模拟的输出结果作为动态力学反应的初始条件,根据混凝土的温度场分布,确定不同温度下有效力学参数,分析其在动态加载下的力学行为。数值计算中,关于高温与高应变率(力学作用)实际上是一种单向耦合。先高温,后高应变率作用,也即是,高温对材料的应变率效应存在影响,而应变率效应对温度场不存在影响。
1.1 细观几何模型
假定骨料为圆形颗粒,采用Monte-Carlo方法投递骨料,生成如图1所示的尺寸为150 mm×150 mm的二维随机骨料混凝土试件。其中,绿色区域代表骨料颗粒,包含两种等效粒径:中石(粒径d=30 mm)颗粒数为6,小石(粒径d=12 mm)颗粒数为56;红色区域表示界面过渡区(interfacial transition zone, ITZ),考虑到计算量的影响,界面厚度设为1 mm;灰色区域代表砂浆基质。在热传导模拟中,采用热传导单元来划分网格;力学模拟中,采用四边形等参单元来划分网格,网格划分尺寸为1 mm。
1.2 本构模型与力学参数
1.2.1 热传导行为及热工参数
经典的热传导方程为:
κ∂2T∂xi2+qV=ρc∂T∂t (1) 式中:κ为导热系数;T为温度,qV为体积热源强度;ρ为介质密度,c为比热容,xi为空间坐标,t为时间。
混凝土的配合比、水分含量及骨料类型等影响着混凝土材料的比热容、导热系数、介质密度等热工性能[14]。这里,参考文献[15-16],考虑热工参数的温度相关性,以温度相关的表达式来确定材料在不同温度下的热工参数:
κTκ0=1.026+0.013(T100)+0.0043(T100)2 (2) κaT=7700.84(350+T)+0.7 (3) cT=790+0.9T (4) caT=(1+T)0.075 (5) ρTρ0={1−0.02(T−115)85115∘C<T≤200∘C0.98−0.03(T−200)200200∘C<T≤400∘C0.95−0.07(T−400)800400∘C<T≤1200∘C (6) 式中:T为温度;κT、cT、ρT分别为高温下砂浆及界面的导热系数、比热容及介质密度;κaT、caT分别为高温下骨料的导热系数、比热容;κ0、ρ0为室温下砂浆及界面的导热系数、介质密度。根据已有研究[17-19]及参数假定,室温下(即20 ℃)混凝土各细观组分热工参数取值详见表1。
1.2.2 力学本构关系及参数
(1) 温度退化效应
相关试验研究[20]表明:高温下混凝土强度受强度等级、骨料类型及温升速度等因素的影响,高温下混凝土的抗压强度、抗拉强度及弹性模量等力学性能均随温度升高而显著下降。
砂浆基质与界面过渡区在高温下关于力学性能的试验研究缺乏,没有直接的计算公式,暂假定高温下砂浆基质和界面相力学性能退化规律与混凝土相同,采用图2 中所示的文献中给出的规范[16]给出的混凝土力学性能随温度退化关系,式(7)为拟合公式:
fcT/fc=1.037+2×103T−8×10−7T2 (7) ftT/ft=1−(T−100)/500 (8) ET/E=1.084+1.384×10−3T (9) 式中:T为温度;fcT、ftT、ET分别为高温下混凝土的抗压强度、抗拉强度、弹性模量;fc、ft、E为室温下混凝土的抗压强度、抗拉强度、弹性模量。
(2) 应变率效应
关于高温下混凝土的应变率效应还没有形成较一致的结论,公认的是存在应变率效应,但试验数据量非常有限。近来,Chen等[21]结合其试验结果,提出了普通混凝土在高温作用下动态强度增大因子的表达式。尽管如此,关于混凝土中各细观组分高温时率效应的相关研究则更为少见,目前尚不能给出一个充分完善的关系式。这里为简便起见,粗糙的处理方式为:暂不考虑高温对混凝土细观组分率效应的影响。
弹性模量、泊松比等力学参数相比于抗压强度和抗拉强度,应变率敏感性较弱。因而本文数值模拟中,仅考虑强度的放大效应,通过强度放大系数δ (动态强度/准静态强度)来表示。CEB规范[16]中用来表征混凝土动态抗压强度提高因子(
δc )及修正的拉伸强度增大因子(δt )公式为:δc=fc,imp,k/fcm=(˙εc/˙εc0)0.014˙εc≤30s−1 (10) δc=fc,imp,k/fcm=0.012(˙εc/˙εc0)1/3˙εc>30s−1 (11) δt=ft,imp,k/ftm=(˙εt/˙εt0)0.018˙εt≤10s−1 (12) δt=ft,imp,k/ftm=0.0062(˙εt/˙εt0)1/3˙εt>10s−1 (13) 式中:fc,imp,k为应变率等于
˙εc 时的混凝土动态抗压强度,fcm为应变率˙εc0=30×10−6s−1 时的准静态抗压强度;ft,imp,k为应变率等于˙εt 时的混凝土动态抗拉强度,ftm为应变率˙εt0=1×10−6s−1 时的准静态抗拉强度。(3) 本构关系模型
试验研究表明[22]:常温静载下,骨料一般不会发生破坏;而高温下,骨料强度、峰值应力、弹性模量均会发生不同程度的降低,经历温度越高,降低幅度越大。同时,在高应变率荷载作用下,骨料很可能破坏。鉴于此,本文中暂采用理想弹塑性模型来描述骨料的力学行为,采用文献[23]中不同温度下的骨料材料参数。Lee等[24]提出了可描述混凝土塑性永久变形及损伤特征的塑性损伤本构模型,适用于单调加载、周期性往复加载及低围压动力加载等多种受力形式,得到了众多学者的广泛应用。由于砂浆基质、界面过渡区与混凝土力学行为类似,同文献[13],本文中采用该塑性损伤本构模型来描述其力学行为,该损伤模型的核心是假定混凝土的破坏形式为拉伸断裂与压缩破坏,由各向同性损伤变量来表征其刚度退化:
σ=(1−d)Del0:(ε−εpl) (14) 式中:
Del0 为各向同性的初始无损伤线弹性张量;εpl 为塑性应变张量;d为刚度退化变量,其值介于0~1之间,其中0表示无任何损伤,1则表示完全失效,而混凝土材料的弹性刚度退化由两个独立的损伤变量来描述,即拉伸损伤因子dt和压缩损伤因子dc。材料在单轴拉伸和压缩条件下的应力-应变关系分别为:σt=(1−dt)E(εt−εplt) (15) σc=(1−dc)E(εc−εplc) (16) 对于高温下混凝土材料动态劈裂拉伸行为的数值模拟,这里采用温度退化效应与应变率效应同时耦合的塑性损伤本构模型来描述砂浆基质与界面过渡区的力学行为。数值模拟中采用的室温下(即20 ℃)混凝土各细观组分的压缩/拉伸强度、弹性模量及泊松比等主要力学参数,详见表2。
1.2.3 热传导分析
本文采用图3所示的国际标准化组织(ISO834)[25]推荐的抗火试验及火灾模拟曲线作为环境温度对混凝土试件进行四面受火分析:
T=T0+345lg(8t+1) (14) 式中:T0表示初始环境温度,℃;t表示加热时间,min;T表示t时刻的温度幅值,℃。
1.2.4 动态加载
《普通混凝土力学性能试验方法标准》[26]规定,劈裂拉伸试验时采用钢制弧形垫块施加集中荷载。为防止混凝土发生局部受压破坏,弧形垫块下部可加设胶合板垫条。高温下,混凝土动态劈拉破坏行为模拟的加载和边界条件为:试件底边与垫条等宽度处采用竖向固定约束;底边中点处采用水平向与竖向约束;两侧为自由边界;试件顶部与垫条等宽度处为荷载输入边界,采用速度v加载控制。试件的名义应变率为
˙ε=v/h (其中,h为混凝土试件高度)。2. 数值计算结果及分析
2.1 温度场分布
图4为混凝土试件在四面受火时不同加热时刻的温度场分布,初始温度为室温,即0 min (20 ℃),加热时间分别为15 min (739 ℃)、30 min (842 ℃)和60 min (945 ℃)。总体来讲,随加热时间的延长混凝土试件靠近加热边界处升温明显,并向内部传递。由于考虑了混凝土的非均质性以及骨料的随机分布,不同时刻下,试件内部的温度场分布是不均匀的,将对混凝土材料内部的力学性能产生不同的影响。
2.2 模拟方法验证
杜敏[27]采用液压式万能试验机对不同尺寸的混凝土立方体试件开展了劈裂拉伸试验。为验证本文混凝土劈裂拉伸行为模拟工作的可靠性,与文献[27]中标准立方体试件的试验结果进行对比分析。这里,对室温下的标准立方体混凝土试件进行准静态劈裂拉伸破坏模拟,模拟及试验获得的劈裂拉伸破坏模式如图5所示。由图5可知,本文中数值模拟获得的拉伸损伤模式与试验结果吻合良好,说明细观数值模型能够很好地模拟混凝土在静态加载、常温作用下劈裂拉伸力学破坏行为。
此外,Jin等[13]采用上述细观数值模拟方法,对混凝土在高温下的动态压缩破坏行为进行了细观尺度数值研究,模拟得到的动态压缩破坏模式及压缩应力-应力关系与已有试验结果吻合良好。综上,本文细观数值分析方法可以用来研究高温作用下混凝土材料的动态劈拉破坏行为。
2.3 高温下混凝土动态劈裂拉伸破坏形态
基于上述已验证的细观数值分析方法,对混凝土试件在高温下的劈裂拉伸破坏行为进行模拟分析,工况为:加热60 min,名义应变率为1 s−1。图6(a)和(b)分别表征了混凝土试件细观结构高温下的竖向应力(σ2)和等效塑性应变(εpe)变化过程,给出了时间t=0.5, 1.5, 2.5, 4.0 ms时刻的应力和应变分布。图6(a)中,正值(红色)代表拉应力,负值(蓝色)代表压应力,可以看出,加载点附近区域处于受压状态,加载点连线两侧区域处于均匀的受拉状态;图6(b)中,由于局部受压以及高温后非均质混凝土的力学损伤分布不均匀,微裂纹首先出现在试件的约束和加载处,随加载时间增大,沿着试件顶部、底部相向延伸并向中心线两侧扩展,损伤越来越明显。综合比较应力、应变分布变化过程,当t =1.5 ms时,达到峰值应力,混凝土试件沿界面过渡区及砂浆产生贯通的劈拉裂缝,随后应力逐渐降低,塑性应变持续增大。
对不同加热时间(0, 15, 30, 60 min)后,名义应变率为10−6 (准静态)、10−4、10−2、10−1、1、10和100 s−1的混凝土试件分别进行劈裂拉伸加载模拟,获得相应的峰值应力处的等效塑性应变如图7所示。从图7可以看出,混凝土试件连接加载及约束处产生贯通的劈拉裂缝。常温下,试件沿破坏面损伤分布均匀连续,随名义应变率升高,裂纹向两侧扩展;高温下,由于不均匀的温度场,试件周边温度损伤尤为严重,主要集中在加载点及约束处附近区域,随名义应变率升高,裂纹变宽并向内部延伸。随加热时间增大,损伤由边界向混凝土试件内部传递扩散,裂纹由细长形延伸至带状。常温下及高温下试件损伤都主要出现在砂浆及界面区域,骨料几乎不受影响,图8给出了高温下高应变率下混凝土各细观组分的损伤分布情况。值得注意的是,名义应变率较高时
(˙ε=100s−1) ,由于冲击速度较大,混凝土试件破坏过程急促,出现局部损伤。常温下,塑性变形较小,属脆性破坏,骨料本身发生破坏(图7中圈记处)。而经历高温后,由于骨料变得松软,塑性变形非常大,破坏荷载低,骨料基本没有破坏,裂缝沿着界面延伸开来,此现象与项凯等[28]、郭金纯等[11]的试验结果较为一致。2.4 峰值应力处竖直应力分布
图9给出了经历不同加热时间后的混凝土试件在不同应变率下达到峰值应力时的竖直应力(σ2)分布。可以看出,由于混凝土试件细观结构的非均质性,其内部应力分布是不连续的,并且呈枣核状由加载点与约束处中心连线向两侧传播扩散。横向观察图9可知,同一温度下,应力随名义应变率升高而逐渐增大,混凝土试件内部的竖直应力分布趋向集中均匀连续;同时,随着名义应变率的升高,由于冲击速度增大,局部受压现象也趋于明显。纵向观察图9,同一名义应变率下,由于混凝土材料温度损伤退化严重,竖直应力较常温下大幅降低,且随加热时间呈明显递减趋势;高温下,由于混凝土试件内部温度场分布不均匀,周边温度较高,塑性损伤显著,应力传播受到抑制,随加热时间增长向试件内部及加载点与约束连线两侧延伸扩散。与高温下混凝土试件相比,常温下混凝土试件的应力分布更加均匀连续。
2.5 拉伸应力-应变关系
图10为不同应变率及不同高温加热下混凝土试件的动态劈裂拉伸应力-应变关系,与压缩曲线相比,受拉曲线更陡,脆性更加明显。由图可知,高温对混凝土有明显的劣化作用,相比于常温下,高温下混凝土材料的劈拉应力显著降低,应力-应变曲线趋于平缓,出现塑性流动现象。同一温度下,随名义应变率增大,峰值应力及峰值应变增大,应力-应变曲线形状不变。在高应变率
(˙ε≥1s−1) 下,由于惯性效应占主导作用,峰值应力增大幅度明显提升,当应变率为100 s−1时,峰值应力突增,应力-应变曲线出现跳跃现象。不难发现,无论常温下还是高温下混凝土试件均表现出显著的应变率效应。除此之外,相比于应变率效应,混凝土强度受温度退化效应影响更为显著,此现象与众多研究者所得结果[29-32]一致。2.6 温度损伤残余因子
通过温度损伤残余因子(即高温下拉伸强度与室温下拉伸强度之比)来反应随温度升高混凝土材料性能的退化程度。图11对比了不同应变率下混凝土试件的温度损伤程度。由图可知,不同应变率下的温度损伤曲线斜率接近,随应变率增大至上而下平行排列,说明应变率对混凝土试件的温度退化规律几乎没有影响。但是同一温度下,逐渐增大的冲击速度加剧了其温度损伤程度,尤其在应变率
˙ε=100s−1 时,局部损伤明显(详见2.3、2.4节),致使其温度损伤残余因子大幅降低,具体数值可参见表3。表 3 不同应变率下的混凝土温度损伤残余因子Table 3. Temperature damage residual factors of concrete at different strain rates and different heating times应变率/s−1 温度损伤残余因子 0 min 15 min 30 min 60 min 1×10−6 1 0.410 0.262 0.242 1×10−4 1 0.401 0.263 0.238 1×10−2 1 0.397 0.262 0.235 1×10−1 1 0.392 0.256 0.229 1×100 1 0.376 0.249 0.228 1×101 1 0.372 0.244 0.221 1×102 1 0.317 0.208 0.194 2.7 动态拉伸强度放大系数
采用同一温度下动态拉伸强度与其准静态拉伸强度(本文中取
˙ε=1×10−6s−1 时强度)之比,即动态拉伸强度增大系数(δt)来描述材料强度随应变率升高而增大的现象。图12给出了在不同加热时间(温度)下混凝土δt与应变率间的散点及其进行拟合关系,比较了常温与高温下δt随应变率的变化关系。由图可知,常温下混凝土较高温下混凝土有更强的应变率敏感性,也进一步说明高温和冲击荷载联合作用下,温度退化效应占主导控制,对应变率效应存在抑制作用。高温作用后,同一应变率下,δt随加热时间增大尚未发现明显变化规律。相比动态压缩强度增大系数(δc),δt增长更为迅速,率敏感性更明显,可参见表4具体数值。表 4 高温时不同应变率下的混凝土动态拉伸强度增大系数Table 4. Dynamic increase factors of tensile strength of concrete under different strain rates and different heating times应变率/s−1 δt δc 0 min 15 min 30 min 60 min CEB规范[16] 1×10−6 1.000 1.000 1.000 1.000 1.000 1.000 1×10−4 1.113 1.087 1.118 1.095 1.086 1.017 1×10−2 1.249 1.207 1.250 1.214 1.180 1.085 1×10−1 1.355 1.295 1.322 1.283 1.230 1.120 1×100 1.544 1.416 1.470 1.455 1.282 1.157 1×101 1.836 1.667 1.709 1.674 1.337 1.195 1×102 3.137 2.423 2.487 2.519 2.878 1.793 3. 结 论
考虑混凝土材料的非均质性,将混凝土视为由骨料、砂浆与界面组成的三相复合材料,综合考虑各细观组分的温度退化效应和应变率增强效应,对高温下混凝土动态劈裂拉伸行为进行了细观数值模拟,获得了高温下混凝土动态劈拉破坏过程及破坏模式、竖直应力分布形态、应力-应变关系等响应结果,总结了应变率对混凝土宏观力学行为温度退化作用的影响,分析了高温下宏观动态拉伸强度增大系数(δt)随温度变化规律,得到以下结论:
(1)本文细观数值方法与试验结果吻合良好,能够很好地模拟混凝土的劈裂拉伸行为。
(2)常温下,试件沿破坏面损伤分布均匀连续;高温下,试件损伤集中在加载点及约束处;随名义应变率升高,裂纹变宽并向内部延伸;随加热时间增长,损伤由边界向混凝土试件内部传递扩散,裂纹由细长形延伸至带状。
(3)名义应变率较高时,混凝土试件破坏过程急促,出现局部损伤;常温下属脆性破坏,骨料本身发生破坏;经历高温后,塑性变形非常大,骨料基本没有破坏,裂缝沿着界面延伸开来。
(4)由于混凝土试件细观结构的非均质性,其内部应力分布不连续,呈枣核状由加载点与约束处连线向两侧传播扩散;随名义应变率升高,局部受压现象明显。
(5)常温下混凝土较高温下混凝土有更强的应变率敏感性,温度退化效应对应变率效应存在抑制作用。
(6)需要说明的是,由于缺乏相关试验数据,本文对高温后混凝土的应变率效应暂时进行粗糙处理,今后应对其修正和完善。
-
表 1 本构模型参数
Table 1. Parameters of constitutive model
A/MPa B/MPa n C m Tm/K ˙ε0/s-1 300 426 0.34 0.015 1.0 775 1.0 CV/(J·kg-1·C-1) D1 D2 D3 D4 D5 875 0.12 0.13 -1.5 0.0175 0 表 2 状态方程参数
Table 2. Parameters of state equation
ρ/(kg·m-3) a b AT/GPa BT/GPa αT/GPa βT/GPa E0/(MJ·kg-1) 2790 0.5 1.63 75 65 5 5 5 表 3 超高速碰撞数值模拟结果对比
Table 3. Comparison of high-velocity impact simulation results
表 4 LY12铝合金本构模型参数
Table 4. Parameters of LY12 aluminium alloy of constitutive model
A/MPa B/MPa n C m Tm/K ˙ε0/s-1 369 684 0.73 0.083 1.7 775 0.1 CV/(J·kg-1·C-1) D1 D2 D3 D4 D5 875 0.13 0.13 -1.5 0.011 0 表 5 采用SPH新方法的数值模拟结果与实验结果比较
Table 5. Comparison of the new SPH algorithm and experimental results
编号 弹丸直径/mm 撞击角度/(°) 撞击速度/(km·s-1) 防护屏弹孔直径/mm 中心损伤区域宽度/mm 舱壁损伤情况 实验04-0090 5 0 5.29 11.5 约50 3处剥落,无穿孔 仿真04-0090 5 0 5.29 12.1 54 无剥落、穿孔 实验04-0080 5 0 6.15 12.6 约50 无剥落、穿孔 仿真04-0080 5 0 6.15 12.6 62 无剥落、穿孔 -
[1] Klinkrad H.Collision risk analysis for LEO[J]. Advances in Space Research, 1993, 13(2):105-112. http://www.sciencedirect.com/science/article/pii/0273117793905883 [2] Whipple F L.Meteorites and space travel[J]. Astronomical Journal, 1947, 52(13):132-137. http://cn.bing.com/academic/profile?id=8269eae916575afbac913486ec4da33f&encoded=0&v=paper_preview&mkt=zh-cn [3] Piekutowski A J.Debris clouds produced by the hypervelocity impact of nonspherical projectiles[J]. International Journal of Impact Engineering, 2001, 26(1-10):613-624. doi: 10.1016/S0734-743X(01)00122-1 [4] 管公顺, 庞宝君, 哈跃, 等.铝合金Whipple防护结构高速撞击实验研究[J].爆炸与冲击, 2005, 25(5):461-466. doi: 10.3321/j.issn:1001-1455.2005.05.012Guan Gongshun, Pang Baojun, Ha Yue, et al.Experimantal investigation of high-velocity impact on aluminum alloy Whipple shield[J]. Explosion and Shock Waves, 2005, 25(5):461-466. doi: 10.3321/j.issn:1001-1455.2005.05.012 [5] Swift H F, Bamford R, Chen R.Designing space vehicle shields for meteoroid protection a new analysis[J]. Advances in Space Research, 1983, 2(12):219-234. http://www.sciencedirect.com/science/article/pii/0273117782903118 [6] Lucy L B.A numerical approach to the testing of the fission hypothesis[J]. Astronomical Journal, 1977, 82:1013-1024. doi: 10.1086/112164 [7] Gingold R A, Monaghan J J.Smoothed particle hydrodynamics:theory and application to non-spherical stars[J]. Monthly Notices of the Royal Astronomical Society, 1977, 181(3):375-389. doi: 10.1093/mnras/181.3.375 [8] Libersky L D, Petscheck A G.Smoothed particle hydrodynamics with strength of materials[C]//Trease H, Fritts J, Crowley W.Proceedings of The Next Free Lagrange Conference.Springer Verlag, New York, 1991, 395: 248-257. [9] Johnson G R, Petersen E H, Stryk R A.Incorporation of an SPH option into the EPIC code for a wide range of high velocity impact computations[J]. International Journal of Impact Engineering, 1993, 14(1-4):385-394. doi: 10.1016/0734-743X(93)90036-7 [10] Hiermaier S, Könke D, Stilp A J, et al.Computational simulation of the hypervelocity impact of Al-spheres on thin plates of different materials[J]. International Journal of Impact Engineering, 1997, 20(1-5):363-374. doi: 10.1016/S0734-743X(97)87507-0 [11] 徐金宏, 汤文辉, 罗永.光滑粒子模拟法在超高速碰撞现象中的应用[J].爆炸与冲击, 2006, 26(1):53-58. doi: 10.3321/j.issn:1001-1455.2006.01.009Xu Jinhong, Tang Wenhui, Luo Yong.Applications of the smoothed particle hydrodynamics method to hypervelocity impact simulations[J]. Explosion and Shock Waves, 2006, 26(1):53-58. doi: 10.3321/j.issn:1001-1455.2006.01.009 [12] Zhou C E, Liu G R, Lou K Y.Three-dimensional penetration simulation using smoothed particle hydrodynamics[J]. International Journal of Computational Methods, 2007, 4(4):671-691. doi: 10.1142/S0219876207000972 [13] Liu G R, Liu M B.Smoothed particle hydrodynamics:a meshfree particle method[M]. Singapore:Word Scientific Publishing Co.Pte.Ltd., 2003. [14] 陈福振, 强洪夫, 高巍然.风沙运动问题的SPH-FVM耦合方法数值模拟[J].物理学报, 2014, 63(13):130202. doi: 10.7498/aps.63.130202Chen Fuzhen, Qiang Hongfu, Gao Weiran.Simulation of aeolian sand transport with SPH-FVM coupled method[J]. Acta Physica Sinica, 2014, 63(13):130202. doi: 10.7498/aps.63.130202 [15] 陈福振, 强洪夫, 高巍然.气-粒两相流传热问题的光滑离散颗粒流体动力学方法数值模拟[J].物理学报, 2014, 63(23):230206. doi: 10.7498/aps.63.230206Chen Fuzhen, Qiang Hongfu, Gao Weiran.Numerical simulation of heat transfer in gas-particle two-phase flow with smoothed discrete particle hydrodynamics[J]. Acta Physica Sinica, 2014, 63(23):230206. doi: 10.7498/aps.63.230206 [16] 陈福振, 强洪夫, 高巍然.燃料爆炸抛撒成雾及其爆轰过程的SDPH方法数值模拟研究[J].物理学报, 2015, 64(11):110202. doi: 10.7498/aps.64.110202Chen Fuzhen, Qiang Hongfu, Gao Weiran.Numerical simulation of fuel dispersal into cloud and its combustion and explosion with smmoothed discrete particle hydrodymamics[J]. Acta Physica Sinica, 2015, 64(11):110202. doi: 10.7498/aps.64.110202 [17] Johnson G R, Cook W H.A constitutive model and data for metals subjected to large strains, high strain rates and high temperatures[C]//Proceedings of the Seventh International Symposium on Ballistics.Hague, Netherlands, 1983: 541-547 https://www.scirp.org/(S(351jmbntvnsjt1aadkposzje))/reference/ReferencesPapers.aspx?ReferenceID=1224970 [18] 李春雷.2A12铝合金本构关系实验研究[M].哈尔滨:哈尔滨工业大学, 2006. [19] 费思聪, 孙秦.铝合金的J-C失效参数标定与仿真分析[J].计算机仿真, 2013, 30(9):46-50. doi: 10.3969/j.issn.1006-9348.2013.09.011Fei Sicong, Sun Qin.Investigation on parameters calibration for the J-C failure model of aluminum alloy[J]. Computer Simulation, 2013, 30(9):46-50. doi: 10.3969/j.issn.1006-9348.2013.09.011 [20] 迟润强.弹丸超高速撞击薄板碎片云建模研究[D].哈尔滨: 哈尔滨工业大学, 2010. [21] 柳森, 李毅, 黄洁, 等.用于验证数值模拟仿真的Whipple屏超高速撞击试验结果[J].宇航学报, 2005, 26(4):505-508. doi: 10.3321/j.issn:1000-1328.2005.04.024Liu Sen, Li Yi, Huang Jie, et al.Hypervelocity impact test results of Whipple shield for the validation of numerical simulations[J]Journal of Astronautics, 2005, 26(4):505-508. doi: 10.3321/j.issn:1000-1328.2005.04.024 [22] 侯日立, 周平, 彭建祥冲击波作用下LY12铝合金结构损伤的数值模拟[J].爆炸与冲击, 2012, 32(5):470-474. doi: 10.3969/j.issn.1001-1455.2012.05.004Hou Rili, Zhou Ping, Peng Jianxiang, Numerical simulation of shock damage of LY12 aluminium alloy structure[J]. Explosion and Shock Waves, 2012, 32(5):470-474. doi: 10.3969/j.issn.1001-1455.2012.05.004 [23] 徐金中.基于SPH方法的空间碎片超高速碰撞特性及其防护结构设计研究[D].长沙: 国防科学技术大学, 2008. 期刊类型引用(8)
1. 陈丁,黄文雄,黄丹. 光滑粒子法中的摩擦接触算法及其在含界面土体变形问题中的应用. 岩土力学. 2024(03): 885-894 . 百度学术
2. 叶纪元,杨扬,徐绯,王逸韬,何宇廷. 基于自适应FEM-SPH耦合算法的飞机典型部位破片冲击战伤的数值研究. 爆炸与冲击. 2024(06): 134-145 . 本站查看
3. 廖祜明,黎波,樊江,焦立新,于帅超,林健宇,裴晓阳. 超高速撞击下碎片云的OTM分析. 爆炸与冲击. 2022(10): 50-60 . 本站查看
4. 樊江,袁圆,廖祜明,袁庆浩,陈高翔,黎波. 基于最优运输无网格法的Whipple屏超高速撞击数值模拟. 爆炸与冲击. 2020(07): 97-107 . 本站查看
5. 朱留宪,孙勇,张永盛,武友德,冯颖珊. 基于SPH方法的钛合金切削仿真分析. 机械研究与应用. 2020(04): 1-3 . 百度学术
6. 强洪夫,孙新亚,王广,黄拳章. 钢箱内部爆炸破坏的SPH数值模拟. 爆炸与冲击. 2019(05): 24-32 . 本站查看
7. 牛伟龙,莫蓉,孙惠斌,韩周鹏. 基于光滑粒子流体动力学方法与TANH本构方程的钛合金切屑形态预测. 上海交通大学学报. 2019(05): 624-632 . 百度学术
8. 王小峰,陶钢,闻鹏,任保祥,庞春桥. SPH方法在超高速撞击问题中的应用研究. 兵器装备工程学报. 2019(09): 7-11+54 . 百度学术
其他类型引用(9)
-