核电站高能管道冲击过程数值模拟研究_彭建.pdf

第 44 卷第 2 期2023 年 2 月哈 尔 滨 工 程 大 学 学 报Journal of Harbin Engineering UniversityVol.44.2Feb.2023核电站高能管道冲击过程数值模拟研究彭建1,2,刘浪1,2,石媛媛1,2(1.中广核工程有限公司 核电安全监控技术与装备国家重点实验室,广东 深圳 518172;2.深圳中广核工程设计有限公司,广东 深圳 518172)摘 要:为了研究管道断裂工况下甩击管道对靶物管道的冲击作用,本文采用有限元软件 ABAQUS 分析了各种结构参数对靶物管道变形的影响。考虑几何非线性、材料非线性和接触非线性建立了数值模型;通过与试验进行对比分析,数值结果与试验结果偏差为 3.5%,验证了模型的有效性;选取核电厂典型材料和结构参数,采用数值模拟方法研究了甩击角度、约束方式、径厚比对靶物管道变形的影响。结果表明:靶物管道的变形包括压缩和弯曲 2 种模式,且压缩和弯曲是相互作用的;甩击角度主要影响靶物管道的压缩变形,对弯曲变形影响较小,靶物管道变形随甩击角度的增加而变大,正交甩击时变形最大;约束方式主要影响靶物管道的弯曲变形,对压缩变形不敏感,且轴向无约束时变形最大;随着径厚比的增加,靶物管道变形越大,冲击截面的直径变化量与应变随时间的变化趋势一致。关键词:高能管道断裂;管道甩动;管道变形;非线性数值模拟;冲击过程;甩击角度;管道约束方式;管道径厚比DOI:10.11990/jheu.202109004网络出版地址:https:/ 文献标志码:A 文章编号:1006-7043(2023)02-0307-07Numerical simulation of the high-energy piping impact process in a nuclear power plantPENG Jian1,2,LIU Lang1,2,SHI Yuanyuan1,2(1.State Key Laboratory of Nuclear Power Safety Monitoring Technology and Equipment,China Nuclear Power Engineering Co.,Ltd.,Shenzhen 518172,China;2.China Nuclear Power Design Company,Co.,Ltd.(Shenzhen),Shenzhen 518172,China)Abstract:To study pipe-to-pipe impact under pipe rupture,the influence of various structural parameters on target pipe deformation is analyzed using the finite element software ABAQUS.First,a numerical model is established that considers geometric nonlinearity,material nonlinearity,and contact nonlinearity.Then,through a comparative analysis with the experiment,the difference between the numerical analysis and the experimental result is only 3.5%,which verifies the validity of the model.Finally,on the basis of the numerical method,the effects of whip-ping angle,restraint mode,and diameter thickness ratio on the target pipe deformation are studied by selecting the typical materials and structural parameters of a nuclear power plant.The results show that target pipe deformation includes two modes:crushing and bending,which interact with each other.The whipping angle mainly affects the compression deformation of the target pipeline and has little effect on the bending deformation.The target pipe de-formation increases with the whipping angle,and the maximum deformation occurs in the orthogonal whipping im-pact.The constraint mode mainly affects the bending deformation of the target pipe,which is insensitive to the crushing deformation,and the maximum deformation occurs under no axial constraint.With an increasing diameter-to-thickness ratio,target pipe deformation increases,and the diameter change of the impact section is consistent with the strain changing with time.Keywords:high-energy pipe rupture;pipe whip;pipe deformation;nonlinear numerical simulation;impact process;whipping angle;restraint mode of piping;diameter thickness ratio of the pipe 收稿日期:2021-09-02.网络出版日期:2022-10-28.基金项目:国家重点研发计划(2018YFB1900200).作者简介:彭建,男,工程师;石媛媛,女,高级工程师.通信作者:石媛媛,E-mail:shiyuanyuan .国际原子能机构(IAEA)和美国国家标准学会(ANSI)在 20 世纪 80 年代对核电站高能管道破裂防护设计开展了大量研究,编制了标准 NS-G-1.111和“ANSI/ANS58.2”2,明确要求核电站设计必须考虑管道甩动效应,提供适当的预防和缓解措施,以保证核安全不受到损害。如果核电厂高能哈 尔 滨 工 程 大 学 学 报第 44 卷管道断裂后,从破口释放出来的流体会对断裂管道产生反作用力,从而引起管道甩动,即管道甩动效应。甩击管道撞击附近的其他管道,它会产生管道对管道的冲击,这可能造成更严重的多米诺效应。管道甩击管道可分为 2 个阶段,即管道自由甩动和管道冲击管道阶段。Reid 等3-5国外学者对管道自由甩动进行了大量的实验研究,设置了由储气罐、试验管道和爆破片组成的试验装置,创建了大挠度动力梁理论模型,并通过试验论证了模型的有效性。Alzheimer 等6国外学者对管道冲击管道也开展了实验研究,首先基于美国某核电厂系统布置模型和压力温度情况,识别出 100 组潜在的管道甩击管道工况;然后采用弹射器装置模拟管道冲击管道试验;最后结论为甩击管道不会引起名义尺寸较大且壁厚也较大的靶物管道发生破裂,可能会引起名义尺寸较小的靶物管道(不论其壁厚如何)发生破裂,会造成名义尺寸相同或较大而壁厚相等或较薄的管道发生贯穿裂纹。上述结论广泛应用于工程设计,但变形机理有待进一步研究。Baum 等7专家学者采用试验方法模拟了管道自由甩动和管道冲击管道的全过程,研究了其变形机理,结果发现靶物管道存在压缩和弯曲两种机理,支撑方式对靶物管道的变形影响较大。但是,由于试验基于 1 寸的小直径管道和理想刚性约束,因此对于大管径和柔性管道系统,还有待进一步研究。国内标准 EJ/T 335 和 NB/T 20516 也提出核电站设计必须考虑管道甩击效应。学者对管道甩击也开展了研究,Yang 等8采用试验和数值模拟的方法研究了冲击位置对靶物管道的变形模式的影响。但是,国内关于管道甩击的研究主要集中在防甩限制件方面,王春霖等9-12研究了防甩限制件的分析方法及影响因素。而管道抗冲击性能的研究主要集中在海洋管道13-20。杨秀娟等13采用非线性动态有限元法模拟了坠物撞击海底管道的过程,分析了物体形状、撞击角度、摩擦、内压、混凝土厚度、埋深及悬空长度与管道凹陷深度的关系。谢丽媛等14利用落锤实验和有限元模型得出海底管道在冲击载荷作用下的破坏模式为整体弯曲变形与冲击凹痕部位的局部弯曲耦合形成,管道钢材的屈服强度的增加可有效减小冲击作用下的局部凹陷变形。张荣等15通过落锤实验证明圆钢管失效形式均为剪切失效,随冲击能量增加圆钢管变形依次为弹性小变形阶段、塑性大变形阶段、剪切失效阶段等 3 个阶段。甩击管道和靶物管道都可能产生较大的变形,并受到很多因素的影响,其中管道的管径比、径厚比、冲击能量、冲击位置、支撑方式已进行了研究。但是,冲击角度的影响暂未见文献报道。工程上,对于甩击管道一般没有要求,主要关注靶物管道的变形。因此,本文采用数值模拟方法,通过模拟管道自由甩动和冲击管道全过程,主要研究冲击角度对靶物管道变形的影响,并分析了约束方式和径厚比的影响。1 管道冲击数值模型1.1 数值算法 ABAQUS 是一套功能强大的工程模拟有限元软件,包括 ABAQUS/Standard 和 ABAQUS/Explicit 分析模块。ABAQUS/Explicit 在处理冲击问题方面非常成熟和方便,用户仅需要提供合理的结构参数、材料参数、边界条件和载荷工况这些工程数据就可以获得精确解。ABAQUS/Explicit 采用中心差分法进行显式时间积分:Mu(t+t)+Cu(t+t)+Ku(t+t)=F(t+t)(1)式中:M 为质量矩阵;C 为阻尼矩阵;K 为刚度矩阵;u(t)为位移矩阵;F(t)为载荷矩阵;t 为时间增量。显式时间积分的求解精度取决于时间积分步长,需要非常小的时间步以保持稳定,而时间步取决于单元尺寸和波速;有稳定的时间增量的限制,完成指定分析一般需要较多的时间增量,不需要求解线性方程组,每个增量步的计算成本较低,适用于冲击等高频动力学问题,能分析各种复杂的接触问题。本文采用 ABAQUS/Explicit 模拟管道甩击管道,通过体积粘性考虑模型阻尼,体积粘性采用软件推荐值,线性参数取 0.06,二次参数取 1.2,最小稳定时间步长和接触算法通过验证分析确定。1.2 材料模型 为描述金属材料的高速冲击效应,利用商业有限元软件 ABAQUS 建立管道甩击管道的数值模型。管道采用壳单元 S4 四边形单元进行网格划分,选择Simpson 积分方法,厚度方向取 5 个积分点,单元应力位置在下表面;材料选用 Mises 弹塑性模型,真应力应变曲线根据 ASME BPVC VIII 的附录,并基于屈服强度、抗拉强度、弹性模型计算获得,本文采用核电站中的典型碳钢材料 P265,其真应力应变曲线见图 1。应变率效应采用经典的 Cowper-Symonds模型:yd=ys1+D()1q()(2)式中:yd为动态屈服应力;ys为静态屈服应力;为应变率;D 为材料常数;q 为材料常数。1.3 结构模型 结构模型由甩击管道和靶物管道组成(见图803第 2 期彭建,等:核电站高能管道冲击过程数值模拟研究2),甩击管道和靶物管道参数选择核电站中的典型参数。甩击管道参数见表 1;靶物管道长度 L3为4 m,外径为 88.9 mm,壁厚、约束方式、角度见表 2。甩击载荷 F 根据标准 ANSI/ANS-58.22计算获得:F=CTPA(3)式中:CT为推力系数;P 为管道正常运行工况压力;A 为管道内截面面积。图 1 真应力应变曲线Fig.1 True Stress-strain curveP 取核电站中典型压力 15.5 MPa,CT取 1.0,则F 为 92 822 N。甩击管道在锚固端采用全约束,在自由端采用壳载荷类型施加载荷 F,该载荷类型的方向伴随单元方向转动。靶物管道的约束方式见表 2。图 2 管道甩击管道模型示意Fig.2 Schematic diagram of pipe whipping pipe model表 1 甩击管道几何参数Table 1 Geometric parameters of whipping pipe外径/mm壁厚/mm弯头半径/mmL1/mmL2/mm/mm/mm114.313.491524 000300500300表 2 靶物管道参数Table 2 Parameters of target pipe工况壁厚/mm约束方式/()111.13两端固定30211.13两端固定45311.13两端固定60411.13两端固定90511.13一端固定,一端简支90611.13一端固定,一端自由90711.13两端简支9087.62两端固定9095.49两端固定90103.05两端固定901111.13两端固定701211.13两端固定802 模型验证 由于试验数据局限性,采用不同的试验分开验证管道自由甩动和冲击阶段的数值模型有效性。2.1 管道自由甩动阶段 1)实验模型文献3开展了管道自由甩动试验,试验装置见图 3。本实验主要由储气罐、试验管道和爆破片3 个部分组成。储气罐模拟设备提供一个稳定的压力源,该储气罐足够大,整个实验过程中压力下降得非常小。爆破片模拟管道的双端剪切断裂,压力达到指定的限值,在 1 ms 内完成爆破,触发试验管道自由甩动,试验参数见表 3。管道的运动使用高速胶片摄像机(Hadlands HYSPEED)进行记录,它提供了一系列管道甩动的瞬时剖面。仿真算法、材料模型和结构模型,基于表 3 的试验参数进行数值模拟。图 3 管道自由甩动装置Fig.3 Schematic diagram of pipe free whip experimental facility表 3 管道自由甩动参数Table 3 Parameters of pipe free whipping外径/mm壁厚/mm长度/mm线质量/(kg m-1)压力/MPa杨氏模量/GPa屈服强度/MPa50.802.603 0003.226.202002952)结果分析通过改变网格尺寸开展最小稳定时间步长的敏感性分析,网格长度与管道周长的比值为 1/30、1/15、1/10、1/6,分别对应最小稳定时间步长为5.010-7、1.110-6、1.710-6、2.210-6。图 4 为不同网格尺寸下管道自由端的变形时程图,从图可知,不同的时间步长,变形最大差异为 2%,因此,后续自由甩动分析的最小稳定时间步长达到 10-7的量级,变形结果的精度是足够的。图 5 为文献3中试验和 ABAQUS 数值模型的自由甩动变形结果对比图,从图可知,变形结果基本吻合。2.2 管道冲击管道阶段 1)实验模型文献8开展了管道冲击管道试验,试验装置见图 6。本实验主要由弹簧弹射器、甩动管道和靶903哈 尔 滨 工 程 大 学 学 报第 44 卷物管道 3 部分组成。弹簧弹射器提供甩动管道的能量,在弹簧作用下,甩动管道对靶物管道产生冲击作用,试验参数见表 4。管道甩动和冲击过程采用高速摄 像 机(Kodak Vidicon)记 录,记 录 速 度 为2 000 帧/s。试验后按轴向和周向方向进行切片,测量了冲击段管道中心线的直径减小和弯曲角。采用第 1 章的仿真算法、材料模型和结构模型,基于表 4的试验参数进行数值模拟。2)结果分析通过改变网格尺寸开展最小稳定时间步长的敏感性分析,网格长度与管道周长的比值为 1/33、1/22、1/13,分别对应最小稳定时间步长为 6.5 10-7、8.510-7、1.410-6,接触约束选择拉格朗日乘子法,接触摩擦采用罚函数,摩擦系数取 0.1。如图 7所示,最小稳定时间步长达 10-7的量级时,冲击截面的变形结果基本无差异。为了分析接触约束算法的影响,保持网格长度为 1/33,分析了接触约束为罚函数的变形情况,如图 7 所示,与拉格朗日乘子法的变形结果基本一致。另外,还分析了接触摩擦的影响,从图 7 可知,摩擦因素对该数值模型也基本无 影响。因此,对于该数值模型,接触约束算法和摩擦的影响可忽略,但最小稳定时间步长需达到 10-7的量级。将文献8中试验数据和 ABAQUS 计算结果进行对比,分别从靶物管道轴向、径向变形进行分析,从图 8 可知,试验结果与数值模拟基本吻合。进一步定量分析发现,冲击截面管道直径变化量 d的计算值和实验值分别为 22.88 mm 和 23.7 mm,差异仅为 3.5%(见图 9)。图 4 管道自由端的变形时程曲线Fig.4 Deformation-time curves at pipe free end图 5 试验和 ABAQUS 的甩动变形结果对比Fig.5 Comparison of whipping deformation between experiment and ABAQUS calculation图 6 管道冲击管道装置Fig.6 Schematic diagram of pipe-to-pipe impact setup表 4 管道冲击管道参数Table 4 Parameters of pipe-to-pipe impact管道外径/mm壁厚/mm长度/mm角速度/(rad s-1)杨氏模量/GPa屈服强度/MPa抗拉强度/MPa甩动管道4531 9007207340490靶物管道441840/207368504图 7 冲击截面直径变化量的敏感性分析Fig.7Sensitivity analysis of diameter reduction at im-pact section013第 2 期彭建,等:核电站高能管道冲击过程数值模拟研究图 8 试验和 ABAQUS 的靶物管道变形结果对比Fig.8Comparison of target pipe deformation between experiment and ABAQUS calculation图 9 靶物管道冲击截面直径变化图Fig.9Diameter reduction of target pipe at the impact section3 工程问题分析3.1 甩击角度效应 采用以上经验证的数值方法和模型,研究不同甩击角度对靶物管道变形的影响。以靶物管道外径88.9 mm、壁厚 11.13 mm、两端固定为例,工况 14、1112 对甩击角度从 3090的 6 种工况进行了数值模拟。从图 10 和图 11 可知,靶物管道的变形可分为 2 种模式,压缩(即局部变形)和弯曲(即整体变形),最后的变形状态是压缩和弯曲的耦合作用。压缩主要发生在冲击截面,可以通过冲击截面的直径变化量 d 表征。弯曲可以通过弯曲角 表征(见图 10)。图 12 展示了冲击截面 d 随时间的变化,冲击瞬间,冲击截面发生了较大的变形,然后随着 的增加而增加,从而进一步说明了 d 和 之间存在一定的相互作用关系。Alzheimer 等6为了综合考虑压缩和弯曲效应,定义了一个变形参数 Td来表征管道的变形行为:Td=dd()2+ddm()2|0.5(4)式中:d 为靶物管道直径;dm为甩击管道直径;为弯曲角。图 10 靶物管道弯曲应变Fig.10 Bending strain of target pipe图 11 靶物管道冲击截面压缩应变Fig.11 Crushing strain of target pipe at the impact sec-tion图 12 直径变化量时程曲线Fig.12 Diameter reduction-time curves从图 13 可知,靶物管道冲击截面的直径变化量d 和应变随着甩击角的变大而变大,这主要是因为随着甩击角度的变大,甩击管道和靶物管道的接触面积变小,从而导致冲击更加局域化。图 14 说明弯曲角 仅仅在甩击角较小的 30时较小,其他甩击角情况下的弯曲角基本一致;变形参数 Td随甩击角的变大而增大,说明其考虑了压缩效应的变化趋势。因此,甩击角度主要影响压缩变形,对弯曲变形不敏感,变形参数 Td综合表征了靶物管道的变形特性,且在甩击角为 90时变形最大。3.2 约束方式影响 为研究不同约束方式对靶物管道变形的影响,以靶物管道外径 88.9 mm、壁厚 11.13 mm、甩击角为 90,工况 47 对不同的约束方式进行了数值模拟。从表 5 中可知,约束方式对压缩基本没有影响,主要影响弯曲。工况 5 和 6 的弯曲角基本一致,而大于工况 4,小于工况 7,这是因为靶物管道轴向约113哈 尔 滨 工 程 大 学 学 报第 44 卷束不同,工况 4 两端轴向都被约束,工况 5 和 6 只有一端轴向被约束,而工况 7 两端轴向都没有约束。变形参数与弯曲角的变化趋势一致。因此,约束方式主要通过管道轴向是否约束影响管道的变形,轴向约束对靶物管道的变形影响较大,且两端简支时变形最大。图 13 直径变化量和应变随甩击角的变化Fig.13Diameter reduction and strain vs.whipping an-gle图 14 弯曲角和变形参数随甩击角的变化Fig.14Bending angle and deformation parameter vs.whipping angle表 5 靶物管道的变形结果Table 5 Deformation results of target pipe工况d/mm/()Td最大应变/%16.80120.1815.9728.95160.2420.86311.74160.2625.12416.65170.3133.77515.47390.5634.07617.86390.5735.22717.53420.6153.65848.84500.9357.85970.86581.1959.261070.87591.2059.271113.29160.2729.411215.59170.3132.103.3 径厚比的影响 为了进一步研究靶物管道的变形特性,选择了更小壁厚的靶物管道进行研究。图 15 给出了工况8 的靶物管道冲击截面的变形过程,随着时间的推移,靶物管道逐渐被压扁。塑性应变主要发生在冲击面,冲击背面的塑性应变非常小,冲击点两侧凸起,形成类似“心形”的截面形状,这可能主要是由于惯性作用导致。工况 9 和 10 的冲击截面变形与工况 8 基本一致,只是变形更加大,最后压成了饼状。图 15 工况 8 的冲击截面变形过程Fig.15 Deformation process of impact section at test 8图 16 展示了冲击点的应变变化情况:冲击瞬间,应变发生突变,随后增长趋势变缓,在 3040 ms时段,应变基本无变化,这主要因为靶物管道两端约束处形成了塑性铰,变形模式由冲击截面压缩变成了整体弯曲模式;4060 ms 时段,由于管道两端轴向约束的作用,整体弯曲刚度大于冲击截面的压缩刚度,变形模式又变成了压缩为主,从而使冲击点的应变随时间的增加而增加。应变变化的趋势与冲击截面直径变化的趋势基本一致。这说明冲击截面直径的变化可以表征冲击截面的应变特性。图 16 工况 8 的直径变化量和应变时程曲线Fig.16 Diameter reduction and strain-time curves at test 84 结论 1)靶物管道的变形包括压缩和弯曲 2 种变形213第 2 期彭建,等:核电站高能管道冲击过程数值模拟研究模式,且压缩和弯曲是相互作用的,综合考虑了压缩和弯曲的变形参数可以很好的表征靶物管道的变形行为。2)靶物管道的变形随甩击角度的增加而增加,正交甩击时变形最大;甩击角度主要影响靶物管道的压缩变形,对弯曲变形影响较小。3)靶物管道的变形与轴向是否约束相关性大,轴向无约束时变形较大;轴向约束主要影响靶物管道的弯曲变形,对压缩变形不敏感。4)随着径厚比的增加,靶物管道变形越大;冲击截面直径变化量与应变随时间的变化趋势一致。参考文献:1 Protection against internal hazards other than fires and ex-plosions in the design of nuclear power plants:IAEA No.NS-G-1.11-2004S.2 Design basis for protection of light water nuclear power plants against the effects of postulated pipe rupture:ANSI/ANS 58.2-1988S.3 REID S R,YU T X,YANG J L,et al.Dynamic elastic-plastic behaviour of whipping pipes:experiments and theo-retical modelJ.International journal of impact engineer-ing,1996,18(7/8):703-733.4 REID S R,YU T X,YANG J L.An elastic-plastic hard-ening-softening cantilever beam subjected to a force pulse at its tip:a model for pipe whipJ.Proc R Soc Lond A,1998,454:997-1029.5 REID S,YANG J L.Pipe Whip:in-plane whipping of bent cantilever pipesJ.Journal of pressure vessel tech-nology,1998,120:170-178.6 ALZHEIMER J M,BAMPTON M C C,FRILEY J R,et al.Pipe-to-pipe impact program.NUREG/CR-3231,PNL-5779R.S.l.:1987.7 BAUM M R.Pipe whip:a summary of the damage ob-served in BNL pipe-on-pipe impact testsJ.Int J Pres Ves Piping,1987,30:217-32.8 YANG J L,LU G Y,YU T X,et al.Experimental study and numerical simulation of pipe-on-pipe impactJ.Inter-national journal of impact engineering,2009,36(10/11):1259-1268.9 王春霖,佘靖策,褚金华.基于 LS-DYNA 的主蒸汽管道甩动仿真分析J.核动力工程,2011,32(S1):93-97.WANG Chunlin,SHE Jingce,CHU Jinhua.Simulation and analysis of main steam pipe whip based on LS-DYNAJ.Nuclear power engineering,2011,32(S1):93-97.10 丁凯,李岗,梁兵兵.核电站高能管道断裂防甩分析方法研究J.核动力工程,2011,32(S1):13-17.DING Kai,LI Gang,LIANG Bingbing.Study on pipe whip analysis of high energy pipe breakJ.Nuclear pow-er engineering,2011,32(S1):13-17.11 袁锋,吕勇波,艾红雷,等.核辅助管道甩击及防甩支架力学分析J.核动力工程,2013,34(6):40-42.YUAN Feng,LYU Yongbo,AI Honglei,et al.Nuclear auxiliary pipe whip and whip restraint stress analysisJ.Nuclear power engineering,2013,34(6):40-42.12 梁向东,冀英杰,何海峰,等.核电站高能管道甩击防护中的能量及力学分析J.压力容器,2019,36(1):41-47,77.LIANG Xiangdong,JI Yingjie,HE Haifeng,et al.Ener-gy and stress analysis on whip restraint of high energy pip-ing for nuclear power plantJ.Pressure vessel technolo-gy,2019,36(1):41-47,77.13 杨秀娟,闫涛,修宗祥,等.海底管道受坠物撞击时的弹塑性有限元分析J.工程力学,2011,28(6):189-194.YANG Xiujuan,YAN Tao,XIU Zongxiang,et al.Elas-tic-plastic finite element analysis of submarine pipeline impacted by dropped objectsJ.Engineering mechanics,2011,28(6):189-194.14 谢丽媛,邵永波,高旭东.单层和双层足尺度海底管道抗冲击性能分析J.振动与冲击,2021,40(1):286-296.XIE Liyuan,SHAO Yongbo,GAO Xudong.Anti-impact performance analysis of single-layer and double-layer full scale submarine pipelinesJ.Journal of vibration and shock,2021,40(1):286-296.15 张荣,支旭东,范峰,等.两端固支圆钢管侧向抗冲击性能试验研究J.天津大学学报(自然科学与工程技术版),2016,49(S1):28-33.ZHANG Rong,ZHI Xudong,FAN Feng,et al.Experi-mental study on the lateral impact of fully clamped circu-lar steel tubesJ.Journal of Tianjin University(science and technology),2016,49(S1):28-33.16 杨政龙,余建星,陈海成,等.深海管道在冲击载荷作用下的局部屈曲特性研究J.天津大学学报(自然科学与工程技术版),2019,52(3):255-261.YANG Zhenglong,YU Jianxing,CHEN Haicheng,et al.Local buckling characteristics of deep-sea pipelines under impact loadingJ.Journal of Tianjin University(science and technology),2019,52(3):255-261.17 赵保磊,余建星,孙震洲,等.深水管道在动力载荷作用下的局部压溃特性研究J.振动与冲击,2017,36(17):104-110,133.ZHAO Baolei,YU Jianxing,SUN Zhenzhou,et al.Local pressure collapse characteristics of a deep water pipeline under dynamic loading J.Journal of vibration and shock,2017,36(17):104-110,133.18 余建星,卞雪航,余杨,等.深水海底管道全尺寸压溃试验及数值模拟J.天津大学学报,2012,45(2):154-159.YU Jianxing,BIAN Xuehang,YU Yang,et al.Full-scale collapse test and numerical simulation of deepwater pipe-lineJ.Journal of Tianjin University,2012,45(2):154-159.19 ZHU Ling,LIU Qingyang,JONES N,et al.Experimental study on the deformation of fully clamped pipes under lat-eral impactJ.International journal of impact engineer-ing,2018,111:94-105.20 QU Hui,HUO Jingsi,XU Chao,et al.Numerical studies on dynamic behavior of tubular T-joint subjected to impact loadingJ.International journal of impact engineering,2014,67:12-26.本文引用格式:彭建,刘浪,石媛媛.核电站高能管道冲击过程数值模拟研究J.哈尔滨工程大学学报,2023,44(2):307-313.PENG Jian,LIU Lang,SHI Yuanyuan.Numerical simulation of the high-energy piping impact process in a nuclear power plantJ.Journal of Harbin Engineering University,2023,44(2):307-313.313