硬质合金与钢电子束焊接接头缺陷及断裂分析

陈国庆,柳峻鹏,张秉刚,冯吉才

doi:10.12073/j.hjxb.20160126001

摘要: 对WC-Co/40Cr进行电子束焊接试验，研究了焊缝和界面组织，同时分析了焊缝中裂纹的产生机理. 结果表明，焊缝主要由马氏体和脆性η相组成，在组织和拉应力作用下生成两种典型裂纹缺陷，一类为冷却过程中尚未凝固的液态薄膜不能填充间隙产生的结晶裂纹；一类为较大残余拉应力使得焊缝中硬脆相开裂形成的淬硬脆化裂纹. 界面由于元素扩散作用产生贫碳环境，在合适的温度梯度下会有η相包络生成. 测试接头的力学性能，发现接头平均抗剪强度为506 MPa，断裂贯穿于硬质合金热影响区、界面、焊缝处，分别呈现沿晶断裂、准解理断裂和解理断裂的特征.

关键词:

0. 序言

焊接结构疲劳评估的主S-N曲线法自提出以来，取得了很好的应用效果，该方法采用与外载荷平衡的结构应力作为参数，解决了焊接结构高周疲劳寿命预测时S-N曲线的选取困难和有限元计算时的网格敏感性问题^[1-4]. 2007年针对低周疲劳问题，在ASME BPVC VIII-2标准中给出了伪弹性算法，即将疲劳试验获得的加载—位移曲线外推插值，深入的研究表明该算法是有局限性的，为解决该局限性提出了结构应变法，并推导了弯曲为主和膜应力为主的理想弹塑性材料的结构应变理论解^[5]. 产品设计时要求应力低于屈服强度，故而很少存在低周疲劳问题，导致结构应变方法的应用较少^[6-9].

2019年，高一迪等人^[10]在国内率先在开展了等效结构应变法在低周疲劳范围的适用性研究，结果表明作为结构应力法的延伸，结构应变法可以应用到低周疲劳范围. 近年来，有些铁路货车产品出现了低周疲劳故障问题，线路测试时发现，原来认为弹性状态的构件，在运用过程中确实存在高于屈服的应力响应. 在设计前端，随着轻量化设计要求的提出，一方面结构要减重，另一方面更要充分发挥构件的潜力，所以设计也开始考虑低周疲劳设计^[11]. 因此详细梳理结构应变法理论，建立分析流程，开展工程领域的适用性研究是目前亟待开展的研究工作.

鉴于此，在前期成果的基础上，推导理想弹塑性模型下结构应变计算公式，并基于理论算法开展程序的设计和验证. 基于平面应变的有限元模型和结构应变理论解进行理论的适用性分析，提出结构应变法的适用范围，实现低周疲劳的结构应变法在铁路货车上的应用，为该方法在工程领域的推广应用奠定基础.

1. 结构应变法理论

结构应变法的理论假设为：对于完全弹塑性各向同性材料，通过弹性和塑性变形后，假定沿着厚度方向上的变形仍然维持在1个平面上；其次，结构应变法理论是以平面应变状态进行考虑，焊接接头模型如图1所示，采用六面体单元进行网格划分，平板承受拉伸载荷，则焊接接头中部为平面应变状态，两侧为平面应力状态，故可将三维问题简化为平面问题，如按最危险部位考核的话，实际疲劳评估应考虑为平面应变状态.

图 1 具有平面应变效应的焊接接头

Figure 1. Welded joints with plane strain effect

下载: 全尺寸图片幻灯片

1.1 结构应变的计算方法

将焊接接头的承载情况分为两种形式：(a)弯曲载荷较大情况，弯曲应力导致板上、下面都发生屈服；(b)膜力载荷较大情况，膜应力导致上表面发生屈服.

1.1.1 上下表面都发生屈服

对于贯穿厚度的1个截面，首先使用弹性状态的结构应力计算的方式，得到膜应力$ {\sigma _{\text{m}}} $和弯曲应力$ {\sigma _{\text{b}}} $，两者求和后显然已经超过了屈服强度$ {S_{\text{Y}}} $，如图2所示.

图 2 上下表面均发生屈服

Figure 2. Yield occurs on both the upper and lower surfaces

下载: 全尺寸图片幻灯片

由于塑性变形原因，弹性轴发生偏移，偏移距离为$ e $，弹性芯变为$ 2c $，内、外表面的应变分别为$ {\varepsilon _{\text{0}}} $和$ {\varepsilon _{\text{i}}} $，内外表面应变对应的结构应力分别为$ {\sigma '_{\text{O}}} $和$ {\sigma '_{\text{i}}} $. 根据弹性状态和塑性状态下的力和弯矩等效，列出

$$ \left\{\begin{array}{l}{\displaystyle \sum F=0,\begin{array}{cc}& {\sigma }_{\text{m}}t={S}_{\text{Y}}\left(\dfrac{t}{2}-(c-e)\right)-{S}_{\text{Y}}\left(\dfrac{t}{2}-(c + e)\right)\end{array}}\\ {\displaystyle \sum {M}_{\text{O}}=0,\begin{array}{cc}& \dfrac{{t}^{2}{\sigma }_{\text{b}}}{6}={S}_{\text{Y}}\left(\dfrac{t}{2}-(c-e)\right)\times \dfrac{1}{2}\left(\dfrac{t}{2} + (c-e)\right) + \dfrac{(2c{)}^{2}{S}_{\text{Y}}}{6} + {S}_{\text{Y}}\left(\dfrac{t}{2}-(c + e)\right)\times \dfrac{1}{2}\left(\dfrac{t}{2} + (c + e)\right)\end{array}}\end{array}\right. $$

(1)

式中：$ {\sigma _{\text{m}}} $和$ {\sigma _{\text{b}}} $为弹性状态下计算的膜应力和弯曲应力；$ {S_{\text{Y}}} $为材料的屈服强度；$ P $为拉伸载荷，$ M $为弯曲载荷；$ t $为板厚；$ e $为弹性轴偏离中性层的距离；$ c $为发生塑性变形后的中性层与塑性区距离；$ {\varepsilon _{\text{i}}} $和$ {\varepsilon _{\text{0}}} $为内、外表面的应变；$ {\sigma '_{\text{O}}} $和$ {\sigma '_{\text{i}}} $为产生内、外表面$ {\varepsilon _{\text{i}}} $和$ {\varepsilon _{\text{0}}} $时的弹性结构应力.求解后可获得式(2)中的弹性轴偏移距离和弹性芯

$$ \left\{ \begin{gathered} e = \dfrac{{{\sigma _{\text{m}}}t}}{{2{S_{\rm{Y}}}}} \\ 2c = t\sqrt {3\left[ {1 - {{\left( {\dfrac{{{\sigma _{\text{m}}}}}{{{S_{\text{Y}}}}}} \right)}^2} - \dfrac{2}{3}\dfrac{{{\sigma _{\text{b}}}}}{{{S_{\text{Y}}}}}} \right]} \\ \end{gathered} \right. $$

(2)

上述计算适用于线性硬化材料或幂硬化材料，弹性芯参数$ 2c $是产生塑性变形后留下的弹性芯大小，中性轴移轴参数$ e $是产生塑性变形以后中性轴的偏移量. 在弹性芯参数$ 2c $和中性轴移轴公式$ e $的计算公式中包含的到$ {\sigma _{\text{m}}} $和$ {\sigma _{\text{b}}} $，两者需要通过数值分析计算得到.假设弹性核的存在占主导地位，根据材料力学的几何关系，弯曲变形的曲率为

$$ \dfrac{1}{R} = \dfrac{M}{{ZI}} = \dfrac{{\dfrac{1}{2}{S_{\text{Y}}} \times c \times \dfrac{2}{3}c \times 2}}{{Z \times \dfrac{{{{(2c)}^3}}}{{12}}}} = \dfrac{{\dfrac{2}{3}{c^2}{S_{\text{Y}}}}}{{\dfrac{2}{3}Z{c^3}}} = \dfrac{{{S_{\text{Y}}}}}{{Zc}} $$

(3)

式中：$ R $为弯曲的曲率半径；$ Z $为弹性模量；$ I $为惯性矩.

发生塑性变性后，中性轴偏移后曲率半径也跟着变化，如图3所示.

图 3 变形关系

Figure 3. The relation of deformation

下载: 全尺寸图片幻灯片

根据图3中的变形关系，可以计算出最外层的应变$ {\varepsilon _{\text{0}}} $和最内层的应变$ {\varepsilon _{\text{i}}} $，即

$$ \begin{gathered} \dfrac{R}{{R + e + \dfrac{t}{2}}} = \dfrac{1}{{1 + 1 \times {\varepsilon _{\text{o}}}}} \Rightarrow {\varepsilon _{\text{o}}} = \dfrac{1}{R}\left(e + \dfrac{t}{2}\right) \\ \dfrac{{\dfrac{t}{2} - e}}{R} = \dfrac{{1 - (1 + {\varepsilon _{\text{i}}})}}{1} \Rightarrow {\varepsilon _{\text{i}}} = \dfrac{1}{R}\left(e - \dfrac{t}{2}\right) \\ \end{gathered} $$

(4)

再根据最外层应变和最内层应变计算膜应变和弯曲应变，即

$$ \begin{gathered} {\varepsilon _{\text{m}}} = \dfrac{{{\varepsilon _{\text{o}}} + {\varepsilon _{\text{i}}}}}{2} = \dfrac{e}{R} \\ {\varepsilon _{\text{b}}} = \dfrac{{{\varepsilon _{\text{o}}} - {\varepsilon _{\text{i}}}}}{2} = \dfrac{t}{{2R}} \\ \end{gathered} $$

(5)

式中：$ {\varepsilon _{\text{m}}} $为膜应变；$ {\varepsilon _{\text{b}}} $为弯曲应变.

获得膜应变$ {\varepsilon _{\text{m}}} $和弯曲应变$ {\varepsilon _{\text{b}}} $后，将两者求和获得结构应变$ {\varepsilon _{\text{s}}} $，同等效结构应力计算方法一样，在循环载荷作用下，等效结构应变的参数为

$$ \Delta {E_{\text{S}}} = \dfrac{{\Delta {\varepsilon _{\text{S}}}}}{{{t^{\tfrac{{2 - m}}{{2m}}}}I{{(r)}^{\tfrac{1}{m}}}}} $$

(6)

式中：$ \Delta {E_{\text{s}}} $等效结构应变的变化范围；$ \Delta {\varepsilon _{\text{s}}} $为结构应变的变化范围；$ m $为应力强度因子的幂指数，取值为3.6；$ r $为弯曲应变$ {\varepsilon _{\text{b}}} $与结构应变$ {\varepsilon _{\text{s}}} $的比值；$ I(r)^{\tfrac{1}{n}} $为载荷比$ r $的无量纲函数^[3].

等效结构应变计算后，结合主$ E - N $曲线，如表1所示，进行寿命$ N $计算，即

表 1 主E-N曲线

Table 1. Main E-N curve

统计依据	试验常数 $ C $	试验常数 $ h $
中值	0.10868	0.3195
$ + 2\sigma $	0.13025
${{ - } }2\sigma$	0.06313
$ + 3\sigma $	0.15610
${{ - } }3\sigma$	0.05268

下载: 导出CSV

| 显示表格

$$ \Delta {E_{\text{S}}} = C \times {N^{ - h}} $$

(7)

式中：$ C $和$ h $为试验常数；$ N $为寿命.

1.1.2 上表面发生屈服

当膜应力较大时，弯曲应力和膜应力只能使板的一侧发生屈服. 对于贯穿厚度的一个截面，使用弹性计算的方式，将得到$ {\sigma _{\text{m}}} $和$ {\sigma _{\text{b}}} $，两者求和后只有外表面超过了屈服强度$ {S_{\text{Y}}} $，如图4所示.

图 4 上表面发生屈服示意图

Figure 4. Schematic diagram of yield on the upper surface relation of deformation

下载: 全尺寸图片幻灯片

根据弹性状态和塑性状态下的力和弯矩等效，得出

$$ \left\{\begin{array}{l}{\displaystyle \sum F=0,\begin{array}{cc}& {\sigma }_{\text{m}}t={S}_{\text{Y}}\left(\dfrac{t}{2}-(c-e)\right)\end{array}}\\ {\displaystyle \sum {M}_{\text{O}}=0,\begin{array}{cc}& \dfrac{{t}^{2}{\sigma }_{\text{b}}}{6}={S}_{\text{Y}}\left(\dfrac{t}{2}-(c-e)\right)\times \dfrac{1}{2}\left(\dfrac{t}{2} + (c\text{ + }e)\right) + \dfrac{(2c{)}^{2}{S}_{\text{Y}}}{6}\end{array}}\end{array}\right. $$

(8)

经计算后，获得中性轴和弹性核，即

$$ \begin{split} & e = \dfrac{1}{8}\dfrac{{{S_{\text{Y}}}t{{\left( { - 3{\sigma _{\text{m}}} - {\sigma _{\text{b}}} + 3{S_{\text{Y}}}} \right)}^2}}}{{{{\left( { - {\sigma _{\text{m}}} + {S_{\text{Y}}}} \right)}^3}}} - \dfrac{1}{2}\dfrac{{\left( { - 2{\sigma _{\text{m}}} - {\sigma _{\text{b}}} + 2{S_{\text{Y}}}} \right)t}}{{ - {\sigma _{\text{m}}} + {S_{\text{Y}}}}} \\& c = \dfrac{t}{2}\left( {\dfrac{{ - 3{\sigma _{\text{m}}} - {\sigma _{\text{b}}} + 3{S_{\text{Y}}}}}{{ - {\sigma _{\text{m}}} + {S_{\text{Y}}}}}} \right) \end{split} $$

(9)

根据材料力学的几何关系，弯曲变形的曲率为

$$ \dfrac{1}{R} = \dfrac{{8{{({S_{\text{Y}}} - {\sigma _{\text{m}}})}^3}}}{{tE{{(3{S_{\text{Y}}} - 3{\sigma _{\text{m}}} - {\sigma _{\text{b}}})}^2}}} $$

(10)

根据曲率公式，计算外表面和内表面的应变为

$$ \begin{gathered} {\varepsilon _{\text{o}}} = \dfrac{{{S_{\text{Y}}}}}{E} + \dfrac{{(t - c)}}{R} \\ {\varepsilon _{\text{i}}} = \dfrac{{{S_{\text{Y}}}}}{E} - \dfrac{c}{R} \\ \end{gathered} $$

(11)

膜应变和弯曲应变、等效结构应变计算公式与式(5)、式(6)一致，寿命计算如式(7)所示. 用有效的屈服应力$ {S'_{\text{Y}}} $替换上述公式中材料的屈服强度$ {S_{\text{Y}}} $为

$$ {S'_{\text{Y}}} = \dfrac{{{S_{\text{Y}}}}}{{\sqrt {1 - \nu + {\nu ^2}} }} $$

(12)

式中：$ \nu $为泊松比.

式(12)将上述公式扩展到平面应变问题. 如果使用冯米塞斯准则，可用下式的$ E' $替换$ E $.

$$ E' = \dfrac{E}{{1 - {\nu ^2}}} $$

(13)

1.1.3 结构应变法的计算实施流程

在弯曲为主的加载式(1)中，将弯矩等效整理后为

$$ {\sigma _{\text{b}}} = \dfrac{{6{S_{\text{Y}}}}}{{{t^2}}}\left[\dfrac{{{t^2}}}{4} - \dfrac{{{c^2}}}{3} - \dfrac{{\sigma _{\text{m}}^2{t^2}}}{{4S_{\text{Y}}^2}}\right]\left( {0 \leqslant c \leqslant \dfrac{t}{2} - e} \right) $$

(14)

如图2所示，当$ c = 0 $时，整个截面均达到塑性，求解式(14)，得

$$ {\sigma _{{\text{bmax}}}} = \dfrac{{6{S_{\text{Y}}}}}{{{t^2}}}\left[\dfrac{{{t^2}}}{4} - \dfrac{{\sigma _{\text{m}}^2{t^2}}}{{4S_{\text{Y}}^2}}\right] = \dfrac{3}{2}{S_{\text{Y}}}\left[ {1 - \dfrac{{\sigma _{\text{m}}^2}}{{S_{\text{Y}}^2}}} \right] $$

(15)

如图4所示，当$ c = 0.5t - e $时，仅截面一侧发生塑性变形，求解式(14)，得

$$ \begin{split} {\sigma _{{\text{bmin}}}}& = \dfrac{{6{S_{\text{Y}}}}}{{{t^2}}}\left[\dfrac{{{t^2}}}{4} - \dfrac{1}{3}{\left(\dfrac{t}{2} - e\right)^2} - \dfrac{{\sigma _{\text{m}}^2{t^2}}}{{4S_{\text{Y}}^2}}\right]\\& = {S_{\text{Y}}}\left[ {1 + \dfrac{{\sigma _{\text{m}}^{}}}{{S_{\text{Y}}^{}}} - \dfrac{{2\sigma _{\text{m}}^{\text{2}}}}{{S_{\text{Y}}^2}}} \right] \end{split}$$

(16)

根据以上计算结果，以弹性计算获得结构应力$ {\sigma _{\text{s}}} = {\sigma _{\text{m}}} + {\sigma _{\text{b}}} $后，整个截面存在以上几种状态：(1)当$ {\sigma _{\text{s}}} \leqslant {S_{\text{Y}}} $时，截面没有发生屈服；(2)当$ {\sigma _{\text{b}}} \leqslant {\sigma _{{\text{bmin}}}} $，且$ {\sigma _{\text{s}}}> {S_{\text{Y}}} $时，外表面都发生屈服；(3)当$ {\sigma _{{\text{bmin}}}} \leqslant {\sigma _{\text{b}}} \leqslant {\sigma _{{\text{bmax}}}} $，且$ {\sigma _{\text{s}}} > {S_{\text{Y}}} $时，内、外表面都发生屈服；(4)当$ {\sigma _{\text{b}}} > > {\sigma _{{\text{bmax}}}} $时，截面发生塑性失效. 基于上述原理程序流程如图5所示.

图 5 焊接结构高、低周疲劳实施流程

Figure 5. Calculation process of high and low cycle fatigue of welded structures

下载: 全尺寸图片幻灯片

2. 试验验证及分析

2.1 试验验证

以Q450NQR1材料的搭接接头试样为研究对象，焊接接头尺寸示意图如图6所示，为保证试验的准确性，使用能够实现塑性应变测量的测试胶水进行应变片粘贴.应变片布置如图7所示，在距离焊趾C截面2 mm、17 mm的B截面、A截面位置的上下表面布置应变片.

图 6 焊接接头宽板试样示意图(mm)

Figure 6. Schematic diagram of wide plate specimen of welded joint

下载: 全尺寸图片幻灯片

图 7 应变片布置位置(mm)

Figure 7. Strain gauge arrangement position

下载: 全尺寸图片幻灯片

值得注意的是由于搭接接头产生了附加弯矩，搭接接头长度为300 mm，两端夹具的夹持长度分别为80 mm，在仿真模型上一定要考虑夹持长度. 另外，由于板的宽度和应变片尺寸限制，导致布置的测点并不在1条直线上，但由于是小试样，不会产生较大影响. 弹性加载下的应变测试结果如表2所示，由于构件焊接过程中存在焊接变形，首次加载后，构件会发生1次变形调整，卸载后再重新加载，测试结果稳定，如图8所示.根据表2中的应变测试结果，采用结构应力的方法进行应力集中系数计算. 第1次拉伸时，试验存在变形微校正，选择拉伸第2次的数据计算应力集中系数为

表 2 弹性加载下的应变测试

Table 2. strain test under elastic loading

载荷类型	加载次数n/次	载荷值P/kN	实测应变 ε(10⁻⁶)
载荷类型	加载次数n/次	载荷值P/kN	1号	2号	3号	4号	5号	6号
拉伸	1	25	714	251	1342	−14	1593	45
	2	25	771	198	802	112	814	122
	3	25	771	197	785	115	811	123
压缩	1	10	330	−61	−320	−37	−336	−35
	2	25	−859	−140	−818	−86	−852	−87
	3	25	−803	−163	−790	−100	−820	−107

下载: 导出CSV

| 显示表格

图 8 焊接接头应变片布置及应变测试

Figure 8. Strain gauge arrangement and strain testing of welded joints

下载: 全尺寸图片幻灯片

$$ \left\{\begin{gathered} {S_1} = 158.8\;\;{\text{MPa}}，\begin{array}{*{20}{c}} {}&{{S_2} = 40.8\;\;{\text{MPa}}} \end{array} \\ {S_3} = 165.2\;\;{\text{MPa}}，\begin{array}{*{20}{c}} {}&{{S_4} = 23.0\;\;{\text{MPa}}} \end{array} \\ \end{gathered}\right. $$

(17)

$$ \left\{\begin{gathered} \sigma _{\text{b}}^{\text{B}} = \dfrac{1}{2}({S_1} - {S_2}) = 59.0\;\;{\text{MPa}} \\ \sigma _{\text{m}}^{\text{B}} = \dfrac{1}{2}({S_1} + {S_2}) = 99.8\;\;{\text{MPa}} \\ \sigma _{\text{b}}^{\text{A}} = \dfrac{1}{2}({S_3} - {S_4}) = 71.1\;\;{\text{MPa}} \\ \end{gathered}\right. $$

(18)

$$\begin{split} \sigma _{\text{s}}^{\text{C}}& = \sigma _{\text{m}}^{\text{C}} + \sigma _{\text{b}}^{\text{C}} = \sigma _{\text{m}}^{\text{A}} + \sigma _{\text{b}}^{\text{A}} + \dfrac{{17}}{{15}}(\sigma _{\text{b}}^{\text{B}} - \sigma _{\text{b}}^{\text{A}})\\& { = {S_1}} + \dfrac{{17}}{{15}}(\sigma _{\text{b}}^{\text{B}} - \sigma _{\text{b}}^{\text{A}}) \\ &{ = 158.8 + } \dfrac{{17}}{{15}}(71.1 - 59.0)\\&{ = 172.54\;\;{\text{MPa}}} \end{split} $$

(19)

通过$ \sigma _{\text{s}}^{\text{C}} $和名义应力比值(25 kN下的名义应力为100 MPa)计算出结构应力集中系数为1.73. 通过试验机采用正弦曲线将试样加载到95 kN载荷(名义应力为380 MPa)，实测的应变数据如图9所示.

图 9 实测的应变曲线

Figure 9. Measured strain curves

下载: 全尺寸图片幻灯片

按实测的应力集中系数1.73，380 MPa载荷下弯矩引起的最大应力为273.6 MPa. 根据上表面屈服的结构应变的式(8) ~ 式(11)可得：上表面焊趾的微应变为2910，焊趾下表面为460. 膜应变和弯曲应变分别为1685和1225. 弹性模量E为206000，$ m $为3.6，板厚$ t $为5 mm. 则弯曲比$ r $，载荷控制下$I{(r)^{\tfrac{1}{m}}} $和等效结构应变$ \Delta {E_{\text{s}}} $分别为

$$ r = \dfrac{{{\varepsilon _{\text{b}}}}}{{{\varepsilon _{\text{s}}}}} = \dfrac{{1\;225}}{{2\;910}} = 0.420\;9 $$

(20)

$$ I{(r)^{\tfrac{1}{m}}} = \dfrac{{1.23 - 0.364r - 0.17{r^2}}}{{1.007 - 0.306r - 0.178{r^2}}} = 1.236\;2 $$

(21)

$$ \Delta {E_s} = \dfrac{{\Delta {\varepsilon _s}}}{{{t^{\tfrac{{2 - m}}{{2m}}}}I{{(r)}^{\tfrac{1}{m}}}}} = \dfrac{{2\;910}}{{0.699\;3 \times 1.236\;2}} = 3\;366 $$

(22)

根据式(7)计算了该试样在不同统计值下的寿命，并与疲劳试验结果进行了对比，结果如表3所示. 可见，中值下的低周疲劳寿命计算结果与试验值吻合.

表 3 寿命计算

Table 3. Life calculation

寿命计算值(次)					寿命试验值 (次)
中值	( + 2σ)	(−2σ)	( + 3σ)	(−3σ)	寿命试验值 (次)
34780	149984	8065	646779	1870	24201

下载: 导出CSV

| 显示表格

2.2 适用性分析

建立的平面应变有限元模型如图10所示，模型一端施加对称约束，另一端施加膜应力和弯曲应力载荷，在靠近对称约束位置，建立路径，路径的起点和终点如图10所示.

图 10 平面应变模型

Figure 10. Plane strain finite element model

下载: 全尺寸图片幻灯片

理论算法和平面应变的有限元模型对比如图11所示. 固定端的膜应力载荷$ {\sigma _{\text{m}}} $取为250 MPa，300 MPa，350 MPa和380 MPa，通过改变端部弯曲应力载荷，计算不同弯曲应力载荷下平面应变模型和结构应变理论解对比如图11(a) ~ 图11(d)所示.

图 11 理论算法和平面应变的有限元模型对比

Figure 11. Comparison between theoretical algorithm and plane strain finite element model. (a) the membrane stress is 250 MPa; (b) the membrane stress is 300 MPa; (c) the membrane stress is 350 MPa; (d) the membrane stress is 380 MPa; (e) the bending stress is 250 MPa; (f) the bending stress is 325 MPa

下载: 全尺寸图片幻灯片

固定端的弯曲应力载荷$ {\sigma _b} $取250 MPa、325 MPa，通过改变端部膜应力载荷，计算不同膜应力载荷. 平面应变模型和结构应变理论解对比如图11(e)、图11(f)所示.计算结果表明：当材料的屈服为550 MPa，当弯曲载荷与膜应力载荷合载荷为700 MPa时，平面应变模型和结构应变理论解结果完全一致，高于700 MPa后，随着载荷的增加，误差也随着增加.

3. 应用实例

以货车车体为研究对象，采用板壳单元建模，单元长度为30 mm，将车体考虑为柔性体，采用固定交界面法完成柔性体计算^[12-13]. 为了能够施加与实际状态更为吻合的载荷，以试验台模型作为车体仿真的边界条件^[14]. 货车虚拟试验模型如图12所示.在虚拟试验台的车钩加载端，施加制动等效压缩载荷，载荷达到峰值后的车体的最大主应力云图如图13所示，主应力最大的位置集中在车体底盘，在车钩和两侧过度位置主应力集中较为严重.一般来讲，主应力与结构正应力往往有正相关的关系，因此我们根据车体最大主应力云图对关心焊缝进行筛选. 提取了如图13所示的A₁和B₁位置的正应力，发现A₁位置明显超过Q450NQR1钢的屈服极限，B₁位置与屈服极限接近，如图14所示.通过计算的节点力先计算结构应力，再根据2.1节中的计算流程计算出结构应变，对比了结构应力法和结构应变的损伤对比结果，如图15所示.

图 12 货车虚拟试验模型

Figure 12. Virtual test model of fast freight car

下载: 全尺寸图片幻灯片

图 13 变形及最大主应力云图

Figure 13. Cloud map of deformation and maximum principal stress. (a) first side of car body; (b) second side of car body

下载: 全尺寸图片幻灯片

图 14 正应力曲线

Figure 14. Normal stress curves

下载: 全尺寸图片幻灯片

图 15 基于结构应力和结构应变所得损伤对比

Figure 15. Comparison of damage based on structural stress and structural strain. (a) The position of A₁; (b) The position of B₁

下载: 全尺寸图片幻灯片

A₁位置 ~ A₄位置由于结构正应力超出材料比例极限较多，因此塑性变形比较明显，这导致了基于结构应变计算的车体疲劳损伤超过基于结构应力计算的车体疲劳损伤，以A₁位置为例，基于结构应变计算的车体疲劳损伤比基于结构应力计算的车体疲劳损伤高2倍.而在车钩附近的B₁位置 ~ B₄位置，由于结构正应力超出材料屈服强度不多，因此基于结构应变计算得到疲劳损伤仅仅略大于基于结构应力计算得到的疲劳损伤.虽然近几年针对低周疲劳问题，针对更为精确的材料应力应变性能曲线，提出了更为准确的计算方法，但由于理想弹塑性材料不需要材料的力学性能曲线，所以开展的结构应变法在工程具有很大的应用价值.

4. 结论

(1)为了解决铁路货车的低周疲劳评估问题，详细梳理低周疲劳的结构应变法理论，完成了程序设计和方法的适用性验证.

(2)平面应变有限元模型和结构应变法理论解的对比结果表明，基于理想弹性性模型的结构应变方法及程序计算的Q450NQR1材料，当结构应力高出屈服强度150 MPa以内时，平面应变模型和结构应变理论解结果完全一致，超过该范围后，随着载荷的增加，误差也增加.

(3)低周疲劳的结构应变法现阶段能够满足铁路货车的设计需求，研究为该方法的工程推广应用提供了良好的技术支撑.

[1]

蒙世合, 罗海辉, 彭宇. WC晶粒离散度对WC-Co硬质合金断裂韧性的影响[J]. 硬质合金, 2017, 34(1): 14-20.Meng Shihe, Luo Haihui, Peng Yu. Effect of WC grain dispersion on fracture toughness of WC-Co cemented carbide[J]. Cemented Carbide, 2017, 34(1): 14-20.[2] 羊建高, 熊继. 稀土硬质合金的研究现状及发展趋势[J]. 稀土, 1992, 13(4): 45-47.Yang Jiangao, Xiong Ji. Research status and development trend of rare earth cemented carbide[J]. Chinese Rare Earths, 1992, 13(4): 45-47.[3] 王悦悦, 殷国涛, 胡小小, 等. 硬质合金/不锈钢激光焊焊接性分析[J]. 热加工工艺, 2017(5): 200-203.Wang Yueyue, Yin Guotao, Hu Xiaoxiao,et al. Laser welding weldability of cemented carbide/stainless steel[J]. Hot Working Technology, 2017(5): 200-203.[4] 李远星, 张晓山, 朱宗涛, 等. Ni元素扩散行为对硬质合金/钢钎焊接头微观组织及力学性能的影响[J]. 稀有金属材料与工程, 2017， 46(4): 1120-1125.Li Yuanxing, Zhang Xiaoshan, Zhu Zongtao,et al. Effect of element Ni diffusion on the mechanical properties and microstructure of joints brazed of cemented carbide and steel[J]. Rare Metal Materials and Engineering, 2017, 46(4): 1120-1125.[5] 付文强, 王海龙. 硬质合金(YG6C)/钢(16Mn)复合螺柱的螺柱焊横焊工艺研究[J]. 热加工工艺, 2017(3): 65-70.Fu Wenqiang, Wang Hailong. Research on horizontal position stud welding of compound stud of hard alloy YG6C/16Mn steel[J]. Hot Working Technology, 2017(3): 65-70.[6] 朱警雷, 黄继华, 张华, 等. 硬质合金与钢异种金属焊接的研究进展[J]. 焊接, 2008(2): 15-19.Zhu Jinglei, Huang Jihua, Zhang Hua,et al. Progress in research on welding of cemented carbide with steel[J]. Welding & Joining, 2008(2): 15-19.[7] Xu P Q, Zhao X J, Yang D X,et al. Study on filler metal (Ni-Fe-C) during GTAW of WC-30Co to 45 carbonsteel[J]. Journal of Materials Science, 2005, 40(24): 6559-6564.[8] 陈国庆, 张秉刚, 吴振中, 等. 硬质合金与钢电子束焊接接头组织及性能[J]. 焊接学报, 2013, 34(6): 9-12.Chen Guoqing, Zhang Binggang, Wu Zhenzhong,et al. Microstructure and mechanical properties of cemented carbide to steel joints by electron beam welding[J]. Transactions of the China Welding Institution, 2013, 34(6): 9-12.[9] 赵秀娟, 王浩, 杨德新, 等. 焊丝成份对YG30/钢焊缝组织和硬度的影响[J]. 硬质合金, 2004, 21(2): 94-98.Zhao Xiujuan, Wang Hao, Yang Dexin,et al. The effect of the welding wire composition on the microstructure and hardness of the welding joints between YG30 and steel[J]. Cemented Carbide, 2004, 21(2): 94-98.[10] Zhao X J, Yang D X. Theηphases and mechanical properties of TIG welded joints of WC-Co cemented carbide and steel[J]. China Welding, 2004, 13(1): 34-39.

期刊类型引用(2)

1.	李传胜，王雷，郭糠，杨祥帆，于文晶，王利忠，王洪宝，张怀强，张英波. Mo对钛合金激光焊接接头组织性能影响研究. 电焊机. 2025(01): 73-79 . 百度学术
2.	李海英，熊淑秋. 船用耐蚀材料激光焊接技术分析. 舰船科学技术. 2024(21): 67-70 . 百度学术