ACEI网站
应用实例
离散元试验

利用MatDEM二次开发建立各类数值模型

隧道和岩爆模型的建立代码Matlab语言,简单易用)

dipD=90;dipA=10;strongT=4;weakT=4;

weakFilter=mfs.getWeakLayerFilter(d.mo.aX,d.mo.aY,d.mo.aZ,dipD,dipA,strongT,weakT);%划分弱层

sampleId=d.getGroupId(‘sample’);%获得组编号

aWFilter=false(size(weakFilter));

aWFilter(sampleId)=true;

sampleWfilter=aWFilter&weakFilter;

d.addGroup(‘WeakLayer’,find(sampleWfilter));%增加弱层组

sX=d.mo.aX(sampleId);sY=d.mo.aY(sampleId);sZ=d.mo.aZ(sampleId);

Rrate=0.2;sampleId=d.getGroupId('sample');

dipD=60;dipA=30;radius=8;height=30;

mX=d.mo.aX(1:d.mNum);mY=d.mo.aY(1:d.mNum);mZ=d.mo.aZ(1:d.mNum);

columnFilter=mfs.getColumnFilter(sX,sY,sZ,dipD,dipA,radius,height);%划分隧道孔

 

d.addGroup('Tunnel',find(columnFilter));

tunnelId=d.getGroupId('Tunnel');

d.delElement(tunnelId);%将隧道孔切出

B.setPlatenFixId();%设定压力板

d.setGroupMat(‘WeakLayer’,‘WeakRock’);%设置各组材料

d.groupMat2Model({‘WeakLayer’},1);%材料性质

d.showFilter(‘SlideY’,0.5,1,’aKN’);%显示正向刚度并做切片


隧道和岩爆模型(切片)


图片1.jpg

水力压裂(切片)

图片2.jpg

Copyright ©2000-2018, Nanjing University, P.R. China. All Rights Reserved 版权所有 ©2000-2018