基于网格映射的非线性损伤累积研究
材料的应变可分为两个部分：其一为弹性应变，其二为塑性应变。在物体的损伤累积过程中，弹性应变存在恢复，而塑性应变将传递到物体下次受载工况中，所以单纯采用有限元的显式算法来预测这种累积效应，只能实现瞬态场量（如瞬态变形几何和应变等）的累积，而不能实现稳态场量（如稳定几何和塑性应变等）的损伤累积，其结果往往不可靠。为了解决瞬态场量不是稳定量这类问题，通常的处理方法是选择物体运动状态转变收敛的情况时，即各能量转换较稳定的状态值，近似认为其处于稳定状态，并将该瞬态时刻作为下次受载工况的起点，从而实现损伤的累积。同时，可以选择只传递塑性应变而不包含弹性应变的方法来避免非真实量的传递。但即使这样，也只能是其瞬态几何累积，而不是物体的稳态变形几何的累积，故存在一定误差。而使用隐式非线性分析理论，能获得物体稳态状态场变量值，但计算资源消耗大，且对于高度非线性问题极难收敛，因此非线性损伤累积在实际中运用较少。
本文通过使用LS-DYNA软件在显式非线性模块进行瞬态计算，获得物体受载条件下的各场量响应，并将这些瞬态场量值（包含应力、应变、厚度），通过PIM[10]（Point Interpolation Method）网格映射，传递到物体的瞬态变形几何，同时运用LS-DYNA软件隐式非线性模块进行回弹分析，从而获得物体在该受载情况下的稳定变形几何和稳态场变量值。然后，通过同样网格映射技术，将稳定场变量值和稳定几何作为第二次瞬态受载工况的初始条件，从而实现其非线性损伤的传递。通过前次相同的方法，可获得物体第二瞬态工况后的稳定场变量值和稳定几何，从而实现其真实损伤传递。依次类推，即可获得物体N次非线性损伤累积。采用这样的处理方法可避免裂纹扩展理论难以计算大物理模型长时间的损伤累积过程和显式有限元计算无法获得稳态场量值的缺陷来进行物体损伤累积过程。本文通过网格映射实现场量传递，通过变形网格的传递考虑物体受载几何变形对损伤的影响，并综合显隐式的计算优势，从而实现非线性损伤累积过程仿真。
本文以某空调外机两次跌落分析为例，获得空调外机连续两次跌落的损伤累积情况。并与试验相互验证，证明该方法的正确性，并最终解决某空调外机跌落变形过大问题。
1 非线性损伤累积方法
1.1 非线性损伤累积分析流程
本文中的物体非线性损伤计算主要包含两个部分：一为部件加载后的部件的瞬态应变、应力和厚度响应;二为部件的稳定几何变形的计算。本方法实现两个部分的损伤累积，其分析流程（以连续2次损伤累积为例，N次损伤累积可以依次类推）如图1所示。
仿真中部件加载后几何形状的传递通过变形的网格来传递。而每次计算后，部件应力、应变、厚度值等场量传递通过网格映射的方法实现。同样在N次瞬态分析后，将部件加载后应力、应变、厚度值映射到物体受载瞬态变形网格中，并进行回弹分析，而获得部件稳态变形及稳态场量值。通过网格映射和变形网格传递，将部件N次受载响应的稳态场量和稳态变形几何，作为第N+1瞬态受载工况的初始条件，从而实现物体的损伤传递和累积。
2 运用实例
本文以某公司某空调外机跌落分析为例。根据该企业运输法规要求，空调外机棱边与地面呈45°，在高度600 mm处连续两次自由跌落，外机零件不能出现明显变形和破坏且外机所有功能正常（如图3所示）。
由于跌落试验是不连续的冲击，并且第二次跌落在第一次的跌落上进行，所以第二次的零件损伤是在第一次跌落后的基础上累积。此外，考虑第一次跌落后，零件存在较大的回弹，特别是空调底脚局部为变形最大区域，同时是承受跌落冲击主要区域。第一次跌落后零件的几何变形将影响跌落冲击载荷的分布，进而影响损伤分布（如图4所示）。因此，必须考虑底脚在第一次跌落后回弹后的稳定几何变形对于第二次跌落损伤影响。在本案例中，采取上文所介绍的基于网格映射的损失累积方法，考虑其影响。
2.1 空调外机有限元模型的建立
根据某空调外机几何模型，将空调外机各钣金件划分成尺寸为10 mm网格，并根据几何特征局部细化，求解使用LS-DYNA软件。模型中包含空调外机各重要质量附件如电机、压缩机等外形并赋予实测质量。有限元模型如图5所示，包括175 487个节点、158 161个单元。空调外机有限元模型总质量为29.66 kg，试验实测质量为30.4 kg。
本次分析中钣金材料和泡沫材料均为实测值，本例中由于初速度为零，冲击加速度仅为1G，故不考虑材料应变率影响，其中主要钣金材料拉伸试验获得的材料曲线如图6所示。
2.2 有限元计算
本次计算案例中，通过显式有限元获得物体跌落变形结果。基于计算时间和精度的考虑，其中沙漏能量控制在4.9%。物体的回弹分析（获得物体稳定状态）开始于瞬态跌落中内能和动能转化较稳定时间点，约30 ms处。其中，第一次跌落的主要能量输出如图7所示。
数值计算中，物体在一次加载工况后的瞬态几何变形，与实际中物体加载后的稳定几何变形存在较大不同。图8为第一次跌落分析后，空调外机底脚的瞬态计算结果和进行回弹处理而获得的稳态计算结果，其應力云图分布、底脚的变形几何均存在较大的差异，前者计算结果明显偏大。直接将瞬态计算结果传递到下次加载工况中将导致过大的误差，因此必须考虑循环载荷作用下，物体在卸载后的几何变形对物体相关局部区域的损伤及其累积的影响（如图9所示）。
回弹约束点的设置对于结果影响较大。对于本文跌落案例中，对于主要关注的部件空调外机底脚，选取其焊点位置为回弹约束点。该处由于为相当较软底脚与相对较硬底盘焊接连接处，约束刚度较大，故为回弹约束点。
3 计算结果与优化
3.1 计算结果与实验
在某空调外机连续两次跌落试验中，外机底脚出现明显变形，实测最大位移偏移约11.2 mm，无法通过企业自身运输法规。底脚的材料为DX52D，厚度为1.2 mm。同时，该企业为降低产品成本，要求将该零件厚度降低为1.0 mm，并通过企业运输法规。在本案例中，通过前文所介绍的基于网格映射的损伤累积方法，计算该空调外机的原始方案，在连续两次跌落后，底脚最大变形量约11.8 mm。不能通过该法规。
根据试验和仿真结果可知：底脚的出现明显变形主要是因为图10所示处，为底脚台阶边缘。针对该处薄弱环节，在底脚两次各增加5 mm长的翻边，同时将底脚的厚度减薄为1.0 mm（如图11所示）。
计算原方案及优化方案的连续两次跌落的最终变形。同时，为研究需要，设计对比算例。该对比算例在第一次跌落分析稳定后，将其塑性应变、厚度及瞬态变形几何直接作为第二次跌落分析的初始条件。计算结果见表1。
根据计算结果可知，优化方案底脚处的最终变形为4.3 mm，底脚减薄厚度并满足企业运输法规的要求，企业最终采用该方案。同时根据对比案例可知：在计算损伤累积时，如果只传递物体的瞬态场量和瞬态变形几何，其并未考虑物体的回弹效应，将导致结果误差为28.6%。而使用本文介绍的基于网格映射的损伤累积方法，将物体的稳态场量和稳态变形几何传递到下次加载工况，误差约5.4%。
4 结语
本文提出在不连续载荷下物体的损伤累积的计算方法，在该过程中，损伤的传递必须是在前一次瞬态工况后物体的稳态场量，其包含物体的应力、应变值，厚度变化及稳定变形几何。此外，提出基于网格映射的损伤累积方法，其利用显式有限元在非线性问题上收敛性和计算效率优势下，来求解物体的加载后瞬态响应。并通PIM网格映射方法，将物体受载后瞬态场量映射到其变形网格上，并进行隐式回弹分析，从而获得物体的稳态场量和稳态变形几何，并将稳定状态值作为下次瞬态载荷工况的初始条件，从而实现物体的损伤传递。这样同时利用显隐式有限元优势，较快捷地完成物体的损伤传递及累积。该方法最终在某机型空调外机的连续两次试验中获得验证，成功解决了不连续冲击下的支架损伤的累积，并提出了优化方案，大大缩短了空调开发的周期，具有重要的意义。
参 考 文 献
[1]Miner M A.Cumulative Damage in Fatigue[J].ournal of Applied Mechanics 1945（12）：159-164.
[2]Kachanov LM.Time of the rupture process under creep condition[J].Izvestiia Akademii Nauk SSSR，1958（8）：26-31.
[3]Lemaitre J，Chaboche J L.A nonliner model of creep-fatigue damage cumulation and Interaction[M].Proc，IUTAM，Symp of Mech-anics of Viscoelastic Media and Bodies. Getheubury：Springer，1974.
[4]Paris PC.The growth of cracks due to variations in loads[C].PhD thesis，Lehigh University，1960.
[5]Suresh S.Fatigue of materials[M].Cambridge：Ca-mbridge University Press，1998.
[6]Kanninen MF，Poplar CH.Advanced fracture mechanics[M].London：Oxford University Press，1985.
[7]Wheeler OE.Spectrum loading and crack growth[J].J Basic Eng 1972，94（1）：181-186.
[8]Willenborg J，Engle RM，Wood H.A crack growth retardation model using an effective stress intensity concept[M].Technical Report TFR 71-701，North American Rockwell，1971.
[9]Christensen RM，Miyano Y.Deterministic and pro-babilistic lifetimes from kinetic crack growth-generalized forms[J].International Journal of Fracture，2007（137）：77-87.
[10]G.R. Liu.，Meshfree methods：Moving Beyond the Finite Element Method[M].2nd Edition，CRC Press：Boca Raton USA，2009.
[責任编辑：钟声贤]