CN106706884B - A kind of method and device of determining rock crackle forming development degree - Google Patents
A kind of method and device of determining rock crackle forming development degree Download PDFInfo
- Publication number
- CN106706884B CN106706884B CN201710017322.0A CN201710017322A CN106706884B CN 106706884 B CN106706884 B CN 106706884B CN 201710017322 A CN201710017322 A CN 201710017322A CN 106706884 B CN106706884 B CN 106706884B
- Authority
- CN
- China
- Prior art keywords
- microfissure
- point
- principal component
- group
- rock
- 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.)
- Active
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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Engineering & Computer Science (AREA)
- Pathology (AREA)
- Immunology (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Medicinal Chemistry (AREA)
- Food Science & Technology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention provides a kind of method and device of determining rock crackle forming development degree, can quantify development degree of each microfissure group on axially extending direction and lateral magnification direction.The described method includes: the GRANULAR FLOW MODEL FOR of building rock sample;According to the GRANULAR FLOW MODEL FOR of building, the coordinate value that microfissure collection and the microfissure in GRANULAR FLOW MODEL FOR destructive process concentrate each microfissure point is obtained;The coordinate value of each microfissure point is concentrated according to the obtained microfissure collection and the microfissure, density based on microfissure point is clustered, and according to the distance of each cluster centre, the microfissure point concentrated to the microfissure is grouped, and obtains multiple microfissure groups;Using Principal Component Analysis, the axially extending direction and lateral magnification direction of each microfissure group are analyzed, and quantifies development degree of each microfissure group on axially extending direction and lateral magnification direction.The present invention is suitable for During Geotechnical Tests and field of engineering technology.
Description
Technical field
The present invention relates to geomaterial Research on Mechanical Properties fields, particularly relate to a kind of determining rock crackle forming development journey
The method and device of degree.
Background technique
In recent years, more research is the rock sample progress failure test for prefabricated crack, probes into and deposits in crack
In case, in rock failure process crack development trend.But be difficult to observe due to microfissure naked eyes, research is most
The information such as the mechanics parameter of rock and the development trend in macroscopical crack are obtained, therefore universal not to research in terms of microcosmic mechanism
Foot.Meanwhile also having relevant research is to combine numerical simulation software, especially discrete meta software carries out the simulation of traditional experiment,
Disclose the thin sight mechanism in rock failure process.Strong supplement of the numerical simulation software as traditional experiment, can effectively catch
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 uniaxial compression becomes
Gesture, but since Cementation failure is regarded as microfissure in simulation process by PFC software, the result finally presented mostly can only
Judge the development trend in macroscopical crack roughly by the aggregation of microfissure, rather than accurate quantification macroscopical crack development is described
Trend and level of breakage.
Summary of the invention
The technical problem to be solved in the present invention is to provide a kind of method and devices of determining rock crackle forming development degree, with solution
The problem of describing macroscopical crack development trend and level of breakage while accurate quantification is certainly unable to present in the prior art.
In order to solve the above technical problems, the embodiment of the present invention provides a kind of method of determining rock crackle forming development degree, packet
It includes:
Construct the GRANULAR FLOW MODEL FOR of rock sample;
According to the GRANULAR FLOW MODEL FOR of building, 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 obtained microfissure collection and the microfissure, is based on microfissure
The density of point is clustered, and according to the distance of each cluster centre, the microfissure point concentrated to the microfissure is grouped, obtains
To multiple microfissure groups;
Using Principal Component Analysis, the axially extending direction and lateral magnification direction of each microfissure group are analyzed, and is quantified each
Development degree of the microfissure group on axially extending direction and lateral magnification direction.
Further, the GRANULAR FLOW MODEL FOR of the building rock sample includes:
S11 using particle and coheres building granule sample, wherein the size of the granule sample and progress are indoor
The size of the rock sample of rock mechanics experiment is identical;
S12 is arranged particle rill evolution and coheres rill evolution;
S13 according to the particle rill evolution of setting and coheres rill evolution, carries out uniaxial compression to the granule sample
Virtual load obtains single compressing stress-strain curve of the granule sample;
S14 by single compressing stress-strain curve of the obtained granule sample and carries out indoor rock mechanics examination
Single compressing stress-the strain curve for the rock sample tested is compared;
S15, if the difference of curve in preset first threshold, according to the particle rill evolution of setting and coheres thin sight
Parameter constructs the GRANULAR FLOW MODEL FOR of rock sample;
Otherwise S16 then adjusts particle rill evolution and coheres rill evolution, and return to S13 and continue to execute.
Further, the GRANULAR FLOW MODEL FOR according to building obtains 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 of building, virtual uniaxial compression is carried out along long axis direction, it is broken to obtain the GRANULAR FLOW MODEL FOR
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 obtains concentrate the coordinate of each microfissure point
Value, the density based on microfissure point are 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 that each microfissure point is concentrated according to the microfissure, obtain microfissure point i and other microfissure points j it
Between distance dij;
If cluster centre is surrounded by surrounding's microfissure point with smaller local density, and cluster centre and other parts are close
It spends between biggish microfissure point there are a certain distance, calculates the local density ρ of microfissure point ii:
ρi=∑jχ(dij-dc);
Microfissure point i is calculated from the distance between other microfissure point j with larger local density δiIf, wherein it is micro-
Crack point i has larger local density, then δi=maxj(dij);Otherwise, distance δiBe by microfissure point i and other have compared with
Minimum range δ between the microfissure point j of big local densityiDetermining,
Wherein, if dij-dc< 0, then χ (dij-dc)=1, otherwise, χ (dij-dc)=0;dcIt is truncation 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 preset second threshold;When
Local density is larger local density when being not less than preset second threshold.
Further, described to use Principal Component Analysis, 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 and include:
Using the method for principal component analysis, each microfissure group is analyzed, obtain first principal component direction, second it is main at
Divide the contribution margin in direction, the contribution margin of first principal component and Second principal component,;
Wherein, the first principal component direction indicates the axially extending direction of corresponding microfissure group, Second principal component, side
To the lateral magnification direction for indicating corresponding microfissure group, first principal component direction is orthogonal with Second principal component, direction;First is main
The contribution margin of ingredient indicates 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 of determining rock crackle forming development degree characterized by comprising
Module is constructed, for constructing the GRANULAR FLOW MODEL FOR of rock sample;
Module is obtained, for the GRANULAR FLOW MODEL FOR according to building, obtains the microfissure collection in GRANULAR FLOW MODEL FOR destructive process
And the microfissure concentrates the coordinate value of each microfissure point;
Cluster module, for concentrating the coordinate of each microfissure point according to obtained the microfissure collection and the microfissure
Value, the density based on microfissure point are 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 analyze axially extending direction and the lateral magnification of each microfissure group for using Principal Component Analysis
Direction, and quantify development degree of each microfissure group on axially extending direction and lateral magnification direction.
Further, the building module includes:
Determination unit, for using particle and cohere building 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 particle rill evolution to be arranged and coheres rill evolution;
Loading unit, for according to the particle rill evolution of setting and cohering rill evolution, to the granule sample into
Row uniaxial compression virtual load obtains single compressing stress-strain curve of the granule sample;
Comparing unit, for single compressing stress-strain curve of the obtained granule sample and progress is indoor
Single compressing stress-strain curve of the rock sample of rock mechanics experiment is compared;
Construction unit, if the difference for curve in preset first threshold, according to the particle rill evolution of setting
With cohere rill evolution, construct the GRANULAR FLOW MODEL FOR of rock sample;
Return unit, if the difference for curve not in preset first threshold, adjusts particle rill evolution and sticks
Rill evolution is tied, and returns to the loading unit and continues to execute.
Further, the acquisition module carries out virtual uniaxial for the GRANULAR FLOW MODEL FOR according to building 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 obtains microfissure point i for concentrating the coordinate value of each microfissure point according to the microfissure
The distance between other microfissure points j dij;
Second computing unit is surrounded by surrounding's microfissure point with smaller local density for setting cluster centre, and gathers
There are a certain distance between the biggish microfissure point in class center and other local densities, calculate the local density of microfissure point i
ρi:
ρi=∑jχ(dij-dc);
Third computing unit, for calculating microfissure point i between other microfissure point j with larger local density
Distance δi, wherein if microfissure point i has larger local density, δi=maxj(dij);Otherwise, distance δiIt is to pass through microfissure
Minimum range δ between point i and other microfissure point j with larger local densityiDetermining,
Wherein, if dij-dc< 0, then χ (dij-dc)=1, otherwise, χ (dij-dc)=0;dcIt is truncation 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 preset second threshold;When
Local density is larger local density when being not less than preset second threshold.
Further, the quantization modules analyze each microfissure group for the method using principal component analysis,
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 indicates the axially extending direction of corresponding microfissure group, Second principal component, side
To the lateral magnification direction for indicating corresponding microfissure group, first principal component direction is orthogonal with Second principal component, direction;First is main
The contribution margin of ingredient indicates 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 advantageous effects of the above technical solutions of the present invention are as follows:
In above scheme, by the GRANULAR FLOW MODEL FOR of building, obtain microfissure collection in GRANULAR FLOW MODEL FOR destructive process and
The microfissure concentrates the coordinate value of each microfissure point;It is concentrated according to the obtained microfissure collection and the microfissure each
The coordinate value of microfissure point, the density based on microfissure point are 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 Principal Component Analysis, the axis of each microfissure group is analyzed
To extending direction and lateral magnification direction, and quantify development of each microfissure group on axially extending direction and lateral magnification direction
Degree.
Detailed description of the invention
Fig. 1 is the flow diagram of the method for determining rock crackle forming development degree provided in an embodiment of the present invention;
Fig. 2 is data analysis process schematic diagram 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 schematic illustration provided in an embodiment of the present invention;
Fig. 8 is clustering provided in an embodiment of the present invention and principal component analysis result schematic diagram;
Fig. 9 is the structural schematic diagram of the device of determining rock crackle forming development degree provided in an embodiment of the present invention.
Specific embodiment
To keep the technical problem to be solved in the present invention, technical solution and advantage clearer, below in conjunction with attached drawing and tool
Body embodiment is described in detail.
The present invention it is existing macroscopical crack development trend and level of breakage are described with being unable to accurate quantification aiming at the problem that, mention
For a kind of method and device of determining rock crackle forming development degree.
Embodiment one
Referring to shown in Fig. 1, a kind of method of determining rock crackle forming development degree provided in an embodiment of the present invention, comprising:
S101 constructs the GRANULAR FLOW MODEL FOR of rock sample;
S102 obtains microfissure collection in GRANULAR FLOW MODEL FOR destructive process and described micro- according to the GRANULAR FLOW MODEL FOR of building
Concentrate the coordinate value of each microfissure point in crack;
S103 concentrates the coordinate value of each microfissure point according to the obtained microfissure collection and the microfissure, is based on
The density of microfissure point is clustered, and according to the distance of each cluster centre, is carried out to the microfissure point that the microfissure is concentrated
Grouping, obtains multiple microfissure groups;
S104 analyzes the axially extending direction and lateral magnification direction of each microfissure group using Principal Component Analysis, and measures
Change development degree of each microfissure group on axially extending direction and lateral magnification direction.
The method of determination rock crackle forming development degree described in the embodiment of the present invention is obtained by the GRANULAR FLOW MODEL FOR of building
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 concentrated to the microfissure is grouped, and obtains multiple microfissures
Group;Using Principal Component Analysis, the axially extending direction and lateral magnification direction of each microfissure group are analyzed, and quantifies each microfissure
Development degree of the group on axially extending direction and lateral magnification direction.
In the specific embodiment of the method for aforementioned determining rock crackle forming development degree, further, the building rock
The GRANULAR FLOW MODEL FOR (S101) of stone sample includes:
S11 using particle and coheres building granule sample, wherein the size of the granule sample and progress are indoor
The size of the rock sample of rock mechanics experiment is identical;
S12 is arranged particle rill evolution and coheres rill evolution;
S13 according to the particle rill evolution of setting and coheres rill evolution, carries out uniaxial compression to the granule sample
Virtual load obtains single compressing stress-strain curve of the granule sample;
S14 by single compressing stress-strain curve of the obtained granule sample and carries out indoor rock mechanics examination
Single compressing stress-the strain curve for the rock sample tested is compared;
S15, if the difference of curve in preset first threshold, according to the particle rill evolution of setting and coheres thin sight
Parameter constructs the GRANULAR FLOW MODEL FOR of rock sample;
Otherwise S16 then adjusts particle rill evolution and coheres rill evolution, and return to S13 and continue to execute.
In the specific embodiment of the method for aforementioned determining rock crackle forming development degree, further, as shown in Fig. 2,
The GRANULAR FLOW MODEL FOR according to building obtains 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 of building, virtual uniaxial compression is carried out along long axis direction, it is broken to obtain the GRANULAR FLOW MODEL FOR
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 procedure 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, and is considered simultaneously
The distance of each cluster centre realizes the grouping 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, the microfissure collection obtained according to S102 can be sat the X axis coordinate of microfissure point and Y-axis
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 grouping 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 that neighbouring microfissure point is divided into a point set using a kind of suitable measurement.
But GRANULAR FLOW MODEL FOR microfissure point is usually 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 can then generate huge workload, and the index classified is also using artificial point-by-point classification
It is difficult to determination, for closing on the microfissure point of two point sets simultaneously, it is difficult to be incorporated into as certain one kind.
In the present embodiment, the cluster centre of the point set can be determined for a denseness of set, to will combine surrounding phase
It closes point and is used as an individual race.On this basis, discrete point can be divided into according to the density and distance of point set
It is a kind of.This cluster mode is mainly applied to field of image processing, and the processing of the algorithm noise and the classification of point set are just applicable in
In the classification of microfissure collection, especially in the larger situation of number of data points, which has certain robustness.This method mentions
Microfissure has been supplied to concentrate a thinking and the direction of microfissure point cluster, and cluster result can be used as the number of principal component analysis
According to source.All microfissure points (microfissure collection) will be generated in GRANULAR FLOW MODEL FOR first as a point set, which are concentrated all
The distance between point is used as input quantity, and executes following steps:
Assuming that cluster centre by with smaller local density surrounding's microfissure point/data point surround, and cluster centre and
There are a certain distance between the biggish point of other local densities;Two parameters are calculated for each data point i, are office respectively
Portion's density piWith from the distance between other microfissure points with larger local density δi.The calculation formula of local density indicates
Are as follows:
ρi=∑jχ(dij-dc)
Wherein, if dij-dc< 0, then χ (dij-dc)=1, otherwise, χ (dij-dc)=0;dcIt is truncation distance, dijIt is two data
The distance between point;This algorithm is only to ρ between different data pointiRelative size it is related, therefore for larger data group, point
Result is analysed relative to dcSelection be robust.
Microfissure point i is calculated from the distance between other microfissure point j with larger local density δiIf, wherein it is micro-
Crack point i has larger local density, then δi=maxj(dij);Otherwise, distance δiBe by microfissure point i and other have compared with
Minimum range δ between the microfissure point j of big local densityiDetermining,Therefore, cluster centre
It is considered as δiBiggish point, wherein be smaller local density when local density is less than preset second threshold;When part is close
It is larger local density when degree is not less than preset second threshold.
In the present embodiment, it regard the distance between all the points in the point set (microfissure collection) as input quantity, passes through above-mentioned calculation
Method seeks multiple cluster centres, is incorporated by the discrete point that cluster centre closes on surrounding as one kind.Herein, mainly pass through tune
The quantity of distance and cluster centre is truncated to obtain the microfissure group result for more meeting the crack rule of development in section.It is quasi- to obtain
Really fine cluster result, can select multiple cluster centres, obtain more microfissure group, while will be far from cluster centre
Microfissure point is considered as discrete point and avoids interference calculation result.Under conditions of analyzing macroscopic cracking development trend, by it is multiple compared with
Small microfissure group is divided into a biggish microfissure group, can exclude the interference of some discrete microfissures in this way, obtain compared with
For accurate principal component analysis result.
In the present embodiment, each microfissure group is carried out using Principal Component Analysis according to the microfissure group result of S103
Analysis, the axially extending direction for obtaining each microfissure group and lateral magnification direction and each microfissure group are in both directions
Development degree.
In the present embodiment, for a better understanding of the present invention, first the key step of Principal Component Analysis is illustrated:
A11 solves the covariance matrix ∑ between variable:
A12 solves the eigen vector of covariance matrix ∑, and acquiring characteristic value is respectively λ1,λ2,…,λp, together
When the corresponding feature vector of characteristic value are as follows:
ui=(u1i,u2i,…,upi) ', i=1,2 ... P
The orthogonal matrix that feature vector is formed are as follows:
y1=u11x1+u21x2+…+up1xp
y2=u12x1+u22x2+…+up2xp
……
yp=u1px1+u2px2+…+uppxp
Following principal component available in this way, respectively with y1,y2,…,ypIt indicates
Y=U 'x, x=(x1,x2,…,xp) ', x indicates variable;
Being write as matrix form is
Wherein, first principal component are as follows:
y1=u11x1+…+up1xp
The information that first principal component can express be in numerous principal components at most, expressing information accounts for the ratio of gross information content
ForThe ratio is bigger, and the information for illustrating that first principal component is included is more, can more express overall feature.Usual situation
Two or three maximum principal components of selected characteristic value just can primary expression general characteristic down.
In the present embodiment, since the coordinate value of the microfissure point is two-dimensional coordinate value (X axis coordinate and Y axis coordinate), because
There are two this variables, using Principal Component Analysis, analyzes each microfissure group, determines that two principal components just can
Express the orientative feature of the microfissure collection.
In the present embodiment, after carrying out clustering 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 indicates corresponding micro-
The axially extending direction of fracture set, Second principal component, direction indicate the lateral magnification direction of corresponding microfissure group, first it is main at
Divide direction orthogonal with Second principal component, direction;The contribution margin of first principal component indicates corresponding microfissure group in axially extending direction
On development degree, the contribution margin of Second principal component, indicates development degree of the corresponding microfissure group on lateral magnification direction.
In the lesser situation of contribution margin of Second principal component, microfissure group can be considered as crackle, it is micro- with being gradually increased for contribution margin
The degrees of expansion of fracture set is also being stepped up, finally when the contribution margin of Second principal component, reaches a certain level, microfissure group
It can be regarded as rupture zone.
In the present embodiment, the present embodiment, sends out the determining rock crackle forming provided in this embodiment in order to better understand
The method for educating degree is described in detail, and the rock sample is by taking diplopore rocks sample as an example, the determining rock crackle forming hair
The specific steps for educating the method for degree may include:
(1) indoor rock mechanics experiment is carried out to rock sample, obtains the rock sample for carrying out indoor rock mechanics experiment
Load-deformation curve, to provide experiment basis to establish GRANULAR FLOW MODEL FOR and matching corresponding rill evolution.
(2) GRANULAR FLOW MODEL FOR of rock sample is constructed: building granule first, to characterize rock sample.Granule by
Grain and cohere composition, cohere similar to the gum material adhered between two particles.In granule, the rill evolution of particle is described
Have: ball density p, smallest particles radius Rmin, it is maximum with minimum grain size ratio Rmax/Rmin, contact Young's modulus Ec, friction factor μ,
Contact normal direction and shear stiffness ratio kn/ks;The contact model used softens contact model for displacement, and displacement softening contact model is thin
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 are displaced sof_uplim;It retouches
Stating the rill evolution cohered has: cohering radius factor λ in parallel, coheres Young's modulusCohere normal direction and shear stiffness ratioCohere normal strength average valueCohere normal strength standard deviationCohere tangential average strength
Cohere tangential tension varianceUsing particle and cohere building granule sample, the size and progress of granule sample
The rock sample of indoor rock mechanics experiment is consistent, as shown in figure 3, carry out the rock sample of indoor rock mechanics experiment
Rock 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 room
The rock bridge inclination angle 3 of the rock sample of interior rock mechanics experiment can be 90 °.Assign the particle and cohere thin that granule sample assumes
Parameter is seen, and uniaxial compression virtual load is carried out to it, obtains single compressing stress-strain curve of the granule sample.
Single compressing stress-the strain curve for the granule sample that will acquire, with the rock for carrying out indoor rock mechanics experiment
Single compressing stress-strain curve of sample is compared.By constantly adjust granule rill evolution (particle rill evolution and
Cohere rill evolution), as shown in figure 4, by the load-deformation curve of the uniaxial compression of the granule sample of acquisition and progress
Single compressing stress-strain curve of the rock sample of indoor rock mechanics experiment is compared, when the granule of acquisition
The single compressing stress-of the rock sample of the load-deformation curve of the uniaxial compression of sample and the indoor rock mechanics experiment of progress
When strain curve more coincide consistent, granule rill evolution can be used as final micro-parameter at this time, according to described final microcosmic
Parameter constructs the GRANULAR FLOW MODEL FOR of rock sample, and in the present embodiment, the uniaxial compression for referring to the granule sample of coincideing is answered
Single compressing stress-strain curve difference of force-strain curve and the rock sample for carrying out indoor rock mechanics experiment is pre-
If first threshold in.
By debugging repeatedly, using granule rill evolution described in Tables 1 and 2, uniaxial fictitious compress recuperation load, mould are carried out
It is quasi- to have obtained single compressing stress-strain curve of granule sample, the curve and the rock for carrying out indoor rock mechanics experiment
Single compressing stress-strain curve of sample more coincide unanimously, it is possible thereby to which granule rill evolution determined by thinking is
Reasonably.
1 granule rill evolution of table
The displacement softening contact model rill evolution of table 2
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 of building.It is illustrated in figure 5 particle body Model failure mode,
Wherein, 4,5 and 6 be microfissure collection, and 7 and 8 pass through particle stream software (Particle after GRANULAR FLOW MODEL FOR is destroyed for circular hole
Flow Code, PFC) included FISH order, export the position coordinates (coordinate value) of each microfissure point.Pass through particle stream mould
Microfissure point coordinate value in type generates the file comprising distance between microfissure point, and this document includes each microfissure
The distance between point and other microfissure points.It can be logical from MATLAB function file is write by the way of in this step
Reading microfissure point coordinate value, then the range data obtained by multiple cycle calculations are crossed, and these range data are exported as
Apart from file.
(4) clustering is carried out by clustering algorithm, determines truncation distance and cluster centre:
As shown in fig. 6, this cluster delimit 8 microfissure groups altogether, 9-16 is the microfissure group that clustering obtains
Number, but it can also be seen that 9,10 and 11 microfissure group of serial number can be divided into a biggish microfissure group, and serial number
It can be divided into another biggish microfissure group for 12,13 and 14 microfissure groups, on the basis of merging microfissure group altogether
Three biggish 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 without principal component analysis.
(5) the three microfissure groups obtained in (4) are subjected to principal component analysis, are illustrated in figure 7 principal component analysis principle
Schematic diagram, 17 be first principal component direction, and 18 be Second principal component, direction, and 19 be data point/microfissure point, obtains three respectively
(as shown in figure 8,20 be the microfissure group of serial number 1,21 be the microfissure group of serial number 2 to microfissure group, and 22 be 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.Since first principal component direction is orthogonal with Second principal component, direction, in table 3
Just repeat no more Second principal component, slope and inclination angle.First principal component direction indicates the axially extending side of corresponding microfissure group
To, Second principal component, direction indicates the lateral magnification direction of corresponding microfissure group, particularly, Second principal component, contribution margin characterization
The degrees of expansion of corresponding microfissure group, the value is bigger, indicate that the microfissure group degrees of expansion is higher.
3 principal component analysis result of table
Embodiment two
The present invention also provides a kind of specific embodiments of the device of determining rock crackle forming development degree, since the present invention mentions
The specific embodiment of the method for the device of the determination rock crackle forming development degree of confession and aforementioned determining rock crackle forming development degree
Corresponding, the device of the determination rock crackle forming development degree can be walked by the process executed in above method specific embodiment
It is rapid to achieve the object of the present invention, therefore explaining in the method specific embodiment of above-mentioned determining rock crackle forming development degree
It is bright, it is also applied for the specific embodiment of the device of determining rock crackle forming development degree provided by the invention, 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 of determining rock crackle forming development degree, comprising:
Module 23 is constructed, for constructing the GRANULAR FLOW MODEL FOR of rock sample;
Module 24 is obtained, for the GRANULAR FLOW MODEL FOR according to building, obtains the microfissure in GRANULAR FLOW MODEL FOR destructive process
Collection and the microfissure concentrate the coordinate value of each microfissure point;
Cluster module 25, for concentrating the seat of each microfissure point according to obtained the microfissure collection and the microfissure
Scale value, the density based on microfissure point are 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 Principal Component Analysis, analyzing the axially extending direction of each microfissure group and laterally expanding
Zhang Fangxiang, and quantify development degree of each microfissure group on axially extending direction and lateral magnification direction.
The device of determination rock crackle forming development degree described in the embodiment of the present invention is obtained by the GRANULAR FLOW MODEL FOR of building
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 concentrated to the microfissure is grouped, and obtains multiple microfissures
Group;Using Principal Component Analysis, the axially extending direction and lateral magnification direction of each microfissure group are analyzed, and quantifies each microfissure
Development degree of the group on axially extending direction and lateral magnification direction.
In the specific embodiment of the device of aforementioned determining rock crackle forming development degree, further, the building mould
Block includes:
Determination unit, for using particle and cohere building 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 particle rill evolution to be arranged and coheres rill evolution;
Loading unit, for according to the particle rill evolution of setting and cohering rill evolution, to the granule sample into
Row uniaxial compression virtual load obtains single compressing stress-strain curve of the granule sample;
Comparing unit, for single compressing stress-strain curve of the obtained granule sample and progress is indoor
Single compressing stress-strain curve of the rock sample of rock mechanics experiment is compared;
Construction unit, if the difference for curve in preset first threshold, according to the particle rill evolution of setting
With cohere rill evolution, construct the GRANULAR FLOW MODEL FOR of rock sample;
Return unit, if the difference for curve not in preset first threshold, adjusts particle rill evolution and sticks
Rill evolution is tied, and returns to the loading unit and continues to execute.
In the specific embodiment of the device of aforementioned determining rock crackle forming development degree, further, the acquisition mould
Block carries out virtual uniaxial compression along long axis direction for the GRANULAR FLOW MODEL FOR according to building, obtains the GRANULAR FLOW MODEL FOR and destroys
Microfissure collection and the microfissure in the process concentrates the coordinate value of each microfissure point.
In the specific embodiment of the device of aforementioned determining rock crackle forming development degree, further, the cluster mould
Block includes:
First computing unit obtains microfissure point i for concentrating the coordinate value of each microfissure point according to the microfissure
The distance between other microfissure points j dij;
Second computing unit is surrounded by surrounding's microfissure point with smaller local density for setting cluster centre, and gathers
There are a certain distance between the biggish microfissure point in class center and other local densities, calculate the local density of microfissure point i
ρi:
ρi=∑jχ(dij-dc);
Third computing unit, for calculating microfissure point i between other microfissure point j with larger local density
Distance δi, wherein if microfissure point i has larger local density, δi=maxj(dij);Otherwise, distance δiIt is to pass through microfissure
Minimum range δ between point i and other microfissure point j with larger local densityiDetermining,
Wherein, if dij-dc< 0, then χ (dij-dc)=1, otherwise, χ (dij-dc)=0;dcIt is truncation 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 preset second threshold;When
Local density is larger local density when being not less than preset second threshold.
In the specific embodiment of the device of aforementioned determining rock crackle forming development degree, further, the quantization mould
Block, for using principal component analysis method, each microfissure group is analyzed, obtain first principal component direction, second it is main at
Divide the contribution margin in direction, the contribution margin of first principal component and Second principal component,;
Wherein, the first principal component direction indicates the axially extending direction of corresponding microfissure group, Second principal component, side
To the lateral magnification direction for indicating corresponding microfissure group, first principal component direction is orthogonal with Second principal component, direction;First is main
The contribution margin of ingredient indicates 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 a preferred embodiment of the present invention, it is noted that for those skilled in the art
For, without departing from the principles of the present invention, it can also make several improvements and retouch, these improvements and modifications
It should be regarded as protection scope of the present invention.
Claims (10)
1. a kind of method of determining rock crackle forming development degree characterized by comprising
Construct the GRANULAR FLOW MODEL FOR of rock sample;
According to the GRANULAR FLOW MODEL FOR of building, the microfissure collection and the microfissure obtained in GRANULAR FLOW MODEL FOR destructive process is concentrated often
The coordinate value of a microfissure point;
The coordinate value that each microfissure point is concentrated according to the obtained microfissure collection and the microfissure, based on microfissure point
Density is clustered, and according to the distance of each cluster centre, the microfissure point concentrated to the microfissure is grouped, obtains more
A microfissure group;
Using Principal Component Analysis, the axially extending direction and lateral magnification direction of each microfissure group are analyzed, and quantifies 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, which is characterized in that the building rock examination
The GRANULAR FLOW MODEL FOR of sample includes:
S11 using particle and coheres building granule sample, wherein the size of the granule sample and the indoor rock of progress
The size of the rock sample of mechanical test is identical;
S12 is arranged particle rill evolution and coheres rill evolution;
S13 according to the particle rill evolution of setting and coheres rill evolution, and it is virtual to carry out uniaxial compression to the granule sample
Load, obtains single compressing stress-strain curve of the granule sample;
S14 by single compressing stress-strain curve of the obtained granule sample and carries out indoor rock mechanics experiment
Single compressing stress-strain curve of rock sample is compared;
S15, if the difference of curve in preset first threshold, according to the particle rill evolution of setting and coheres thin sight ginseng
Number, constructs the GRANULAR FLOW MODEL FOR of rock sample;
Otherwise S16 then adjusts particle rill evolution and coheres rill evolution, and return to S13 and continue to execute.
3. the method for determining rock crackle forming development degree according to claim 1, which is characterized in that described according to building
GRANULAR FLOW MODEL FOR, the microfissure collection and the microfissure obtained in the GRANULAR FLOW MODEL FOR destructive process concentrate each microfissure point
Coordinate value include:
According to the GRANULAR FLOW MODEL FOR of building, virtual uniaxial compression is carried out along long axis direction, the GRANULAR FLOW MODEL FOR is obtained and destroyed
Microfissure collection and the microfissure in journey concentrate the coordinate value of each microfissure point.
4. the method for determining rock crackle forming development degree according to claim 1, which is characterized in that the basis 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 concentrated to the microfissure is grouped, and obtains multiple microfissure group packets
It includes:
The coordinate value that each microfissure point is concentrated according to the microfissure obtains between microfissure point i and other microfissure points j
Distance dij;
If cluster centre is surrounded by surrounding's microfissure point with smaller local density, cluster centre and other local densities are larger
Microfissure point between there are a certain distance, calculate the local density ρ of microfissure point ii:
ρi=∑jχ(dij-dc);
Microfissure point i is calculated from the distance between other microfissure point j with larger local density δiIf, wherein microfissure point
I has larger local density, then δi=maxj(dij);Otherwise, distance δiIt is with other by microfissure point i with larger part
Minimum range δ between the microfissure point j of densityiDetermining, δi=minJ: ρ j > ρ i(dij);
Wherein, if dij-dc< 0, then χ (dij-dc)=1, otherwise, χ (dij-dc)=0;dcIt is truncation distance, dijIt is microfissure point i
The distance between other microfissure points j;It is smaller local density when local density is less than preset second threshold;Work as part
Density is larger local density when being not less than preset second threshold.
5. the method for determining rock crackle forming development degree according to claim 1 or 4, which is characterized in that described using master
Componential analysis, analyzes the axially extending direction and lateral magnification direction of each microfissure group, and quantifies each microfissure group in axial direction
Development degree on extending direction 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 indicates the axially extending direction of corresponding microfissure group, Second principal component, direction table
Show that the lateral magnification direction of corresponding microfissure group, first principal component direction are orthogonal with Second principal component, direction;First principal component
Contribution margin indicate development degree of the corresponding microfissure group on axially extending direction, the contribution margin of Second principal component, indicates phase
Development degree of the microfissure group answered on lateral magnification direction.
6. a kind of device of determining rock crackle forming development degree characterized by comprising
Module is constructed, for constructing the GRANULAR FLOW MODEL FOR of rock sample;
Module is obtained, for the GRANULAR FLOW MODEL FOR according to building, obtains the microfissure collection in GRANULAR FLOW MODEL FOR destructive process and institute
State the coordinate value that microfissure concentrates each microfissure point;
Cluster module, for concentrating the coordinate value of each microfissure point according to obtained the microfissure collection and the microfissure,
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, obtains multiple microfissure groups;
Quantization modules, for analyzing the axially extending direction and lateral magnification direction of each microfissure group using Principal Component Analysis,
And quantify development degree of each microfissure group on axially extending direction and lateral magnification direction.
7. the device of determining rock crackle forming development degree according to claim 6, which is characterized in that the building module packet
It includes:
Determination unit, for using particle and cohering building granule sample, wherein the size and progress of the granule sample
The size of the rock sample of indoor rock mechanics experiment is identical;
Setting unit, for particle rill evolution to be arranged and coheres rill evolution;
Loading unit carries out the granule sample single for according to the particle rill evolution of setting and cohering rill evolution
Axis compresses virtual load, obtains single compressing stress-strain curve of the granule sample;
Comparing unit, for by single compressing stress-strain curve of the obtained granule sample and carrying out indoor rock
Single compressing stress-strain curve of the rock sample of mechanical test is compared;
Construction unit, if the difference for curve in preset first threshold, according to the particle rill evolution of setting and sticks
Rill evolution is tied, the GRANULAR FLOW MODEL FOR of rock sample is constructed;
Return unit, if the difference for curve not in preset first threshold, adjusts particle rill evolution and coheres thin
Parameter is seen, and returns to the loading unit and continues to execute.
8. the device of determining rock crackle forming development degree according to claim 6, which is characterized in that the acquisition module,
For the GRANULAR FLOW MODEL FOR according to building, virtual uniaxial compression is carried out along long axis direction, the GRANULAR FLOW MODEL FOR is obtained and destroyed
Microfissure collection and the microfissure in journey concentrate the coordinate value of each microfissure point.
9. the device of determining rock crackle forming development degree according to claim 6, which is characterized in that the cluster module packet
It includes:
First computing unit obtains microfissure point i and its for concentrating the coordinate value of each microfissure point according to the microfissure
The distance between his microfissure point j dij;
Second computing unit is surrounded for setting cluster centre by surrounding's microfissure point with smaller local density, and in cluster
There are a certain distance between the biggish microfissure point of the heart and other local densities, calculate the local density ρ of microfissure point ii:
ρi=∑jχ(dij-dc);
Third computing unit, for calculating microfissure point i from other the distance between microfissure point j with larger local density
δi, wherein if microfissure point i has larger local density, δi=maxj(dij);Otherwise, distance δiIt is by microfissure point i
Minimum range δ between other microfissure point j with larger local densityiDetermining, δi=minJ: ρ j > ρ i(dij);
Wherein, if dij-dc< 0, then χ (dij-dc)=1, otherwise, χ (dij-dc)=0;dcIt is truncation distance, dijIt is microfissure point i
The distance between other microfissure points j;It is smaller local density when local density is less than preset second threshold;Work as part
Density is larger local density when being not less than preset second threshold.
10. the device of determining rock crackle forming development degree according to claim 6 or 9, which is characterized in that the quantization mould
Block, for using principal component analysis method, each microfissure group is analyzed, obtain first principal component direction, second it is main at
Divide the contribution margin in direction, the contribution margin of first principal component and Second principal component,;
Wherein, the first principal component direction indicates the axially extending direction of corresponding microfissure group, Second principal component, direction table
Show that the lateral magnification direction of corresponding microfissure group, first principal component direction are orthogonal with Second principal component, direction;First principal component
Contribution margin indicate development degree of the corresponding microfissure group on axially extending direction, the contribution margin of Second principal component, indicates 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 CN106706884A (en) | 2017-05-24 |
CN106706884B true 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) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109085321B (en) * | 2018-07-25 | 2021-06-01 | 石家庄铁道大学 | Calibration method of rock mesomechanics parameters and terminal equipment |
CN113177302B (en) * | 2021-04-14 | 2023-07-14 | 北京科技大学 | Crack rock particle flow model construction and new crack occurrence analysis method |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2670181C (en) * | 2006-11-29 | 2013-03-19 | Baker Hughes Incorporated | Discrete element modeling of rock destruction under high pressure conditions |
CN101915764B (en) * | 2010-08-10 | 2012-02-01 | 武汉武大卓越科技有限责任公司 | 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 |
CN105093264B (en) * | 2014-05-12 | 2018-10-19 | 中国石油化工股份有限公司 | Rock stratum effective fracture development belt prediction technique and system |
CN105844032B (en) * | 2016-03-30 | 2020-07-24 | 李新卫 | Numerical simulation research method for irregular particles based on particle flow microwave induced damage |
-
2017
- 2017-01-11 CN CN201710017322.0A patent/CN106706884B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN106706884A (en) | 2017-05-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Tran et al. | A benchmark of batch-effect correction methods for single-cell RNA sequencing data | |
Patterson-Cross et al. | Selecting single cell clustering parameter values using subsampling-based robustness metrics | |
Ioviţă | Comparing stone tool resharpening trajectories with the aid of elliptical Fourier analysis | |
Chung et al. | Evaluation of the anisotropy of the void distribution and the stiffness of lightweight aggregates using CT imaging | |
CN106706884B (en) | A kind of method and device of determining rock crackle forming development degree | |
CN104866863B (en) | A kind of biomarker screening technique | |
CN110473597A (en) | Conglomerate mechanical property evaluation analysis method and system | |
CN103890560A (en) | Component fracture evaluation device, component fracture evaluation method, and computer program | |
CN102841985B (en) | Method for identifying key proteins based on characteristics of structural domain | |
Lehký et al. | Solving inverse structural reliability problem using artificial neural networks and small-sample simulation | |
CN105389471A (en) | Method for reducing training set of machine learning | |
Leichsenring et al. | Numerical simulation of wooden structures with polymorphic uncertainty in material properties | |
Ghosh et al. | Evaluating the morphology of erythrocyte population: An approach based on atomic force microscopy and flow cytometry | |
CN110472363B (en) | Surrounding rock deformation grade prediction method and system suitable for high-speed railway tunnel | |
Pietruszczak et al. | Description of induced anisotropy in microstructure of granular soils | |
Saputri et al. | A study of cross-national differences in Happiness factors using machine learning approach | |
Garg et al. | A new simulation approach of genetic programming in modelling of soil water retention property of unsaturated soil | |
Hafeez et al. | GPU-based real-time detection and analysis of biological targets using solid-state nanopores | |
Shafayet Zamil et al. | A multiscale FEA framework for bridging cell-wall to tissue-scale mechanical properties: the contributions of middle lamella interface and cell shape | |
EP2948885B1 (en) | Methods, systems, and computer readable media for data analysis and inference of particle diffusion in target materials and target material simulants | |
Yu et al. | A synchronized statistical characterization of size dependence and random variation of breakage strength of individual brittle particles | |
Ganusov et al. | Negative correlation between average speed and average turning angle naturally arises for sparsely sampled cell movement data | |
Jia et al. | Accounting for technical noise in single-cell RNA sequencing analysis | |
Jung et al. | Cryo-EM tomography and automatic segmentation delineate modular structures in the postsynaptic density | |
Wang et al. | Origins of brain tissue elasticity under multiple loading modes by analyzing the microstructure-based models |
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 |