一个井筒建模出问题 请大家帮忙 谢谢

chenhonglei

具体步骤:
1.先模拟土的初始地应力
2.开挖并改变材料参数进行支护.
但是在第二步中总是出现单元压碎的现象 我把子步弄的很大 也是不行 请高手帮忙看一下程序 是不是合理 谢谢
/COM,  Structural   
/PREP7
!******材料模型********  
ET,1,PLANE182   
KEYOPT,1,1,0
KEYOPT,1,3,1
KEYOPT,1,6,0
KEYOPT,1,10,0   
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,2e-3
MPDATA,EX,1,,370
MPDATA,PRXY,1,,0.28
TB,CREE,1,1,4,6
TBTEMP,0
TBDATA,,-1.853e-4,0.725,-1.829,0,,   
MPTEMP,,,,,,,,  
MPTEMP,1,0
MPDATA,DENS,2,,2e-3
MPDATA,EX,2,,276
MPDATA,PRXY,2,,0.33
TB,CREE,2,1,4,6
TBTEMP,0
TBDATA,,-8.550e-4,0.732,-1.8648,0,,   
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,3,,2e-3
MPDATA,EX,3,,183
MPDATA,PRXY,3,,0.38
TB,CREE,3,1,4,6
TBTEMP,0
TBDATA,,-2.86714e-4,2.102,-1.9079,0,,   
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,4,,0.05e-3
MPDATA,EX,4,,0.421  
MPDATA,PRXY,4,,0.3
TB,CREE,4,1,4,6
TBTEMP,0
TBDATA,,1E-60,0,0,0,,
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,5,,0.05e-3
MPDATA,EX,5,,2.4
MPDATA,PRXY,5,,0.3
TB,CREE,5,1,4,6
TBTEMP,0
TBDATA,,1E-60,0,0,0,,
MPTEMP,,,,,,,,  
MPTEMP,1,0
MPDATA,DENS,6,,0.05e-3  
MPDATA,EX,6,,40000  
MPDATA,PRXY,6,,0.3
TB,CREE,6,1,4,6
TBTEMP,0
TBDATA,,1E-60,0,0,0,,
MPTEMP,,,,,,,,  
MPTEMP,1,0
MPDATA,DENS,7,,2e-3  
MPDATA,EX,7,,40000  
MPDATA,PRXY,7,,0.2
TB,CREE,7,1,4,6
TBTEMP,0
TBDATA,,1E-60,0,0,0,,  
MPTEMP,,,,,,,,  
MPTEMP,1,0
MPDATA,DENS,8,,2e-3
MPDATA,EX,8,,20.4   
MPDATA,PRXY,8,,0.28
TB,CREE,8,1,4,6
TBTEMP,0
TBDATA,,1E-60,0,0,0,,  
!**************************
RECTNG,0,4.275,0,25,  
*do,I,25,32.5,0.5  
RECTNG,0,3.3,I,I+0.5,   
RECTNG,3.3,4.2,I,I+0.5,
RECTNG,4.2,4.275,I,I+0.5,
*enddo
RECTNG,0,3.3,33,60,
RECTNG,3.3,4.2,33,60,
RECTNG,4.2,4.275,33,60,  
RECTNG,4.275,5.775,0,60,
RECTNG,5.775,7.275,0,60,
RECTNG,7.275,9.275,0,60,
RECTNG,9.275,10.775,0,60,
RECTNG,10.775,12.275,0,60,
RECTNG,12.275,100,0,60,
ALLSEL
AGLUE,ALL
ASEL,ALL
AATT,8,,,,
!***********************************
LESIZE,all,0.5,,,,,,,
MSHKEY,0
MSHAPE,0,2D
AMESH,ALL
!*****************
/SOL
ANTYPE,TRANS
TRNOPT,FULL
LUMPM,OFF
!******边界条件的定义*******
LSEL,S,LOC,Y,-0.1,0.1
DL,ALL,,UY,,
LSEL,S,LOC,X,99.9,100.1
DL,ALL,,UX,,
LSEL,S,LOC,Y,59.9,60.1
SFL,ALL,PRES,9,
ALLSEL,ALL
RATE,OFF
TIMINT,OFF
NROPT,FULL
NLGEOM,1
DELTIM,0.005
OUTRES,ALL,1
TIME,0.01
ALLSEL,ALL
ISWRITE,1
SOLVE
ASEL,S,LOC,X,-0.1,3.3001
ASEL,R,LOC,Y,32.9,60.1
ESLA
EKILL,ALL
ASEL,S,LOC,X,3.299,4.2001
ASEL,R,LOC,Y,32.9,60.1
ESLA
MPCHG,7,ALL,
ASEL,S,LOC,X,4.199,4.2751
ASEL,R,LOC,Y,32.9,60.1
ESLA
MPCHG,5,ALL,
ASEL,S,LOC,X,4.274,5.776
ASEL,A,LOC,X,10.774,12.276
ESLA
MPCHG,3,ALL,
ASEL,S,LOC,X,5.774,7.276
ASEL,A,LOC,X,9.274,10.776
ESLA
MPCHG,2,ALL,
ASEL,S,LOC,X,7.274,9.276
ESLA
MPCHG,1,ALL,
ALLSEL
ISFILE,READ,sheji1,ist, ,0
LSEL,S,LOC,X,-0.1,4.276
LSEL,R,LOC,Y,59.9,60.1
SFLDELE,ALL,PRES
LSEL,S,LOC,X,3.29,3.31
LSEL,R,LOC,Y,32.9,60.1
DL,ALL,,ALL,,
ASEL,S,loc,x,0,4.276
ASEL,r,loc,y,32.4,33.1
ESLA,S
EKILL,ALL
DELTIM,0.002
TIME,0.02
ALLSEL
SOLVE

chenhonglei

怎么没人回复啊 顶起

chenhonglei

没人帮忙吗