谢谢解答!也不是这个问题,我是平面问题,当然没有z方向的位移!我贴出命令流,请高手指点下
COM
fini
/cle
*set,d,4
*set,a,3
*set,b,4
*set,s1,0.4
*set,x,20
*set,y,20
*set,H,2.9e7!垂直应力
*set,V,8.7e7!水平应力
*set,k,2.9e3!密度
*set,e,50e9!弹性模量
*set,u,0.3!泊松比
*set,c,0.1e6!粘结力
*set,j,30!内摩擦角
/prep7
ET,1,PLANE182
MPTEMP,1,0
MPDATA,DENS,1,,k!密度
!dp模型
MPDATA,EX,1,,e
MPDATA,PRXY,1,,u
!TB,DP,1,,,
!TBMODIF,1,1,c
!TBMODIF,1,2,j
!TBMODIF,1,3,
!建模
WPCSYS,-1,0
RECTNG,,x,,y,
WPCSYS,-1,0
CYL4,x/2,y/2,0,0,d/2,180
wpoff,x/2,y/2
RECTNG,-b/2,b/2,0,-a,
ASEL,S, , , 2
ASEL,A, , , 3
CM,m,AREA
ALLSEL,ALL
ASBA, 1,m
ALLSEL,ALL
!切分
WPCSYS,-1,0
wpoff,x/2,y/2
wpoff,-d/2-1
wpro,,,90
ASBW, all
WPCSYS,-1,0
wpoff,x/2,y/2
wpro,,,90
ASBW, all
WPCSYS,-1,0
wpoff,x/2,y/2
wpoff,d/2+1
wpro,,,90
ASBW, all
WPCSYS,-1,0
wpoff,x/2,y/2
wpro,,90
ASBW, all
WPCSYS,-1,0
wpoff,x/2,y/2
wpoff,,d/2+1
wpro,,90
ASBW, all
WPCSYS,-1,0
wpoff,x/2,y/2-a
wpro,,90
ASBW, all
WPCSYS,-1,0
wpoff,x/2-b/2,y/2-a
wpro,,,90
ASBW, 20
WPCSYS,-1,0
wpoff,x/2+b/2,y/2-a
wpro,,,90
ASBW, 19
NUMCMP,KP
NUMCMP,LINE
NUMCMP,AREA
!网格划分
TYPE, 1
MAT, 1
ESIZE,s1,0,
AMESH,15!巷道拱顶部分
AMESH,16!
AMESH,4
AMESH,5
AMESH,1
AMESH,2
AMESH,13
AMESH,14
ESIZE,s1,0,
AMESH,9!巷道直墙部分
AMESH,10!
ESIZE,s1,0,
AMESH,7!巷道直墙部分地板
AMESH,12!
AMESH,3
AMESH,11
AMESH,6
AMESH,8
AMESH,17
AMESH,18
WPCSYS,-1,0
!约束边界条件
NSEL,S,LOC,X,0
NSEL,A,LOC,X,x
CM,V,NODE
ALLSEL,ALL
NSEL,S,LOC,Y,0
NSEL,A,LOC,Y,Y
CM,H,NODE
ALLSEL,ALL
!施加荷载
SF,H,PRES,H
SF,V,PRES,V
NROPT,full,,
/SOL
ANTYPE,0
NLGEOM,1
|