焊接应力与变形数值模拟领域的若干关键问题.docVIP

焊接应力与变形数值的关键问题 摘要为了进一步提升焊接应力与变形数值模拟在工业技术领域的应用，本文提出了目前该领域所面临的主要问题，并针对材料性能参数的缺失、边界条件施加的随意性、大型焊接结构数值模拟效率低下、连续加工过程实现困难等几方面进行详细的阐述。提出了解决该类问题的一下思想方法。 关键词 焊接应力变形数值模拟 0前言 数值模拟技术已经广泛应用于焊接科学与技术的各个领域，依据拉达伊的分类方式[12]，大致将焊接领域的数值模拟分为三类：焊接过程模拟、材料模拟以及结构的模拟。图1中则表明了与这三个方面相关的热力过程。相应的今年来很多学者对于这三方面进行了研究。目前焊接数值模拟技术已经应用于几乎所有的焊接方法，如电弧焊[3]、激光焊[4]、扩散连接[5]、钎焊[6]、搅拌摩擦焊[7]等；同时数值模拟技术也在对某些焊接过程中的一些现象的机理进行了研究，如冷裂纹形成机理、应力腐蚀开裂机理[89]、熔池内流体的流动[10]及断裂力学[11]等；在焊接残余应力与变形领域[12-14]，该技术也得到了广泛的应用。 图1 与焊接数值模拟相关的物理领域关系原理图(修改自拉达伊) 其中与工程联系紧密的焊接应力与变形问题一直是焊接模拟领域的研究热点，经过三十余年的发展，目前数值模拟技术已经实现了焊接应力与变形从简单到复杂、从二维到三维、从定性到定量的分析。 虽然在焊接应力与变形数值模拟过程中，研究者考虑的焊接因素越来越多，结果也越来越精确，但同时实际焊接结构的数值模拟也面临很大的挑战，本文将就这些困难与问题进行论述，并提出一些相应的解决思路。期待该综述能够为未来的焊接结构焊接应力与变形的数值模拟提供有益的思路。 1 材料高温性能的确定原则 1.1 材料高温性能的缺失 材料的性能参数是数值模拟过程中必不可少的物理量之一，其选取的正确与否，对模拟结果的准确性起决定的作用。从目前的文献来看，对材料的物理性能研究主要集中于相对较低的温度区间，除个别常用材料外，对其他材料的高温物理性能（接近熔点及熔化后）参数研究还比较少，这也是目前数值模拟工作的难点之一。 目前对高温下的材料性能参数值一般是通过外延的方法得到的。对于弹性模量及屈服强度，在材料达到熔点附近及以上温度时，都被设置为一个很小的值；热导率、比热及线膨胀系数在各个温度下变化不是很大，高温参数值与低温部分基本保持不变。 这里可以看出很多参数的设置还带有一定的随意性，一般都是以计算能够收敛作为基本的要求。如弹性模量及屈服强度通常取室温下参数的千分之一。另外诸如热导率在熔化前后实际上差别似乎非常大的，这主要是由于熔化前后改变了传热的基本模式，在固态下材料主要的传热模式为热传导，而熔化后，液态金属的对流将对于传热起到非常大的作用。这也是面状热源仅适合于熔深较浅、电弧吹力比较小的场合的原因。 另外一个影响数值模拟准确性的因素是模拟过程中材料高温参数通常仅仅设置为温度的函数，而对于应变大小及应变速率对于材料性能的影响则较少涉及，而实际上高温性能参数尤其是力学性能应该为温度、应变量及应变速率的函数。然而目前的通用有限元软件很难综合考虑以上的因素，这一点将是焊接结构与力学数值模拟领域的研究的难点。 目前的问题不仅仅是对于以上的综合考虑温度、应变量及应变速率的材料本构模型尚不完善，即使对于单个因素的影响的研究也存在一定的缺陷，如大部分的研究者所采用的材料性能参数随温度变化的数据是通过高温拉伸试验机，获得的相应恒定温度下的结果，事实上，一定的材料的在一定的温度下的性能会随着这个材料的加热或者冷却历程而存在差异，也就是说这对特定温度焊接加热过程中材料参数与冷却过程中的参数会存在一定的差异，这种差异会直接影响计算结果的准确性。图2为10Ni5CrMoV钢屈服强度与温度的关系[15]，从中可以看出在室温与700℃的温度区间内, 冷却条件下测定的各温度的屈服强度要远小于稳态条件下的测定值。 图2 10Ni5CrMoV钢屈服强度与温度的关系[15]1.2 材料参数的处理方法 针对材料参数在熔点附近会出现变化的问题，其可能的解决方法主要包括以下两方面，其一是熔化后采用液态下的相关参数，当然如前所述，这部分的性能比较难于测定。尤其是力学性能，当材料变成液态后，屈服及弹性模量的具体定义都会出现问题。所以如果采用这种方式还是没办法脱离传统的处理方法局限性（即收敛判据），也只能对于热性能参数进行适当的优化。另外一种方法是综合采用弹塑性力学与流体力学，对于处于固态的热影响区及母材采用弹塑性本构方程，而熔池采用流体力学的计算模型。目前流体力学在焊接领域的应用主要在熔池流体流动行为研究以及搅拌摩擦焊等领域的应用，将该方法应用于熔化焊接还需要进行深入的研究。 采用与加热及冷却过程直接相关的随温度变化的力学参数进行数值模拟，有望大幅度提