【分享】焊接的相变分析（补上次发过的）
这是ansys.net上唯一的关于welding的inp文件，大家共同学习：!注意开始文件名是test，不然会有问题
!FEAmodelofaweld
!Themodelconsistsoftwopartsthatareweldedtogetherthrough2welds.
!Thisisaphasechangeanalysis(transient).
!Thefirstpartconsistsofsimulatingthephasechangeandsolidificationofoneweld.
!Consequently,the2ndweldisactivatedandthetempdistributioniscalculated.
!
FINISH
/CLEAR
/TITLE,Solidificationof2welds
!
!*********************
/PREP7
!*********************
ET,1,SOLID70!thermalsolid
MP,KXX,1,.5e-3!weldonright
MP,C,1,.2
MP,DENS,1,.2833
MPTEMP,1,0,2643,2750,2875,3000!enthalpyasfunctionoftime
MPDATA,ENTH,1,1,0,128.1,163.8,174.2,184.6
MP,KXX,2,.5e-3!solidmaterial
MP,C,2,.2
MP,DENS,2,.2833
MP,KXX,3,0.5e-3!weldonleft
MP,C,3,.2!fort=:onlyrightweld,leftweldhasnothermconduct
MP,DENS,3,.2833!fort>1000:leftweldisactivated
W=1.20
TL=.34
TWELD=.2!weldthickness
TPIVOT=TL!thicknessofsupportbar
TPHALF=TPIVOT/2
WIDTH=W!widthofsupportbar
WHALF=WIDTH/2
LCYL=2.5!lengthofcylinder
LHALF=LCYL/2
RID=.50!insideradiusofcylinder
ROD=.75!outsideradiusofcylinder
*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!vol16(usedtotrimwelds)
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!righthandweld(makemat1)
ESIZE,.1
VMESH,ALL
MAT,3
vsel,s,,,25,27,2!lefthandweld(makemat3-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
!Initializationoftemperatures
!----------------------------
TIMINT,OFF
ESEL,S,MAT,,1!weldonright
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
!Phasechangeforweldonright(timevariesfrom1to100andfrom100to1000)
!----------------------------------------------------------------------------------
OUTRES,NSOL,1
DDELE,ALL,TEMP
TIMINT,ON
TINTP,,,,1
TIME,100!phasechange
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
!Phasechangeforweldonright(timevariesfrom1e3to1100andfrom1100to2000)
!---------------------------------------------------------------------------------------
/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
/COPY,test,rth,,temp,rth
FINI
!*********************
/POST1
!*********************
FILE,temp,rth
SET,,,1,,1,,
PLNS,TEMP,
ANTIME,20,0.5,,1,2,1001,1005
!
!
回复： 【分享】焊接的相变分析（补上次发过的）
inp就是input文件在工作目录下，新建一个记事本，把名字改为abc.inp，inp是扩展名。然后用记事本编辑这个inp文件，把命令流帖上去，保存。在ANSYS里，命令栏输入/input,abc,inp
回复： 【分享】焊接的相变分析（补上次发过的）
非常好a回复： 【分享】焊接的相变分析（补上次发过的）
用txt文件就可以不一定非得inp
页:
[1]