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 PDF

Info

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
Application number
CN201710017322.0A
Other languages
Chinese (zh)
Other versions
CN106706884A (en
Inventor
周喻
刘冰
王雪
吴顺川
周建新
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Science and Technology Beijing USTB
Original Assignee
University of Science and Technology Beijing USTB
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Science and Technology Beijing USTB filed Critical University of Science and Technology Beijing USTB
Priority to CN201710017322.0A priority Critical patent/CN106706884B/en
Publication of CN106706884A publication Critical patent/CN106706884A/en
Application granted granted Critical
Publication of CN106706884B publication Critical patent/CN106706884B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/24Earth materials
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/08Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0001Type of application of the stress
    • G01N2203/0003Steady
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0014Type of force applied
    • G01N2203/0016Tensile or compressive
    • G01N2203/0019Compressive
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0058Kind of property studied
    • G01N2203/006Crack, flaws, fracture or rupture
    • G01N2203/0062Crack or flaws
    • G01N2203/0066Propagation 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

A kind of method and device of determining rock crackle forming development degree
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 λ12,…,λ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.
CN201710017322.0A 2017-01-11 2017-01-11 A kind of method and device of determining rock crackle forming development degree Active CN106706884B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

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