CN113609558A - Slope influence factor determination method and device based on uniform design - Google Patents
Slope influence factor determination method and device based on uniform design Download PDFInfo
- Publication number
- CN113609558A CN113609558A CN202110845202.6A CN202110845202A CN113609558A CN 113609558 A CN113609558 A CN 113609558A CN 202110845202 A CN202110845202 A CN 202110845202A CN 113609558 A CN113609558 A CN 113609558A
- Authority
- CN
- China
- Prior art keywords
- test
- factors
- slope
- determining
- uniform design
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Evolutionary Computation (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Structural Engineering (AREA)
- Civil Engineering (AREA)
- Architecture (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Probability & Statistics with Applications (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The application provides a slope influence factor determination method and device based on uniform design. The method comprises the following steps: determining a test factor; determining the value range of each test factor based on the test factors and the slope structure model to be tested; constructing a uniform design table based on the test factors and the value range of each test factor; testing the test factors through the uniform design table to obtain a stability coefficient corresponding to each test number; and determining slope influence sensitive factors in the test factors based on the stability coefficient corresponding to each test number. The uniform design can well carry out sensitivity analysis on all factors influencing the slope stability, is particularly suitable for multi-factor and multi-level conditions, can greatly reduce the test times and the calculated amount, and also enables test points to be more uniformly distributed in the test range, thereby achieving the ideal effect.
Description
Technical Field
The application relates to the technical field of data processing, in particular to a slope influence factor determination method and device based on uniform design.
Background
The construction of infrastructures such as roads, railways and water conservancy forms a large number of side slopes, and the increasing number of side slope accidents cause huge loss to people. The slope damage is formed under the multi-factor composite action, the influence degree of each factor on the slope stability determines a sliding mechanism and a damage mode, and the sensitivity analysis of the slope stability influence factors has important significance on slope management, monitoring and forecasting.
Slope stability sensitivity analysis is a current hotspot problem, and many experts and scholars carry out deeper research on the slope stability sensitivity analysis, and the slope stability sensitivity analysis is summarized into a single-factor analysis method and a multi-factor analysis method.
The single-factor analysis method essentially selects an index value, changes one factor, simultaneously assumes that other factors are kept unchanged, compares the change of the reference value with the factor, and then intuitively reflects the sensitivity of each factor in a reference value-factor change curve. The method can intuitively reflect the influence of each factor on the reference value, but needs certain assumption premise and is not consistent with the actual situation.
The common characteristics of the methods are that the probability statistical principle and the electronic computer technology are adopted to reduce the test times and the calculation workload to a certain extent, but when the factors and the levels are more, the test times and the calculation workload are still large, and the difficulty is brought to the test and the numerical analysis.
Disclosure of Invention
The embodiment of the application aims to provide a slope influence factor determining method and device based on uniform design so as to solve the problem that the existing test times and calculation amount are large and difficulty is brought to test and numerical analysis.
The invention is realized by the following steps:
in a first aspect, an embodiment of the present application provides a slope influence factor determination method based on uniform design, including: determining a test factor; wherein the test factors are derived from a slope structure model to be tested; determining the value range of each test factor based on the test factors and the slope structure model to be tested; constructing a uniform design table based on the test factors and the value range of each test factor; the uniform design table comprises test numbers and numerical values of test factors in each test number; testing the test factors through the uniform design table to obtain a stability coefficient corresponding to each test number; and determining slope influence sensitive factors in the test factors based on the stability coefficient corresponding to each test number.
In the embodiment of the application, after the test factors and the value range of each test factor are determined, a uniform design table is constructed; the uniform design can well carry out sensitivity analysis on all factors influencing the slope stability, is particularly suitable for multi-factor and multi-level conditions, can greatly reduce the test times and the calculated amount, and also enables test points to be more uniformly distributed in the test range, thereby achieving the ideal effect. Meanwhile, the value range of each test factor is determined by the test factors and the slope structure model to be tested, so that the correctness and the rationality of a later analysis conclusion can be improved.
With reference to the technical solution provided by the first aspect, in some possible implementations, the uniform design table is constructed through the following steps; determining the test times n; obtain vector set h ═ h (h)1,h2,...hm) (ii) a Wherein h is1,h2,...hmAre all numbers less than n, and the greatest common divisor with n is 1; generating a jth column of the uniform design table until the uniform design table is obtained; wherein, the expression in the jth column is: u. ofij＝(i×hj)[mod n](ii) a Wherein, [ modn]Representing a congruence operation; if i × hjIf n is exceeded, a preset multiple of n is subtracted to make the fall within [1, n ]]To (1); u. ofijGenerating through a recursion algorithm; the uniform design is denoted as Un(nm) (ii) a The number of m corresponding to the determined n and h is obtained by an Euler function E (n); when n is a prime number, e (n) n-1; when n is a prime power, i.e. n ═ plWhen the temperature of the water is higher than the set temperature,
p is prime number, l is positive integer; when n is neither a prime nor a prime power, n is expressed as a square power product of different prime numbers; namely, it is
At this time, the process of the present invention,
p1,p2,L pmare different prime numbers,/1,l2,L lmAre integers.
In the embodiment of the application, the number n of tests is determined first, then a set of numbers which are smaller than n and have a greatest common divisor with n of 1 is obtained, and then the formula u is usedij＝(i×hj)[mod n]And generating each column in the uniform design table, and finally obtaining the uniform design table. The uniform design table obtained in this way allows one and only one test to be done for each level of each factor; the test points for any two factors can be marked on a planar grid with a column for each row and a columnThe uniformity of a test scheme consisting of any two columns of a uniform design table with only one test point is generally not equivalent; when the level of the factor increases, the number of trials increases as the number of levels increases.
With reference to the technical solution provided by the first aspect, in some possible implementations, the test factors include: volume weight, internal friction angle, cohesion, slope angle and seismic coefficient.
With reference to the technical solution provided by the first aspect, in some possible implementation manners, the number of times of the test n is twice as large as the number of the test factors.
In the embodiment of the application, the number of the test times n is twice of the number of the test factors, so that the test points are distributed more uniformly in the test range while the reduction of the test times and the calculation amount is ensured, and a more ideal effect is achieved.
With reference to the technical solution provided by the first aspect, in some possible implementation manners, the testing the test factors through the uniform design table to obtain the stability coefficient corresponding to each test number includes: testing the test factors through the uniform design table; obtaining a stability coefficient corresponding to each test number by a simplified Bischu method; wherein the stability factor is expressed as:
The stability factor; c. CiThe cohesion corresponding to the ith test; phi is aiThe internal friction angle corresponding to the ith test; wiThe block self weight is the ith test; u. ofiRepresents the pore water pressure; biThe width of the block for the ith test; beta is aiThe slope angle corresponding to the ith test; qHiHorizontal seismic force; kHIs a horizontal seismic coefficient; cZIs a comprehensive influence coefficient; alpha is alphaiIs the seismic acceleration distribution coefficient.
In the embodiment of the application, the stability coefficient corresponding to each test number is obtained by simplifying the Bishou method, so that the stability coefficient is more accurate, and the slope stability under each test factor can be more intuitively and reliably shown.
With reference to the technical solution provided by the first aspect, in some possible implementation manners, the determining, based on the stability coefficient corresponding to each test number, a slope influence sensitive factor in the test factors includes: regression analysis is carried out on the test factors by adopting a backspacing method, and a regression equation is established; carrying out significance test on the regression equation; and determining the slope influence sensitive factors in the test factors based on the result of the significance test.
In the embodiment of the application, regression analysis is carried out on the test factors by adopting a backspacing method, and a regression equation is established; carrying out significance test on the regression equation; based on the result of the significance test, slope influence sensitive factors in the test factors are determined, and by means of the method, slope stability under each test factor can be effectively analyzed, so that accurate sensitive factors influencing slope stability sensitivity are determined subsequently.
With reference to the technical solution provided by the first aspect, in some possible implementation manners, the performing significance test on the regression equation includes: carrying out significance test on the regression equation to determine the insignificant factors in the test factors; removing the non-significant factors and reconstructing a new regression equation; performing significance test on the new regression equation; correspondingly, the determining of the slope influence sensitive factors in the test factors based on the result of the significance test comprises: and determining slope influence sensitive factors in the test factors based on the result of the significance test on the new regression equation.
In the embodiment of the application, the non-significant factors in the test factors are determined by performing significance test on the regression equation; removing the non-significant factors and reconstructing a new regression equation; and the significance test is carried out on the new regression equation, so that the reliability of the significance test is further improved.
With reference to the technical solution provided by the first aspect, in some possible implementation manners, the determining the test factors includes: establishing a slope numerical value to construct a test slope structure model; wherein the model is a homogeneous soil slope; and determining the test factors based on all numerical values of the slope structure model to be tested.
In the embodiment of the application, a slope numerical value is established to construct a test slope structure model, and then the test slope structure model passes through, and then test factors are determined based on various numerical values of the slope structure model to be tested; the model is a homogeneous soil slope, the test is carried out by constructing the model and determining test factors based on the model, the method has important significance for the disaster mechanism and prevention and control of the side slope, and especially has important theoretical value and practical value for the research and application of the prevention and control of the bank side slope in the three gorges reservoir area.
In a second aspect, an embodiment of the present application provides a slope influence factor determination device based on a uniform design, including: the first determining module is used for determining a test factor; wherein the test factors are derived from a slope structure model to be tested; the second determination module is used for determining the value range of each test factor based on the test factors and the slope structure to be tested; the construction module is used for constructing a uniform design table based on the test factors and the value range of each test factor; the uniform design table comprises test numbers and numerical values of test factors in each test number; the processing module is used for testing the test factors through the uniform design table to obtain a stability coefficient corresponding to each test number; and the third determining module is used for determining the slope influence sensitive factors in the test factors based on the stability coefficient corresponding to each test number.
In a third aspect, an embodiment of the present application provides an electronic device, including: a processor and a memory, the processor and the memory connected; the memory is used for storing programs; the processor is configured to invoke a program stored in the memory to perform a method as provided in the above-described first aspect embodiment and/or in combination with some possible implementations of the above-described first aspect embodiment.
In a fourth aspect, embodiments of the present application provide a storage medium having stored thereon a computer program, which, when executed by a processor, performs a method as provided in the above-described first aspect embodiment and/or in connection with some possible implementations of the above-described first aspect embodiment.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are required to be used in the embodiments of the present application will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present application and therefore should not be considered as limiting the scope, and that those skilled in the art can also obtain other related drawings based on the drawings without inventive efforts.
Fig. 1 is a block diagram of an electronic device according to an embodiment of the present disclosure.
Fig. 2 is a flowchart of steps of a slope influence factor determination method based on uniform design according to an embodiment of the present application.
Fig. 3 is a schematic diagram of a slope structure model to be tested according to an embodiment of the present application.
Fig. 4 is a block diagram of a slope influence factor determination device based on uniform design according to an embodiment of the present application.
Icon: 100-an electronic device; 10-a processor; 11-a memory; 200-a slope influencing factor determining device based on uniform design; 201-a first determination module; 202-a second determination module; 203-a building block; 204-a processing module; 205-third determination module.
Detailed Description
The technical solutions in the embodiments of the present application will be described below with reference to the drawings in the embodiments of the present application.
In view of the problems that the conventional slope stability sensitivity analysis has a large number of factors and levels, the number of tests and the calculation amount are still large, and the tests and the numerical analysis are difficult, the inventor of the present application has made research and proposes the following embodiments to solve the problems.
Referring to fig. 1, a schematic structural block diagram of an electronic device 100 for a slope influence factor determination method and apparatus based on uniform design according to an embodiment of the present application is provided. In the embodiment of the present application, the electronic Device 100 may be, but is not limited to, a Personal Computer (PC), a smart phone, a tablet Computer, a Personal Digital Assistant (PDA), a Mobile Internet Device (MID), and the like. Structurally, electronic device 100 may include a processor 10 and a memory 11.
The processor 10 and the memory 11 are electrically connected directly or indirectly to enable data transmission or interaction, for example, the components may be electrically connected to each other via one or more communication buses or signal lines. The slope influencing factor determining means based on the uniform design comprises at least one software module which can be stored in the form of software or Firmware (Firmware) in the memory 11 or solidified in an Operating System (OS) of the electronic device 100. The processor 10 is configured to execute executable modules stored in the memory 11, such as software functional modules and computer programs included in the slope influencing factor determining apparatus based on uniform design, so as to implement the slope influencing factor determining method based on uniform design. The processor 10 may execute the computer program upon receiving the execution instruction.
The processor 10 may be an integrated circuit chip having signal processing capabilities. The Processor 10 may also be a general-purpose Processor, for example, a Central Processing Unit (CPU), a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a discrete gate or transistor logic device, or a discrete hardware component, which may implement or execute the methods, steps, and logic blocks disclosed in the embodiments of the present Application. Further, a general purpose processor may be a microprocessor or any conventional processor or the like.
The Memory 11 may be, but is not limited to, a Random Access Memory (RAM), a Read Only Memory (ROM), a Programmable Read-Only Memory (PROM), an Erasable Programmable Read-Only Memory (EPROM), and an electrically Erasable Programmable Read-Only Memory (EEPROM). The memory 11 is used for storing a program, and the processor 10 executes the program after receiving the execution instruction.
It should be noted that the structure shown in fig. 1 is only an illustration, and the electronic device 100 provided in the embodiment of the present application may also have fewer or more components than those shown in fig. 1, or have a different configuration than that shown in fig. 1. Further, the components shown in fig. 1 may be implemented by software, hardware, or a combination thereof.
Referring to fig. 2, fig. 2 is a flowchart illustrating steps of a slope influence factor determining method based on uniform design according to an embodiment of the present application, where the method is applied to the electronic device 100 shown in fig. 1. It should be noted that, the slope influence factor determining method based on uniform design provided in the embodiment of the present application is not limited by the sequence shown in fig. 2 and the following, and the method includes: step S101-step S105.
Step S101: determining a test factor; wherein, the test factors come from the slope structure model to be tested.
Step S102: and determining the value range of each test factor based on the test factors and the slope structure model to be tested.
Step S103: constructing a uniform design table based on the test factors and the value range of each test factor; wherein, the uniform design table comprises the test numbers and the numerical values of the test factors in each test number.
Step S104: and testing the test factors through a uniform design table to obtain the stability coefficient corresponding to each test number.
Step S105: and determining the slope influence sensitive factors in the test factors based on the stability coefficient corresponding to each test number.
In the embodiment of the application, after the test factors and the value range of each test factor are determined, a uniform design table is constructed; the uniform design can well carry out sensitivity analysis on all factors influencing the slope stability, is particularly suitable for multi-factor and multi-level conditions, can greatly reduce the test times and the calculated amount, and also enables test points to be more uniformly distributed in the test range, thereby achieving the ideal effect. Meanwhile, the value range of each test factor is determined by the test factors and the slope structure model to be tested, so that the correctness and the rationality of a later analysis conclusion can be improved.
For ease of understanding, the uniform design table provided in the embodiments of the present application will be described first. Specifically, the uniform design table is constructed by the following steps: determining the test times n; obtain vector set h ═ h (h)1,h2,...hm) (ii) a Wherein h is1,h2,...hmAre all numbers less than n, and the greatest common divisor with n is 1; and generating a jth column of the uniform design table until the uniform design table is obtained.
Wherein, the expression in the jth column is:
uij＝(i×hj)[mod n] (1)
in equation (1), [ modn]Representing a congruence operation; if i × hjIf n is exceeded, a preset multiple of n is subtracted to make the fall within [1, n ]]Among them.
uijGenerating through a recursion algorithm; the uniform design is expressed as Un(nm) (ii) a The number of m corresponding to n and h is determined by an Euler function E (n). Specifically, the euler function e (n) is obtained under the following conditions:
when n is a prime number, e (n) is n-1.
When n is a prime power, i.e. n ═ plWhen the temperature of the water is higher than the set temperature,
p is prime number and l is positive integer.
When n is neither a prime nor a prime power, n is expressed as a square power product of different prime numbers; namely, it is
At this time, the process of the present invention,
p1,p2,L pmare different prime numbers,/1,l2,L lmAre integers.
It should be noted that, in the above three cases, n is the prime number, which is the best case, at most n-1 columns can be obtained, i.e. n-1 factors can be arranged for testing, and when n is a non-prime number, n-1 columns are not possible in the structure of the uniform design table, and even the possible number of columns is less. Therefore, U can be replacedn(nm) Is constructed by removing the last row of
The method of (1). Due to the fact that
Than Un(nm) Have better uniformity and are therefore generally considered important
In addition, each uniform table corresponds to a usage table. Using representation to guarantee slave Un(nm) Or
The s columns selected in (1) have good uniformity, and uniformity is measured by deviation. Let x1,x2,L xnIs CmN points in (c), any vector x ═ x1,x2,L xm)∈CmIs denoted by v (x) ═ x1L xmIs a rectangle [0, x]Volume of (1), nxIs x1,x2,L xnFalls into [0, x ]]The number of points of
Called point set { x1,x2,L xnAt CmDeviation in (2).
The above steps are exemplified below.
In step S101, a test factor is determined, which specifically includes: establishing a slope numerical value to construct a test slope structure model; wherein the model is a homogeneous soil slope; and determining test factors based on various numerical values of the slope structure model to be tested.
In the embodiment of the present application, the test factors may include any one or more of volume weight, internal friction angle, cohesive force, slope angle, and seismic coefficient, and of course, in other embodiments, the test factors may also be any other factors, such as water content, and the like, and the present application is not limited thereto.
It can be understood that, in this step, a corresponding test slope structure model is constructed according to the slope to be tested, and then a corresponding numerical value is given. And then, determining test factors according to various numerical values of the constructed slope structure model to be tested. That is, the factors to be tested can be screened out according to various numerical values of the constructed slope structure model to be tested, or the factors corresponding to all the numerical values can be included, and the method is not limited in the application.
In summary, in the embodiment of the application, a slope numerical value is established to construct a test slope structure model, and then the test factors are determined based on various numerical values of the slope structure model to be tested; the model is a homogeneous soil slope, the test is carried out by constructing the model and determining test factors based on the model, the method has important significance for the disaster mechanism and prevention and control of the side slope, and especially has important theoretical value and practical value for the research and application of the prevention and control of the bank side slope in the three gorges reservoir area.
Optionally, in the embodiment of the present application, the number of test times n is twice of the number of test factors, so that while the reduction of the number of test times and the reduction of the calculated amount are ensured, the test points are distributed more uniformly in the test range, and a more ideal effect is achieved.
It is understood that in other embodiments, the number of tests n is three times the number of test factors, and the application is not limited thereto.
In step S102, a value range of each test factor is determined based on the test factors and the slope structure model to be tested, and the value range may be determined by floating up and down a preset range respectively with a value given to the test slope structure model as a reference value.
In step S104, the test factors are tested through the uniform design table to obtain the stability coefficient corresponding to each test number, including: testing the test factors through a uniform design table; obtaining a stability coefficient corresponding to each test number by a simplified Bisshop method (Bisshop);
wherein the stability coefficient is expressed as:
in the formula (2), K is a stability coefficient; c. CiThe cohesion corresponding to the ith test; phi is aiThe internal friction angle corresponding to the ith test; wiThe block self weight is the ith test; u. ofiRepresents the pore water pressure; biThe width of the block for the ith test; beta is aiThe slope angle corresponding to the ith test; qHiHorizontal seismic force; kHIs a horizontal seismic coefficient; cZIs a comprehensive influence coefficient; alpha is alphaiIs the seismic acceleration distribution coefficient.
The biotransformation method is an arc sliding analysis method in which the interaction force between soil strips is considered in the soil slope stability analysis. Still based on the extreme balance principle, the slip-crack soil body is used as a rigid body to rotate around the circle center, the sliding force and the anti-sliding force of the slip-crack soil body are calculated in strips, finally, the stability safety coefficient is calculated, the interaction force among the soil strips is considered during calculation, and the method is an improved arc sliding method.
In the embodiment of the application, the stability coefficient corresponding to each test number is obtained by simplifying the Bishou method, so that the stability coefficient is more accurate, and the slope stability under each test factor can be more intuitively and reliably shown.
In step S105, based on the stability coefficient corresponding to each test number, determining a slope influence sensitive factor in the test factors, which specifically includes: regression analysis is carried out on the test factors by adopting a backspacing method, and a regression equation is established; carrying out significance test on the regression equation; and determining the slope influence sensitive factors in the test factors based on the result of the significance test.
In the embodiment of the application, regression analysis is carried out on the test factors by adopting a backspacing method, and a regression equation is established; carrying out significance test on the regression equation; based on the result of the significance test, slope influence sensitive factors in the test factors are determined, and by means of the method, slope stability under each test factor can be effectively analyzed, so that accurate sensitive factors influencing slope stability sensitivity are determined subsequently.
Optionally, the performing the significance test on the regression equation specifically includes: carrying out significance test on the regression equation to determine insignificant factors in the test factors; removing the non-significant factors and reconstructing a new regression equation; carrying out significance test on the new regression equation; correspondingly, the slope influence sensitive factors in the test factors are determined based on the results of the significance test, and include: and determining slope influence sensitive factors in the test factors based on the result of the significance test on the new regression equation.
The expression of the regression equation above is:
y＝b0+b1×γ+b2×φ+b3×c+b4×β+b5×α (3)
in the formula (3), b0＝1.75、b1＝-2.5×10-2、b2＝2.14×10-2、b3＝4.62×10-2、b4＝-4.34×10-2、b5＝-3.18×10-4(ii) a Alpha is significance level and takes 0.05.
In the embodiment of the application, the non-significant factors in the test factors are determined by performing significance test on the regression equation; removing the non-significant factors and reconstructing a new regression equation; and the significance test is carried out on the new regression equation, so that the reliability of the significance test is further improved.
For the convenience of understanding, a slope influence factor determination method based on uniform design provided by the embodiment of the present application is described below with a specific example.
Referring to fig. 3, fig. 3 is a schematic diagram illustrating a test slope structure model constructed by establishing slope values. Wherein the slope height h is 20m, the soil volume weight gamma is 20kN/m3(ii) a The internal friction angle phi is 17 degrees, the cohesive force c is 42kPa, and the seismic coefficient alpha is 0.05 g.
And then, determining the value range of each test factor based on the test factors and the slope structure model to be tested. Assuming that the above parameters are used as reference values, the range is a preset value which is floated up and down.
Illustratively, the value range of each test factor is:
the soil volume weight gamma is 18-27 (kN/m)3)；
The range of the internal friction angle phi is 15-33 (°);
the cohesive force c is in the range of 27 to 54 (kPa);
the range of the slope angle beta is 31-49 (°); the reference value of the slope angle can be calculated according to the graph;
the range of the seismic coefficient alpha is 10-100 (cm/s)2)。
After the test factors and the value range of each test factor are determined, a uniform design table is constructed
The number of times n of the test is twice of the number of the test factors, namely the number of times of the test is 10, and then the test factors are tested through a uniform design table to obtain the stability coefficient corresponding to each test number. (see tables 1 to 3 for details)
TABLE 1 Uniform design Table
TABLE 2 Uniform design Table
Using table
Table 3 uniform design calculation results
And finally, determining the slope influence sensitive factors in the test factors based on the stability coefficient corresponding to each test number. In the embodiment of the application, the stability coefficient corresponding to each test number can be obtained by simplifying the Bishou method, and then regression analysis is performed on the test factors by adopting a backspacing method to establish a regression equation; and (4) carrying out significance test on the regression equation, and further determining the slope influence sensitive factors in the test factors.
Specifically, regression analysis was performed by a regression method, and the significance level α was taken to be 0.05.
The expression of the regression equation is:
y＝b0+b1×γ+b2×φ+b3×c+b4×β+b5×α (4)
in the formula (4), b0＝1.75、b1＝-2.5×10-2、b2＝2.14×10-2、b3＝4.62×10-2、b4＝-4.34×10-2、b5＝-3.18×10-4(ii) a The formula (4) is the same as the formula (3).
The regression equation significance test can be referred to table 4.
TABLE 4 analysis of variables
The regression equation is significant with the sample volume N being 10, the test value Ft being 52.20, the threshold value F (0.05,5,4) being 6.265, Ft > F (0.05,5, 4).
F, checking value: f (1) ═ 5.992, F (2) ═ 17.50, F (3) ═ 184.3, F (4) ═ 72.27, F (5) ═ 9.707 × 10-2, the contribution of each equation term to regression is c, β, Φ, γ, α in turn, the contribution of equation term α to regression is minimal, the significance test is performed on the contribution, the critical value F (0.05,1,4) ═ 7.709, F (5) ≦ F (0.05,1,4), the equation term is not significant and needs to be eliminated.
Eliminating an insignificant equation item, and establishing a new regression equation:
y＝b0+b1×γ+b2×φ+b3×c+b4×β (5)
in the formula (5), b0＝1.73、b1＝-2.42×10-2、b2＝2.10×10-2、b3＝4.65×10-2、b4＝-4.38×10-2；
The regression equation significance test can be referred to table 5.
TABLE 5 analysis of variables
The regression equation is significant with the sample volume N being 10, the test value Ft being 79.61, the threshold value F (0.05,4,5) being 5.192, Ft > F (0.05,4, 5).
F, checking value: f (1) ═ 7.312, F (2) ═ 21.95, F (3) ═ 242.6, F (4) ═ 95.81, and the contribution of each equation term to regression is in turn: c. β, Φ, γ, the equation term γ, which has the smallest contribution to the regression, was tested for significance, with the cutoff value F (0.05,1,5) ═ 6.608, F (1) > F (0.05,1,5), this equation term being significant.
The residual analysis can be referred to table 6.
Table 6 residual error analysis table
According to the regression analysis result, the earthquake action alpha is an insignificant factor in the areas with the earthquake fortification intensity below 7 degrees, the soil volume weight gamma, the internal friction angle phi, the cohesive force c and the slope angle beta are significant factors, and the significance is as follows: cohesive force c, slope angle beta, internal friction angle phi and soil volume weight gamma.
Referring to fig. 4, based on the same inventive concept, an embodiment of the present application further provides a slope influence factor determining apparatus 200 based on uniform design, the apparatus including: a first determination module 201, a second determination module 202, a construction module 203, a processing module 204, and a third determination module.
A first determining module 201, configured to determine a test factor; wherein the test factors are derived from the slope structure model to be tested.
A second determining module 202, configured to determine a value range of each test factor based on the test factor and the slope structure to be tested.
A construction module 203, configured to construct a uniform design table based on the test factors and the value range of each test factor; wherein, the uniform design table comprises test numbers and numerical values of test factors in each test number.
And the processing module 204 is configured to test the test factors through the uniform design table to obtain a stability coefficient corresponding to each test number.
And a third determining module 205, configured to determine, based on the stability coefficient corresponding to each test number, a slope influence sensitive factor in the test factors.
Optionally, the constructing module 203 is specifically configured to determine the number n of tests; obtain vector set h ═ h (h)1,h2,...hm) (ii) a Wherein h is1,h2,...hmAre all numbers less than n, and the greatest common divisor with n is 1; generating a jth column of the uniform design TableUntil the uniform design table is obtained; wherein, the expression in the jth column is: u. ofij＝(i×hj)[modn](ii) a Wherein, [ modn]Representing a congruence operation; if i × hjIf n is exceeded, a preset multiple of n is subtracted to make the fall within [1, n ]]To (1); u. ofijGenerating through a recursion algorithm; the uniform design is denoted as Un(nm) (ii) a The number of m corresponding to the determined n and h is obtained by an Euler function E (n); when n is a prime number, e (n) n-1; when n is a prime power, i.e. n ═ plWhen the temperature of the water is higher than the set temperature,
p is prime number, l is positive integer; when n is neither a prime nor a prime power, n is expressed as a square power product of different prime numbers; namely, it is
At this time, the process of the present invention,
p1,p2,L pmare different prime numbers,/1,l2,L lmAre integers.
Optionally, the processing module 204 is specifically configured to test the test factors through the uniform design table; and obtaining the stability coefficient corresponding to each test number by a simplified Bischu method.
Optionally, the third determining module 205 is specifically configured to perform regression analysis on the test factors by using a regression method to establish a regression equation; carrying out significance test on the regression equation; and determining the slope influence sensitive factors in the test factors based on the result of the significance test.
Optionally, the third determining module 205 is further specifically configured to perform significance test on the regression equation to determine insignificant factors in the test factors; removing the non-significant factors and reconstructing a new regression equation; performing significance test on the new regression equation; and determining slope influence sensitive factors in the test factors based on the result of the significance test on the new regression equation.
Optionally, the first determining module 201 is specifically configured to establish a slope numerical value to construct a test slope structural model; wherein the model is a homogeneous soil slope; and determining the test factors based on all numerical values of the slope structure model to be tested.
It should be noted that, as those skilled in the art can clearly understand, for convenience and brevity of description, the specific working processes of the above-described systems, apparatuses and units may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.
Based on the same inventive concept, embodiments of the present application further provide a computer-readable storage medium, on which a computer program is stored, and when the computer program is executed, the computer program performs the methods provided in the above embodiments.
The storage medium may be any available medium that can be accessed by a computer or a data storage device including one or more integrated servers, data centers, and the like. The usable medium may be a magnetic medium (e.g., floppy Disk, hard Disk, magnetic tape), an optical medium (e.g., DVD), or a semiconductor medium (e.g., Solid State Disk (SSD)), among others.
In the embodiments provided in the present application, it should be understood that the disclosed apparatus and method may be implemented in other ways. The above-described embodiments of the apparatus are merely illustrative, and for example, the division of the units is only one logical division, and there may be other divisions when actually implemented, and for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection of devices or units through some communication interfaces, and may be in an electrical, mechanical or other form.
In addition, units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
Furthermore, the functional modules in the embodiments of the present application may be integrated together to form an independent part, or each module may exist separately, or two or more modules may be integrated to form an independent part.
In this document, relational terms such as first and second, and the like may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions.
The above description is only an example of the present application and is not intended to limit the scope of the present application, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, improvement and the like made within the spirit and principle of the present application shall be included in the protection scope of the present application.
Claims (10)
1. A slope influence factor determination method based on uniform design is characterized by comprising the following steps:
determining a test factor; wherein the test factors are derived from a slope structure model to be tested;
determining the value range of each test factor based on the test factors and the slope structure model to be tested;
constructing a uniform design table based on the test factors and the value range of each test factor; the uniform design table comprises test numbers and numerical values of test factors in each test number;
testing the test factors through the uniform design table to obtain a stability coefficient corresponding to each test number;
and determining slope influence sensitive factors in the test factors based on the stability coefficient corresponding to each test number.
2. The method of claim 1, wherein the uniform design table is constructed by;
determining the test times n;
obtain vector set h ═ h (h)1,h2,...hm) (ii) a Wherein h is1,h2,...hmAre all numbers less than n, and the greatest common divisor with n is 1;
generating a jth column of the uniform design table until the uniform design table is obtained;
wherein, the expression in the jth column is: u. ofij＝(i×hj)[mod n](ii) a Wherein, [ mod n]Representing a congruence operation; if i × hjIf n is exceeded, a preset multiple of n is subtracted to make the fall within [1, n ]]To (1); u. ofijGenerating through a recursion algorithm; the uniform design is denoted as Un(nm) (ii) a The number of m corresponding to the determined n and h is obtained by an Euler function E (n);
when n is a prime number, e (n) n-1;
when n is a prime power, i.e. n ═ plWhen the temperature of the water is higher than the set temperature,
p is prime number, l is positive integer;
when n is neither a prime nor a prime power, n is expressed as a square power product of different prime numbers; namely, it is
At this time, the process of the present invention,
p1,p2,L pmare different prime numbers,/1,l2,L lmAre integers.
3. The method of claim 2, wherein the test factors comprise: volume weight, internal friction angle, cohesion, slope angle and seismic coefficient.
4. The method of claim 3, wherein the number of trials n is twice the number of trial factors.
5. The method of claim 4, wherein said testing said test factors through said uniform design table to obtain a stability factor corresponding to each of said test numbers comprises:
testing the test factors through the uniform design table;
obtaining a stability coefficient corresponding to each test number by a simplified Bischu method;
wherein the stability factor is expressed as:
wherein K is the stability coefficient; c. CiThe cohesion corresponding to the ith test; phi is aiThe internal friction angle corresponding to the ith test; wiThe block self weight is the ith test; u. ofiRepresents the pore water pressure; biThe width of the block for the ith test; beta is aiThe slope angle corresponding to the ith test; qHiHorizontal seismic force; kHIs a horizontal seismic coefficient; cZIs a comprehensive influence coefficient; alpha is alphaiIs the seismic acceleration distribution coefficient.
6. The method according to claim 1, wherein the determining the slope influence sensitive factors of the test factors based on the stability coefficient corresponding to each test number comprises:
regression analysis is carried out on the test factors by adopting a backspacing method, and a regression equation is established;
carrying out significance test on the regression equation;
and determining the slope influence sensitive factors in the test factors based on the result of the significance test.
7. The method of claim 1, wherein the performing a significance test on the regression equation comprises:
carrying out significance test on the regression equation to determine the insignificant factors in the test factors;
removing the non-significant factors and reconstructing a new regression equation;
performing significance test on the new regression equation;
correspondingly, the determining of the slope influence sensitive factors in the test factors based on the result of the significance test comprises:
and determining slope influence sensitive factors in the test factors based on the result of the significance test on the new regression equation.
8. The method of claim 1, wherein determining the test factor comprises:
establishing a slope numerical value to construct a test slope structure model; wherein the model is a homogeneous soil slope;
and determining the test factors based on all numerical values of the slope structure model to be tested.
9. A slope influence factor determination device based on uniform design, comprising:
the first determining module is used for determining a test factor; wherein the test factors are derived from a slope structure model to be tested;
the second determination module is used for determining the value range of each test factor based on the test factors and the slope structure to be tested;
the construction module is used for constructing a uniform design table based on the test factors and the value range of each test factor; the uniform design table comprises test numbers and numerical values of test factors in each test number;
the processing module is used for testing the test factors through the uniform design table to obtain a stability coefficient corresponding to each test number;
and the third determining module is used for determining the slope influence sensitive factors in the test factors based on the stability coefficient corresponding to each test number.
10. An electronic device, comprising: a processor and a memory, the processor and the memory connected;
the memory is used for storing programs;
the processor is configured to execute a program stored in the memory to perform the method of any of claims 1-7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110845202.6A CN113609558A (en) | 2021-07-26 | 2021-07-26 | Slope influence factor determination method and device based on uniform design |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110845202.6A CN113609558A (en) | 2021-07-26 | 2021-07-26 | Slope influence factor determination method and device based on uniform design |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113609558A true CN113609558A (en) | 2021-11-05 |
Family
ID=78338351
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110845202.6A Pending CN113609558A (en) | 2021-07-26 | 2021-07-26 | Slope influence factor determination method and device based on uniform design |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113609558A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114357731A (en) * | 2021-12-16 | 2022-04-15 | 中国人民解放军国防科技大学 | Method for acquiring basic model under influence of mixed multi-type factors |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111914457A (en) * | 2020-08-18 | 2020-11-10 | 中国能源建设集团江苏省电力设计院有限公司 | Power transmission line tower footing slope stability judging method and device and storage medium |
| CN113033108A (en) * | 2021-04-19 | 2021-06-25 | 昆明理工大学 | Side slope reliability judgment method based on AdaBoost algorithm |
-
2021
- 2021-07-26 CN CN202110845202.6A patent/CN113609558A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111914457A (en) * | 2020-08-18 | 2020-11-10 | 中国能源建设集团江苏省电力设计院有限公司 | Power transmission line tower footing slope stability judging method and device and storage medium |
| CN113033108A (en) * | 2021-04-19 | 2021-06-25 | 昆明理工大学 | Side slope reliability judgment method based on AdaBoost algorithm |
Non-Patent Citations (4)
| Title |
|---|
| FEI KANG 等: "System probabilistic stability analysis of soil slopes using Gaussian process regression with Latin hypercube sampling", 《COMPUTERS AND GEOTECHNICS》 * |
| 庄楚强 等: "《应用数理统计基础 第4版本》", 31 August 2013 * |
| 苟海军: "基于均匀设计的边坡稳定性分析", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》 * |
| 陈高峰 等: "基于均匀设计的边坡敏感性灰色关联分析", 《地质与勘测》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114357731A (en) * | 2021-12-16 | 2022-04-15 | 中国人民解放军国防科技大学 | Method for acquiring basic model under influence of mixed multi-type factors |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11301029B2 (en) | Method, apparatus, and system for allocating power to graphics processing unit | |
| CN111324471B (en) | Service adjustment method, device, equipment and storage medium | |
| CN105183852A (en) | Database migration method and device | |
| CN113609558A (en) | Slope influence factor determination method and device based on uniform design | |
| Ghosh et al. | System resiliency quantification using non-state-space and state-space analytic models | |
| CN109582649A (en) | A kind of metadata storing method, device, equipment and readable storage medium storing program for executing | |
| CN113792995A (en) | Method, device, equipment and storage medium for determining power resource dominance degree | |
| CN108696397B (en) | Power grid information security assessment method and device based on AHP and big data | |
| CN108599774B (en) | Compression method, system, device and computer readable storage medium | |
| CN116953478B (en) | Ultralow failure rate analysis method and device for integrated circuit and computing equipment | |
| CN111028931B (en) | Medical data processing method and device, electronic equipment and storage medium | |
| KR102280298B1 (en) | Memory management system and method considering application usage patterns analysis | |
| CN117236693A (en) | Storage tank region dynamic risk evaluation method and device | |
| CN116911223A (en) | Ultra-low failure rate upper bound estimation method, device and medium for integrated circuit | |
| CN116340889A (en) | Data acquisition method and system based on big data platform | |
| CN114389953B (en) | Method and system for dynamic expansion and contraction of Kubernetes container based on flow prediction | |
| JP4934660B2 (en) | Communication bandwidth calculation method, apparatus, and traffic management method | |
| CN115128476A (en) | Lithium ion battery pre-lithium amount estimation method, device, equipment and storage medium | |
| CN114385368A (en) | Method and device for flashing pre-written log | |
| CN112187564A (en) | vSAN performance test method, apparatus, computer device and storage medium | |
| US11438408B2 (en) | Transferring applications from overutilized arrays of computer systems to underutilized arrays of computer systems | |
| CN117851729A (en) | Method, device, equipment and medium for calculating soil performance index in basic design | |
| CN107608782A (en) | The method and device of memory compression, operating system, electronic equipment | |
| CN117792800B (en) | Information verification method and system based on Internet of things security evaluation system | |
| Arroyo et al. | Probabilistic model to predict soil-water characteristics on unsaturated soils |
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 |