资源描述：
《第30例单元生死应用实例焊接模拟》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、第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,
2、1,7859,7770,7659,7600,6750,6750MPDATA,C,1,1,450,514,967,1370,2800,2000!比热MPDATA,C,2,1,450,514,967,1370,2800,2000MPDATA,KXX,1,1,67,53,39,31,20,20!传导系数MPDATA,KXX,2,1,67,53,39,31,20,20RECT,0,0.10,0,0.01!因为对称，创建1/2实体模型K,5,0.002K,6,0.002,0.002K,7,0.00662,0.01K,8,0.00339,0.01112K,9,0,0.0115LAR
3、C,7,9,8A,1,5,6,7,9AGEN,2,2ASBA,1,2AGLUE,ALLK,20,0,-0.00385CIRCLE,20,0.008,,,90CIRCLE,20,0.012,,,90ASBL,3,4ASBL,2,5ESIZE,0.00075SMRTSIZE,7MSHAPE,1AMESH,3,5,1ESIZE,0.003AMESH,1!划分单元EXTOPT,ESIZE,L_NUMBEREXTOPT,ACLEAR,1MAT,1VEXT,3,5,1,,,LENGTH!挤出MAT,2VEXT,1,,,,,LENGTHNUMMRG,NODE!合并两次挤出产生的位置重
4、合的节点ALLSFINISH!热分析求解!作稳态热分析施加初始温度/SOLU!进入求解器AUTOTS,ON!打开自动时间步长OUTPR,NSOL,ALL!输出控制OUTRES,NSOL,ALLKBC,1!阶跃载荷ANTYPE,TRANS!瞬态分析TIMINT,OFF!关闭时间积分NSUBST,4!子步数目D,ALL,TEMP,20!施加初始温度TIME,0.01!时间SOLVE!解答!杀死焊缝单元ESEL,S,MAT,,1!选择焊缝单元EKILL,ALL!杀死单元ESEL,ALLTIME,0.02SOLVE!开始激活焊缝单元，进行瞬态热分析TIMINT,ON,THER!
5、打开时间积分，进行瞬态分析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,
6、ARRAY,NODE_MAX!定义数组，存储节点相邻单元数量*DO,J,NODE_MIN,NODE_MAX!计算节点相邻单元数量，循环开始ALLS*IF,NSEL(J),EQ,1,THENNSEL,S,,,J!选择节点JESLN,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!输入热量QM
7、AX=3*Q/3.1415926/R/R!变量XC=0.11190E-02!第一层焊缝质心X坐标YC=0.21451E-02!第一层焊缝质心Y坐标V_NUM=2（这些根据我自己的模型修改的）!第一层焊缝的体号TIME_INC=LENGTH/L_NUMBER/V!第一层焊缝时间增量*ENDIF*IF,KK,EQ,2,THEN!第二层焊缝的基本参数Q=4500*0.7QMAX=3*Q/3.1415926/R/RXC=0.20373E-02YC=0.60454E-02（根据自己模型）V_NUM=1TIME_INC=LENGTH/L_NUMBER
