Modeling of MatDEM

Tunnel and rockburst model establishment code (Matlab language, easy to use)

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

weakFilter=mfs.getWeakLayerFilter(d.mo.aX,d.mo.aY,d.mo.aZ,dipD,dipA,strongT,weakT);%Divided into weak layer

sampleId=d.getGroupId(‘sample’);%Get group number

aWFilter=false(size(weakFilter));

aWFilter(sampleId)=true;

sampleWfilter=aWFilter&weakFilter;

d.addGroup(‘WeakLayer’,find(sampleWfilter));%Add weak layer groups

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);%Split tunnel hole

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

tunnelId=d.getGroupId('Tunnel');

d.delElement(tunnelId);%Cut out the tunnel hole

B.setPlatenFixId();%Set pressure plate

d.setGroupMat(‘WeakLayer’,‘WeakRock’);%Set up each group of materials

d.groupMat2Model({‘WeakLayer’},1);%Material property

d.showFilter(‘SlideY’,0.5,1,’aKN’);%Show forward stiffness and slice

Tunnel and rock burst model (section)

图片1.jpg

Hydraulic fracturing (section)

图片2.jpg