焊接.txt

ansys.com上的ansys命令流,主要为热,流体,接触分析方面的内容

这个例子的看点: 
1.建模. 
2.死活单元. 
3.相变. 
4.局部坐标系

/PREP7 
!********************* 
ET,1,SOLID70 !thermal solid 

MP,KXX,1,.5e-3 !weld on right 
MP,C,1,.2 
MP,DENS,1,.2833 
MPTEMP,1,0,2643,2750,2875,3000 !enthalpy as function of time 
MPDATA,ENTH,1,1,0,128.1,163.8,174.2,184.6 

MP,KXX,2,.5e-3 !solid material 
MP,C,2,.2 
MP,DENS,2,.2833 

MP,KXX,3,0.5e-3 !weld on left 
MP,C,3,.2 !for t=[0,1000]: only right weld, left weld has no therm conduct 
MP,DENS,3,.2833 !for t>1000: left weld is activated 

W=1.20 
TL=.34 
TWELD=.2 ! weld thickness 
TPIVOT=TL ! thickness of support bar 
TPHALF=TPIVOT/2 
WIDTH=W ! width of support bar 
WHALF=WIDTH/2 
LCYL=2.5 ! length of cylinder 
LHALF=LCYL/2 
RID=.50 ! inside radius of cylinder 
ROD=.75 ! outside radius of cylinder 
*AFUN,DEG 
ANG1=50-ASIN(TPHALF/ROD)-2*ASIN(TWELD/(2*ROD)) 
ANG2=50-ASIN(TPHALF/ROD) 
ANG3=50+ASIN(TPHALF/ROD) 
ANG4=50+ASIN(TPHALF/ROD)+2*ASIN(TWELD/(2*ROD)) 
WPST,-1 ,1 ,-10 ,10 ,.003 ,3 

/PNUM,VOLU,1 
/VIEW,1,-.2,.2,1 
CYLIND,RID,ROD,-LHALF,-WHALF,-90,ANG1 
CYLIND,RID,ROD,-WHALF,+WHALF,-90,ANG1 
CYLIND,RID,ROD,+WHALF,+LHALF,-90,ANG1 
CYLIND,RID,ROD,-LHALF,-WHALF,ANG1,ANG2 
CYLIND,RID,ROD,-WHALF,+WHALF,ANG1,ANG2 
CYLIND,RID,ROD,+WHALF,+LHALF,ANG1,ANG2 
CYLIND,RID,ROD,-LHALF,-WHALF,ANG2,ANG3 
CYLIND,RID,ROD,-WHALF,+WHALF,ANG2,ANG3 
CYLIND,RID,ROD,+WHALF,+LHALF,ANG2,ANG3 
CYLIND,RID,ROD,-LHALF,-WHALF,ANG3,ANG4 
CYLIND,RID,ROD,-WHALF,+WHALF,ANG3,ANG4 
CYLIND,RID,ROD,+WHALF,+LHALF,ANG3,ANG4 
CYLIND,RID,ROD,-LHALF,-WHALF,ANG4,90 
CYLIND,RID,ROD,-WHALF,+WHALF,ANG4,90 
CYLIND,RID,ROD,+WHALF,+LHALF,ANG4,90 

/TYPE,1,4 
WPROT,50 ! vol 16 (used to trim welds) 
BLOCK,.6,rod*cos(asin(tphalf/rod))+tweld,-tphalf,tphalf,-whalf,whalf 
CSWPLA,11,0 
KSEL,S,LOC,X,.75*COS(ASIN(TPHALF/ROD)) 
KSEL,R,LOC,Y,-TPHALF 
KSEL,R,LOC,Z,LHALF 
*GET,KNUM1,KP,,NUM,MAX 
KSEL,S,LOC,X,.75*COS(ASIN(TPHALF/ROD)) 
KSEL,R,LOC,Y,TPHALF 
KSEL,R,LOC,Z,LHALF 
*GET,KNUM2,KP,,NUM,MAX 
KSEL,ALL 
WPRO,-50 
CSYS,0 
LOCAL,12,0,KX(KNUM1),KY(KNUM1) 
CSYS,0 
LOCAL,13,0,KX(KNUM2),KY(KNUM2) 
WPCSYS,-1,12 
VSEL,U,,,1,16 
CYLIND,,TWELD,-WHALF,WHALF,(ANG1+ANG2)/2-90,70 
/DIST,1,.6 
/FOCUS,1,.47,.59 
WPCSYS,-1,13 
CYLIND,,TWELD,-WHALF,WHALF,35,90+(ANG3+ANG4)/2 
VSEL,ALL 
BOPT,keep,yes, 
bopt,version,rv52 
VSEL,S,,,5,17,12 
VSBV,17,5 
VSEL,A,,,16,19,3 
VSBV,19,16 
VSEL,A,,,11,18,7 
VSBV,18,11 
VSEL,A,,,16,21,5 
VSBV,21,16 
VSEL,ALL 
WPSTYLE,default 
/FOCUS,,.5,.5 
/DIST,1,.6 
CYLIN,.6,SQRT((ROD+.01)**2+TPHALF**2),-2,+2,10,90 
BOPT,keep,no 
VSEL,S,,,20,22,2 
VSEL,A,,,23 
VSEL,ALL 
VOVLAP,20,22,23 
VDELE,16,21,1,1 
VDELE,28,,,1 
/TYPE,1,4 

CSYS,0 
WPCSYS,-1 
WPROT,50 

NUMSTR,VOLU,29 
block,ROD+.01,ROD*COS(ASIN(TPHALF/ROD))+TWELD,-TPHALF,TPHALF,-WHALF,WHALF 
NUMSTR,VOLU,30 
BLOCK,ROD*COS(ASIN(TPHALF/ROD))+TWELD,3,-TPHALF,TPHALF,-WHALF,WHALF 
NUMMRG,KP 
VSEL,S,,,24,26,2 ! right hand weld (make mat 1) 
ESIZE,.1 
VMESH,ALL 
MAT,3 
vsel,s,,,25,27,2 ! left hand weld (make mat 3 - dead) 
vmesh,all 
ALLSEL,ALL 
VSEL,S,,,24,27,1 
VSEL,A,,,30 
VSEL,INVE 
ESIZE,.1 
MAT,2 
VMESH,ALL 
ESIZE,.2 
LSEL,S,,,198,199 
LESIZE,ALL,,,8,5,1 
LSEL,S,,,196,201,5 
LESIZE,ALL,,,8,.2,1 
ALLSEL,ALL 
VMESH,30 
NUMCMP,all 
WPSTYLE,,,,,,,,0 
esel,s,mat,,3 
ekill,all 
allsel,all 
SAVE,mesh,DB 
FINI 
! 
! 
!********************* 
/SOLU 
!********************* 
ANTYPE,TRANS 

!Initialization of temperatures 
!---------------------------- 
TIMINT,OFF 
ESEL,S,MAT,,1 !weld on right 
NSLE 
D,ALL,TEMP,3000 !To(weld)=3000 
NSEL,INVE 
D,ALL,TEMP,70 !To(steel)=70 
TIME,1 
KBC,0 
NSUBST,2,2,2 
OUTRES,BASIC,2 
ALLSEL,all 
SOLVE 
SAVE,therm,db 

!Phase change for weld on right (time varies from 1 to 100 and from 100 to 1000) 
!---------------------------------------------------------------------------------- 
OUTRES,NSOL,1 
DDELE,ALL,TEMP 
TIMINT,ON 
TINTP,,,,1 
TIME,100 !phase change 
DELTIME,1,.5,10 
AUTOTS,ON 
ASEL,S,EXT 
ASEL,U,LOC,X,0 
SFA,ALL,,CONV,5E-5,70 
ALLSEL,ALL 
SOLVE 

TIME,1000 !solidification 
DELTIME,50,10,500 
AUTOTS,ON 
SOLVE 
FINI 

!Phase change for weld on right (time varies from 1e3 to 1100 and from 1100 to 2000) 
!--------------------------------------------------------------------------------------- 
/PREP7 
EALIVE,ALL 
ALLSEL,all 
FINI 

/SOLU 
ANTYPE,TRANS,REST 
vsel,s,,,17 
vsel,a,,,19 
nslv,,1 
d,all,temp,3000 
time,1001 
nsubst,2,2,2 
allsel,all 
solve 

ANTYPE,TRANS,REST 
DDELE,ALL,TEMP 
TIME,1100 
DELTIME,1,.5,10 
SOLVE 
TIME,2000 
DELTIME,100,10,200 
SOLVE 
FINI