本人新人,在计算过程中遇到一个问题,求各位前辈帮忙!激光焊模拟,在计算过程中出现 value=-833.33 is outside function range EXP,请问什么意思,该怎么解决,谢谢!
命令流如下:
/config,nres,60000
/prep7
ET,1,SOLID70
/prep7
mp,dens,1,4440
mptemp,1,20,100,200,300,400,500,
mptemp,7,600,700,800,1000,1200,1500,,,,
mpdata,kxx,1,1,6.8,7.4,8.7,9.8,10.3,11.8,
mpdata,kxx,1,7,13.3,14.8,16.3,19.3,22.3,26.8,,,,
mpdata,c,1,1,611,624,653,674,691,703,
mpdata,c,1,7,711,719.6,726.8,739.1,749.3,762,,,,
MPTEMP,,,,,,,,
MPTEMP,1,20
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,600
MPTEMP,8,700
MPTEMP,9,800
MPDE,EX,1
MPDE,EY,1
MPDE,EZ,1
MPDE,NUXY,1
MPDE,NUYZ,1
MPDE,NUXZ,1
MPDE,PRXY,1
MPDE,PRYZ,1
MPDE,PRXZ,1
MPDE,GXY,1
MPDE,GYZ,1
MPDE,GXZ,1
MPDATA,EX,1,,109E9
MPDATA,EX,1,,102E9
MPDATA,EX,1,,95E9
MPDATA,EX,1,,87E9
MPDATA,EX,1,,80E9
MPDATA,EX,1,,72E9
MPDATA,EX,1,,65E9
MPDATA,EX,1,,57E9
MPDATA,EX,1,,49E9
MPDATA,PRXY,1,,0.34
MPDATA,PRXY,1,,0.34
MPDATA,PRXY,1,,0.34
MPDATA,PRXY,1,,0.34
MPDATA,PRXY,1,,0.34
MPDATA,PRXY,1,,0.34
MPDATA,PRXY,1,,0.34
MPDATA,PRXY,1,,0.34
MPDATA,PRXY,1,,0.34
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!model
BLOCK,0,0.001,0,0.12,-0.0015,0.0015
BLOCK,0.001,0.007,0,0.12,-0.0015,0.0015
block,0.007,0.015,0,0.12,-0.0015,0.0015
VGLUE,all
!***mesh weld line(2mm)
!y-direction along the weld direction
LESIZE,1, , ,480, , , , ,1
LESIZE,3, , ,480, , , , ,1
LESIZE,6, , ,480, , , , ,1
LESIZE,8, , ,480, , , , ,1
!x-direction
LESIZE,2, , ,4, , , , ,1
LESIZE,4, , ,4,, , , ,1
LESIZE,5, , ,4, , , , ,1
LESIZE,7, , ,4, , , , ,1
!z-direction
LESIZE,10, , ,10, , , , ,1
LESIZE,9, , ,10, , , , ,1
LESIZE,11, , ,10, , , , ,1
LESIZE,12, , ,10, , , , ,1
!heat affect zone
!x-direction
LESIZE,38, , ,12, , , , ,1
LESIZE,39, , ,12, , , , ,1
LESIZE,37, , ,12, , , , ,1
LESIZE,40, , ,12, , , , ,1
!y-direction
LESIZE,3, , ,480, , , , ,1
LESIZE,6, , ,480, , , , ,1
LESIZE,15, , ,480, , , , ,1
LESIZE,18, , ,480, , , , ,1
!z-direction
LESIZE,10, , ,10, , , , ,1
LESIZE,11, , ,10, , , , ,1
LESIZE,22, , ,10, , , , ,1
LESIZE,23, , ,10, , , , ,1
! the third part
!x-direction
LESIZE,42, , ,4, , , , ,1
LESIZE,43, , ,4, , , , ,1
LESIZE,41, , ,4, , , , ,1
LESIZE,44, , ,4, , , , ,1
!y-direction
LESIZE,15, , ,480,, , , ,1
LESIZE,18, , ,480,, , , ,1
LESIZE,30, , ,480,, , , ,1
LESIZE,27, , ,480,, , , ,1
!z-direction
LESIZE,34, , ,10,, , , ,1
LESIZE,35, , ,10,, , , ,1
!mesh
MSHAPE,0,3D
MSHKEY,1
vmesh,1
vmesh,4
vmesh,5
!solve
!sol
*DIM,conv,TABLE,9,1,1,temp, ,
!*
*SET,CONV(0,1,1) , 0
*SET,CONV(1,0,1) , 20
*SET,CONV(1,1,1) , 50
*SET,CONV(2,0,1) , 300
*SET,CONV(2,1,1) , 50
*SET,CONV(3,0,1) , 600
*SET,CONV(3,1,1) , 50
*SET,CONV(4,0,1) , 750
*SET,CONV(4,1,1) , 50
*SET,CONV(5,0,1) , 900
*SET,CONV(5,1,1) , 50
*SET,CONV(6,0,1) , 1200
*SET,CONV(6,1,1) , 50
*SET,CONV(7,0,1) , 1500
*SET,CONV(7,1,1) , 50
*SET,CONV(8,0,1) , 1800
*SET,CONV(8,1,1) , 50
*SET,CONV(9,0,1) , 2000
*SET,CONV(9,1,1) , 50
!***********
*DIM,conv1,TABLE,9,1,1,temp, ,
!*
*SET,CONV1(0,1,1) , 0
*SET,CONV1(1,0,1) , 20
*SET,CONV1(1,1,1) , 50
*SET,CONV1(2,0,1) , 300
*SET,CONV1(2,1,1) , 50
*SET,CONV1(3,0,1) , 600
*SET,CONV1(3,1,1) , 50
*SET,CONV1(4,0,1) , 750
*SET,CONV1(4,1,1) , 50
*SET,CONV1(5,0,1) , 900
*SET,CONV1(5,1,1) , 50
*SET,CONV1(6,0,1) , 1200
*SET,CONV1(6,1,1) , 50
*SET,CONV1(7,0,1) , 1500
*SET,CONV1(7,1,1) , 50
*SET,CONV1(8,0,1) , 1800
*SET,CONV1(8,1,1) , 50
*SET,CONV1(9,0,1) , 2000
*SET,CONV1(9,1,1) , 50
!*************
a=10e-3
a1=5e-3
b=0.3e-3
c=4.5e-3
c1=4e-3
e=2.72
p=4000
eff=0.85
pi=3.14
h=3e-3
r0=0.3e-3
f1=0.4
f2=0.25
f3=0.2
!*******************************
/sol
ANTYPE,4
TUNIF,20,
TREF,20,
IC,all,TEMP,20,
NSEL,s,LOC,z,0.0015
SF,all,CONV,%CONV% ,20
allsel,all
NSEL,s,LOC,z,-0.0015
SF,all,CONV,%CONV% ,20
allsel,all
j=0
*do,i,0.01,0.11,0.00025
j=j+0.015
local,11,1,0,i,
nsel,s,loc,x,0,0.005
nsel,r,loc,z,-0.0015,0.0015
nplot
csys,0
*get,nmax,node,,num,max
*get,nmin,node,,num,min
*do,ei,nmin,nmax
*if,nsel(ei),eq,1,then
*get,elx,node,ei,loc,x
*get,ely,node,ei,loc,y
*get,elz,node,ei,loc,z
*if,elz,ge,-0.0015,and,elz,le,-0.0012,then
ely1=ely-i
elz1=elz+0.0015
!_____________________
qa1=6*sqrt(3.0)*p
qa2=qa1/(a*b*c*pi*sqrt(pi))
qa3=exp(-3*elx**2/c**2)
qa4=qa3*exp(-3*ely1**2/a**2)
qa5=qa4*exp(-3*elz1**2/b**2)
qm1=eff*qa2*qa5*f1
q10=3.15*p
q20=q10/(pi*h*r0**2)
q30=exp(-3*elx**2/r0**2)
q40=q30*exp(-3*ely1**2/r0**2)
qm2=eff*q20*q40*f2
qm3=qm1+qm2
BF,ei,HGEN,qm3
*elseif,elz,ge,-0.0012,and,elz,le,0.0012,then
ely1=ely-i
q10=3.15*p
q20=q10/(pi*h*r0**2)
q30=exp(-3*elx**2/r0**2)
q40=q30*exp(-3*ely1**2/r0**2)
qm2=eff*q20*q40*f2
qm3=qm2
BF,ei,HGEN,qm3
*elseif,elz,ge,0.0012,and,elz,le,0.0015,then
ely1=ely-i
elz1=elz-0.0015
!_____________________
qa1=6*sqrt(3.0)*p
qa2=qa1/(a1*b*c1*pi*sqrt(pi))
qa3=exp(-3*elx**2/c1**2)
qa4=qa3*exp(-3*ely1**2/a1**2)
qa5=qa4*exp(-3*elz1**2/b**2)
qm1=eff*qa2*qa5*f3
q10=3.15*p
q20=q10/(pi*h*r0**2)
q30=exp(-3*elx**2/r0**2)
q40=q30*exp(-3*ely1**2/r0**2)
qm2=eff*q20*q40*f2
qm3=qm1+qm2
BF,ei,HGEN,qm3
*endif
*endif
*enddo
TIME,j
LNSRCH,1
TIMINT,1
!TINTP,0.007, , ,1, , ,
AUTOTS,1
DELTIM,0.03
outres,nsol,all
KBC,1
!*
TSRES,ERASE
allsel,all
solve
!*********
BFDELE,ALL,HGEN
BF,ALL,HGEN,0
allsel,all
!*********
*ENDDO
!!cooling
allsel,all
BFDELE,ALL,HGEN
BF,ALL,HGEN,0
SFDELE,all,CONV
NSEL,s,LOC,z,0.0015
SF,all,CONV,%CONV1% ,20
allsel,all
NSEL,s,LOC,z,-0.0015
sf,all,conv,%CONV1% ,20
allsel,all
NSEL,S,LOC,y,0
sf,all,conv,%CONV1% ,20
allsel,all
NSEL,S,LOC,y,0.12
sf,all,conv,%CONV1% ,20
allsel,all
NSEL,s,LOC,x,0.015
SF,all,CONV,%CONV1% ,20
allsel,all
time,9
deltim,0.2
!outres,all,last
OUTRES,NSOL,all,
autots,on
!kbc,0
kbc,1
solve
time,200
deltim,2
!outres,all,last
OUTRES,NSOL,all,
autots,on
kbc,1
solve
time,300
deltim,5
!outres,all,last
OUTRES,NSOL,all,
autots,on
kbc,1
solve
time,500
deltim,10
!outres,all,last
OUTRES,NSOL,all,
autots,on
kbc,1
solve
save
finish |