焊接模擬ansys實例

1、焊接模擬ansys實例!下面的命令流進行的是一個簡單的二維焊接分析, 利用ANSYS單元生死和熱-結(jié)構(gòu)耦合分析功能進 !行焊接過程仿真, 計算焊接過程中的溫度分布和應(yīng)力分布以及冷卻后的焊縫殘余應(yīng)力。finish/clear/filnam,1-2D element birth and death /title,Weld Analysis by Element Birth and Death/prep7/unit,si !采用國際單位制!*et,1,13,4!13號二維耦合單元, 同時具有溫度和位移自由度et,2,13,4!1號材料是鋼!2號材料是鋁!3號材料是銅!鋁是本次分析中的焊料, 它將鋼結(jié)

2、構(gòu)部分和銅結(jié)構(gòu)部分焊接起來!下面是在幾個溫度點下, 各材料的彈性模量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!假設(shè)各材料都是雙線性隨動硬化彈塑性本構(gòu)關(guān)系!下面是各材料在各溫度點下的屈服應(yīng)力和屈服后的彈性模量tb,bkin,1,5tbtemp,20,1tb

3、data,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,bkin,2,5tbtemp,20,1tbdata,1,800e6,0.102e11tbtemp,500,2tbdata,1,400e6,0.050e11tbtemp,1000,3tbdata,1, 70e6,0.008e11tbtemp,1500,4tbda

4、ta,1, 1e6,0.0001e11tbtemp,2000,5tbdata,1,0.1e6,0.00001e11!tb,bkin,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!材料密度(假設(shè)為常值)mp,dens,1,8030mp,dens,2,4850mp,dens,3,8900

5、! 熱膨脹系數(shù)(假設(shè)為常值)mp,alpx,1,1.78e-5mp,alpx,2,9.36e-6mp,alpx,3,1.66e-5!泊松比(假設(shè)為常值)mp,nuxy,1,0.29mp,nuxy,2,0.30mp,nuxy,3,0.30!熱傳導(dǎo)系數(shù)(假設(shè)為常值)mp,kxx,1,16.3mp,kxx,2,7.44mp,kxx,3,393!比熱(假設(shè)為常值)mp,c,1,502mp,c,2,544mp,c,3,385!熱膨脹系數(shù)(假設(shè)為常值)!由于該13號單元還有磁自由度, 此處假設(shè)一磁特性, 但并不影響我們所關(guān)心的結(jié)果mp,murx,1,1mp,murx,2,1mp,murx,3,1!假設(shè)焊料

6、(鋁)焊上去后的初始溫度是1500mp,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!劃分網(wǎng)格!esize,0.025type,2mat,2amesh,2 ! esize,0.05 !網(wǎng)格劃分出現(xiàn)問題type,1mat,1amesh,1!mat,3amesh,3eplot!/soluanty

7、pe,4! 瞬態(tài)分析trnopt,full!在模型的左邊界加位移約束!nsel,all*get,minx,node,mnloc,xnsel,s,loc,x,minxd,all,ux,0*get,miny,node,mnloc,ynsel,r,loc,y,minyd,all,uy,0!*假設(shè)模型的左右邊界處溫度始終保持在20攝氏度左右*!其他邊界條件如對流和輻射等均可施加，此處因為只是示意而已，故只施加恒溫邊界條件nsel,all*get,minx,node,mnloc,xnsel,s,loc,x,minxd,all,temp,20nsel,all*get,maxx,node,mxloc,xns

8、el,s,loc,x,maxxd,all,temp,20!由于第2個面是焊接所在區(qū)域，因此首先將該區(qū)域的單元“殺死”!nna=2esel,all*get,emax,elem,num,maxasel,s,area,nnaesla*get,nse,elem,count*dim,ne,nse*dim,nex,nse*dim,ney,nse*dim,neorder,nsemine=0!* !下面的do循環(huán)用于將焊料區(qū)的單元按其形心y坐標排序!以便后面模擬焊料由下向上逐步“生長”過程*do,i1,1,nseesel,u,elem,mine*get,nse1,elem,countii=0*do,i,1,e

9、max*if,esel(i),eq,1,thenii=ii+1ne(ii)=i*endif*enddo*do,i,1,nse1*get,ney(i),elem,ne(i),cent,y*get,nex(i),elem,ne(i),cent,x*enddominy=1e20minx=1e20*do,i,1,nse1*if,ney(i),lt,miny,thenminy=ney(i)minx=nex(i)mine=ne(i)*else*if,ney(i),eq,miny,then*if,nex(i),lt,minx,thenminy=ney(i)minx=nex(i)mine=ne(i)*endi

10、f*endif*endif*enddoneorder(i1)=mine*enddo!* max_tem=1500 !按照前面假設(shè)，焊料的初始溫度為1500dt1=1e-3 !用于建立初始條件的一個很小的時間段dt=5 !焊接一個單元所需的時間t=0!起始時間esel,alleplot/auto,1/replot*do,i,1,nseekill,neorder(i)esel,s,liveeplot*enddoallsel,alloutres,all,allic,all,temp,20kbc,1timint,0,structtimint,1,thermtimint,0,magtintp,0.005

11、,1,0.5,0.2!nsub1=2nsub2=40 !* *do,i,1,nseealive,neorder(i)esel,s,liveeplotesel,all!*下面的求解用于建立溫度的初始條件*t=t+dt1time,tnsubst,1*do,j,1,4d,nelem(neorder(i),j),temp,max_tem*enddosolve!*下面的求解用于保證初始的升溫速度為零*t=t+dt1time,tsolve!*下面的步驟用于求解溫度分布* *do,j,1,4ddele,nelem(neorder(i),j),temp*enddot=t+dt-2*dt1time,tnsubs

12、t,nsub1solve*enddot=t+50000 !*下面的步驟用于冷卻過程求解* time,tnsubst,nsub2solvesavefinish!后處理過程! /post1!*下面的一系列命令用于生成應(yīng)力的動畫文件*/seg,dele/cont,1,15,0,1200e6/16,1200e6/dscale,1,1.0avprin,0,0avres,1/seg,multi,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,1p

13、lnsol,s,eqv*do,j,1,nsub1set,(i-1)*3+3,jplnsol,s,eqv*enddo*enddo*do,i,1,nsub2set,(nse-1)*3+4,iplnsol,s,eqv*enddo/seg,off,stress1,0.1/anfile,save,stress1,avi!*下面的一系列命令用于生成溫度的動畫文件*/seg,dele/cont,1,15,0,1500/16,1500/dscale,1,1.0avprin,0,0avres,1/seg,multi,temp1,0.1esel,all*do,i,1,nseesel,u,elem,neorder(i)*end