| 查看: 381 | 回复: 0 | ||
[求助]
堆焊热分析 ——总是提示无法收敛，请高手指教
|
|
由于项目需要，想做一个堆焊过程仿真，观察基体温度的变化过程，但在第二个子步及t=0.4s时提示无法收敛，不知道是什么原因，麻烦高人帮忙看一下，指出问题所在，谢谢了！ /filname,weldtem /title,WELDING TEMPERATIUTE FIELD SIMULATION /units,si /prep7 et,1,solid70 keyopt,1,2,1 et,2,surf152 r,2 keyopt,2,4,1 keyopt,2,5,0 keyopt,2,8,3 mptemp,1,0,100,200,300,400,500 !定义材料1属性温度范围 mptemp,,1200,1250,1350,1400,1450,1500 mptemp,,1505,1705,1905,2105,2305,2500 mpdata,dens,1,1,7820,7800,7800,7800,7800,7800 !定义材料1密度(kg/m3) mpdata,dens,1,,7800,7800,7800,7800,7800,7800 mpdata,dens,1,,7800,7800,7800,7800,7800,7800 mpdata,kxx,1,1,52,50.7,48.6,46.1,42.3,38.9 !热传导系数(W/(m*K)) mpdata,kxx,1,,30,30,30,30,30,30 mpdata,kxx,1,,30,30,30,30,30,30 mpdata,c,1,1,450,469,481,508.5,536,569 !考虑相变潜热的比热容(J/kg*K) mpdata,c,1,,700,2172,5116,6589,8061,9533 mpdata,c,1,,9533,7757,5982,4206,2431,700 !基体材料单元类型 mpdata,dens,2,1,7850,7850,7800,7800,7800,7800 !定义材料2密度(kg/m3) mpdata,dens,2,,7800,7800,7800,7800,7800,7800 mpdata,dens,2,,7800,7800,7800,7800,7800,7800 mpdata,kxx,2,1,52,50.7,48.6,46.1,42.3,38.9 !热传导系数(W/(m*K)) mpdata,kxx,2,,30,30,30,30,30,30 mpdata,kxx,2,,30,30,30,30,30,30 mpdata,c,2,1,450,469,481,508.5,536,569 !考虑相变潜热的比热容(J/kg*K) mpdata,c,2,,700,2172,5116,6589,8061,9533 mpdata,c,2,,9533,7757,5982,4206,2431,700 pi=3.1415926 I=160 U=23 efficient=0.78 r=0.0045 Qm=3*I*U*efficient/(pi*r**2) kk=0.85 Qmax=Qm*kk ly=0.02 lx=0.02 lz=0.01 lsize=0.001 !单元尺寸 v=0.0025 tm=lx/v !移动的总时间 tinc=lsize/v !在每个单元的移动时间 max_time=1+tm/tinc max_x=1+lx/lsize !x方向节点个数 max_y=1+ly/lsize !y方向节点个数 /view,1,1,1,1 block,0,lx,-ly/2,ly/2,0.0003,lz+0.0003 !建模 block,0,lx,-ly/2,ly/2,0,0.0003 vglue,all numcmp,all save esize,lsize !单元的大小 mshkey,1 !Use mapped meshing type,1 mat,1 vmesh,1 esize,lsize/4 type,2 mat,2 mshkey,1 vmesh,2 /pnum,line,off /pnum,defa eplot finish --------------- *dim,flux2,table,max_x,max_y,max_time,x,y,time *do,k,1,max_time,1 *do,i,1,max_x,1 flux2(i,0,k)=(i-1)*lsize !经过x方向的距离 *enddo *do,j,1,max_y,1 flux2(0,j,k)=(j-1)*lsize !经过y方向的距离 *enddo *do,i,1,max_x,1 *do,j,1,max_y,1 xcenter=v*(k-1)*tinc ycenter=0 distance=sqrt(((i-1)*lsize-xcenter)**2+((j-1)*lsize-ycenter)**2) *if,distance,le,r,then flux2(i,j,k)=Qmax*exp(-3*distance**2/(r**2)) *else flux2(i,j,k)=0 *endif *enddo *enddo flux2(0,0,k)=(k-1)*tinc *enddo /solu antype,trans tunif,15 kbc,1 timint,on autots,off esel,s,type,,2 nsle,s sf,all,conv,-1 ! -1; current flow into the element face (IN ) allsel *do,t,1e-6,tm,tinc time,t autos,0 deltim,,,,1 kbc,0 outres,all,all tsres,erase sf,all,hflux,%flux2% solve *enddo Finish /seg,off,WJ,0.2 anim,1,1 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有5人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有5人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有5人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
