Earthquake

Step 1Step 2Step 3

%build the geometrical model
clear;
fs.randSeed(2);%build random model
B=obj_Box;%build a box object
B.name='BoxEarthquake';
B.GPUstatus='auto';
B.ballR=1;%default value is 0.5, use 1 to increase speed
B.isClump=0;
B.distriRate=0.2;
B.sampleW=200;
B.sampleL=0;
B.sampleH=160;
%B.BexpandRate=4;
%B.PexpandRate=4;
B.type='topPlaten';
%B.type='TriaxialCompression';
B.setType();
B.buildInitialModel();%B.show();

B.setUIoutput();

d=B.d;%d.breakGroup('sample');d.breakGroup('lefPlaten');
%--------------end initial model------------
B.gravitySediment();
B.compactSample(1);%input is compaction time
%------------return and save result--------------
d.status.dispEnergy();%display the energy of the model

d.clearData(1);%clear dependent data
d.recordCalHour('Box1Finish');
save(['TempModel/' B.name '1.mat'],'B','d');
save(['TempModel/' B.name '1R' num2str(B.ballR) '-distri' num2str(B.distriRate)  'aNum' num2str(d.aNum) '.mat']);
d.calculateData();
d.show('aR');

%initilization
clear
load('TempModel/BoxEarthquake1.mat');
B.setUIoutput();%set output of message
d=B.d;
d.calculateData();
d.mo.setGPU('off');
d.getModel();%get xyz from d.mo

%---------set material of model
matTxt=load('Mats\Soil3.txt');%material of soil
Mats{1,1}=material('Soil1',matTxt,B.ballR);
Mats{1,1}.Id=1;
matTxt2=load('Mats\Rock1.txt');
Mats{2,1}=material('Rock1',matTxt2,B.ballR);%material of rock
Mats{2,1}.Id=2;
d.Mats=Mats;
%---------end set material of model
%---------cut the model and make layers
C=Tool_Cut(d);%cut the model
lSurf=load('slope/Earthquake.txt');%load the surface data
C.addSurf(lSurf);%add the surfaces to the cut
C.setLayer({'sample'},[1,2,3,4]);%set layers according geometrical data
gNames={'lefPlaten';'rigPlaten';'botPlaten';'layer1';'layer2';'layer3'};
d.makeModelByGroups(gNames);%build new model
%---------end cut the model and make layers
%---------set material of group
d.setGroupMat('layer1','Rock1');
d.setGroupMat('layer2','Soil1');
d.setGroupMat('layer3','Rock1');
d.groupMat2Model({'layer1','layer2','layer3'},2);
%---------end set material of group

%---------balance the model
d.balanceBondedModel0();%balance the bonded model without friction

%---------define a block on left side of the model to generate wave
mX=d.mo.aX(1:d.mNum);
leftBlockId=find(mX<0.05*max(mX));%choose element Id of block
d.addGroup('LeftBlock',leftBlockId);%add a new group
d.setClump('LeftBlock');%set the block clump
d.mo.zeroBalance();
%---------end define a block on left side of the model to generate wave

%---------balance the model
d.mo.bFilter(:)=true;
d.mo.dT=d.mo.dT*4;
d.balance('Standard',1);
d.mo.bFilter(:)=true;
d.balance('Standard',0.5);
d.mo.dT=d.mo.dT/4;
%---------end balance the model

d.show('aMatId');%show material of model
%---------save data
d.mo.setGPU('off');
d.clearData(1);
d.recordCalHour('BoxModel2Finish');
save(['TempModel/' B.name '2.mat'],'B','d');
save(['TempModel/' B.name '2R' num2str(B.ballR) '-distri' num2str(B.distriRate)  'aNum' num2str(d.aNum) '.mat']);
d.calculateData();

clear;
load('TempModel/BoxEarthquake2.mat');
B.setUIoutput();
d=B.d;
d.calculateData();
d.mo.setGPU('off');
d.getModel();%reset the initial status of the model
d.resetStatus();%initialize model status, which records running information

d.mo.isHeat=1;%calculate heat in the model
d.mo.isCrack=1;%record in the model
visRate=0.000001;
d.mo.mVis=d.mo.mVis*visRate;
gpuStatus=d.mo.setGPU('auto');
d.setStandarddT();

%------------move left boundary to generate compressive wave
d.moveBoundary('left',0.01,0,0);

totalCircle=40;
d.tic(totalCircle);
fName=['data/step/' B.name  num2str(B.ballR) '-' num2str(B.distriRate) 'loopNum'];
save([fName '0.mat']);%return;
for i=1:totalCircle
    d.mo.setGPU(gpuStatus);
    d.balance('Standard',0.01);
    d.clearData(1);
    save([fName num2str(i) '.mat']);
    d.calculateData();
    d.toc();%show the note of time
end

d.show('mV');

%-----------save data
d.mo.setGPU('off');
d.clearData(1);
d.recordCalHour('BoxModel3Finish');
save(['TempModel/' B.name '3.mat'],'B','d');
save(['TempModel/' B.name '3R' num2str(B.ballR) '-distri' num2str(B.distriRate)  'aNum' num2str(d.aNum) '.mat']);
d.calculateData();