一个采用双椭球热源焊接温度场的例子
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
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.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,,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