CN106706884A - Method and apparatus for determining development degree of rock cracks - Google Patents
Method and apparatus for determining development degree of rock cracks Download PDFInfo
- Publication number
- CN106706884A CN106706884A CN201710017322.0A CN201710017322A CN106706884A CN 106706884 A CN106706884 A CN 106706884A CN 201710017322 A CN201710017322 A CN 201710017322A CN 106706884 A CN106706884 A CN 106706884A
- Authority
- CN
- China
- Prior art keywords
- microfissure
- point
- principal component
- rock
- group
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0001—Type of application of the stress
- G01N2203/0003—Steady
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0014—Type of force applied
- G01N2203/0016—Tensile or compressive
- G01N2203/0019—Compressive
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/006—Crack, flaws, fracture or rupture
- G01N2203/0062—Crack or flaws
- G01N2203/0066—Propagation of crack
Abstract
The invention provides a method and an apparatus for determining the development degree of rock cracks. The development degree of all micro-crack groups in an axial extension direction and in a transverse expansion direction is quantified. The method comprises the following steps: constructing a grain flow model of a rock sample; acquiring a micro-crack set in the destroy process of the grain flow model and the coordinate value of every micro-crack point in the micro-crack set according to the constructed grain flow model; carrying out clustering based on the density of the micro-crack points according to the obtained micro-crack set and the coordinate value of every micro-crack point in the micro-crack set, and grouping the micro-crack points in the micro-crack set according to the distance of every clustering center to obtain multiple micro-crack groups; and analyzing the axial extension direction and the transverse expansion direction of every micro-crack group by adopting a main component analysis technology, and quantifying the development degree of every micro-crack group in the axial extension direction and in the transverse expansion direction. The method and the apparatus are suitable for the technical fields of rock soil test and engineering.
Description
Technical field
The present invention relates to geomaterial Research on Mechanical Properties field, a kind of determination rock crackle forming development journey is particularly related to
The method and device of degree.
Background technology
In recent years, the rock sample for having more research to be directed to prefabricated crack carries out failure test, probes into and is deposited in crack
In case, in rock failure process crack development trend.But because microfissure naked eyes are difficult to what is observed, research is most
The information such as the mechanics parameter of rock and the development trend in macroscopical crack is obtained, thus it is universal not to microcosmic mechanism aspect research
Foot.Meanwhile, the research for also having correlation is to combine numerical simulation software, and particularly discrete meta software carries out the simulation of traditional experiment,
Disclose the thin sight mechanism in rock failure process.Numerical simulation software can effectively catch as the strong supplement of traditional experiment
Grasp the thin generation development for seeing microfissure and the process for ultimately forming Macroscopic band.Related scholar combine laboratory test and
Crack development of the distinct element methods such as grain stream (Particle Flow Code, PFC) the research sample in the case of uniaxial compression becomes
Gesture, but because Cementation failure is regarded as microfissure by PFC softwares in simulation process, therefore the final result for presenting mostly can only
Judge the development trend in macroscopical crack roughly by the aggregation of microfissure, macroscopical crack development is described rather than accurate quantification
Trend and level of breakage.
The content of the invention
The technical problem to be solved in the present invention is to provide a kind of method and device for determining rock crackle forming development degree, to solve
Certainly the describe macroscopical crack development trend and level of breakage problem with being unable to accurate quantification existing for prior art.
In order to solve the above technical problems, the embodiment of the present invention provides a kind of method for determining rock crackle forming development degree, bag
Include:
Build the GRANULAR FLOW MODEL FOR of rock sample;
According to the GRANULAR FLOW MODEL FOR for building, the microfissure collection and the microfissure collection in GRANULAR FLOW MODEL FOR destructive process are obtained
In each microfissure point coordinate value;
The coordinate value of each microfissure point is concentrated according to the microfissure collection for obtaining and the microfissure, based on microfissure
The density of point is clustered, and according to the distance of each cluster centre, the microfissure point that the microfissure is concentrated is grouped, and is obtained
To multiple microfissure groups;
Using PCA, the axially extending direction and lateral magnification direction of each microfissure group are analyzed, and quantify each
Development degree of the microfissure group on axially extending direction and lateral magnification direction.
Further, the GRANULAR FLOW MODEL FOR for building rock sample includes:
S11, using particle and cohere structure granule sample, wherein, the size of the granule sample with carry out interior
The size of the rock sample of rock mechanics experiment is identical;
S12, sets particle rill evolution and coheres rill evolution;
S13, according to set particle rill evolution and cohere rill evolution, uniaxial compression is carried out to the granule sample
Virtual load, obtains the single compressing stress-strain curve of the granule sample;
S14, the single compressing stress-strain curve of the granule sample that will be obtained with carry out indoor rock mechanics examination
Single compressing stress-the strain curve of the rock sample tested is compared;
S15, if the difference of curve is in default first threshold, according to the particle rill evolution for setting and coheres thin sight
Parameter, builds the GRANULAR FLOW MODEL FOR of rock sample;
S16, otherwise, then adjusts and particle rill evolution and coheres rill evolution, and return to S13 and continue executing with.
Further, it is described according to the GRANULAR FLOW MODEL FOR for building, obtain the fine fisssure in the GRANULAR FLOW MODEL FOR destructive process
Gap collection and the microfissure concentrate the coordinate value of each microfissure point to include:
According to the GRANULAR FLOW MODEL FOR for building, virtual uniaxial compression is carried out along long axis direction, obtain the GRANULAR FLOW MODEL FOR and break
Microfissure collection and the microfissure during bad concentrate the coordinate value of each microfissure point.
Further, the microfissure collection and the microfissure that the basis is obtained concentrate the coordinate of each microfissure point
Value, the density based on microfissure point is clustered, and according to the distance of each cluster centre, the microfissure concentrated to the microfissure
Point is grouped, and obtaining multiple microfissure groups includes:
The coordinate value of each microfissure point is concentrated according to the microfissure, obtain microfissure point i and other microfissure points j it
Between apart from dij;
If cluster centre is surrounded by the surrounding's microfissure point with smaller local density, and cluster centre and other parts are close
There is a certain distance between the larger microfissure point of degree, calculate the local density ρ of microfissure point ii:
ρi=∑jχ(dij-dc);
Calculate microfissure point i has the distance between the microfissure point j of larger local density δ from otheri, wherein, if micro-
Crack point i has larger local density, then δi=maxj(dij);Otherwise, apart from δiBe by microfissure point i and other have compared with
Minimum range δ between the microfissure point j of big local densityiDetermine,
Wherein, if dij-dc<0, then χ (dij-dc)=1, otherwise, χ (dij-dc)=0;dcIt is to block distance, dijIt is microfissure
The distance between point i and other microfissure points j;It is smaller local density when local density is less than default Second Threshold;When
Local density is larger local density when being not less than default Second Threshold.
Further, the use PCA, analyzes axially extending direction and the lateral magnification of each microfissure group
Direction, and quantify development degree of each microfissure group on axially extending direction and lateral magnification direction and include:
Using the method for principal component analysis, each microfissure group is analyzed, obtain first principal component direction, second it is main into
Divide the contribution margin in direction, the contribution margin of first principal component and Second principal component,;
Wherein, the first principal component direction represents the axially extending direction of corresponding microfissure group, Second principal component, side
To the lateral magnification direction for representing corresponding microfissure group, first principal component direction is orthogonal with Second principal component, direction;First master
The contribution margin of composition represents development degree of the corresponding microfissure group on axially extending direction, the contribution margin table of Second principal component,
Show development degree of the corresponding microfissure group on lateral magnification direction.
The embodiment of the present invention also provides a kind of device for determining rock crackle forming development degree, it is characterised in that including:
Build module, the GRANULAR FLOW MODEL FOR for building rock sample;
Acquisition module, for according to the GRANULAR FLOW MODEL FOR for building, obtaining the microfissure collection in GRANULAR FLOW MODEL FOR destructive process
And the microfissure concentrates the coordinate value of each microfissure point;
Cluster module, the coordinate for concentrating each microfissure point according to the microfissure collection and the microfissure that obtain
Value, the density based on microfissure point is clustered, and according to the distance of each cluster centre, the microfissure concentrated to the microfissure
Point is grouped, and obtains multiple microfissure groups;
Quantization modules, for using PCA, analyze axially extending direction and the lateral magnification of each microfissure group
Direction, and quantify development degree of each microfissure group on axially extending direction and lateral magnification direction.
Further, the structure module includes:
Determining unit, for using particle and cohere structure granule sample, wherein, the size of the granule sample with
The size for carrying out the rock sample of indoor rock mechanics experiment is identical;
Setting unit, for setting particle rill evolution and cohering rill evolution;
Loading unit, for according to set particle rill evolution and cohere rill evolution, the granule sample is entered
Row uniaxial compression virtual load, obtains the single compressing stress-strain curve of the granule sample;
Comparing unit, for the single compressing stress-strain curve of the granule sample that will obtain with carry out interior
Single compressing stress-the strain curve of the rock sample of rock mechanics experiment is compared;
Construction unit, if for curve difference in default first threshold, according to set particle rill evolution
With cohere rill evolution, build the GRANULAR FLOW MODEL FOR of rock sample;
Returning unit, if for curve difference not in default first threshold, adjustment particle rill evolution and glutinous
Rill evolution is tied, and returns to the loading unit and continued executing with.
Further, the acquisition module, for according to the GRANULAR FLOW MODEL FOR for building, virtual single shaft being carried out along long axis direction
Compression, obtains the coordinate that microfissure collection and the microfissure in the GRANULAR FLOW MODEL FOR destructive process concentrate each microfissure point
Value.
Further, the cluster module includes:
First computing unit, the coordinate value for concentrating each microfissure point according to the microfissure obtains microfissure point i
With the distance between other microfissure points j dij;
Second computing unit, is surrounded, and gather for setting cluster centre by the surrounding's microfissure point with smaller local density
There is a certain distance between the larger microfissure point in class center and other local densities, calculate the local density of microfissure point i
ρi:
ρi=∑jχ(dij-dc);
3rd computing unit, for calculate microfissure point i have from other between microfissure point j of larger local density away from
From δi, wherein, if microfissure point i has larger local density, δi=maxj(dij);Otherwise, apart from δiIt is by microfissure point i
And other minimum range δ between there is the microfissure point j of larger local densityiDetermine,
Wherein, if dij-dc<0, then χ (dij-dc)=1, otherwise, χ (dij-dc)=0;dcIt is to block distance, dijIt is microfissure
The distance between point i and other microfissure points j;It is smaller local density when local density is less than default Second Threshold;When
Local density is larger local density when being not less than default Second Threshold.
Further, the quantization modules, for the method using principal component analysis, are analyzed to each microfissure group,
Obtain the contribution margin of first principal component direction, Second principal component, direction, the contribution margin of first principal component and Second principal component,;
Wherein, the first principal component direction represents the axially extending direction of corresponding microfissure group, Second principal component, side
To the lateral magnification direction for representing corresponding microfissure group, first principal component direction is orthogonal with Second principal component, direction;First master
The contribution margin of composition represents development degree of the corresponding microfissure group on axially extending direction, the contribution margin table of Second principal component,
Show development degree of the corresponding microfissure group on lateral magnification direction.
Above-mentioned technical proposal of the invention has the beneficial effect that:
In such scheme, by build GRANULAR FLOW MODEL FOR, obtain GRANULAR FLOW MODEL FOR destructive process in microfissure collection and
The microfissure concentrates the coordinate value of each microfissure point;Each is concentrated according to the microfissure collection for obtaining and the microfissure
The coordinate value of microfissure point, the density based on microfissure point is clustered, and according to the distance of each cluster centre, to the fine fisssure
The microfissure point that gap is concentrated is grouped, and obtains multiple microfissure groups;Using PCA, the axle of each microfissure group is analyzed
To bearing of trend and lateral magnification direction, and quantify development of each microfissure group on axially extending direction and lateral magnification direction
Degree.
Brief description of the drawings
Fig. 1 is the schematic flow sheet of the method for determination rock crackle forming development degree provided in an embodiment of the present invention;
Fig. 2 is data analysis schematic flow sheet provided in an embodiment of the present invention;
Fig. 3 is rock sample scale diagrams provided in an embodiment of the present invention;
Fig. 4 is the load-deformation curve under rock sample provided in an embodiment of the present invention and granule sample uniaxial loading
Schematic diagram;
Fig. 5 is particle body Model failure mode schematic diagram provided in an embodiment of the present invention;
Fig. 6 is particle body Model microfissure collection group result schematic diagram provided in an embodiment of the present invention;
Fig. 7 is principal component analysis principle schematic provided in an embodiment of the present invention;
Fig. 8 is cluster analysis provided in an embodiment of the present invention and principal component analysis result schematic diagram;
Fig. 9 is the structural representation of the device of determination rock crackle forming development degree provided in an embodiment of the present invention.
Specific embodiment
To make the technical problem to be solved in the present invention, technical scheme and advantage clearer, below in conjunction with accompanying drawing and tool
Body embodiment is described in detail.
The present invention is carried for existing describe macroscopical crack development trend and level of breakage problem with being unable to accurate quantification
For a kind of method and device for determining rock crackle forming development degree.
Embodiment one
Referring to shown in Fig. 1, a kind of method for determining rock crackle forming development degree provided in an embodiment of the present invention, including:
S101, builds the GRANULAR FLOW MODEL FOR of rock sample;
S102, according to the GRANULAR FLOW MODEL FOR for building, obtains the microfissure collection and described micro- in GRANULAR FLOW MODEL FOR destructive process
Concentrate the coordinate value of each microfissure point in crack;
S103, the coordinate value of each microfissure point is concentrated according to the microfissure collection for obtaining and the microfissure, is based on
The density of microfissure point is clustered, and according to the distance of each cluster centre, the microfissure point that the microfissure is concentrated is carried out
Packet, obtains multiple microfissure groups;
S104, using PCA, analyzes the axially extending direction and lateral magnification direction of each microfissure group, and measure
Change development degree of each microfissure group on axially extending direction and lateral magnification direction.
The method of the determination rock crackle forming development degree described in the embodiment of the present invention, by the GRANULAR FLOW MODEL FOR for building, obtains
Obtain the coordinate value that microfissure collection and the microfissure in GRANULAR FLOW MODEL FOR destructive process concentrate each microfissure point;According to obtaining
The microfissure collection and the microfissure concentrate the coordinate value of each microfissure point, the density based on microfissure point is gathered
Class, and according to the distance of each cluster centre, the microfissure point that the microfissure is concentrated is grouped, obtain multiple microfissures
Group;Using PCA, the axially extending direction and lateral magnification direction of each microfissure group are analyzed, and quantify each microfissure
Development degree of the group on axially extending direction and lateral magnification direction.
In the specific embodiment of the method for foregoing determination rock crackle forming development degree, further, the structure rock
The GRANULAR FLOW MODEL FOR (S101) of stone sample includes:
S11, using particle and cohere structure granule sample, wherein, the size of the granule sample with carry out interior
The size of the rock sample of rock mechanics experiment is identical;
S12, sets particle rill evolution and coheres rill evolution;
S13, according to set particle rill evolution and cohere rill evolution, uniaxial compression is carried out to the granule sample
Virtual load, obtains the single compressing stress-strain curve of the granule sample;
S14, the single compressing stress-strain curve of the granule sample that will be obtained with carry out indoor rock mechanics examination
Single compressing stress-the strain curve of the rock sample tested is compared;
S15, if the difference of curve is in default first threshold, according to the particle rill evolution for setting and coheres thin sight
Parameter, builds the GRANULAR FLOW MODEL FOR of rock sample;
S16, otherwise, then adjusts and particle rill evolution and coheres rill evolution, and return to S13 and continue executing with.
In the specific embodiment of the method for foregoing determination rock crackle forming development degree, further, as shown in Fig. 2
The GRANULAR FLOW MODEL FOR according to structure, obtains the microfissure collection and the microfissure collection in the GRANULAR FLOW MODEL FOR destructive process
In the coordinate value (S103) of each microfissure point include:
According to the GRANULAR FLOW MODEL FOR for building, virtual uniaxial compression is carried out along long axis direction, obtain the GRANULAR FLOW MODEL FOR and break
Microfissure collection and the microfissure during bad concentrate the coordinate value of each microfissure point.
In the present embodiment, the microfissure collection is that the coordinate value of the microfissure point is two-dimensional coordinate value, the microfissure
The coordinate value of point can be expressed as X-axis coordinate and Y-axis coordinate.
In the present embodiment, in S103, the clustering method of use is different from common k-means clustering procedures and hierarchical clustering
The classics clustering method such as method.In S103, the clustering method of use is that the density based on data point is clustered, while considering
The distance of each cluster centre, realizes the packet of data point, obtains multiple microfissure groups;Wherein, the data point is microfissure collection
The coordinate value of middle microfissure point, that is to say, that in data analysis process, the coordinate value of microfissure point is considered as data point.
In the present embodiment, according to the microfissure collection that S102 is obtained, the X-axis coordinate of microfissure point and Y-axis can be sat
Mark imported into MATLAB, using the clustering method of the density based on data point, according to microfissure point aggregation extent and
The distance of each cluster centre, realizes that microfissure concentrates the packet of microfissure point.
In the present embodiment, using the clustering method of the density based on data point to the microfissure point in GRANULAR FLOW MODEL FOR
During being classified, maximum difficulty is, using a kind of suitable measurement, neighbouring microfissure point to be divided into a point set.
But, GRANULAR FLOW MODEL FOR microfissure point is generally more, for the GRANULAR FLOW MODEL FOR comprising amounts of particles close to 20,000, microfissure point
Quantity can be up to thousands of, if classifying using artificial pointwise, then can produce huge workload, and the index classified is also
It is difficult to what is determined, for closing on two microfissure points of point set simultaneously, it is difficult to which it is a certain class to be incorporated into.
In the present embodiment, the cluster centre that point denseness of set determines the point set can be directed to, so as to the phase around it will be combined
Pass o'clock is used as a single race.On this basis, just discrete point can be divided into according to the density and distance of point set
One class.This cluster mode is mainly and is applied to image processing field, and the treatment of the algorithm noise and the classification of point set are just applicable
In the classification of microfissure collection, especially in the case of number of data points is larger, the algorithm has certain robustness.The method is carried
Microfissure has been supplied to concentrate a thinking and the direction of microfissure point cluster, and cluster result can be as the number of principal component analysis
According to source.All microfissure points (microfissure collection) an as point set will be produced in GRANULAR FLOW MODEL FOR first, the point is concentrated all
The distance between point performs following steps as input quantity:
Assuming that cluster centre is surrounded by the surrounding microfissure point/data point with smaller local density, and cluster centre and
There is a certain distance between the larger point of other local densities;Two parameters are calculated for each data point i, is respectively office
Portion's density piWith the distance between the microfissure point from other with larger local density δi.The computing formula of local density is represented
For:
ρi=∑jχ(dij-dc)
Wherein, if dij-dc<0, then χ (dij-dc)=1, otherwise, χ (dij-dc)=0;dcIt is to block distance, dijIt is two data
The distance between point;This algorithm is only to ρ between different pieces of information pointiRelative size it is relevant, therefore for larger data group, point
Analysis result is relative to dcSelection be robust.
Calculate microfissure point i has the distance between the microfissure point j of larger local density δ from otheri, wherein, if micro-
Crack point i has larger local density, then δi=maxj(dij);Otherwise, apart from δiBe by microfissure point i and other have compared with
Minimum range δ between the microfissure point j of big local densityiDetermine,Therefore, cluster centre
It is considered as δiLarger point, wherein, it is smaller local density when local density is less than default Second Threshold;When local close
It is larger local density when degree is not less than default Second Threshold.
In the present embodiment, using institute in the point set (microfissure collection) the distance between a little as input quantity, by above-mentioned calculation
Method asks for multiple cluster centres, and it is a class to be incorporated into the discrete point that surrounding is closed on by cluster centre.Herein, it is main by adjusting
Section blocks the microfissure group result that the quantity of distance and cluster centre more meets the crack rule of development to obtain.It is accurate to obtain
Really fine cluster result, can select multiple cluster centres, more microfissure group be obtained, while will be far from cluster centre
Microfissure point is considered as discrete point and avoids interference result of calculation.Analyze macroscopic cracking development trend under conditions of, by multiple compared with
Small microfissure group is divided into a larger microfissure group, can so exclude the interference of some discrete microfissures, obtain compared with
It is accurate principal component analysis result.
In the present embodiment, according to the microfissure group result of S103, using PCA, each microfissure group is carried out
Analysis, obtains the axially extending direction and lateral magnification direction of each microfissure group, and each microfissure group is in both directions
Development degree.
In the present embodiment, for a better understanding of the present invention, first the key step to PCA is illustrated:
A11, solves the covariance matrix ∑ between variable:
A12, solves the eigen vector of covariance matrix ∑, tries to achieve characteristic value respectively λ1,λ2,…,λp, together
When the corresponding characteristic vector of characteristic value be:
ui=(u1i,u2i,…,upi) ', i=1,2 ... P
Characteristic vector formed orthogonal matrix be:
y1=u11x1+u21x2+…+up1xp
y2=u12x1+u22x2+…+up2xp
……
yp=u1px1+u2px2+…+uppxp
Following principal component can be so obtained, respectively with y1,y2,…,ypRepresent
Y=U 'x, x=(x1,x2,…,xp) ', x represents variable;
Being write as matrix form is
Wherein, first principal component is:
y1=u11x1+…+up1xp
The information that first principal component can be expressed is most in numerous principal components, and expressing information accounts for the ratio of gross information content
ForThe information that the bigger explanation first principal component of the ratio is included is more, can more express overall feature.Normal conditions
Two or three maximum principal components of selected characteristic value just can primary expression general characteristic down.
In the present embodiment, because the coordinate value of the microfissure point is two-dimensional coordinate value (X-axis coordinate and Y-axis coordinate), because
This variable only has two, using PCA, each microfissure group is analyzed, and determines that two principal components just can
Express the orientative feature of the microfissure collection.
In the present embodiment, after carrying out cluster analysis and principal component analysis to microfissure collection, first principal component direction is obtained
With the contribution margin on Second principal component, direction, and two principal component directions;Wherein, first principal component direction represents corresponding micro-
The axially extending direction of fracture set, Second principal component, direction represents the lateral magnification direction of corresponding microfissure group, first it is main into
Divide direction orthogonal with Second principal component, direction;The contribution margin of first principal component represents corresponding microfissure group in axially extending direction
On development degree, the contribution margin of Second principal component, represents development degree of the corresponding microfissure group on lateral magnification direction.
In the case of the contribution margin of Second principal component, is less, microfissure group can be considered as crackle, it is micro- with the gradually increase of contribution margin
The degrees of expansion of fracture set is also being stepped up, finally when the contribution margin of Second principal component, reaches to a certain degree, microfissure group
Rupture zone can be regarded as.
In the present embodiment, in order to more fully understand the present embodiment, the determination rock crackle forming hair provided the present embodiment
The method for educating degree is described in detail, and the rock sample by taking diplopore rocks sample as an example, send out by the determination rock crackle forming
The specific steps for educating the method for degree can include:
(1) indoor rock mechanics experiment is carried out to rock sample, acquisition carries out the rock sample of indoor rock mechanics experiment
Load-deformation curve, so as to provide experiment basis to set up GRANULAR FLOW MODEL FOR and the corresponding rill evolution of matching.
(2) GRANULAR FLOW MODEL FOR of rock sample is built:Granule is built first, to characterize rock sample.Granule by
Grain and cohere compositions, cohere similar between two particles attachment gum material.In granule, the rill evolution of particle is described
Have:Ball density p, smallest particles radius Rmin, maximum compare R with minimum grain sizemax/Rmin, contact Young's modulus Ec, friction factor μ,
Contact normal direction and shear stiffness compare kn/ks;The contact model for using softens contact model for displacement, and it is thin that displacement softens contact model
Sight parameter be coefficient of friction sof_fric, shear strength sof_fsmax, tensile strength sof_ftmax, compression stiffness sof_knc,
Tensible rigidity sof_knt, shearing rigidity sof_ks, remaining coefficient of friction sof_rfric, Limited Plasticity displacement sof_uplim;Retouch
Stating the rill evolution for cohering has:It is parallel to cohere radius factor λ, cohere Young's modulusCohere normal direction and shear stiffness ratioCohere normal strength average valueCohere normal strength standard deviationCohere tangential average strengthIt is glutinous
Tie tangential tension varianceUsing particle and cohere structure granule sample, the size of granule sample with carry out room
The rock sample of interior rock mechanics experiment is consistent, as shown in figure 3, carrying out the rock of the rock sample of indoor rock mechanics experiment
Bridge length 1 can be 60mm, and the Circularhole diameter 2 for carrying out the rock sample of indoor rock mechanics experiment can be 8mm, carry out interior
The rock bridge inclination angle 3 of the rock sample of rock mechanics experiment can be 90 °.Assign the particle of granule sample hypothesis and cohere thin sight
Parameter, and uniaxial compression virtual load is carried out to it, obtain the single compressing stress-strain curve of the granule sample.Will
Single compressing stress-the strain curve of the granule sample for getting, tries with the rock for carrying out indoor rock mechanics experiment
Single compressing stress-the strain curve of sample is compared.By constantly adjustment granule rill evolution (particle rill evolution and glutinous
Knot rill evolution), as shown in figure 4, by obtain the granule sample uniaxial compression load-deformation curve with carry out room
Single compressing stress-the strain curve of the rock sample of interior rock mechanics experiment is compared, when the granule examination for obtaining
The single compressing stress of the load-deformation curve of the uniaxial compression of sample and the rock sample for carrying out indoor rock mechanics experiment-should
When varied curve more coincide consistent, now granule rill evolution can be as final micro-parameter, according to the final microcosmic ginseng
Number build rock samples GRANULAR FLOW MODEL FORs, in the present embodiment, coincide refer to the granule sample uniaxial compression stress-
The difference of strain curve and the single compressing stress-strain curve of the rock sample for carrying out indoor rock mechanics experiment is default
In first threshold.
By debugging repeatedly, using the granule rill evolution described in Tables 1 and 2, single shaft fictitious compress recuperation loading, mould are carried out
Plan has obtained the single compressing stress-strain curve of granule sample, the curve and the rock for carrying out indoor rock mechanics experiment
Single compressing stress-the strain curve of sample more coincide unanimously, it is possible thereby to think that identified granule rill evolution is
Reasonably.
The granule rill evolution of table 1
The displacement of table 2 softens contact model rill evolution
Coefficient of friction | Shear strength/Pa | Tensile strength/Pa | Compression stiffness/GPa |
0 | 55.40 | 0.17 | 16.38 |
Remaining coefficient of friction | Limited Plasticity displacement | Shearing rigidity/GPa | Tensible rigidity/GPa |
0 | 14×10-6 | 7 | 15 |
(3) uniaxial compression test is carried out to the GRANULAR FLOW MODEL FOR for building.Particle body Model failure mode is illustrated in figure 5,
Wherein, 4,5 and 6 is microfissure collection, and 7 and 8 is circular hole, after GRANULAR FLOW MODEL FOR is destroyed, by particle stream software (Particle
Flow Code, PFC) the FISH orders that carry, derive the position coordinates (coordinate value) of each microfissure point.By particle stream mould
Microfissure point coordinates value file of the generation comprising distance between microfissure point in type, this document includes each microfissure
The distance between point and other microfissure points.Can be logical from MATLAB function files are write by the way of in this step
Reading microfissure point coordinates value, then the range data obtained by multiple cycle calculations are crossed, and these range data are exported as
Apart from file.
(4) cluster analysis is carried out by clustering algorithm, it is determined that blocking distance and cluster centre:
As shown in fig. 6, this cluster delimit 8 microfissure groups altogether, 9-16 is the microfissure group that cluster analysis is obtained
Numbering, but it can also be seen that the microfissure group of serial number 9,10 and 11 can be divided into a larger microfissure group, and sequence number
For 12,13 and 14 microfissure groups can be divided into another larger microfissure group, on the basis of microfissure group is merged altogether
Three larger microfissure groups can be obtained, microfissure group is redistributed into serial number 1,2 and 3.And No. 15 microfissure groups due to
Discreteness is larger, therefore does not carry out principal component analysis.
(5) the three microfissure groups obtained in (4) are carried out into principal component analysis, is illustrated in figure 7 principal component analysis principle
Schematic diagram, 17 is first principal component direction, and 18 is Second principal component, direction, and 19 is data point/microfissure point, and three are obtained respectively
(as shown in figure 8,20 is the microfissure group of serial number 1,21 is the microfissure group of serial number 2 to microfissure group, and 22 is serial number 3
Microfissure group) first principal component and Second principal component, relevant information, 3 fracture sets can distinguish in different colors, each microfissure
Group principal component direction and contribution margin are shown in Table 3.Because first principal component direction is orthogonal with Second principal component, direction, therefore in table 3
Just repeat no more Second principal component, slope and inclination angle.First principal component direction represents the axially extending side of corresponding microfissure group
To Second principal component, direction represents the lateral magnification direction of corresponding microfissure group, and especially, Second principal component, contribution margin is characterized
The degrees of expansion of corresponding microfissure group, the value is bigger, represents that the microfissure group degrees of expansion is higher.
The principal component analysis result of table 3
Embodiment two
The present invention also provides a kind of specific embodiment of the device for determining rock crackle forming development degree, because the present invention is carried
The device of the determination rock crackle forming development degree of confession with it is foregoing determination rock crackle forming development degree method specific embodiment
Corresponding, the device of the determination rock crackle forming development degree can be walked by performing the flow in above method specific embodiment
It is rapid to realize explaining in the purpose of the present invention, therefore the method specific embodiment of above-mentioned determination rock crackle forming development degree
It is bright, the specific embodiment of the device of the determination rock crackle forming development degree of present invention offer is also applied for, below the present invention
Specific embodiment in will not be described in great detail.
As shown in figure 9, the embodiment of the present invention also provides a kind of device for determining rock crackle forming development degree, including:
Build module 23, the GRANULAR FLOW MODEL FOR for building rock sample;
Acquisition module 24, for according to the GRANULAR FLOW MODEL FOR for building, obtaining the microfissure in GRANULAR FLOW MODEL FOR destructive process
Collection and the microfissure concentrate the coordinate value of each microfissure point;
Cluster module 25, the seat for concentrating each microfissure point according to the microfissure collection and the microfissure that obtain
Scale value, the density based on microfissure point is clustered, and according to the distance of each cluster centre, the fine fisssure concentrated to the microfissure
Gap point is grouped, and obtains multiple microfissure groups;
Quantization modules 26, for using PCA, analyze the axially extending direction of each microfissure group and laterally expand
Zhang Fangxiang, and quantify development degree of each microfissure group on axially extending direction and lateral magnification direction.
The device of the determination rock crackle forming development degree described in the embodiment of the present invention, by the GRANULAR FLOW MODEL FOR for building, obtains
Obtain the coordinate value that microfissure collection and the microfissure in GRANULAR FLOW MODEL FOR destructive process concentrate each microfissure point;According to obtaining
The microfissure collection and the microfissure concentrate the coordinate value of each microfissure point, the density based on microfissure point is gathered
Class, and according to the distance of each cluster centre, the microfissure point that the microfissure is concentrated is grouped, obtain multiple microfissures
Group;Using PCA, the axially extending direction and lateral magnification direction of each microfissure group are analyzed, and quantify each microfissure
Development degree of the group on axially extending direction and lateral magnification direction.
In the specific embodiment of the device of foregoing determination rock crackle forming development degree, further, the structure mould
Block includes:
Determining unit, for using particle and cohere structure granule sample, wherein, the size of the granule sample with
The size for carrying out the rock sample of indoor rock mechanics experiment is identical;
Setting unit, for setting particle rill evolution and cohering rill evolution;
Loading unit, for according to set particle rill evolution and cohere rill evolution, the granule sample is entered
Row uniaxial compression virtual load, obtains the single compressing stress-strain curve of the granule sample;
Comparing unit, for the single compressing stress-strain curve of the granule sample that will obtain with carry out interior
Single compressing stress-the strain curve of the rock sample of rock mechanics experiment is compared;
Construction unit, if for curve difference in default first threshold, according to set particle rill evolution
With cohere rill evolution, build the GRANULAR FLOW MODEL FOR of rock sample;
Returning unit, if for curve difference not in default first threshold, adjustment particle rill evolution and glutinous
Rill evolution is tied, and returns to the loading unit and continued executing with.
In the specific embodiment of the device of foregoing determination rock crackle forming development degree, further, the acquisition mould
Block, for according to the GRANULAR FLOW MODEL FOR for building, virtual uniaxial compression being carried out along long axis direction, obtains the GRANULAR FLOW MODEL FOR destruction
During microfissure collection and the microfissure concentrate the coordinate value of each microfissure point.
In the specific embodiment of the device of foregoing determination rock crackle forming development degree, further, the cluster mould
Block includes:
First computing unit, the coordinate value for concentrating each microfissure point according to the microfissure obtains microfissure point i
With the distance between other microfissure points j dij;
Second computing unit, is surrounded, and gather for setting cluster centre by the surrounding's microfissure point with smaller local density
There is a certain distance between the larger microfissure point in class center and other local densities, calculate the local density of microfissure point i
ρi:
ρi=∑jχ(dij-dc);
3rd computing unit, for calculate microfissure point i have from other between microfissure point j of larger local density away from
From δi, wherein, if microfissure point i has larger local density, δi=maxj(dij);Otherwise, apart from δiIt is by microfissure point
I and other minimum range δ between there is the microfissure point j of larger local densityiDetermine,
Wherein, if dij-dc<0, then χ (dij-dc)=1, otherwise, χ (dij-dc)=0;dcIt is to block distance, dijIt is microfissure
The distance between point i and other microfissure points j;It is smaller local density when local density is less than default Second Threshold;When
Local density is larger local density when being not less than default Second Threshold.
In the specific embodiment of the device of foregoing determination rock crackle forming development degree, further, the quantization mould
Block, for the method using principal component analysis, is analyzed to each microfissure group, obtain first principal component direction, second it is main into
Divide the contribution margin in direction, the contribution margin of first principal component and Second principal component,;
Wherein, the first principal component direction represents the axially extending direction of corresponding microfissure group, Second principal component, side
To the lateral magnification direction for representing corresponding microfissure group, first principal component direction is orthogonal with Second principal component, direction;First master
The contribution margin of composition represents development degree of the corresponding microfissure group on axially extending direction, the contribution margin table of Second principal component,
Show development degree of the corresponding microfissure group on lateral magnification direction.
The above is the preferred embodiment of the present invention, it is noted that for those skilled in the art
For, on the premise of principle of the present invention is not departed from, some improvements and modifications can also be made, these improvements and modifications
Should be regarded as protection scope of the present invention.
Claims (10)
1. it is a kind of determine rock crackle forming development degree method, it is characterised in that including:
Build the GRANULAR FLOW MODEL FOR of rock sample;
According to the GRANULAR FLOW MODEL FOR for building, obtain microfissure collection and the microfissure in GRANULAR FLOW MODEL FOR destructive process and concentrate every
The coordinate value of individual microfissure point;
The coordinate value of each microfissure point is concentrated according to the microfissure collection for obtaining and the microfissure, based on microfissure point
Density is clustered, and according to the distance of each cluster centre, the microfissure point that the microfissure is concentrated is grouped, and obtains many
Individual microfissure group;
Using PCA, the axially extending direction and lateral magnification direction of each microfissure group are analyzed, and quantify each fine fisssure
Development degree of the gap group on axially extending direction and lateral magnification direction.
2. the method for determining rock crackle forming development degree according to claim 1, it is characterised in that structure rock examination
The GRANULAR FLOW MODEL FOR of sample includes:
S11, using particle and cohere structure granule sample, wherein, the size of the granule sample with carry out indoor rock
The size of the rock sample of mechanical test is identical;
S12, sets particle rill evolution and coheres rill evolution;
S13, according to set particle rill evolution and cohere rill evolution, it is virtual to carry out uniaxial compression to the granule sample
Loading, obtains the single compressing stress-strain curve of the granule sample;
S14, the single compressing stress-strain curve of the granule sample that will be obtained with carry out indoor rock mechanics experiment
Single compressing stress-the strain curve of rock sample is compared;
S15, if the difference of curve is in default first threshold, joins according to the particle rill evolution for setting with thin sight is cohered
Number, builds the GRANULAR FLOW MODEL FOR of rock sample;
S16, otherwise, then adjusts and particle rill evolution and coheres rill evolution, and return to S13 and continue executing with.
3. it is according to claim 1 determine rock crackle forming development degree method, it is characterised in that it is described according to build
GRANULAR FLOW MODEL FOR, the microfissure collection and the microfissure obtained in the GRANULAR FLOW MODEL FOR destructive process concentrates each microfissure point
Coordinate value include:
According to the GRANULAR FLOW MODEL FOR for building, virtual uniaxial compression is carried out along long axis direction, obtain the GRANULAR FLOW MODEL FOR and destroyed
Microfissure collection and the microfissure in journey concentrate the coordinate value of each microfissure point.
4. it is according to claim 1 determine rock crackle forming development degree method, it is characterised in that what the basis was obtained
The microfissure collection and the microfissure concentrate the coordinate value of each microfissure point, and the density based on microfissure point is clustered,
And according to the distance of each cluster centre, the microfissure point that the microfissure is concentrated is grouped, obtain multiple microfissure group bags
Include:
The coordinate value of each microfissure point is concentrated according to the microfissure, is obtained between microfissure point i and other microfissure points j
Apart from dij;
If cluster centre is surrounded by the surrounding microfissure point with smaller local density, and cluster centre and other local densities compared with
There is a certain distance between big microfissure point, calculate the local density ρ of microfissure point ii:
ρi=∑jχ(dij-dc);
Calculate microfissure point i has the distance between the microfissure point j of larger local density δ from otheri, wherein, if microfissure point
I has larger local density, then δi=maxj(dij);Otherwise, apart from δiIt is that there is larger part with other by microfissure point i
Minimum range δ between the microfissure point j of densityiDetermine,
Wherein, if dij-dc<0, then χ (dij-dc)=1, otherwise, χ (dij-dc)=0;dcIt is to block distance, dijIt is microfissure point i
The distance between with other microfissure points j;It is smaller local density when local density is less than default Second Threshold;Work as part
Density is larger local density when being not less than default Second Threshold.
5. the method for the determination rock crackle forming development degree according to claim 1 or 4, it is characterised in that described using master
Componential analysis, analyze the axially extending direction and lateral magnification direction of each microfissure group, and quantify each microfissure group in axial direction
Development degree on bearing of trend and lateral magnification direction includes:
Using the method for principal component analysis, each microfissure group is analyzed, obtains first principal component direction, Second principal component, side
To, the contribution margin of first principal component and the contribution margin of Second principal component,;
Wherein, the first principal component direction represents the axially extending direction of corresponding microfissure group, Second principal component, direction table
Show the lateral magnification direction of corresponding microfissure group, first principal component direction is orthogonal with Second principal component, direction;First principal component
Contribution margin represent development degree of the corresponding microfissure group on axially extending direction, the contribution margin of Second principal component, represents phase
Development degree of the microfissure group answered on lateral magnification direction.
6. it is a kind of determine rock crackle forming development degree device, it is characterised in that including:
Build module, the GRANULAR FLOW MODEL FOR for building rock sample;
Acquisition module, for according to the GRANULAR FLOW MODEL FOR for building, obtaining the microfissure collection in GRANULAR FLOW MODEL FOR destructive process and institute
State the coordinate value that microfissure concentrates each microfissure point;
Cluster module, the coordinate value for concentrating each microfissure point according to the microfissure collection and the microfissure that obtain,
Density based on microfissure point is clustered, and according to the distance of each cluster centre, the microfissure point concentrated to the microfissure
It is grouped, is obtained multiple microfissure groups;
Quantization modules, for using PCA, analyze the axially extending direction and lateral magnification direction of each microfissure group,
And quantify development degree of each microfissure group on axially extending direction and lateral magnification direction.
7. it is according to claim 6 determine rock crackle forming development degree device, it is characterised in that the structure module bag
Include:
Determining unit, for using particle and cohere structure granule sample, wherein, the size of the granule sample with carry out
The size of the rock sample of indoor rock mechanics experiment is identical;
Setting unit, for setting particle rill evolution and cohering rill evolution;
Loading unit, for according to set particle rill evolution and cohere rill evolution, list is carried out to the granule sample
Axle compresses virtual load, obtains the single compressing stress-strain curve of the granule sample;
Comparing unit, for the single compressing stress-strain curve of the granule sample that will obtain with carry out indoor rock
Single compressing stress-the strain curve of the rock sample of mechanical test is compared;
Construction unit, if for curve difference in default first threshold, according to the particle rill evolution for setting and glutinous
Knot rill evolution, builds the GRANULAR FLOW MODEL FOR of rock sample;
Returning unit, if for curve difference not in default first threshold, adjustment particle rill evolution and cohere thin
Parameter is seen, and returns to the loading unit and continued executing with.
8. it is according to claim 6 determine rock crackle forming development degree device, it is characterised in that the acquisition module,
For according to the GRANULAR FLOW MODEL FOR for building, virtual uniaxial compression being carried out along long axis direction, obtain the GRANULAR FLOW MODEL FOR and destroyed
Microfissure collection and the microfissure in journey concentrate the coordinate value of each microfissure point.
9. it is according to claim 6 determine rock crackle forming development degree device, it is characterised in that the cluster module bag
Include:
First computing unit, the coordinate value for concentrating each microfissure point according to the microfissure, obtain microfissure point i and its
The distance between his microfissure point j dij;
Second computing unit, is surrounded for setting cluster centre by the surrounding's microfissure point with smaller local density, and in cluster
There is a certain distance between the larger microfissure point of the heart and other local densities, calculate the local density ρ of microfissure point ii:
ρi=∑jχ(dij-dc);
3rd computing unit, has the distance between microfissure point j of larger local density for calculating microfissure point i from other
δi, wherein, if microfissure point i has larger local density, δi=maxj(dij);Otherwise, apart from δiIt is by microfissure point i
And other minimum range δ between there is the microfissure point j of larger local densityiDetermine,
Wherein, if dij-dc<0, then χ (dij-dc)=1, otherwise, χ (dij-dc)=0;dcIt is to block distance, dijIt is microfissure point i
The distance between with other microfissure points j;It is smaller local density when local density is less than default Second Threshold;Work as part
Density is larger local density when being not less than default Second Threshold.
10. the device of the determination rock crackle forming development degree according to claim 6 or 9, it is characterised in that the quantization mould
Block, for the method using principal component analysis, is analyzed to each microfissure group, obtain first principal component direction, second it is main into
Divide the contribution margin in direction, the contribution margin of first principal component and Second principal component,;
Wherein, the first principal component direction represents the axially extending direction of corresponding microfissure group, Second principal component, direction table
Show the lateral magnification direction of corresponding microfissure group, first principal component direction is orthogonal with Second principal component, direction;First principal component
Contribution margin represent development degree of the corresponding microfissure group on axially extending direction, the contribution margin of Second principal component, represents phase
Development degree of the microfissure group answered on lateral magnification direction.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710017322.0A CN106706884B (en) | 2017-01-11 | 2017-01-11 | A kind of method and device of determining rock crackle forming development degree |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710017322.0A CN106706884B (en) | 2017-01-11 | 2017-01-11 | A kind of method and device of determining rock crackle forming development degree |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106706884A true CN106706884A (en) | 2017-05-24 |
CN106706884B CN106706884B (en) | 2018-12-25 |
Family
ID=58908174
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710017322.0A Active CN106706884B (en) | 2017-01-11 | 2017-01-11 | A kind of method and device of determining rock crackle forming development degree |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106706884B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109085321A (en) * | 2018-07-25 | 2018-12-25 | 石家庄铁道大学 | The scaling method and terminal device of Micromechanical Parameters of Rock |
CN113177302A (en) * | 2021-04-14 | 2021-07-27 | 北京科技大学 | Fractured rock particle flow model construction and new crack occurrence analysis method |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090132218A1 (en) * | 2006-11-29 | 2009-05-21 | Ledgerwood Iii Leroy W | Discrete element modeling of rock destruction under high pressure conditions |
CN101915764A (en) * | 2010-08-10 | 2010-12-15 | 武汉武大卓越科技有限责任公司 | Road surface crack detection method based on dynamic programming |
CN102621009A (en) * | 2012-03-21 | 2012-08-01 | 武汉大学 | Test method for simulating long-term deformation of rockfill |
CN103940666A (en) * | 2014-03-18 | 2014-07-23 | 中国矿业大学 | Determination method for mesoscopic parameters simulating mechanical properties of intermittent crack rock |
CN105093264A (en) * | 2014-05-12 | 2015-11-25 | 中国石油化工股份有限公司 | Rock stratum effective fracture-developed zone predicting method and system |
CN105844032A (en) * | 2016-03-30 | 2016-08-10 | 李新卫 | Particle flow based irregular particle numerical simulation research method for microwave induced damage |
-
2017
- 2017-01-11 CN CN201710017322.0A patent/CN106706884B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090132218A1 (en) * | 2006-11-29 | 2009-05-21 | Ledgerwood Iii Leroy W | Discrete element modeling of rock destruction under high pressure conditions |
CN101915764A (en) * | 2010-08-10 | 2010-12-15 | 武汉武大卓越科技有限责任公司 | Road surface crack detection method based on dynamic programming |
CN102621009A (en) * | 2012-03-21 | 2012-08-01 | 武汉大学 | Test method for simulating long-term deformation of rockfill |
CN103940666A (en) * | 2014-03-18 | 2014-07-23 | 中国矿业大学 | Determination method for mesoscopic parameters simulating mechanical properties of intermittent crack rock |
CN105093264A (en) * | 2014-05-12 | 2015-11-25 | 中国石油化工股份有限公司 | Rock stratum effective fracture-developed zone predicting method and system |
CN105844032A (en) * | 2016-03-30 | 2016-08-10 | 李新卫 | Particle flow based irregular particle numerical simulation research method for microwave induced damage |
Non-Patent Citations (3)
Title |
---|
吴顺川: "基于颗粒元模拟的含软弱夹层类土质边坡变形破坏过程分析", 《岩土力学》 * |
杨国勇: "基于层次分析-模糊聚类分析法的导水裂隙带发育高度研究", 《采矿与安全工程学报》 * |
杨圣奇: "断续双裂隙岩石抗拉强度与裂纹扩展颗粒流分析", 《中国矿业大学学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109085321A (en) * | 2018-07-25 | 2018-12-25 | 石家庄铁道大学 | The scaling method and terminal device of Micromechanical Parameters of Rock |
CN109085321B (en) * | 2018-07-25 | 2021-06-01 | 石家庄铁道大学 | Calibration method of rock mesomechanics parameters and terminal equipment |
CN113177302A (en) * | 2021-04-14 | 2021-07-27 | 北京科技大学 | Fractured rock particle flow model construction and new crack occurrence analysis method |
CN113177302B (en) * | 2021-04-14 | 2023-07-14 | 北京科技大学 | Crack rock particle flow model construction and new crack occurrence analysis method |
Also Published As
Publication number | Publication date |
---|---|
CN106706884B (en) | 2018-12-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Pathan et al. | Predictions of the mechanical properties of unidirectional fibre composites by supervised machine learning | |
Tran et al. | A benchmark of batch-effect correction methods for single-cell RNA sequencing data | |
Lin et al. | Experimental design and experimental inference in stone artifact archaeology | |
Huang et al. | Nanoscience and technology publications and patents: a review of social science studies and search strategies | |
López‐Angulo et al. | The role of root community attributes in predicting soil fungal and bacterial community patterns | |
Ioviţă | Comparing stone tool resharpening trajectories with the aid of elliptical Fourier analysis | |
US10373708B2 (en) | Systems and methods for generating biomarker signatures with integrated dual ensemble and generalized simulated annealing techniques | |
CN1741036A (en) | Heuristic method of classification belongs to different state bio-samples model | |
Das | Pattern Recognition using the Fuzzy c-means Technique | |
CN103116713A (en) | Method of predicting interaction between chemical compounds and proteins based on random forest | |
WO2017194431A1 (en) | System for predicting efficacy of a target-directed drug to treat a disease | |
CN105046080B (en) | A kind of rock mass quality classification | |
Clarkson | Regional diversity within the core technology of the Howiesons Poort techno-complex | |
CN106706884A (en) | Method and apparatus for determining development degree of rock cracks | |
CN107729555A (en) | A kind of magnanimity big data Distributed Predictive method and system | |
Sholts et al. | Tracing social interactions in Pleistocene North America via 3D model analysis of stone tool asymmetry | |
Villa et al. | Taxonomic and functional beta diversity of woody communities along Amazon forest succession: The relative importance of stand age, soil properties and spatial factor | |
CN103898890B (en) | Soil layer quantization layering method based on double-bridge static sounding data of BP neural network | |
Melendez-Melendez et al. | An improved algorithm for partial clustering | |
Massalov et al. | Application of soft computing techniques to estimate cutter life index using mechanical properties of rocks | |
Wu et al. | DeepRetention: a deep learning approach for intron retention detection | |
Garg et al. | A new simulation approach of genetic programming in modelling of soil water retention property of unsaturated soil | |
Chowdhury | Effective clustering of scRNA-seq data to identify biomarkers without user input | |
EP2948885B1 (en) | Methods, systems, and computer readable media for data analysis and inference of particle diffusion in target materials and target material simulants | |
CN103970973B (en) | A kind of simulation RCT analysis methods based on True Data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
CB03 | Change of inventor or designer information |
Inventor after: Zhou Yu Inventor after: Liu Bing Inventor after: Wang Xue Inventor after: Wu Shunchuan Inventor after: Zhou Jianxin Inventor before: Liu Bing Inventor before: Zhou Yu Inventor before: Wang Xue Inventor before: Wu Shunchuan Inventor before: Zhou Jianxin |
|
CB03 | Change of inventor or designer information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |