Discontinuous impact fatigue failure model and microscopic mechanism of pure titanium under high strain-rate loading
-
摘要: 基于传统的分离式霍普金森拉杆系统,设计了应变控制的冲击疲劳寿命测试实验,研究了冲击疲劳加载下纯钛的微观演化机制及冲击疲劳对材料宏观力学行为的影响。通过对不同冲击疲劳试验阶段的试样开展准静态力学性能测试,借助扫描电子显微镜 (scanning electron microscope, SEM) 和电子背散射衍射 (electron backscatter diffraction, EBSD) 技术表征试样在不同阶段的微观组织以及冲击疲劳失效后的断口形貌,研究纯钛在冲击疲劳失效过程中的循环硬化/软化规律及其微观演化机制。结果表明:通过改变子弹长度可以实现应变控制的冲击疲劳寿命测试;Manson-Coffin疲劳寿命模型可以较好地反映纯钛的冲击疲劳寿命与应变幅值之间的关系;纯钛在冲击疲劳失效过程中表现出循环硬化的现象,这主要是疲劳过程中孪生变形引起的细晶强化和塑性变形引起的应变硬化共同作用的结果,纯钛的冲击疲劳损伤主要表现为变形能力的损失。
-
关键词:
- 冲击疲劳 /
- 纯钛 /
- Manson-Coffin模型 /
- 微观组织
Abstract: The fatigue failure behavior of structural materials under repeated impact loads has always attracted much attention. Mastering its damage accumulation process and evolution mechanism at the micro-scale is the fundamental way to understand the impact fatigue failure mechanism. Due to the complexity of the impact fatigue load itself and the limitations of the current experimental equipment, there are still major problems in the study of impact fatigue failure of materials. Therefore, pure titanium was used as the research object and a strain-controlled impact fatigue life test was designed based on the traditional split Hopkinson tension bar system. The strain-controlled impact fatigue life test was achieved by changing the length of the striker, and the amplitude of the incident wave needed to be kept at the same level when using different striker tests. The relationship between strain amplitude and impact fatigue life was analyzed. The impact fatigue interruption experiments of 5 times, 10 times and 20 times were carried out with 100 mm bullets. The microstructure of the samples after different impact times were characterized by electron backscatter diffraction (EBSD) and then the quasi-static mechanical properties were tested. The fracture morphology after impact fatigue failure was observed by scanning electron microscope (SEM). The cyclic hardening/softening law and its microscopic evolution mechanism of pure titanium during impact fatigue failure were studied. The results show that the strain-controlled impact fatigue life test can be realized by changing the striker length. The Manson-Coffin fatigue life model can better reflect the relationship between impact fatigue life and strain amplitude of pure titanium. Moreover, pure titanium exhibits cyclic hardening during impact fatigue failure, which is mainly due to the combined effect of fine grain strengthening caused by twin deformation and strain hardening caused by plastic deformation during fatigue. Finally, the impact fatigue damage of pure titanium is mainly manifested as the loss of deformation ability.-
Key words:
- impact fatigue /
- pure titanium /
- Manson-Coffin model /
- microstructure
-
大型商用飞机撞击核电厂是一种潜在的威胁, 9·11事件以后尤其得到重视。对核电厂而言, 一旦发生破坏, 造成的危害不堪设想, 目前在核岛结构的设计中考虑大型商用飞机撞击的影响已经纳入规范。美国核管会在2009年将评价核设施结构防御飞机撞击能力作为联邦法规的新条款颁布[1], 同年发布NEI 07-13[2], 提供评价大型商用飞机撞击下核设施结构的响应的方法, 并在2011年正式发布了RG 1.217[3], 要求核设施申请者应提交大型商用飞机撞击下核设施结构响应的评价报告。在我国, 《核电厂厂址选择的外部人为事件》[4]中尚未强制要求评价核电厂结构防御大型商用飞机撞击能力, 但福岛核事故后, 对于一些超设计基准事件, 包括大型商用飞机撞击对核电厂的影响, 也越来越受到重视[5]。新项目审批重新开启时要求新项目须符合三代标准, 提高准入门槛, 按照国际最高安全要求新建核电项目。这说明, 新建的项目除应满足我国规范体系外, 还应按较高的国际标准进行校核验算, 这就要求新建的三代核电厂须评价核岛结构防御飞机撞击的能力。
近年来, 人们开展了在飞机撞击下混凝土结构响应的研究[6-7, 11-13], 这些研究多借鉴了其他冲击荷载(如弹体、弹丸及爆炸等)下混凝土结构响应的研究成果[6-10], 而飞机撞击混凝土结构的数值模拟及动力学分析仍需进一步深入。由于大型商用飞机撞击核电厂结构会产生两种效应, 局部破坏和整体振动, 即由飞机发动机引起的局部破坏(穿入、散裂、碎甲和穿透)和由整个飞机冲击引起的整体结构的反应分析, 因此, 需要进行两种典型的结构破坏模式的评估。由于飞机撞击混凝土结构具有短时强载的特征, 其破坏机制非常复杂[12-13]。整个撞击过程是一个瞬态、高应变率、复杂的破坏过程, 包含着几何、材料非线性方面的诸多问题, 也涉及接触面上伴随接触-撞击过程而产生的撞击能量的传递与交换、接触面的滑移、分离等非线性接触问题。在进行飞机撞击下混凝土结构响应的评价时, 通常先基于局部破坏的数值模拟, 验证选用的结构材料参数及失效准则的合理性, 再进行结构的整体响应分析。对于混凝土结构在飞机撞击下的局部破坏, 由于飞机属于硬碰撞物, 局部侵彻需要考虑[2, 10-11]。局部破坏由连续的3个阶段组成[2, 8], 飞射物穿入标靶, 标靶的散裂和碎甲, 飞射物穿透标靶。各状态定义如下:穿入(Penetration)为飞射物进入标靶; 散裂(Spalling)为标靶正前面材料的散裂(飞射物冲击的正面); 碎甲(Scabbing)为标靶后面材料的碎片(冲击面的背面); 穿透(Perforation)为飞射物完全穿透标靶。穿透速度是飞射物恰好完全穿透标靶并且没有穿出时的初始速度。剩余速度是当飞射物的初始速度大于穿透速度时, 飞射物从标靶中穿出的速度。在飞机撞击下混凝土结构响应数值模拟的精确性, 取决于混凝土材料本构方程及其参数选取的正确性, 因此数值模拟中如何选用合适的本构模型以及模型中参数, 已经成为数值模拟的关键。
本文中, 基于ANSYS/LS-DYNA的几种常用的材料参数模型, 进行在飞机撞击下混凝土结构局部破坏的数值分析, 并将数值分析结果与实验结果比较, 给出各个材料参数在模拟飞机撞击混凝土结构的适用性和准确性, 为整体评价飞机撞击下核电厂结构的防御能力提供可信的模拟方法和基础参数。
1. 基于ANSYS/LS-DYNA的数值模拟分析
1.1 实验模型
对于混凝土结构在飞机撞击下的局部破坏评价, 多采用侵彻动力分析方法。由于侵彻动力分析的准确性受多个参数的影响, 有很大的不确定性, 在进行侵彻动力分析中, 常辅以真实的实验过程及实验数据校正数值模拟结果的可靠性。比较典型的飞机撞击混凝土结构实验有, H.Tsubota等[16-18]描述的飞机撞击钢筋混凝土及钢板混凝土结构的实验。为评价飞机撞击混凝土结构的数值模拟结果, 本文中选取比较典型的飞机撞击钢筋混凝土响应的实验。飞机撞击钢筋混凝土实验中, 飞机垂直撞击钢筋混凝土靶板; 钢筋混凝土板的尺寸为1 500 mm×1 500 mm×600 mm, 钢筋间距为D3 mm@25 mm; 取1/7.5全尺寸飞机模型, 飞机模型由4种材料组成, 总重量为245.6 N。表 1为实验中飞机模型中4种材料的主要材料参数和重量; 图 1为实验中飞机撞击钢板混凝土结构的模型剖面示意图。
表 1 飞机模型的材料参数和重量Table 1. The parameters and the weight of aircraft model部件 材料 E/GPa Y/MPa εf W/N 机身外壳 玻璃 6.37 82.3 0.013 81.3 机身填充物 高密度泡沫 2.28 0.6 0.10 129.4 引擎 钢 206 797.4 0.20 28.1 赫氏支撑轴 碳纤维板 0.168 12.9 ∞ 6.8 1.2 有限元模型
有限元模型建立过程中, 钢筋混凝土的位移协调分离式模型, 即钢筋采用空间梁单元beam 161划分网格, 混凝土采用三维实体solid 164单元, 钢筋与混凝土单元通过Constrained_Lagrange_In_Solid连接, 飞机模型中机身外壳和支撑轴采用壳shell 163单元, 机身填充物和引擎用实体solid 164单元。为保证求解过程的高效及计算结果的收敛性, 混凝土结构模型采用映射法划分网格, 飞机有限元模型均采用与混凝土结构模型一致的网格, 整个有限元模型如图 2所示。对于材料单元中无失效准则的, 通过Mat_Add_Erosion定义失效, 飞机模型中材料的失效准则均与实验模型中提供的参数一致, 采用混凝土材料的失效准则。飞机撞击混凝土结构的接触算法通过Contact_Eroding_Surface_To_Surface和Contact_Surfact_To_Surfact定义。计算过程中不同混凝土材料模型时, 飞机模型与混凝土结构间的接触算法不变, 同时为防止出现沙漏现象, 还加入了Control_Energy和Control_Hourglass。
1.3 混凝土本构模型
ANSYS/LS-DYNA广泛应用于爆炸和冲击荷载作用下结构响应的计算, 提供了丰富的混凝土材料模型[14-15], 如:Mat_Pseudo_Tensor(MAT016);Mat_Concrete_Damage_R3(MAT072R3);Mat_Winfrith_Concrete(MAT084);Mat_Johson_Holmquist_Concrete(MAT111);Mat_Soil_Concrete(MAT078);Mat_Brittle_Damage(MAT096);Mat_CSCM_Concrete(MAT159)。每个模型都有自己的优缺点, 需要进行分析才能选到合适的混凝土材料模型。本文中选取4种模型进行飞机撞击混凝土结构的模拟, 用以分析各个混凝土材料参数模拟飞机撞击混凝土结构的反应规律。
1.3.1损伤混凝土材料模型
MAT072R3为损伤混凝土材料模型, 是原72号混凝土材料模型的改进版, 采用3个剪切失效面, 同时考虑损伤和应变率的影响。由于材料模型的参数众多且复杂, 最后又对材料模型添加了参数自动生成功能, 该模型一般只需要混凝土抗压强度等较少的信息, 所以该材料模型适合于模拟混凝土爆炸以及碰撞载荷下的响应。
1.3.2连续帽盖混凝土模型
MAT159为连续帽盖混凝土模型, 该模型针对用于安全设施的混凝土模型, 如高速路旁的混凝土护栏、桥梁护栏等在受到碰撞时, 其变形和破坏能够很好地用该材料模型模拟, 它是地帽模型的扩展, 包括各向同性连续方程、屈服及硬化强度面、破坏方程, 同时也考虑了应变率效应。该材料模型主要用于模拟埋有钢筋的加强混凝土, 可用钢筋、混凝土耦合算法共同模拟钢筋混凝土。
1.3.3混凝土累积损伤模型
MAT111是混凝土累积损伤模型, Johnson-Holmquist-Concrete本构模型及损伤模型是针对混凝土在大应变、高应变率及高围压条件下提出的一种混凝土材料本构关系, 该模型考虑了材料损伤、应变率效应以及静水压力对于屈服应力的影响。它能较好地描述混凝土在高速撞击与侵彻下所产生的损伤、破碎及断裂(或层裂), 在数值模拟中得到了广泛的使用。
1.3.4混凝土Ottosen塑性模型
MAT084是基于剪切破坏提出的混凝土的塑性模型, 主要用于解决冲击荷载作用下钢筋混凝土的响应, 该模型允许每个单元的三向裂纹开展, 而且可在Ls-Prepost查看裂纹。该模型主要应用在模拟冲击荷载作用下混凝土结构的响应。
这些材料模型的理论公式及详细介绍可参见ANSYS/LSDYNA使用手册, 参见文献[14-15]。在第三代核电机组设计中评价飞机撞击混凝土安全壳结构时, 这些模型有很好的应用, 如AP1000核电机组设计中评价屏蔽厂房在大型商用飞机撞击下的响应时, 采用的是混凝土Ottosen塑性模型; EPR核电机组设计中评价混凝土安全壳结构在大型商用分析撞击下的响应时, 采用连续帽盖混凝土模型。
2. 数值分析结果
2.1 计算结果
飞机撞击混凝土钢板的数值计算过程为15 ms。5 ms左右前, 飞机引擎未到达混凝土靶板, 混凝土几乎没有破坏; 5 ms左右时, 飞机引擎开始撞击混凝土靶板; 9 ms左右, 整个飞机模型即将穿透混凝土靶板; 11 ms后, 飞机引擎穿透混凝土结构后飞离混凝土靶板。图 3为9 ms时靶板的正面破坏状态; 图 4为9 ms时靶板的背面的破坏状态; 图 5为9 ms时飞机穿透靶板的侧视图。由图可知, MAT084、MAT159、MAT111和MAT72R3均能模拟飞机撞击混凝土结构的局部效应, 飞机引擎撞击混凝土结构时均出现混凝土的破坏失效和钢筋的断裂的现象; 混凝土结构的背面破坏面积大于正面破坏面积; 相同撞击条件下, 根据混凝土结构破坏面积得到MAT072R3和MAT084两种模型计算结果较接近, MAT111和MAT159两种模型计算结果较接近, 且MAT072R3和MAT084模型比MAT111和MAT159模型破坏面积大。
2.2 与实验结果的比较
为更精确地说明4种混凝土材料模型在模拟飞机撞击混凝土结构的使用性, 图 6为飞机引擎和穿透机身速度的变化过程; 表 2列出了飞机撞击混凝土的数值分析结果与实验结果的对比。表中, dp为正面破坏面积直径, Ab为背面破坏面积。从图 6中可以看到, MAT084、MAT111和MAT72R3变化趋势一致, 穿透机身外壳的剩余速度在110 m/s左右, 引擎的剩余速度在60 m/s左右。表 2列出的实验结果机身外壳的剩余速度在110 m/s, 引擎的剩余速度在80 m/s。MAT084、MAT111和MAT72R3材料模型计算的机身外壳的剩余速度均与实验结果较接近; 对于引擎剩余速度, MAT72R3计算的68m/s较接近于实验结果。
表 2 数值分析结果与实验结果的对比Table 2. Comparison of numerical and the experimental results模型 v0/(m·s-1) vr/(m·s-1) dp/cm Ab/(cm·cm) 机身 引擎 实验 142 110 80 33 50×55 MAT084 142 108 52 30 39×42 MAT159 142 10 5 24 18×27 MAT111 142 113 58 18 30×33 MAT072R3 142 113 68 30 42×42 3. 结论
利用ANSYS/LSDYNA提供的可用于模拟冲击作用下混凝土结构性能的4种材料参数, 对日本Kobori综合研究公司的飞机撞击混凝土结构的实验进行了数值模拟分析, 比较了不同混凝土材料模型在飞机撞击下混凝土结构动力分析时, 混凝土结构靶板的正、背面破坏面积及飞机剩余速度等计算结果。得到:
(1) ANSYS/LSDYNA中混凝土材料模型较多, 各自有着其适用的领域, 在相同分析计算领域中, 不同的混凝土材料模型计算结果也不尽一致。选择合适的混凝土材料模型本构, 有利于提高数值分析计算结果的准确度。
(2) 飞机撞击混凝土结构数值模拟中, 采用的4种混凝土模型参数均能模拟飞机引擎穿入、散裂、碎甲等混凝土结构的局部破坏效应, 但计算的飞机撞击混凝土的正、背面破坏面积及剩余速度不同, MAT072R3和MAT084材料模型的计算结果与实验结果较接近, MAT111材料模型次之, 而MAT159材料模型有较大的差异, 正、背面破坏面积偏小, 剩余速度几乎降低为零。
(3) 由于数值分析计算结果存在一定的不确定性, 所以采用的MAT072R3、MAT084和MAT111混凝土材料模型参数均可为后续评估混凝土结构安全壳抵抗飞机撞击能力时提供基础参数, 也可应用于飞机撞击混凝土结构最不利撞击角度、飞机撞击混凝土结构的穿透公式及飞机撞击混凝土结构的安全壳的整体响应分析等方面的进一步研究中。
-
表 1 不同冲击应变幅值下纯钛的冲击疲劳寿命
Table 1. Impact fatigue life of pure titanium under different impact strain amplitudes
Δεt/2 冲击疲劳寿命/次 0.017 58 62 0.023 21 23 25 0.026 12 16 17 0.035 10 8 11 0.051 5 6 7 0.070 4 4 4 -
[1] 刘正, 胡冶昌, 魏志芳. 复进簧冲击疲劳应力响应及其寿命预测 [C] //首届兵器工程大会论文集. 2017: 52–57. [2] 希弦. 微观航母之舰载机拦阻钩 [J]. 兵器知识, 2015(3): 72–75. DOI: 10.19437/j.cnki.11-1470/tj.2015.03.016.XI X. Arresting hook of carrier-based aircraft on aircraft carrier [J]. Ordnance Knowledge, 2015(3): 72–75. DOI: 10.19437/j.cnki.11-1470/tj.2015.03.016. [3] JOHNSON A A, STOREY R J. The impact fatigue properties of iron and steel [J]. Journal of Sound and Vibration, 2007, 308(3/4/5): 458–466. DOI: 10.1016/j.jsv.2007.06.044. [4] YANG S S, BAI C Y, YANG Q, et al. Review on impart fatigue of metallic materials and structures [J]. Aeronautical Science & Technology, 2021, 32(2): 1–13. DOI: 10.19452/j.issn1007-5453.2021.02.001. [5] JGUCHI K T H, TAIRA S. Failure mechanisms in impact fatigue of metals [J]. Fatigue and Fracture of Engineering Materials and Srrucfures, 1979, 2(2): 165–176. DOI: 10.1111/j.1460-2695.1979.tb01352.x. [6] NAKAYAMA H, TANAKA T. Impact fatigue crack growth behaviors of high strength low alloy steel [J]. International Journal of Fracture, 1984, 26(9): 19–24. DOI: 10.1007/BF01152319. [7] TANAKA T, KINOSHITA K, NAKAYAMA H. Fatigue crack growth and microscopic crack opening behaviour under impact fatigue load [J]. International Journal of Fatigue, 1989, 11(2): 117–123. DOI: 10.1016/0142-1123(89)90006-6. [8] YANG P, LIAO X, ZHU J, et al. High strain-rate low-cycle impact fatigue of a medium-carbon alloy steel [J]. International Journal of Fatigue, 1994, 16(5): 327–330. DOI: 10.1016/0142-1123(94)90270-4. [9] ZHANG M, YANG P S, TAN Y X, et al. An observation of crack initiation and early crack growth under impact fatigue loading [J]. Materials Science and Engineering: A, 1999, 271(1/2): 390–394. DOI: 10.1016/S0921-5093(99)00264-6. [10] 李会会, 易丹青, 刘会群, 等. 硬质合金冲击疲劳行为的研究 [J]. 硬质合金, 2014, 31(2): 100–111. DOI: 10.3969/j.issn.1003-7292.2014.02.006.LI H H, YI D Q, LIU H Q, et al. Research on impact fatigue behaviour of cemented carbide [J]. Cemented Carbide, 2014, 31(2): 100–111. DOI: 10.3969/j.issn.1003-7292.2014.02.006. [11] 陈鼎, 姚亮, 陈振华, 等. WC-Co类硬质合金的低周冲击疲劳性能研究 [J]. 稀有金属与硬质合金, 2017, 45(3): 71–76. DOI: CNKI:SUN:XYJY.0.2017-03-014.CHEN D, YAO L, CHEN Z H, et al. Study on low cycle impact fatigue performance of WC-Co cemented carbides [J]. Rare Metals and Cemented Carbides, 2017, 45(3): 71–76. DOI: CNKI:SUN:XYJY.0.2017-03-014. [12] STANTON L B. The resistance of materials to impact [J]. Proceedings of the Institute of Mechanical Engineers, 1908: 889–919. DOI: 10.1243/PIME_PROC_1908_075_019_02. [13] JOHNSON D J. The impact fatigue properties of pearlitic plain carbon steels [J]. Fatigue of and Fracture Engineering Materials and Structures, 1981, 4(3): 279–285. DOI: 10.1111/j.1460-2695.1981.tb01125.x. [14] JOHNSON A A. The low cycle impact fatigue properties of pearlitic plain carbon steels [J]. Fatigue and Fracture of Engineering Mateirals and Structures, 1985, 8(3): 287–294. DOI: 10.1111/j.1460-2695.1985.tb00428.x. [15] 张遥辉. 钢铁材料冲击疲劳行为综述 [J]. 中国设备工程, 2020(6): 211–216. DOI: 10.3969/j.issn.1671-0711.2020.06.131.ZHANG Y H. Overview of impact fatigue behavior of steel materials [J]. China Plant Engineering, 2020(6): 211–216. DOI: 10.3969/j.issn.1671-0711.2020.06.131. [16] SUN Q, LIU X R, LIANG K. Impact fatigue life prediction for notched specimen of steel AerMet100 subjected to high strain rate loading [J]. International Journal of Applied Mechanics, 2018, 10(3): 1850030. DOI: 10.1142/s1758825118500308. [17] WANG B W, QIAN C C, BAI C Y, et al. Study on impact fatigue test and life prediction method of TC18 titanium alloy [J]. International Journal of Fatigue, 2023, 168: 1–17. DOI: 10.1016/j.ijfatigue.2022.107391. [18] GAO P F, LEI Z N, LI Y K, et al. Low-cycle fatigue behavior and property of TA15 titanium alloy with tri-modal microstructure [J]. Materials Science and Engineering: A, 2018, 736: 1–11. DOI: 10.1016/j.msea.2018.08.080. [19] 张欠欠, 刘晓燕, 雷罗, 等. 工业纯钛的室温低周疲劳行为 [J]. 塑性工程学报, 2019, 26(2): 219–224. DOI: 10.3969/j.issn.1007-2012.2019.02.029.ZHANG Q Q, LIU X Y, LUO L, et al. Low-cycle fatigue behavior of commercially pure titanium at room temperature [J]. Journal of Plasticity Engineering, 2019, 26(2): 219–224. DOI: 10.3969/j.issn.1007-2012.2019.02.029. [20] CHANG L, LV C, KITAMURA T, et al. Slip dominated planar anisotropy of low cycle fatigue behavior of commercially pure titanium [J]. Materials Science and Engineering: A, 2022, 854(9):143807. DOI: 10.1016/j.msea.2022.143807. [21] SONG X P, CHEN G L, GU H C. Low cycle fatigue behavior of commercial purity titanium in liquid nitrogen [J]. International Journal of Fatigue, 2002, 24: 49–56. DOI: 10.1016/S0142-1123(01)00047-0. [22] ABDUL-LATIF A. Continuum damage model for low-cycle fatigue of metals: an overview [J]. International Journal of Damage Mechanics, 2021, 30(7): 1036–1078. DOI: 10.1177/1056789521991620. 期刊类型引用(2)
1. 王亚磊, 郑立刚, 于水军, 朱小超, 李刚, 杜德朋, 窦增果. 约束端面对管内甲烷爆炸特性的影响. 爆炸与冲击. 2019(09): 139-148 . 本站查看
2. 郑立刚, 王亚磊, 于水军, 朱小超, 李刚, 杜德朋, 窦增果. NaHCO_3抑制瓦斯爆炸火焰与压力的耦合分析. 化工学报. 2018(09): 4129-4136 . 百度学术
其他类型引用(0)
-