一条焊缝的命令流程序，请高手指教！
大家好我编写了一条焊缝的焊接命令流，但是只是加热阶段的程序，冷却阶段的老是处错误，希望有高手指教下，本人将不胜感激！补充一点，我用的是直接法耦合的，下面是我焊接加热阶段的命令流程序，运行没有错误。（注意记得glue一下）/PreP7mPtemP,l,20,250,500,750,1000,1500,1700,2500mpdate,kxx,l,l,50,47,40,27,30,35,140,142 mpdata,dens,l,l,7820,7700,7610,7550,7490,7350,7300,7090 mpdata,c,1,l,460,480,530,675,670,660,780,820 MPDATA,PRXY,l,l,0.28,0.29,0.31,0.35,0.4,0.45,0.48,0.5 mPdata,alPx,l,l,1.10E-05,1.22E-05,1.39E-05,l.48E-05,l.34E-05,l.33E-05,1.32E-05,1.3lE-05 mPdata,ex,l,l,2.05E+11,1.87E+11,1.50E+11,7.00E+10,2.00E+10,1.90E+10,1.80E+10,l.20E+10 MPTEMP,,,,,,,,MPTEMP,l,20MPTEMP,2,1533MPTEMP,3,1590MPTEMP,4,1670MPDATA,ENTH,l,,0MPDATA,ENTH,l,,7.5E+009MPDATA,ENTH,l,,9.6E+009MPDATA,ENTH,l,,1.1E+010MPDATA,MURX,1,,1MPTEMP,,,,,,,,MPTEMP,1,20 MPTEMP,2,250MPTEMP,3,500MPTEMP,4,750MPTEMP,5,1000
MPTEMP,6,1500
MPTEMP,7,1700
MPTEMP,8,2500
MPDATA,HF,1,,1
MPDATA,HF,1,,3.5MPDATA,HF,1,,5.2MPDATA,HF,1,,10 MPDATA,HF,1,,15 MPDATA,HF,1,,30 MPDATA,HF,1,,31 MPDATA,HF,1,,35 !下面在各温度点下的屈服应力和屈服后的弹性模量TB,BISO,l,6,2,TBTEMP,20TBDATA,,2.20E+08,2.05E+10,,,,TBTEMP,500TBDATA,,1.30E+08,l.50E+10,,,,TBTEMP,750TBDATA,,4.00E+07,7.00E+09,,,,TBTEMP,1000TBDATA,,2.50E+07,2.00E+09,,,,TBTEMP,1700TBDATA,,1.00E+06,1.80E+09,,,,TBTEMP,2500TBDATA,,1.00E+05,1.20E+09,,,,SAVE /prep7
/units,si
Rb=0.002 Length=0.12width=0.03 height=0.003et,l,solid5et,2,solid98block,0,length,0,width,0,height block,0,length,0,3*Rb,0,heightblock,0,length,0,6*Rb,0,height block,0,length,0,12*Rb,0,heightVOVLAP,ALL /view,1,1,1,1 !下面定义焊缝区的单元尺寸Lesize,4,rbLesize,5,rbLesize,9,rbLesize,14,rbLesize,50,rbLesize,10,rbLesize,52,rbLesize,23,rbLesize,51,rbLesize,9,rbLesize,49,rbLesize,24,rb!定义过渡区的单元尺寸Lesize,31,2*rbLesize,26,2*rb!定义远离焊缝区的单元尺寸Lesize,43,4*rbLesize,38,4*rbLesize,2,4*rbLesize,7,4*rbLesize,53,4*rbLesize,54,4*rbLesize,55,4*rbLesize,56,4*rbLesize,11,1.5*rbLesize,12,1.5*rbLesize,47,1.5*rbLesize,48,1.5*rbSave记得glue一下 Type,1Mshkey,1 Mshape,0Vmesh,5 Vmesh,6Type,2Mshkey,0Mshape,1Vmesh,7 Vmesh,8Save
/SOLUDA,3, /GODL,53,,ux,DL,61,,UX,DL,57,,UX,DL,49,,UX,DL,4,,uz,
DL,9,,uy, /Solu!时间步长设置Antype,trans,new!Nlgeom,on!Timint,1,struct!Timint,1,therm! Timint,0,mag!Timint,0,electTintP,0.005,,,-1,0.5,0.2Tref, 25Nropt,auto,,Autos,onOutres,,allv=0.01tinc=rb/v!移动热源设置U=25I=180ef =0.75P=25*180*efr=0.006 qmax=p/(3.14159*r**2) Max_time=length/rb+1max_x=length/rb+1max_y=width/rb+1*dim,flux2,table,max_x,max_y,max_time,x,y,time*do,k,1,max_time,1 *do,i,1,max_x,1Flux2(i,0,k)=(i-1)*rb*enddo*do,j,1,max_y,1Flux2(0,j,k)= (j-1)*rb*enddo*do,i,1,max_x,1*do,j,1,max_y,1!设置热源中心位置坐标xcenter=v*(k-1)*tinc ycenter=0distance=sqrt(((i-1)*rb-xcenter)**2+((j-1)*rb-ycenter)**2) *if,distance,le,r,then Flux2(i,j,k)=qmax/EXP(3*DISTANCE**2/(R**2)) *elseFlux2(i,j,k)=0*endif*enddo*enddoFlux2(0,0,k)=(k-1)*tincNsel,s,loc,z,height Sf,all,hflux,%flux2%Allsel,allNsel,s,loc,x,0 Sf,all,conv,10,20Allsel,allNsel,s,loc,x,length Sf,all,conv,10,20Allsel,allNsel,s,loc,y,0Sf,all,conv,10,20Allsel,allNsel,s,loc,y,widthSf,all,conv,10,20Allsel,allNsel,s,loc,z,0Sf,all,conv,10,20Allsel,allTime,k*rb/vNsubst,5 solve*enddo
save下面冷却阶段的命令流哪位高手能给指教一下，谢谢了！
页:
[1]