问一个简单问题

laoyin2

为什么这个命令流我拉伸后不对称呢

FINISH  ! Make sure we are at BEGIN level   
/CLEAR,NOSTART  ! Clear model since no SAVE found   
/NOPR   
/PMETH,OFF,0
KEYW,PR_SET,1   
KEYW,PR_STRUC,1



/PREP7  
!*********************面
ET,1,plane42


!**********************混凝土
ET,2,SOLID45

!********************8围岩
ET,3,SOLID65




!******************************混凝土mp1
!18.96mpa(兆帕)
MP,EX,1,20000e06
MP,PRXY,1,0.2  


!**************************************岩土mp2


MP,EX,2,300e6
MP,NUXY,2,0.4
mp,dens,2,2e3



!*******************************点  
K, ,0,0,,   
K, ,-4,-4,,
K, ,-4,4,,  
K, ,4,-4,,  
K, ,4,4,,   


k,,2,2,
k,,2,-2
k,,-2,-2
k,,-2,2


!************************线***************************************


LSTR,       2,       4  
LSTR,       4,       5  
LSTR,       5,       3  
LSTR,       3,       2


LSTR,       3,       9  
LSTR,       6,       5  
LSTR,       7,       4  
LSTR,       8,       2  

larc,6,7,1,2.828
larc,7,8,1,2.828
larc,8,9,1,2.828
larc,9,6,1,2.828


al,9,10,11,12


al,3,5,6,12
al,4,5,8,11
al,1,7,8,10
al,2,6,7,9

lsel,s,line,,1,12,1

lesize,all,,,4


asel,s,area,,2,5,1
mshkey,3
amesh,all

asel,s,area,,1
mshkey,3
amesh,all




!******************************************拉伸**********************************
mat,1                    
! 混凝土
type,2
!
esize,,5
! Anzahl der Elemente in L鋘gsrichtung
asel,s,area,,1           
! 选面从几到几
vext,all,,,,,20         
! 拉伸50m
allsel


!拉伸
mat,2                    
! 岩土
type,3
!砌衬
esize,,5
! Anzahl der Elemente in L鋘gsrichtung
asel,s,area,,2,5,1
! 选面从几到几
vext,all,,,,,20         
! 拉伸50m
ALLSEL,ALL  


!**********************清除面
asel,s,type,,1
aclear,all

!重力荷载
acel,,10 !加速度


!边界条件

nsel,s,loc,y,-4
! Markiert alle Punkte in y-Richtung auf H鰄e null
D,All,UY,0
! Freiheitsgrad UY festgehalten ...

allsel
nsel,s,loc,x,4
d,all,ux,0

allsel
nsel,s,loc,x,-4
d,all,ux,0

allsel
nsel,s,loc,z,0
d,all,uz,0

allsel
nsel,s,loc,z,20
d,all,uz,0



!last********************************************************last
/solu
!under face!可以不考虑  以后加荷载







!solution*********************************


/SOL  
!nsel,all
alls
eplo
ALLSEL



/solu                !进入求解器
antype,static            !指定为静力学求解器
pred,on                !打开时间步长预测器
lnsrch,on            !打开线性搜索
nlgeom,on            !打开大变形求解
!pcg 法
!EQSLV,PCG, ,1,
!nropt,full            !设定全N-R求解
OUTRES,ALL,all,      !最后一步，all，所有项
!NSUBST,subs,100,10
!Substepps
NEQIT,1000   
!迭代次数
!CNVTOL,F, ,0.01,2, ,   !精度是0.01

laoyin2

我举个简单例子 我在做隧道模拟但是为什么拉伸后y轴在重力下不对称呢

laoyin2

模型很简单 就是2种材料 一个圆形 一个方形 在重力下变形为什么沉降不一样呢 报酬50元 帮我看出来