跳转至

SlopeNet

user_BoxSlopeNet1.m
 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
%make a box model, which will be put on a slope
clear;
fs.randSeed(1);%build random model
B=obj_Box;%build a box object
B.name='BoxSlopeNet';
B.GPUstatus='auto';
B.ballR=0.05;
B.isClump=0;
B.distriRate=0.2;
B.sampleW=1;
B.sampleL=1;
B.sampleH=1.2;
%B.BexpandRate=4;
%B.PexpandRate=4;
B.type='GeneralSlope';
%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(2);%input is compaction time
mfs.reduceGravity(d,10);%reduce the gravity of element

%------------return and save result--------------
d.status.dispEnergy();%display the energy of the model
d.show('aR');

d.mo.setGPU('off');
d.clearData(1);%clear dependent data
d.recordCalHour('Step1Finish');
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();%because data is clear, it will be re-calculated
user_BoxSlopeNet2.m
 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
%step1: packing the elements
clear;
load('TempModel/BoxSlopeNet1.mat');
B.d.mo.setGPU('off');
packBallR=B.ballR;%record the ballR of the small box block
packBoxObj=B.d.group2Obj('sample');%make struct object from a group
packBoxObj=mfs.moveObj2Origin(packBoxObj);
sphereObj=mfs.cutSphereObj(packBoxObj,B.sampleW/2-B.ballR);
sphereObj2=mfs.cutSphereObj(packBoxObj,B.sampleW/3-B.ballR);
boxObj=mfs.cutBoxObj(packBoxObj,B.sampleW/2,B.sampleW/2.5,B.sampleW/3);

%make a big box for the simulation
fs.randSeed(1);%random model seed, 1,2,3...
B=obj_Box;%declare a box object
B.name='BoxSlopeNet';
%--------------initial model------------
B.GPUstatus='auto';%program will test the CPU and GPU speed, and choose the quicker one
B.ballR=packBallR;
B.isShear=0;
B.isClump=0;%if isClump=1, particles are composed of several balls
B.distriRate=0.2;%define distribution of ball radius, 
B.sampleW=7;%width, length, height, average radius
B.sampleL=4;%when L is zero, it is a 2-dimensional model
B.sampleH=4;
B.BexpandRate=2;%boundary is 4-ball wider than 
B.PexpandRate=0;
B.type='botPlaten';%add a top platen to compact model
B.isSample=0;
%B.type='TriaxialCompression';
B.setType();
B.buildInitialModel();
d=B.d;
d.mo.setGPU('off');
d.showB=3;
cellW=B.sampleL/10;
cellH=B.sampleL/10;
netObj=mfs.denseModel(0.8,@mfs.makeNet,B.sampleL,B.sampleL/3,cellW,cellH,B.ballR);
fs.showObj(netObj);
netId=d.addElement(1,netObj);%add a slope boundary
d.addGroup('net',netId);%add a new group
d.rotateGroup('net','XY',90,0,0,0);%rotate the group along XZ plane
d.rotateGroup('net','XZ',180);%rotate the group along XZ plane
%@@@@@@@@@@@@@@@@@add dZ here 
d.moveGroup('net',B.sampleW*0.8,0,0);
d.setClump('net');
%find the boundary of the net and fix them
netX=d.mo.aX(netId);netY=d.mo.aY(netId);netZ=d.mo.aZ(netId);
yFilter=(netY==min(netY)|netY==max(netY));
zFilter=(netZ==min(netZ));
netBoundaryFilter=yFilter|zFilter;
boundaryNetId=netId(netBoundaryFilter);
d.addFixId('X',boundaryNetId);%fix the X-coordinate
d.addFixId('Y',boundaryNetId);
d.addFixId('Z',boundaryNetId);

%d.show('aR');return;
slopeW=6;slopeL=B.sampleL;slopeH=B.ballR;
slope=mfs.denseModel(0.6,@mfs.makeBox,slopeW,slopeL,slopeH,B.ballR);%make a pile struct
slopeId=d.addElement(1,slope,'boundary');%add a slope boundary
d.addGroup('slope',slopeId);%add a new group
rotateX=max(d.aX(d.GROUP.slope));rotateZ=max(d.aZ(d.GROUP.slope));
d.rotateGroup('slope','XZ',-30,rotateX,0,rotateZ);%rotate the group along XZ plane
d.moveGroup('slope',-min(d.aX(d.GROUP.slope)),0,-min(d.aZ(d.GROUP.slope)));%move the slope

% slopeW=6/4;slopeL=B.sampleL;slopeH=B.ballR;
% slope=mfs.denseModel(0.6,@mfs.makeBox,slopeW,slopeL,slopeH,B.ballR);%make a pile struct
% slopeId=d.addElement(1,slope,'boundary');%add a slope boundary
% d.addGroup('step',slopeId);%add a new group
% rotateX=max(d.aX(d.GROUP.slope));rotateZ=max(d.aZ(d.GROUP.slope));
% d.moveGroup('step',B.sampleW*0.4/2,0,B.sampleH*0.8/2);%move the slope

%import the packed box (from step 1) on the top left side of the model
sphereObjId=d.addElement(1,sphereObj);%add regular model elements (not boundary)
d.addGroup('sphere',sphereObjId);%add a new group
d.moveGroup('sphere',B.sampleW/8,B.sampleL*1/2,B.sampleH*3.05/4);
d.setClump('sphere');

sphereObjId=d.addElement(1,sphereObj2);%add regular model elements (not boundary)
d.addGroup('sphere2',sphereObjId);%add a new group
d.moveGroup('sphere2',B.sampleW*1.7/8,B.sampleL*1.38/4,B.sampleH*3.05/4);
d.setClump('sphere2');

boxObjId=d.addElement(1,boxObj);%add regular model elements (not boundary)
d.addGroup('box',boxObjId);%add a new group
%d.rotateGroup('box','XZ',10);
d.moveGroup('box',B.sampleW/8,B.sampleL*1/4,B.sampleH*2.84/4);
d.setClump('box');

d.delElement('botPlaten');%remove bottom platen
d.show('StressZZ');
%------------record data
d.mo.setGPU('off');
d.clearData(1);%clear dependent data in d
d.recordCalHour('BoxSlopeNet2Finish');
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();
user_BoxSlopeNet3.m
 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
clear;
load('TempModel/BoxSlopeNet2.mat');
B.setUIoutput();
d=B.d;
d.calculateData();
d.mo.setGPU('off');
d.getModel();%d.setModel();%reset the initial status of the model
d.status=modelStatus(d);%initialize model status, which records running information

d.mo.isHeat=1;%calculate heat in the model
visRate=0.0001;
d.mo.mVis=d.mo.mVis*visRate;%use low viscosity for dynamic simulation
d.mo.setGPU('auto');
d.setStandarddT();
%d.mo.isShear=0;
totalCircle=20;
d.tic(totalCircle);
%.mat files will be saved in the folder data/step
fName=['data/step/' B.name num2str(B.ballR) '-' num2str(B.distriRate) 'loopNum'];
save([fName '0.mat']);%return;
for i=1:totalCircle
    d.balance(50,d.SET.packNum*1.5);
    d.clearData(1);%clear data in d.mo
    save([fName num2str(i) '.mat']);
    d.calculateData();
    d.toc();%show the note of time
end

d.mo.setGPU('off');
d.clearData(1);
d.recordCalHour('BoxSlope3Finish');
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();