CN115879325B - Sliding surface shear strength parameter inversion method, device, equipment and readable storage medium - Google Patents
Sliding surface shear strength parameter inversion method, device, equipment and readable storage medium Download PDFInfo
- Publication number
- CN115879325B CN115879325B CN202310150425.XA CN202310150425A CN115879325B CN 115879325 B CN115879325 B CN 115879325B CN 202310150425 A CN202310150425 A CN 202310150425A CN 115879325 B CN115879325 B CN 115879325B
- Authority
- CN
- China
- Prior art keywords
- shear strength
- bar
- strength parameter
- area
- point safety
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/23—Dune restoration or creation; Cliff stabilisation
Landscapes
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
Abstract
The invention provides a sliding surface shear strength parameter inversion method, a device, equipment and a readable storage medium, which relate to the technical field of parameter inversion and comprise the steps of vertically dividing a sliding body in a sliding slope to obtain a plurality of strips; dividing a first area, a second area and a third area from the plurality of strips; obtaining vertex coordinates of each bar in the first area, the second area and the third area, and calculating geometric parameters of each bar; determining a first shear strength parameter, a second shear strength parameter, a first threshold value of a first regional point safety coefficient and a second threshold value of a second regional point safety coefficient of the landslide; and calculating a first point safety coefficient of the third area, and obtaining the sliding surface shear strength parameter by changing the second shear strength parameter to enable the first point safety coefficient to be equal to the second threshold value. The method is used for solving the technical problems that in the prior art, the calculation accuracy is far less than the engineering requirement, and the calculated landslide thrust is inaccurate, so that engineering control measures are affected.
Description
Technical Field
The invention relates to the technical field of parameter inversion, in particular to a sliding surface shear strength parameter inversion method, a device, equipment and a readable storage medium.
Background
Regarding the landslide parameter inversion study, landslide can be divided into two types, one is that the main section of the landslide before sliding can be recovered, and the other is that the main section of the landslide before sliding cannot be recovered. And for the main section of the landslide before sliding can be recovered, a limit balance equation is listed through the integral stability coefficient to solve the cohesive force or the internal friction angle. And for the landslide of the main section of the landslide before unrecoverable sliding, selecting a stability coefficient according to the state of the landslide, and substituting the stability coefficient into a formula to perform back calculation. The two calculation modes are only considered by a strong deformation area or the whole body alone, and for complex landslide, the calculation precision of the two methods can not reach the engineering requirement, and can cause a large gap between cohesive force or an internal friction angle and reality, so that the calculated landslide thrust is inaccurate, thereby influencing engineering prevention measures.
Disclosure of Invention
The invention aims to provide a sliding surface shear strength parameter inversion method, a device, equipment and a readable storage medium, so as to solve the problems. In order to achieve the above purpose, the technical scheme adopted by the invention is as follows:
in a first aspect, the present application provides a slip plane shear strength parameter inversion method comprising:
vertically dividing a sliding body in a landslide to obtain a plurality of strips;
dividing a first area, a second area and a third area from the plurality of strips, wherein the second area is a strong deformation area;
obtaining the vertex coordinates of each bar in the first area, the second area and the third area, calculating the geometric parameters of each bar according to the vertex coordinates, and forming a first matrix by the geometric parameters of each bar;
determining a first shear strength parameter, a second shear strength parameter, a first threshold value of a first regional point safety coefficient and a second threshold value of a second regional point safety coefficient of the landslide;
and calculating a first point safety coefficient of the third region according to the first shear strength parameter, the second shear strength parameter, the first matrix and the first threshold value, enabling the first point safety coefficient to be equal to the second threshold value by changing the second shear strength parameter, and taking the second shear strength parameter when the first point safety coefficient is equal to the second threshold value as the sliding surface shear strength parameter.
In a second aspect, the present application further provides a slip plane shear strength parameter inversion apparatus, comprising:
a first dividing module: vertically dividing a sliding body in a landslide to obtain a plurality of strips;
a second dividing module: dividing a first area, a second area and a third area from the plurality of strips, wherein the second area is a strong deformation area;
a first acquisition module: obtaining the vertex coordinates of each bar in the first area, the second area and the third area, calculating the geometric parameters of each bar according to the vertex coordinates, and forming a first matrix by the geometric parameters of each bar;
and a second acquisition module: determining a first shear strength parameter, a second shear strength parameter, a first threshold value of a first regional point safety coefficient and a second threshold value of a second regional point safety coefficient of the landslide;
and an inversion module: and calculating a first point safety coefficient of the third region according to the first shear strength parameter, the second shear strength parameter, the first matrix and the first threshold value, enabling the first point safety coefficient to be equal to the second threshold value by changing the second shear strength parameter, and taking the second shear strength parameter when the first point safety coefficient is equal to the second threshold value as the sliding surface shear strength parameter.
In a third aspect, the present application further provides a sliding surface shear strength parameter inversion apparatus, including:
a memory for storing a computer program;
and the processor is used for realizing the step of the sliding surface shear strength parameter inversion method when executing the computer program.
In a fourth aspect, the present application also provides a readable storage medium having stored thereon a computer program which, when executed by a processor, implements the steps of the sliding surface shear strength parameter inversion method described above.
The beneficial effects of the invention are as follows:
according to the invention, the landslide is identified and partitioned by combining image analysis, and the landslide is recovered and judged to judge the parameter inversion region of the landslide, so that compared with the prior art, the overall stability and the local stability of the landslide are combined and analyzed, and the obtained result is more reasonable and accurate. Meanwhile, parameters to be inverted are determined based on the influence degree of cohesive force and an internal friction angle on landslide stability, the parameters obtained through inversion are more accurate, and engineering calculation errors are reduced.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be apparent from the description, or may be learned by practice of the embodiments of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims thereof as well as the appended drawings.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings that are needed in the embodiments will be briefly described below, it being understood that the following drawings only illustrate some embodiments of the present invention and therefore should not be considered as limiting the scope, and other related drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a schematic flow chart of a method for inverting the shear strength parameters of the sliding surface according to embodiment 1 of the present invention;
fig. 2 is a schematic diagram of the structure of a landslide in embodiment 2 of the invention;
FIG. 3 is a schematic diagram of the division of bars of a landslide in embodiment 2 of the invention;
FIG. 4 is a graph showing the safety factor of all bars in example 2 of the present invention;
fig. 5 is a schematic diagram of the structure of a landslide in embodiment 3 of the invention;
FIG. 6 is a schematic diagram showing the division of bars of a landslide in embodiment 3 of the invention;
FIG. 7 is a graph showing the safety factor of all bars in example 3 of the present invention;
FIG. 8 is a schematic diagram of a device for inverting the shear strength parameters of the sliding surface according to embodiment 4 of the present invention;
FIG. 9 is a schematic diagram of a sliding surface shear strength parameter inversion apparatus according to example 5 of the present invention.
The marks in the figure:
01. a first dividing module; 02. a second dividing module; 021. a first sorting unit; 022. a first collection unit; 023. a second collection unit; 024. a third collection unit; 03. a first acquisition module; 031. a first acquisition unit; 032. a second sorting unit; 033. a first calculation unit; 034. a second matrix constructing unit; 035. a first matrix constructing unit; 04. a second acquisition module; 041. a second calculation unit; 042. a first judgment unit; 043. a second acquisition unit; 05. an inversion module; 051. a third calculation unit; 0511. a third acquisition unit; 0512. a fourth calculation unit; 0513. a fifth calculation unit; 05131. a third judgment unit; 05132. a fourth judgment unit; 0514. a sixth calculation unit; 0515. a seventh calculation unit; 0516. an eighth calculation unit; 052. a second judgment unit;
800. a sliding surface shear strength parameter inversion device; 801. a processor; 802. a memory; 803. a multimedia component; 804. an I/O interface; 805. a communication component.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments of the present invention. The components of the embodiments of the present invention generally described and illustrated in the figures herein may be arranged and designed in a wide variety of different configurations. Thus, the following detailed description of the embodiments of the invention, as presented in the figures, is not intended to limit the scope of the invention, as claimed, but is merely representative of selected embodiments of the invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
It should be noted that: like reference numerals and letters denote like items in the following figures, and thus once an item is defined in one figure, no further definition or explanation thereof is necessary in the following figures. Meanwhile, in the description of the present invention, the terms "first", "second", and the like are used only to distinguish the description, and are not to be construed as indicating or implying relative importance.
Example 1:
the embodiment provides a sliding surface shear strength parameter inversion method.
Referring to fig. 1, the method is shown to include:
s1, vertically dividing a sliding body in a landslide to obtain a plurality of strips;
specifically, a history image before sliding of a landslide, an image after sliding, and survey data are acquired. Judging the landslide, outputting the type of the landslide, wherein the type of the landslide comprises a push type landslide, a traction type landslide and an integral type landslide, judging a strong deformation area of the landslide according to the type of the landslide, recovering the landslide, acquiring a typical sectional view before the landslide slides, and vertically dividing the typical sectional view. S2, dividing a first area, a second area and a third area from the plurality of strips, wherein the second area is a strong deformation area;
specifically, the step S2 includes:
s21, sorting a plurality of strips: sequencing all the bars from the top of the slope to the bottom of the slope, wherein the bar at the top of the slope is the first bar, and the bar at the bottom of the slope is the last bar;
s22, classifying the strip blocks in the strong deformation area into a second area B:
wherein b i Representing the ith bar, b in the second region n Representing the nth bar in the second region.
S23, classifying the bar blocks before the second region is ordered into a first region A:
wherein a is j Represents the j-th bar in the first region, a m Representing the mth bar in the first region.
S24, classifying the first bar block after the second area is ordered into a third area C, wherein the third area only comprises one bar block C, and the bar block C is not necessarily the last bar block.
S3, obtaining the vertex coordinates of each bar in the first area, the second area and the third area, calculating the geometric parameters of each bar according to the vertex coordinates, and forming a first matrix by the geometric parameters of each bar;
specifically, the step S3 includes:
s31, acquiring coordinates of four vertexes of any block in the first area, the second area and the third area, and sequentially storing the coordinates of the four vertexes in a matrix according to the sequence of an upper right corner, an upper left corner, a lower right corner and a lower left corner;
wherein the set of vertex coordinates of the bar of the first region is represented as:
the set of vertex coordinates of the bar of the second region is expressed as:
the vertex coordinate set of the bar of the third region is expressed as:
s32, sorting the vertex coordinates according to the magnitude relation of the ordinate of the vertex coordinates;
specifically, the ordinate is obtained according to the order from big to small:
wherein a is jy A set of vertex ordinates representing a bar of the first region, b iy A set of vertex ordinates representing a bar of the second region, c y And a set of vertex ordinate representing the bar of the third region.
S33, calculating the area, the bottom surface length and the bottom surface dip angle of the bar blocks based on the ordering of the vertex coordinates;
wherein S is k Represents the area of the kth bar, x 4 Representing the left lower corner abscissa of the bar, x 3 Representing the lower right-hand abscissa of the bar, y max Representing the maximum ordinate, y, of the bar 1m Representing the ordinate, y, of the ordered second of the bar 2m Representing the third ordinate, y, of the ordering of the bars min Representing the smallest ordinate of the bar.
Wherein, I k Represents the bottom surface length of the kth bar, y 4 Representing the lower left-hand ordinate, y, of the bar 3 Indicating the lower right vertical axis of the bar.
Wherein a is k Indicating the inclination of the bottom surface of the kth bar.
S34, forming a second matrix by the area of one bar, the length of the bottom surface and the inclination angle of the bottom surface, wherein the second matrix can be expressed as;
s35, forming a first matrix by all the second matrices.
S4, determining a first shear strength parameter, a second shear strength parameter, a first threshold value of a first regional point safety coefficient and a second threshold value of a second regional point safety coefficient of the landslide;
specifically, the step S4 includes:
s41, calculating the influence of cohesive force on the landslide and the influence of an internal friction angle on the landslide according to an engineering analogy method, wherein the engineering analogy method is a prior art known to a person skilled in the art, and is not repeated here;
s42, if the influence of the cohesive force on the landslide is larger than that of the internal friction angle on the landslide, taking the internal friction angle as a first shear strength parameter and taking the cohesive force as a second shear strength parameter;
otherwise, the cohesive force is taken as a first shear strength parameter, and the internal friction angle is taken as a second shear strength parameter.
In this embodiment, taking the case that the influence of the internal friction angle is larger as an example, the cohesive force is taken as a first shear strength parameter, the internal friction angle is taken as a second shear strength parameter, and the internal friction angle is inverted by using the cohesive force.
S43, determining the value of the first shear strength parameter according to the survey report.
S5, calculating a first point safety coefficient of the third region according to the first shear strength parameter, the second shear strength parameter, the first matrix and the first threshold value, enabling the first point safety coefficient to be equal to the second threshold value through changing the second shear strength parameter, and taking the second shear strength parameter when the first point safety coefficient is equal to the second threshold value as the sliding surface shear strength parameter.
Specifically, the step S5 includes:
s51, calculating a first point safety coefficient of the third area according to the second shear strength parameter, the first shear strength parameter and the first matrix;
specifically, the step S51 includes:
s511, acquiring geometric parameters of a current bar in the first matrix, wherein the current bar is any one of a first area and a second area;
s512, calculating a first point safety coefficient of the current bar according to the geometric parameter of the current bar, the current value of the second shear strength parameter and the first shear strength parameter;
wherein F is k1 A first point safety factor representing the kth bar, R k Representing the anti-slip force, T, of the kth bar k Indicating the sliding down force of the kth bar.
Specifically, the anti-slip force R k The calculation formula of (2) is as follows:
wherein c represents cohesion, V k-1 Representing the supporting force provided by the kth bar to the kth-1 bar, alpha k-1 Represents the inclination angle of the bottom surface of the kth-1 bar,
represents the internal friction angle, W k Represents the dead weight stress of the kth bar, said W k According to S k And calculating the density of the bar blocks.
Specifically, the sliding force T k The calculation formula of (2) is as follows:
s513, calculating a second point safety coefficient of the current bar and a supporting force required by the current bar based on the first point safety coefficient;
specifically, the step S513 includes:
s5131, judging the area to which the current bar belongs:
if the current bar belongs to the first region, the preset threshold is set to be a first threshold, in this embodiment, the first threshold
；/>
If the current bar belongs to the second area, the preset threshold is set to be a second threshold, the value of the second threshold is 1.05 under the natural working condition, and the value of the second threshold is 1.0 under the heavy rain working condition, in this embodiment, taking the natural working condition as an example, the second threshold is set to be
；
S5132, judging whether the first point safety coefficient is greater than or equal to a preset threshold value:
if yes, the second point safety coefficient is equal to the first point safety coefficient, and the supporting force required by the current bar is 0;
if it is
Let->
，/>
，/>
Representing a second point safety factor, said V k Representing the holding force required for the kth bar, i.e. the holding force provided by the kth+1th bar for the kth bar, fs represents a preset threshold.
If not, the second point safety coefficient is enabled to be equal to a preset threshold value, and the supporting force required by the current bar block is calculated under the condition of meeting the static balance.
If it is
Let->
The supporting force V required by the current bar k The calculation formula of (2) is as follows:
s514, calculating a second point safety coefficient and required supporting force of the next bar of the current bar under the action of the required supporting force of the current bar;
s515, sequentially calculating a second point safety coefficient and required supporting force of the rest bar under the action of the required supporting force of the last bar until the required supporting force of the last bar of the second area is calculated;
in this embodiment, the steps S511-S513 are repeated to calculate the required supporting force of the 1 st to the (m+n) th bars;
s516, calculating a first point safety coefficient of the third area based on the supporting force required by the last bar block of the second area.
Based on the supporting force required by the m+n-th bar
Calculating the first point safety factor of the c-th bar, wherein the first point safety factor of the bar is the third area because the third area has only one barFirst point security factor of domain.
Wherein F is c A first point safety factor representing a third region or c-th bar, R c Representing the anti-slip force, T, of the c-th bar c Indicating the sliding down force of the c-th bar.
S52, judging whether the first point safety coefficient of the third area is equal to the second threshold value:
specifically, judge
And->
Whether or not to equal:
if the values are equal, the current value of the second shear strength parameter is taken as the sliding surface shear strength parameter, namely the current internal friction angle phi 0 The internal friction angle of the slip surface as a reversal of the method;
if not, changing the value of the second shear strength parameter, recalculating the first point safety factor of the third region, specifically changing the value of the internal friction angle, and repeating step S51.
Example 2:
as shown in FIGS. 2-4, the landslide has a weight of 20KN/m 3 Dividing the landslide into 21 pieces, wherein 1-20 pieces are the second region, 21 pieces are the third region, and the first threshold of the second region
According to engineering analogy, the cohesion c is taken as the first shear strength parameter, the internal friction angle +.>
As a second shear strength parameter c=15 kPa was determined, let initial value +.>
=7.93°, noBreak change input +.>
Value of>
At 8.94 DEG, the first point safety factor of the third region approaches 1.05, whereby the sliding surface shear strength parameter is obtained>
= 7.942 ° as shown in table 1.
TABLE 1
In Table 1, the point safety coefficients of the 1 st to 20 th bars are the second point safety coefficient, the point safety coefficient of the 21 st bar is the first point safety coefficient, and it can be seen from Table 1 that when
At= 7.942 °, the first point safety factor of the third region is 1.052, very close to 1.05.
Example 3
As shown in FIGS. 5-7, the landslide has a weight of 20KN/m 3 Dividing the landslide into 21 strips, wherein the 1 st strip is a first area, the 2 nd-12 th strips are a second area, the 13 th strip is a third area, and the first threshold value of the second area is set
According to engineering analogy, using the cohesion c as the first shear strength parameter and the internal friction angle phi as the second shear strength parameter, determining c=13.70 kPa, let the initial value +.>
=14°, constantly changing the input +.>
Value of>
At 13.91 DEG the first point safety factor of the third zone approaches 1.05, whereby the slip shear strength parameter is obtained>
=13.9° as shown in table 2.
TABLE 2
Table 2, point safety coefficients of the 1 st to 12 th blocks are second point safety coefficients, point safety coefficient of the 13 th block is first point safety coefficient, as can be seen from table 2, when
At=13.9°, the first point safety factor of the third region is 1.052, very close to 1.05.
Example 4
As shown in fig. 8, the embodiment provides a sliding surface shear strength parameter inversion device, which includes:
first division module 01: vertically dividing a sliding body in a landslide to obtain a plurality of strips;
the second division module 02: dividing a first area, a second area and a third area from the plurality of strips, wherein the second area is a strong deformation area;
the first acquisition module 03: obtaining the vertex coordinates of each bar in the first area, the second area and the third area, calculating the geometric parameters of each bar according to the vertex coordinates, and forming a first matrix by the geometric parameters of each bar;
the second acquisition module 04: determining a first shear strength parameter, a second shear strength parameter, a first threshold value of a first regional point safety coefficient and a second threshold value of a second regional point safety coefficient of the landslide;
inversion module 05: and calculating a first point safety coefficient of the third region according to the first shear strength parameter, the second shear strength parameter, the first matrix and the first threshold value, enabling the first point safety coefficient to be equal to the second threshold value by changing the second shear strength parameter, and taking the second shear strength parameter when the first point safety coefficient is equal to the second threshold value as the sliding surface shear strength parameter.
Based on the above embodiments, the second dividing module 02 specifically includes:
the first sorting unit 021: ordering the plurality of bars: the bar at the top of the landslide is made to be the first bar, and the bar at the bottom of the landslide is made to be the last bar;
the first collection unit 022: grouping the strips of the strong deformation area into a second area;
a second collection unit 023: grouping the bar blocks before the second region is ordered into a first region;
third aggregation unit 024: the first bar after the second region is ordered is grouped into a third region.
Based on the above embodiments, the first obtaining module 03 specifically includes:
first acquisition unit 031: obtaining vertex coordinates of any one block in the first area, the second area and the third area;
a second sorting unit 032: sorting the vertex coordinates according to the magnitude relation of the ordinate of the vertex coordinates;
first calculation unit 033: calculating the area, the bottom surface length and the bottom surface dip angle of the bar blocks based on the ordering of the vertex coordinates;
second matrix constructing unit 034: forming a second matrix by the area of one bar, the length of the bottom surface and the inclination angle of the bottom surface;
first matrix constructing unit 035: the first matrix is formed by all the second matrices.
Based on the above embodiments, the second obtaining module 04 specifically includes:
the second calculation unit 041: calculating the influence of the cohesive force on the landslide and the influence of the internal friction angle on the landslide according to an engineering analogy method;
first judgment unit 042: if the influence of the cohesive force on the landslide is larger than the influence of the internal friction angle on the landslide, taking the internal friction angle as a first shear strength parameter and taking the cohesive force as a second shear strength parameter;
otherwise, the cohesive force is taken as a first shear strength parameter, and the internal friction angle is taken as a second shear strength parameter.
Second acquisition unit 043: and determining the value of the first shear strength parameter according to the survey report.
Based on the above embodiments, the inversion module 05 specifically includes:
third computing unit 051: calculating a first point safety coefficient of the third area according to the second shear strength parameter, the first shear strength parameter and a first matrix;
the second judgment unit 052: judging whether the first point safety coefficient of the third area is equal to the second threshold value or not:
if the two parameters are equal, taking the current value of the second shear strength parameter as the sliding surface shear strength parameter;
if the first and second points are not equal, changing the value of the second shear strength parameter, and recalculating the first point safety coefficient of the third area.
Based on the above embodiments, the third computing unit specifically includes:
third acquisition unit 0511: obtaining geometric parameters of a current bar in a first matrix, wherein the current bar is any one of a first area and a second area;
fourth calculation unit 0512: calculating a first point safety coefficient of the current bar according to the geometric parameter of the current bar, the current value of the second shear strength parameter and the first shear strength parameter;
fifth calculation unit 0513: calculating a second point safety coefficient of the current bar and a supporting force required by the current bar based on the first point safety coefficient;
sixth calculation unit 0514: calculating a second point safety coefficient and a required supporting force of a next bar of the current bar under the action of the required supporting force of the current bar;
seventh calculation unit 0515: sequentially calculating a second point safety coefficient and a required supporting force of the rest bar blocks under the action of the required supporting force of the last bar block until the required supporting force of the last bar block of the second area is calculated;
eighth calculation unit 0516: the first point security factor of the third region is calculated based on the holding force required for the last bar of the second region.
Based on the above embodiments, the fifth computing unit 0513 specifically includes:
third judgment unit 05131: judging the area to which the current bar belongs:
if the current bar block belongs to the first area, enabling a preset threshold value to be a first threshold value; if the current bar block belongs to the second area, enabling the preset threshold value to be a second threshold value;
fourth judgment unit 05132: judging whether the first point safety coefficient is larger than or equal to a preset threshold value or not:
if yes, the second point safety coefficient is equal to the first point safety coefficient, and the supporting force required by the current bar is 0;
if not, the second point safety coefficient is enabled to be equal to a preset threshold value, and the supporting force required by the current bar block is calculated under the condition of meeting the static balance.
It should be noted that, regarding the apparatus in the above embodiments, the specific manner in which the respective modules perform the operations has been described in detail in the embodiments regarding the method, and will not be described in detail herein.
Example 5:
corresponding to the above method embodiment, a sliding surface shear strength parameter inversion apparatus is further provided in this embodiment, and a sliding surface shear strength parameter inversion apparatus described below and a sliding surface shear strength parameter inversion method described above may be referred to correspondingly.
FIG. 9 is a block diagram illustrating a slip plane shear strength parameter inversion apparatus 800, according to an example embodiment. As shown in fig. 9, the slip plane shear strength parameter inversion apparatus 800 may include: a processor 801, a memory 802. The slide shear strength parameter inversion apparatus 800 may also include one or more of a multimedia component 803, an I/O interface 804, and a communication component 805.
The processor 801 is configured to control the overall operation of the sliding surface shear strength parameter inversion apparatus 800 to perform all or part of the steps of the sliding surface shear strength parameter inversion method described above. The memory 802 is used to store various types of data to support the operation of the slide shear strength parameter inversion device 800, which may include, for example, instructions for any application or method operating on the slide shear strength parameter inversion device 800, as well as application-related data, such as contact data, messaging, pictures, audio, video, and the like. The Memory 802 may be implemented by any type or combination of volatile or non-volatile Memory devices, such as static random access Memory (Static Random Access Memory, SRAM for short), electrically erasable programmable Read-Only Memory (ElectricallyErasable Programmable Read-Only Memory, EEPROM for short), erasable programmable Read-Only Memory (Erasable Programmable Read-Only Memory, EPROM for short), programmable Read-Only Memory (Programmable Read-Only Memory, PROM for short), read-Only Memory (ROM for short), magnetic Memory, flash Memory, magnetic disk, or optical disk. The multimedia component 803 may include a screen and an audio component. Wherein the screen may be, for example, a touch screen, the audio component being for outputting and/or inputting audio signals. For example, the audio component may include a microphone for receiving external audio signals. The received audio signals may be further stored in the memory 802 or transmitted through the communication component 805. The audio assembly further comprises at least one speaker for outputting audio signals. The I/O interface 804 provides an interface between the processor 801 and other interface modules, which may be a keyboard, mouse, buttons, etc. These buttons may be virtual buttons or physical buttons. The communication component 805 is configured to provide wired or wireless communication between the slide shear strength parameter inversion apparatus 800 and other apparatus. Wireless communication, such as Wi-Fi, bluetooth, near Field Communication (NFC), 2G, 3G, or 4G, or a combination of one or more thereof, the corresponding communication component 805 may include: wi-Fi module, bluetooth module, NFC module.
In an exemplary embodiment, the slip plane shear strength parameter inversion apparatus 800 may be implemented by one or more application specific integrated circuits (Application Specific Integrated Circuit, ASIC), digital signal processors (DigitalSignal Processor, DSP), digital signal processing apparatus (Digital SignalProcessing Device, DSPD), programmable logic devices (Programmable Logic Device, PLD), field programmable gate arrays (Field Programmable Gate Array, FPGA), controllers, microcontrollers, microprocessors, or other electronic components for performing the slip plane shear strength parameter inversion methods described above.
In another exemplary embodiment, a computer readable storage medium is also provided, comprising program instructions which, when executed by a processor, implement the steps of the above-described slide shear strength parameter inversion method. For example, the computer readable storage medium may be the memory 802 described above that includes program instructions executable by the processor 801 of the slide shear strength parameter inversion apparatus 800 to perform the slide shear strength parameter inversion method described above.
Example 6:
corresponding to the above method embodiment, a readable storage medium is also provided in this embodiment, and a readable storage medium described below and a sliding surface shear strength parameter inversion method described above may be referred to correspondingly.
A readable storage medium having stored thereon a computer program which when executed by a processor performs the steps of the sliding surface shear strength parameter inversion method of the above method embodiments.
The readable storage medium may be a usb disk, a removable hard disk, a Read-Only Memory (ROM), a random access Memory (RandomAccess Memory, RAM), a magnetic disk, or an optical disk, and the like.
The above description is only of the preferred embodiments of the present invention and is not intended to limit the present invention, but various modifications and variations can be made to the present invention by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
The foregoing is merely illustrative of the present invention, and the present invention is not limited thereto, and any person skilled in the art will readily recognize that variations or substitutions are within the scope of the present invention. Therefore, the protection scope of the invention is subject to the protection scope of the claims.
Claims (8)
1. The method for inverting the shear strength parameter of the sliding surface is characterized by comprising the following steps of:
vertically dividing a sliding body in a landslide to obtain a plurality of strips;
dividing a first area, a second area and a third area from the plurality of strips, wherein the second area is a strong deformation area;
obtaining the vertex coordinates of each bar in the first area, the second area and the third area, calculating the geometric parameters of each bar according to the vertex coordinates, and forming a first matrix by the geometric parameters of each bar;
determining a first shear strength parameter, a second shear strength parameter, a first threshold value of a first regional point safety coefficient and a second threshold value of a second regional point safety coefficient of the landslide;
calculating a first point safety factor for the third region based on the first shear strength parameter, the second shear strength parameter, the first matrix, and the first threshold value, comprising:
obtaining geometric parameters of a current bar in a first matrix, wherein the current bar is any one of a first area and a second area;
calculating a first point safety coefficient of the current bar according to the geometric parameter of the current bar, the current value of the second shear strength parameter and the first shear strength parameter:
wherein F is k1 A first point safety factor representing the kth bar, R k Representing the anti-slip force, T, of the kth bar k Representing the sliding down force of the kth bar;
calculating a second point safety coefficient of the current bar and a supporting force required by the current bar based on the first point safety coefficient;
calculating a second point safety coefficient and a required supporting force of a next bar of the current bar under the action of the required supporting force of the current bar;
sequentially calculating a second point safety coefficient and a required supporting force of the rest bar blocks under the action of the required supporting force of the last bar block until the required supporting force of the last bar block of the second area is calculated;
calculating a first point safety factor for the third region based on the holding power required for the last bar of the second region:
wherein F is c A first point safety factor representing a third region or c-th bar, R c Representing the anti-slip force, T, of the c-th bar c Representing the sliding force of the c-th bar;
and (3) changing the second shear strength parameter to enable the first point safety coefficient to be equal to a second threshold value, and taking the second shear strength parameter when the first point safety coefficient is equal to the second threshold value as the sliding surface shear strength parameter.
2. The method for inverting the sliding surface shear strength parameter according to claim 1, wherein the calculating the first point safety factor of the third region according to the first shear strength parameter, the second shear strength parameter, the first matrix and the first threshold value, and the changing the second shear strength parameter to make the first point safety factor equal to the second threshold value, takes the second shear strength parameter when the first point safety factor is equal to the second threshold value as the sliding surface shear strength parameter, specifically comprises:
calculating a first point safety coefficient of the third area according to the second shear strength parameter, the first shear strength parameter and a first matrix;
judging whether the first point safety coefficient of the third area is equal to the second threshold value or not:
if the two parameters are equal, taking the current value of the second shear strength parameter as the sliding surface shear strength parameter;
if the first and second points are not equal, changing the value of the second shear strength parameter, and recalculating the first point safety coefficient of the third area.
3. The sliding surface shear strength parameter inversion method according to claim 1, wherein the calculating based on the first point safety factor obtains a second point safety factor of the current bar and a supporting force required by the current bar, specifically comprising:
judging the area to which the current bar belongs:
if the current bar block belongs to the first area, enabling a preset threshold value to be a first threshold value; if the current bar block belongs to the second area, enabling the preset threshold value to be a second threshold value;
judging whether the first point safety coefficient is larger than or equal to a preset threshold value or not:
if yes, the second point safety coefficient is equal to the first point safety coefficient, and the supporting force required by the current bar is 0;
if not, the second point safety coefficient is enabled to be equal to a preset threshold value, and the supporting force required by the current bar block is calculated under the condition of meeting the static balance.
4. A slip plane shear strength parameter inversion apparatus, comprising:
a first dividing module: vertically dividing a sliding body in a landslide to obtain a plurality of strips;
a second dividing module: dividing a first area, a second area and a third area from the plurality of strips, wherein the second area is a strong deformation area;
a first acquisition module: obtaining the vertex coordinates of each bar in the first area, the second area and the third area, calculating the geometric parameters of each bar according to the vertex coordinates, and forming a first matrix by the geometric parameters of each bar;
and a second acquisition module: determining a first shear strength parameter, a second shear strength parameter, a first threshold value of a first regional point safety coefficient and a second threshold value of a second regional point safety coefficient of the landslide;
and an inversion module: calculating a first point safety factor for the third region based on the first shear strength parameter, the second shear strength parameter, the first matrix, and the first threshold value, comprising:
obtaining geometric parameters of a current bar in a first matrix, wherein the current bar is any one of a first area and a second area;
calculating a first point safety coefficient of the current bar according to the geometric parameter of the current bar, the current value of the second shear strength parameter and the first shear strength parameter:
wherein F is k1 A first point safety factor representing the kth bar, R k Representing the anti-slip force, T, of the kth bar k Representing the sliding down force of the kth bar;
calculating a second point safety coefficient of the current bar and a supporting force required by the current bar based on the first point safety coefficient;
calculating a second point safety coefficient and a required supporting force of a next bar of the current bar under the action of the required supporting force of the current bar;
sequentially calculating a second point safety coefficient and a required supporting force of the rest bar blocks under the action of the required supporting force of the last bar block until the required supporting force of the last bar block of the second area is calculated;
calculating a first point safety factor for the third region based on the holding power required for the last bar of the second region:
wherein F is c A first point safety factor representing a third region or c-th bar, R c Representing the anti-slip force, T, of the c-th bar c Representing the sliding force of the c-th bar;
and (3) changing the second shear strength parameter to enable the first point safety coefficient to be equal to a second threshold value, and taking the second shear strength parameter when the first point safety coefficient is equal to the second threshold value as the sliding surface shear strength parameter.
5. The slip shear strength parameter inversion device of claim 4, wherein the inversion module specifically comprises:
a third calculation unit: calculating a first point safety coefficient of the third area according to the second shear strength parameter, the first shear strength parameter and a first matrix;
a second judgment unit: judging whether the first point safety coefficient of the third area is equal to the second threshold value or not:
if the two parameters are equal, taking the current value of the second shear strength parameter as the sliding surface shear strength parameter;
if the first and second points are not equal, changing the value of the second shear strength parameter, and recalculating the first point safety coefficient of the third area.
6. The slip shear strength parameter inversion device of claim 4, wherein the fifth calculation unit comprises:
a third judgment unit: judging the area to which the current bar belongs:
if the current bar block belongs to the first area, enabling a preset threshold value to be a first threshold value; if the current bar block belongs to the second area, enabling the preset threshold value to be a second threshold value;
fourth judgment unit: judging whether the first point safety coefficient is larger than or equal to a preset threshold value or not:
if yes, the second point safety coefficient is equal to the first point safety coefficient, and the supporting force required by the current bar is 0;
if not, the second point safety coefficient is enabled to be equal to a preset threshold value, and the supporting force required by the current bar block is calculated under the condition of meeting the static balance.
7. A slip plane shear strength parameter inversion apparatus comprising:
a memory for storing a computer program;
a processor for carrying out the steps of the slide shear strength parameter inversion method according to any one of claims 1 to 3 when executing the computer program.
8. A computer-readable storage medium, characterized by: a computer program stored on a computer readable storage medium, which when executed by a processor, implements the steps of the sliding surface shear strength parameter inversion method according to any one of claims 1 to 3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310150425.XA CN115879325B (en) | 2023-02-22 | 2023-02-22 | Sliding surface shear strength parameter inversion method, device, equipment and readable storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310150425.XA CN115879325B (en) | 2023-02-22 | 2023-02-22 | Sliding surface shear strength parameter inversion method, device, equipment and readable storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115879325A CN115879325A (en) | 2023-03-31 |
| CN115879325B true CN115879325B (en) | 2023-06-09 |
Family
ID=85761531
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310150425.XA Active CN115879325B (en) | 2023-02-22 | 2023-02-22 | Sliding surface shear strength parameter inversion method, device, equipment and readable storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115879325B (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103306672A (en) * | 2013-05-24 | 2013-09-18 | 中国石油大学(北京) | Method for predicting abrasiveness of shaly stratum in different drilling directions |
| CN115374528A (en) * | 2022-10-24 | 2022-11-22 | 西南交通大学 | Slope safety analysis method, system and equipment and readable storage medium |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109726444A (en) * | 2018-12-07 | 2019-05-07 | 青岛理工大学 | The inverting measuring method of rain-induced landslide shear strength parameter |
| CN114004055B (en) * | 2021-07-23 | 2023-02-17 | 中国电建集团西北勘测设计研究院有限公司 | Slope shear strength parameter inversion analysis method based on equivalent soil pressure effect |
| CN114692441B (en) * | 2022-02-17 | 2023-05-09 | 成都理工大学 | Loess landslide stability prediction method, electronic equipment and storage medium |
| CN115546433B (en) * | 2022-11-23 | 2023-04-07 | 西南交通大学 | Slope softening area prediction method, system and equipment and readable storage medium |
-
2023
- 2023-02-22 CN CN202310150425.XA patent/CN115879325B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103306672A (en) * | 2013-05-24 | 2013-09-18 | 中国石油大学(北京) | Method for predicting abrasiveness of shaly stratum in different drilling directions |
| CN115374528A (en) * | 2022-10-24 | 2022-11-22 | 西南交通大学 | Slope safety analysis method, system and equipment and readable storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115879325A (en) | 2023-03-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Dyson et al. | Prediction and classification for finite element slope stability analysis by random field comparison | |
| JP6784239B2 (en) | Face evaluation support system, face evaluation support method and face evaluation support program | |
| Ducea et al. | High-volume magmatic events in subduction systems | |
| CN104424648B (en) | Method for tracing object and equipment | |
| CN109382822B (en) | Simulation device and simulation method for robot system | |
| Sarkar et al. | Rock mass classification and slope stability assessment of road cut slopes in Garhwal Himalaya, India | |
| CN105069248B (en) | Large landslide polylith sliding scale | |
| JP2019012037A (en) | Material characteristic estimation device and material characteristic estimation method | |
| CN115879325B (en) | Sliding surface shear strength parameter inversion method, device, equipment and readable storage medium | |
| CN106021770A (en) | Stability analysis method and device for airplane structure rectangular flat plate | |
| CN115546433B (en) | Slope softening area prediction method, system and equipment and readable storage medium | |
| CN113221371A (en) | Method and device for determining critical sliding surface of side slope and terminal equipment | |
| JP7348575B2 (en) | Deterioration detection device, deterioration detection system, deterioration detection method, and program | |
| JP6438169B1 (en) | Earthquake prediction system and earthquake prediction method | |
| CN115858996B (en) | Safety coefficient calculation method, device, equipment and medium based on sectional landslide | |
| JP2016106684A (en) | Skeleton model creating device, method and program | |
| US20150084889A1 (en) | Stroke processing device, stroke processing method, and computer program product | |
| JP5429886B2 (en) | Interface particle determination method and apparatus in particle method, and interface particle determination program | |
| Quiroz‐Ramírez et al. | Evaluation of the intensity measure approach in performance‐based earthquake engineering with simulated ground motions | |
| CN115935490B (en) | Anti-skid resistance calculation method, device, equipment and medium based on point safety coefficient | |
| Bahari et al. | Nonlinear estimation model of minimum void ratio for sand–silt mixtures | |
| CN113642071A (en) | Stability determination method and device for rock drilling chamber and computer readable storage medium | |
| WO2016142965A1 (en) | Video processing device, video processing method, and storage medium storing video processing program | |
| US10545262B2 (en) | Method for determining a proportion cube | |
| Zhang et al. | Elimination of boundary effect for laminated overburden model in pillar stability analysis |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |