/filn,a11
/units,si
/prep7
et,1,solid5
et,2,solid5
mptemp,1,20,500,1000,1500,2000
mpdata,ex,1,1,2.06e11,1.7e11,1e11,0.5e11,0.1e11
mpdata,dens,1,1,7800,7800,7800,7800,7800
mpdata,nuxy,1,1,0.3,0.3,0.3,0.3,0.3
mpdata,kxx,1,1,16.3,16.3,16.3,16.3,16.3
mpdata,alpx,1,1,1.06e-5,1.06e-5,1.06e-5,1.06e-5,1.06e-5
mpdata,c,1,1,502,502,502,502,502
tb,bkin,1,5
tbtemp,20,1
tbdata,1,1.4e9,2.06e10
tbtemp,500,2
tbdata,1,1e9,1.7e10
tbtemp,1000,3
tbdata,1,0.5e9,0.9e10
tbtemp,1500,4
tbdata,1,0.07e9,0.2e10
tbtemp,2000,5
tbdata,1,0.007e9,0.01e10
mpdata,ex,2,1,2.07e11,1.8e11,0.9e11,0.3e11,0.1e11
mpdata,dens,2,1,7811,7811,7811,7811,7811
mpdata,nuxy,2,1,0.3,0.3,0.3,0.3,0.3
mpdata,kxx,2,1,14.9,14.9,14.9,14.9,14.9
mpdata,alpx,2,1,1.08e-5,1.08e-5,1.08e-5,1.08e-5,1.08e-5
mpdata,c,2,1,522,522,522,522,522
tb,bkin,2,5
tbtemp,20,1
tbdata,1,1.5e9,2.07e10
tbtemp,500,2
tbdata,1,1.1e9,1.8e10
tbtemp,1000,3
tbdata,1,0.6e9,1e10
tbtemp,1500,4
tbdata,1,0.08e9,0.3e10
tbtemp,2000,5
tbdata,1,0.008e9,0.02e10
k,1,0,0,0
k,2,0.15
k,3,0.137,0.025
k,4,0,0.025
k,5,0.152
k,6,0.3
k,7,0.3,0.025
k,8,0.152,0.025
k,9,0.145,0.013
a,1,2,3,4
a,5,6,7,8,9
k,14,0.145,0.0225
larc,10,11,14
l,10,13
l,12,13
l,12,9
l,11,9
k,20,0,0,0.01
l,1,20
vglue,all
nummrg,all
numcmg,all
allsel
/image,save,vmesh,bmp
/solu
antype,4
nsel,s,loc,x,0
d,all,ux,0
nsel,r,loc,y,0
d,all,uy,0
nsel,r,loc,z,0
d,all,uz,0
nsel,all
nsel,s,loc,x,0
d,all,temp,20
nsel,r,loc,x,0.3
d,all,temp,20
esel,all
*get,emax,elem,,num,max
vsel,s,volu,,3
eslv
*get,nse,elem,,count
*dim,ne,,nse
*dim,nex,,nse
*dim,ney,,nse
*dim,nez,,nse
*dim,neorder,,nse
mine=0
*do,i1,1,nse
esel,u,elem,,mine
*get,nse1,elem,,count
ii=0
*do,i,1,emax
*if,esel(i),eq,1,then
ii=ii+1
ne(ii)=i
*endif
*enddo
*do,i,1,nse1
*get,nez(i),elem,ne(i),cent,z
*get,ney(i),elem,ne(i),cent,y
*get,nex(i),elem,ne(i),cent,x
*enddo
minz=1e20
miny=1e20
minx=1e20
*do,i,1,nse1
*if,nez(i),lt,minz,then
minz=nez(i)
miny=ney(i)
minx=nex(i)
mine=ne(i)
*else
*if,nez(i),eq,minz,then
*if,ney(i),lt,miny,then
minz=nez(i)
miny=ney(i)
minx=nex(i)
mine=ne(i)
*else
*if,ney(i),eq,miny,then
*if,nex(i),lt,minx,then
minz=nez(i)
miny=mey(i)
minx=nex(i)
mine=ne(i)
*endif
*endif
*endif
*endif
*endif
*enddo
neorder(i1)=mine
*enddo
max_tem=1600
dt1=1e-8
dt=0.1111
t=0
esel,all
eplot
*do,i,1,nse
ekill,neorder(i)
esel,s,live
*enddo
allsel,all
outres,all,all
ic,all,temp,20
kbc,1
timint,0,struct
timint,1,therm
timint,0,mag
tintp,0.005,,,1,0.5,0.2
nsub1=2
nsub2=40
*do,i,1,nse
ealive,neorder(i)
esel,s,live
esel,all
t=t+dt1
time,t
nsubst,1
*do,j,1,8
d,nelem(neorder(i),j),temp,max_tem
*enddo
solve
t=t+dt1
time,t
solve
*do,j,1,8
ddele,nelem(neorder(i),j),temp
*enddo
t=t+dt-2*dt1
time,t
nsubst,nsub1
solve
*enddo
/image,save,1.bmp
t=t+10000
time,t
nsubst,nsub2
solve
finish
/post1
set,last
plnsol,temp,,0,1.0
/image,save,temp1,bmp
plnsol,s,eqv,2,1.0
/image,save,eqv1,bmp
/post26
/axlab,x,Time(s)
/axlab,y,Seqv(Mpa)
ansol,2,1,s,eqv,n1
ansol,3,2,s,eqv,n2
ansol,4,16,s,eqv,n16
ansol,5,23,s,eqv,n23
ansol,6,36,s,eqv,n36
xvar,1
plvar,2,3,4,5,6,
/image,save,eqv2,bmp
/axlab,x,time(s)
/axlab,y,temp(c)
nsol,7,1,temp,,n1
nsol,8,2,temp,,n2
nsol,9,23,temp,,n23
nsol,10,16,temp,,n16
nsol,11,36,temp,,n36
plvar,7,8,9,10,
/image,save,temp2,bmp
finish |