Static cone penetration test parameter inversion method, system, equipment and medium

Abstract

The invention discloses a static cone penetration test parameter inversion method, system, equipment and medium, and the method comprises the steps: firstly obtaining basic parameters of target region soil according to set variables, and carrying out the numerical simulation of a static cone penetration test based on a material point method according to the basic parameters of the target region soil to obtain a static cone penetration test simulation result; carrying out regression analysis by adopting the static sounding test simulation result to obtain a regression equation of the static sounding test simulation result, and carrying out regression analysis by adopting multiple groups of static sounding test actual measurement results to obtain a regression equation of the static sounding test actual measurement results; and finally, performing inversion according to the regression equation of the simulation result of the static penetration test and the regression equation of the actual measurement result of the static penetration test to obtain static penetration test inversion parameters of the target area soil. By means of the method, unnecessary repeated static sounding tests can be avoided, resource waste caused by the repeated tests is reduced, and therefore the efficiency of engineering geological exploration is improved.

Description

Technical Field

[0001] The present invention relates to the field of novel low-carbon geotechnical engineering survey technology, and in particular to a static penetration test parameter inversion method, system, equipment and medium. Background Art

[0002] The static penetration test uses static pressure to uniformly press a conical probe into the soil at a certain rate to measure its penetration resistance, and divide the soil layers according to the resistance to determine the engineering properties of the soil. The static penetration test is widely used in engineering geological surveys to determine the spatial distribution of various types of soil and their engineering characteristics, such as dividing soil layers and identifying soil types, calculating the engineering property indicators of soil layers, and determining the bearing layer at the pile end and the bearing capacity of a single pile.

[0003] Among them, the method for obtaining the engineering property index of the soil layer (such as penetration resistance) is generally for engineers to study theoretical models for static penetration tests, and obtain values ​​such as cone tip resistance and specific penetration resistance through static penetration tests, and then establish empirical formulas applicable to different regions or different soil properties by analyzing the relationship between these values ​​and the friction angle and cohesion of the soil. Although these formulas are still applicable in actual engineering, they are all summarized by engineers based on a lot of experience and can only be applied under specific conditions. In actual static penetration tests, repeated tests may be necessary due to factors such as site selection and terrain, resulting in unnecessary resource loss. Summary of the invention

[0004] The technical problem to be solved by the present invention is to provide a static penetration test parameter inversion method, system, equipment and medium, which can avoid unnecessary repeated tests, reduce resource waste and improve the efficiency of engineering geological survey.

[0005] In order to solve the above technical problems, an embodiment of the present invention provides a static penetration test parameter inversion method, comprising:

[0006] Obtain the basic parameters of the soil in the target area according to the set variables;

[0007] According to the basic parameters of the soil in the target area, a numerical simulation of the static penetration test is performed based on the material point method to obtain a simulation result of the static penetration test;

[0008] Using the static penetration test simulation results to perform regression analysis, a regression equation of the static penetration test simulation results is obtained;

[0009] The regression analysis was carried out using multiple groups of static penetration test results to obtain the regression equation of the static penetration test results.

[0010] Inversion is performed based on the regression equation of the static penetration test simulation results and the regression equation of the static penetration test actual measurement results to obtain the static penetration test inversion parameters of the soil in the target area.

[0011] Preferably, the basic parameters of the soil include soil moisture content, density, specific gravity, internal friction angle and cohesion.

[0012] Furthermore, the set variables include a cohesive soil moisture content variable, a cohesive soil moisture content variable, and a density variable;

[0013] The variable of the moisture content of the clay soil is sampled at ten gradients between 2% and the liquid limit of the moisture content of the clay soil;

[0014] The moisture content variable of the cohesionless soil is between 2% and 10% of the moisture content of the cohesionless soil, and the samples are taken in five gradients;

[0015] The density variable is sampled in five gradients between the maximum dry density and the minimum dry density.

[0016] Preferably, the numerical simulation of the static penetration test is performed based on the material point method according to the basic parameters of the soil in the target area to obtain the simulation result of the static penetration test, including:

[0017] According to the basic parameters of the soil in the target area, the interaction force between the particles is calculated using the three-phase single-point formula in the material point method to simulate the three-phase structure of the soil in soil mechanics;

[0018] Based on the three-phase structure of the soil, a two-dimensional axisymmetric method is used to model the static penetration test;

[0019] The corresponding parameters and values ​​are set for the established model, and the simulation results of the static penetration test are output.

[0020] Preferably, the established model includes a probe model, a probe rod material model and a soil model;

[0021] Wherein, the setting parameters of the probe model include the width, height and top angle of the probe model;

[0022] The probe rod material model adopts a linear elastic model, and the setting parameters of the linear elastic model include Young's modulus and Poisson's ratio;

[0023] The soil model adopts the Mohr-Coulomb model, and the setting parameters of the Mohr-Coulomb model include basic parameters of the soil.

[0024] Preferably, the static penetration test parameter inversion method further includes:

[0025] The specific penetration resistance value and the specific penetration resistance value error boundary are calculated according to the inversion parameters of the static penetration test.

[0026] In order to solve the above technical problems, an embodiment of the present invention further provides a static penetration test parameter inversion system, comprising:

[0027] The basic soil parameter acquisition module is used to obtain the basic soil parameters of the target area according to the set variables;

[0028] A static penetration test simulation module is used to perform numerical simulation of the static penetration test based on the material point method according to the basic parameters of the soil in the target area, and obtain the simulation results of the static penetration test;

[0029] A first regression analysis module is used to perform regression analysis using the static penetration test simulation results to obtain a regression equation of the static penetration test simulation results;

[0030] The second regression analysis module is used to perform regression analysis using multiple groups of static penetration test measured results to obtain a regression equation of the static penetration test measured results;

[0031] The parameter inversion module is used to perform inversion according to the regression equation of the static penetration test simulation results and the regression equation of the static penetration test actual measurement results to obtain the static penetration test inversion parameters of the soil in the target area.

[0032] Preferably, the static penetration test simulation module comprises:

[0033] The first simulation unit is used to calculate the interaction force between the mass points according to the basic parameters of the soil in the target area by using the three-phase single-point formula in the material point method to simulate the three-phase structure of the soil in soil mechanics;

[0034] A modeling unit, used for modeling a static penetration test using a two-dimensional axisymmetric method based on the three-phase structure of the soil;

[0035] The second simulation unit is used to set corresponding parameters and values ​​for the established model and output the simulation results of the static penetration test.

[0036] Furthermore, the established model includes a probe model, a probe rod material model and a soil model;

[0037] Wherein, the setting parameters of the probe model include the width, height and top angle of the probe model;

[0038] The probe rod material model adopts a linear elastic model, and the setting parameters of the linear elastic model include Young's modulus and Poisson's ratio;

[0039] The soil model adopts the Mohr-Coulomb model, and the setting parameters of the Mohr-Coulomb model include basic parameters of the soil.

[0040] Preferably, the static penetration test parameter inversion system further includes:

[0041] The soil engineering property index calculation module is used to calculate the specific penetration resistance value and the specific penetration resistance value error boundary according to the static penetration test inversion parameters.

[0042] To solve the above technical problems, an embodiment of the present invention further provides a terminal device, comprising a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, wherein when the processor executes the computer program, the static penetration test parameter inversion method as described in any one of the above items is implemented.

[0043] To solve the above technical problems, an embodiment of the present invention further provides a computer-readable storage medium, wherein the computer-readable storage medium includes a stored computer program, wherein when the computer program is running, the device where the computer-readable storage medium is located is controlled to execute the static penetration test parameter inversion method as described in any one of the above.

[0044] Compared with the prior art, the embodiments of the present invention provide a method, system, device and medium for inversion of static penetration test parameters. First, the basic parameters of the soil in the target area are obtained according to the set variables, and according to the basic parameters of the soil in the target area, the numerical simulation of the static penetration test is performed based on the material point method to obtain the simulation results of the static penetration test, and then the simulation results of the static penetration test are used for regression analysis to obtain the regression equation of the simulation results of the static penetration test, and multiple groups of measured results of the static penetration test are used for regression analysis to obtain the regression equation of the measured results of the static penetration test, and finally, the regression equation of the simulation results of the static penetration test and the regression equation of the measured results of the static penetration test are used for inversion to obtain the static penetration test inversion parameters of the soil in the target area. The present invention simulates the static penetration process by combining the basic parameters of the soil with the material point method, and performs regression analysis on the numerical simulation results and the field measured results to obtain the inversion parameters. When survey and design parameters outside the site are needed and actual static penetration tests are not required according to specifications, the present invention can be used to obtain engineering property indicators of soil, thereby avoiding unnecessary repeated tests, reducing the waste of resources caused by repeated tests, and thus improving the efficiency of engineering geological surveys. BRIEF DESCRIPTION OF THE DRAWINGS

[0045] In order to more clearly illustrate the technical solution of the present invention, the drawings used in the implementation mode will be briefly introduced below. Obviously, the drawings described below are only some implementation modes of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0046] Figure 1 It is a schematic flow chart of a static penetration test parameter inversion method provided by an embodiment of the present invention;

[0047] Figure 2 It is a static penetration curve and its regression analysis diagram provided by an embodiment of the present invention;

[0048] Figure 3 The static penetration curve and regression analysis diagram of the geological survey report provided by the embodiment of the present invention;

[0049] Figure 4 It is a structural block diagram of a static penetration test parameter inversion system provided by an embodiment of the present invention;

[0050] Figure 5 It is a structural block diagram of a terminal device provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0051] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0052] See also Figure 1 , Figure 1 1 is a flow chart of a static penetration test parameter inversion method provided by an embodiment of the present invention, wherein the static penetration test parameter inversion method comprises steps S1 to S5:

[0053] S1. Obtain the basic parameters of the soil in the target area according to the set variables;

[0054] In an optional embodiment, the basic parameters of the soil include the moisture content, density, specific gravity, internal friction angle and cohesion of the soil.

[0055] Furthermore, the set variables include a cohesive soil moisture content variable, a cohesive soil moisture content variable, and a density variable;

[0056] The variable of the moisture content of the clay soil is sampled at ten gradients between 2% and the liquid limit of the moisture content of the clay soil;

[0057] The moisture content variable of the cohesionless soil is between 2% and 10% of the moisture content of the cohesionless soil, and the samples are taken in five gradients;

[0058] The density variable is sampled in five gradients between the maximum dry density and the minimum dry density.

[0059] Specifically, soil samples are retrieved from the target area where the static penetration test is required, and then the moisture content, density, and specific gravity of the original soil are determined according to tests such as moisture content test, specific gravity test, and direct shear test. Among them, the variable of density can be set to 5 gradients between the maximum dry density and the minimum dry density, the variable of moisture content of cohesive soil can be set to 10 gradients between 2% and liquid limit, and the variable of moisture content of cohesive soil can be set to 5 gradients between 2% and 10%, and then the internal friction angle and cohesion can be determined respectively.

[0060] It should be noted that the moisture content test in the embodiment of the present invention adopts the drying method. Considering that the climate in some areas is dry and the temperature difference is large, taking the soil sample back to the laboratory for moisture content test will lead to excessive errors. Therefore, the embodiment of the present invention puts the soil sample of the target area into the weighing box, covers the box lid, and immediately weighs the mass of the humidified soil in the weighing box. At the same time, in order to ensure that the error is not very large, three samples can be taken from each sampling point and put into the weighing box, and then the average value is calculated. The density test adopts the ring knife method. After the ring knife sampling is carried out according to the requirements of the geotechnical test, the sample is taken back to the laboratory to weigh the mass and then the density is calculated. In order to minimize the error, four ring knife samples can be taken at each location. Specific gravity test: The test is carried out using the specific gravity bottle method, and the test steps are carried out according to the requirements of the geotechnical test. Direct shear test: A quick shear test was used to prepare the dry density, maximum dry density and minimum dry density of the three soil samples in their natural state. The moisture content of the cohesive soil was varied in 10 gradients from 2% to the liquid limit, and the moisture content of the non-cohesive soil was varied in 5 gradients from 2% to 10%. Then, quick shear tests were carried out in accordance with the requirements of geotechnical tests to measure the internal friction angle and cohesion.

[0061] It is worth mentioning that, through the above method, more accurate basic parameters of the soil can be obtained, thereby further improving the accuracy of the numerical simulation of the static penetration test.

[0062] S2. According to the basic parameters of the soil in the target area, numerical simulation of the static penetration test is performed based on the material point method to obtain simulation results of the static penetration test;

[0063] In an optional embodiment, the numerical simulation of the static penetration test is performed based on the material point method according to the basic parameters of the soil in the target area to obtain the simulation result of the static penetration test, including:

[0064] According to the basic parameters of the soil in the target area, the interaction force between the particles is calculated using the three-phase single-point formula in the material point method to simulate the three-phase structure of the soil in soil mechanics;

[0065] Based on the three-phase structure of the soil, a two-dimensional axisymmetric method is used to model the static penetration test;

[0066] The corresponding parameters and values ​​are set for the established model, and the simulation results of the static penetration test are output.

[0067] Furthermore, the established model includes a probe model, a probe rod material model and a soil model;

[0068] Wherein, the setting parameters of the probe model include the width, height and top angle of the probe model;

[0069] The probe rod material model adopts a linear elastic model, and the setting parameters of the linear elastic model include Young's modulus and Poisson's ratio;

[0070] The soil model adopts the Mohr-Coulomb model, and the setting parameters of the Mohr-Coulomb model include basic parameters of the soil.

[0071] Specifically, the calculation formula for numerical simulation based on the material point method adopts the three-phase single-point formula, and the principle is shown in formulas (1)(2)(3):

[0072]

[0073] In formula (1), (2), and (3), is the node acceleration and velocity vector of each phase, M G 、M L , and is the lumped mass matrix of gas, liquid and solid, is the internal and external nodal force vector of gas, liquid phase and mixture, Q G and Q L is the resistance matrix.

[0074] It is worth noting that the three-phase single-point formula assumes the same conception of soil structure as the three-phase structure of soil in soil mechanics. In the three-phase MPM (Material Point Method) framework, the solid mass in each MPM remains unchanged, so the solid mass is automatically conserved. However, liquids and gases can leave the solid skeleton, so the conservation of fluid mass depends entirely on the momentum balance of liquids and gases. In contrast to the two-phase single-point method, the three-phase single-point formula does not ignore the flux caused by the spatial variation of liquid and gas mass, and takes into account all terms of mass balance.

[0075] It is worth noting that the material point method is a complete Lagrangian particle method. In each step, there is no relative motion between the particle and the computational grid, which avoids the numerical difficulties caused by the nonlinear convection term of the Euler method and makes it easy to track the material interface. For the computational grid of the material point method, the discrete momentum balance equation is usually solved at the nodes of the computational grid, while the mass conservation equation and the constitutive equation are usually solved at the material points of the computational grid. The information associated with the grid is not needed for the next step of analysis, so it can be discarded, so any grid deformation can be avoided. That is, it not only does not need to deal with the convection term during the simulation process, but also avoids the problem of grid distortion, making it very suitable for the simulation of static penetration test, and the parameters are dynamically changed during the simulation process, so the simulation can know the range of change.

[0076] Furthermore, the static penetration test model was constructed based on the three-phase results of the soil.

[0077] Specifically, two-dimensional axisymmetric modeling is adopted, and type I-1 single-bridge probe is taken as an example. After modeling, the probe diameter, width, height and top angle of the probe model are set; the probe rod material model is set to a linear elastic model, and the parameters Young's modulus, Poisson's ratio and penetration velocity are input; the soil model adopts the Mohr-Coulomb model, and the parameters are selected according to the basic parameters of different soil layers. For example, the basic parameters of the input soil layer include porosity, density, static earth pressure coefficient, Poisson's ratio, Young's modulus and internal friction angle.

[0078] Among them, the void ratio, density, and static earth pressure coefficient can be calculated from the parameters obtained from indoor geotechnical tests. The calculation formulas are shown in (4)(5)(6):

[0079]

[0080] ρ s =G s ρ w (5);

[0081]

[0082] In formula (4), (5), and (6), G s is the specific gravity, is the internal friction angle, e is the porosity, n is the porosity ratio, ρ s is the density of soil, ρ w is the density of water, γ d is the dry soil weight, γ w is the density of water, and K0 is the static earth pressure coefficient.

[0083] Poisson's ratio can be obtained by consulting the standard; because the compression modulus value of radially isotropic materials is often approximately equal to their Young's modulus value, the compression modulus value can be approximately taken as the Young's modulus value; the internal friction angle is the measured value in the geological survey report.

[0084] After the above parameters are input into the numerical simulation software, a surface is defined in the software to record the reaction force (the recorded reaction force is the output result), namely the cone tip resistance and the cone side resistance. The numerical simulation software calculates the result through its built-in formula and finally outputs the result. The specific penetration resistance is converted from the cone tip resistance according to the formula. The formula is shown in (7). Considering that the measured result is only the specific penetration resistance and depth curve, and the specific penetration resistance can be converted to the cone tip resistance, the cone side resistance and depth curve will not be analyzed in the following.

[0085] P s =1.093q c +0.358(MPa)(7);

[0086] In formula (7), P s is the relative penetration resistance, q c is the cone tip resistance;

[0087] The final static penetration test simulation result is specifically a static penetration test curve.

[0088] S3, using the static penetration test simulation results to perform regression analysis to obtain a regression equation for the static penetration test simulation results;

[0089] S4. Perform regression analysis using multiple groups of measured results of static penetration tests to obtain a regression equation for the measured results of static penetration tests;

[0090] S5. Invert the static penetration test inversion parameters of the target area soil according to the regression equation of the static penetration test simulation results and the regression equation of the static penetration test actual measurement results.

[0091] In an optional embodiment, the static penetration test parameter inversion method further includes:

[0092] The specific penetration resistance value and the specific penetration resistance value error boundary are calculated according to the inversion parameters of the static penetration test.

[0093] For example, the outputted static penetration test curve graph is subjected to regression analysis, and at the same time, no less than 14 groups of static penetration test measured curve graphs are subjected to regression analysis, and the regression equation of the static penetration test curve and the regression equation of the static penetration test measured curve are obtained respectively. 模拟 and f 实际 and through Finally, the inversion parameters K1K2K3 are obtained, and the specific penetration resistance value is q s实际=K1q s模拟 , the error margin of the penetration resistance value is K2q s模拟 ~K3q s模拟 .

[0094] In summary, a method for inverting parameters of a static penetration test provided by an embodiment of the present invention first obtains the basic parameters of the soil in the target area according to the set variables, and performs numerical simulation of the static penetration test based on the material point method according to the basic parameters of the soil in the target area to obtain the simulation results of the static penetration test, and then uses the simulation results of the static penetration test to perform regression analysis to obtain the regression equation of the simulation results of the static penetration test, and uses multiple groups of measured results of static penetration test to perform regression analysis to obtain the regression equation of the measured results of the static penetration test, and finally performs inversion according to the regression equation of the simulation results of the static penetration test and the regression equation of the measured results of the static penetration test to obtain the static penetration test inversion parameters of the soil in the target area. The present invention simulates the static penetration test process by combining the basic parameters of the soil with the material point method, and performs regression analysis on the numerical simulation results and the field measured results to obtain the inversion parameters. When survey and design parameters outside the site are needed and actual static penetration tests are not required according to specifications, the present invention can be used to obtain engineering property indicators of soil, thereby avoiding unnecessary repeated tests, reducing the waste of resources caused by repeated tests, and thus improving the efficiency of engineering geological surveys.

[0095] In order to make those skilled in the art more aware of the static penetration test parameter inversion method provided by the present embodiment, the static penetration test parameter inversion method provided by the present embodiment is described in detail below in combination with data during specific implementation.

[0096] Specifically, taking the static penetration test of a photovoltaic project in a wind-blown sand field in a desert area as an example, the static penetration test parameter inversion is carried out.

[0097] Step 1: Obtain the basic parameters of the soil in the target desert area according to the set variables:

[0098] Specifically, after retrieving soil samples in the area where the static penetration test is conducted, the moisture content, density, specific gravity, internal friction angle and cohesion of the undisturbed soil are measured, and the density variable is set to three gradients of maximum dry density, minimum dry density and natural dry density, and the moisture content is respectively proportioned to 4%, 6%, 8% and 10%. Then, quick shear tests are carried out in accordance with the requirements of geotechnical tests to measure the internal friction angle and cohesion, respectively.

[0099] The moisture content of the soil was obtained by a moisture content test using a drying method. During the test, considering that the climate in the desert is dry and the temperature difference is large, taking the soil samples back to the laboratory for moisture content test will result in large errors. Therefore, the inventor selected soil samples from three sampling locations and put them into a weighing box, closed the box lid, and immediately weighed the mass of the humidified soil in the weighing box. Three samples were taken from each sampling point and put into the weighing box. Finally, the average value was calculated, and the moisture contents of the three sampling points were finally obtained to be 0.0280172, 0.0243789, and 0.0471074, respectively.

[0100] The density of the soil was obtained by the density test using the knife ring method. Specifically, after knife ring sampling according to the geotechnical test requirements, the samples were taken back to the laboratory to weigh the mass and then the density was calculated. In order to minimize the error, four knife ring samples were taken at each sampling location. The final density of the soil at the three sampling points was 1.35, 1.47, and 1.48 (g / cm 3 ).

[0101] The specific gravity of the soil is obtained by a specific gravity test using a pycnometer method. Specifically, the test steps are carried out according to the requirements of geotechnical tests, and the specific gravity of the soil at the three sampling points in the embodiment of the present invention is 2.66, 2.67, and 2.56 respectively.

[0102] Step 2: According to the basic parameters of the desert soil, numerical simulation of the static penetration test is carried out based on the material point method:

[0103] The built-in calculation formula in the numerical simulation software based on the material point method adopts the three-phase single-point formula. The assumption of soil structure in this formula is the same as the three-phase structure of soil in soil mechanics, and can be used to simulate the three-phase structure of soil.

[0104] Furthermore, the static penetration test model was constructed based on the three-phase results of the soil.

[0105] Specifically, two-dimensional axisymmetric modeling is adopted. Taking the type Ⅰ-1 single-bridge probe as an example, after modeling, the probe diameter in the probe model is set to 0.0357m, the width of the probe model is 0.0179m, the height is 0.0309m, and the top angle is 60°; the probe rod material model is set to a linear elastic model, and the input parameters are Young's modulus of 50000, Poisson's ratio of 0.33, and the penetration speed is a uniform penetration speed of 0.02m / s; the soil model adopts the Mohr-Coulomb model, and the parameters are selected according to the basic parameters of different soil layers. By way of example, the basic parameters of the input soil layer include internal friction angle, porosity, density, static earth pressure coefficient, Poisson's ratio, and Young's modulus.

[0106] Specifically, the internal friction angle, void ratio, density, and static earth pressure coefficient can be calculated from the parameters obtained from indoor geotechnical tests. The calculation formula is as follows (4)(5)(6). The internal friction angle finally calculated is is 27°, the void ratio n is 0.51, and the density of the soil is ρ s 2660kg / m 3 , the coefficient of static earth pressure K0 is 0.55; Poisson's ratio can be found in the standard and is 0.15; Young's modulus is set to 5000 kPa.

[0107] Then, the reaction force is calculated by defining the cone tip resistance and the cone side resistance, and then these two values ​​are output through numerical simulation, and the specific penetration resistance is converted according to formula (7), and finally the static penetration curve is output. Figure 2 , Figure 2 It is a static penetration test curve and its regression analysis diagram provided by the embodiment of the present invention. Figure 3 Analysis shows that the trends of the static penetration curve output by the numerical simulation based on the material point method and the curve of the static penetration in the geological survey report are similar in most aspects, and the specific penetration resistance obtained by the numerical simulation is very close to the standard value of the specific penetration resistance in the geological survey report. The specific penetration resistance of the first layer of soil in the numerical simulation in the embodiment of the present invention is between 0 and 1.2 MPa, and the specific penetration resistance of the second layer of soil in the numerical simulation is between 0.8 and 1.2 MPa, while the standard value in the geological survey report is 1 MPa. The specific penetration resistance of the third layer of soil in the numerical simulation is between 0.9 and 1.5 MPa, and the standard value in the geological survey report is 1.4 MPa. The specific penetration resistance of the fourth layer of soil in the numerical simulation is between 0.8 and 1.2 MPa, and the standard value in the geological survey report is 1.8 MPa. That is, the static penetration curve output by the numerical simulation is consistent with the law of the measured static penetration curve, and the specific penetration resistance of the same soil layer is very close to the measured standard value.

[0108] Step 3: Perform regression analysis on the simulation results of the static penetration test and multiple groups of measured data of the static penetration test to obtain the regression equation of the simulation results of the static penetration test and the regression equation of the measured results of the static penetration test:

[0109] Specifically, the static penetration test curve output in step 2 is subjected to regression analysis, and no less than 14 sets of static penetration test measured data are subjected to regression analysis to obtain the regression equation of the static penetration test curve and the regression equation of the static penetration test measured curve f respectively. 模拟 and f 实际 and through Finally, the inversion parameters K1K2K3 are obtained.

[0110] Among them, Figure 2 As shown in the figure, the regression analysis of the numerical simulation static penetration test curve adopts the exponential form fitting, and the R 2 The value is 0.99866>0.7. According to the statistical requirement, when it is not less than 0.7, the fitting effect is good. 2R-squared, also known as the coefficient of determination, is an important concept in statistics. It is mainly used in regression analysis to explain the proportion of the variation of the response variable explained by the variation of the independent variable. 2 The values ​​range from 0 to 1, with higher values ​​indicating a better fit to the data.

[0111] Further, see Figure 3 , Figure 3 The following is a static penetration curve and regression analysis diagram of the geological survey report provided by the embodiment of the present invention. Figure 3 As shown in the figure, the measured static penetration test selected 14 drilling points (ZK01, ZK04, ZK07, ZK15, ZK18, ZK21, ZK29, ZK32, ZK35, ZK36, ZK41, ZK44, ZK48, and ZK55), and put their static penetration test curves in the same curve graph for regression analysis. The linear fitting with boundaries was used, and the R of the regression analysis fitting function was 0. 2 The value 0.82725>0.7. Similarly, according to the statistical requirement, it should be no less than 0.7 to indicate a good fitting effect.

[0112] Step 4: When no actual static penetration test is conducted in the aeolian sand area, the specific penetration resistance value can be calculated based on the obtained inversion parameters K1K2K3, that is, the specific penetration resistance value q s实际 =K1q s模拟 , the error margin of the penetration resistance value is K2q s模拟 ~K3q s模拟 .

[0113] Based on the above method embodiments, the embodiments of the present invention provide corresponding system embodiments.

[0114] See also Figure 4 , Figure 4 1 is a structural block diagram of a static penetration test parameter inversion system provided by an embodiment of the present invention, wherein the static penetration test parameter inversion system comprises:

[0115] The basic soil parameter acquisition module 21 is used to acquire the basic soil parameters of the target area according to the set variables;

[0116] The static penetration test simulation module 22 is used to perform numerical simulation of the static penetration test based on the material point method according to the basic parameters of the soil in the target area, and obtain the simulation result of the static penetration test;

[0117] A first regression analysis module 23, used to perform regression analysis using the static penetration test simulation results to obtain a regression equation of the static penetration test simulation results;

[0118] A second regression analysis module 24 is used to perform regression analysis using multiple groups of static penetration test measured results to obtain a regression equation of the static penetration test measured results;

[0119] The parameter inversion module 25 is used to perform inversion according to the regression equation of the static penetration test simulation result and the regression equation of the static penetration test actual measurement result to obtain the static penetration test inversion parameters of the soil in the target area.

[0120] In an optional embodiment, the static penetration test simulation module includes:

[0121] The first simulation unit is used to calculate the interaction force between the mass points according to the basic parameters of the soil in the target area by using the three-phase single-point formula in the material point method to simulate the three-phase structure of the soil in soil mechanics;

[0122] A modeling unit, used for modeling a static penetration test using a two-dimensional axisymmetric method based on the three-phase structure of the soil;

[0123] The second simulation unit is used to set corresponding parameters and values ​​for the established model and output the simulation results of the static penetration test.

[0124] Furthermore, the established model includes a probe model, a probe rod material model and a soil model;

[0125] Wherein, the setting parameters of the probe model include the width, height and top angle of the probe model;

[0126] The probe rod material model adopts a linear elastic model, and the setting parameters of the linear elastic model include Young's modulus and Poisson's ratio;

[0127] The soil model adopts the Mohr-Coulomb model, and the setting parameters of the Mohr-Coulomb model include basic parameters of the soil.

[0128] In an optional embodiment, the static penetration test parameter inversion system further includes:

[0129] The soil engineering property index calculation module is used to calculate the specific penetration resistance value and the specific penetration resistance value error boundary according to the static penetration test inversion parameters.

[0130] It should be noted that a static penetration test parameter inversion system provided in an embodiment of the present invention is used to execute all process steps of a static penetration test parameter inversion method in the above embodiment. The working principles and beneficial effects of the two correspond one to one, and thus will not be described in detail.

[0131] The embodiment of the present invention also provides a terminal device, such as Figure 5, which is a block diagram of a terminal device provided by an embodiment of the present invention. The terminal device includes a processor 31, a memory 32, and a computer program stored in the memory 32 and configured to be executed by the processor 31. When the processor 31 executes the computer program, the static penetration test parameter inversion method described in any of the above embodiments is implemented.

[0132] In addition, an embodiment of the present invention also provides a computer-readable storage medium, which includes a stored computer program, wherein when the computer program is running, the device where the computer-readable storage medium is located is controlled to execute the static penetration test parameter inversion method described in any of the above embodiments.

[0133] When the processor 31 executes the computer program, the steps in the above-mentioned static penetration test parameter inversion method embodiment are implemented, for example: Figure 1 Alternatively, when the processor 31 executes the computer program, the functions of each module in the above-mentioned static penetration test parameter inversion system embodiment are realized, for example Figure 4 The functions of each module of the static penetration test parameter inversion system are shown.

[0134] Preferably, the computer program can be divided into one or more modules / units, which are stored in the memory 32 and executed by the processor 31 to implement the present invention. The one or more modules / units can be a series of computer program instruction segments that can implement specific functions, and the instruction segments are used to describe the execution process of the computer program in the terminal device.

[0135] The processor 31 can be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor can be a microprocessor, or the processor 31 can be any conventional processor. The processor 31 is the control center of the terminal device, and various parts of the terminal device are connected using various interfaces and lines.

[0136] The memory 32 mainly includes a program storage area and a data storage area, wherein the program storage area can store an operating system, an application required for at least one function, etc., and the data storage area can store related data, etc. In addition, the memory 32 can be a high-speed random access memory, or a non-volatile memory, such as a plug-in hard disk, a smart memory card (Smart Media Card, SMC), a secure digital (Secure Digital, SD) card, and a flash card (Flash Card), etc., or the memory 32 can also be other volatile solid-state storage devices.

[0137] It should be noted that the above terminal device may include, but is not limited to, a processor and a memory. Those skilled in the art will understand that Figure 5 The structural block diagram shown is only a structural example of the above-mentioned terminal device and does not constitute a structural limitation of the above-mentioned terminal device. The above-mentioned terminal device may include more or fewer components than shown in the figure, or combine certain components, or different components.

[0138] The above is a preferred embodiment of the present invention. It should be pointed out that a person skilled in the art can make several improvements and modifications without departing from the principle of the present invention. These improvements and modifications are also considered to be within the scope of protection of the present invention.

Claims

1. A static penetration test parameter inversion method, characterized in that: include: Obtain the basic parameters of the soil in the target area according to the set variables; According to the basic parameters of the soil in the target area, a numerical simulation of the static penetration test is performed based on the material point method to obtain a simulation result of the static penetration test; Using the static penetration test simulation results to perform regression analysis, a regression equation of the static penetration test simulation results is obtained; The regression analysis was carried out using multiple groups of static penetration test results to obtain the regression equation of the static penetration test results. Inversion is performed based on the regression equation of the static penetration test simulation results and the regression equation of the static penetration test actual measurement results to obtain the static penetration test inversion parameters of the soil in the target area.

2. The static penetration test parameter inversion method according to claim 1, characterized in that: The basic parameters of the soil include soil moisture content, density, specific gravity, internal friction angle and cohesion.

3. The static penetration test parameter inversion method according to claim 1, characterized in that: The set variables include a cohesive soil moisture content variable, a cohesive soil moisture content variable, and a density variable; The variable of the moisture content of the clay soil is sampled at ten gradients between 2% and the liquid limit of the moisture content of the clay soil; The moisture content variable of the cohesionless soil is between 2% and 10% of the moisture content of the cohesionless soil, and the samples are taken in five gradients; The density variable is sampled in five gradients between the maximum dry density and the minimum dry density.

4. The static penetration test parameter inversion method according to claim 1, characterized in that: The numerical simulation of the static penetration test is performed based on the material point method according to the basic parameters of the soil in the target area to obtain the simulation results of the static penetration test, including: According to the basic parameters of the soil in the target area, the interaction force between the mass points is calculated using the three-phase single-point formula in the material point method to simulate the three-phase structure of the soil in soil mechanics; Based on the three-phase structure of the soil, a two-dimensional axisymmetric method is used to model the static penetration test; The corresponding parameters and values ​​are set for the established model, and the simulation results of the static penetration test are output.

5. The static penetration test parameter inversion method according to claim 4, characterized in that: The established models include probe model, probe rod material model and soil model; Wherein, the setting parameters of the probe model include the width, height and top angle of the probe model; The probe rod material model adopts a linear elastic model, and the setting parameters of the linear elastic model include Young's modulus and Poisson's ratio; The soil model adopts the Mohr-Coulomb model, and the setting parameters of the Mohr-Coulomb model include basic parameters of the soil.

6. The static penetration test parameter inversion method according to claim 1, characterized in that: Also includes: The specific penetration resistance value and the specific penetration resistance value error boundary are calculated according to the inversion parameters of the static penetration test.

7. A static penetration test parameter inversion system, characterized in that: include: The basic soil parameter acquisition module is used to obtain the basic soil parameters of the target area according to the set variables; A static penetration test simulation module is used to perform numerical simulation of the static penetration test based on the material point method according to the basic parameters of the soil in the target area, and obtain the simulation results of the static penetration test; A first regression analysis module is used to perform regression analysis using the static penetration test simulation results to obtain a regression equation of the static penetration test simulation results; The second regression analysis module is used to perform regression analysis using multiple groups of static penetration test measured results to obtain a regression equation of the static penetration test measured results; The parameter inversion module is used to perform inversion according to the regression equation of the static penetration test simulation results and the regression equation of the static penetration test actual measurement results to obtain the static penetration test inversion parameters of the soil in the target area.

8. The static penetration test parameter inversion system according to claim 7, characterized in that: The static penetration test simulation module comprises: The first simulation unit is used to calculate the interaction force between the mass points according to the basic parameters of the soil in the target area by using the three-phase single-point formula in the material point method to simulate the three-phase structure of the soil in soil mechanics; A modeling unit, used for modeling a static penetration test using a two-dimensional axisymmetric method based on the three-phase structure of the soil; The second simulation unit is used to set corresponding parameters and values ​​for the established model and output the simulation results of the static penetration test.

9. The static penetration test parameter inversion system according to claim 8, characterized in that: The established models include probe model, probe rod material model and soil model; Wherein, the setting parameters of the probe model include the width, height and top angle of the probe model; The probe rod material model adopts a linear elastic model, and the setting parameters of the linear elastic model include Young's modulus and Poisson's ratio; The soil model adopts the Mohr-Coulomb model, and the setting parameters of the Mohr-Coulomb model include basic parameters of the soil.

10. The static penetration test parameter inversion system according to claim 7, characterized in that: Also includes: The soil engineering property index calculation module is used to calculate the specific penetration resistance value and the specific penetration resistance value error boundary according to the static penetration test inversion parameters.

11. A terminal device, characterized in that: The method comprises a processor, a memory and a computer program stored in the memory and configured to be executed by the processor, wherein when the processor executes the computer program, the static penetration test parameter inversion method as described in any one of claims 1 to 6 is implemented.

12. A computer-readable storage medium, characterized in that: The computer-readable storage medium includes a stored computer program, wherein when the computer program is executed, the device where the computer-readable storage medium is located is controlled to execute the static penetration test parameter inversion method according to any one of claims 1 to 6.

Images