求助大神 我做了一个温度场的模拟 但是结果一直都是20度

757817641

以下是程序流  按照一个书上的程序改的 但是结果一直不对
温度始终保持在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,  !开始求解

小七工作室

边界加载有问题哦