【讨论】一个焊接的例子

htjp

下面是龚曙光书中的一个例子,好像有问题,请各位高手看一看:
/BATCH
/TITLE,WELD ANALYSIS
/PREP7
/UNIT,SI
ET,1,13,4
ET,2,13,4
MPTEMP,1,20,500,800,1200,1500
MPDATA,EX,1,1,2.09E11,1.72E11,1.33E11,0.84E11,0.45E11
MPDATA,EX,2,1,2.16E11,1.78E11,1.42E11,1.04E11,0.68E11
MPDATA,EX,3,1,2.12E11,1.75E11,1.39E11,1.07E11,0.83E11
TB,BKIN,1,5
TBTEMP,20,1
TBDATA,1,315E6,0.209E11
TBTEMP,500,2
TBDATA,1,212E6,0.172E11
TBTEMP,800,3
TBDATA,1,153E6,0.133E11
TBTEMP,1200,4
TBDATA,1,93E6,0.084E11
TBTEMP,1500,5
TBDATA,1,33E6,0.045E11
TB,BKIN,2,5
TBTEMP,20,1
TBDATA,1,314E6,0.216E11
TBTEMP,500,2
TBDATA,1,203E6,0.178E11
TBTEMP,800,3
TBDATA,1,153E6,0.142E11
TBTEMP,1200,4
TBDATA,1,63E6,0.104E11
TBTEMP,1500,5
TBDATA,1,13E6,0.068E11
TB,BKIN,3,5
TBTEMP,20,1
TBDATA,1,330E6,0.212E11
TBTEMP,500,2
TBDATA,1,213E6,0.175E11
TBTEMP,800,3
TBDATA,1,153E6,0.139E11
TBTEMP,1200,4
TBDATA,1,73E6,0.107E11
TBTEMP,1500,5
TBDATA,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,0.30
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
CSYS,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
D,ALL,UX,0
*GET,miny,NODE,,MNLOC,Y
NSEL,S,LOC,Y,miny
D,ALL,UY,0
NSEL,ALL
*GET,minx,NODE,,MNLOC,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),EQ,1,THEN
ii=ii+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
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
n1(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
mine=0
*DIM,ne2,,d
*DIM,nex2,,d
*DIM,ney2,,d
*DO,i1,1,d
ESEL,U,ELEM,,mine
*GET,nse2,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(i1)=mine
*ENDDO
maxe=0
ESEL,NONE
*DO,i,1,c
*IF,n1(i),NE,0,THEN
ESEL,A,ELEM,,n1(i)
*ENDIF
*ENDDO
EPLOT
*DIM,nel,,c
*DIM,nes1,,c
*DIM,ney1,,c
*DO,i1,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
ne1(ii)=i
*ENDIF
*ENDDO
*DO,i,1,nse3
*GET,ney1(i),ELEM,ne1(i),CENT,Y
*GET,nex1(i),ELEM,ne1(i),CENT,X
*ENDDO
maxy=-1e10
maxx=-1e10
*DO,i,1,nse3
*IF,ney1(i),GT,maxy,THEN
maxy=ney1(i)
maxx=nex1(i)
maxe=ne1(i)
*ELSE
*IF,ney1(i),EQ,maxy,THEN
*IF,nex1(i),LT,maxx,THEN
maxy=ney1(i)
maxx=nex1(i)
maxe=ne1(i)
*ENDIF
*ENDIF
*ENDIF
*ENDDO
neorder(i1+d)=maxe
*ENDDO
max_tem=1500
dt1=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+dt1
TIME,t
NSUBST,1
*DO,j,1,4
D,NELEM(neorder(i),j),TEMP,max_tem
*ENDDO
SOLVE
t=t+dt1
TIME,t
SOLVE
*DO,j,1,4
DDELE,NELEM(neorder(i),j),TEMP
*ENDDO
t=t+dt-2*dt1
TIME,t
NSUNST,nsub1
SOLVE
*ENDDO
t=t+5000
TIME,t
NSUBST,nsub2
SOLVE
FINISH

jnjl

怎么了？哪个地方有问题直接标出来啊！

yayasmiling

这么长，不好看啊
还是直接说出哪里有毛病比较方便一点
呵呵

anzi888

1把ESEL,U,ELEM,,mine 放到neorder(i1)=mine 的后边
在
2在mine=01后边加:
asel,s,area,,nna
esla

htjp

谢谢,anzi888师兄

fangdian

好厉害！
但是在这个程序中还是会出现warning：
The mesh of area 2 contains PLANE13 triangles, which are much too stiff
in bending.  Use quadratic (6- or 8-node) elements if possible.

怎么办？

谢谢