我下下来了,不要花钱了,这个内容是:
/title,weld analysis by "element birth and death"
/prep7
/unit,si
et,1,13,4
et,2,13,4
mptemp,1,20,500,800,1200,1500
mpdate,ex,1,1,2.09e11,1.72e11,1.33e11,0.84e11,0.45e11
mpdate,ex,2,1,2.16e11,1.78e11,1.42e11,1.04e11,0.68e11
mpdate,ex,3,1,2.12e11,1.75e11,1.39e11,1.07e11,0.83e11
tb,bkin,1,5 !以表格的方式输入经典双曲线性随动强化参数
tbtemp,20,1
tbdate,1,315e6,0.209e11
tbtemp,500,2
tbdate,1,212e6,0.172e11
tbtemp,800,3
tbdate,1,153e6,0.133e11
tbtemp,1200,4
tbdate,1,93e6,0.084e11
tbtemp,1500,5
tbdate,1,33e6,0.045e11
tb,bkin,2,5 !材料2
tbtemp,20,1
tbdate,1,314e6,0.216e11
tbtemp,500,2
tbdate,1,203e6,0.178e11
tbtemp,800,3
tbdate,1,153e6,0.142e11
tbtemp,1200,4
tbdate,1,63e6,0.104e11
tbtemp,1500,5
tbdate,1,13e6,0.068e11
tb,bkin,3,5 !材料3
tbtemp,20,1
tbdate,1,330e6,0.212e11
tbtemp,500,2
tbdate,1,213e6,0.175e11
tbtemp,800,3
tbdate,1,153e6,0.139e11
tbtemp,1200,4
tbdate,1,73e6,0.107e11
tbtemp,1500,5
tbdate,1,13e6,0.083e11
mp,dens,1,7850
mp,dens,2,7770
mp,dens,3,7860
mp,alpx,1,1.23e-5
mp,alpx,2,1.38e-5
mp,alpx,3,1.48e-5
mp,nuxy,1,030
mp,nuxy,2,0.29
mp,nuxy,3,0.29
mp,kxx,1,33
mp,kxx,2,32
mp,kxx,3,34
mp,c,1,561
mp,c,2,683
mp,c,3,983
mp,murx,1,1
mp,murx,2,1
mp,murx,3,1
mp,reft,1,20
mp,reft,2,1500
mp,reft,3,20
!下面将生成分析模型和网格
*afun,deg
csy,0
width=1
height=0.3
ar=0.14*tan(30)
ex=0.01*tan(60)
df=ar+ex
cf=df/cos(60)
k,1,0,0,0
k,2,width/2-0.01-ar,0,0
k,3,width/2+0.01-ar,0,0
k,4,width,0,0,
k,5,width/2-0.01,height/2-0.01,0
k,6,width/2+0.01,height/2-0.01,0
k,7,width/2-0.01,height/2+0.01,0
k,8,width/2+0.01,height/2+0.01,0
k,9,0,height,0
k,10,width/2-0.01-ar,height,0
k,11,width/2+0.01+ar,height,0
k,12,width,height,0
n,1,width/2+0.01+ar,height,0
n,2,width,height,0
n,3,width/2+0.01+ar,0,0
larc,10,11,7,cf
larc,2,3,5,cf
a,1,2,5,7,10,9
a,10,7,5,2,3,6,8,11
a,3,4,12,11,8,6
/triad,off
/plots,info,0
/devices,vector,1
esize,0.03
type,2
mat,2
amesh,2
esize,0.03
type,1
mat,1
amesh,1
mat,3
amesh,3
nummrg,all,,,,low
numcmp,all
eplot
/image,save,vmesh,bmp
/devices,vector,0
finish
!下面开始求解分析
/solu
antype,4
trnopt,full
nsel,all
*get,minx,node,,mnloc,x
nsel,s,loc,x,minx
*get,miny,node,,mnloc,y
nsel,r,loc,y,miny
d,all,uy,0
nsel,all
*get,minx,node,,minloc,x
nsel,s,loc,x,minx
d,all,temp,20
nsel,all
*get,maxx,node,,mxloc,x
nsel,s,loc,x,maxx
d,all ,temp,20
*set,nna,2
esel,all
*get,emax,elem,,num,max
asel,s,area,,nna
esla
*get,nse,elem,count
*dim,ne,,nse
*dim,n1,,nse
*dim,n2,,nse
*dim,nex,,nse
*dim,ney,,nse
*dim,neorder,,nse
ii=0
*do,i,1,emax
*if,esel(i),1,then
ii=i+1
ne(ii)=i
*endif
enddo
*do,i,1,nse
*get,ney(i),elem,ne(i),cent,y
*get,nex(i),elem,ne(i),cent,x
*enddo
!将 x型焊缝区得单元分成两部分
d=0
*do,i,1,nse
*if,ney(i),ge,height/2,then
d=d+1
n2(d)=ne(i)
*endif
*enddo
c=0
*do,i,1,nse
*if,ney(i),lt,height/2,then
c=c+1
nl(c)=ne(i)
*endif
*enddo
esel,none
*do,i,1,d
*if,n2(i),ne,0,then
esel,a,elem,,n2(i)
*endif
*enddo
eplot
!对焊缝区上的部分单元按质心y坐标排序
mine=0
*dim,ne2,,d
*dim,nex2,,d
*dim,ney2,,d
*do,il,1,d
esel,u,elem,,count
ii=0
*do,i,1,emax
*if,esel(i),eq,1,then
ii=ii+1
ne2(ii)=i
*endif
*enddo
*do,i,1,nse2
*get,ney2(i),elem,ne2(i),cent,y
*get,nex2(i),elem,ne2(i),cent,x
*enddo
miny=1e20
minx=1e20
*do,i,1,nse2
if,ney2(i),lt,miny,then
miny=ney2(i)
minx=nex2(i)
mine=ne2(i)
*else
*if,ney2(i),eq,miny,then
*if,nex2(i),lt,minx,then
miny=ney2(i)
minx=nex2(i)
mine=ne2(i)
*endif
*endif
*endif
*enddo
neorder(il)=mine
*enddo
!坐标排序 y,与上类似
maxe=0
esel,none
ii=0
*do,i,1,c
*if,n1(i),ne,0,then
esel,a,elem,,n1(i)
*endif
*enddo
eplot
*dim,nel,,c
*dim,nexl,,c
*dim,ney1,,c
*do,il,1,c
esel,u,elem,,maxe
*get,nse3,elem,,count
ii=0
*do,i,1,emax
*if,esel(i),eq,1,then
ii=ii+1
nel(ii)=i
*endif
*enddo
*do,i,1,nse3
*get,neyl(i),elem,nel(i),cent,y
*get,neyl(i),elem,nel(i),cent,x
*enddo
maxy=-1e10
maxx=-1e10
*do,i,1,nse3
*if,neyl(i),gt,maxy,then
maxy=neyl(i)
maxx=nexl(i)
maxe=nel(i)
*else
*if,neyl(i),eq,maxy,then
!428
*if,nex(i),lt,maxx,then
maxy=neyl(i)
maxx=nexl(i)
maxe=nel(i)
*endif
*endif
*endif
*enddo
neorder(il+d)=maxe
*enddo
max_tem=1500
dtl=1e-3
dt=5
t=0
esel,all
eplot
/auto,1
/replot
*do,i,1,nse
ekill,neorder(i)
esel,s,live
eplot
*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
eplot
esel,all
t=t+dtl
time,t
nsubst,1
*do,j,1,4
d,nelem( neorder(i),j),temp,max_tem
*enddo
solve
t=t+dtl
time,t
solve
*do,j,1,4
ddele,nelem(neorder(i),j),temp
*enddo
t=t+dt-2*dtl
time,t
nsubst,nsub1
solve
*enddo
t=t+50000
time,t
nsubst,nsub2
solve
finish
/post26
nsol,2,379,temp,temp379
/gropt,axnsc,1.5
plvar,2
/image,save,temp379,bmp
finish
/post1
/seg,dele
/cont,1,15,0,1200e6/16,1200e6
/dscale,1,1.0
avprin,0,0
avres,1
/seg,multi,stress,0.1
esel,all
*do,i,1,nse
esel,u,elem,,neorder(i)
*enddo
*do,i,1,nse
esel,u,elem,,neorder(i)
*enddo
*do,i,1,nse
esel,a,elem,,neorder(i)
set,(i-1)*3+1,1
plnsol,s,eqv
*do,j,1,nsub1
set,(i-1)*3+3,j
plnsol,s,eqv
*enddo
*enddo
*do,i,1,nsub2
set,(nse-1)*3+4,i
plnsol,s,eqv
*enddo
!430
/seg,off,stress,0.1
/image,save,stress,bmp
finish
/post1
/seg,dele
/cont,1,15,0,1500/16,1500
/dscale,1,1.0
avprin,0,0
avres,1
/seg,multi,temp,0.1
esel,all
*do,i,1,nse
esel,u,elem,neorder(i)
*enddo
*do,i,1,nse
esel,a,elem,,neorder(i)
set,(i-1)*3+1,1
plnsol,temp
*do,i,1,nsub1
set,(i-1)*3+3,j
plnsol,temp
*enddo
*enddo
!/image,save,tem1,bmp
*do,i,1,nsub2
set,(nse-1)*3+4,i
plnsol,temp
*enddo
/seg,off,temp,0.1
anim,1,1
!现实路径结果
set,last
csys,0
path,xdirection,2,10
ppath,1,1
ppath,2,2
pdef,,s,eqv
/gropt,axnsc,1.5,
plpath,seqv
/image,save,xpath,bmp
path,ydirection,2,10
ppath,1,1
ppath,2,3
pdef,,s,eqv
/gropt,axnsc,1.5,
plpath,seqv |