fini
/cle
!/batch,list
/prep7
!title,weld analysisusing"birth and death"
et,1,13,4
et,2,13,4
et,3,5,
!UX, UY, TEMP, AZ if KEYOPT (1) = 4,plane13
!dimension mm,n,g,Mpa
!热传导系数
mp,kxx,1,70
mp,kxx,2,70
!比热
mp,c,1,460
mp,c,2,460
!密度
mp,dens,1,7850
mp,dens,2,7850
!热膨胀系数
mp,alpx,1,1.2e-5
mp,alpx,2,1.2e-5
!泊松比
mp,nuxy,1,0.30
mp,nuxy,2,0.30
!磁特性
mp,murx,1,1
mp,murx,2,1
!t为隔板厚d为顶板厚
t=0.03
d=0.06
k,1,-(5*d+t/2),0
k,2,-(5*d+t/2),d
k,3,-(t/2+t/3),d
k,4,-(t/2+t/3),d-t/3
k,5,-t/6,d
k,6,-t/2,d+t/3
k,7,-t/2,d+5*t
k,8,-(t/2+t/3),0
k,9,0,0
k,10,0,d
k,11,0,d+5*t
a,3,4,5,6
a,1,2,3,4,8
a,8,9,10,5,4
a,6,5,10,11,7
!1为母材2为焊缝
type,2
mat,2
asel,s,,,1
lesize,4,,,6
lesize,2,,,8
lesize,3,,,6
lesize,1,,,4
amesh,all,,
type,1
mat,1
asel,s,,,2,4
lesize,6,,,20
lesize,8,,,20
lesize,7,,,6
lesize,10,,,8
lesize,5,,,4
lesize,9,,,4
lesize,11,,,2
lesize,13,,,2
lesize,14,,,18
lesize,12,,,20
amesh,all,,
asel,all
arsym,x,all
nummrg,kp,,,
nummrg,node,,,
nummrg,kp,,,
nummrg,node,,,
save,birth,db
resum,birth,db
mptemp,1,0,100,600,800,1000,1500
mpdata,ex,1,1,2.1e11,2.0e11,1.5e11,10e10,5e10,2e9
mpdata,ex,2,1,2.1e11,2.0e11,1.5e11,10e10,5e10,2e9
tb,bkin,1,6
tbtemp,0
tbdata,1,370e6,1e11
tbtemp,100
tbdata,1,370e6,1e11
tbtemp,600
tbdata,1,200e6,5e10
tbtemp,800
tbdata,1,120e6,1e10
tbtemp,1000
tbdata,1,70e6,5e9
tbtemp,1500
tbdata,1,50e6,1e8
tbcopy,bkin,1,2
ASEL,S,MAT,,2
TYPE,3
MAT,2
EXTOPT,ESIZE,10,0,
EXTOPT,ACLEAR,1
vext,all,,,,,-0.1
ASEL,S,MAT,,1
TYPE,3
MAT,1
EXTOPT,ESIZE,10,0,
EXTOPT,ACLEAR,1
vext,all,,,,,-0.1
allsel
nummrg,kp,,,
nummrg,node,,,
nummrg,kp,,,
nummrg,node,,,
VSEL,S,MAT,,2
eslv,s
nsle,,all
cm,wnode,node
cm,welem,elem
nsel,s,loc,y,d+5*t
d,all,ux,0
d,all,uy,0
d,all,uz,0
nsel,all
esel,all
asel,s,,,20,
asel,a,,,40
lsla,s,
nsll,s,1
sf,all,conv,13.9,5
allsel
finish
/solu
antype,trans
trnopt,full
nropt,auto,,
eqslv, !select the eqation solver
solcontrol,on
autots,on !specified whether to use automatic load step or time step.
kbc,1 !specified ramped or stepped loading within load step
outres,all,all, ! control the solution data into database
!load step 1, initial conditions 25 c
time,1
deltim,0.2,0.2,0.2
timint,on !time integration = off
tunif,25,
solve
!load step 2..... , apply moving heat flux
timint,on !time integration = on
j=2
allsel,all
*do,i,0,4
time,j
deltim,0.11,0.11,0.11
sfdele,all,hflux
cmsel,s,wnode,node
nsel,r,loc,z,-0.02*i-0.02,-0.02*i
sf,all,hflux,167.44e5
allsel,all
nsubst,1
outres,all,all
solve
j=j+1
*enddo
time,7
deltim,0.11,0.11,0.11
allsel,all
sfdele,all,hflux
allsel,all
nsubst,1
outres,all,all
solve
time,8
deltim,0.11,0.11,0.11
allsel,all
ddele,all,temp
nsubst,1
outres,all,all
solve
finish
为什么我运行到time8的温度场结果并没有全部删去,还有温度呢?另外求教sf,all,hflux,%,%可不可以用热流密度除以节点数求得?谢谢! |