引言

TC4钛合金因密度低、强度高、韧性好蹬着点，宽泛利用于航空航天、国防军工等关键领域[1]，但该资料对应力集中敏感，且切削加工难度大[2]。金属增材造作技术的发展为钛合金复杂构件的造作提供了便捷蹊径，进一步拓展了其在航空航天、工业设备等领域的利用领域[3]。然而，增材造作过程中逐层堆积的工艺个性易导致构件内部产生孔隙、未熔合等冶金缺点，显著降低钛合金的委顿机能[4]。增材造作钛合金构件的委顿机能普遍低于传统铸造件，而委顿失效占机械结构故障总数的90%以上。委顿断裂通常无显著宏观塑性变形预兆，易引发苦难性变乱，造成沉大经济损失与社会风险[5]。其中，委顿裂纹扩大行为的评价与预测是预防和节造委顿失效的关键环节[6]。裂纹尖端的塑性变形是决定委顿裂纹扩大速度的主题成分。因而，发展裂尖左近应力应变场的分析拥有沉要的工程意思。

成形方向是影响增材造作构件力学机能各向异性的关键成分，其对裂尖委顿机能的影响尤为显著。现有钻研多聚焦于成形方向对增材造作钛合金委顿裂纹扩大速度及宏观委顿机能的影响。RANS等[7]钻研了5种倾斜成形方向对增材造作TC4薄板委顿机能的影响，发现柱状晶粒取向对委顿裂纹扩大机能无显著影响。白澄岩等[8]系统分析了成形方向、冶金缺点、显微组织及表表处置对增材造作Ti-6Al-4V低周委顿机能的影响，阐了然低周委顿危险机造。了局批注，水平成形试样的委顿机能优于垂直成形试样。CAIN等[9]对比了3种成形方向及热处置工艺对增材造作钛合金力学机能的影响，发现热处置前成形方向对断裂韧性及委顿裂纹扩大速度的影响最为显著。XU等[10]通过原位扫描电镜观察了分歧成形方向试样的委顿裂纹萌生与扩大过程，发现90°成形试样的抗裂纹扩大能力显著降低。POULIN等[11]探索了委顿裂纹扩大行为与成形方向的关联，了局批注，裂纹扩大门槛值及近门槛区扩大速度均受成形方向的显著影响。

裂尖委顿机能由部门循环变形行为决定，发展成形方向对增材造作资料裂尖循环塑性行为的钻研，有助于深刻揭示委顿危险演化法规。然而，目前关于成形方向对委顿裂纹尖端循环塑性区应力、应变影响的钻研较少。针对上述问题，选取有限元步骤分析裂尖左近应力、应变场的演化法规，钻研激光选区溶解增材造作TC4钛合金在L1(0°)、L2(45°)、L3(90°)3种成形方向下的裂尖循环塑性行为，并探求载荷迸纂载荷幅对其的影响法规。

1、资料与试样

钻研选取的资料为TC4钛合金，是一种中等强度的α-β型两相钛合金，含有6%的α不变元素Al和4%的β不变元素V，具体化学成分如表1所示。原资料为球形粉末，资料粒径散布为15~53 μm。激光选区溶解成形工艺参数：激光功率为340 W；扫描间距为0.12 mm；层厚为60 μm；扫描速度为1250 mm/s。为解除成形过程中的残存应力，打印实现后，选取800 ℃的温度保温2 h，以随炉冷却的方式进行热处置。

表1 增材造作TC4钛合金化学成分

金相组织如图1所示，组织结构重要为少量的β晶（图中玄色）、部门细长针状α晶（图中白色针状）以及部门α晶（图中白色板条状）组成，与张海英等[12]的钻研了局一致。图2为L1、L2、L3这3种成形方向示意图。除此之表，资料内部存在少量未熔合及孔隙缺点。设计并造备与打印方向呈L1(0°)、L2(45°)、L3(90°)的3种成形方向的紧凑拉伸（Compact Tension, CT）试样，如图2所示。具体尺寸参数为：试样厚度B=5 mm，试样宽度W=25 mm；裂纹长度a=13 mm，如图3所示。

2、有限元分析

金属资料的委顿裂纹扩大速度由裂纹前沿的应力、应变状态决定，裂纹尖端的塑性变形对裂纹扩大拥有决定性影响。忽略裂尖循环塑性变形将无法正确诠释委顿裂纹扩大的诸多景象[13]。在循环载荷作用下，裂尖左近会产生塑性应变堆集[14-15]，即便整体应力场处于弹性领域，裂尖部门仍会产生塑性变形[16]。对于无数韧性金属资料，裂尖常超出幼领域屈服领域，线弹性断裂力学不再合用。因而，成形方向对增材造作钛合金裂尖部门循环塑性行为的影响，直接决定了分歧方向的委顿裂纹扩大速度。成形方向引起的资料宏观力学机能各向异性，会进一步扭转裂尖左近的部门力学行为。

选取有限元分析步骤可能从宏观方面分析裂尖左近的力学行为和裂尖左近部门委顿危险演变法规的整体情况。通过对3种成形方向的资料作力学机能测试获取资料参数。使用Abaqus仿真仿照软件成立带裂纹的CT试样模型，将3种分歧成形方向得到的资料参数输入到资料属性中，从而实现3种分歧成形方向的CT试样裂尖力学行为分析。其中，有限元模型使用Chaboche资料本构模型进行分析，该模型结合各向同性强化与随动强化机造，可正确描述资料在循环加载过程中的包辛格效应、棘轮效应及硬化、软化行为[17]。

表2 增材造作TC4钛合金Chaboche本构模型参数

混合强化机造使其可能精确仿照资料在反复加载卸载过程中的应力-应变响应，合用于循环塑性行为分析；诘4强度理论，屈服面表白式为

式中，σ为应力张量；α为非线性随动强化变量；R_k为等向强化变量；k为初始的屈服面半径；J为偏应力空间的Mises距离，界说为

式中，α^{dev}为背应力的偏应力张量；S为偏应力张量。背应力α能够暗示成很多背应力分量的组合，其中，第i个背应力分量α_i的演化模型界说为

式中，σ?为屈服面的尺寸；C_i、γ_i均为随动强化参数； ε ˉ pl为等效塑性应变率；α_i为背应力的第i个分量。背应力α的表白式为

式中，N为背应力分量的个数。通过单轴拉伸试验拟合得到TC4钛合金的Chaboche模型参数，如表2所示。

凭据试验数据，成立等效塑性应变与屈服面尺寸的对应关系，实现各向同性强化的描述。资料的等向强化曲线如图4所示。由图4可知，随着循环塑性应变的累积，资料出现显著的循环硬化个性。

图5所示为划分的有限元网格，选取陆续的四节点平面应变单元（CPE4）进行离散，对裂纹尖端区域进行网格细化以正确捉拿应力-应变梯度。选取加强的沙漏节造，以预防沙漏模式[18]。针对每种成形方向，设置分歧的载荷比R和载荷幅ΔF，共设计9种分析规划，如表3所示。

表3 增材造作CT试样有限元分析工况规划

3、了局与会商

图6和图7所示为距离裂纹尖端0.94 μm处加载方向塑性应变随循环周次的演化曲线。由图6、图7可知，循环加载过程中裂尖塑性应变持续演化，并因塑性应变累积产生棘轮变形。在每个加载周期内，塑性应变在加载峰值处达到最大值，卸载过程中塑性应变产生部门回复。由图6可知，在一样载荷迸纂载荷幅前提下，L?方向累积的塑性应变总是最大，L?方向次之，L?方向的塑性应变最幼。且成形方向对塑性应变的影响随循环周次的增长逐步显著。由图7(a)可知，高载荷比会加快塑性应变累积，初次加载时高载荷比下的塑性应变显著高于低载荷比工况。由图7(b)可知，载荷幅增大导致应力幅升高，进而加快塑性应变累积速度；低载荷幅下塑性应变累积速度显著低于高载荷幅工况。

图8所示为第15次循环峰值时刻在分歧成形方向的"加载方向应力-与裂尖的距离"曲线。由图8可知，裂尖区域存在显著应力集中，应力值随与裂尖距离的增长呈衰减趋向，应力集中是导致构件委顿失效的重要诱因。L?方向的应力集中水平最高，应力峰值最大；L?方向次之；L?方向下的应力峰值最幼。但随着与裂尖距离的增长，成形方向对应力散布的影响逐步减弱。图9所示为第15次循环峰值时刻在分歧载荷比或载荷幅下的"加载方向应力-与裂尖的距离"曲线。由图9可知，载荷迸纂载荷幅对应力散布的影响显著：低载荷比、低载荷幅工况下应力峰值较幼，应力衰减梯度较大，高应力区域领域较；高载荷比、高载荷幅工况则相反，批注高载荷迸纂高载荷幅会加剧裂尖应力集中。

图10所示为第15次循环峰值时刻在分歧成形方向的"加载方向塑性应变-与裂尖的距离"曲线。由图10可知，裂尖左近因应力集中产生较大塑性应变，塑性应变值随与裂尖距离的增长逐步减幼，在特定地位出现拐点后，最终趋于0。在裂尖左近区域，L?方向的塑性应变最大，L?方向的塑性应变最幼。但随着与裂尖距离的增长，3种成形方向的塑性应变差值逐步缩幼，最终趋于一致。图11所示为第15次循环峰值时刻在分歧载荷比或载荷幅下的"加载方向塑性应变-与裂尖的距离"曲线。由图11可知，载荷迸纂载荷幅对塑性应变散布的影响显著，各循环周期内高载荷比、高载荷幅工况下的塑性应变均高于低载荷比、低载荷幅工况。

图12所示为第15次循环峰值时刻在分歧成形方向的"部门应力比-与裂尖的距离"曲线。由图12可知，裂尖部门应力迸纂宏观载荷比存在差距，随与裂尖距离的增长，部门应力比总体出现由负转正、先增大后减幼的趋向，最终趋于宏观载荷比并维持不变，批注在循环加载谷值点，裂尖周围存在肯定领域的压应力区。成形方向对部门应力比散布的影响较幼，3种成形方向的曲线险些沉合。图13所示为第15次循环峰值时刻在分歧载荷比或载荷幅下的"部门应力比-与裂尖的距离"曲线。由图13可知，高载荷比下裂尖均匀应力水平较高；低载荷比下裂尖均匀应力水平较低，但裂尖左近压应力区领域更大。载荷幅对均匀应力比的影响较幼，但高载荷幅下拐点地位更靠后，批注其应力集中区域更大。

图14所示为第15次循环峰值时刻在分歧成形方向的"Mises应力-与裂尖的距离"曲线，反映了裂尖的综合应力状态。由图14可知，Mises应力总体随与裂尖距离的增长而减幼，裂尖地位因应力集中出现颠簸。3种成形方向的Mises应力散布根基一致。图15所示为第15次循环峰值时刻在分歧载荷比或载荷幅下的"Mises应力-与裂尖的距离"曲线，反映了裂尖的综合应力状态。由图15可知，高载荷比、高载荷幅工况下Mises应力值更高。图16和图17所示别离为等效塑性应变在分歧成形方向及分歧载荷比、载荷幅下随与裂尖的距离的变动，其法规与Mises应力类似。

图18~图21所示为CT试样整体弹性应变能与塑性应变能的分析了局。对比弹性应变能与塑性应变能的数值可知，只管塑性应变能随循环周次逐步增长，但弹性应变能始终占主导职位。由图18可知，成形方向对弹性应变能的影响可忽略不计。由图19可知，载荷迸纂载荷幅对弹性应变能影响显著，高载荷比、高载荷幅工况下的均匀弹性应变能更高。由图20可知，塑性应变能随循环周次的增长逐步增大，成形方向对塑性应变能的影响显著：L?方向的塑性应变能最高，L?方向次之，L?方向最低，且该差距随循环周次的增长逐步扩大。由图21可知，高载荷比、高载荷幅工况下的塑性应变能显著高于低载荷比、低载荷幅工况，且低载荷工况下塑性应变能的增长速度更慢。

4、试验验证

为验证有限元分析了局的靠得住性，选取数字图像有关（Digital Image Correlation, DIC）技术实时丈量委顿加载过程中的部门变形，分析成形方向对应变场的影响。该技术通过度析试件表表散斑的灰度变动获取位移场与应变场，拥有非接触、实时丈量蹬着点，宽泛利用于委顿试验钻研[19-20]。试验装置如图22所示。拔取裂纹耽搁线作为分析蹊径L′，由图22(c)可知，蹊径L′上裂尖根部区域存在显著的应力集中。试验加载前提为：载荷比R=0.1，载荷幅ΔF=2076.3 N，试样蕴含3种成形方向。

由于DIC技术可直接获取总应变，塑性应变与应力可通过Ramberg-Osgood（R-O）本构方程换算得到，其推算式为

式中，ε为总应变；ε_e、ε_p别离为弹、塑性应变分量；σ为应力；E为弹性模量；H为强度系数；n为应变硬化指数。涉及到的资料力学机能参数如表4所示。

I型裂纹在蹊径L′上的应力场解析表白式为

式中，σ?、σ?、σ?别离为第1、第2、第3主应力；σ_eq为Mises等效应力；K?为I型裂纹应力强度因子；ν为泊松比，取0.34；r为与裂尖的水平距离。

表4 分歧成形方向增材造作TC4钛合金根基力学机能

在DIC分析中，应力取加载峰值时刻的Mises等效应力，应变取峰值处的等效塑性应变。图23为3种成形方向试样加载峰值时刻的应变云图。由图23可知，L?方向的裂尖变形最高，L?方向最低。图24所示为第15次加载峰值时刻塑性应变与Mises应力沿蹊径L′的散布曲线。由图24可知，应变与应力均随与裂尖距离的增长而减；3种成形方向的应变差值随与裂尖距离的增长而逐步减幼，应力散布差距则不显著。

图25、图26所示为DIC试验了局与有限元分析了局的对比情况。由图25可知，3种成形方向的Mises应力DIC丈量值均高于有限元推算值，且介于平面应力与平面应变解析解之间，与平面应力解析解更为靠近。在距离裂尖约莫0.02 mm以内，有限元了局显著低于解析解，在距离裂尖0.02 mm以表，有限元了局与解析解根基一致。产生该差距的重要原因是，DIC丈量的是试样表表的变形，更靠近平面应力状态；而有限元分析选取平面应变模型。裂尖左近区域因塑性变形的存在，不再满足线弹性断裂力学的合用前提；同时塑性变形导致裂尖钝化，引起应力松弛，因而距裂尖约0.02 mm以内的区域不受应力强度因子K?节造。由图26可知，3种成形方向塑性应变的DIC丈量值与有限元推算值随距裂尖距离的变动趋向一致，但裂尖左近DIC丈量值普遍高于有限元推算值。该差距重要源于平面应力状态下的裂尖塑性区尺寸显著大于平面应变状态。试验了局与有限元分析了局整体吻合，验证了有限元模型的靠得住性。

5、结论

选取有限元步骤分析了成形方向对增材造作钛合金裂尖循环塑性行为的影响，得到如下重要结论：

在整个循环过程中，L?方向裂尖塑性应变始终最大，L?方向的塑性应变均最幼，且该差距随循环周次的增长逐步扩大；成形方向对分歧地位处的应力及部门应力迸装响较幼。

成形方向对弹性应变能的影响可忽略，但对塑性应变能影响显著，L?成形方向塑性应变能高，L?方向次之，L?方向最幼。一样前提下，高载荷比、高载荷幅工况下的均匀弹性应变能与塑性应变能均更高。

DIC试验了局与有限元分析了局吻合优良，验证了成形方向对裂尖应力、应变场影响法规的正确性。

参考文件

[1] 詹奇云，靳刚，韩进，等. TC4钛合金低温拉伸行为与本构建模[J]. 塑性工程学报，2024，31(1)：204-215.

[2] 张天任，蒋睿嵩，王栋辉，等. TC4钛合金加工表表的超声冲击强化表表齐全性尝试[J]. 机械科学与技术，2025，44(11)：1955-1965.

[3] 高韬. 增材造作钛合金TC4的变形及失效机理钻研[D]. 十堰：十堰大学，2023：14.

[4] 王华明. 高机能大型金属构件激光增材造作：若干资料基础问题[J]. 航空学报，2014，35(10)：2690-2698.

[5] 李军. 弹塑性状态下复合型委顿裂纹扩大钻研[D]. 北京：北京理工大学，2014：9.

[6] PAUL S K, TARAFDER S. Cyclic plastic deformation response at fatigue crack tips[J]. International Journal of Pressure Vessels and Piping, 2013, 101: 81-90.

[7] RANS C, MICHIELSSEN J, WALKER M, et al. Beyond the orthogonal: on the influence of build orientation on fatigue crack growth in SLM Ti-6Al-4V[J]. International Journal of Fatigue, 2018, 116: 344-354.

[8] 白澄岩，兰亮，辛如意，等. 增材造作Ti-6Al-4V钛合金低周委顿机能钻研进展[J]. 资料科学与工艺，2023，31(1)：79-90.

[9] CAIN V, THIJS L, VAN HUMBEECK J, et al. Crack propagation and fracture toughness of Ti6Al4V alloy produced by selective laser melting[J]. Additive Manufacturing, 2015, 5: 68-76.

[10] XU Z W, LIU A, WANG X S. The influence of building direction on the fatigue crack propagation behavior of Ti6Al4V alloy produced by selective laser melting[J]. Materials Science and Engineering: A, 2019, 767: 138409.

[11] POULIN J R, BRAILOVSKI V, TERRIAULT P. Long fatigue crack propagation behavior of Inconel 625 processed by laser powder bed fusion: influence of build orientation and post-processing conditions[J]. International Journal of Fatigue, 2018, 116: 634-647.

[12] 张海英，李刚，臧伟锋，等. 激光选区溶解Ti-6Al-4V钛合金委顿裂纹扩大机能钻研[J]. 机械强度，2023，45(6)：1355-1360.

[13] VORMWALD M. Effect of cyclic plastic strain on fatigue crack growth[J]. International Journal of Fatigue, 2016, 82: 80-88.

[14] MCCLUNG R C. Crack closure and plastic zone sizes in fatigue[J]. Fatigue & Fracture of Engineering Materials & Structures, 1991, 14(4): 455-468.

[15] IINO Y. Local fatigue damage accumulation around Notch attending crack initiation[J]. Metallurgical and Materials Transactions A, 1995, 26(6): 1419-1430.

[16] JAMES M N, CHRISTOPHER C J, LU Y W, et al. Local crack plasticity and its influences on the global elastic stress field[J]. International Journal of Fatigue, 2013, 46: 4-15.

[17] 宋建，阎琨. 基于Chaboche模型的3D打印GH3536资料本构参数拟合钻研[J]. 塑性工程学报，2023，30(10)：167-173.

[18] 齐威. ABAQUS 6.14超等进建手册[M]. 北京：人民邮电出版社，2016：109.

[19] 代巧，何爵亨，周金宇，等. 基于DIC与Irwin模型的工业纯钛委顿裂纹尖端循环应变场表征[J]. 罕见金属资料与工程，2021，50(8)：2815-2822.

[20] GOMEZ GONZALES G L, ORTIZ GONZALEZ J A, ANTUNES F V, et al. Experimental determination of the reversed plastic zone size around fatigue crack using digital image correlation[J]. Theoretical and Applied Fracture Mechanics, 2023, 125: 103901.

（注，原文标题：成形方向对增材造作钛合金裂尖循环塑性行为的影响_刘庆）