Dynamic Response Characteristics of Double-Row Suspended Steel Sheet Piles on Nearshore Slope under the Impact of Pile Hammer
-
摘要: 岸坡锤击管桩施工过程中,为保证岸坡的安全稳定性,准确评估桩锤的动力荷载作用对岸坡支护钢板桩结构的影响非常重要。针对华容煤炭某码头一期工程,采用动力有限元数值模拟方法建立双排悬挂钢板桩结构计算模型;结合现场振动测试验证,分析桩锤振动作用下双排钢板桩的动力响应特性,实现桩锤振动影响下双排钢板桩的安全性评价。结果表明:钢板桩的位移变形与锤击位点相对位置的关系不大,上、下排钢板桩的最大位移分别为2.51和3.14 cm;最大主应力与锤击位点相对位置的关系不大,上、下排钢板桩的最大主应力均小于20 MPa;Mises应力最大值出现在钢板桩与锤击点垂直连线处,上、下排钢板桩的最大Mises应力分别为20.85和25.40 MPa。双排钢板桩的位移变形与应力处于安全控制范围内。Abstract: During the construction process of pile driving on nearshore slope, it is crucial to accurately assess the impact of dynamic loading from the pile hammer on the stability and safety of the slope’s supporting steel sheet pile structure. In the case of the first phase of the dock project at Huarong, a dynamic finite element numerical simulation method was employed to establish a computational model for the double-row suspended steel sheet pile structure. Through field vibration testing verification, the dynamic response characteristics of the double-row steel sheet piles under the influence of pile hammer impact were analyzed, enabling the evaluation of the safety of the double-row steel sheet piles. The research findings revealed that the displacement deformation of the steel sheet piles had little correlation with the impact location of the hammer. The maximum displacement of the lower row of steel sheet piles was 3.14 cm, while the maximum displacement of the upper row was 2.51 cm. The maximum principal stress had no significant relationship with the impact location, and the maximum principal stress of both the upper and lower rows of steel sheet piles was less than 20 MPa. The maximum von Mises stress occurred at the intersection between the steel sheet pile and the vertical line passing through the impact point of the hammer. The maximum von Mises stress for the upper row of steel sheet piles was 20.85 MPa, while for the lower row it was 25.40 MPa. The displacement deformation and stress of the double-row steel sheet piles remained within the safe control range.
-
Key words:
- pile hammer vibration /
- double-row steel sheet piles /
- dynamic response /
- safety control
-
随着全球经济与科技日新月异的发展,越来越多的乘客选择飞机作为出行工具,越来越多的航空器参与到医疗救援、物流运输、环境监测等活动中。但是,飞机在起飞、巡航及降落过程中不可避免地要面临冰雹冲击的威胁。韩登安等[1]在模拟冰雹冲击航空碳纤维增强复合材料时发现,不同冲击能量的冰雹会对复合材料造成不同程度的损伤。鉴于冰雹形成过程的特殊性,相比其他材料,冰雹冲击载荷的加载过程为短时急剧线性加载,表现出强烈的弹脆性材料属性,而卸载过程为一个时间相对较长的非线性卸载,具有明显的流体性质[2-4]。并且,该加载过程中没有平台期,其动态压缩过程受应变率影响显著,因而,冰为典型的应变率相关材料。在低应变率加载下,冰展现出很强的韧性,但随着应变率的增加,脆性属性更加明显[3]。
自然界中的冰雹一般由透明与不透明的冰层相间交替组成[5-6]。对于冰结构的力学性能,国内外进行了大量研究,取得了一定成果。20世纪70年代末,美国科研人员Haynes[7]就对冰雹冲击进行了研究,通过对冰材料开展应变率恒定的准静态压缩实验,探究了温度对压缩强度的影响。之后,National Aeronautics and Space Administra研究员Schulson [8]和Dempsey等[9]进一步探究了温度和晶体结构对冰力学性能和裂纹扩展的影响,通过拉伸、压缩准静态与静态实验,发现冰的压缩强度远大于拉伸强度,当压缩应变率从10−8 s−1变化到10−3 s−1时,压缩强度从0.5 MPa增加至10 MPa,当应变率为10−1 s−1时,压缩强度降低至6 MPa。而Cole[10]指出,在10−2~10−1 s−1的应变率范围内冰的失效应力与应变率无关。随着实验测试技术的快速发展,在20世纪90年代末到21世纪初,冰力学性能的动态实验研究得以进一步展开[3, 11-12]。研究发现,冰为典型的应变率敏感材料,动态压缩性能受应变率的影响很大,而拉伸状态下材料性能受应变率影响较不敏感,且应力强度相对较低[13]。Schulson[14]率先对淡水冰和海水冰进行应变率在10−8~10−3 s−1范围内的动态压缩力学性能测试,首次观察到随应变率的增加冰由韧性向脆性转变的过程。然而,由于测试结果在应变率10−3 s−1时表现出很强的离散性,因此,其实验所能达到的有效应变率受限于此。近年来,Kim等[15]对冰开展了400~2600 s−1应变率范围内的霍普金森杆动态压缩实验,结果显示,随着应变率的提高,压缩强度离散性增强,可能与试件尺寸(直径12.5 mm)有关。目前,我国航空航天发展水平较发达国家还有很大差距,关于冰的动态力学性能研究工作较少。宋振华[16]在前人研究的基础上开展了冰柱在37 mm霍普金森杆作用下的动态力学性能测试,获得了更加真实准确且更高应变率下的强度极限。
然而,自然界中形成的冰雹与固体冰的力学性能有着很大的不同。ASTM F320-21指出天然冰雹比普通的固体透明冰更轻更坚硬,可将质量分数为12%的棉纤维填充至透明冰中以模拟自然界中的冰雹[17]。Swift[6]对透明冰和含棉纤维冰球进行应变率分别为0.052、0.01、0.1、0.4 s−1的压缩试验,发现棉纤维对冰材料承载能力的增强机制与复合材料承载能力的增强机制基本相同。张丽芬等[18]利用低温试验机对不同种类、不同含量的棉纤维冰球进行应变率为0.16 s−1的单轴压缩试验,结果表明添加棉纤维可以提高冰的抗压强度和韧性。
探究含棉纤维冰材料的结构强度、失效模式以及力学性能的影响因素对适航审定具有一定的应用价值,虽然前人对冰雹的力学性能开展了相关研究,但不同应变率下含棉纤维冰材料的力学性能数据尚不完整,需要进一步研究。为此,本研究采用万能试验机对棉纤维质量分数为0、3%、6%和12%的含棉冰试样开展不同应变率下的压缩实验,探究棉纤维质量分数与应变率对冰试样力学性能的影响,并利用应变能相关理论分析其失效模式。
1. 实 验
1.1 试件制备
ASTM F320-21标准为美国材料与实验协会针对航空器上的透明外壳(如挡风玻璃、透镜盖等)的抗冰雹冲击性能研究而制定的标准。其对撞击的冰雹尺寸、填充物材料、填充物质量分数等参数做出了规定,如表1所示。为了获取冰材料的力学性能,采用的冰材料为圆柱试样,其长径比采用文献[17]的试样参数,即长径比为
√3/2 ,并制备了棉纤维质量分数分别为0、3%、6%、12%的圆柱状冰,具体参数如表2所示。表 1 ASTM F320-21标准冰雹参数Table 1. Standardization of hail in ASTM F320-21Hail diameter/mm Hail mass/g Cotton fiber mass fraction/% Cotton fiber mass/g Error/% Temperature/℃ 13 1.0 12 0.14 ±30 −18 25 8.2 12 1.00 ±30 −18 50 66.4 12 8.00 ±30 −18 表 2 人工制备圆柱冰尺寸Table 2. Parameters of artificial iciclesIcicle diameter/mm Icicle height/mm Icicle mass/g Cotton fiber mass fraction/% Error/% Temperature/℃ 18 8 2.0 0, 3, 6, 12 ±30 −18 16 7 1.4 0, 3, 6, 12 ±30 −18 为制备无气泡、无裂纹的高质量冰柱,采用蒸馏水制备试件。将蒸馏水煮沸约15 min,排净水中溶解的气体后,密封放置冷却,至40~50 ℃后,将其缓慢平稳的注入提前放好棉纤维的模具中,使棉纤维质量分数分别为0、3%、6%、12%。实验所用模具由聚甲基丙烯酸甲酯板材(有机玻璃)制成,如图1所示。有机玻璃材料的导热性能较低,能很好地降低冰块结晶速度,使水中的溶解氧排出,并在冷冻过程中使整体充分且均匀地制冷。根据ASTM F320-21标准,棉纤维采用具有良好亲水性的医用脱脂棉,棉纤维平均长度不小于10 mm,如图2所示。图3显示了制备的棉纤维质量分数不同的含棉冰试件,其中,含棉纤维质量分数为零的试件为透明冰。
1.2 实验方法
准静态压缩实验在带有温度箱的万能试验机上进行,如图4所示。冰是温度敏感性材料,且冰柱试件制作温度为−18 ℃,因此,在开展实验前先将温度箱内温度调节至−18 ℃,待温度稳定10 min后再将冰柱放入夹具,然后等温度稳定至−18 ℃后开始压缩实验。通过控制万能试验机压头的速度,对直径18 mm、高度8 mm,以及直径16 mm、高度7 mm的圆柱体冰试件进行应变率分别为10−4、10−3、10−2 s−1的轴向压缩实验。应变率为10−4和10−3 s−1时,试验机的压头速度能快速上升至相应速度,因此,直接调节压头与试件上表面密实接触;当应变率为10−2 s−1时,压头达到相应速度需要一定时间,因此,令上压头下表面距试件上表面有15 mm空隙,使其在接触试件时达到相应速度。
2. 实验结果与讨论
透明冰的离散性较大,因此,每种应变率下保证至少有5组有效数据;含棉冰在每种应变率下保证至少有3组有效数据。
2.1 棉纤维质量分数对冰强度的影响
图5为直径18 mm、高度8 mm的不同含棉纤维质量分数的圆柱冰试件在10−4、10−3和10−2 s−1应变率下的真应力-应变曲线。实验结果表明,在10−4~10−2 s−1应变率范围内,随着棉纤维质量分数的增加,冰柱的强度呈现逐渐增大的趋势。在该应变率范围内,透明冰在10−4 s−1应变率下表现出韧性性质,达到屈服应力后出现应力软化行为,破坏形式为塑性破坏;当应变率为10−3 s−1时,透明冰在达到屈服应力后应力迅速降低,但没有完全失去承载能力,应力反复增减,是因为该应变率下冰破裂后碎片不断脱落,但由于应变率较小,试件并没有完全破裂,此时,透明冰处于韧性向脆性材料转变的阶段;应变率为10−2 s−1时,透明冰表现出脆性材料的特质,加载段近似于直线加载,在达到屈服强度后,直接破裂,失去承载能力,应力瞬间减小,没有塑性阶段。与透明冰不同,压缩应变率在10−4~10−2 s−1范围内时,含棉冰表现出韧性材料的性质,在达到屈服应力后出现应力软化进而塑性破坏,并且,随着棉纤维的增多,应力-应变曲线的下降段逐渐变缓。
图6显示了含棉纤维质量分数不同的各冰柱的屈服强度随应变率的变化。图6显示,在10−4 s−1应变率下,随着棉纤维质量分数的增加,冰柱的屈服强度逐渐增大,因为棉纤维增强了冰晶之间的黏结性,该应变率下透明冰和含棉冰均表现出韧性性质,均经历从开始压缩到逐渐压实的过程,并伴随表面冰碎片掉落;应变率为10−3 s−1时,相比于质量分数为3%的含棉冰,透明冰的屈服强度很小,质量分数为6%的含棉冰的屈服强度提高3.5%,质量分数为12%的含棉冰屈服强度更大,但低于质量分数为6%的含棉冰的屈服强度,因为将冰晶粘结在一起的棉纤维在一定程度上起到加强压缩强度的作用,但当棉纤维质量分数继续增加,冰柱自身的强度降低,且棉纤维的加入会增加冰界面的薄弱层,因此,质量分数为12%的含棉冰的屈服强度较质量分数为6%的含棉冰的屈服强度在一定程度上有所降低。当应变率为10−2 s−1时,随着棉纤维质量分数的增加,冰柱的屈服强度增加,这是因为冰柱含棉纤维质量分数低时,与棉纤维没有粘连的冰掉落,降低了冰柱的屈服强度。这也说明在应变率为10−2 s−1、棉纤维质量分数处于3%~12%时,棉纤维质量分数与冰柱强度呈正相关。
图7为直径18 mm、高8 mm的透明冰在10−4 s−1应变率下单轴压缩时的应力-应变曲线。图8显示了压缩过程中不同时刻的冰柱状态。图8(a)显示了试件在初始时刻的状态;从图8(b)可以看出,试件在压力的作用下出现裂纹,裂纹从中间开始发展,说明试件内部出现了损伤,此时应力-应变曲线中的应力出现骤降,对应图7中的点A;试件继续受压,如图8(c)所示,试件中间的裂纹向下向外扩展至最外侧冰面,裂纹贯穿整个冰柱,对应图7中的点B,AB段曲线的骤降是由裂纹的产生和强度的降低造成的,但此时冰柱仍具有承载力,对其继续加载,强度持续增加;图8(d)显示,试件出现多条纵向扩展的裂纹,并伴随冰块表面有少量碎块掉落;图8(e)显示,有碎片掉落,冰柱向侧面鼓起,说明冰块处于压实过程;图8(f)为压缩完成后冰块的形态,可以看出该应变率下冰块较为完整,周围有掉落的碎片,但并没有向四周飞溅。随着应变率的增大,透明冰在压缩过程中出现越来越多的碎块,并发出崩裂声。其中,在10−2 s−1应变率下,透明冰碎裂较分散,没有压实过程。图9显示了棉纤维质量分数不同的各个冰柱受压后的破坏形态,可以看出,棉纤维冰在压缩破坏时完整性保持良好,随着棉纤维质量分数的增加,冰柱出现“裂而不碎”的现象,破坏形式趋于塑性破坏,说明棉纤维对冰晶的粘结具有重要作用。
2.2 应变率对冰强度的影响
表3为含棉纤维质量分数不同的各个含棉冰柱在不同应变率下的压缩屈服强度。表3显示,在10−4~10−2 s−1应变率范围内,随着应变率的增大,透明冰的压缩屈服强度逐渐降低,屈服强度在10−3 s−1应变率下较10−4 s−1应变率下降低58.20%,在10−2 s−1应变率下较10−3 s−1应变率下降低了50.75%;质量分数为3%、6%和12%的含棉冰的屈服强度均随应变率的升高逐渐增大,与应变率为10−4 s−1时的屈服强度相比,应变率为10−3 s−1时分别提高36.21%、57.51%、22.08%,与应变率为10−3 s−1时的屈服强度相比,应变率为10−2 s−1时分别提高2.84%、9.79%、58.49%。当应变率从10−4 s−1增加到10−3 s−1时,含棉冰的屈服强度增加幅度较大。其中,质量分数为3%的含棉冰在10−3和10−2 s−1应变率下的屈服强度十分相近,这是由于含量较少的棉纤维不能将冰晶全部粘连起来,在应变率较大的压缩过程中,透明冰层的屈服强度较小,会率先掉落冰晶而产生破坏,从而屈服强度没有大幅度提高。图10为含棉纤维质量分数不同的试件在不同应变率下的应力-应变曲线。可以看出,含棉冰达到屈服强度后,随着应变率的增大,曲线下降阶段变陡,说明含棉冰向脆性阶段转变。
表 3 试件在不同应变率下的压缩屈服强度Table 3. Compressive yield strength of specimens at different strain ratesCotton fiber mass fraction/% Strain rates/s−1 Compressive yield strength/MPa Change/% 0 10−4 3.744 10−3 1.563 −58.20 10−2 0.770 −50.75 3 10−4 8.920 10−3 12.150 36.21 10−2 12.495 2.84 6 10−4 8.920 10−3 14.050 57.51 10−2 15.425 9.79 12 10−4 11.200 10−3 13.673 22.08 10−2 21.670 58.49 2.3 压缩能量与破坏形式
能量转换存在于试件发生物理变化的整个过程中,材料的破坏是一种能量驱动下发生的材料失稳现象[19]。在材料破坏过程中,材料的物理性质发生不可逆的变化,对应产生各种形式的能量转换和耗散。研究不同工况下冰柱破坏过程中能量特征与破坏形式之间的联系,更利于反映外载作用下冰柱的破坏机理。图11为不同工况下各个试件从受外载作用到失效整个过程中单位体积所吸收的应变能,即临界应变能密度,其公式为
dWdV=∫εij0σijdεij ,其中,dεij 为试件变形过程中的应变增量,σij 代表试件受压过程中的应力。结果表明,棉质量分数为6%、12%的含棉冰的临界应变能密度均随应变率的升高而增大,棉质量分数为3%的含棉冰的临界应变能密度在应变率为10−3 s−1时达到最大,透明冰的临界应变能密度随着应变率的升高而减小。这种现象的产生是由不同应变率下冰柱破坏过程中能量转换特性决定的。在准静态下,压缩所产生的能量主要转化成弹性能、塑性能以及裂纹表面能。棉纤维质量分数为6%、12%的试件在10−4~10−2 s−1应变率范围内的韧性破坏明显,试件内的破坏以微小裂纹缺陷为主,微小缺陷在发展和相互限制过程中产生塑性变形,能量转化为弹性能、塑性能以及较少的裂纹表面能。棉纤维质量分数为6%、12%的试件随着应变率的升高屈服强度增大,应变能密度增大;棉纤维质量分数为3%的试件在10−4和10−3 s−1应变率下主要呈现与棉纤维质量分数为6%、12%的试件相同的能量转换形式,在10−2 s−1应变率下,试件呈现出一定的脆性破坏趋势,此时,试件内能量除转化为弹性能和塑性能外,还转化为黏结在一起的冰晶的裂纹表面能。透明冰的破坏形貌随应变率的增大呈现从周向膨胀到劈裂的形态,在10−4 s−1应变率时为韧性破坏,能量转化形式与棉纤维质量分数高的试件相同,在10−3和10−2 s−1应变率下,试件向脆性破坏转换,试件几乎未发生塑性变形就直接断裂,能量主要转化为弹性势能和裂纹表面能。从临界能量密度的数值来看,在准静态条件下,转化为裂纹表面能比转化为塑性能所需的能量更少。
2.4 试件尺寸的影响
透明冰的屈服强度受温度、应变率、晶粒尺寸、晶粒方向等多种因素的影响。为探究试样尺寸对冰材料屈服强度测试结果的影响,将直径18 mm、高8 mm的含棉冰试件和直径16 mm、高7 mm的含棉冰试件进行10−4~10−2 s−1应变率的单轴压缩实验,得到两种试件的屈服强度,如表4所示。
表 4 试件尺寸及屈服强度Table 4. Dimensions and compression strengths of specimensStrain rates/s−1 Cotton fiber mass fraction/% Specimen size Mean yield strength/MPa Diameter/mm Height/mm 10−4 0 16 7 4.47 0 18 8 3.74 3 16 7 8.43 3 18 8 8.62 6 16 7 9.70 6 18 8 8.92 12 16 7 12.75 12 18 8 11.53 10−3 0 16 7 1.20 0 18 8 1.56 3 16 7 12.24 3 18 8 12.15 6 16 7 15.40 6 18 8 14.05 12 16 7 16.80 12 18 8 13.80 10−2 0 16 7 0.76 0 18 8 0.77 3 16 7 13.79 3 18 8 12.50 6 16 7 20.00 6 18 8 15.43 12 16 7 20.10 12 18 8 21.67 表4显示,当应变率较小,试件表现出韧性性能且破坏形貌为周向膨胀时,含棉试件中直径为16 mm的试件的屈服强度略大于直径为18 mm的试件,可能是由于端部效应的影响,即当试件受到压缩作用时,试件与压缩圆盘接触的两端受摩擦力影响不能自由地周向膨胀而导致强度增大。当应变率增大,试件表现为脆性破坏且破坏形貌为劈裂时,试件的屈服强度不受尺寸影响。因此,可以得出结论:直径为16与18 mm的试件在受压时,尺寸对压缩强度没有影响,18 mm直径的试件满足实验要求。
3. 结 论
通过万能试验机对棉纤维质量分数不同的冰柱进行压缩试验,分析了10−4~10−2 s−1应变率下冰柱试样的应力-应变曲线,得到如下结论。
(1) 在10−4~10−2 s−1应变率范围内,透明冰在应变率约为10−3 s−1时由韧性材料转化为脆性材料,含棉冰则无明显转化阶段,表现出韧性性质。这是因为棉纤维将冰晶黏结在一起,在压缩的过程中表现出“裂而不碎”的现象,其应力-应变曲线呈现出弹性加载、应力软化、平台区3个阶段。
(2) 在应变率为10−4 与10−2 s−1的压缩载荷下,随着棉纤维质量分数的增加,冰的屈服强度增加;当压缩载荷应变率为10−3 s−1,棉纤维质量分数不大于6%时,随棉纤维质量分数的增加,含棉冰的屈服强度逐渐升高,棉纤维质量分数大于6%时,含棉冰的屈服强度减小。
(3) 随着压缩载荷应变率的增加,透明冰的屈服强度降低,含棉冰的屈服强度呈增大趋势。其中,质量分数为3%的含棉冰在10−3和10−2 s−1应变率下的屈服强度较为接近。
(4) 由临界应变能密度可以看出,在准静态条件下,相比于塑性能,转化为裂纹表面能的能量更少。
-
表 1 现场振动测试数据
Table 1. Field vibration test data
No. Test No. PPV/(cm·s−1) N Distance/m Ⅰ 1 2.000 7 13 2 0.632 3 21 3 0.158 6 38 Ⅱ 1 2.045 7 13 2 0.644 2 21 3 0.160 6 38 Ⅲ 1 2.025 7 13 2 0.639 3 21 3 0.160 3 38 表 2 数值模拟参数
Table 2. Numerical simulation parameters
Material ρ/(kg·m−3) E0/GPa μ A/MPa B/MPa n C σ0/GPa Steel 7830 205 0.30 792 510 0.26 0.014 Clay 1900 11 0.20 0.19, 0.15 Mudstone 2230 15 0.22 0.30 表 3 质点峰值振动速度对比
Table 3. Comparison of PPVs
No. Burst distance/m Test PPV/(cm·s−1) Simulated PPV/(cm·s−1) Relative error/% 1 13 2.023 1.867 7.7 2 21 0.638 0.576 9.7 3 38 0.159 0.205 28.9 -
[1] 吕鹏, 董胜, 焉振, 等. U型钢板桩横向承载性能模型试验研究[J]. 海洋湖沼通报, 2023, 45(1): 1−9.LYU P, DONG S, YAN Z, et al. Experimental study on lateral bearing performance model of U-shaped steel sheet piles [J]. Marine and Lake Sciences Bulletin, 2023, 45(1): 1–9. [2] 焉振, 王元战. 考虑软基不排水强度循环弱化的格型钢板桩防波堤动力有限元分析 [J]. 岩土力学, 2017, 38(5): 1454–1462.YAN Z, WANG Y Z. Dynamic finite element analysis of geogrid steel sheet pile breakwater considering undrained strength degradation of soft foundation [J]. Rock and Soil Mechanics, 2017, 38(5): 1454–1462. [3] 郑国兵, 黄朝煊, 袁文喜, 等. 深厚淤泥中双排钢板桩结构安全稳定性研究 [J]. 水利水电科技进展, 2020, 40(4): 71–76.ZHENG G B, HUANG C X, YUAN W X, et al. Research on the safety and stability of double-row steel sheet pile structure in deep and thick silt [J]. Advances in Science and Technology of Water Resources and Hydropower Engineering, 2020, 40(4): 71–76. [4] 杨熠. 双排钢板桩围堰结构受力变形及稳定性分析 [D]. 南昌: 南昌大学, 2020.YANG Y. Analysis of stress, deformation, and stability of double-row steel sheet pile caisson structure [D]. Nanchang: Nanchang University, 2020. [5] 曹胜敏. 高桩码头桩竖向荷载下静动力学行为研究 [D]. 成都: 西南交通大学, 2008.CAO S M. Study on the static and dynamic behavior of high pile wharf piles under vertical loads [D]. Chengdu: Southwest Jiaotong University, 2008. [6] 王娜娜. 打桩下沉过程中的动力响应研究 [D]. 天津: 天津大学, 2012.WANG N N. Study on the dynamic response during pile driving process [D]. Tianjin: Tianjin University, 2012. [7] XIA Y, JIANG N, ZHOU C, et al. Dynamic behaviors of buried reinforced concrete pipelines with gasketed bell-and-spigot joints subjected to tunnel blasting vibration [J]. Tunnelling and Underground Space Technology, 2021, 118: 104172. doi: 10.1016/j.tust.2021.104172 [8] JIANG N, ZHU B, ZHOU C, et al. Blasting vibration effect on the buried pipeline: a brief overview [J]. Engineering Failure Analysis, 2021, 129: 105709. doi: 10.1016/j.engfailanal.2021.105709 [9] HALLQUIST J O. LS-DYNA theory manual [M]. Livermore, USA: Livermore Software Technology Corporation, 2006. [10] HALLQUIST J O. LS-DYNA keyword user’s manual version 971 [M]. Livermore, USA: Livermore Software Technology Corporation, 2007. [11] ZHOU Y, ZHANG S. Rupture and perforation responses of pressurized tubular members subjected to medium-velocity transverse impact loading [J]. Engineering Failure Analysis, 2021, 127: 105387. doi: 10.1016/j.engfailanal.2021.105387 [12] YAO S, ZHANG D, LU Z, et al. Experimental and numerical investigation on the dynamic response of steel chamber under internal blast [J]. Engineering Structures, 2018, 168: 877–888. doi: 10.1016/j.engstruct.2018.03.067 [13] 邓灵敏. 基于LS-DYNA的动力沉桩全过程分析 [D]. 武汉: 武汉理工大学, 2013.DENG L M. Dynamic analysis of the entire process of pile driving based on LS-DYNA [D]. Wuhan: Wuhan University of Technology, 2013. [14] JIANG N, ZHU B, HE X, et al. Safety assessment of buried pressurized gas pipelines subject to blasting vibrations induced by metro foundation pit excavation [J]. Tunnelling and Underground Space Technology, 2020, 102: 103448. doi: 10.1016/j.tust.2020.103448 [15] 杜闯, 丁红岩, 张浦阳, 等. 钢板桩围堰有限元分析 [J]. 岩土工程学报, 2014, 36(Suppl 2): 159–164.DU C, DING H Y, ZHANG P Y, et al. Finite element analysis of steel sheet pile caisson [J]. Chinese Journal of Geotechnical Engineering, 2014, 36(Suppl 2): 159–164. -