焊接模拟ansys实例(2020年整理).doc

1、学海无涯焊接模拟ansys实例!下面的命令流进行的是一个简单的二维焊接分析，利用ANSYS单元生死和热-结构耦合分析功能进!行焊接过程仿真，计算焊接过程中的温度分布和应力分布以及冷却后的焊缝残余应 力。fin ish/clear/filnam,1-2D element birth and death/title,Weld Analysis by "Element Birth and Death"/prep7/u nit,si!采用国际单位制!*et,1,13,4!13号二维耦合单元，同时具有温度和位移自由度et,2,13,4!1号材料是钢!2号材料是铝!3号材料是铜!铝是本

2、次分析中的焊料，它将钢结构部分和铜结构部分焊接起来!下面是在几个温度点下，各材料的弹性模量mptemp,1,20,500,1000,1500,2000mpdata,ex,1,1,1.93e11,1.50e11,0.70e11,0.10e11,0.01e11mpdata,ex,2,1,1.02e11,0.50e11,0.08e11,0.001e11,0.0001e11mpdata,ex,3,1,1.17e11,0.90e11,0.30e11,0.05e11,0.005e11!假设各材料都是双线性随动硬化弹塑性本构关系!下面是各材料在各温度点下的屈服应力和屈服后的弹性模量tb,bk in ,1,5

3、tbtemp,20,1tbdata,1,1200e6,0.193e11tbtemp,500,2tbdata,1, 933e6,0.150e11tbtemp,1000,3tbdata,1, 435e6,0.070e11tbtemp,1500,4tbdata,1, 70e6,0.010e11tbtemp,2000,5tbdata,1, 7e6,0.001e11!tb,bk in ,2,5tbtemp,20,1tbdata,1,800e6,0.102e11tbtemp,500,2tbdata,1,400e6,0.050e11tbtemp,1000,3tbdata,1, 70e6,0.008e11tb

4、temp,1500,4tbdata,1, 1e6,0.0001e11tbtemp,2000,5tbdata,1,0.1e6,0.00001e11!tb,bk in ,3,5tbtemp,20,1tbdata,1,900e6,0.117e11tbtemp,500,2tbdata,1,700e6,0.090e11tbtemp,1000,3tbdata,1,230e6,0.030e11tbtemp,1500,4tbdata,1, 40e6,0.005e11tbtemp,2000,5tbdata,1, 4e6,0.0005e11!材料密度（假设为常值）mp,de ns,1,8030mp,de ns,2

5、,4850mp,de ns,3,8900!热膨胀系数（假设为常值）mp,alpx,1,1.78e-5mp,alpx,2,9.36e-6mp,alpx,3,1.66e-5!泊松比（假设为常值）mp ,nu xy,1,0.29mp,nu xy,2,0.30mp,nu xy,3,0.30!热传导系数（假设为常值）mp,kxx,1,16.3mp,kxx,2,7.44mp,kxx,3,393止匕热（假设为常值）mp,c,1,502mp,c,2,544mp,c,3,385!热膨胀系数（假设为常值）!由于该13号单元还有磁自由度，此处假设一磁特性，但并不影响我们所关心的结果mp,murx,1,1mp,mur

6、x,2,1mp,murx,3,1!假设焊料（铝）焊上去后的初始温度是1500 Cmp,reft,1,20mp,reft,2,1500mp,reft,3,20!*!下面建立几何模型csys,0k,1,0,0,0k,2,0.5,0,0k,3,1,0,0!长 1 米k,4,0,0.3,0!厚度0.3米（二维中叫做宽度）k,5,0.35,0.3,0k,6,0.65,0.3,0k,7,1,0.3,0a,1,2,5,4a,2,6,5a,2,3,7,6!戈寸分网格!esize,0.025type,2mat,2amesh,2!esize,0.05网格划分出现问题type,1mat,1amesh,1!mat,3

7、amesh,3eplot!/soluan type,4! 瞬态分析trno pt,full!在模型的左边界加位移约束!nsel,all*get,m inx,no de,m nl oc,xn sel,s,loc,x,mi nxd,all,ux,0*get,m iny,no de,m nl oc,yn sel,r,loc,y,mi nyd,all,uy,0!*假设模型的左右边界处温度始终保持在20摄氏度左右*!其他边界条件如对流和辐射等均可施加，此处因为只是示意而已，故只施加恒温边 界条件nsel,all*get,m inx,no de,m nl oc,xn sel,s,loc,x,mi nxd,

8、all,temp,20nsel,all*get,maxx ,no de,mxloc,xn sel,s,loc,x,maxxd,all,temp,20 首先将该区域的单元杀死” !nn a=2esel,all*get,emax,elem, num ,maxasel,s,area, nnaesla*get ,n se,elem,co unt*di m,n e, nse*di m,n ex, nse*di m,n ey, nse*di m,n eorder, nsemin e=0!*!下面的do循环用于将焊料区的单元按其形心y坐标排序!以便后面模拟焊料由下向上逐步生长”过程*do,i1,1, nse

9、esel,u,elem,mi ne*get ,n se1,elem,co untii=0*do,i,1,emax*if,esel(i),eq,1,the nii=ii+1n e(ii)=i*en dif*enddo*do,i,1, nse1*get, ney(i),elem ,n e(i),ce nt,y*get, nex(i),elem ,n e(i),ce nt,x*enddomin y=1e20min x=1e20*do,i,1, nse1*if, ney(i),lt,m in y,the nminy二n ey(i)minx二n ex(i)mine二n e(i)*else*if, ney

10、(i),eq, min y,the n*if, nex(i),lt,m in x,the nminy二n ey(i)minx二n ex(i)mine二n e(i)*en dif*en dif*en dif*enddon eorder(i1)=mine*enddo1500 C!* max_tem=1500 !按照前面假设，焊料的初始温度为 dt仁1e-3!用于建立初始条件的一个很小的时间段dt=5!焊接一个单元所需的时间t=0!起始时间esel,alleplot/auto,1/replot*do,i,1, nseekill, neorder(i)esel,s,liveeplot*enddoall

11、sel,alloutres,all,allic,all,temp,20kbc,1timi nt,0,structtimi nt,1,thermtimi nt,0,magtin tp,0.005,1,0.5,0.2!nsub 仁2nsub2=40!*do,i,1, nseealive ,n eorder(i)esel,s,liveeplotesel,all!*下面的求解用于建立温度的初始条件*t=t+dt1time,tn subst,1*do,j,1,4d,n elem( neorder(i),j),temp,max_tem*enddosolve!*下面的求解用于保证初始的升温速度为零*t=t+

12、dt1time,tsolve!*下面的步骤用于求解温度分布*do,j,1,4ddele ,n elem( neorder(i),j),temp*enddot=t+dt-2*dt1time,tn subst, nsublsolve*enddot二t+50000!*下面的步骤用于冷却过程求解 *time,tn subst, nsub2solvesavefin ish!后处理过程！!！/post1!*下面的一系列命令用于生成应力的动画文件*/seg,dele/con t,1,15,0,1200e6/16,1200e6/dscale,1,1.0avpri n,0,0avres,1/seg,multi,

13、stress1,0.1esel,all*do,i,1, nseesel,u,elem, neorder(i)*enddo*do,i,1, nseesel,a,elem, neorder(i)set,(i-1)*3+1,1pin sol,s,eqv*do,j,1, nsublset,(i-1)*3+3,jpin sol,s,eqv*enddo*enddo*do,i,1, nsub2set,( nse-1)*3+4,ipin sol,s,eqv*enddo/seg,off,stress1,0.1/an file,save,stress1,avi!*下面的一系列命令用于生成温度的动画文件*/seg,dele/con t,1,15,0,1500/16,1500/dscale,1,1.0avpri n,0,0avres,1/seg,multi,temp1,0.1esel,all*do,i,1, nseesel,u,elem, neorder(i)*enddo*