激光焊接ANSYS命令流脉冲热源Word文档下载推荐.docx
《激光焊接ANSYS命令流脉冲热源Word文档下载推荐.docx》由会员分享，可在线阅读，更多相关《激光焊接ANSYS命令流脉冲热源Word文档下载推荐.docx（16页珍藏版）》请在冰豆网上搜索。
MPTEMP,,,,,,,,
MPTEMP,1,293
MPTEMP,2,373
MPTEMP,3,473
MPTEMP,4,573
MPTEMP,5,673
MPTEMP,6,773
MPTEMP,7,873
MPTEMP,8,973
MPTEMP,9,1073
MPTEMP,10,1173
MPTEMP,11,1273
MPTEMP,12,1373
MPTEMP,13,1500
MPTEMP,14,1600
MPTEMP,15,1700
MPTEMP,16,1800
MPTEMP,17,1900
MPTEMP,18,2000
MPTEMP,19,2200
MPTEMP,20,2400
MPTEMP,21,2500
MPTEMP,22,2600
MPTEMP,23,2700
MPTEMP,24,2800
MPTEMP,25,2900
MPTEMP,26,3000
MPDATA,KXX,1,,10.8
MPDATA,KXX,1,,11.9
MPDATA,KXX,1,,12.9
MPDATA,KXX,1,,14.2
MPDATA,KXX,1,,15.8
MPDATA,KXX,1,,17.6
MPDATA,KXX,1,,19.5
MPDATA,KXX,1,,21.3
MPDATA,KXX,1,,23.2
MPDATA,KXX,1,,24.8
MPDATA,KXX,1,,26.3
MPDATA,KXX,1,,27.4
MPDATA,KXX,1,,28.2
MPDATA,KXX,1,,41.1
定义密度
MPTEMP,1,0
MPDATA,DENS,1,,8890
定义比热容
MPDATA,C,1,,423
MPDATA,C,1,,454
MPDATA,C,1,,472
MPDATA,C,1,,489
MPDATA,C,1,,505
MPDATA,C,1,,522
MPDATA,C,1,,538
MPDATA,C,1,,553
MPDATA,C,1,,569
MPDATA,C,1,,584
MPDATA,C,1,,598
MPDATA,C,1,,613
MPDATA,C,1,,830
MPDATA,C,1,,1650
MPDATA,C,1,,660
MPDATA,C,1,,670
MPDATA,C,1,,680
MPDATA,C,1,,688
MPDATA,C,1,,695
MPDATA,C,1,,699
MPDATA,C,1,,703
MPDATA,C,1,,705
建立模型
S_DIS=0.0005
RECTNG,S_DIS,0.050+S_DIS,0,0.0008
RECTNG,S_DIS,0.050+S_DIS,0.0008,0.0018
RECTNG,S_DIS,0.050+S_DIS,0.0018,0.005
RECTNG,S_DIS,0.050+S_DIS,0.005,0.050
aglue,all
numcmp,all
划分网格
------------------------面4
LESIZE,12,,,20
LESIZE,13,,,20
LESIZE,7,,,20
MSHKEY,1
MSHAPE,0,2D
AMESH,4
-------------------------面3
LESIZE,10,,,2
LESIZE,11,,,2
AMESH,3
LREFINE,5,,,1,1,OFF
-------------------------面2
LESIZE,8,,,2
LESIZE,9,,,2
AMESH,2
LREFINE,3,,,1,1,OFF
------------------------面1
LESIZE,2,,,4,
LESIZE,4,,,4,
AMESH,1
ALLSEL,ALL
NUMCMP,ALL
/REPLOT
FINISH
/SOLU
ANTYPE,TRANS
TRNOPT,FULL
TUNIF,293
nlgeom,off
nropt,full
KK=3
N_ONCE=TT/TINC
MAX_TIME=N_ONCE+1
MAX_X=（LENGTH+2*S_DIS）/LSIZE+1
MAX_Y=2*L_ZONE/LSIZE+1
H_MAX_Y=L_ZONE/LSIZE+1
------------------------定义载荷并求解-------------------------------
SFA,2,,CONV,400,293
SFA,3,,CONV,250,293
SFA,4,,CONV,50,293
*DIM,FLUX2,TABLE,MAX_X,MAX_Y,MAX_TIME,X,Y,TIME
*DO,K,1,MAX_TIME,1
*DO,I,1,MAX_X,1
FLUX2（I,0,K）=（I-1）*LSIZE
*ENDDO
*DO,J,1,MAX_Y,1
FLUX2（0,J,K）=-（H_MAX_Y-J）*LSIZE
XCENTER=（K-1）*Vlaser/F
YCENTER=0
XCOORD=（I-1）*LSIZE
YCOORD=-（H_MAX_Y-J）*LSIZE
DISTANCE=SQRT（（XCOORD-XCENTER）**2+（YCOORD-YCENTER）**2）
*IF,DISTANCE,LE,Rlaser,THEN
FLUX2（I,J,K）=QMAX/EXP（KK*（DISTANCE**2）/（Rlaser**2））
*ELSE
FLUX2（I,J,K）=0
*ENDIF
FLUX2（0,0,K）=（K-1）*TINC
OUTRES,BASIC,ALL
KBC,1
TSRES,ERASE
TIMINT,ON
LNSRCH,ON
TM1=0.005+（K-1）*TINC
TIME,TM1
SFA,1,,HFLUX,%FLUX2%
AUTOTS,-1
NSUBST,4
/GST,ON
/REP,FAST
/PSF,HFLUX,,1!
表面上热流加载用带颜色填充曲面表示
SOLVE
SFADELE,1,,HFLUX
TM2=K*TINC
TIME,TM2
SFA,1,,CONV,400,293
SFADELE,1,,CONV
SAVE
***********************************
冷却阶段
*DO,k,1,30,1
TIME,TT+5*k!
载荷步结束时间
NSUBST,25
outres,all,all
finish
/FILNAME,stress3,1
/PREP7
et,1,solid45!
把热单元转化成结构单元
R,1,0.0005
------------------------设置结构分析的材料属性-----------------------
MPTEMP,4,513
MPTEMP,5,573
MPTEMP,6,673
MPTEMP,7,773
MPTEMP,8,873
MPTEMP,9,973
MPTEMP,10,1073
MPTEMP,11,1173
MPTEMP,12,1273
MPTEMP,13,1373
MPTEMP,14,1423
MPTEMP,15,1473
MPTEMP,16,1596
MPTEMP,17,1673
MPTEMP,18,1873
MPTEMP,19,2073
MPTEMP,20,2273
MPTEMP,21,2473
MPTEMP,22,2673
MPTEMP,23,2873
MPTEMP,24,3073
MPTEMP,25,3273
MPDATA,ALPX,1,,0.184E-5
MPDATA,ALPX,1,,1.17E-5
MPDATA,ALPX,1,,1.21E-5
MPDATA,ALPX,1,,1.28E-5
MPDATA,ALPX,1,,1.31E-5
MPDATA,ALPX,1,,1.35E-5
MPDATA,ALPX,1,,1.40E-5
MPDATA,ALPX,1,,1.47E-5
MPDATA,ALPX,1,,1.55E-5
MPDATA,ALPX,1,,1.60E-5
MPDATA,ALPX,1,,1.65E-5
MPDATA,ALPX,1,,1.8E-5
MPTEMP,13,1473
MPTEMP,14,1596
MPTEMP,15,1673
MPTEMP,16,1873
MPTEMP,17,2073
MPTEMP,18,2273
MPTEMP,19,2473
MPTEMP,20,2673
MPTEMP,21,2873
MPTEMP,22,3073
MPTEMP,23,3273
MPDATA,EX,1,,2.05E11
MPDATA,EX,1,,2.03E11
MPDATA,EX,1,,1.98E11
MPDATA,EX,1,,1.92E11
MPDATA,EX,1,,1.86E11
MPDATA,EX,1,,1.80E11
MPDATA,EX,1,,1.78E11
MPDATA,EX,1,,1.67E11
MPDATA,EX,1,,1.59E11
MPDATA,EX,1,,1.50E11
MPDATA,EX,1,,1.41E11
MPDATA,EX,1,,1.2E11
MPDATA,EX,1,,0.8E11
MPDATA,EX,1,,0.08E11
MPDATA,PRXY,1,,0.33
TB,BISO,1,6,2,
TBTEMP,298
TBDATA,,354e6,5.89e9,,,,
TBTEMP,873
TBDATA,,212e6,1.7e9,,,,
TBTEMP,1073
TBDATA,,199e6,1.4e9,,,,
TBTEMP,1473
TBDATA,,100e6,1e8,,,,
TBTEMP,1596
TBDATA,,50e6,1e7,,,,
TBTEMP,2173
TBDATA,,5e6,1e7,,,,
save
antype,trans
trnopt,full!
指定为完全瞬态分析
nropt,full,,on!
完全牛顿拉夫逊法自适应下降
nlgeom,1!
大变形
autots,on!
激活时间步自动阶跃
timint,on,struct!
打开结构瞬态分析
tintp,0.005,,,1,0.5,0.2
kbc,0
CNVTOL,F,,0.005,2,1.0e-6
CNVTOL,U,,0.005,2,1.0e-6
OUTRES,ERASE
OUTRES,NSOL,LAST
nsel,s,loc,y,0.015,0.05
d,all,all,0
allsel,all
nsel,s,loc,y,0
dsym,symm,y,,
*SET,ii,1
*DO,T,1E-6,1,TINC/2
TM1=TINC*（ii-1）/2+1E-6
ldread,temp,,,TM1,,Thermal3,rth
nsubst,20
solve
*SET,ii,ii+1
*enddo
*SET,aa,1
*DO,T,1,30.5,TINC/2
TM2=TINC*（aa-1）/2+1
ldread,temp,,,TM2,,Thermal3,rth
nsubst,6
*SET,aa,aa+1
*SET,k,1
*DO,T,30.51,55.51,0.1
TM3=（k-1）*0.1+30.51
TIME,TM3
ldread,temp,,,TM3,,Thermal3,rth
*SET,k,k+1
nsel,s,loc,y,0.03
d,all,ux,0
d,all,uz,0
*SET,bb,1
*DO,T,55.52,57.52,0.02
TM4=（bb-1）*0.02+55.52
TIME,TM4
ldread,temp,,,TM4,,Thermal3,rth
nsubst,30
*SET,bb,bb+1
*SET,jj,1
*DO,T,57.53,180,2
TM5=（jj-1）*2+57.53
TIME,TM5
ldread,temp,,,TM5,,Thermal3,rth
*SET,jj,jj+1
升级会员 