高级检索

选区激光熔化制备多孔结构的成形偏差及力学性能与压缩失效分析

徐榕蔚, 张振杰, 刘清原, 张光辉, 龙芋宏

徐榕蔚, 张振杰, 刘清原, 张光辉, 龙芋宏. 选区激光熔化制备多孔结构的成形偏差及力学性能与压缩失效分析[J]. 焊接学报, 2022, 43(10): 49-56. DOI: 10.12073/j.hjxb.20211005001
引用本文: 徐榕蔚, 张振杰, 刘清原, 张光辉, 龙芋宏. 选区激光熔化制备多孔结构的成形偏差及力学性能与压缩失效分析[J]. 焊接学报, 2022, 43(10): 49-56. DOI: 10.12073/j.hjxb.20211005001
XU Rongwei, ZHANG Zhenjie, LIU Qingyuan, ZHANG Guanghui, LONG Yuhong. The forming deviation, mechanical properties and compression failure of porous structures fabricated by laser melting were analyzed[J]. TRANSACTIONS OF THE CHINA WELDING INSTITUTION, 2022, 43(10): 49-56. DOI: 10.12073/j.hjxb.20211005001
Citation: XU Rongwei, ZHANG Zhenjie, LIU Qingyuan, ZHANG Guanghui, LONG Yuhong. The forming deviation, mechanical properties and compression failure of porous structures fabricated by laser melting were analyzed[J]. TRANSACTIONS OF THE CHINA WELDING INSTITUTION, 2022, 43(10): 49-56. DOI: 10.12073/j.hjxb.20211005001

选区激光熔化制备多孔结构的成形偏差及力学性能与压缩失效分析

基金项目: 国家自然科学基金青年科学基金资助项目(62004050, 52165056);广西自然科学基金重点项目(2019JJD160010);广西研究生教育创新项目(2020YCXS010, 2021YCXS001);2020年全国大学生创新创业训练计划(No.202010595032).
详细信息
    作者简介:

    徐榕蔚,硕士;主要研究方向为金属增材制造与结构优化设计. Email: 974613097@qq.com

    通讯作者:

    龙芋宏,博士,教授,博士研究生导师;Email: longyuhong@guet.edu.cn

  • 中图分类号: TG495,TF124

The forming deviation, mechanical properties and compression failure of porous structures fabricated by laser melting were analyzed

  • 摘要: 因多孔结构轻质高强度、力学性能可调节的特点,被广泛用于骨骼医疗、航空航天等领域. 为了探索多孔结构选区激光熔化(Selective Laser Melting, SLM)成形误差与压缩失效性能,以钻石型晶格和六孔开口球形两种多孔结构为例,采用理论预测与试验测试研究SLM制造多孔结构的压缩力学行为,使用ANSYS软件对所研究的多孔结构进行准静态压缩模拟,并对SLM成形的多孔结构进行单轴压缩试验,最后结合仿真和试验,观测和分析它们的变形过程和失效机制. 对比后发现数值设计的多孔结构尺寸与最终制造的结构存在偏差,导致力学性能理论值与试验值存在一定差异,但应力应变场变化规律一致. 试验结果表明,在孔隙率50% ~ 80%时,钻石型晶格结构屈服强度为31.85 ~ 182.13 MPa,弹性模量为1.45 ~ 2.30 GPa;六孔开口球形结构屈服强度为35.19 ~ 130.64 MPa,弹性模量为1.59 ~ 2.90 GPa,不同多孔结构随孔隙率的增大,力学性能变化趋势不一致.
    Abstract: Due to the characteristics of light, high strength and adjustable mechanical properties of porous structure, it is widely used in bone medicine, aerospace and other fields. In order to explore the forming error and compression failure performance of porous structure with selective laser melting (SLM), this paper takes two kinds of porous structure with diamond lattice and spherical six-hole opening as examples to study the compressive mechanical behavior of porous structure manufactured by SLM by theoretical prediction and experimental test. ANSYS software was used to simulate the quasi-static compression of the studied porous structure, and the uniaxial compression experiment of the SLM formed porous structure was carried out. Finally, the deformation process and failure mechanism of the SLM formed porous structure were observed and analyzed combined with the simulation and experiment. After comparison, it is found that the size of the numerical design porous structure deviates from that of the final manufactured structure, resulting in a certain difference between the theoretical value of mechanical properties and the experimental value, but the variation law of stress and strain field is consistent. The experimental results show that when the porosity is 50% ~ 80%, the yield strength and elastic modulus of diamond lattice structure are 31.85 ~ 182.13 MPa and 1.45 ~ 2.30 GPa respectively. The yield strength and elastic modulus of six-hole spherical structure are 35.19 ~ 130.64 MPa and 1.59 ~ 2.90 GPa respectively. The mechanical properties of different porous structures vary with the increase of porosity.
  • 不锈钢轨道车辆以其高耐腐蚀性能、较小的车体自身质量、较低的维修费用和运营成本等显著优点成为国内外轨道交通轻量化车体的主流,而不锈钢车体,通常使用电阻点焊连接,每辆车的点焊焊点多达3万个[1]. 点焊作为一种经济、快速的连接方式,具有高效率、低能耗、高产出等优点,已成为地铁车辆、城市轻轨车辆、甚至高速动车组不锈钢车体结构金属板件之间的主要连接形式[2]. 由于点焊焊核处在封闭状态,而疲劳裂纹常常萌生于板的内侧,使得对点焊进行无损探伤以及裂纹长度和深度的监测比较困难,而通过疲劳试验得到的一般都是载荷—寿命(ΔF-N)曲线,因载荷模式、板厚组合、焊核直径不同而不同,不利于工程上的应用[3]. 因此,建立一种合理表征点焊疲劳特性的方法,对预测和评估点焊疲劳寿命有重要意义.

    点焊疲劳寿命的预测方法主要包括名义应力法、热点应力法、缺口应力法、断裂力学法、局部法和结构应力法等. 其中大多数方法都存在对有限元网格敏感性大、适用范围小等局限,而结构应力法由于对网格不敏感,得到了较广泛应用. 现行基于结构应力的点焊疲劳寿命预测方法主要有Radaj法[4]、Rupp法[5]、Sheppard法[6-7]、Kang法[8]和Dong法[9]等. 早期,Radaj[4]根据点焊试验的外载荷、焊核直径和板厚来计算拉剪载荷下点焊试件焊核边缘处的最大结构应力. Rupp[5]提出基于结构应力方法预测点焊寿命,并通过建立壳单元和梁单元有限元模型仿真得到点焊结构应力,结果与真实应力较为接近. Dong[9]将承受复杂载荷的焊接结构应力等效地转化为用膜应力和弯曲应力表示的简单结构应力,建立了一条考虑各种影响因素的主S-N曲线. 闫坤[10]结合Dong法和Rupp法,提出了一种新的点焊疲劳特征参量和简化的结构应力,并评估了基于Kang结构应力、Sheppard结构应力和等效应力强度因子等6种特征参量下的点焊疲劳性能. 吴骁[11]提出了一种对Rupp结构应力计算式中的经验参数项进行优化的方法,并基于该方法得到了相关性较高的点焊S-N曲线. 谢素明等人[12]基于结构应力法分析了不锈钢车体中主横梁与底架边梁区域的点焊接头结构应力分布规律,并对焊点进行了寿命预测.

    尽管围绕点焊的结构应力计算、有限元模型简化、疲劳寿命预测方法、结构应力优化方法的研究已取得一系列进展,但工程中对点焊结构进行疲劳分析时,仍采用ΔF-N曲线或通过建立复杂的有限元模型来进行寿命预测. 然而,不锈钢车体等实际结构往往以不同板厚、焊核直径、焊点数的组合点焊形式呈现,如果对每种形式都进行疲劳试验或有限元建模,过程繁琐且费时费力. 因此,需要简化点焊模型,对不同ΔF-N曲线进行归一化处理,建立综合考虑载荷比和几何参数等因素的点焊S-N曲线,以及适合工程应用的P-S-N曲线.

    文中基于广泛使用的Rupp结构应力法,研究点焊接头的结构应力,并根据疲劳试验寿命拟合得到点焊结构的S-N曲线,然后采用Goodman修正法将非对称循环下的结构应力当量为对称循环下的结构应力来消除应力比的影响,最终建立点焊的P-S-N曲线,使之更具工程参考价值.

    文章针对轨道车辆不锈钢车体结构中常见的点焊连接类型,设计点焊试件进行拉剪疲劳试验. 其中点焊试件的结构如图1所示,为了防止试验过程中因试件板厚不同而导致加载偏心,在夹持试件的端部焊接了辅助板,保证了夹持厚度一致和焊接接头受力均匀. 所有点焊试件的几何参数及载荷比见表1.

    图  1  不锈钢点焊试件结构简图
    Figure  1.  Schematic diagram of stainless steel spot welded sample
    表  1  点焊试件几何参数及载荷比
    Table  1.  Geometric parameters and load ratios of spot welded samples
    试件板厚组合s/mm焊核直径d/mm总长度L/mm重叠长度l/mm板宽D/mm载荷比r有效试样材料
    s1s2
    1.52.06.025025500.120SUS301L,1/8H + 1/8H
    1.53.06.025025500.120SUS301L,1/8H + 1/16H
    2.02.07.025025500.120SUS301L,1/4H + 1/4H
    2.02.07.025025500.120SUS301L,1/8H + 1/8H
    0.81.24.531025300.59SUS301L + C850
    0.82.54.531025300.58
    2.02.06.031025300.59
    1.52.06.031025300.57
    1.51.56.031025300.19
    1.51.56.031025300.510
    下载: 导出CSV 
    | 显示表格

    试件在RUMUL高频疲劳试验机上进行拉剪疲劳试验,以正弦波加载,试验有效频率范围为50~100 Hz,由于从观测到点焊试件产生裂纹到试件完全断裂的时间很短,为方便记录试验结果,以点焊试件断裂作为失效评判依据,最终试验ΔF-N数据如图2所示. 由图可知,由于受板厚、焊核直径、载荷形式等的影响,ΔF-N数据分散性较大,拟合曲线线性相关系数平方R2仅为0.695 3. 如果针对不同的板厚、开孔直径组合给出不同的预测曲线,选取过程又将变得繁琐,因此绘制适用于不同板厚、焊核直径和载荷类型的点焊S-N曲线,更具有工程意义.

    图  2  点焊试件ΔF-N散点图
    Figure  2.  ΔF-N scatter diagram of spot welded samples

    点焊结构应力计算方法中最为经典的是Rupp提出的等效结构应力法,如图3所示,Rupp评估法[5]对于点焊疲劳强度的预测主要包括对焊核和焊核周围母材区这两部分,对于厚板($d > 3.5\sqrt s $)和小直径焊核点焊结构易在焊核处产生裂纹[13],焊核处的结构应力主要由正应力和切应力构成. 对于薄板($d < 3.5\sqrt s $)和大焊核直径的点焊结构易在母材区产生裂纹,母材区的结构应力主要由径向应力决定.

    图  3  典型的点焊结构模型
    Figure  3.  Typical model of spot welded structure. (a) base material crack; (b) weld nugget crack

    根据结构应力计算模型,当板厚满足$d > 3.5\sqrt s $时,Rupp得到母材区的结构应力为[5]

    $$\begin{aligned} {\sigma _{v1}} =& - {\sigma _{\max }}\left( {{F_{x1}}} \right)\cos \theta - {\sigma _{\max }}\left( {{F_{y1}}} \right)\sin \theta + {\sigma _{\max }}\left( {{F_{{\textit{z}}1}}} \right)+ \\& {\sigma _{\max }}\left( {{M_{x1}}} \right)\sin \theta - {\sigma _{\max }}\left( {{M_{y1}}} \right)\cos \theta \\ \end{aligned} $$ (1)

    其中,

    $${\sigma _{\max }}\left( {{F_{x1}}} \right){\rm{ = }}\frac{{{F_{x1}}}}{{{\text{π}} d{s_1}}}$$ (2)
    $${\sigma _{\max }}\left( {{F_{y1}}} \right){\rm{ = }}\frac{{{F_{y1}}}}{{{\text{π}} d{s_1}}}$$ (3)
    $${\sigma _{\max }}\left( {{M_{x1}}} \right){\rm{ = }}{K_1}\left( {\frac{{1.872{M_{x1}}}}{{ds_1^2}}} \right)$$ (4)
    $${\sigma _{\max }}\left( {{M_{y1}}} \right){\rm{ = }}{K_1}\left( {\frac{{1.872{M_{y1}}}}{{ds_1^2}}} \right)$$ (5)

    式中:${F_{{\textit{z}}1}} > 0$为焊核直径(mm);${F_{{\textit{z}}1}} > 0$为板1的厚度(mm);${F_{{\textit{z}}1}} > 0$为经验修正系数;${F_{{\textit{z}}1}} > 0$为所计算应力的角度.

    当焊核所受轴向力为拉力时,会导致疲劳破坏,需考虑轴向力的影响,即当${F_{{\textit{z}}1}} > 0$

    $${\sigma _{\max }}\left( {{M_{y1}}} \right){\rm{ = }}{K_1}\left( {\frac{{1.872{M_{y1}}}}{{ds_1^2}}} \right)$$ (6)

    当焊核所受轴向力为压力时,不会产生疲劳,此时忽略轴向力的作用,即当${F_{{\textit{z}}1}} \leqslant 0$

    $$\sigma \left( {{F_{{\textit{z}}1}}} \right){\rm{ = }}0$$ (7)

    当点焊接头满足$d < 3.5\sqrt s $时焊核处的结构应力计算公式为[4, 14]

    $$ \tau = {\tau _{\max }}\left( {{F_{x3}}} \right)\sin \theta + {\tau _{\max }}\left( {{F_{y3}}} \right)\cos \theta $$ (8)
    $$ \sigma = \sigma \left( {{F_{{\textit{z}}3}}} \right) + {\sigma _{\max }}\left( {{M_{x3}}} \right)\sin \theta - {\sigma _{\max }}\left( {{M_{y3}}} \right)\cos \theta $$ (9)

    其中,

    $${\tau _{\max }}\left( {{F_{x3}}} \right) = \frac{{16{F_{x3}}}}{{3{\text{π}} {d^2}}}$$ (10)
    $${\tau _{\max }}\left( {{F_{y3}}} \right) = \frac{{16{F_{y3}}}}{{3{\text{π}} {d^2}}}$$ (11)

    $ {F_{{\textit{z}}3}} > 0$

    $$\sigma \left( {{F_{z3}}} \right) = \frac{{4{F_{z3}}}}{{{\text{π}} {d^2}}}$$ (12)

    ${F_{{\textit{z}}3}} \leqslant 0$

    $$\sigma \left( {{F_{{\textit{z}}3}}} \right) = 0$$ (13)
    $${\sigma _{\max }}\left( {{M_{x3}}} \right) = \frac{{32{M_{x3}}}}{{{\text{π}} {d^3}}}$$ (14)
    $${\sigma _{\max }}\left( {{M_{y3}}} \right) = \frac{{32{M_{y3}}}}{{{\text{π}} {d^3}}}$$ (15)

    然后根据应力分量,可以计算出焊核的绝对最大主应力${\sigma _{\rm{AMP}}}$

    $\sigma \geqslant 0$

    $${\sigma _{\rm{AMP}}} = \frac{\sigma }{2} + \sqrt {{{\left( {\frac{\sigma }{2}} \right)}^2} + {\tau ^2}} $$ (16)

    $ \sigma < 0$

    $${\sigma _{\rm{AMP}}} = \frac{\sigma }{2} - \sqrt {{{\left( {\frac{\sigma }{2}} \right)}^2} + {\tau ^2}} $$ (17)

    分别建立如图4所示的点焊试件体单元和壳单元有限元模型,焊核分别采用体单元和CBAR单元模拟. 通过准静态拉伸试验和有限元仿真对比来验证有限元模型的准确性,准静态试验载荷为Fmin = 611 N,Fmax = 6 111 N,则r为0.1. 应变片贴片位置如图5所示,在距变形环外缘2 mm处[15].

    图  4  点焊试件有限元模型
    Figure  4.  Finite element model of spot welded sample. (a) solid element model; (b) shell element model
    图  5  准静态试验中应变片贴片位置
    Figure  5.  Position of the strain gauge on the sample in the quasi-static test

    试验与仿真结果对比如表2所示[15],可以看出体单元模型仿真所得应变结果与试验测得应变值十分接近,而壳单元模型中梁单元的节点力与体单元模型焊核横截面的力非常接近,证明采用两种有限元模型均可以较为准确地模拟点焊,为了方便提取点焊结构的受力,文中采用点焊壳单元模型来模拟点焊试件.

    表  2  体单元模型仿真结果与试验结果对比
    Table  2.  Comparison of simulation result of solid element model and test result
    提取结果应变剪切力Fσ/N轴向力FN/N
    试验518
    体单元模型5156111139.2
    壳单元模型6110.2140.47
    下载: 导出CSV 
    | 显示表格

    将壳单元有限元模型加载后,导入Nastran软件进行仿真,得到各类型试件在最大载荷Fmax和最小载荷Fmin下梁单元节点的力和力矩,然后带入结构应力公式计算结构应力. 文中点焊试件板厚满足$d > 3.5\sqrt s $,故将得到的梁单元节点的力和力矩带入式(1) ~ 式(7),得到拉剪载荷下的点焊试件的等效结构应力.

    在双对数坐标系下,以试验的疲劳寿命N为横坐标,以等效结构应力幅值Δσs为纵坐标,基于最小二乘法对表1中材料为SUS301L + C850的点焊试验数据进行线性回归拟合,得到点焊试件不区分及区分应力比r(载荷比)的S-N曲线如图6所示. 在疲劳数据处理中,一般以5倍寿命范围来评价S-N曲线的好坏,如果大多数疲劳数据位于5倍寿命范围之内,则认为该S-N曲线关联性较强,可以用来预测疲劳寿命[11, 16].

    图  6  材料为SUS301L + C850的点焊试件的S-N曲线
    Figure  6.  S-N curves of spot welded samples made of SUS301L + C850. (a) S-N curve of spot welded samples under all stress ratios; (b) S-N curve of spot welded samples under stress ratio r = 0.1; (c) S-N curve of spot welded samples under stress ratio r = 0.5

    由图可知,考虑所有应力比的S-N曲线R2 = 0.658 2,相关性较差,个别数据点位于5倍寿命范围之外. 应力比r = 0.1和0.5下的S-N曲线R2分别为0.824 5和0.848 2,所有的数据点都位于5倍寿命范围之内,相关性较好. 可见,应力比对点焊试件的S-N曲线影响较大,因此分析点焊试件在拉剪状态下的疲劳特性时,不可忽略应力比的影响. 在实际的应用过程中应当对不同应力比的S-N曲线加以区分.

    将试验数据根据应力均值Sm大小划分为三级,得到不同应力均值水平下点焊的S-N曲线如图7所示. 可以看出应力均值水平越高,S-N曲线越低,点焊的疲劳强度越小. 对比图6b图6c可知,应力比r越大,S-N曲线越低,即Sm越大,点焊的疲劳强度越小,因为rSm相关联且满足[17]

    $${S_{\rm{m}}} = (1 + r){S_{\rm{a}}}/(1 - r)$$ (18)

    式中:Sa为应力幅值.

    图  7  不同应力均值水平下点焊试件的S-N曲线
    Figure  7.  S-N curves of spot welded samples at different stress mean levels

    为了消除应力比对点焊试件S-N曲线的影响,需要对同种材料不同应力比下点焊试件的结构应力进行修正,从而得点焊相关性更高的基本S-N曲线(r = −1). 文中点焊试验得到的Sa-Sm曲线(应力幅—应力均值曲线)趋于Goodman直线,故采用Goodman修正法将非对称循环下的结构应力当量为对称循环下的结构应力. 以材料为SUS301L + C850的点焊试件为例,将应力比为r = 0.1和0.5的结构应力修正为r = −1时的结构应力,最后得到点焊试件的基本S-N曲线如图8所示. 根据同样的方法可以将其他应力比下的结构应力通过Goodman修正法等效为r = −1时的结构应力,从而消除应力比对点焊试件S-N曲线影响. 对比图6a图8可以看出,修正后点焊的基本S-N曲线线性相关程度(R2 = 0.729 6)较未修正的所有应力比下的S-N曲线(R2 = 0.658 2)高.

    图  8  基于Goodman修正法修正后点焊试件的基本S-N曲线
    Figure  8.  Basic S-N curve of spot welded samples modified by Goodman correction method

    点焊S-N曲线是对点焊结构进行疲劳寿命预测的前提和基础. 在实际疲劳分析中,尤其是构件的可靠性设计中,往往需要使用不同存活率的S-N曲线,即P-S-N曲线. 文中基于表1中前4组点焊试样的S-N曲线,采用3种方法绘制其P-S-N曲线,并确认最优的点焊P-S-N曲线,具体步骤如下.

    (1)通过点焊拉剪疲劳试验得到寿命N,然后对点焊有限元模型进行仿真,提取梁单元的力和力矩并基于Rupp评估法得到点焊的等效结构应力.

    (2)计算试件的破坏率.

    $${p_i} = \frac{i}{{1 + n}}$$ (19)

    (3)将pi作为概率值,在标准正态分布表中查得对应的概率密度函数值ui.

    (4)令yi = lg(Ni),xi = uiP-S-N曲线方程为A + B × up= lg(N). 根据xiyi计算系数LxxLyy,和Lxy,均值$\overline X $$\overline Y $[17].

    $${L_{xx}}{\rm{ = }}\sum {(x - \overline X )^2} = \sum {x_i}^2 - {\left(\sum {x_i}\right)^2}/n$$ (20)
    $${L_{yy}}{\rm{ = }}\sum {(y - \overline Y )^2} = \sum {y_i}^2 - {\left(\sum {y_i}\right)^2}/n$$ (21)
    $${L_{xy}}{\rm{ = }}\sum (x - \overline X )(y - \overline Y ) = \sum {x_i}{y_i} - \sum {x_i}\sum {y_i}/n$$ (22)

    (5)计算P-S-N曲线方程系数AB[17].

    $${{B = }}{L_{xy}}{\rm{/}}{L_{xx}}$$ (23)
    $${{A = }}\overline Y - {\rm{B}}\overline X $$ (24)

    (6)根据需要的存活率ps,查标准正态分布表中对应的概率密度函数值up, 例如:ps = 50%,则up = 0;ps = 99.9%,则up = −3.09.

    (7)将应力Si下存活率psi对应的upi带入P-S-N曲线方程A + Bup = lg(N),得到存活率为psi,应力为Si的寿命的估计值Npi.

    $$ {{A}} + {{B}}\times{u_{pi}} = \lg ({N_{pi}}) $$ (25)

    (8)在对数坐标系中以结构应力幅值Δσs为纵坐标,估算寿命Npi为横坐标,基于最小二乘法拟合得到不同存活率psi下点焊试件的P-S-N曲线如图9所示.

    图  9  点焊试件的P-S-N曲线(方法一)
    Figure  9.  P-S-N curve of spot welded samples (method 1)

    (1)通过试验得到点焊试件的寿命N,然后通过仿真并基于Rupp评估法得到点焊的等效结构应力S.

    (2)求试件在每组载荷下对数寿命的均值$\overline N $,对均值寿命$\overline N $和结构应力S使用最小二乘法拟合,得到均值下的S-N曲线(即50%存活率曲线).

    (3)求试件在每组载荷下对数寿命的标准差$\sigma $,然后求出所有试件对数寿命标准差的均值$\overline \sigma $.

    (4)求不同存活率下的点焊试件的寿命值$\lg {N_p}_i$.

    $$\lg {N_{pi}} = \overline N + \overline \sigma \times {u_{pi}}$$ (26)

    式中:upi为根据需要的存活率psi查标准正态分布表得到对应的概率密度函数值.

    (5)根据结构应力幅值Δσs和不同存活率下的寿命${N_p}_i$,基于最小二乘法拟合得到点焊的P-S-N曲线如图10所示.

    图  10  点焊试件的P-S-N曲线(方法二)
    Figure  10.  P-S-N curve of spot welded samples (method 2)

    (1)通过点焊拉剪疲劳试验得到试件的寿命N,然后通过仿真并基于Rupp评估法得到点焊的等效结构应力S.

    (2)对等效结构应力S和试验寿命N取对数,并采用最小二乘法拟合得到均值下的S-N曲线(50%存活率曲线).

    (3)将所有的点焊试验寿命值N通过S-N曲线等效至5 × 106次循环时的疲劳寿命Ne[18].

    (4)求所有对数等效寿命$\lg {N_e}$的标准差${\sigma _e}$.

    (5)根据需要的存活率psi,在标准正态分布表中可以查得psi值对应的概率密度函数值upi.

    (6)求不同存活率psi对应的点焊估算对数寿命lgNpi.

    $$\lg {N_{pi}} = N + {\sigma _e}\times {u_{pi}}$$ (27)

    (7)在对数坐标系中以等效结构应力幅值Δσs为纵坐标,以估算对数寿命Npi为横坐标,基于最小二乘法拟合得到点焊的P-S-N曲线如图11所示.

    图  11  点焊试件的P-S-N曲线(方法三)
    Figure  11.  P-S-N curve of spot welded samples (method 3)

    通过对比3种方法得到的点焊P-S-N曲线可知,由方法二得到的点焊P-S-N曲线图中一些数据点落在99.9%存活率曲线以外,数据点相对曲线较分散,拟合效果不佳;由方法一和方法三得到的点焊P-S-N曲线图中大部分数据点都在99.9%存活率曲线以内,数据相对曲线较紧凑,拟合效果较好. 综上可得,采用方法一和方法三拟合得到的P-S-N曲线更适合于点焊的疲劳寿命预测,为实际工程应用提供参考.

    (1) 通过拉剪疲劳试验得到点焊试件在不同板厚、焊核直径和载荷比下的疲劳寿命. 根据点焊壳单元有限元模型仿真得到梁单元的力和力矩,并基于Rupp评估法得到点焊试件的等效结构应力.

    (2) 基于点焊试件的疲劳寿命和等效结构应力,采用最小二乘法拟合得到点焊试件在不同应力比以及所有应力比下的S-N曲线.

    (3)对比点焊S-N曲线得出应力比对点焊的S-N曲线影响很大,而且应力比与应力均值成正比,应力均值越大,点焊的疲劳强度越低. 基于Goodman修正法将非对称循环下的结构应力当量为对称循环下的结构应力后得到了点焊的基本S-N曲线,消除了应力比对S-N曲线的影响.

    (4) 基于点焊的S-N曲线,通过三种方法得到了点焊的P-S-N曲线,对比三种P-S-N曲线的拟合效果认为基于破坏率的方法一和基于等效寿命的方法三法拟合效果较好,更适合实际工程应用.

  • 图  1   多孔胞元模型

    Figure  1.   Porous cell model. (a) diamond lattice; (b) spherical six-hole opening structure

    图  2   有限元分析模型

    Figure  2.   Finite element analysis model

    图  3   316L粉末形貌

    Figure  3.   Morphology of 316L powder

    图  4   SLM成形的多孔结构

    Figure  4.   Porous structure formed by SLM

    图  5   不同孔隙率钻石型晶格多孔件表面形貌

    Figure  5.   Surface morphology of diamond lattice porous parts with different porosity. (a) porosity 50%; (b) porosity 60%; (c) porosity 70%; (d) porosity 80%

    图  6   50% ~ 80%孔隙率钻石型晶格结构试验及模拟的应力—应变曲线

    Figure  6.   Experimental and simulated stress-strain curves of diamond lattice structure with 50% ~ 80% porosity. (a) porosity 50%; (b) porosity 60%; (c) porosity 70%; (d) porosity 80%

    图  7   50% ~ 80%孔隙率六孔开口球形多孔结构试验及模拟的应力—应变曲线

    Figure  7.   Experimental and simulated stress-strain curves of spherical porous structures with six-hole openings with 50% ~ 80% porosity. (a) porosity 50%; (b) porosity 60%; (c) porosity 70%; (d) porosity 80%

    图  8   孔隙率对多孔结构力学性能的影响

    Figure  8.   Effect of porosity on mechanical properties of porous structures. (a) change trend of yield strength-porosity; (b) change trend of modulus of elasticity-porosity

    图  9   多孔结构压缩变形过程

    Figure  9.   Compression deformation process of porous structure. (a) diamond lattice with 80% porosity; (b) spherical six-hole opening structure with 80% porosity

    图  10   压缩仿真(左)和试验(右)各阶段变形

    Figure  10.   Deformation at different stages of compression simulation (left) and experiment (right). (a) diamond lattice; (b) spherical six-hole opening structure

    表  1   多孔结构设计参数

    Table  1   Design parameters of porous structure

    孔隙率A(%)钻石型晶格-杆直径d/mm六孔开口球形-壁厚t/mm
    501.27060.7498
    601.09180.5001
    700.91060.3355
    800.71620.2070
    下载: 导出CSV

    表  2   316L粉末的化学成分(质量分数,%)

    Table  2   Chemical composition of 316L powder

    SiCrNiMnMoCSPFe
    0.6416.7911.070.682.530.0270.00560.022余量
    下载: 导出CSV

    表  3   钻石型晶格多孔试样实测参数与设计参数对比

    Table  3   Comparison between measured parameters and design parameters of diamond lattice porous sample

    孔隙率A(%)杆件直径d/mm误差∆d/mm相对误差B(%)
    理论值实际值
    501.2701.3110.0413.24
    601.0911.1660.0756.86
    700.9100.9760.0667.27
    800.7160.7890.07310.16
    下载: 导出CSV
  • [1]

    Wallach J C, Gibson L J. Mechanical behavior of a three-dimensional truss material[J]. International Journal of Solids and Structures, 2001, 38(40): 7181 − 7196.

    [2]

    Wadley H N G. Cellular metals manufacturing[J]. Advanced engineering materials, 2002, 4(10): 726 − 733. doi: 10.1002/1527-2648(20021014)4:10<726::AID-ADEM726>3.0.CO;2-Y

    [3] 常帅. 不锈钢阵列结构选区激光熔化制备与电化学抛光技术研究[D]. 哈尔滨: 哈尔滨工业大学, 2019.

    Chang Shuai. Research on selective laser melting and electrochemical polishing of stainless steel array structures[D]. Harbin: Harbin Institute on Technology, 2019.

    [4] 张钱城, 卢天健, 闻婷. 轻质高强点阵金属材料的制备及其力学性能强化的研究进展[J]. 力学进展, 2010, 40(2): 157 − 169. doi: 10.6052/1000-0992-2010-2-J2008-152

    Zhang Qiancheng, Lu Tianjian, Wen Ting. Processes in the study on enhanced mechanical properties of high-performance lightweight lattice metallic materials[J]. Advances in Mechanics, 2010, 40(2): 157 − 169. doi: 10.6052/1000-0992-2010-2-J2008-152

    [5] 吴林志, 熊健, 马力, 等. 新型复合材料点阵结构的研究进展[J]. 力学进展, 2012, 42(1): 41 − 67. doi: 10.6052/1000-0992-2012-1-lxjzJ2011-095

    Wu Zhilin, Xiong Jian, Ma Li, et al. Processes in the study on novel composite sandwich panels with lattice truss cores[J]. Advances in Mechanics, 2012, 42(1): 41 − 67. doi: 10.6052/1000-0992-2012-1-lxjzJ2011-095

    [6]

    Chantarapanich N, Laohaprapanon A, Wisutmethangoon S, et al. Fabrication of three-dimensional honeycomb structure for aeronautical applications using selective laser melting: a preliminary investigation[J]. Rapid Prototyping Journal, 2014, 20(6): 551 − 558. doi: 10.1108/RPJ-08-2011-0086

    [7]

    Mullen L, Stamp R C, Fox P, et al. Selective laser melting: A unit cell approach for the manufacture of porous, titanium, bone in-growth constructs, suitable for orthopedic applications. II. Randomized structures[J]. Journal of Biomedical Materials Research Part B:Applied Biomaterials, 2010, 92B(1): 178 − 188. doi: 10.1002/jbm.b.31504

    [8]

    Xiao Z, Yang Y, Xiao R, et al. Evaluation of topology-optimized lattice structures manufactured via selective laser melting[J]. Materials & design, 2018, 143: 27 − 37.

    [9] 郑权, 冀宾, 李昊, 等. 基于增材制造的多层金字塔点阵夹芯板抗压缩性能[J]. 航空材料学报, 2018, 38(3): 77 − 82. doi: 10.11868/j.issn.1005-5053.2017.000036

    Zheng Quan, Ji Bin, Li Hao, et al. Compressive behavior of sandwich panels with multilayer pyramidal truss cores by additive manufacturing[J]. Journal of Aeronautical Materials, 2018, 38(3): 77 − 82. doi: 10.11868/j.issn.1005-5053.2017.000036

    [10] 曾寿金, 吴启锐, 叶建华. 选区激光熔化成型316L不锈钢多孔结构的力学性能[J]. 红外与激光工程, 2020(8): 67 − 75.

    Zeng Shoujin, Wu Qirui, Ye Jianhua. Mechanical properties of 316L stainless steel porous structure formed by selective laser melting[J]. Infrared and Laser Engineering, 2020(8): 67 − 75.

    [11]

    Shi C, Lu N, Qin Y, et al. Study on mechanical properties and permeability of elliptical porous scaffold based on the SLM manufactured medical Ti6Al4V[J]. PLOS ONE, 2021, 16(3): e247764.

    [12]

    Hasan R. Progressive collapse of titanium alloy micro-lattice structures manufactured using selective laser melting[D]. UK: University of Liverpool, 2013.

    [13]

    Huo P, Zhao Z, Bai P, et al. Deformation evolution and fracture mechanism of porous TC4 alloy scaffolds fabricated using selective laser melting under uniaxial compression[J]. Journal of Alloys and Compounds, 2021, 861: 158529. doi: 10.1016/j.jallcom.2020.158529

    [14]

    Yang K, Wang J, Jia L, et al. Additive manufacturing of Ti-6Al-4V lattice structures with high structural integrity under large compressive deformation[J]. Journal of Materials Science and Technology, 2019(2): 303 − 308.

    [15]

    Feng Q, Tang Q, Liu Y, et al. Quasi-static analysis of mechanical properties of Ti6Al4V lattice structures manufactured using selective laser melting[J]. The International Journal of Advanced Manufacturing Technology, 2018, 94(5-8): 2301 − 2313. doi: 10.1007/s00170-017-0932-7

    [16] 陈继民, 王颖, 曹玄扬, 等. 选区激光熔融技术制备多孔支架及其单元结构的拓扑优化[J]. 北京工业大学学报, 2017, 43(4): 489 − 495.

    Chen Jimin, Wang Ying, Cao Xuanyang, et al. Topology optimization of microstructure and selective laser melting manufacture for porous scaffolds[J]. Journal of Beijing University of Technology, 2017, 43(4): 489 − 495.

  • 期刊类型引用(14)

    1. 韩彬,高建章,周聪,贾彦杰,刘雪光,谢斐,牛盛源,李立英. 坡度X70管线钢全自动外焊接头组织与性能. 中国石油大学学报(自然科学版). 2024(03): 145-153 . 百度学术
    2. 靳红星. D1219 mm管道山地自动焊系统研制及应用. 压力容器. 2024(06): 36-42 . 百度学术
    3. 高建章,韩彬,贾彦杰,刘雪光,先云泽,秦小省,李立英. 大坡度X70管线钢全自动焊内根焊接头组织与性能. 焊接. 2024(08): 61-68 . 百度学术
    4. 刘铎. 金属管道焊接接头延性断裂行为及显微表征. 信息记录材料. 2023(02): 49-51 . 百度学术
    5. 贾彦杰,周聪,李立英,高建章,杨健,吴易哲,吴彦鹏,韩彬. 管线钢焊接接头韧性影响因素及评价方法. 电焊机. 2023(05): 43-53 . 百度学术
    6. 薛永盛,刘鹏鹏,王峰,宋泽宇. 大管径X80级管线钢焊接弯头热处理工艺研究. 焊接技术. 2023(05): 75-78 . 百度学术
    7. 王鹏宇,闫臣. 从油气管道工程建设的发展看焊接技术的进步. 焊接. 2023(06): 44-51+64 . 百度学术
    8. 丰帆,马明,高彦伟,叶陈峣,李立英,韩彬. 天然气场站用1Cr18Ni9Ti管道全自动横焊接头组织与性能. 焊管. 2023(12): 40-46 . 百度学术
    9. 何小东,高雄雄,David Han,池强,霍春勇,José B.Bacalhau. 高铌X80管道不同强度匹配环焊接头性能及应变演化. 焊接学报. 2022(02): 34-39+115 . 本站查看
    10. 文学,钱建康,汪宏辉,王鹏宇,郭祥,付伟杰,雷正龙. 装配条件对X80钢外根焊接头成分与性能的影响. 焊接. 2022(05): 24-28+35 . 百度学术
    11. 汪宏辉,王鹏宇. 水网地区铜衬垫外根焊全自动焊接技术适应性. 焊接. 2022(07): 60-64 . 百度学术
    12. 武建朝,王晓娥. 自动焊在现代煤化工项目上的研究与应用. 机械制造文摘(焊接分册). 2022(05): 38-43 . 百度学术
    13. 陈自振,程义远,李天伟,任建志,马志军,韩卫亮. X80管线钢药芯焊丝气保护全自动焊焊缝性能研究. 材料导报. 2022(S2): 352-355 . 百度学术
    14. 王汉奎,商学欣,柳旺,宋明,何仁洋. 高钢级油气管线环焊缝应变时效脆化研究. 压力容器. 2021(05): 8-13 . 百度学术

    其他类型引用(1)

图(10)  /  表(3)
计量
  • 文章访问数:  328
  • HTML全文浏览量:  39
  • PDF下载量:  78
  • 被引次数: 15
出版历程
  • 收稿日期:  2021-10-04
  • 网络出版日期:  2022-07-05
  • 刊出日期:  2022-10-30

目录

/

返回文章
返回