/FILNAME,weldingstress,0
/PREP7
!定义焊接参数
L=1E-1 !焊件的长度
H=6E-3 !焊件的高度
W1=360 !整个圆360度
w=5 !焊接角速度
R=0.007 !电弧有效加热半径
Qm=4e7 !加热斑点中心最大热流密度
!*
ET,1,PLANE55
ET,2,SOLID70
MPTEMP,,,,,,,,
MPTEMP,1,20
MPTEMP,2,200
MPTEMP,3,500
MPTEMP,4,750
MPTEMP,5,1000
MPTEMP,6,1500
MPTEMP,7,1700
MPTEMP,8,2500
MPDATA,KXX,1,,50
MPDATA,KXX,1,,47
MPDATA,KXX,1,,40
MPDATA,KXX,1,,27
MPDATA,KXX,1,,30
MPDATA,KXX,1,,35
MPDATA,KXX,1,,40
MPDATA,KXX,1,,55
MPDATA,DENS,1,,7820
MPDATA,DENS,1,,7700
MPDATA,DENS,1,,7610
MPDATA,DENS,1,,7550
MPDATA,DENS,1,,7490
MPDATA,DENS,1,,7350
MPDATA,DENS,1,,7300
MPDATA,DENS,1,,7090
MPDATA,C,1,,460
MPDATA,C,1,,480
MPDATA,C,1,,530
MPDATA,C,1,,675
MPDATA,C,1,,670
MPDATA,C,1,,660
MPDATA,C,1,,780
MPDATA,C,1,,820
MPDATA,EX,1,,2.05E11
MPDATA,EX,1,,1.87E11
MPDATA,EX,1,,1.5E11
MPDATA,EX,1,,0.7E11
MPDATA,EX,1,,0.2E11
MPDATA,EX,1,,0.19E2
MPDATA,EX,1,,0.18E2
MPDATA,EX,1,,0.12e2
MPDATA,PRXY,1,,0.28
MPDATA,PRXY,1,,0.29
MPDATA,PRXY,1,,0.31
MPDATA,PRXY,1,,0.35
MPDATA,PRXY,1,,0.4
MPDATA,PRXY,1,,0.45
MPDATA,PRXY,1,,0.48
MPDATA,PRXY,1,,0.5
UIMP,1,REFT,,,20
MPDATA,ALPX,1,,1.1e-5
MPDATA,ALPX,1,,1.22e-5
MPDATA,ALPX,1,,1.39e-5
MPDATA,ALPX,1,,1.48e-5
MPDATA,ALPX,1,,1.34e-5
MPDATA,ALPX,1,,1.33e-5
MPDATA,ALPX,1,,1.32e-5
MPDATA,ALPX,1,,1.31e-5
TB,BISO,1,6,2,
TBTEMP,20
TBDATA,,220e6,0,,,,
TBTEMP,250
TBDATA,,175e6,0,,,,
TBTEMP,500
TBDATA,,80e6,0,,,,
TBTEMP,750
TBDATA,,40E6,0,,,,
TBTEMP,1000
TBDATA,,10E6,0,,,,
TBTEMP,1500
TBDATA,,1E-5,0,,,,
K,1,L,L,0
K,2,L,-L,0
K,3,-L,-L,0
K,4,-L,L,0
A,1,2,3,4
ESIZE,0.005
AMESH,1
TYPE, 2
EXTOPT,ESIZE,2,0,
EXTOPT,ACLEAR,1
!*
EXTOPT,ATTR,1,0,0
REAL,_Z4
ESYS,0
!*
VOFFST,1,H, ,
NUMMRG,ALL, , , ,LOW
/SOL
!*
ANTYPE,4
!*
TRNOPT,FULL
LUMPM,0
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNCCSYS
*SET,_FNCNAME,'zhu1yzzb'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*SET,_FNC_C1(1),5
*SET,_FNC_C2(1),0.007
*SET,_FNCCSYS,1
! /INPUT,H:\Desktop\csys1zhugui.func,,,1
*DIM,%_FNCNAME%,TABLE,6,20,1,,,,%_FNCCSYS%
! Begin of equation: 4e7*exp(-3*((0.5)^2+(w*{time})^2)/r^2)
*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), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,3,1), 1, -3, 0, 1, -1, 2, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 3, 0, 0, -3
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 0.5, 0, 0, 0
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 2, 0, 0, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -1, 17, -3
*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 1, 17, 3, 1
*SET,%_FNCNAME%(0,10,1), 0.0, -3, 0, 2, 0, 0, -1
*SET,%_FNCNAME%(0,11,1), 0.0, -5, 0, 1, -1, 17, -3
*SET,%_FNCNAME%(0,12,1), 0.0, -1, 0, 1, -4, 1, -5
*SET,%_FNCNAME%(0,13,1), 0.0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,14,1), 0.0, -1, 0, 2, 0, 0, 18
*SET,%_FNCNAME%(0,15,1), 0.0, -2, 0, 1, 18, 17, -1
*SET,%_FNCNAME%(0,16,1), 0.0, -1, 0, 1, -3, 4, -2
*SET,%_FNCNAME%(0,17,1), 0.0, -1, 7, 1, -1, 0, 0
*SET,%_FNCNAME%(0,18,1), 0.0, -2, 0, 4e7, 0, 0, -1
*SET,%_FNCNAME%(0,19,1), 0.0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,20,1), 0.0, 99, 0, 1, -3, 0, 0
! End of equation: 4e7*exp(-3*((0.5)^2+(w*{time})^2)/r^2)
TUNIF,20, !定义初始温度
SFA,2,1,CONV,30,20
SFA,3,1,CONV,30,20
SFA,3,1,CONV,30,20
SFA,4,1,CONV,30,20
SFA,5,1,CONV,30,20
!施加高斯热源
SFA,1,1,HFLUX, %zhu1yzzb%
OUTRES,ALL,ALL,
TIME,W1/W !设置求解时间
AUTOTS,-1
NSUBST,50,50,50
KBC,0
TSRES,ERASE
LSWRITE,1, !写入载荷文件为1
!*
TIME,80
AUTOTS,1
NSUBST,20,20,20
KBC,0
!*
TSRES,ERASE
LSWRITE,2,
TIME,100
AUTOTS,1
NSUBST,30,30,30
KBC,0
!*
TSRES,ERASE
LSWRITE,3,
TIME,1100
AUTOTS,1
NSUBST,105,105,105
KBC,0
!*
TSRES,ERASE
LSWRITE,4,
LSSOLVE,1,4,1, !开始求解
我是用柱坐标csys1,w为角速度,w1是360圆 |