为什么我算出来的焊接温度场比实际的偏低？

lijm114

/PREP7
/UNITS,SI
ET,1,SOLID70
MPTEMP,,,,,,,,

MPTEMP,1,20
MPTEMP,2,250
MPTEMP,3,500
MPTEMP,4,750
MPTEMP,5,1000

MPTEMP,6,1500

MPTEMP,7,1700

MPTEMP,8,2500

MPDATA,EX,1,,2.05E+011

MPDATA,EX,1,,1.87E+011

MPDATA,EX,1,,1.50E+011

MPDATA,EX,1,,7.00E+010

MPDATA,EX,1,,2.00E+010

MPDATA,EX,1,,1.90E+010

MPDATA,EX,1,,1.80E+010

MPDATA,EX,1,,1.20E+010

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

MPTEMP,,,,,,,,

MPTEMP,1,20
MPTEMP,2,250
MPTEMP,3,500
MPTEMP,4,750
MPTEMP,5,1000

MPTEMP,6,1500

MPTEMP,7,1700

MPTEMP,8,2500

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
MPTEMP,,,,,,,,

MPTEMP,1,20
MPTEMP,2,250
MPTEMP,3,500
MPTEMP,4,750
MPTEMP,5,1000

MPTEMP,6,1500

MPTEMP,7,1700

MPTEMP,8,2500

UIMP,1,REFT,,,

MPDATA,ALPX,1,,1.10E-005
MPDATA,ALPX,1,,1.22E-005
MPDATA,ALPX,1,,1.39E-005
MPDATA,ALPX,1,,1.48E-005
MPDATA,ALPX,1,,1.34E-005
MPDATA,ALPX,1,,1.33E-005
MPDATA,ALPX,1,,1.32E-005
MPDATA,ALPX,1,,1.31E-005
MPTEMP,,,,,,,,

MPTEMP,1,20
MPTEMP,2,250
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,,140

MPDATA,KXX,1,,142

MPTEMP,,,,,,,,

MPTEMP,1,20
MPTEMP,2,250
MPTEMP,3,500
MPTEMP,4,750
MPTEMP,5,1000

MPTEMP,6,1500

MPTEMP,7,1700

MPTEMP,8,2500

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
TB,BISO,1,6,2,

TBTEMP,20

TBDATA,,0.220E+009,,,,,
TBTEMP,250

TBDATA,,0.175E+009,,,,,
TBTEMP,500

TBDATA,,0.130E+009,,,,,
TBTEMP,750

TBDATA,,0.040E+009,,,,,
TBTEMP,1000
TBDATA,,0.025E+009,,,,,
TBTEMP,1500
TBDATA,,0.002E+009,,,,,
MPTEMP,,,,,,,,

MPTEMP,1,20
MPTEMP,2,250
MPTEMP,3,500
MPTEMP,4,750
MPTEMP,5,1000

MPTEMP,6,1500

MPTEMP,7,1700

MPTEMP,8,2500

MPDATA,HF,1,,100
MPDATA,HF,1,,350
MPDATA,HF,1,,520
MPDATA,HF,1,,1000

MPDATA,HF,1,,1500

MPDATA,HF,1,,3000

MPDATA,HF,1,,3100

MPDATA,HF,1,,3500
MPTEMP,,,,,,,,
MPTEMP,1,20
MPTEMP,2,1533
MPTEMP,3,1590
MPTEMP,4,1670
MPDE,ENTH,1
MPDATA,ENTH,1,,0
MPDATA,ENTH,1,,7.5E+009
MPDATA,ENTH,1,,9.6E+009
MPDATA,ENTH,1,,1.1E+010
*SET,RB,0.002

*SET,LENGTH,0.12
*SET,WIDTH,0.06
*SET,HEIGHT,0.006

BLOCK,0,LENGTH,0,WIDTH,0,HEIGHT,
BLOCK,0,LENGTH,0,3*RB,0,HEIGHT,
BLOCK,0,LENGTH,0,6*RB,0,HEIGHT,
BLOCK,0,LENGTH,0,12*RB,0,HEIGHT,
FLST,2,4,6,ORDE,2

FITEM,2,1

FITEM,2,-4

VPTN,P51X

SAVE
/PNUM,KP,0

/PNUM,LINE,1
/PNUM,AREA,0
/PNUM,VOLU,1
/PNUM,NODE,0
/PNUM,TABN,0
/PNUM,SVAL,0
/NUMBER,0

!*

/PNUM,ELEM,0
/REPLOT
!*

FLST,5,12,4,ORDE,10
FITEM,5,4

FITEM,5,-5

FITEM,5,9

FITEM,5,-10
FITEM,5,14

FITEM,5,19

FITEM,5,23

FITEM,5,-24
FITEM,5,49

FITEM,5,-52
CM,_Y,LINE

LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y

!*

LESIZE,_Y1,RB, , , , , , ,1
!*

FLST,5,2,4,ORDE,2

FITEM,5,26

FITEM,5,31

CM,_Y,LINE

LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y

!*

LESIZE,_Y1,2*RB, , , , , , ,1

!*

FLST,5,8,4,ORDE,6

FITEM,5,2

FITEM,5,7

FITEM,5,38

FITEM,5,43

FITEM,5,53

FITEM,5,-56
CM,_Y,LINE

LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y

!*

LESIZE,_Y1,4*RB, , , , , , ,1

!*

FLST,5,4,4,ORDE,4

FITEM,5,11

FITEM,5,-12
FITEM,5,47

FITEM,5,-48
CM,_Y,LINE

LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y

!*

LESIZE,_Y1,1.5*RB, , , , , , ,1
!*

MSHAPE,0,3D
MSHKEY,1
!*

FLST,5,2,6,ORDE,2

FITEM,5,5

FITEM,5,-6

CM,_Y,VOLU

VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'

CMSEL,S,_Y

!*

VMESH,_Y1

!*

CMDELE,_Y

CMDELE,_Y1

CMDELE,_Y2

!*

MSHAPE,1,3D
MSHKEY,0
!*

FLST,5,2,6,ORDE,2

FITEM,5,7

FITEM,5,-8

CM,_Y,VOLU

VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'

CMSEL,S,_Y

!*

VMESH,_Y1

!*

CMDELE,_Y

CMDELE,_Y1

CMDELE,_Y2

!*

/DIST,1,1.08222638492,1
/REP,FAST

/DIST,1,0.924021086472,1
/REP,FAST

SAVE
/solu
antype,trans
trnopt,full
nropt,auto, ,
autots,on
kbc,1
tintp,0.005,,,1,0.5,0.2
tref,25
outres,all
*set,v,0.01
*set,tinc,rb/v
*set,U,25
*set,I,180
*set,ef,0.75
*set,p,25*180*ef
*set,r,0.006
*set,qmax,p/(3.14159*r**2)
*set,max_time,length/rb
*set,max_x,length/rb+1
*set,max_y,width/rb+1
*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)*0.002
*enddo
*do,j,1,max_y,1
flux2(0,j,k)=(j-1)*0.002
*enddo
*do,i,1,max_x,1
*do,j,1,max_y,1
xcenter=v*(k-1)*tinc
ycenter=0
distance=sqrt(((i-1)*0.002-xcenter)**2+((j-1)*0.002-ycenter)**2)
*if,distance,le,r,then
flux2(i,j,k)=qmax/EXP(3*distance**2/(r**2))
*else
flux2(i,j,k)=0
*endif
*enddo
*enddo
flux2(0,0,k)=(k-1)*tinc
nsel,s,loc,z,height
sf,all,hflux,%flux2%
allsel,all
nsel,s,loc,x,0
sf,all,conv,10,20
allsel,all
nsel,s,loc,x,length
sf,all,conv,10,20
allsel,all
nsel,s,loc,y,width
sf,all,conv,10,20
allsel,all
nsel,s,loc,z,0
sf,all,conv,10,20
allsel,all
time,k*0.002/v
nsubst,5
solve
*enddo
save