镁合金焊缝数值模拟.docVIP

镁合金的残余应力及温度场分析（直接法） 一、问题的描述 如图1所示为焊缝的剖面示意图，为镁合金的双面焊接，求在焊接后的焊接残余应力以及在焊接过程中的温度变化。假定没有对流和辐射，结构左面固定约束，左右两边界均给定20℃的温度约束。整个材料都是镁合金，其中焊缝左右的初始温度是环境温度20℃，焊缝是1500℃。 镁合金材料的属性如表1所示，其中材料1和3，分别是焊缝的左边和右边金属的属性，材料2是焊缝金属的属性。 材料1、2 摄氏 温度℃ 弹性模量 Pa 屈服强度 Pa 切变模量 Pa 材料密度 Kg/m3 泊松比 传热系数 W/(m·K) 线胀系数 1/℃ 比热容 J/(Kg℃) 镁 合 金 20 1.02e11 0.8e9 1.02e10 1740 0.35 51 26.1e-6 102 500 0.5e11 0.4e9 0.5e10 1000 0.08e11 0.07e9 0.08e10 1500 0.001e11 0.001e9 0.001e10 2000 0.0001e11 0.0001e9 0.0001e10 材料3 摄氏 温度℃ 弹性模量 Pa 屈服强度 Pa 切变模量 Pa 材料密度 Kg/m3 泊松比 传热系数 W/(m·K) 线胀系数 1/℃ 比热容 J/(Kg℃) 镁 合 金 20 2.02e11 0.8e9 1.02e10 1740 0.35 51 26.1e-6 102 500 0.6e11 0.4e9 0.5e10 1000 0.09e11 0.07e9 0.08e10 1500 0.002e11 0.001e9 0.001e10 2000 0.0002e11 0.0001e9 0.0001e10 表1 二 前处理 （1）定义工作文件名：Utility MenuFileChange Jobname,在弹出的对话框中输入文件名为Weld，单击OK按钮。 （2）定义工作标题：Utility MenuFileChange Title，在弹出的对话框中输入weld，单击OK按钮。 （3）关闭坐标符号的显示：Utility MenuPlotCtrlsWindow ControlsWindow Options,弹出一个对话框，在Location of triad后面的下拉式选择栏中选择No Shown，单击OK按钮。 （4）定义单元类型：Main MenuPreprocessorElement TypeAdd/Edit/Delete,出现Element Type对话框，单击Add按钮，出现Library of Element Type对话框，如图2所示，在该对话框左边的列表中单击Coupled Field，在其右面的列表栏中单击Vector Qaud13，单击OK按钮。 图2 （5）设置单元选项：单击Element Types对话框的Options按钮，弹出PLANE13 element type options对话框，如图3所示，在K1后面下拉列表中选择UX UY TEMP AZ，在K3后面的下拉列表中选择Plane strain，单击OK按钮。最后单击Element Types对话框的Close按钮，关闭该对话框。 图3 (6)定义材料的弹性模量：Main MenuPreprocessorMaterial PropsMaterial Models,出现Difine Meterial Model Behavior对话框，如图4所示，在Material Model Available栏中，连续双击StructuralLinearElasticIsotropic,出现Linear Isotropic Material Properties For Material Number 1对话框，连四次单击Add Temperature，将生成如图5所示的表格图，在Temperature后面分别填入200，500，1000，1500，2000，在EX后面分别填入镁合金在各温度时刻对应的弹性模量，在PRXY后面填入对应的泊松比，单击Graph按钮可绘出弹性模量的曲线图，如图6，然后单击OK。 图4 图5 图6 （7）设定材料的线膨胀系数：连续双击Difine Meterial Model Behavior对话框中的StructuralThermal Expansion Secant CofficientIsotropic，出现Thermal Expansion Secant Cofficient For Material Number 1对话框，在Reference Temperature后面输