An engineering-scale soil-rock mixed slope numerical modeling and critical size determination method, system, terminal and storage medium

Abstract

The application belongs to the technical field of numerical simulation research of slope engineering, and discloses an engineering scale soil-rock mixed slope numerical modeling and critical size determination method, system, terminal and storage medium, which comprises the following steps: generating the geometric parameters of the stone blocks by using the discrete element method, and establishing a soil-rock mixed slope finite element model according to the geometric parameters; calculating the safety factor of the soil-rock mixed slope finite element model, and determining the critical size of the soil-rock mixed slope based on the hypothesis testing method to obtain the numerical results required by the calibration soil-rock mixed slope critical size calculation model; calibrating the model parameters of the critical size calculation model based on a preset number of numerical simulation results; obtaining the stone block information of the soil-rock mixed slope in the engineering scale test; substituting the obtained stone block information into the soil-rock mixed slope critical size calculation model, solving and outputting the critical size of the engineering scale soil-rock mixed slope. The application realizes a fast and accurate determination method of the representative unit size of the soil-rock mixed slope under the engineering scale.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of numerical simulation of slope engineering, and in particular to an engineering scale soil-rock mixed slope numerical modeling and critical size determination method, system, terminal and storage medium. BACKGROUND

[0002] The soil-rock mixed slope refers to a slope body composed of gravel and soil particles in a certain proportion, which has significant heterogeneity. Due to the large difference in strength between the soil and the stone and the strong randomness of the spatial distribution of the stone, the soil-rock mixed slope is prone to instability disasters under adverse conditions such as weathering erosion, rainfall infiltration and engineering disturbance, causing casualties and economic losses. And building a slope model that can truly reflect the mechanical response of engineering scale under the condition of containing stone is the primary prerequisite for reliable evaluation of the safety and stability of the soil-rock mixed slope.

[0003] Because the soil-rock mixed slope is formed by weathering and deposition process, the content, gradation and shape of the stone in different regions differ significantly, showing strong spatial variability. The extreme heterogeneity of the stone makes the mechanical response reflected by the soil-rock mixed slope calculation model of different sizes fluctuate, and the failure slip surface of the slope when it fails has development modes such as "stone around" and "stone wrapping". However, the traditional modeling method is difficult to establish an engineering scale high-stone-content soil-rock mixed slope model, and the spatial distribution characteristics of the stone in the engineering scale slope are not fully quantitatively considered when determining the model size, therefore, a large amount of trial and comparison tests are needed to verify the rationality of the model size, which consumes a lot of test resources or computing power, seriously restricting the efficiency and accuracy of the high-stone-content slope stability evaluation.

[0004] The prior art also has the problem of low accuracy of numerical modeling and critical size determination of engineering scale soil-rock mixed slope, so the prior art still needs to be improved. SUMMARY

[0005] The technical problem to be solved by the present application is that, in view of the defects of the prior art, the present application provides an engineering scale soil-rock mixed slope numerical modeling and critical size determination method, system, terminal and storage medium, to solve the problem of low efficiency and poor accuracy of engineering scale slope stability evaluation due to the difficulty of traditional modeling method in constructing high-stone-content slope and the failure to quantitatively consider the influence of stone characteristics on model critical size.

[0006] The technical solution adopted by the present application to solve the technical problem is as follows:

[0007] In a first aspect, the present application provides an engineering scale soil-rock mixed slope numerical modeling and critical size determination method, comprising:

[0008] Based on the preset slope size and the stone information, a discrete element method is used to generate geometric parameters of the stones, and a soil-rock mixed slope finite element model is established according to the geometric parameters;

[0009] A safety factor of the soil-rock mixed slope finite element model is calculated by using a strength reduction method, and a critical size of the soil-rock mixed slope is determined based on a hypothesis testing method, so as to obtain a numerical result required by a calibrated soil-rock mixed slope critical size calculation model;

[0010] Based on a preset number of numerical simulation results, model parameters of the soil-rock mixed slope critical size calculation model are calibrated;

[0011] Stone information of a soil-rock mixed slope in an engineering scale test is obtained;

[0012] The obtained stone information is substituted into the soil-rock mixed slope critical size calculation model, and an engineering scale soil-rock mixed slope critical size is solved and output.

[0013] In an implementation manner, based on the preset slope size and the stone information, a discrete element method is used to generate geometric parameters of the stones, and a soil-rock mixed slope finite element model is established according to the geometric parameters, including:

[0014] Based on the preset slope size and the stone information, a stone boundary of the slope size shape is set in a discrete element software;

[0015] Initial stones with random coordinates are generated in the stone boundary, and initial sizes and velocity vectors of the initial stones are set;

[0016] The initial stones are expanded based on the initial sizes, when the stones collide with each other or with a preset boundary, the stones will start to move based on Newton's second law, so as to realize random distribution of the stones in the slope boundary, and geometric parameters of the random distribution of the initial stones in the stone boundary are derived;

[0017] A finite element analysis software is called, the geometric parameters are imported into the finite element analysis software, and a model instance is generated, so as to establish the soil-rock mixed slope finite element model.

[0018] In an implementation manner, the safety factor of the soil-rock mixed slope finite element model is calculated by using the strength reduction method, including:

[0019] Based on the strength reduction method, a soil body strength parameter is reduced by using the following formula:

[0020] ;

[0021] ;

[0022] wherein,c and φ are respectively the soil cohesion and internal friction angle; c r and φ r are respectively the reduced soil strength parameters, F r is the material performance reduction factor;

[0023] According to the reduced soil strength parameters, numerical calculation non-convergence is selected as the judgment basis for slope instability, and the reduction factor at the time of instability is taken as the safety factor of the soil-rock mixed slope finite element model.

[0024] In an implementation, the assumption testing method is used to determine the critical size of the soil-rock mixed slope, and numerical results required for a critical size calculation model of the soil-rock mixed slope are obtained, including:

[0025] The soil-rock mixed slope model is calculated 10 times for each group of representative unit sizes based on the assumption testing method, and the critical size calculation model of the soil-rock mixed slope is obtained:

[0026] ;

[0027] ;

[0028] ;

[0029] wherein, n is the sample size; X n is the safety factor of each calculation of the slope model; and S are respectively the average value and standard deviation of the safety factor of each group of soil-rock mixed slope models; ε is the relative error in the assumption testing method t- ; t α (n-1) represents the degree of freedom n -1; t is the theoretical error critical value under a given confidence interval α ;

[0030] The size ratio of the soil-rock mixed slope model when the relative error ε is less than 1% is taken as the critical size of the soil-rock mixed slope.

[0031] In an implementation, the model parameters of the critical size calculation model of the soil-rock mixed slope are calibrated based on a preset number of numerical simulation results, including:

[0032] According to the critical size calculation model of the soil-rock mixed slope, the critical size is determined H / d a According to the linear relationship between the rock content and the non-uniformity coefficient of gradation C u

[0033] ;

[0034] wherein, k 、 b 、 c are model parameters;

[0035] Based on the linear relationship, the model parameters of the critical size calculation model of the soil-rock mixed slope are calibrated by using a preset number of rock contents and non-uniformity coefficients of gradation C u .

[0036] In an implementation manner, the rock information of the soil-rock mixed slope in the engineering scale test is obtained, including:

[0037] The size of the rock in the soil-rock mixed slope in the engineering scale test is obtained d a , the rock content and the non-uniformity coefficient C u , to obtain the rock information.

[0038] In an implementation manner, the rock information obtained is substituted into the critical size calculation model of the soil-rock mixed slope, and the critical size of the engineering scale soil-rock mixed slope is solved and output, including:

[0039] The model parameters k , the model parameters b , the model parameters c , the rock content and the non-uniformity coefficient C u are substituted into the critical size calculation model of the soil-rock mixed slope, to obtain the critical size of the engineering scale soil-rock mixed slope;

[0040] The critical size of the engineering scale soil-rock mixed slope is output.

[0041] In a second aspect, the present application provides a numerical modeling and critical size determination system for an engineering scale soil-rock mixed slope, including:

[0042] ​The finite element model establishing module is configured to generate geometric parameters of the stone blocks by using a discrete element method based on preset sizes of the slope and the stone block information, and to establish a finite element model of the soil and stone mixed slope according to the geometric parameters.

[0043] The critical size calculation model generating module is configured to calculate a safety factor of the finite element model of the soil and stone mixed slope by using a strength reduction method, and to determine a critical size of the soil and stone mixed slope based on a hypothesis testing method to obtain numerical results required for calibrating the critical size calculation model of the soil and stone mixed slope.

[0044] The model parameter calibrating module is configured to calibrate model parameters of the critical size calculation model of the soil and stone mixed slope based on a preset number of numerical simulation results.

[0045] The stone block information acquiring module is configured to acquire stone block information of the soil and stone mixed slope in an engineering scale test.

[0046] The critical size calculation and output module is configured to substitute the acquired stone block information into the critical size calculation model of the soil and stone mixed slope, to solve and output the critical size of the soil and stone mixed slope in the engineering scale.

[0047] In a third aspect, the present application provides a terminal, comprising a processor and a memory, wherein the memory stores an engineering scale soil and stone mixed slope numerical modeling and critical size determination program, and the engineering scale soil and stone mixed slope numerical modeling and critical size determination program is configured to implement operations of the engineering scale soil and stone mixed slope numerical modeling and critical size determination method according to the first aspect when executed by the processor.

[0048] In a fourth aspect, the present application further provides a computer readable storage medium, wherein the computer readable storage medium stores an engineering scale soil and stone mixed slope numerical modeling and critical size determination program, and the engineering scale soil and stone mixed slope numerical modeling and critical size determination program is configured to implement operations of the engineering scale soil and stone mixed slope numerical modeling and critical size determination method according to the first aspect when executed by a processor.

[0049] The technical solution of the present application has the following effects:

[0050] The application is based on presetting the size of the slope and the information of the stone blocks, using the discrete element method to generate the geometric parameters of the stone blocks, and establishing a finite element model of the soil and stone mixed slope according to the geometric parameters; using the strength reduction method to calculate the safety factor of the finite element model of the soil and stone mixed slope, and determining the critical size of the soil and stone mixed slope based on the hypothesis testing method to obtain the numerical results required by the calibration of the critical size calculation model of the soil and stone mixed slope; based on a preset number of numerical simulation results, the model parameters of the critical size calculation model of the soil and stone mixed slope are calibrated; the information of the stone blocks in the soil and stone mixed slope in the engineering scale test is obtained; the obtained stone block information is substituted into the critical size calculation model of the soil and stone mixed slope, and the critical size of the soil and stone mixed slope in the engineering scale is solved and output; the application provides a simple and efficient numerical modeling and critical size determination method for the soil and stone mixed slope in the engineering scale, which can quickly determine the representative unit size of the soil and stone mixed slope in the engineering scale according to only 4 groups of numerical simulation results, provides a theoretical basis for slope disaster analysis, and ensures the safety of engineering infrastructure. BRIEF DESCRIPTION OF DRAWINGS

[0051] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the drawings needed to be used in the embodiments or the prior art description will be briefly introduced as follows. Obviously, the drawings in the following description only constitute some embodiments of the present application, and other drawings can also be obtained according to the structures shown in the drawings without creative labor for those skilled in the art.

[0052] Figure 1 It is a flow chart of the numerical modeling and critical size determination method for the soil and stone mixed slope in the engineering scale in the present application.

[0053] Figure 2 It is a modeling flow chart of the high-content soil and stone mixed slope in the present application.

[0054] Figure 3 It is a comparison chart of the critical size prediction value and the simulation value of the soil and stone mixed slope in the present application. t- Test relative error ε Schematic diagram of the change law with the increase of the slope model size ratio.

[0055] Figure 4 It is a comparison chart of the critical size prediction value and the simulation value of the soil and stone mixed slope in the present application.

[0056] Figure 5 It is a functional principle diagram of the terminal in one implementation manner of the present application.

[0057] The purpose implementation, functional characteristics and advantages of the present application will be further described with reference to the embodiments and the drawings. DETAILED DESCRIPTION

[0058] In order to make the objects, technical solutions and advantages of the present application clearer and more apparent, the present application will be further described in detail below with reference to the drawings and examples. It should be understood that the specific examples described herein are only used to explain the present application and not to limit the present application.

[0059] Exemplary method

[0060] At present, the soil-rock mixed slope is mostly formed by weathering and deposition process, and the content, gradation and shape of the stone block in different regions are significantly different, showing strong spatial variability. The extreme heterogeneity of the stone block makes the mechanical response reflected by the calculation model of the soil-rock mixed slope of different sizes fluctuate, and the failure sliding surface of the slope when the slope is unstable has development modes such as "stone wrapping" and "stone wrapping". However, the traditional modeling method is difficult to establish a soil-rock mixed slope model of high stone content at an engineering scale, and the spatial distribution characteristics of the stone block in the slope at an engineering scale are not fully quantitatively considered when determining the model size, therefore, a large amount of trial and comparison tests are required to verify the rationality of the model size, which consumes a large amount of test resources or computing power, and seriously restricts the efficiency and accuracy of the stability evaluation of the high stone content slope.

[0061] In view of the above technical problems, the present application provides an engineering scale soil-rock mixed slope numerical modeling and critical size determination method, comprising: generating the geometric parameters of the stone block by using the discrete element method, and establishing a soil-rock mixed slope finite element model according to the geometric parameters; calculating the safety factor of the soil-rock mixed slope finite element model by using the strength reduction method, and determining the critical size of the soil-rock mixed slope based on the hypothesis testing method to obtain the numerical results required for calibrating the critical size calculation model of the soil-rock mixed slope; calibrating the model parameters of the critical size calculation model based on a preset number of numerical simulation results; obtaining the stone block information of the soil-rock mixed slope in the engineering scale test; substituting the obtained stone block information into the critical size calculation model of the soil-rock mixed slope, solving and outputting the critical size of the engineering scale soil-rock mixed slope. The present application provides a simple and efficient large-scale soil-rock mixed slope modeling and critical size determination method, which only needs to be based on 4 groups of numerical simulation results to quickly determine the representative unit size of the soil-rock mixed slope at an engineering scale, provides a theoretical basis for slope disaster analysis, and ensures the safety of slope engineering infrastructure.

[0062] As shown in Figure 1 The present application provides an engineering scale soil-rock mixed slope numerical modeling and critical size determination method, comprising the following steps:

[0063] Step S100, based on the preset slope size and stone block information, the geometric parameters of the stone block are generated by using the discrete element method, and a soil-rock mixed slope finite element model is established according to the geometric parameters.

[0064] In this embodiment, based on the preset slope dimensions and stone information, the geometric parameters of the stones are generated using the discrete element method. The geometric parameters of the stones are then imported into Abaqus (a finite element analysis software) using Python (a general programming language) to establish a finite element model of the soil-rock mixed slope.

[0065] Specifically, in one implementation of this embodiment, step S100 includes the following steps:

[0066] Step S101: Based on the preset slope dimensions and stone information, set the stone boundary of the slope size and shape in the discrete element software;

[0067] Step S102: Generate an initial stone block with random coordinates within the boundary of the stone block, and set an initial size for the initial stone block;

[0068] Step S103: Based on the initial size, the initial stone is expanded. When the stone collides with each other or with the preset boundary, the stone will start to move based on Newton's second law, realizing the random distribution of the stone within the slope boundary, and deriving the geometric parameters of the random distribution of the initial stone within the stone boundary.

[0069] Step S104: Call the finite element analysis software, import the geometric parameters into the finite element analysis software and generate a model instance to establish the finite element model of the soil-rock mixed slope.

[0070] In this embodiment, establishing the finite element model of the soil-rock mixed slope mainly includes:

[0071] First, in the discrete element method (DEM) software, the boundary of the stone blocks with the desired slope size and shape is set according to the preset slope dimensions. Then, initial stone blocks with random coordinates are generated within this boundary, and they are assigned small initial sizes and velocity vectors. Next, the motion of these initial stone blocks is simulated. During the simulation, the size of the stone blocks continuously expands until they reach a specified stone content. and gradation unevenness coefficient C u When stones collide with each other or with the boundary, their direction of movement changes accordingly, ultimately resulting in a random distribution of stones within the boundary. Finally, the geometric parameters of the stones are recorded and output, including size, coordinates, and rotation information.

[0072] After obtaining the geometric parameters of the stones, Python code is used to call the Abaqus Part module (a core functional module for finite element analysis preprocessing, mainly used to create and manage geometric parts). Relying on the Scripting Interface, the geometric parameters of the stones are imported into model instances, creating batch finite element models of soil-rock mixed slopes containing the randomly distributed stones. In the soil-rock mixed slope finite element model, H / d a Defined as the size ratio of a soil-rock mixed slope. d a The equivalent diameter of the rock or the polydisperse ligand d 50 .

[0073] As an example, the modeling process for a high-content soil-rock mixed slope model at the engineering scale in this embodiment is as follows: Figure 2 As shown, the size, coordinates, and rotation information of the stone were generated using the discrete element method software Yade to obtain the geometric parameters of the stone. Then, a Python program was written to call the Part module of Abaqus to import the geometric parameters of the stone into a model instance, batch-constructing models with stone contents of 10% and 55% and different coefficients of uniformity. C u The simulation schemes for the soil-rock mixed slope models 1 and 4 (both with a rock content of 40%) are shown in Table 1.

[0074] Table 1. Simulation scheme for soil-rock mixed slope model:

[0075]

[0076] The boundary conditions of the model are as follows Figure 2 As shown in (c), displacement in the x and y directions is restricted at the bottom of the slope, and displacement in the x direction is restricted on both sides. The soil and rocks are meshed using six-node triangular elements with plane strain, and the mesh is locally refined at the contact surface between the rocks and the soil and on the slope surface. The material parameters of the soil-rock mixture are shown in Table 2.

[0077] Table 2. Material parameters of soil-rock mixture:

[0078]

[0079] In this embodiment, a finite element model of a soil-rock mixed slope is established in the manner described above. The safety factor of the model can be calculated using the strength reduction method, and the critical dimensions of the soil-rock mixed slope can be determined by the hypothesis testing method. This results in a critical dimension calculation model for the soil-rock mixed slope that can quickly calculate the representative unit dimensions of the soil-rock mixed slope at the engineering scale.

[0080] like Figure 1 As shown in the figure, this invention provides a method for numerical modeling and critical dimension determination of engineering-scale soil-rock mixed slopes, including the following steps:

[0081] Step S200: The safety factor of the finite element model of the soil-rock mixed slope is calculated using the strength reduction method, and the critical dimensions of the soil-rock mixed slope are determined based on the hypothesis testing method, so as to obtain the numerical results required for calibrating the critical dimension calculation model of the soil-rock mixed slope.

[0082] Specifically, in one implementation of this embodiment, step S200 includes the following steps:

[0083] Step S201: Based on the strength reduction method, reduce the soil strength parameters using the following formula;

[0084] Step S202: Based on the reduced soil strength parameters, the non-convergence of numerical calculation is selected as the criterion for judging slope instability, and the reduction coefficient at the time of instability is used as the safety factor of the finite element model of the soil-rock mixed slope.

[0085] In this embodiment, based on the batch-established finite element models of soil-rock mixed slopes in Table 1 above, the slope safety factor is calculated using the strength reduction method:

[0086] ;

[0087] ;

[0088] in, c and φ These are soil cohesion and internal friction angle, respectively. c r and φ r These are the reduced soil strength parameters. F r This is the material performance reduction factor.

[0089] Then, the non-convergence of numerical calculations was used as the criterion for judging slope instability, and the reduction factor at the time of instability was used as the safety factor of the slope.

[0090] Step S203: Based on the hypothesis testing method, determine the soil-rock mixed slope model by calculating the distribution of stone positions 10 times for each group of representative unit sizes, and obtain the critical size calculation model of the soil-rock mixed slope.

[0091] Step S204, relative error ε The size ratio of the soil-rock mixed slope model when it is less than 1% is used as the critical size of the soil-rock mixed slope.

[0092] In this embodiment, based on the safety factor of the finite element model of the soil-rock mixed slope, the slope model corresponding to the size of each group of representative elements, calculated 10 times for the distribution of rock positions, is determined, thus obtaining the critical size calculation model of the soil-rock mixed slope; the calculation method of the critical size calculation model of the soil-rock mixed slope is shown in the following formula:

[0093] ;

[0094] ;

[0095] ;

[0096] in, n For sample size; X n The safety factor is calculated for each slope model; and S These represent the average and standard deviation of the safety factor for each group of soil-rock mixed slope models; ε Hypothesis testing method t- Relative error in testing; t α (n-1) represents the number of degrees of freedom. n -1 t Distributed within a given confidence interval α The theoretical error threshold can be determined by directly looking up... t The distribution table was obtained.

[0097] In this embodiment, in relative error ε When the relative error is less than 1%, the size ratio of the soil-rock mixed slope model is considered to be at the critical size. ε The variation pattern with increasing size ratio is as follows: Figure 3 As shown, Figure 3 (a) shows the critical dimensions of a soil-rock mixed slope with a 10% rock content. Figure 3 (b) shows the critical dimensions of a soil-rock mixed slope with 55% rock content. Figure 3 (c) is C u =1、 C u A schematic diagram of the critical dimensions of a soil-rock mixed slope with a rock content of 4; where the rock content is... Soil-rock mixed slopes of 10% and 55% respectively ( C u 1) The critical dimensions are 4 and 12 respectively. C u The critical dimensions for soil-rock mixed slopes (with a rock content of 40%) with values ​​of 1 and 4 are 10 and 11, respectively.

[0098] like Figure 1 As shown in the figure, this invention provides a method for numerical modeling and critical dimension determination of engineering-scale soil-rock mixed slopes, including the following steps:

[0099] Step S300: Based on a preset number of numerical simulation results, calibrate the model parameters of the critical dimension calculation model for the soil-rock mixed slope.

[0100] Specifically, in one implementation of this embodiment, step S300 includes the following steps:

[0101] Step S301: Determine the critical dimensions based on the critical dimension calculation model for the soil-rock mixed slope. H / d a With stone content and gradation unevenness coefficient C u A linear relationship;

[0102] Step S302: Based on the linear relationship, use a preset amount of stone. and gradation unevenness coefficient C u The model parameters of the calculation model for the critical dimensions of the soil-rock mixed slope were calibrated.

[0103] In this embodiment, the critical dimensions are determined based on the aforementioned calculation model for the critical dimensions of soil-rock mixed slopes. H / d a With stone content and gradation unevenness coefficient C u The linear relationship is as follows:

[0104] ;

[0105] in, k , b , c These are all model parameters; only two types of stone content and two types of... C u The critical dimensions of the soil-rock mixed slope model can be calibrated; the model parameters are obtained through fitting. k , b and c The values ​​are 1.60, 19.67, and 0.25, respectively.

[0106] like Figure 1 As shown in the figure, this invention provides a method for numerical modeling and critical dimension determination of engineering-scale soil-rock mixed slopes, including the following steps:

[0107] Step S400: Obtain information on the rocks in the soil-rock mixed slope during the engineering scale test.

[0108] Specifically, in one implementation of this embodiment, step S400 includes the following steps:

[0109] Step S401: Obtain the dimensions of the stones in the soil-rock mixed slope during the engineering scale test. d a Stone content and coefficient of uniformity C u The information about the stone was obtained.

[0110] In this embodiment, the dimensions of the stones in the soil-rock mixed slope are obtained during engineering-scale tests. d a Stone content and coefficient of uniformity C u The information on the rocks in the soil-rock mixed slope is obtained.

[0111] As an example, in a sample of a soil-rock mixed slope, the rock content of the slope... The coefficients of non-uniformity are 40% and 50%, respectively. C u They are 2 and 3 respectively.

[0112] like Figure 1 As shown in the figure, this invention provides a method for numerical modeling and critical dimension determination of engineering-scale soil-rock mixed slopes, including the following steps:

[0113] Step S500: Substitute the obtained stone information into the critical dimension calculation model of the soil-rock mixed slope, solve and output the critical dimension of the soil-rock mixed slope in engineering scale.

[0114] Specifically, in one implementation of this embodiment, step S500 includes the following steps:

[0115] Step S501, set the model parameters k Model parameters b Model parameters c The stone content and the non-uniformity coefficient C u Substituting the values ​​into the critical dimension calculation model for the soil-rock mixed slope, the critical dimensions of the soil-rock mixed slope at the engineering scale are obtained;

[0116] Step S502: Output the critical dimensions of the engineering-scale soil-rock mixed slope.

[0117] As an example, in this embodiment, the model parameters are... k Model parametersb Model parameters c The stone content and the non-uniformity coefficient C u Substituting into the critical dimension calculation model for soil-rock mixed slopes, calculate the representative unit dimensions for the following four groups of soil-rock mixed slopes:

[0118] 1) Coefficient of non-uniformity C u The value is 2, and the stone content is 40%.

[0119] 2) Coefficient of non-uniformity C u It has a content of 3% and a stone content of 40%.

[0120] 3) Coefficient of non-uniformity C u 2. Stone content is 50%;

[0121] 4) Coefficient of non-uniformity C u It has a value of 3 and a stone content of 50%.

[0122] After substitution and calculation, the representative unit sizes of the above-mentioned soil-rock mixed slopes are 10.0, 10.2, 11.9, and 12.2, respectively.

[0123] To verify the effectiveness of the engineering-scale soil-rock mixed slope critical dimension calculation model proposed in this embodiment, more than 500 additional sets of soil-rock mixed slope critical dimension simulation tests were carried out, as shown in Table 3.

[0124] Table 3. Simulation scheme for verifying the soil-rock mixed slope model:

[0125]

[0126] The predicted and simulated values ​​of the critical dimension calculation model for engineering-scale soil-rock mixed slopes proposed in this embodiment are as follows: Figure 4 As shown, Figure 4 The results show that the critical size calculation model proposed in this embodiment has good prediction accuracy.

[0127] This embodiment achieves the following technical effects through the above technical solution:

[0128] This embodiment proposes a simple and efficient numerical modeling and critical size determination method for soil-rock mixed slopes at the engineering scale. Based on only four sets of numerical simulation results, the representative unit size of soil-rock mixed slopes at the engineering scale can be quickly determined, providing a theoretical basis for slope disaster analysis and ensuring the safety of engineering infrastructure.

[0129] Exemplary device

[0130] Based on the above embodiments, the present invention also provides a numerical modeling and critical dimension determination system for engineering-scale soil-rock mixed slopes, comprising:

[0131] The finite element model building module is used to generate the geometric parameters of the stones based on the preset slope dimensions and stone information using the discrete element method, and to build a finite element model of the soil-rock mixed slope based on the geometric parameters.

[0132] The critical dimension calculation model generation module is used to calculate the safety factor of the finite element model of the soil-rock mixed slope using the strength reduction method, and to determine the critical dimension of the soil-rock mixed slope based on the hypothesis testing method, so as to obtain the numerical results required for calibrating the critical dimension calculation model of the soil-rock mixed slope.

[0133] The model parameter calibration module is used to calibrate the model parameters of the critical dimension calculation model of the soil-rock mixed slope based on a preset number of numerical simulation results;

[0134] The stone information acquisition module is used to acquire stone information of soil-rock mixed slopes in engineering scale tests;

[0135] The critical dimension calculation and output module is used to substitute the acquired stone information into the critical dimension calculation model of the soil-rock mixed slope, solve and output the critical dimensions of the soil-rock mixed slope in engineering scale.

[0136] This embodiment achieves the following technical effects through the above technical solution:

[0137] The numerical modeling and critical size determination system for soil-rock mixed slopes at the engineering scale provided in this embodiment can quickly determine the representative unit size of soil-rock mixed slopes at the engineering scale based on only 4 sets of numerical simulation results, providing a theoretical basis for slope disaster analysis and ensuring the safety of engineering infrastructure.

[0138] Based on the above embodiments, the present invention also provides a terminal, the principle block diagram of which can be as follows: Figure 5 As shown.

[0139] The terminal includes: a processor, a memory, an interface, a display screen, and a communication module connected via a system bus; wherein, the processor of the terminal provides computing and control capabilities; the memory of the terminal includes a computer-readable storage medium and internal memory; the computer-readable storage medium stores an operating system and computer programs; the internal memory provides an environment for the operation of the operating system and computer programs in the computer-readable storage medium; the interface is used to connect to external devices; the display screen is used to display relevant information; and the communication module is used to communicate with a cloud server or other devices.

[0140] When executed by a processor, this computer program is used to implement methods for numerical modeling and critical dimension determination of engineering-scale soil-rock mixed slopes.

[0141] It will be understood by those skilled in the art that Figure 5 The schematic diagram shown is merely a partial structural diagram related to the present invention and does not constitute a limitation on the terminal to which the present invention is applied. A specific terminal may include more or fewer components than those shown in the figure, or combine certain components, or have different component arrangements.

[0142] In one embodiment, a terminal is provided, comprising: a processor and a memory, the memory storing a program for numerical modeling and critical dimension determination of engineering-scale soil-rock mixed slopes, the program being executed by the processor to implement the above-described method for numerical modeling and critical dimension determination of engineering-scale soil-rock mixed slopes.

[0143] In one embodiment, a computer-readable storage medium is provided, wherein the computer-readable storage medium stores a program for numerical modeling and critical dimension determination of engineering-scale soil-rock mixed slopes, the program being executed by a processor to implement the above-described method for numerical modeling and critical dimension determination of engineering-scale soil-rock mixed slopes.

[0144] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile storage medium, and when executed, it can include the processes of the embodiments of the methods described above. Any references to memory, storage, database, or other media used in the embodiments provided by this invention can include both non-volatile and volatile memory.

[0145] In summary, this invention provides a method, system, terminal, and storage medium for numerical modeling and critical dimension determination of soil-rock mixed slopes at the engineering scale. The method includes: generating geometric parameters of the rocks using the discrete element method and establishing a finite element model of the soil-rock mixed slope based on these parameters; calculating the safety factor of the finite element model of the soil-rock mixed slope using the strength reduction method and determining the critical dimensions of the soil-rock mixed slope based on hypothesis testing, obtaining the numerical results required for calibrating the critical dimension calculation model of the soil-rock mixed slope; calibrating the model parameters of the critical dimension calculation model based on a preset number of numerical simulation results; acquiring rock information of the soil-rock mixed slope in engineering-scale tests; substituting the acquired rock information into the critical dimension calculation model of the soil-rock mixed slope, solving for and outputting the critical dimensions of the soil-rock mixed slope at the engineering scale. This invention achieves a rapid and accurate method for determining the representative unit dimensions of soil-rock mixed slopes at the engineering scale.

[0146] It should be understood that the application of the present invention is not limited to the examples above. Those skilled in the art can make improvements or modifications based on the above description, and all such improvements and modifications should fall within the protection scope of the appended claims.

Claims

1. A method for numerical modeling and critical dimension determination of soil-rock mixed slopes at engineering scale, characterized in that, include: Based on the preset slope dimensions and stone information, the discrete element method is used to generate the geometric parameters of the stones, and a finite element model of the soil-rock mixed slope is established based on the geometric parameters. The safety factor of the finite element model of the soil-rock mixed slope was calculated using the strength reduction method, and the critical dimensions of the soil-rock mixed slope were determined based on the hypothesis testing method, thus obtaining the numerical results required for calibrating the critical dimension calculation model of the soil-rock mixed slope. The model parameters of the critical dimension calculation model for the soil-rock mixed slope are calibrated based on a preset number of numerical simulation results; Obtain information on the rocks in a soil-rock mixed slope during engineering-scale testing; Substitute the acquired stone information into the critical dimension calculation model of the soil-rock mixed slope to solve and output the critical dimension of the soil-rock mixed slope in engineering scale. Based on preset slope dimensions and stone information, the discrete element method is used to generate the geometric parameters of the stones, and a finite element model of the soil-rock mixed slope is established based on the geometric parameters, including: Based on the preset slope dimensions and stone information, the stone boundaries of the slope dimensions and shapes are set in the discrete element software. An initial stone with random coordinates is generated within the boundary of the stone block, and an initial size is set for the initial stone block; Based on the initial size, the initial stone is expanded. When the stone collides with each other or with the preset boundary, the stone will start to move based on Newton's second law, realizing the random distribution of the stone within the slope boundary, and deriving the geometric parameters of the random distribution of the initial stone within the stone boundary. The geometric parameters are imported into the finite element analysis software and a model instance is generated to establish the finite element model of the soil-rock mixed slope. The method of determining the critical dimensions of soil-rock mixed slopes based on hypothesis testing yields the numerical results required for calibrating the critical dimension calculation model of soil-rock mixed slopes, including: Based on the hypothesis testing method, the size of each representative unit is determined, and the soil-rock mixed slope model with 10 rock block location distributions is calculated to obtain the critical size calculation model of the soil-rock mixed slope: ； ； ； in, n For sample size; X n The safety factor is calculated for each slope model; and S These represent the average and standard deviation of the safety factor for each group of soil-rock mixed slope models; ε Hypothesis testing method t- Relative error in testing; t α (n-1) represents the number of degrees of freedom. n -1 t Distributed within a given confidence interval α The theoretical critical error value; relative error ε The size ratio of the soil-rock mixed slope model when it is less than 1% is taken as the critical size of the soil-rock mixed slope; The model parameters of the critical dimension calculation model for the soil-rock mixed slope are calibrated based on a preset number of numerical simulation results, including: Based on the aforementioned critical dimension calculation model for soil-rock mixed slopes, the critical dimensions are determined. H / d a With stone content and gradation unevenness coefficient C u Linear relationship: ； in, k , b , c All are model parameters; Based on the aforementioned linear relationship, using a preset quantity of stones and gradation unevenness coefficient C u The model parameters of the calculation model for the critical dimensions of the soil-rock mixed slope were calibrated.

2. The method for numerical modeling and critical dimension determination of engineering-scale soil-rock mixed slopes according to claim 1, characterized in that, The calculation of the safety factor of the finite element model of the soil-rock mixed slope using the strength reduction method includes: Based on the strength reduction method, the soil strength parameters are reduced using the following formula: ； ； in, c and φ These are soil cohesion and internal friction angle, respectively. c r and φ r These are the reduced soil strength parameters. F r This is the material property reduction factor; Based on the reduced soil strength parameters, the non-convergence of numerical calculations was selected as the criterion for judging slope instability, and the reduction factor at the time of instability was used as the safety factor of the finite element model of the soil-rock mixed slope.

3. The method for numerical modeling and critical dimension determination of engineering-scale soil-rock mixed slopes according to claim 1, characterized in that, The acquisition of rock information from soil-rock mixed slopes in engineering-scale tests includes: Obtain the dimensions of the stones in the soil-rock mixed slope during the engineering scale test. d a Stone content and coefficient of uniformity C u The information about the stone was obtained.

4. The method for numerical modeling and critical dimension determination of engineering-scale soil-rock mixed slopes according to claim 3, characterized in that, The process of substituting the acquired stone information into the critical dimension calculation model for the soil-rock mixed slope, solving for and outputting the critical dimensions of the soil-rock mixed slope at engineering scale, includes: Model parameters k Model parameters b Model parameters c The stone content and the non-uniformity coefficient C u Substituting the values ​​into the critical dimension calculation model for the soil-rock mixed slope, the critical dimensions of the soil-rock mixed slope at the engineering scale are obtained; Output the critical dimensions of the soil-rock mixed slope at the engineering scale.

5. A numerical modeling and critical dimension determination system for engineering-scale soil-rock mixed slopes, used to implement the numerical modeling and critical dimension determination method for engineering-scale soil-rock mixed slopes as described in any one of claims 1-4, characterized in that, include: The finite element model building module is used to generate the geometric parameters of the stones based on the preset slope dimensions and stone information using the discrete element method, and to build a finite element model of the soil-rock mixed slope based on the geometric parameters. The critical dimension calculation model generation module is used to calculate the safety factor of the finite element model of the soil-rock mixed slope using the strength reduction method, and to determine the critical dimension of the soil-rock mixed slope based on the hypothesis testing method, so as to obtain the numerical results required for calibrating the critical dimension calculation model of the soil-rock mixed slope. The model parameter calibration module is used to calibrate the model parameters of the critical dimension calculation model of the soil-rock mixed slope based on a preset number of numerical simulation results; The stone information acquisition module is used to acquire stone information of soil-rock mixed slopes in engineering scale tests; The critical dimension calculation and output module is used to substitute the acquired stone information into the critical dimension calculation model of the soil-rock mixed slope, solve and output the critical dimensions of the soil-rock mixed slope in engineering scale.

6. A terminal, characterized in that, include: The processor and memory, wherein the memory stores a program for numerical modeling and critical dimension determination of engineering-scale soil-rock mixed slopes, wherein the program for numerical modeling and critical dimension determination of engineering-scale soil-rock mixed slopes, when executed by the processor, is used to implement the operation of the method for numerical modeling and critical dimension determination of engineering-scale soil-rock mixed slopes as described in any one of claims 1-4.

7. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a program for numerical modeling and critical dimension determination of engineering-scale soil-rock mixed slopes. When executed by a processor, the program is used to implement the operation of the method for numerical modeling and critical dimension determination of engineering-scale soil-rock mixed slopes as described in any one of claims 1-4.

Images