CN110243666A - A kind of soil body plane of fracture based on shearing strain determines method and system - Google Patents
A kind of soil body plane of fracture based on shearing strain determines method and system Download PDFInfo
- Publication number
- CN110243666A CN110243666A CN201910598198.0A CN201910598198A CN110243666A CN 110243666 A CN110243666 A CN 110243666A CN 201910598198 A CN201910598198 A CN 201910598198A CN 110243666 A CN110243666 A CN 110243666A
- Authority
- CN
- China
- Prior art keywords
- strain
- soil body
- node
- coordinate
- displacement
- 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.)
- Pending
Links
- 239000002689 soil Substances 0.000 title claims abstract description 132
- 238000010008 shearing Methods 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000006073 displacement reaction Methods 0.000 claims abstract description 121
- 238000004364 calculation method Methods 0.000 claims abstract description 53
- 239000000700 radioactive tracer Substances 0.000 claims abstract description 34
- 238000005516 engineering process Methods 0.000 claims description 7
- 238000012163 sequencing technique Methods 0.000 claims description 6
- 238000000917 particle-image velocimetry Methods 0.000 description 8
- 230000008859 change Effects 0.000 description 5
- 239000005321 cobalt glass Substances 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000037237 body shape Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 238000007493 shaping process Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 239000005341 toughened glass Substances 0.000 description 1
Classifications
-
- 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
-
- 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/0067—Fracture or rupture
-
- 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/02—Details not specific for a particular testing method
- G01N2203/025—Geometry of the test
- G01N2203/0258—Non axial, i.e. the forces not being applied along an axis of symmetry of the specimen
-
- 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/02—Details not specific for a particular testing method
- G01N2203/06—Indicating or recording means; Sensing means
- G01N2203/0641—Indicating or recording means; Sensing means using optical, X-ray, ultraviolet, infrared or similar detectors
-
- 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/02—Details not specific for a particular testing method
- G01N2203/06—Indicating or recording means; Sensing means
- G01N2203/067—Parameter measured for estimating the property
- G01N2203/0682—Spatial dimension, e.g. length, area, angle
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The present invention relates to a kind of soil body planes of fracture based on shearing strain to determine method and system, is related to soil survey field.The following steps are included: S1: obtaining one equipped with the soil body as model casing, and the soil body is deformed, the time of the deformation of the soil body is divided at least two period, observe the deformation of the soil body and obtains the displacement field of soil body tracer point corresponding with each period;S2: a period is obtained as calculation interval, and constructs three node strain calculation units in the model casing of the calculation interval;S3: obtaining tracer point displacement field interpolation according to the displacement field of the soil body tracer point, and the modal displacement coordinate of the three nodes strain calculation unit is calculated according to the tracer point displacement field interpolation, the three nodes strain calculation unit.This programme solves the technical issues of form for how accurately determining the soil body plane of fracture, suitable for determining the form of the soil body plane of fracture.
Description
Technical field
The present invention relates to soil survey field, in particular to a kind of soil body plane of fracture based on shearing strain determines method and is
System.
Background technique
It in the design of the supporting constructions such as side slope, barricade, pattern foundation pit supporting structure and abutment, needs to carry out stability analysis, and props up
The determination of the soil body plane of fracture is the important prerequisite condition of stable calculation after gear works.The judgement side of the existing soil body plane of fracture
Method mainly have the 2 kinds: 1st kind of method be using layering smalt carry out model test, by the damage -form of the naked-eye observation soil body come
It is fitted the plane of fracture of the soil body;Since the observation of people has certain subjectivity and randomness, and plane of fracture during observation
In dynamic change, therefore, this method is difficult to effectively determine the actual form of the plane of fracture.2nd kind of method is based on particle picture
Velocity measuring technique (Particle Image Velocimetry, abbreviation PIV) obtains the displacement field in soil deformation destructive process, so
It is artificial to draw a circle to approve the soil body plane of fracture afterwards according to the soil body for generating displacement and its lower borderline region for stablizing the soil body, by human factor shadow
It rings random larger.Although this method can obtain the displacement field during damage of soil body, do not have between the plane of fracture and displacement field
Direct corresponding relationship is only capable of determining approximately plane of fracture form.
Summary of the invention
The technical problem to be solved by the present invention is to how accurately determine the form of the soil body plane of fracture.
The technical scheme to solve the above technical problems is that a kind of soil body plane of fracture determination side based on shearing strain
Method, comprising the following steps:
S1: obtaining the model casing that the soil body is housed, and deform to the soil body, by the deformation of the soil body when
Between be divided at least three period, observe the deformation of the soil body and obtain soil body tracer corresponding with each period
The displacement field of point;
S2: one period of acquisition constructs as calculation interval, and in the model casing of the calculation interval
Three node strain calculation units, the three nodes strain calculation unit include first node, second node and third node;
S3: obtaining tracer point displacement field interpolation according to the displacement field of the soil body tracer point, is displaced according to the tracer point
Field interpolation, the first node, the second node and the third node calculate the section of the three nodes strain calculation unit
Point displacement coordinate, the modal displacement coordinate includes the first displacement coordinate corresponding with the first node, with second section
The corresponding second displacement coordinate of point, third displacement coordinate corresponding with the third node;
S4: described three are calculated according to first displacement coordinate, the second displacement coordinate and the third displacement coordinate
The normal strain and shearing strain of node strain calculation unit;
S5: principal strain and maximum shear strain are obtained according to the normal strain and the shearing strain;
S6: obtaining the maximum shear strain of the calculation interval previous period as original shearing strain, according to it is described most
Big shearing strain and the original shearing strain obtain maximum shear strain increment;
S7: repeating step S2- step S6 until traversing all periods;
S8: a period is obtained as objective time interval, the model casing in the objective time interval is uniformly divided
At the vertical stripe region of at least two;
S9: obtaining the maximum value in all maximum shear strain increments in the vertical strip, and by the maximum
The position coordinates of value are recorded as maximum coordinates;
S10: curve is fitted to according to least square method and the maximum coordinates, and is recorded as the soil body in the objective time interval
The geometric shape of plastic zone;
S11: repeating step S8- step S10 until traversing all periods;
S12: according to the sequencing of the period, the geometric shape for recording first perforation model casing is
Soil body plane of fracture geometric shape.
The beneficial effects of the present invention are: the process that this programme can be destroyed according to soil deformation, the shearing strain of analytical soil sample,
To directly determine plane of fracture geometric shape;With it is existing have smalt determining method and displacement field determining method compared with, this programme is not required to
It not to be needed using there is smalt to carry out layered shaping using manual measurement yet, to reduce error, how accurately be solved
The technical issues of determining the form of the soil body plane of fracture.
Based on the above technical solution, the present invention can also be improved as follows.
Further, step S1 is deformed equipped with the soil body as model casing, and to the soil body specifically, obtaining one,
The time of the deformation of the soil body is divided at least two period, the deformation of the soil body is observed based on PIV technology and is obtained
The displacement field of soil body tracer point corresponding with each period.
Further, step S4 specifically:
S41: obtaining the coordinate of three nodes of three node strain unit, and the coordinate of three nodes is respectively the
One node coordinate (xi,yi), second node coordinate (xj,yj) and third node coordinate (xm,ym), and calculated according to following equation
First node displacement coordinate (u corresponding with the first node coordinatei,vi), the second section corresponding with the second node coordinate
Point displacement coordinate (ui,vi) and third modal displacement coordinate (u corresponding with the third node coordinatei,vi):
Wherein, the area that A is surrounded by 3 nodes:
Further, step S4 specifically: according to following equation, first displacement coordinate, the second displacement coordinate and
The third displacement coordinate calculates the normal strain ε of the three nodes strain calculation unitx,εyWith shearing strain γmax:
WhereinIndicate that u seeks partial derivative to x,Indicate that v seeks partial derivative to y,Indicate that u asks inclined to y
Derivative,Indicate that v seeks partial derivative to x.
Further, step S5 specifically: according to following equation, the normal strain εx,εyWith the shearing strain γxyIt is led
Strain stress1,ε2With maximum shear strain γmax:
Another technical solution that the present invention solves above-mentioned technical problem is as follows:
A kind of soil body plane of fracture based on shearing strain determines system, including observation module, computing module and geometric form morphotype
Block, the observation module is deformed equipped with the soil body as model casing, and to the soil body for obtaining one, by the soil body
Time of deformation be divided at least two period, observe the deformation of the soil body and obtain corresponding with each period
Soil body tracer point displacement field;The computing module is used to obtain a period as calculation interval, and described
Three node strain calculation units are constructed in the model casing of calculation interval, the three nodes strain calculation unit includes first segment
Point, second node and third node;The computing module is also used to obtain tracer point according to the displacement field of the soil body tracer point
Displacement field interpolation, according to the tracer point displacement field interpolation, the first node, the second node and the third node meter
The modal displacement coordinate of the three nodes strain calculation unit is calculated, the modal displacement coordinate includes corresponding with the first node
The first displacement coordinate, second displacement coordinate corresponding with the second node, the displacement of corresponding with third node third
Coordinate;The computing module is also used to be displaced according to first displacement coordinate, the second displacement coordinate and the third and sit
Mark calculates the normal strain and shearing strain of the three nodes strain calculation unit;The computing module is also used to according to the normal strain
Principal strain and maximum shear strain are obtained with the shearing strain;The computing module be also used to obtain the calculation interval it is previous when
Between section maximum shear strain as original shearing strain, obtain maximum according to the maximum shear strain and the original shearing strain and cut to answer
Become increment;The computing module be also used to repeat to obtain the period until traverse all periods and obtain with each
Period corresponding maximum shear strain increment;The geometric shape module is for obtaining a period as target
The model casing in the objective time interval is uniformly divided into the vertical stripe region of at least two by the period;The geometric shape module
It is also used to obtain the maximum value in all maximum shear strain increments in the vertical strip, and by the position of the maximum value
Setting coordinate record is maximum coordinates;The geometric shape module is also used to be fitted to according to least square method and the maximum coordinates
Curve, and it is recorded as the geometric shape in Plastic flow of soil area in the objective time interval;The geometric shape module is also used to repeat to obtain
The period is taken until traversing all periods and obtaining the several of Plastic flow of soil area corresponding with each period
What form;The geometric shape module is also used to the sequencing according to the period, records first perforation model
The geometric shape of case is soil body plane of fracture geometric shape.
Further, the observation module becomes equipped with the soil body as model casing, and to the soil body for obtaining one
Shape, is divided at least two period for the time of the deformation of the soil body, observes the deformation of the soil body simultaneously based on PIV technology
Obtain the displacement field of soil body tracer point corresponding with each period.
Further, the coordinate of three node of the computing module for obtaining the three nodes strain unit, described three
The coordinate of a node is respectively first node coordinate (xi,yi), second node coordinate (xj,yj) and third node coordinate (xm,ym),
And first node displacement coordinate (u corresponding with the first node coordinate is calculated according to following equationi,vi) and described second
The corresponding second node displacement coordinate (u of node coordinatei,vi) and third modal displacement corresponding with third node coordinate seat
Mark (ui,vi):
Wherein, the area that A is surrounded by 3 nodes:
Further, the computing module is also used to be sat according to following equation, first displacement coordinate, the second displacement
Mark and the third displacement coordinate calculate the normal strain ε of the three nodes strain calculation unitx,εyWith shearing strain γmax:
WhereinIndicate that u seeks partial derivative to x,Indicate that v seeks partial derivative to y,Indicate that u asks inclined to y
Derivative,Indicate that v seeks partial derivative to x.
Further, the computing module is also used to according to following equation, the normal strain εx,εyWith the shearing strain γxy
Obtain principal strain ε1,ε2With maximum shear strain γmax:
The advantages of additional aspect of the invention, will be set forth in part in the description, and will partially become from the following description
It obtains obviously, or practice is recognized through the invention.
Detailed description of the invention
Fig. 1 is the method flow diagram that the embodiment of method is determined the present invention is based on the soil body plane of fracture of shearing strain;
Fig. 2 is the system structure diagram that the embodiment of system is determined the present invention is based on the soil body plane of fracture of shearing strain.
Specific embodiment
The principle and features of the present invention will be described below with reference to the accompanying drawings, and the given examples are served only to explain the present invention, and
It is non-to be used to limit the scope of the invention.
Embodiment is substantially as shown in Fig. 1:
The soil body plane of fracture in the present embodiment based on shearing strain determines method, comprising the following steps:
S1: it obtains one and is deformed equipped with the soil body as model casing, and to the soil body, the deformation in the present embodiment can be
It is squeezed and deformed, shear-deformable or stretcher strain, the time of the deformation of the soil body is divided at least three period, observes the soil body
Deform and obtain the displacement field of soil body tracer point corresponding with each period, in the present embodiment can by the deformation of the soil body when
Between be divided into 3 periods;
S2: a period is obtained as calculation interval, and constructs three node strain gauges in the model casing of calculation interval
Unit is calculated, three node strain calculation units include first node, second node and third node;
S3: obtaining tracer point displacement field interpolation according to the displacement field of soil body tracer point, according to tracer point displacement field interpolation,
One node, second node and third node calculate the modal displacement coordinate of three node strain calculation units, modal displacement coordinate packet
The first displacement coordinate corresponding with first node is included, second displacement coordinate corresponding with second node is corresponding with third node
Third displacement coordinate;
S4: three node strain calculation units are calculated according to the first displacement coordinate, second displacement coordinate and third displacement coordinate
Normal strain and shearing strain;
S5: principal strain and maximum shear strain are obtained according to normal strain and shearing strain;
S6: the maximum shear strain of calculation interval previous period is obtained as original shearing strain, according to maximum shear strain
Maximum shear strain increment is obtained with original shearing strain;
S7: repetition step S2- step S6 up to needing to be repeated 2 times in traversal all periods, the present embodiment, ability time
Go through 3 periods;
S8: a period is obtained as objective time interval, it is vertical that the model casing in objective time interval is uniformly divided at least two
Stripe region, the vertical stripe region in the present embodiment can be 3;
S9: obtaining the maximum value in all maximum shear strain increments in vertical strip, and by the position coordinates of maximum value
It is recorded as maximum coordinates;
S10: being fitted to curve according to least square method and maximum coordinates, and is recorded as Plastic flow of soil area in objective time interval
Geometric shape;
S11: step S8- step S10 is repeated up to traversing all periods, is repeated twice in the present embodiment;
S12: according to the sequencing of period, the geometric shape for recording first failure model case is that the soil body plane of fracture is several
What form.
The beneficial effects of the present invention are: the process that this programme can be destroyed according to soil deformation, the shearing strain of analytical soil sample,
To directly determine plane of fracture geometric shape;With it is existing have smalt determining method and displacement field determining method compared with, this programme is not required to
It not to be needed using there is smalt to carry out layered shaping using manual measurement yet, to reduce error, how accurately be solved
The technical issues of determining the form of the soil body plane of fracture;
In the present embodiment, while loading beginning, tempered glass is penetrated using high definition camera, at a certain time interval
(0.5s-1.0s) photographs to record the soil body of the different moments under load action, and the photograph of the adjacent moment soil body is analyzed by PIV technology
Piece situation of change, deformation of the available test soil body particle in the period, carries out PIV points to the photo of all periods
The displacement field variation of the soil body in entire loading procedure can be obtained in analysis.
Based on the above technical solution, the present invention can also be improved as follows.
Optionally, in some other embodiments, step S1 is equipped with the soil body as model casing specifically, obtaining one, and
The soil body is deformed, the time of the deformation of the soil body is divided at least two period, the change based on the PIV technology observation soil body
Shape and the displacement field for obtaining soil body tracer point corresponding with each period.
Optionally, in some other embodiments, step S4 specifically:
S41: obtaining the coordinate of three nodes of three nodes strain unit, and the coordinate of three nodes is respectively that first node is sat
Mark (xi,yi), second node coordinate (xj,yj) and third node coordinate (xm,ym), and calculated according to following equation and first segment
The corresponding first node displacement coordinate (u of point coordinatei,vi), second node displacement coordinate (u corresponding with second node coordinatei,
vi) and third modal displacement coordinate (u corresponding with third node coordinatei,vi):
Wherein, the area that A is surrounded by 3 nodes:
Optionally, in some other embodiments, step S5 specifically: according to following equation, the first displacement coordinate, second
Displacement coordinate and third displacement coordinate calculate the normal strain ε of three node strain calculation unitsx,εyWith shearing strain γmax:
WhereinIndicate that u seeks partial derivative to x,Indicate that v seeks partial derivative to y,Indicate that u asks inclined to y
Derivative,Indicate that v seeks partial derivative to x.
Optionally, in some other embodiments, step S6 specifically: according to following equation, normal strain εx,εyIt is answered with cutting
Become γxyObtain principal strain ε1,ε2With maximum shear strain γmax:
The another embodiment of above-mentioned technical proposal are as follows:
As shown in Fig. 2, a kind of soil body plane of fracture based on shearing strain determines system, including observation module 1, computing module
2 and geometric shape module 3, observation module 1 is deformed equipped with the soil body as model casing, and to the soil body for obtaining one, by
The time of the deformation of the soil body is divided at least three period, observes the deformation of the soil body and obtains soil corresponding with each period
The displacement field of body tracer point;Computing module 2 is used to obtain a period as calculation interval, and in the model casing of calculation interval
Three node strain calculation units of interior building, three node strain calculation units include first node, second node and third node;Meter
Calculate module 2 be also used to obtain tracer point displacement field interpolation according to the displacement field of soil body tracer point, according to tracer point displacement field interpolation,
First node, second node and third node calculate the modal displacement coordinate of three node strain calculation units, modal displacement coordinate
Including the first displacement coordinate corresponding with first node, second displacement coordinate corresponding with second node is corresponding with third node
Third displacement coordinate;Computing module 2 is also used to be calculated according to the first displacement coordinate, second displacement coordinate and third displacement coordinate
The normal strain and shearing strain of three node strain calculation units;Computing module 2 is also used to obtain master according to normal strain and shearing strain and answers
Change and maximum shear strain;The maximum shear strain that computing module 2 is also used to obtain the calculation interval previous period is cut as original
Strain, obtains maximum shear strain increment according to maximum shear strain and original shearing strain;Computing module 2 is also used to repeat acquisition time
Duan Zhizhi traverses all periods and obtains maximum shear strain increment corresponding with each period;Geometric shape module 3 is used for
A period is obtained as objective time interval, the model casing in objective time interval is uniformly divided into the vertical stripe region of at least two;It is several
What configuration module 3 is also used to obtain the maximum value in all maximum shear strain increments in vertical strip, and by the position of maximum value
Setting coordinate record is maximum coordinates;Geometric shape module 3 is also used to be fitted to curve according to least square method and maximum coordinates, and
It is recorded as the geometric shape in Plastic flow of soil area in objective time interval;Geometric shape module 3 be also used to repeat acquisition time section until time
It goes through all periods and obtains the geometric shape in Plastic flow of soil corresponding with each period area;Geometric shape module 3 is also used to
According to the sequencing of period, the geometric shape for recording first failure model case is soil body plane of fracture geometric shape.
Optionally, in some other embodiments, observation module 1 is equipped with the soil body as model casing for obtaining one, and
The soil body is deformed, the time of the deformation of the soil body is divided at least two period, the change based on the PIV technology observation soil body
Shape and the displacement field for obtaining soil body tracer point corresponding with each period.
Optionally, in some other embodiments, computing module 2 is used to obtain three nodes of three nodes strain unit
Coordinate, the coordinate of three nodes are respectively first node coordinate (xi,yi), second node coordinate (xj,yj) and third node coordinate
(xm,ym), and first node displacement coordinate (u corresponding with first node coordinate is calculated according to following equationi,vi) and second
The corresponding second node displacement coordinate (u of node coordinatei,vi) and third modal displacement coordinate corresponding with third node coordinate
(ui,vi):
Wherein, the area that A is surrounded by 3 nodes:
Optionally, in some other embodiments, computing module 2 is also used to according to following equation, the first displacement coordinate,
Two displacement coordinates and third displacement coordinate calculate the normal strain ε of three node strain calculation unitsx,εyWith shearing strain γmax:
WhereinIndicate that u seeks partial derivative to x,Indicate that v seeks partial derivative to y,Indicate that u asks inclined to y
Derivative,Indicate that v seeks partial derivative to x.
Optionally, in some other embodiments, computing module 2 is also used to according to following equation, normal strain εx,εyWith cut
Strain γxyObtain principal strain ε1,ε2With maximum shear strain γmax:
It should be noted that the various embodiments described above are product embodiments corresponding with above-mentioned each method embodiment, for this
In embodiment the explanation of each constructional device and optional embodiment can with reference in above-mentioned each method embodiment pair it should be noted that
This is repeated no more.
Reader should be understood that in the description of this specification reference term " one embodiment ", " is shown " some embodiments "
The description of example ", " specific example " or " some examples " etc. mean specific features described in conjunction with this embodiment or example, structure,
Material or feature are included at least one embodiment or example of the invention.In the present specification, above-mentioned term is shown
The statement of meaning property need not be directed to identical embodiment or example.Moreover, particular features, structures, materials, or characteristics described
It may be combined in any suitable manner in any one or more of the embodiments or examples.In addition, without conflicting with each other, this
The technical staff in field can be by the spy of different embodiments or examples described in this specification and different embodiments or examples
Sign is combined.
It is apparent to those skilled in the art that for convenience of description and succinctly, the dress of foregoing description
The specific work process with unit is set, can refer to corresponding processes in the foregoing method embodiment, details are not described herein.
More than, only a specific embodiment of the invention, but scope of protection of the present invention is not limited thereto, and it is any to be familiar with
Those skilled in the art in the technical scope disclosed by the present invention, can readily occur in various equivalent modifications or substitutions,
These modifications or substitutions should be covered by the protection scope of the present invention.Therefore, protection scope of the present invention should be wanted with right
Subject to the protection scope asked.
Claims (10)
1. a kind of soil body plane of fracture based on shearing strain determines method, which comprises the following steps:
S1: it obtains one and is deformed equipped with the soil body as model casing, and to the soil body, by the time of the deformation of the soil body
It is divided at least three period, observe the deformation of the soil body and obtains soil body tracer point corresponding with each period
Displacement field;
S2: one period of acquisition constructs three sections as calculation interval, and in the model casing of the calculation interval
Point strain calculation unit, the three nodes strain calculation unit includes first node, second node and third node;
S3: obtaining tracer point displacement field interpolation according to the displacement field of the soil body tracer point, is inserted according to the tracer point displacement field
Value, the first node, the second node and the third node calculate the node position of the three nodes strain calculation unit
Coordinate is moved, the modal displacement coordinate includes the first displacement coordinate corresponding with the first node, with the second node pair
The second displacement coordinate answered, third displacement coordinate corresponding with the third node;
S4: three node is calculated according to first displacement coordinate, the second displacement coordinate and the third displacement coordinate
The normal strain and shearing strain of strain calculation unit;
S5: principal strain and maximum shear strain are obtained according to the normal strain and the shearing strain;
S6: the maximum shear strain of the calculation interval previous period is obtained as original shearing strain, is cut according to the maximum
Strain and the original shearing strain obtain maximum shear strain increment;
S7: repeating step S2- step S6 until traversing all periods;
S8: obtaining a period as objective time interval, by the model casing in the objective time interval be uniformly divided into
Few 2 vertical stripe regions;
S9: obtaining the maximum value in all maximum shear strain increments in the vertical strip, and by the maximum value
Position coordinates are recorded as maximum coordinates;
S10: curve is fitted to according to least square method and the maximum coordinates, and is recorded as Plastic flow of soil in the objective time interval
The geometric shape in area;
S11: repeating step S8- step S10 until traversing all periods;
S12: according to the sequencing of the period, the geometric shape for recording first perforation model casing is the soil body
Plane of fracture geometric shape.
2. the soil body plane of fracture according to claim 1 based on shearing strain determines method, it is characterised in that: step S1 is specific
To obtain one and being deformed equipped with the soil body as model casing, and to the soil body, the time of the deformation of the soil body is divided
For at least two period, the deformation of the soil body is observed based on PIV technology and obtains the soil body corresponding with each period
The displacement field of tracer point.
3. the soil body plane of fracture according to claim 1 based on shearing strain determines method, which is characterized in that step S4 is specific
Are as follows:
S41: the coordinate of three nodes of the three nodes strain unit is obtained, the coordinate of three nodes is respectively first segment
Point coordinate (xi,yi), second node coordinate (xj,yj) and third node coordinate (xm,ym), and calculated according to following equation and institute
State the corresponding first node displacement coordinate (u of first node coordinatei,vi), second node position corresponding with the second node coordinate
Move coordinate (ui,vi) and third modal displacement coordinate (u corresponding with the third node coordinatei,vi):
Wherein, the area that A is surrounded by 3 nodes:
4. the soil body plane of fracture according to claim 3 based on shearing strain determines method, which is characterized in that step S5 is specific
Are as follows: described three are calculated according to following equation, first displacement coordinate, the second displacement coordinate and the third displacement coordinate
The normal strain ε of node strain calculation unitx,εyWith shearing strain γmax:
WhereinIndicate that u seeks partial derivative to x,Indicate that v seeks partial derivative to y,Indicate that u seeks partial derivative to y,Indicate that v seeks partial derivative to x.
5. the soil body plane of fracture according to claim 4 based on shearing strain determines method, which is characterized in that step S6 is specific
Are as follows: according to following equation, the normal strain εx,εyWith the shearing strain γxyObtain principal strain ε1,ε2With maximum shear strain γmax:
6. a kind of soil body plane of fracture based on shearing strain determines system, it is characterised in that: including observation module, computing module and several
What configuration module, the observation module are deformed equipped with the soil body as model casing, and to the soil body for obtaining one, will
The time of the deformation of the soil body is divided at least two period, observe the deformation of the soil body and obtain with it is each described when
Between the corresponding soil body tracer point of section displacement field;The computing module for obtaining a period as calculation interval,
And three node strain calculation units are constructed in the model casing of the calculation interval, the three nodes strain calculation unit packet
Include first node, second node and third node;The computing module is also used to be obtained according to the displacement field of the soil body tracer point
To tracer point displacement field interpolation, according to the tracer point displacement field interpolation, the first node, the second node and described
Three nodes calculate the modal displacement coordinate of the three nodes strain calculation unit, and the modal displacement coordinate includes and described first
Corresponding first displacement coordinate of node, second displacement coordinate corresponding with the second node are corresponding with the third node
Third displacement coordinate;The computing module is also used to according to first displacement coordinate, the second displacement coordinate and described
Triple motion coordinate calculates the normal strain and shearing strain of the three nodes strain calculation unit;The computing module is also used to according to institute
It states normal strain and the shearing strain obtains principal strain and maximum shear strain;The computing module is also used to obtain the calculation interval
The maximum shear strain of previous period is obtained as original shearing strain according to the maximum shear strain and the original shearing strain
Maximum shear strain increment;The computing module is also used to repeat to obtain the period until traversing all periods and obtaining
To maximum shear strain increment corresponding with each period;The geometric shape module is for obtaining a period
As objective time interval, the model casing in the objective time interval is uniformly divided into the vertical stripe region of at least two;The geometry
Configuration module is also used to obtain the maximum value in all maximum shear strain increments in the vertical strip, and by described in most
The position coordinates being worth greatly are recorded as maximum coordinates;The geometric shape module is also used to according to least square method and the maximum seat
Mark is fitted to curve, and is recorded as the geometric shape in Plastic flow of soil area in the objective time interval;The geometric shape module is also used
In repeating to obtain the period until traversing all periods and obtaining the soil body corresponding with each period and mould
The geometric shape in property area;The geometric shape module is also used to the sequencing according to the period, records first perforation
The geometric shape of the model casing is soil body plane of fracture geometric shape.
7. the soil body plane of fracture according to claim 6 based on shearing strain determines system, it is characterised in that: the observation mould
Block is deformed equipped with the soil body as model casing, and to the soil body for obtaining one, by the time of the deformation of the soil body
It is divided at least two period, the deformation of the soil body is observed based on PIV technology and is obtained corresponding with each period
The displacement field of soil body tracer point.
8. the soil body plane of fracture according to claim 6 based on shearing strain determines system, it is characterised in that: the calculating mould
Block is used to obtain the coordinate of three nodes of the three nodes strain unit, and the coordinate of three nodes is respectively first node
Coordinate (xi,yi), second node coordinate (xj,yj) and third node coordinate (xm,ym), and according to following equation calculate with it is described
The corresponding first node displacement coordinate (u of first node coordinatei,vi), second node corresponding with the second node coordinate displacement
Coordinate (ui,vi) and third modal displacement coordinate (u corresponding with the third node coordinatei,vi):
Wherein, the area that A is surrounded by 3 nodes:
9. the soil body plane of fracture according to claim 8 based on shearing strain determines system, it is characterised in that: the calculating mould
Block is also used to be calculated according to following equation, first displacement coordinate, the second displacement coordinate and the third displacement coordinate
The normal strain ε of the three nodes strain calculation unitx,εyWith shearing strain γmax:
WhereinIndicate that u seeks partial derivative to x,Indicate that v seeks partial derivative to y,Indicate that u seeks local derviation to y
Number,Indicate that v seeks partial derivative to x.
10. the soil body plane of fracture according to claim 9 based on shearing strain determines system, it is characterised in that: the calculating
Module is also used to according to following equation, the normal strain εx,εyWith the shearing strain γxyObtain principal strain ε1,ε2It cuts and answers with maximum
Become γmax:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910598198.0A CN110243666A (en) | 2019-07-04 | 2019-07-04 | A kind of soil body plane of fracture based on shearing strain determines method and system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910598198.0A CN110243666A (en) | 2019-07-04 | 2019-07-04 | A kind of soil body plane of fracture based on shearing strain determines method and system |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110243666A true CN110243666A (en) | 2019-09-17 |
Family
ID=67890917
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910598198.0A Pending CN110243666A (en) | 2019-07-04 | 2019-07-04 | A kind of soil body plane of fracture based on shearing strain determines method and system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110243666A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110851983A (en) * | 2019-11-12 | 2020-02-28 | 湘潭大学 | Novel method for searching large-range unstable fracture dislocation surface of stope in complex environment |
CN110987335A (en) * | 2019-12-20 | 2020-04-10 | 防灾科技学院 | Method and device for determining position of potential fracture surface of two-stage reinforced soil structure |
CN113008675A (en) * | 2021-03-01 | 2021-06-22 | 浙江科技学院 | Detection method for bonding failure of undisturbed sand |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1139242A (en) * | 1996-02-02 | 1997-01-01 | 曾攀 | Numerical composite unit method for obtaining static and dynamic characteristics of engineering structure |
CN102081030A (en) * | 2010-04-08 | 2011-06-01 | 上海海事大学 | Geotechnical model test system based on macro-micro mechanics as well as elaboration test method |
KR101769952B1 (en) * | 2016-05-30 | 2017-08-21 | 울산과학기술원 | Triaxial stress analyzing method |
CN109030199A (en) * | 2018-05-22 | 2018-12-18 | 湖南理工学院 | Flexible retaining wall Earth Pressure Models test method |
CN109446731A (en) * | 2018-12-11 | 2019-03-08 | 哈工大机器人(合肥)国际创新研究院 | A kind of Geotechnical Engineering method based on ABAQUS |
-
2019
- 2019-07-04 CN CN201910598198.0A patent/CN110243666A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1139242A (en) * | 1996-02-02 | 1997-01-01 | 曾攀 | Numerical composite unit method for obtaining static and dynamic characteristics of engineering structure |
CN102081030A (en) * | 2010-04-08 | 2011-06-01 | 上海海事大学 | Geotechnical model test system based on macro-micro mechanics as well as elaboration test method |
KR101769952B1 (en) * | 2016-05-30 | 2017-08-21 | 울산과학기술원 | Triaxial stress analyzing method |
CN109030199A (en) * | 2018-05-22 | 2018-12-18 | 湖南理工学院 | Flexible retaining wall Earth Pressure Models test method |
CN109446731A (en) * | 2018-12-11 | 2019-03-08 | 哈工大机器人(合肥)国际创新研究院 | A kind of Geotechnical Engineering method based on ABAQUS |
Non-Patent Citations (5)
Title |
---|
刘金操: "非均质土体中的地基承载力试验研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
戴宏亮等: "《工程有限元及数值计算》", 30 June 2019 * |
曹耿等: "均布条形荷载下地基塑性区分析", 《低温建筑技术》 * |
李勇军: "一种确定均布条形荷载下塑性区的方法", 《辽宁工学院学报》 * |
王西彬等: "《精密制造工学基础》", 31 January 2018 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110851983A (en) * | 2019-11-12 | 2020-02-28 | 湘潭大学 | Novel method for searching large-range unstable fracture dislocation surface of stope in complex environment |
CN110987335A (en) * | 2019-12-20 | 2020-04-10 | 防灾科技学院 | Method and device for determining position of potential fracture surface of two-stage reinforced soil structure |
CN113008675A (en) * | 2021-03-01 | 2021-06-22 | 浙江科技学院 | Detection method for bonding failure of undisturbed sand |
CN113008675B (en) * | 2021-03-01 | 2024-04-26 | 浙江科技学院 | Method for detecting failure of original sand cementation |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110243666A (en) | A kind of soil body plane of fracture based on shearing strain determines method and system | |
Korumaz et al. | An integrated Terrestrial Laser Scanner (TLS), Deviation Analysis (DA) and Finite Element (FE) approach for health assessment of historical structures. A minaret case study | |
CN108824816B (en) | High-altitude long-span net frame sliding, positioning, installing and monitoring method | |
CN107423462A (en) | Workpiece considers the Prediction method for fatigue life and system of three-dimensional rough surface morphology | |
DE112020000229T5 (en) | Method and system for detecting a flatness of a road surface | |
CN105954796B (en) | A kind of method and apparatus for the hypocentral location for determining microseism | |
CN106294969A (en) | A kind of method of steel structure net rack hydraulic pressure lift precision controlling based on BIM | |
JP6201059B2 (en) | Ground shape estimation program, ground shape estimation device, and ground shape estimation method | |
CN111581836A (en) | Method for calculating stability of power transmission line landslide mass | |
Costa-Jover et al. | Using the terrestrial laser scanner and simple methodologies for geometrically assessing complex masonry vaults | |
CN110390800B (en) | Network type disaster monitoring and early warning system | |
CN111141473A (en) | Real-time display method and system for product vibration test data | |
JP2011040055A (en) | Simulation device and program | |
CN111830561B (en) | Method for constructing fault three-dimensional structure based on seismic distribution characteristics | |
DE102010063818A1 (en) | Method for automatic acquisition of mobility of patient during predetermined mobility test in clinic, involves evaluating sensor data by evaluating device to obtain mobility parameter related to mobility of person | |
CN111507015A (en) | In-service slope technical condition evaluation method | |
CN115115813B (en) | Intelligent construction method for human skeleton standard body position | |
CN106525192A (en) | Real-time river section measurement device and method | |
RU2374697C2 (en) | Experimental method for parametrisation of minimum surfaces with complex contour | |
KR101255022B1 (en) | Detecting method of building crack using point group | |
CN109308373A (en) | Consider that seismic response analysis method and device is built in place-urban effect region | |
CN106040750B (en) | Position compensation device and its position compensation method in the processing of looping mill rolling part | |
JP3642481B2 (en) | Surface stability evaluation data creation system | |
Alfaro et al. | Bioimplantable bone stress sensor | |
CN113216922A (en) | Platform and method for researching expansion rule of hydraulic fracture in layered medium |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190917 |