以下是程序流 按照一个书上的程序改的 但是结果一直不对
温度始终保持在20度
求大神指教!!!!
/FILNAME,weldingstress,0
/PREP7
!定义参数
W=1E-1
L=4E-2
H=6E-3
P=30
A=0.5
R=0.002
V=0.02
Qm=2*P*A/3.14/R/R
!*
ET,1,PLANE55
ET,2,SOLID70
MPTEMP,,,,,,,,
MPTEMP,1,0
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,600
MPTEMP,7,800
MPTEMP,8,1200
MPTEMP,9,1300
MPTEMP,10,1500
MPDATA,KXX,1,,14.6
MPDATA,KXX,1,,15.1
MPDATA,KXX,1,,16.1
MPDATA,KXX,1,,17.9
MPDATA,KXX,1,,18
MPDATA,KXX,1,,20.8
MPDATA,KXX,1,,23.9
MPDATA,KXX,1,,32.2
MPDATA,KXX,1,,33.7
MPDATA,KXX,1,,120
MPDATA,DENS,1,,7900
MPDATA,DENS,1,,7880
MPDATA,DENS,1,,7830
MPDATA,DENS,1,,7790
MPDATA,DENS,1,,7750
MPDATA,DENS,1,,7660
MPDATA,DENS,1,,7560
MPDATA,DENS,1,,7370
MPDATA,DENS,1,,7320
MPDATA,DENS,1,,7320
MPDATA,C,1,,462
MPDATA,C,1,,496
MPDATA,C,1,,512
MPDATA,C,1,,525
MPDATA,C,1,,540
MPDATA,C,1,,577
MPDATA,C,1,,604
MPDATA,C,1,,676
MPDATA,C,1,,692
MPDATA,C,1,,700
MPDATA,EX,1,,1.985E11
MPDATA,EX,1,,1.93E11
MPDATA,EX,1,,1.85E11
MPDATA,EX,1,,1.76E11
MPDATA,EX,1,,1.67E11
MPDATA,EX,1,,1.59E11
MPDATA,EX,1,,1.51E11
MPDATA,EX,1,,0.6e11
MPDATA,EX,1,,0.2e11
MPDATA,EX,1,,0.1e11
MPDATA,PRXY,1,,0.294
MPDATA,PRXY,1,,0.295
MPDATA,PRXY,1,,0.301
MPDATA,PRXY,1,,0.31
MPDATA,PRXY,1,,0.318
MPDATA,PRXY,1,,0.326
MPDATA,PRXY,1,,0.333
MPDATA,PRXY,1,,0.339
MPDATA,PRXY,1,,0.342
MPDATA,PRXY,1,,0.388
UIMP,1,REFT,,,20
MPDATA,ALPX,1,,1.7e-5
MPDATA,ALPX,1,,1.74e-5
MPDATA,ALPX,1,,1.80e-5
MPDATA,ALPX,1,,1.86e-5
MPDATA,ALPX,1,,1.91e-5
MPDATA,ALPX,1,,1.96e-5
MPDATA,ALPX,1,,2.02e-5
MPDATA,ALPX,1,,2.07e-5
MPDATA,ALPX,1,,2.11e-5
MPDATA,ALPX,1,,2.16e-5
TB,BISO,1,6,2,
TBTEMP,100
TBDATA,,218e6,0,,,,
TBTEMP,300
TBDATA,,170E6,0,,,,
TBTEMP,600
TBDATA,,149E6,0,,,,
TBTEMP,800
TBDATA,,91E6,0,,,,
TBTEMP,1200
TBDATA,,25E6,0,,,,
TBTEMP,1500
TBDATA,,10E6,0,,,,
K,1,0,0,0
K,2,W/2*0.5,0,0
K,3,W/2*0.7,0,0
K,4,W/2*0.85,0,0
K,5,W/2+W/2*0.15,0,0
K,6,W/2+W/2*0.3,0,0
K,7,W/2+W/2*0.5,0,0
K,8,W,0,0
K,9,W,L,0
K,10,W/2+W/2*0.5,L,0
K,11,W/2+W/2*0.3,L,0
K,12,W/2+W/2*0.15,L,0
K,13,W/2*0.85,L,0
K,14,W/2*0.7,L,0
K,15,W/2*0.5,L,0
K,16,0,L,0
K,17,0,0,H
A,1,2,15,16
A,2,3,14,15
A,3,4,13,14
A,4,5,12,13
A,5,6,11,12
A,6,7,10,11
A,7,8,9,10
ESIZE,0.0012
AMESH,4
ESIZE,0.0025
AMESH,3,5
ESIZE,0.005
AMESH,2,6
ESIZE,0.0065
AMESH,1,7
TYPE, 2
EXTOPT,ESIZE,2,0,
EXTOPT,ACLEAR,1
!*
EXTOPT,ATTR,1,0,0
REAL,_Z4
ESYS,0
!*
VOFFST,4,H, ,
VOFFST,3,H, ,
VOFFST,2,H, ,
VOFFST,1,H, ,
VOFFST,5,H, ,
VOFFST,6,H, ,
VOFFST,7,H, ,
EPLOT
NUMMRG,ALL, , , ,LOW
/SOL
!*
ANTYPE,4
!*
TRNOPT,FULL
LUMPM,0
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNCCSYS
*SET,_FNCNAME,'GAOSI'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*SET,_FNC_C1(1),Qm
*SET,_FNC_C2(1),V
*SET,_FNC_C3(1),R
*SET,_FNCCSYS,0
! /INPUT,gaosi.func,,,1
*DIM,%_FNCNAME%,TABLE,6,19,1,,,,%_FNCCSYS%
!
! Begin of equation: Qm*exp(-2*({X}^2+({Y}-V*{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), %_FNC_C3(1)%
*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), 0, -3, 0, 1, -1, 2, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 2, 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, 2, 0, 0, 2
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, 2, 17, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 1, 18, 3, 1
*SET,%_FNCNAME%(0,9,1), 0.0, -4, 0, 1, 3, 2, -1
*SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 2, 0, 0, -4
*SET,%_FNCNAME%(0,11,1), 0.0, -5, 0, 1, -4, 17, -1
*SET,%_FNCNAME%(0,12,1), 0.0, -1, 0, 1, -3, 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, 19
*SET,%_FNCNAME%(0,15,1), 0.0, -2, 0, 1, 19, 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, 1, 17, 3, -1
*SET,%_FNCNAME%(0,19,1), 0.0, 99, 0, 1, -2, 0, 0
! End of equation: Qm*exp(-2*({X}^2+({Y}-V*{TIME})^2)/R^2)
TUNIF,20, !定义初始温度
!定义对流换热边界
SFA,23,1,CONV,30,20
SFA,18,1,CONV,30,20
SFA,33,1,CONV,30,20
SFA,38,1,CONV,30,20
SFA,40,1,CONV,30,20
SFA,24,1,CONV,30,20
SFA,19,1,CONV,30,20
SFA,14,1,CONV,30,20
SFA,9,1,CONV,30,20
SFA,29,1,CONV,30,20
SFA,34,1,CONV,30,20
SFA,39,1,CONV,30,20
SFA,41,1,CONV,30,20
SFA,36,1,CONV,30,20
SFA,31,1,CONV,30,20
SFA,11,1,CONV,30,20
SFA,16,1,CONV,30,20
SFA,21,1,CONV,30,20
SFA,26,1,CONV,30,20
SFA,27,1,CONV,30,20
SFA,1,1,CONV,30,20
SFA,2,1,CONV,30,20
SFA,3,1,CONV,30,20
SFA,4,1,CONV,30,20
SFA,5,1,CONV,30,20
SFA,6,1,CONV,30,20
SFA,7,1,CONV,30,20
SFA,13,1,CONV,30,20
SFA,28,1,CONV,30,20
!施加高斯热源
SFA,8,1,HFLUX, %GAOSI%
OUTRES,ALL,ALL,
TIME,L/V !设置求解时间
AUTOTS,-1
NSUBST,50,50,50
KBC,0
TSRES,ERASE
LSWRITE,1, !写入载荷文件为1
!*
TIME,4
AUTOTS,1
NSUBST,20,20,20
KBC,0
!*
TSRES,ERASE
LSWRITE,2,
TIME,10
AUTOTS,1
NSUBST,20,20,20
KBC,0
!*
TSRES,ERASE
LSWRITE,3,
TIME,50
AUTOTS,1
NSUBST,50,50,50
KBC,0
!*
TSRES,ERASE
LSWRITE,4,
LSSOLVE,1,4,1, !开始求解
|