Modeling of 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’);%显示正向刚度并做切片
隧道和岩爆模型(切片)
水力压裂(切片)