为什么我算出来的焊接温度场比实际的偏低？
/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
页:
[1]