一个采用双椭球热源焊接温度场的例子.txt

我做的一些ansys的练习,命令流,对学习ansys很有帮助

finish
/clear
/BATCH    
WPSTYLE,,,,,,,,0
/FILNAME,Laser Beam Welding,0  
/TITLE,laser Beam Welding  
/units,si 
!!!参数定义
L=0.02
B1=0.003
B2=0.006
B3=0.02
H=0.0032 

D=0.001                !焊缝宽度 
ah=D/2
chf=2*ah
chb=8*ah
bh=0.005
rf=2*chf/(chf+chb)
rb=2*chb/(chf+chb)
N=0.27                     !功率有效系数 
W0=2200                    !功率
W=W0*N                     !有效功率 

/PREP7  
!定义导热系数
ET,1,SOLID70 
MPTEMP,1,20
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,600
MPTEMP,8,700
MPTEMP,9,800
MPTEMP,10,900
MPTEMP,11,1000
MPTEMP,12,1200   
MPTEMP,13,1540
MPTEMP,14,1700
MPTEMP,15,2000
MPTEMP,16,3000
MPDATA,KXX,1,,6 
MPDATA,KXX,1,,6.5  
MPDATA,KXX,1,,7
MPDATA,KXX,1,,8
MPDATA,KXX,1,,9
MPDATA,KXX,1,,10
MPDATA,KXX,1,,11 
MPDATA,KXX,1,,12
MPDATA,KXX,1,,13
MPDATA,KXX,1,,13.5
MPDATA,KXX,1,,13.8
MPDATA,KXX,1,,13.9
MPDATA,KXX,1,,14
MPDATA,KXX,1,,14.2
MPDATA,KXX,1,,14.4
MPDATA,KXX,1,,14.4
!定义比热容
MPTEMP,1,20
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,600
MPTEMP,8,700
MPTEMP,9,800
MPTEMP,10,900
MPTEMP,11,1000
MPTEMP,12,1200   
MPTEMP,13,1540
MPTEMP,14,1700
MPTEMP,15,2000
MPTEMP,16,3000
MPDATA,C,1,,600
MPDATA,C,1,,620
MPDATA,C,1,,640
MPDATA,C,1,,660
MPDATA,C,1,,680
MPDATA,C,1,,700
MPDATA,C,1,,725
MPDATA,C,1,,750
MPDATA,C,1,,775
MPDATA,C,1,,800
MPDATA,C,1,,805
MPDATA,C,1,,810
MPDATA,C,1,,830
MPDATA,C,1,,850
MPDATA,C,1,,900
MPDATA,C,1,,900 
!定义密度
MP,DENS,1,4450 

!!!建立几何模型 
BLOCK,0,B1,0,L,0,H,
BLOCK,B1,B2,0,L,0,H,
BLOCK,B2,B3,0,L,0,H,  
VGLUE,all 
!!!划分网格

esize,h/16,0
VSWEEP,1 
esize,h/2,0
VSWEEP,5 

MSHAPE,1,3D 
MSHKEY,0
VMESH,4   
  
FINISH  
!!!施加载荷 
STEP=0.002                !
N=L/STEP
V=0.004                    !速度 
*DO,I,0,N 
YY=I*STEP
T=(YY+STEP)/V 

/SOLU 
ANTYPE,4
TRNOPT,FULL 
TUNIF,293, 
OUTRES,NSOL,LAST
TIME,T 
NROPT,FULL,,ON
AUTOTS,1
NSUBST,5, , ,1   !!!
KBC,1 
lnsrch,1     !!!
TSRES,ERASE 
TINTP,,,,1,,,,,!!!
TIMINT,0,struct
TIMINT,1,THERM
TIMINT,0,MAG 

!!!添加第一个载荷步（焊接阶段，有对流和热源） 
NSEL,S,EXT  
NSEL,U,LOC,X,0 
ASEL,S, , ,5  
ASEL,U,LOC,X,0   
LSLA,S    
NSLL,A,1
SF,ALL,CONV,200,293
ALLSEL,ALL
!NSEL,ALL
!*GET,CENTER,NODE,,MNLOC,X
!NSEL,S,LOC,X,CENTER
!SF,ALL,HFLUX,0 

!!!双椭球热源的加载 
!!!前半球节点的选取 

LOCAL,11,2,0,YY,0, , , ,3,8,  
NSEL,S,LOC,X,,D  
NSEL,R,LOC,Y,,90
NSEL,R,LOC,Z,,90
CM,FRONT,NODE 

!!!!!!!!!!!!!!!!!!!!!!!!!!!
LOCAL,13,0,0,YY,0, , , ,1,1, 
!!!前半球函数、表格定义 
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNC_C4
*DEL,_FNC_C5
*DEL,_FNCCSYS
*DEL,'TFRONT'    !!!!!!!! 
*SET,_FNCNAME,'TFRONT'  
*DIM,_FNC_C1,,1 
*DIM,_FNC_C2,,1 
*DIM,_FNC_C3,,1 
*DIM,_FNC_C4,,1 
*DIM,_FNC_C5,,1  
*SET,_FNC_C1(1),rf           !!!rf
*SET,_FNC_C2(1),w            !!!Q
*SET,_FNC_C3(1),chf          !!!chf
*SET,_FNC_C4(1),ah           !!!ah
*SET,_FNC_C5(1),bh           !!!bh
*SET,_FNCCSYS,13            
*DIM,%_FNCNAME%,TABLE,7,38,1,,,,%_FNCCSYS%  
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0 
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
*SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)%
*SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)%
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 16, 3, 0, 0, 0 
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 6, 0, 0, -1 
*SET,%_FNCNAME%(0,3,1),   0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, -3, 3, 17
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -1, 3, 18
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 0, 0, 0, 0  
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, 0, 0, -1 
*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -1, 2, -3
*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 2, 0, 0, 2  
*SET,%_FNCNAME%(0,10,1), 0.0, -3, 0, 1, 2, 17, -1   
*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 2, 0, 0, 19
*SET,%_FNCNAME%(0,12,1), 0.0, -5, 0, 1, 19, 17, -1  
*SET,%_FNCNAME%(0,13,1), 0.0, -1, 0, 1, -3, 4, -5   
*SET,%_FNCNAME%(0,14,1), 0.0, -3, 0, 2, 0, 0, 3 
*SET,%_FNCNAME%(0,15,1), 0.0, -5, 0, 1, 3, 17, -3   
*SET,%_FNCNAME%(0,16,1), 0.0, -3, 0, 2, 0, 0, 20
*SET,%_FNCNAME%(0,17,1), 0.0, -6, 0, 1, 20, 17, -3  
*SET,%_FNCNAME%(0,18,1), 0.0, -3, 0, 1, -5, 4, -6   
*SET,%_FNCNAME%(0,19,1), 0.0, -5, 0, 1, -1, 1, -3   
*SET,%_FNCNAME%(0,20,1), 0.0, -1, 0, 2, 0, 0, 4 
*SET,%_FNCNAME%(0,21,1), 0.0, -3, 0, 1, 4, 17, -1   
*SET,%_FNCNAME%(0,22,1), 0.0, -1, 0, 2, 0, 0, 21
*SET,%_FNCNAME%(0,23,1), 0.0, -6, 0, 1, 21, 17, -1   
*SET,%_FNCNAME%(0,24,1), 0.0, -1, 0, 1, -3, 4, -6   
*SET,%_FNCNAME%(0,25,1), 0.0, -3, 0, 1, -5, 1, -1   
*SET,%_FNCNAME%(0,26,1), 0.0, -1, 0, 3, 0, 0, -3
*SET,%_FNCNAME%(0,27,1), 0.0, -5, 0, 1, -1, 3, -3   
*SET,%_FNCNAME%(0,28,1), 0.0, -1, 0, 1, -4, 3, -5   
*SET,%_FNCNAME%(0,29,1), 0.0, -1, 7, 1, -1, 0, 0
*SET,%_FNCNAME%(0,30,1), 0.0, -3, 0, 1, -2, 3, -1   
*SET,%_FNCNAME%(0,31,1), 0.0, -1, 0, 1, 20, 3, 21   
*SET,%_FNCNAME%(0,32,1), 0.0, -2, 0, 1, -1, 3, 19   
*SET,%_FNCNAME%(0,33,1), 0.0, -1, 0, 3.14159265358979310, 0, 0, -2  
*SET,%_FNCNAME%(0,34,1), 0.0, -4, 0, 1, -2, 3, -1   
*SET,%_FNCNAME%(0,35,1), 0.0, -1, 16, 3.14159265358979310, 0, 0, 0  
*SET,%_FNCNAME%(0,36,1), 0.0, -2, 0, 1, -4, 3, -1   
*SET,%_FNCNAME%(0,37,1), 0.0, -1, 0, 1, -3, 4, -2   
*SET,%_FNCNAME%(0,38,1), 0.0, 99, 0, 1, -1, 0, 0 
!!!!!!!!!!!!!!!!!! 
BF,FRONT,HGEN, %TFRONT%     
CMDELE,FRONT   
ALLSEL,ALL 
!!!后半球节点的选取 
LOCAL,12,2,0,YY,0, , , ,6,8,   
NSEL,S,LOC,X,,D   
NSEL,R,LOC,Y,-90,0
NSEL,R,LOC,Z,,90
CM,BEHIND,NODE 

!!!后半球加载 
!!!!!!!!!!!!!!!!!!!
LOCAL,14,0,0,YY,0, , , ,1,1, 
!!!后半球函数、表格定义 
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNC_C4
*DEL,_FNC_C5
*DEL,_FNCCSYS 
*DEL,'TBEHIND'  !!!!!!!!! 
*SET,_FNCNAME,'TBEHIND' 
*DIM,_FNC_C1,,1 
*DIM,_FNC_C2,,1 
*DIM,_FNC_C3,,1 
*DIM,_FNC_C4,,1 
*DIM,_FNC_C5,,1 
*SET,_FNC_C1(1),rb                   !!!rb
*SET,_FNC_C2(1),w                   !!!Q
*SET,_FNC_C3(1),chb                !!!chb
*SET,_FNC_C4(1),ah                 !!!ah
*SET,_FNC_C5(1),bh                 !!!bh
*SET,_FNCCSYS,14                        !!!   0
*DIM,%_FNCNAME%,TABLE,7,38,1,,,,%_FNCCSYS%  
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0 
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
*SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)%
*SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)%
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 16, 3, 0, 0, 0 
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 6, 0, 0, -1 
*SET,%_FNCNAME%(0,3,1),   0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, -3, 3, 17
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -1, 3, 18
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 0, 0, 0, 0  
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, 0, 0, -1 
*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -1, 2, -3
*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 2, 0, 0, 2  
*SET,%_FNCNAME%(0,10,1), 0.0, -3, 0, 1, 2, 17, -1   
*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 2, 0, 0, 19
*SET,%_FNCNAME%(0,12,1), 0.0, -5, 0, 1, 19, 17, -1  
*SET,%_FNCNAME%(0,13,1), 0.0, -1, 0, 1, -3, 4, -5   
*SET,%_FNCNAME%(0,14,1), 0.0, -3, 0, 2, 0, 0, 3 
*SET,%_FNCNAME%(0,15,1), 0.0, -5, 0, 1, 3, 17, -3   
*SET,%_FNCNAME%(0,16,1), 0.0, -3, 0, 2, 0, 0, 20
*SET,%_FNCNAME%(0,17,1), 0.0, -6, 0, 1, 20, 17, -3  
*SET,%_FNCNAME%(0,18,1), 0.0, -3, 0, 1, -5, 4, -6   
*SET,%_FNCNAME%(0,19,1), 0.0, -5, 0, 1, -1, 1, -3   
*SET,%_FNCNAME%(0,20,1), 0.0, -1, 0, 2, 0, 0, 4 
*SET,%_FNCNAME%(0,21,1), 0.0, -3, 0, 1, 4, 17, -1   
*SET,%_FNCNAME%(0,22,1), 0.0, -1, 0, 2, 0, 0, 21
*SET,%_FNCNAME%(0,23,1), 0.0, -6, 0, 1, 21, 17, -1  
*SET,%_FNCNAME%(0,24,1), 0.0, -1, 0, 1, -3, 4, -6   
*SET,%_FNCNAME%(0,25,1), 0.0, -3, 0, 1, -5, 1, -1   
*SET,%_FNCNAME%(0,26,1), 0.0, -1, 0, 3, 0, 0, -3
*SET,%_FNCNAME%(0,27,1), 0.0, -5, 0, 1, -1, 3, -3   
*SET,%_FNCNAME%(0,28,1), 0.0, -1, 0, 1, -4, 3, -5   
*SET,%_FNCNAME%(0,29,1), 0.0, -1, 7, 1, -1, 0, 0
*SET,%_FNCNAME%(0,30,1), 0.0, -3, 0, 1, -2, 3, -1   
*SET,%_FNCNAME%(0,31,1), 0.0, -1, 0, 1, 20, 3, 21   
*SET,%_FNCNAME%(0,32,1), 0.0, -2, 0, 1, -1, 3, 19   
*SET,%_FNCNAME%(0,33,1), 0.0, -1, 0, 3.14159265358979310, 0, 0, -2  
*SET,%_FNCNAME%(0,34,1), 0.0, -4, 0, 1, -2, 3, -1   
*SET,%_FNCNAME%(0,35,1), 0.0, -1, 16, 3.14159265358979310, 0, 0, 0  
*SET,%_FNCNAME%(0,36,1), 0.0, -2, 0, 1, -4, 3, -1   
*SET,%_FNCNAME%(0,37,1), 0.0, -1, 0, 1, -3, 4, -2   
*SET,%_FNCNAME%(0,38,1), 0.0, 99, 0, 1, -1, 0, 0 
BF,BEHIND,HGEN, %TBEHIND%    
CMDELE,BEHIND
ALLSEL,ALL 
SOLVE 
BFDELE,all,HGEN 
*ENDDO 
SAVE
!!!添加冷却载荷第一步
CSYS,0 
NSEL,S,EXT  
SF,all,CONV,200,293 
ALLSEL,ALL  
TIME,L/V+1
NSUBST,5, , ,1 
!!!DELTIM,0.05,0.01,0.1,1  
SOLVE  
!!!添加冷却载荷第二步
CSYS,0 
TIME,L/V+10 
NSUBST,5, , ,1 
!!!DELTIM,0.2,0.1,0.4,1
SOLVE  

/post26 
nsol,6,node(0.000,5*L/10,0.000),temp,,origin
!!!在表面沿着远离焊缝中心（+x）方向取点 
nsol,7,node(0.0005,2*L/4,0),temp,,top1
nsol,8,node(0.001,2*L/4,0),temp,,top2
nsol,9,node(0.0012,2*L/4,0),temp,,top3
nsol,10,node(0.0015,2*L/4,0),temp,,top4
plvar,6,7,8,9,10,
save
finish