基于ANSYS的中厚板焊接有限元三维数值模拟.docVIP

收稿日期:2002206203. 作者简介:罗金华(19622 , 男, 博士研究生; 武汉, 华中科技大学材料科学与工程学院(430074 . 现工作单位:武昌 造船厂(430070 . 基金项目:国防科技预研基金资助项目. 基于ANS YS 的中厚板焊接有限元三维数值模拟 罗金华　王晓熙　胡伦骥 (华中科技大学材料科学与工程学院 摘要:基于大型有限元分析软件ANSYS , 对中厚板焊接的温度场和应力应变场进行三维数值动态模拟, 并将计算量控制在可接受的范围内. 建模时采用两种单元结合以获得焊缝处细密、, 热载荷施加过程中采用余量控制法, . 却过程中角变形沿纵向并不总是线性分布的. 关　键　词:焊接; 有限元分析; 数值模拟; ANSYS ; 中图分类号:TG 44　　文献标识码:A 　　 204 　　, (, 考, 对中厚板焊接的温度场和应力应变场进行三维数值动态模拟, 并控制计算量在可接受的范围内, 以期能将研究成果运用到更为复杂构件焊接过程的数值模拟中. 1　建模 研究10. 7mm ×200mm ×200mm 的D32钢板. 焊接试验采用SH5073. 2mm 焊条沿正面中心线进行焊接. 焊接参数为:焊接电流I =120A , 焊接电压U =25～26V , 焊接速度V =3mm/s. 试板正面经加工刨平, 反面沿焊缝线钻一个4. 5mm 不同深度的平底孔并焊上铂2铑热电偶用以 测量温度. 焊前、焊后通过测量事先布置好的测量杆顶部的跨距变化来测定焊接产生的角变形. 焊接过程是个局部快速加热到高温并随后冷却的过程, 随着热源的移动, 整个焊件的温度随时间和空间急剧变化, 材料热物理性能也随温度剧烈变化, 同时还存在熔化和相变时的潜热现象[1]. 由于高度非线性及计算量等因素, 因此在对平板堆焊的温度场和应力应变场分析中作以下假设:a . 不考虑工件与实验台之间的热传导, 假设工件的所有外边界仅与空气发生对流换热, 将辐射换热的影响折算到对流换热中, 而不单独考 虑; b . 忽略熔池内部的化学反应和搅拌、对流等现象; c . 忽略焊条与母材材料的不一致性, 采用随温度变化的热物性参数; d . 工件的初始温度为 室温(20℃ ; e . 焊接以恒定速度V 进行, 电弧的能量密度服从高斯分布. 有限元模型是真实系统的数值模型. 确定单元类型、材料特性、几何模型以及网格划分都是温度场分析和应力应变场分析中的关键环节. 本研究采用文献[2]中的材料热物理性能参数, 分析中应用间接耦合场法, 下面只讨论建模中单元类型和网格划分方案的确定. 1. 1　定义单元类型 焊接过程是个热力耦合的过程, 在热弹塑性分析的过程中, 既需要计算温度场, 又需要计算应力应变场, 选择的单元必须是自由度为温度的热 单元, 具有热传导、对流能力(有的单元仅有辐射、对流功能 , 而且必须能够进行热力耦合分析. 在本课题中, 研究的是中厚板焊接过程的应力应变场, 必须考察板内厚度方向的温度和应力应变分布, 因此采用三维单元. 此外, 考察的是整个焊接过程的温度场和应力应变场的发展过程, 因此采用能够进行瞬态非线性分析的单元类型. 在此研究的是平板堆焊, 由以上几条原则, 最终可把选用单元的范围缩小为ANSYS 单元类型库中的8节点三维热实体单元solid70和20节点的三维热实体单元solid90[3]. 为了保证焊接熔池 第30卷第11期　　　华　中　科　技　大　学　学　报(自然科学版 　　　　Vol. 30　No. 112002年　11月　　J. Huazhong Univ. of Sci. Tech. (Nature Science Edition 　　　Nov. 　2002 处的计算精度并控制适当的计算工作量, 模型的不同区域采用不同的单元, 在焊缝和远离焊缝的区域用solid70, 中间部分用solid90单元过渡. 整个模型单元数为6355, 节点数为8104, 不会出现网格过渡错误, 计算量也可以接受. 另一项节省计算时间的措施是转变退化四面体单元. 用solid90的目的是为了生成过渡的金字塔单元, 在生成过渡的金字塔后, 可将模型中的20节点退化四面体单元solid90转化为相应的10节点非退化单元solid87, 这可使每个单元所需的随机存储单元(RAM 更少[4]. 1. 2　网格划分方案自由网格对于单元形状无限制, 但排列不规则. 而映射网格对包含的单元形状有限制, 须满足特定的规则. 角形单元, 且, , 是排的单元[5] , 的, . 出于提高计算效率的考虑, 在远离焊缝处应采用较焊缝处粗略的网格. 若全部采用六面体网格, 则由于六面体网格的局限性, 焊缝处网格的细化只能发生在平板的宽度和厚度方向, 而长度方向(焊接方