第30例单元生死应用实例—焊接模拟.doc

第30例 单元生死应用实例—焊接模拟 /CLEAR !清除数据库，新建文件 /FILNAME, EXAMPLE30 !改变工作名 !开始热分析 !建立有限元模型 /PREP7 !进入前处理器 ET,1,PLANE77 !平面热单元 ET,2,SOLID90 !空间热单元 V=0.005 !焊接线速度 LENGTH=0.05 !钢板沿焊缝方向长度 L_NUMBER=LENGTH/V !沿焊缝方向单元段数 MPTEMP,1,0,300,600,1000,1300,2000 !定义材料模型，温度表 MPDATA,DENS,1,1,7859,7770,7659,7600,6750,6750 !密度 MPDATA,DENS,2,1,7859,7770,7659,7600,6750,6750 MPDATA,C,1,1,450,514,967,1370,2800,2000 !比热 MPDATA,C,2,1,450,514,967,1370,2800,2000 MPDATA,KXX,1,1,67,53,39,31,20,20 !传导系数 MPDATA,KXX,2,1,67,53,39,31,20,20 RECT,0,0.10,0,0.01 !因为对称，创建1/2实体模型 K,5,0.002 K,6,0.002,0.002 K,7,0.00662,0.01 K,8,0.00339,0.01112 K,9,0,0.0115 LARC,7,9,8 A,1,5,6,7,9 AGEN,2,2 ASBA,1,2 AGLUE, ALL K,20,0,-0.00385 CIRCLE,20,0.008,,,90 CIRCLE,20,0.012,,,90 ASBL,3,4 ASBL,2,5 ESIZE, 0.00075 SMRTSIZE,7 MSHAPE,1 AMESH,3,5,1 ESIZE, 0.003 AMESH,1 !划分单元 EXTOPT,ESIZE,L_NUMBER EXTOPT,ACLEAR,1 MAT,1 VEXT,3,5,1,,,LENGTH !挤出 MAT,2 VEXT,1,,,,,LENGTH NUMMRG,NODE !合并两次挤出产生的位置重合的节点 ALLS FINISH !热分析求解 !作稳态热分析施加初始温度 /SOLU !进入求解器 AUTOTS,ON !打开自动时间步长 OUTPR,NSOL,ALL !输出控制 OUTRES,NSOL,ALL KBC,1 !阶跃载荷 ANTYPE,TRANS !瞬态分析 TIMINT,OFF !关闭时间积分 NSUBST,4 !子步数目 D,ALL,TEMP,20 !施加初始温度 TIME,0.01 !时间 SOLVE !解答 !杀死焊缝单元 ESEL,S,MAT,,1 !选择焊缝单元 EKILL,ALL !杀死单元 ESEL,ALL TIME,0.02 SOLVE !开始激活焊缝单元，进行瞬态热分析 TIMINT,ON,THER DDELE,ALL,TEMP !删除前面用D命令施加的温度载荷 TIME=0.02 !变量TIME的初始值 *DIM,QQ,ARRAY,8 !定义数组，存储单元节点热流密度 *DIM,E_NODE,ARRAY,20 !定义数组，存储单元所属节点编号 *GET,ELEM_MAX, ELEM, 0, NUM,MAX !获取最大单元号 *GET,ELEM_MIN, ELEM, 0,NUM,MIN !获取最小单元号 *GET,NODE_MAX, NODE, 0, NUM,MAX !获取最大节点号 *GET,NODE_MIN, NODE, 0,NUM,MIN !获取最小节点号 *DIM,NODE_AVE,ARRAY,NODE_MAX !定义数组，存储节点相邻单元数量 *DO,J,NODE_MIN,NODE_MAX !计算节点相邻单元数量，循环开始 ALLS *IF,NSEL(J),EQ,1,THEN NSEL,S,,,J !选择节点J ESLN,S !选择节点J相邻单元 *GET,NODE_AVE(J), ELEM, 0,COUNT !获取单元数量 *ENDIF *ENDDO !计算节点相邻单元数量，循环结束 ALLS *DO,KK,1,3(我用的是两层焊缝) !计算三层焊缝温度场，循环开始 R=0.02 !电弧有效加热半径 *IF,KK,EQ,1,THEN !第一层焊缝的基本参数 Q=3200*0.7 !输入热量 QMAX=3*Q/3.1415926/R/R !变量 XC= 0.11190E-02 !第一层焊缝质心X坐标 YC= 0.21451