Dynamic crack propagation and arrest investigated with a cracked eccentrically-holed flattened disc of rock
-
摘要: 针对平台圆环构型的优点, 提出偏心圆孔单裂纹平台圆盘(cracked eccentrically holed flattened disc, CEHFD), 该试样具有更长的断裂路径。利用霍普金森压杆加载系统, 径向冲击CEHFD试样, 完成Ⅰ型动态断裂实验。砂岩试样表面粘贴应变片和裂纹扩展计, 用于监测裂纹动态起裂、扩展和止裂的全过程。实验表明, 在整个断裂过程中, 裂纹非匀速扩展, 裂纹扩展速度在裂纹起裂后加速上升, 在裂纹止裂前有明显的减速, 与地震时断层的动态破裂全过程完全吻合。采用实验-数值-解析法得到动态应力强度因子, 其时间历程呈现先增大后减小的趋势。根据断裂过程不同时刻, 得到相应的动态起裂韧度、扩展韧度及止裂韧度。在动态断裂全过程中, 动态扩展韧度为速度的函数, 变化趋势与速度一致, 随着时间先增大后减小; 动态起裂韧度大于动态止裂韧度, 止裂韧度随着裂纹最大扩展速度的增大而降低, 并且有较大的离散性。Abstract: Expanding the advantages of flattened ring geometry, a new test specimen geometry with a longer fracture path was advised to investigate the complete dynamic fracture process.This geometry was cracked eccentrically holed flattened disc (CEHFD) which was more conducive to study the unstable dynamic crack propagation and dynamic crack arrest.In order to study the mode Ⅰ (opening mode) dynamic fracture process of a rock, CEHFD specimens with strain gauges and crack propagation gauge glued on the specimen's surface were diametrically impacted by a split Hopkinson pressure bar.The experiment investigated crack initiation, rapid crack propagation, and crack arrest, all in one specimen.In a single complete fracture event, the crack accelerated after initiating and decelerated before arresting.The process monitored by the laboratory is completely consistent with the whole process of the dynamic rupture of faults in earthquake.A hybrid experimental-numerical-analytical method was used to determine the dynamic stress intensity factor that increased with increasing time.The dynamic fracture toughness was obtained with the fracture time information.Crack propagation velocity was first increased and then decreased with the time, and the dynamic propagation toughness had the same variation tendency with time as it was the function of velocity.The dynamic arrest toughness was smaller than the dynamic initiation toughness, and decreased when the maximum crack propagation speed increased.
-
液舱结构在弹丸侵彻下的动态响应问题近年来受到工程界广泛关注。战斗部接触爆炸会产生大量高速破片,对舰船舷侧防护液舱结构的侵彻会造成严重毁伤。针对高速弹体侵彻液舱结构的破损问题已开展了大量研究工作。J.H.McMillen等[1]早在20世纪40年代就针对弹体入水冲击波传播等问题通过实验进行了分析。D.Townsend等[2]则结合实验对高速破片侵彻薄壁水箱结构引起的水锤效应做了细致的研究,并提出了衰减冲击波压力的初步构想。随后,R.J.Disimile等[3]通过实验分析了钢质和铝质球形弹丸穿透背水薄铝板后在水中引起的水锤效应。D.Varas等[4]则从压力变化和管壁变形等方面对球形弹丸侵彻注水方管引起的水锤效应进行了研究,并分析了注水液面的影响。E.Deletombe等[5]比较了密闭和敞开水箱两种情形下水锤效应引起的冲击波压力差异。沈晓乐等[6]、徐双喜等[7]、李典等[8]针对高速破片侵彻液舱问题也开展了大量研究。通过分析可知,目前针对破片侵彻液舱问题的研究大多集中于穿透背水钢板后引起的水锤效应,对于破片穿透后的剩余特性研究较少,且研究对象大多为垂直背水薄板。本文中通过弹道实验,分析破片高速侵彻中厚背水钢板后的剩余特性,从弹体和靶板的破坏模式以及破片穿透背水钢板后的瞬时余速等方面,探讨破片高速侵彻倾斜和垂直背水钢板后的剩余特性的差异,力图揭示中厚背水钢板的抗高速侵彻机理,为后续的理论分析和数值模拟研究提供数据支撑。
1. 实验设计
实验采用15 mm口径滑膛弹道枪,通过火药发射。破片初速采用激光测速系统(1#和2#光幕靶)得到。实验弹丸采用质量3.3 g、边长7.5 mm的立方体破片。破片材料为经淬火处理的45钢。实验靶板为5 mm左右的Q235船用低碳钢。为模拟破片垂直和斜侵彻背水钢板的情形,制作了垂直90°和倾斜60°两种水箱结构,分别如图 1~2所示。
2. 弹道实验结果
采用高速摄影机对破片穿透背水钢板后在水中的运动过程进行了观测,得到背水钢板情形下破片余速即弹丸穿透钢板入水瞬间的速度。表 1中给出了弹道实验结果及主要参数。
表 1 实验结果及主要参数Table 1. Experimental results and main parameters工况 靶板厚度/mm 靶后介质 入射角/(°) 破片初速/(m·s-1) 破片余速/(m·s-1) 破坏情况 1 4.98 水(水箱Ⅰ) 90 1 105.0 229.3 穿透 2 4.96 水(水箱Ⅰ) 90 1 231.2 347.1 穿透 3 4.98 水(水箱Ⅱ) 60 1 058.1 167.8 穿透 4 4.94 水(水箱Ⅱ) 60 1 290.3 310.5 穿透 5 4.92 空气 90 1 029.0 384.0 穿透 6 4.96 空气 60 1 029.9 224.1 穿透 3. 实验结果及分析
3.1 弹体破坏模式分析
图 3给出了实验后收集的破片变形破坏形貌。由图可知,立方体破片不管是侵彻背水钢板还是背空钢板,头部均发生了严重的镦粗变形。对于垂直侵彻情形,破片头部近似呈圆形;而对于60°斜侵彻情形,破片头部虽然也近似呈圆形,但圆形中心向入射速度反方向偏移。进一步观察可知,立方体破片基本上是以面接触的形式侵彻靶板,且能够较好地保持侵彻姿态。
3.2 靶板破坏模式分析
由于钢板背后有液体存在,使得背液钢板的抗侵彻过程及耗能机制与背空钢板[9]差异较大。图 4中给出了中厚背水钢板的典型破坏形貌。由图 4可知,高速侵彻下,中厚背水钢板的穿甲破坏模式主要是剪切冲塞破坏。对于垂直和60°斜侵彻两种侵彻情形,主要区别在于穿孔形状。进一步观察图 5可知,背水钢板在厚度方向存在一定的绝热剪切现象。这说明高速侵彻背水钢板过程中,绝热剪切效应不容忽视。而未考虑绝热剪切效应的弹体侵彻背水钢板的能量耗散理论计算[7]是不符合实际情况的。
3.3 破片在水中运动轨迹分析
破片在水中运动的过程主要包括拖曳和空化两个阶段[3-4]。破片在水中运动过程中,垂直方向则受到水动升力和重力及低压效应的共同作用[10]。
结合图 6~9可知,破片初速较大时,即穿透背水钢板后的瞬时余速较大时,破片在水中的运动轨迹向上偏,如工况2(图 7)和工况4(图 9);相反,则向下偏,如工况3(图 8);若破片速度能使水动升力和重力及低压效应保持平衡,则破片的水中运动轨迹近似成水平,如工况1(图 6)。进一步比较可看出,破片在水中运动轨迹的偏转主要还是受入水初速的影响,钢板倾斜角度对轨迹偏转的影响较小。
结合图 6~7可看出,破片入水初期即产生了明显射流,射流方向与水平方向夹角均为约45°。这与弹体低速入水深闭合所产生的铅垂方向的射流[10]有较大差异。进一步分析图 8~9可看出,虽然破片在穿透倾斜背水钢板后也产生了上下方向的射流,但由于倾斜钢板的影响,向上射流的方向平行于钢板平面,与水平方向成60°夹角;向下射流的方向仍与水平方向成45°夹角。进一步结合图 6~9可知,破片入水初期形成的射流强度与破片入水初速有关,破片入水初速越大,射流越强。
进一步分析图 6~9中空泡的相关特性可知,破片穿透垂直背水钢板后的空泡形状比较对称,而穿透60°倾斜背水钢板后形成的空泡则往下偏。从空泡的大小来看,主要还是取决于破片穿透背水钢板后的余速即破片入水初速。
3.4 破片穿透背水钢板后的余速分析
通过分析高速摄影图片得到破片穿透背水钢板后的瞬时余速。通过拟合求导得到各工况破片穿透背水钢板后的瞬时余速,如表 1所示。由表 1通过计算可知,随着破片初速提高,背水钢板单位面密度吸能值相应提高。破片初速越高,提高值越大。相对于背空钢板而言,相同破片初速下,背水钢板的单位面密度吸能要稍大。
3.5 初始冲击波压力特性分析
大量研究结果表明[2, 5, 8],弹体高速侵彻水介质时形成的冲击波压力峰值较大。因此,本节主要分析初始冲击波压力特性。假设稀疏波以对称于自由液面虚拟镜像波源为起点向外传播,则初始冲击波和稀疏波的传播过程如图 10所示。以初始冲击波产生时刻为起点,得到初始冲击波传播至压力测点的时间t1为
t1=√L2h+B2+(D+h−H)2/cw (1) 式中:Lh为压力测点离撞击点的水平距离,B为压力测点离轨迹线的水平横向距离,D为撞击点离自由液面的高度,h为压力测点离水底高度,H为自由液面高度,cw为水中声速。
稀疏波到达压力测点所需的时间t2为
t2=√L2h+B2+(D+H−h)2/cw (2) 由于H>h,因此t2始终要大于t1,但随着水平距离Lh的增大,t2与t1的差值变小。
考虑到初始冲击波有一定的持续作用时间Δt,若稀疏波在初始冲击波持续作用时间内到达,则需考虑稀疏波对初始冲击波的影响;反之,稀疏波对初始冲击波无影响。由此得到以下判据:
{t2<t1+Δt稀疏波对初始冲击波有影响t2≥t1+Δt稀疏波对初始冲击波无影响 (3) 对于线性可压缩介质而言,波传播具有频率恒定的特性[11]。因此,初始冲击波在水中传播时,距离上的几何衰减仅对冲击波载荷幅值即峰值压力大小产生影响,对载荷的持续作用时间没有影响[4]。因而某一压力测点上,初始冲击波持续作用时间Δt可按下式计算[12]:
Δt=3Rv0w/(2cw) (4) 式中:R为被压缩水介质的区域半径,v0w为弹丸(破片)入水初速。
若压力测点位于水平弹道轨迹上,则稀疏波和初始冲击波传播至该点的时间差ΔT为
ΔT=t2−t1=(√L2h+4D2−Lh)/cw (5) 假设该压力测点的初始冲击波峰值压力与距离的平方成反比,且近似呈三角形脉冲形式,则考虑稀疏波的影响下,该压力测点的比冲量I近似等于:
I≈12P(Lh)ΔT=P(L0)(L0Lh)2√L2h+4D2−Lh2cw (6) 式中:P(L0)为距离撞击点L0(L0 < Lh)处的峰值压力值。
将上式对D求导得到:
∂I∂D=P(L0)(L0Lh)4Dcw√L2h+4D2 (7) 由上式可以看出,撞击点水深D对压力测点比冲量I的影响与该点离撞击点距离Lh的3次方成反比,即稀疏波对初始冲击波的影响程度随距离的3次方倒数下降,说明距离越远,稀疏波的影响越小。
图 11中给出了典型工况(工况2)中压力测点1~2的初始冲击波压力时程曲线。采用式(1)~(4)计算得到:工况2中初始冲击波到达测点2和测点1的时间差约为49 μs,而实验得到两测点峰值时间差约为51 μs,两者吻合较好。结合式(3)的判据可得,对于测点1,稀疏波无影响;而对于测点2,稀疏波有影响。而实际情况是,测点2受到了稀疏波影响,测点1未见明显影响。因此,考虑破片后续运动和空化效应影响,本文中对测点1~2中稀疏波是否影响的预测判别还是比较符合实验结果的。
表 2中给出了各工况下压力测点1~2各自的初始冲击波峰值压力大小。由表 2可知,各工况中测点1和测点2处的初始冲击波峰值压力大小顺序是相同的,即工况2>工况4>工况1>工况3。由于初始冲击波峰值压力大小与破片的入水初速近似呈线性增大[2, 4, 8, 12],因而由峰值压力大小顺序可反推得到破片入水初速即破片穿透背水钢板后的瞬时余速的大小顺序应该相同,这与上一节的分析结果是一致的。
表 2 破片侵彻背水钢板各工况测点1~2的峰值压力测量结果Table 2. Measured peak pressures of test points 1 and 2 for cases of water-bakced steel plates penetrated by fragments工况 破片初速v0/(m·s-1) 入射角/(°) 破片余速vr1/(m·s-1) 峰值压力P1max/MPa 峰值压力P2max/MPa 1 1 105.0 90 286.0 6.184 4.580 2 1 231.2 90 462.7 7.707 6.132 3 1 058.1 60 202.9 4.322 3.894 4 1 290.3 60 385.9 6.432 5.563 4. 结论
通过弹道实验,从弹体和靶板的破坏模式、破片穿透背水钢板后的瞬时余速等方面,分析了破片高速侵彻中厚背水钢板的剩余特性,得出主要结论如下:
(1) 破片高速侵彻中厚背水钢板过程中,头部会发生严重镦粗变形,而镦粗变形主要发生在侵彻背水钢板阶段,但背部水介质的“动支撑”作用不可忽视;
(2) 破片高速侵彻下,中厚背水钢板的主要破坏模式为剪切冲塞,迎弹面出现靶材的反向飞溅,而背弹面穿孔则存在绝热剪切现象;
(3) 破片穿透背水钢板初期,会产生空泡和射流,空泡大小和射流强度主要与破片入水初速有关,而空泡形状和射流方向则受背水钢板的倾斜角度影响;
(4) 破片穿透背水钢板后,在水介质中的运动轨迹会发生偏转,偏转方向主要与破片入水初速有关,破片初速较大时,运动轨迹向上偏,反之则向下偏;
(5) 由于水介质的“动支撑”作用和动能耗散效应,破片穿透背水钢板后的瞬时余速较穿透背空钢板后的余速要低,动能损失要大;
(6) 由于有自由液面的存在,水介质中初始冲击波压力特性需考虑稀疏波的影响,而稀疏波影响的可能性随离撞击点的距离的增大而提高,但影响程度会大大降低;
(7) 相同初速条件下,破片侵彻垂直背水钢板引起的初始冲击波的峰值压力大小较侵彻倾斜背水钢板要大,但峰值压力的提高幅度随破片初速的增大或离撞击点距离的增大均相应减弱。
-
表 1 动态实验数据
Table 1. Dynamic experimental data
试样编号 t0/μs t4/μs tf/μs ta/μs tf′/μs (tf-tf′)/μs CEHFD-1 840.9 916.9 76.0 253.6 1 781.8 1 705.8 CEHFD-2 846.4 928.2 81.8 CEHFD-3 846.3 929.1 82.8 309.8 1 797.0 1 714.2 CEHFD-4 838.0 923.7 85.7 327.5 1 804.2 1 718.5 CEHFD-5 844.1 916.9 72.8 282.8 1 793.0 1 720.2 CEHFD-6 844.0 923.5 79.5 252.6 1 774.6 1 695.1 表 2 砂岩的动态起裂和止裂断裂韧度
Table 2. Dynamic initiation and propagation toughness of sandstone
试样编号 ˙KⅠ/(GPa·m1/2·s-1) KⅠC(d)/(MPa·m1/2) vmax/cR KⅠC(a)/(MPa·m1/2) CEHFD-1 27.4 2.08 0.44 0.72 CEHFD-2 36.1 2.95 0.51 CEHFD-3 32.0 2.65 0.47 0.69 CEHFD-4 41.5 3.56 0.35 1.52 CEHFD-5 24.3 1.77 0.40 1.07 CEHFD-6 26.9 2.14 0.51 0.61 -
[1] BRIAN C, YANG Qingda.In quest of virtual test for structural composite[J].Science, 2006, 314(5802):1102-1107.DOI: 10.1126/science.1131624. [2] 余寿文.断裂力学的历史发展和思考[J].力学与实践, 2015, 39(3):390-394.DOI: 10.6052/1000-0879-15-164.YU Shouwen.The historical development of fracture mechanics and rethinking[J].Mechanics in Engineering, 2015, 39(3):390-394.DOI: 10.6052/1000-0879-15-164. [3] DAS S.Dynamic fracture mechanics in the study of the earthquake rupturing process:Theory and observation[J].Journal of the Mechanics and Physics of Solids, 2003, 51(11/12):1939-1955.DOI: 10.1016/j.jmps.2003.09.025. [4] DAS S.The need to study speed[J].Science, 2007, 317(5840):905-906.DOI: 10.1126/science.1142143. [5] ROBINSON D P, BROUGH C, DAS S.The Mw 7.8, 2001 Kunlunshan earthquake:Extreme rupture speed variability and effect of fault geometry[J].Journal of Geophysical Research, 2006(111):B08303.DOI: 10.1029/2005JB004137. [6] KALTHOFF J F, BEINERT J, WINCKLER S.Experimental analysis of dynamic effects in different crack arrest test specimens[J].ASTM STP 1980, 711:109-127.DOI: 10.1520/STP27443S. [7] KOBAYASHI T, DALLY J W.In fast fracture and crack arrest[J].ASTM STP 1977, 627:257-273. doi: 10.1016-0013-7944(86)90191-8/ [8] SCHARDIN H.Velocity effects in fracture[M].1959: 297-330. [9] RAVI-CHANDAR K, KNAUSS W G.An experimental investigation into dynamic fracture:Ⅰ.Crack initiation and arrest[J].International Journal of Fracture, 1984, 25(4):247-262. doi: 10.1007/BF00963460 [10] FREUND L B.Dynamic fracture mechanics[M].Cambridge:Cambridge University Press, 1990. [11] DALLY J W.Dynamic photo elastic studies of fracture[J].Experimental Mechanics, 1979, 19(10), 349-361. doi: 10.1007/BF02324250 [12] KALTHOFF J F.On some current problems in experimental fracture dynamics[C]//Workshop on Dynamic Fracture.Pasadena: California Institute of Technology, 1983: 11-35. [13] RAVI-CHANDAR K, KNAUSS W G.An experimental investigation into dynamic fracture:Ⅲ.On steady state crack propagation and crack branching[J].International Journal of Fracture, 1984, 26(2):141-154.DOI: 10.1007/BF01157550. [14] FINEBERG J, BOUCHBINDER E.Recent developments in dynamic fracture:some perspectives[J].International Journal of Fracture, 2015, 196(1):33-57.DOI: 10.1007/s10704-015-0038-x. [15] ARAKAWA K, MADA T.Unsteady dynamic crack propagation in a brittle polymer[J].Experimental Mechanics, 2007, 47(5):609-615.DOI: 10.1007/s11340-006-9020-x. [16] SINGH R P, KAVATURU M, SHUKLA A.Initiation, propagation and arrest of an interface crack subjected to controlled stress wave loading[J].International Journal of Fracture, 1997, 83(3):291-304. doi: 10.1023/A:1007358901588 [17] 王自强, 陈少华.高等断裂力学[M].北京:科学出版社, 2009. [18] LEVENT T, CIGDEM K.Effects of geometric factors on mode Ⅰ fracture toughness for modified ring tests[J].International Journal of Rock Mechanics and Mining Sciences, 2012, 51(1):149-161. https://www.deepdyve.com/lp/elsevier/effects-of-geometric-factors-on-mode-i-fracture-toughness-for-modified-xPK0koo9Zh [19] IUNG T, PINEAU A.Dynamic crack propagation and crack arrest investigated with a new specimen geometry:Part Ⅰ:Experimental and numerical calculations[J].Fatigue and Fracture of Engineering Materials and Structures, 1996, 19(11):1357-1367. doi: 10.1111/ffe.1996.19.issue-11 [20] 杨井瑞, 张财贵, 周妍, 等.用CSTBD试样确定砂岩的动态起裂和扩展韧度[J].爆炸与冲击, 2014, 34(3):264-271. http://www.bzycj.cn/CN/abstract/abstract8838.shtmlYANG Jingrui, ZHANG Caigui, ZHOU Yan, et al.Determination of dynamic initiation toughness and propagation toughness of sandstone using CSTBD specimens[J].Explosion and Shock Waves, 2014, 34(3):264-271. http://www.bzycj.cn/CN/abstract/abstract8838.shtml [21] 张盛, 王启智, 谢和平.岩石动态断裂韧度的尺寸效应[J].爆炸与冲击, 2008, 28(6):544-551.DOI: 10.11883/1001-1455(2008)06-0544-08.ZHANG Sheng, WANG Qizhi, XIE Heping.Size effect of rock dynamic fracture toughness[J].Explosion and Shock Waves, 2008, 28(6):544-551.DOI: 10.11883/1001-1455(2008)06-0544-08. [22] 周妍, 张财贵, 杨井瑞, 等.圆孔内单边(或双边)裂纹平台巴西圆盘应力强度因子的全面标定[J].应用数学和力学, 2015, 36(1):16-30. http://d.old.wanfangdata.com.cn/Periodical/yysxhlx201501002ZHOU Yan, ZHANG Caigui, YANG Jingrui, et al.Comprehensive calibration of stress intensity factor for flattened Brazilian disc with holed inner single or double cracks[J].Applied Mathematics and Mechanics, 2015, 36(1):16-30. http://d.old.wanfangdata.com.cn/Periodical/yysxhlx201501002 [23] 周妍, 张财贵, 王启智.用圆孔内单边裂纹平台巴西圆盘和实验-数值-解析法确定砂岩的动态起裂和扩展韧度[J].振动与冲击, 2017, 36(5):37-47. http://d.old.wanfangdata.com.cn/Periodical/zdycj201705007ZHOU Yan, ZHANG Caigui, WANG Qizhi.Determination of dynamic initiation toughness and dynamic propagation toughness of sandstone[J].Journal of Vibration and Shock, 2017, 36(5):37-47. http://d.old.wanfangdata.com.cn/Periodical/zdycj201705007 [24] 李炼, 杨丽萍, 曹富, 等.冲击加载下的砂岩动态断裂全过程的实验和分析[J].煤炭学报, 2016, 41(8):1912-1922. http://d.old.wanfangdata.com.cn/Periodical/mtxb201608006LI Lian, YANG Liping, CAO Fu, et al.Complete dynamic fracture process of sandstone under impact loading:Experiment and analysis[J].Journal of China Coal Society, 2016, 41(8):1912-1922. http://d.old.wanfangdata.com.cn/Periodical/mtxb201608006 [25] 唐志平, 王礼立.SHPB实验的电脑化实验处理系统[J].爆炸与冲击, 1986, 6(4):320-327. http://www.bzycj.cn/CN/abstract/abstract11099.shtmlTANG Zhiping, WANG Lili.A computerized system of data processing used in SHPB experiments[J].Explosion and Shock Waves, 1986, 6(4):320-327. http://www.bzycj.cn/CN/abstract/abstract11099.shtml [26] 李夕兵, 古德生, 赖海辉.冲击荷载下岩石动态应力-应变全图测试中的合理加载波形[J].爆炸与冲击, 1993, 13(2):125-130. http://www.bzycj.cn/CN/abstract/abstract10667.shtmlLI Xibing, GU Desheng, LAI Haihui.On the reasonable loading stress waveforms determined by dynamic stress-strain curves of rocks by SHPB[J].Explosion and Shock Waves, 1993, 13(2):125-130. http://www.bzycj.cn/CN/abstract/abstract10667.shtml [27] GUO W G, LI Y L, LIU Y Y.Analytical and experimental determination of dynamic impact stress intensity factor for 40Cr steel[J].Theoretical and Applied Fracture Mechanics, 1997, 26(1):29-34. doi: 10.1016/S0167-8442(96)00031-6 [28] WEISBROD G, RITTLE D.A method for dynamic fracture toughness determination using short beams[J].International Journal of Fracture, 2000, 104(1):89-103. doi: 10.1023/A:1007673528573 [29] WANG Q Z, YANG J R, ZHANG C G, et al.Determination of dynamic crack initiation and propagation toughness of a rock using a hybrid experimental-numerical approach[J].Journal of Engineering Mechanics, 2016, 142(12):04016097.DOI: 10.1061/(ASCE)EM.1943-7889.0001155. [30] 杨井瑞, 张财贵, 周妍, 等.用SCDC试样测试岩石动态断裂韧度的新方法[J].岩石力学与工程学报, 2015, 34(2):279-292. http://d.old.wanfangdata.com.cn/Periodical/yslxygcxb201502007YANG Jingrui, ZHANG Caigui, ZHOU Yan, et al.A new method for determining dynamic fracture toughness of rock using SCDC specimens[J].Chinese Journal of Rock Mechanics and Engineering, 2015, 34(2):279-292. http://d.old.wanfangdata.com.cn/Periodical/yslxygcxb201502007 [31] CHEN Y M, WILKINS M L.Stress analysis of problems with a three-dimensional, time-dependent computer program[J].International Journal of Fracture, 1976, 12(4):607-617.DOI: 10.1007/BF00034646. [32] CHEN L S, KYANG J H.A modified linear extrapolation formula for determination of stress intensity factors[J].International Journal of Fracture, 1992, 54(1):R3-R8.DOI: 10.1007/BF00040859. [33] RAVI-CHANDAR K K.Dynamic fracture[M].Elsevier, 2004. 期刊类型引用(2)
1. 周磊,朱哲明,董玉清,牛草原,邓帅,姜亚成. 动载荷作用下含倒U型孔洞裂纹试件的动力损伤行为研究. 岩石力学与工程学报. 2021(01): 113-125 . 百度学术
2. 刘宁,胡梦凡,周飞. 基于键基近场动力学理论的单裂纹圆孔板冲击破坏研究. 工程力学. 2020(12): 9-17 . 百度学术
其他类型引用(4)
-