基于固体力学的搅拌摩擦焊接数值模拟.docxVIP

基于固体力学的搅拌摩擦焊接数值模拟 实验研究现状 随着新材料、新工艺、新设备的不断涌现，材料过程研究和加工技术以及材料行为之间的联系变得越来越重要。同时，由于金属接合件和其他材料的使用，材料行为的研究提供了新的方向。搅拌摩擦焊接作为一种新型的焊接工艺,已经开始广泛应用于航空航天以及船舶制造等领域,其加工过程如图1所示。搅拌摩擦焊接技术能够较好的解决异种金属焊接这一难题以及不同厚度板类构件连接问题,并可以有效提高焊接接头的性能。Chen等研究了AISI 1018钢和Al 6061之间的搅拌摩擦焊接过程,分析了焊接过程中的温度场的分布情况,并发现在搅拌区内能够形成致密的多相焊接微结构。Li等对2024铝合金和6061铝合金之间的搅拌摩擦焊接工艺进行了研究,在材料流动中观察到有螺旋以及涡旋的产生,并且根据温度的不同,微结构颗粒的大小也有一些差异。Reynolds等将小的Al5454标示物放入到2195焊接构件的焊缝中心线附近的切口内,焊接完成之后对焊接构件进行抛光可以观察到标示物的位置,以研究材料在搅拌摩擦焊接中的流动行为;Guerra等的研究表明,在搅拌摩擦焊接过程中,材料的流动有两种情况,一方面,前进侧的材料将形成一个旋转区,在这个区域内的材料围绕搅拌头旋转的同时伴随在搅拌头轴线方向上的材料的运动,并在旋转数周之后在尾迹当中沉积下来;另一方面,后退侧的材料不会随搅拌头旋转,而仅仅是填充由于前进侧材料的沉积所形成的空间。Li等的研究显示,焊接搅拌区是由一系列相等的半椭圆形区域组成的。Uzun等的工作说明搅拌摩擦焊接可以实现铝合金以及不锈钢之间的连接。 搅拌摩擦焊接作为一种固态焊接工艺,在焊接过程中最高温度也只能达到基体金属焊接熔点的80%～90%。Chao等对搅拌摩擦焊接过程中的温度场做了试验和数值模拟方面的研究,认为在搅拌摩擦焊接工艺中,产生的热仅有大约5%传到搅拌头上,因此,在搅拌摩擦焊接中的热效率可以达到大约95%,这个效率远高于传统熔化焊接工艺当中的60%～80%。Song等用有限元方法建立了搅拌摩擦焊接的传热三维模型,认为预热可以对焊接构件产生有益的影响。汪建华等建立了一个搅拌摩擦焊接传热模型,并尝试使用温度循环的方法获得焊接构件的近似残余应力分布。栾国红等使用Navier-Stocks方程建立了搅拌摩擦焊接流变区域的数值模型,并认为搅拌摩擦焊接构件前进面的塑性流体结构中存在不稳定区域,并在试验中观察到搅拌头周围初始熔化现象的产生。Zhang等采用数值手段对搅拌摩擦焊接过程中材料的流动形式以及材料的一些力学行为进行了分析。张昭等采用数值方法研究了搅拌摩擦焊接构件的残余应力分布。 搅拌摩擦焊接技术可以被用来焊接非直线焊缝,在焊接非直线焊缝的过程中,焊接方向变化的瞬间,搅拌头-焊接构件接触压力的变化以及材料流动行为的变化对深入理解搅拌摩擦焊接工艺的机理具有重要意义。因此,在本文的研究工作中,采用基于固体力学的有限元方法,结合ALE网格技术以及自适应网格重剖分技术,建立了搅拌摩擦焊接过程的三维模型,并着重研究了使用搅拌摩擦焊接技术连接非直线焊缝时材料的流动行为以及一些力学特征,以期可以加深对搅拌摩擦焊接机理的认识,并对工业应用提供有益的帮助。 1 计算模型 1.1 流动率的模拟过程 两块焊接平板的尺寸均为100mm×40mm×8mm,搅拌头的半径R=5mm,肩台的半径为22.7mm。搅拌头在焊接过程中确认为刚体,平板的单元为八节点的六面体单元,采用减缩积分以及沙漏控制以避免单元自锁现象的产生。平板共划分了31105个节点和22785个单元,其几何模型如图1所示,RθZ是一个固定在搅拌头上的柱状坐标系,Z=0mm和Z=8mm分别表示焊接构件的下表面和上表面。在模拟计算中,平移速度被施加到焊接构件的边界上,大小为2.363mm/s,等效于实际加工过程中搅拌头在相反方向的移动,而搅拌头以240rpm的角速度旋转。 焊接构件的材料采用AL6061-T6,作为非率相关的材料处理,材料性质是依赖于温度的,具体如表1所示。 对于理想弹塑性,屈服条件为 Fd(σij?ˉε)=12σijσ′ij-13σ2s(ˉε)(1)Fd(σij?εˉ)=12σijσ′ij?13σ2s(εˉ)(1) 式中σs——屈服应力 ˉεεˉ——等效塑性应变 σij——Kirchhoff应力 σ′ij——Kirchhoff应力的偏量部分 对于非率相关材料的相关联流动,材料的流动方向与屈服面的法线方向ns一致,有 ns=?Fd/?σ∥?Fd/?σ∥=Dp∥Dp∥(2)ns=?Fd/?σ∥?Fd/?σ∥=Dp∥Dp∥(2) 式中Dp——塑性变形率 材料的流动法则为 Dp=˙λ3σ′ˉσ(3)Dp=λ˙3σ′σˉ(3) 式中˙λλ˙——用来表示衡量塑性流动率的一个标量 ˉσ