【原创】用强度折减法计算岩质节理边坡的稳定系数
本意是用强度折减法计算岩质节理边坡的稳定系数，可就是得不出收敛值啊，请高手多多指教。这里没有用到接触，不知道是不是这个原因fini
/clea
/filnam,rock,1
/TITLE,the 2D analysis of rock mass slope
/PREP7
l_1=105 !底边长
l_2=30 !坡底长
l_3=55 !坡顶长
h_1=20 !坡底高
h_2=40 !坡全高
hh=5 !岩石被切割宽度
ll=15 !岩石被切割长度
alfa=30 !第一组节理倾角
beta=60 !第二组节理倾角
width=0.8 !节理宽度
f_trial=1.23
cohe1=1e6
ifrag1=38
cohe2=0.12e6
ifrag2=24
K,1,0,0,0,
K,2,l_1,0,0,
K,3,l_1,h_2,0,
K,4,(l_1-l_3),h_2,0,
K,5,l_2,h_1,0,
K,6,0,h_1,0,
A,1,2,3,4,5,6
*afun,DEG
ET,1,plane82 !定义单元1为平面4节点单元（岩石）
KEYOPT,1,3,2 !定义为平面应变
MP,EX,1,1e11
MP,PRXY,1,0.2
MP,dens,1,2500
TB,dp,1
TBDATA,1,cohe1/f_trial,atan(tan(ifrag1)/f_trial),0
ET,2,plane82 !定义单元2为平面4节点单元（节理）
KEYOPT,2,3,2 !定义为平面应变
MP,EX,2,1e8
MP,PRXY,2,0.3
MP,dens,2,1700
TB,dp,2
TBDATA,1,cohe2/f_trial,atan(tan(ifrag2)/f_trial),0
nn=NINT((l_1*sin(alfa)+h_1*cos(alfa))/hh)!第一组节理的条数
*do,i,1,nn
k,2000+i,-h_1/tan(alfa)+(i-1)*hh/sin(alfa),0,0
k,3000+i,(h_2-h_1)/tan(alfa)+(i-1)*hh/sin(alfa),h_2,0
l,2000+i,3000+i
*enddo
LGEN,2,7,20,1,width/sin(alfa),,,, !复制线7－20并左移产生第一组节理
ASBL,1,all !将面1由线划分
/PNUM,AREA,1 !显示面的编号
APLOT
mm=NINT((l_1*cos(alfa)+h_2*sin(alfa))/ll)+1 !第二组节理的条数
*do,i,1,mm
k,4000+i,-h_2/tan(beta)+(i-1)*ll/sin(beta),h_2,0
k,5000+i,(i-1)*ll/sin(beta),0,0
l,4000+i,5000+i
*enddo
LGEN,2,1,8,1,width/sin(beta),,,, !复制线1－8并左移产生第二组节理
ASBL,28,all,,,keep
ASBL,3,all,,,keep
ASBL,30,all,,,keep
ASBL,8,all,,,keep
ASBL,14,all,,,keep
ASBL,18,all,,,keep
ASBL,22,all,,,keep
ASBL,2,all
*do,i,1,mm-1
k,6000+i,-h_2/tan(beta)+(i-0.5)*ll/sin(beta),h_2,0
k,7000+i,(i-0.5)*ll/sin(beta),0,0
l,6000+i,7000+i
*enddo
LGEN,2,1,7,1,width/sin(beta),,,, !复制线1－8并左移产生第二组节理
ASBL,26,all,,,keep
ASBL,12,all,,,keep
ASBL,11,all,,,keep
ASBL,6,all,,,keep
ASBL,10,all,,,keep
ASBL,16,all,,,keep
ASBL,20,all,,,keep
ASBL,24,all
aadd,27,34,29,85,87,32,90,92,13,36,89,38,40,42,44,46,5,104,102
aadd,24,100,98,96,95,7,58,56,54,52,50,49,9,107,109,111,113,115,31
aadd,5,61,63,65,67,69,15,117,119,121,123,17,71,73,75,19,125,127,129
aadd,7,21,77,79,23,131,25,81,82
aesize,all,1 !定义单元大小为5
mshape,1,2d !定义网格划分为映射划分
mshkey,0
type,2 !用单元2
mat,2 !用材料2
real,2
ASEL,S,,,3,5,1
AMESH,all
type,1 !用单元1
mat,1 !用材料1
real,1
ASEL,INVERT !反选面积
AMESH,ALL
!LREFINE,ALL,,,1 !在节理面上加深网格
/SOLU
ANTYPE,static
NLGEOM,on
NSUBST,12,0,0
ARCLEN,1,0,0
AUTOTS,on
OUTRES,ALL,all
PRED,ON
ACEL,0,9.8,0
NSEL,S,LOC,Y,0
NSEL,A,LOC,X,0
NSEL,A,LOC,X,105
D,ALL,ALL
NSEL,ALL
SOLVE
FINISH
/POST1
/EFACE,1
AVPRIN,0, ,
PLNSOL,S,1,0,1
回复： 【原创】用强度折减法计算岩质节理边坡的稳定系数
请问是建模有问题呢还是/solu有问题啊？小弟才疏学浅，贴出命令流抛砖引玉啊回复： 【原创】用强度折减法计算岩质节理边坡的稳定系数
*afun,DEG是设定角度以度为单位，默认的是弧度nn=NINT((l_1*sin(alfa)+h_1*cos(alfa))/hh) /NINT计算一组节理的条数
循环是画出一组节理的线，在后面的程序再来对面进行分割，生成节理
上面的结尾约束还是有问题，应将
NSEL,S,LOC,Y,0
NSEL,A,LOC,X,0
NSEL,A,LOC,X,105
D,ALL,ALL
改为
NSEL,S,LOC,Y,0
D,ALL,ALL
NSEL,s,LOC,X,0
d，all，ux，0
NSEL,s,LOC,X,l_1
d，all，ux，0
这样约束更合理些，不过计算结果还是不理想，应该把节理中尖角的地方处理下，避免计算时应力集中，使线性方程发生畸变
Re:【原创】用强度折减法计算岩质节理边坡的稳定系数
已经添加到收藏夹了.Re:【原创】用强度折减法计算岩质节理边坡的稳定系数
楼主有没有强度折减计算土质变坡安全系数，或者搜索临界破坏面的例子，能不能共享一下。谢谢了 强:victory: :这个问题很复杂呀
页:
[1]