【讨论】请大牛讨论和指正：ansys+岩土＋接触＋单元生死

flac-ansys

ANSYS到底能不能用在岩土工程？下为一课题的ANSYS命令流，请大牛跑后将问题指正：
finish
/cle
/units,si         
a_1=35   
/prep7
！建模
wprota,,-a_1
csys,4
K,1,,-5.25,,      
K,2,,,,            
K,3,,5.65,,   
K,4,-9.8/2,8.08-5.25,,   
K,5,-9.8/2,-5.25,,  
  
K,11,,-7,-40,      
K,12,,,-40,           
K,13,,7,-40,   
K,14,-7,,-40,   
K,15,-7,-7,-40,  
  
l,1,3
larc,3,4,2,5.65
l,4,5
l,5,1
l,11,13
larc,13,14,12,7
l,14,15
l,15,11
l,1,11
l,3,13
l,4,14
l,5,15
  
al,1,2,3,4
al,5,6,7,8
al,2,11,6,10
al,3,12,7,11
al,4,9,8,12
al,1,9,5,10
va,1,2,3,4,5,6
  
csys,0   
KGEN,2,11, , , ,-20, ,100,0  
KGEN,2,111, , , , ,60,1,0
KGEN,2,111, , , , ,-30,2,0
KGEN,2,112,113,1, ,80,,2,0
l,11,111
l,111,112
l,111,113
l,114,112
lang,14,1,90,
lang,16,1,60,
VDRAG,2, , , , , ,15
KWPAVE,113  
wprota,,a_1
VSBW,2   
VDELE,3, , ,1
VDRAG,1, , , , , ,20
KWPAVE,2
wprota,,-a_1
csys,4
KGEN,2,3,10,7, ,2,,200,0
KGEN,2,13,,,,2,,200,0
l,114,115  
lang,37,213,90,
lang,37,203,90,
lang,37,210,90,
LDELE,16,37,21,1
l,17,115
lang,16,213,90,
a,1,6,112,9
a,1,6,111,11
a,11,111,113,8
a,9,1,203,210
a,1,11,213,203
a,11,8,24,213
a,203,22,23,210
a,203,213,21,22
a,213,24,115,21
KGEN,2,115, , ,-30, , , ,0  
l,25,115
ASEL,S,,,13,14
asel,a,,,19,26
VDRAG,all, , , , , ,53
VPTN,all
NUMMRG,ALL, , , ,LOW
NUMCMP,ALL  
ALLSEL
VGEN,2,5           
VDELE,5, , ,1
  
  
EHNT=3.5e10         
UHNT=0.1667      
DENSHNT=2600        
EWEIYAN=7.5e9      
UWEIYAN=0.21         
DENSWEIYAN=2400      
ET,1,SOLID45         
R,1                  
MP,EX,1,EHNT         
MP,PRXY,1,UHNT      
MP,DENS,1,DENSHNT   
mp,mu,1,0.40         
  
ET,2,SOLID45              
R,2                     
MP,EX,2,EWEIYAN         
MP,PRXY,2,UWEIYAN        
MP,DENS,2,DENSWEIYAN     
mp,mu,2,0.40
tb,dp,2
tbdata,1,1.35e6,47.7,0      
vsel,s,,,13
VATT,1,1,1,
vsel,inve
VATT,2,2,2,
allsel
LESIZE,90, , ,10
LESIZE,91, , ,7
LESIZE,92, , ,8
LESIZE,93, , ,5
LESIZE,99, , ,20
LESIZE,101, , ,20
VSWEEP,13
LESIZE,87, , ,2
LESIZE,2, , ,7
LESIZE,3, , ,8
LESIZE,83, , ,11
LESIZE,62, , ,9
LESIZE,63, , ,18
LESIZE,12, , ,20
LESIZE,10, , ,20
VSWEEP,10
LESIZE,1, , ,10
LESIZE,4, , ,5
LESIZE,23, , ,10
LESIZE,20, , ,10
VSWEEP,1
LESIZE,8, , ,5
LESIZE,33, , ,8
LESIZE,22, , ,8
VSWEEP,6
VSWEEP,11
VSWEEP,12
VSWEEP,3
VSWEEP,2
VSWEEP,4
LESIZE,18, , ,8
VSWEEP,7
VSWEEP,9
VSWEEP,8
allsel
  
/prep7
ALLSEL
R,3  
REAL,3   
ET,3,170
ET,4,174
keyopt,4,9,0         
keyopt,4,12,2         
asel,s,,,61,63,1
cm,_target,area
type,3
nsla,s,1
esurf,all
allsel
asel,s,,,2,3,1
asel,a,,,42
cm,_contact,area
type,4
nsla,s,1
esurf,all
allsel
  
vsel,s,,,1
eslv,r
cm,death1,elem
allsel
vsel,s,,,13
eslv,r
cm,death2,elem
allsel
esel,s,ename,,170
esel,a,ename,,174
cm,death3,elem
allsel
csys,0
acel,,9.8   
!加边界条件
ASEL,S,LOC,X,0      
DA,all,SYMM
allsel
ASEL,S,LOC,X,-30     
DA,all,ux
allsel
ASEL,S,LOC,y,-48.67712     
DA,all,uy
allsel
ASEL,S,LOC,z,-58.75105      
ASEL,A,LOC,z,31.24895      
DA,all,uz
allsel
  
!求解
/solu
antype,static
lnsrch,on  
pred,on  
nropt,full,,on
nlgeom,on
autots,on
nsub,10,20,10
  
time,1
allsel
esel,s,,,death1
EKILL,ALL  
esel,s,live
nsle,s
nsel,inve
d,all,all,0
allsel
solve
SAVE
  
time,2
allsel
ddele,all,all
esel,s,,,death2
esel,a,,,death3
EKILL,ALL  
esel,s,live
nsle,s
nsel,inve
d,all,all,0
allsel
ASEL,S,LOC,X,0     
DA,all,SYMM
allsel
ASEL,S,LOC,X,-30     
DA,all,ux
allsel
ASEL,S,LOC,y,-48.67712      
DA,all,uy
allsel
ASEL,S,LOC,z,-58.75105      
ASEL,A,LOC,z,31.24895      
DA,all,uz
allsel
solve
SAVE
  
time,3
allsel
ddele,all,all
esel,s,,,death2
esel,a,,,death3
ealive,all
allsel
esel,s,live
nsle,s
nsel,inve
d,all,all,0
allsel
ASEL,S,LOC,X,0      
DA,all,SYMM
allsel
ASEL,S,LOC,X,-30      
DA,all,ux
allsel
ASEL,S,LOC,y,-48.67712   
DA,all,uy
allsel
ASEL,S,LOC,z,-58.75105      
ASEL,A,LOC,z,31.24895      
DA,all,uz
allsel
solve
SAVE
time,4
SFA,60,1,PRES,1022466.7   
solve
SAVE

flac-ansys

不要认为问题简单！其实岩土的问题能在软件中实现，关键是要看仿真后的数据分析能力！请大牛和斑竹回复！！！我谢谢先！