A calculation method for dynamic residual deformation of saturated soil and its implementation method

Abstract

The present invention discloses a method for calculating the dynamic residual deformation of saturated soil and an implementation method thereof; it relates to the field of dynamic characteristics of geotechnical engineering, and the calculation method comprises: S1, establishing a calculation model for the dynamic residual deformation of saturated soil; S2, obtaining existing original test data; S3, processing the test data to obtain fitting data; S4, training and fitting the calculation model to obtain a trained calculation model; S5, calculating the dynamic residual deformation of saturated soil according to the calculation model; the implementation method comprises: S1, establishing a worksheet in Excel; S2, inputting the calculation model for the dynamic residual deformation of saturated soil; S3, calculating the parameters of the sample after consolidation, the axial dynamic stress and axial dynamic strain; S4. Filter the data in step S3 and output them to the data processing worksheet; S5. In the data processing worksheet, calculate the residual shear strain and residual volume strain; S6. Use the regression analysis calculation formula to calculate the model parameters and the determination coefficient R to obtain the calculation model; S7. Use the calculation model to calculate the dynamic residual deformation of saturated soil; This application first improves the original saturated soil dynamic residual deformation calculation model, and then uses a large amount of data to fit the model so that the model conforms to the data obtained from the actual experiment. The model can be used to measure the relationship between the vibration cycle of the soil and the dynamic residual deformation.

Description

Technical Field

[0001] The present invention relates to the field of geotechnical engineering dynamic characteristics, and in particular to a calculation method for dynamic residual deformation of saturated soil and an implementation method thereof. Background Art

[0002] Dynamic triaxial testing involves applying a dynamic load along the axial direction of a cylindrical soil sample under a constant confining pressure. The dynamic stress, dynamic strain, and pore water pressure are measured over time to determine the soil's dynamic characteristic parameters and dynamic residual deformation. In civil engineering structures such as buildings, bridges, and roads, dynamic residual deformation of the soil directly affects foundation stability, leading to uneven settlement and, in turn, compromising the safety and service life of the structure. Calculating dynamic residual deformation of the soil can be used to assess foundation stability and ensure the safety of civil engineering structures.

[0003] Previous studies have investigated dynamic residual deformation and proposed models to describe the relationship between soil vibration cycles and dynamic residual deformation. The Shen Zhujiang model, however, suggests that the residual deformation of earth-rock dam materials increases linearly on a semi-logarithmic scale with increasing vibration cycles. However, experiments have shown that the dynamic residual deformation of saturated soil exhibits a nonlinear relationship on a semi-logarithmic scale with increasing vibration cycles, significantly different from the Shen Zhujiang model. Consequently, using this model to calculate the dynamic residual deformation of saturated soil can result in certain errors. Summary of the Invention

[0004] In order to calculate the dynamic residual deformation of saturated soil more accurately, the present application provides a method for calculating the dynamic residual deformation of saturated soil and an implementation method thereof.

[0005] In a first aspect, the present application provides a method for calculating the dynamic residual deformation of saturated soil, which adopts the following technical solution:

[0006] A method for calculating the dynamic residual deformation of saturated soil includes the following steps:

[0007] S1. Establish a calculation model for dynamic residual deformation of saturated soil;

[0008] S2. Obtain existing raw test data;

[0009] S3, processing the test data obtained in S2 to obtain fitting data for training the model;

[0010] S4, using the fitting data obtained in S3, the computational model in S1 is trained and fitted to obtain a trained computational model;

[0011] S5. Calculate the dynamic residual deformation of saturated soil according to the trained calculation model obtained in S4.

[0012] Optionally, the calculation model for the dynamic residual deformation of saturated soil in step S1 is:

[0013] γ r =a·log b (N+1)

[0014] ε vr = c·log d (N+1);

[0015] where γ r is the residual shear strain, ε vr is the residual volume strain, a, b, c, d are the model fitting parameters to be determined, and N is the number of vibrations.

[0016] Optionally, the test data to be acquired in step S2 include: sample starting height h0, sample starting volume V0, sample starting area A0, sample height after consolidation h c , consolidation drainage ΔV, sample area after consolidation A c ;

[0017] Among them, the height of the sample after consolidation h c

[0018]

[0019] Sample area after consolidation A c

[0020]

[0021] Sample volume after consolidation V c

[0022] V c =h c A c .

[0023] Optionally, the fitting data in step S3 includes:

[0024] Axial dynamic strain ε d

[0025]

[0026] Where Δh d is the axial dynamic deformation, in millimeters, h c is the height of the specimen after consolidation, in millimeters;

[0027] Dynamic volume strain ε V

[0028]

[0029] Dynamic shear strain γ d

[0030] γ d =εd (1+μ)

[0031] Where μ is Poisson's ratio;

[0032] Residual axial strain ε dr

[0033] ε dr =ε di -ε d0

[0034] Where ε di is the axial strain of the i-th vibration, ε d0 is the initial axial strain;

[0035] Residual shear strain γ r

[0036] γ r =γ i -γ0

[0037] Where, γ i is the shear strain of the ith vibration, γ0 is the initial shear strain;

[0038] Residual volume strain ε vr

[0039] ε vr =ε vi -ε v0

[0040] Where, ε vi is the volume strain of the ith vibration, ε v0 is the initial volume strain.

[0041] Optionally, in step S4, the data obtained in S3 are used to create an X dataset, a Y dataset, and a Z dataset, where the data in the X dataset is log(log(N+1)), and the data in the Y dataset is logγ r , the data in the Z data set is logε vr ,but:

[0042]

[0043]

[0044] Where x is the sample in the X data set, y is the sample in the Y data set, z is the sample in the Z data set, and n is the sample size. is the sample mean of the X data set, is the sample mean of the Y data set; is the sample mean of the Z dataset.

[0045] Optionally, in step S4, R1 and R2 are the determination coefficients of the saturated soil dynamic residual deformation calculation model. When |R1| ≥ 0.9, the fitting of a and b is stopped, and when |R2| ≥ 0.9, the fitting of c and d is stopped, thereby obtaining a mature calculation model, where

[0046]

[0047] In a second aspect, the present application provides a method for calculating the dynamic residual deformation of saturated soil, which adopts the following technical solution:

[0048] A method for calculating the dynamic residual deformation of saturated soil includes the following steps:

[0049] S1. Create three worksheets in Excel: original data, data processing, and dynamic residual deformation calculation model;

[0050] S2. Input the saturated soil dynamic residual deformation calculation model into the dynamic residual deformation calculation model worksheet;

[0051] S3. Input the dynamic triaxial test data into an Excel raw data worksheet and filter the data to obtain the test data corresponding to the minimum axial force at each vibration frequency; and calculate the post-consolidation specimen parameters, axial dynamic stress, and axial dynamic strain;

[0052] S4. Calculate the specimen parameters, axial dynamic stress, and axial dynamic strain after consolidation in the original data worksheet;

[0053] S5. In the data processing worksheet, calculate the residual shear strain and the residual volume strain based on the test data selected in step S4;

[0054] S6. In the dynamic residual deformation calculation model worksheet, use the regression analysis formula to calculate the parameters and determination coefficient R of the saturated soil dynamic residual deformation calculation model to obtain a mature calculation model;

[0055] S7. Calculate the dynamic residual deformation of saturated soil using the mature calculation model obtained in S6.

[0056] In summary, this application has the following beneficial technical effects:

[0057] This application first improves the original calculation model of saturated soil dynamic residual deformation, and then uses a large amount of data to fit the model to make the model consistent with the data obtained from actual experiments. The model can then be used to measure the relationship between the vibration cycles and dynamic residual deformation of the soil. BRIEF DESCRIPTION OF THE DRAWINGS

[0058] Figure 1 This is the overall flow chart of a method for calculating the dynamic residual deformation of saturated soil in this application;

[0059] Figure 2 This is an overall flow chart of a method for calculating the dynamic residual deformation of saturated soil according to the present application;

[0060] Figure 3 Schematic diagram of basic parameter input in the embodiment;

[0061] Figure 4 Schematic diagram of the original test data import in the embodiment;

[0062] Figure 5 Schematic diagram of screening data in the embodiment

[0063] Figure 6 This is a schematic diagram of the calculation of dynamic residual deformation in the embodiment.

[0064] Figure 7 Schematic diagram of model parameter calculation and model calculation results in the embodiment

[0065] Figure 8-10 is a schematic diagram of the drawing process in the embodiment;

[0066] Figure 11 Schematic diagram of the calculation of dynamic residual deformation model parameters in the embodiment. DETAILED DESCRIPTION

[0067] The following is combined with Figure 1-11 This application is described in further detail.

[0068] The present invention discloses a method for calculating the dynamic residual deformation of saturated soil, comprising the following steps:

[0069] S1. Establish a calculation model for the dynamic residual deformation of saturated soil, where the calculation model for the dynamic residual deformation of saturated soil is:

[0070] γ r =a·log b (N+1)

[0071] ε vr = c·log d (N+1);

[0072] where γ r is the residual shear strain, ε vr is the residual volume strain, a, b, c, d are the model fitting parameters to be determined, and N is the number of vibrations;

[0073] S2. Obtain the existing original test data, including: sample starting height h0, sample starting volume V0, sample starting area A0, sample height after consolidation h c , consolidation drainage ΔV, sample area after consolidation A c ;

[0074] Among them, the height of the sample after consolidation h c

[0075]

[0076] Sample area after consolidation A c

[0077]

[0078] Sample volume after consolidation V c

[0079] V c =h c A c ;

[0080] S3. Process the test data obtained in S2 to obtain fitting data for training the model. The fitting data includes:

[0081] Axial dynamic stress σ d

[0082]

[0083] Where F is the dynamic load;

[0084] Axial dynamic strain ε d

[0085]

[0086] Where Δh d is the axial dynamic deformation, in millimeters, h c is the height of the specimen after consolidation, in millimeters;

[0087] Dynamic volume strain ε V

[0088]

[0089] Dynamic shear strain γ d

[0090] γ d =ε d (1+μ)

[0091] Where μ is Poisson's ratio;

[0092] Residual axial strain ε dr

[0093] ε dr =ε di -ε d0

[0094] Where εdi is the axial strain of the i-th vibration, ε d0 is the initial axial strain;

[0095] Residual shear strain γ r

[0096] γ r =γ i -γ0

[0097] Where, γ i is the shear strain of the ith vibration, γ0 is the initial shear strain;

[0098] Residual volume strain ε vr

[0099] ε vr =ε vi -ε v0

[0100] Where, ε vi is the volume strain of the ith vibration, ε v0 is the initial volume strain;

[0101] S4. Using the fitted data obtained in S3, the data obtained in S3 are used to establish the X data set, the Y data set, and the Z data set. The data in the X data set is log(log(N+1)), and the data in the Y data set is logγ r , the data in the Z data set is logε vr ,but:

[0102]

[0103] Where x is the sample in the X data set, y is the sample in the Y data set, z is the sample in the Z data set, and n is the sample size. is the sample mean of the X data set, is the sample mean of the Y data set; is the sample mean of the Z data set;

[0104] R1 and R2 are the determination coefficients of the calculation model for the dynamic residual deformation of saturated soil. When |R1| ≥ 0.9, the fitting of a and b is stopped. When |R2| ≥ 0.9, the fitting of c and d is stopped. Thus, a mature calculation model is obtained.

[0105]

[0106] S5. Calculate the dynamic residual deformation of saturated soil based on the mature calculation model obtained in S4.

[0107] This application also provides a method for calculating the dynamic residual deformation of saturated soil. The above calculation method can be implemented in Excel, which includes the following steps:

[0108] The first step is to create a new Excel file and create five worksheets: "Instructions", "Raw Data", "Data Processing", "Drawing", and "Dynamic Residual Deformation Calculation Model". The "Instructions" worksheet introduces the specific use of the present invention; the "Raw Data" worksheet is used to import the basic parameters and raw data related to the dynamic triaxial test; the "Data Processing" worksheet filters the data in the "Raw Data" and calculates dynamic stress, dynamic strain, and dynamic residual deformation; the "Drawing" worksheet plots the relevant data calculated in the "Data Processing" worksheet; the "Dynamic Residual Deformation Calculation Model" worksheet proposes a calculation model for the dynamic residual deformation of saturated soil and calculates the model parameters and determination coefficient, where the calculation model is:

[0109] γ r =a·log b (N+1)

[0110] ε vr = c·log d (N+1);

[0111] where γ r is the residual shear strain, ε vr is the residual volume strain, a, b, c, d are the model fitting parameters to be determined, and N is the number of vibrations;

[0112] The second step is Figure 3 As shown, enter the basic parameters of the dynamic triaxial test in the "Raw Data" worksheet. Specifically, enter the specimen starting diameter, specimen starting height h0, consolidation displacement ΔV, number of sampling points per vibration, and Poisson's ratio μ in cells B3, D3, F3, H3, and J3 respectively;

[0113] The third step is to process the basic parameters of the dynamic triaxial test. Specifically, insert the function = 0.25*PI()*B3^2 in cell B4 to obtain the sample starting area A0; insert the function = B4*D3 in cell D4 to obtain the sample starting volume V0; insert the function = D3*(1-F3 / D4)^(1 / 3) in cell F4 to obtain the sample height h after consolidation. c ;Insert the function =B4*(1-F3 / D4)^(2 / 3) in cell B5 to obtain the area A of the sample after consolidation. c ; Insert the function =F4*B5 in cell D5 to obtain the volume V of the sample after consolidation c ;

[0114] The fourth step is to find the minimum axial force of the first vibration based on the number of sampling points per vibration entered in the second step, insert the function =MIN(INDIRECT("B9:B"&(8+H3))) in cell H4, and insert the function =VLOOKUP(H4,B:H,7,FALSE) in cell J4 to output the corresponding serial number of the sampling point with the minimum axial force of the first vibration; it should be noted that the screening process here is a conventional means of obtaining experimental data. In the experimental process, the hysteresis curve is like a spring. A test vibrates 50 times, and each vibration has 50 sampling points, that is, 50 connected circles (one vibration corresponds to one circle, and each circle has 50 data points), a total of 2500 test data. The screening here is to find the point corresponding to the minimum axial force in each vibration, that is, in each circle, and calculate the residual deformation.

[0115] The explanation in yellow above can be replaced with the following:

[0116] Here, the original data of the dynamic triaxial test are screened to determine the test data corresponding to the minimum axial force under each vibration (that is, the test data corresponding to the unloading state of each vibration), and to calculate the dynamic axial strain, dynamic body strain, and dynamic shear strain, and then calculate the dynamic residual body strain and dynamic residual shear strain.

[0117] Combine Figure 4 , perform the following steps in the "Raw Data" worksheet:

[0118] Step 5: Import the dynamic triaxial test data. Specifically, import the dynamic load, dynamic axial deformation, dynamic pore water pressure, confining pressure, dynamic volume change, and sampling time in columns B to G, and import each data item from row 8 to row n (n = number of sampling points per vibration * total vibrations + 8).

[0119] Step 6. Insert the function =B8*10 / $B$5 in cell I8 (*10 here is for unit conversion, i.e. N / cm2 / 10=103N / m2=kPa). Then select cell I8 and drag the mouse to copy it to the cell in column I, row n (n=number of sampling points per vibration*total vibrations+8) to obtain the axial dynamic stress σ. d ;

[0120] Step 7. Insert the function =C8*10 / $F$4 in cell J8 (*10 here is for unit conversion, i.e. N / cm2 / 10=103N / m2=kPa). Then select cell J8 and drag the mouse to copy it to the cell in column J and row n (n=number of sampling points per vibration*total vibrations+8) to obtain the axial dynamic strain ε. d ;

[0121] Combine Figure 5 , perform the following steps in the "Data Processing" worksheet:

[0122] In the eighth step, after finding the minimum axial force of the first vibration in the fourth step and finding the corresponding sampling point number, output the entire row of data under the corresponding sampling point at intervals of n (the number of sampling points per vibration) to obtain the minimum axial force data group under each vibration and output it to the data processing part. Specifically, there are the following steps:

[0123] Insert the function =OriginalData!J4 in cell I4;

[0124] Insert the function =I4+original data!$H$3 in cell I5, then select cell I5 and drag the mouse to copy it to the cell in column I, row m (m = total vibrations + 4). Determine the number of sampling points corresponding to the minimum axial force at each vibration.

[0125] Step 9: After determining the number of sampling points corresponding to the minimum axial force at each vibration frequency in step 8, find and output the entire row of data corresponding to the specified number of sampling points. Specifically, there are the following steps:

[0126] Insert the function =VLOOKUP(I4,original data!$A:$G,2,FALSE) in cell C4, then select cell C4 and drag the mouse to copy it to the cell in column C, row m (m=vibration times + 4) to filter out the dynamic load of each vibration time.

[0127] Insert the function =VLOOKUP(original data!$H4,original data!$B:$G,2,FALSE) in cell D4, then select cell D4 and drag the mouse to copy it to the cell in row m of column D (m=vibration times+4) to filter out the dynamic load of each vibration time.

[0128] Insert the function =VLOOKUP(original data!$H4,original data!$B:$G,3,FALSE) in cell E4, then select cell E4 and drag the mouse to copy it to the cell in the mth row of column E (m=vibration times + 4), and filter out the dynamic axial deformation of each vibration time;

[0129] Insert the function =VLOOKUP(original data!$H4,original data!$B:$G,4,FALSE) in cell F4, then select cell F4 and drag the mouse to copy it to the cell in row m of column F (m=vibration times + 4), and filter out the confining pressure of each vibration time;

[0130] Insert the function =VLOOKUP(original data!$H4,original data!$B:$G,5,FALSE) in cell G4, then select cell G4 and drag the mouse to copy it to the cell in column G, row m (m=vibration times + 4), and filter out the dynamic volume change of each vibration time;

[0131] Insert the function =VLOOKUP(original data!$H4,original data!$B:$G,6,FALSE) in cell H4, then select cell H4 and drag the mouse to copy it to the cell in the Hth column and the mth row (m=vibration times + 4) to filter out the sampling time of each vibration time;

[0132] Combine Figure 6 , perform the following steps in the "Data Processing" worksheet:

[0133] Step 10. After filtering out the required data in step 9, use Excel's calculation function to calculate the axial dynamic stress, axial dynamic strain, and dynamic volume strain. Specifically, there are the following steps:

[0134] Insert the function =C4*10 / original data!$B$5 in cell K4 (*10 here is for unit conversion, i.e. N / cm2 / 10=103N / m2=kPa). Then select cell K4 and drag the mouse to copy it to the cell in column K and row m (m=number of vibrations + 4). Calculate the axial dynamic stress σ at each number of vibrations. d ;

[0135] Insert the function =D4*10 / original data!$F$4 in cell L4 (*10 here is for unit conversion, i.e. N / cm2 / 10=103N / m2=kPa). Then select cell L4 and drag the mouse to copy it to the cell in column L and row m (m=number of vibrations + 4). Calculate the axial dynamic strain ε at each vibration. d ;

[0136] Insert the function = 100 * G4 / original data! $D$5 in cell M4, then select cell M4 and drag the mouse to copy it to the cell in column M and row m (m = vibration + 4), and calculate the dynamic volume strain ε at each vibration. V ;

[0137] Step 11: Based on the dynamic strain and dynamic body strain calculated in step 10, use Excel's calculation function to further calculate the residual shear strain γ r and residual volume strain ε vr Specifically, there are the following steps:

[0138] Insert the function = (L4-$L$4)*(1+original data!$J$3) in cell N4, then select cell N4 and drag the mouse to copy it to the cell in the Nth column and the mth row (m = vibration times + 4), and calculate the residual shear strain γ at each vibration time. r ;

[0139] Insert the function =M4-$M$4 in cell O4, then select cell O4 and drag the mouse to copy it to the cell in column O and row m (m = vibration number + 4), and calculate the residual volume strain at each vibration number.

[0140] Combine Figure 7 , perform the following steps in the "Data Processing" worksheet:

[0141] The twelfth step is to process the existing data for subsequent model parameter calculation and to calculate γ r , ε vr And log(N+1) takes the base 10 logarithm. Specifically, there are the following steps:

[0142] Insert the function =LOG10(N5:N53) into cell Q5, then select cell Q5 and drag the mouse to copy it to the cell in the Qth column and the mth row (m = vibrations + 4) to calculate logγ r , and similarly find logε vr , log(log(N+1));

[0143] Step 13: Based on the model calculation results, output the γ predicted by the model at each vibration r and ε vr , specifically the following steps:

[0144] Insert the function = dynamic residual deformation calculation model! $C$8*A5^dynamic residual deformation calculation model! $C$9 in cell T5, then select cell T5 and drag the mouse to copy it to the cell in the Tth column and the mth row (m = vibration times + 4) to calculate γ r , we can find ε by the same logic vr ;

[0145] Combine Figure 8-10 , do the following steps in the Drawing worksheet:

[0146] Step 14: Insert an XY smooth curve chart in the "Chart Type" option bar of the Excel software. Select the data columns corresponding to cells H4 and K4 in the "Data Processing" worksheet, drag the mouse to the end of the column data, set the two columns of data as X and Y respectively, and draw the dynamic stress time history curve.

[0147] In the fifteenth step, modify the corresponding data of the Y series into the data columns corresponding to the "dynamic strain" and "pore water pressure" cells in the "data processing" worksheet, that is, the corresponding columns of cell L4 and cell M4, and draw the dynamic strain time history curve and pore water pressure time history curve in the same way as drawing the dynamic stress time history curve.

[0148] Step 16: Insert an XY smooth curve chart in the "Chart Type" option bar of the "Insert" tab of the Excel software. Select the data columns corresponding to cells K4 and L4 in the "Data Processing" worksheet, and drag the mouse to the end of the column data. Set the two columns of data as X and Y respectively, and draw the stress-strain curve.

[0149] Step 17. Insert an XY smooth curve chart in the "Chart Type" of the "Insert" option bar in Excel software, select the data columns corresponding to cells B4 and N4 in the "Data Processing" worksheet, use the mouse to drag to the end of the column data, set the two columns of data as X and Y respectively, and draw the residual shear strain versus vibration number curve; modify the corresponding data of the Y series to the data column corresponding to cell O4 in the "Data Processing" worksheet, and draw the residual volume strain versus vibration number curve in the same way as drawing the residual shear strain versus vibration number curve.

[0150] Combine Figure 11 , perform the following steps in the "Dynamic Residual Deformation Calculation Model" worksheet:

[0151] Step 18. Insert the function =10^INTERCEPT(data processing!Q5#,data processing!S5#) in cell C8 to obtain the model parameter a.

[0152] Step 19. Insert the function =SLOPE(data processing!Q5#, data processing!S5#) in cell C9 to obtain the model parameter b.

[0153] Step 20. Insert the function =CORREL(data processing!S5:S53, data processing!Q5:Q53) in cell C10 to obtain the coefficient of determination R1.

[0154] Similarly

[0155] Step 21. Insert the function =10^INTERCEPT(data processing!Q5#,data processing!S5#) in cell G8 to obtain the model parameter c.

[0156] Step 22. Insert the function =SLOPE(data processing!Q5#, data processing!S5#) in cell G9 to obtain the model parameter d.

[0157] Step 23. Insert the function =CORREL(data processing!S5:S53, data processing!Q5:Q53) in cell G10 to obtain the coefficient of determination R2.

[0158] Step 24: If the dynamic residual deformation of saturated soil needs to be calculated later, the basic data can be automatically calculated by simply entering it into the original data table. The results will be directly output in the "Dynamic Residual Deformation Calculation Model" worksheet; at the same time, relevant curve graphs will also be generated.

[0159] The implementation method uses the function calculation function of Excel software to automatically solve the dynamic residual shear strain and dynamic residual volume strain, calculate the calculation model parameters, and uses the drawing function to automatically draw the dynamic stress time history curve, dynamic strain time history curve, pore water pressure time history curve, stress strain curve, and the residual shear strain, residual volume strain and vibration frequency change relationship curve and model prediction curve. The operation is carried out on a visual interface, which is simple to operate, with accurate calculation results, automatic drawing and automatic fitting, and beautiful graphics.

[0160] The above are all preferred embodiments of the present application, and are not intended to limit the scope of protection of the present application. Therefore, any equivalent changes made based on the structure, shape, and principle of the present application should be included in the scope of protection of the present application.

Claims

1. A method for calculating dynamic residual deformation of saturated soil, characterized in that The steps include: S1. Establish a calculation model for dynamic residual deformation of saturated soil; S2. Obtain existing raw test data; S3, processing the test data obtained in S2 to obtain fitting data for training the model; S4, using the fitting data obtained in S3, the computational model in S1 is trained and fitted to obtain a trained computational model; S5, calculating the dynamic residual deformation of saturated soil according to the calculation model trained in S4; The calculation model of saturated soil dynamic residual deformation in step S1 is: c r =a·log b (N+1) e vr =c log d (N+1); where γ r is the residual shear strain, ε vr is the residual volume strain, a, b, c, d are the model fitting parameters to be determined, and N is the number of vibrations; In step S4, the data obtained in S3 are used to create X dataset, Y dataset and Z dataset. The data in X dataset is log(log(N+1)), and the data in Y dataset is logγ r , the data in the Z data set is logε vr ,but: Where x is the sample in the X data set, y is the sample in the Y data set, z is the sample in the Z data set, and n is the sample size. is the sample mean of the X data set, is the sample mean of the Y data set; is the sample mean of the Z data set; In step S4, R1 and R2 are the determination coefficients of the saturated soil dynamic residual deformation calculation model. When |R1| ≥ 0.9, the fitting of a and b is stopped, and when |R2| ≥ 0.9, the fitting of c and d is stopped, thereby obtaining a mature calculation model, where 2. A method for calculating dynamic residual deformation of saturated soil according to claim 1, characterized in that: The test data to be acquired in step S2 include: sample starting height h0, sample starting volume V0, sample starting area A0, sample height after consolidation h c , consolidation drainage ΔV, sample area after consolidation A c ; Among them, the height of the sample after consolidation h c Sample area after consolidation A c Sample volume after consolidation V c In c =h c And c 。 3. The method for calculating dynamic residual deformation of saturated soil according to claim 2, characterized in that: The fitting data in step S3 includes: Axial dynamic strain ε d Where Δh d is the axial dynamic deformation, in millimeters, h c is the height of the specimen after consolidation, in millimeters; Dynamic volume strain ε V Dynamic shear strain γ d c d =e d (1+m) Where μ is Poisson's ratio; Residual axial strain ε dr e dr =e di -e d0 Where ε di is the axial strain of the i-th vibration, ε d0 is the initial axial strain; Residual shear strain γ r c r =c i -γ0 Where, γ i is the shear strain of the ith vibration, γ0 is the initial shear strain; Residual volume strain ε vr e vr =e vi -e v0 Where, ε vi is the volume strain of the ith vibration, ε v0 is the initial volume strain.

4. A method for calculating the dynamic residual deformation of saturated soil, which is used to implement the method for calculating the dynamic residual deformation of saturated soil described in claim 3 in an Excel spreadsheet, characterized in that The steps include: S1. Create three worksheets in Excel: original data, data processing, and dynamic residual deformation calculation model; S2. Input the saturated soil dynamic residual deformation calculation model into the dynamic residual deformation calculation model worksheet; S3. Input the dynamic triaxial test data into an Excel raw data worksheet and filter the data to obtain the test data corresponding to the minimum axial force at each vibration frequency; And the parameters of the sample after consolidation, axial dynamic stress and axial dynamic strain are calculated; S4. Calculate the specimen parameters, axial dynamic stress, and axial dynamic strain after consolidation in the original data worksheet; S5. In the data processing worksheet, calculate the residual shear strain and the residual volume strain based on the test data selected in step S4; S6. In the dynamic residual deformation calculation model worksheet, use the regression analysis formula to calculate the parameters and determination coefficient R of the saturated soil dynamic residual deformation calculation model to obtain a mature calculation model; S7. Calculate the dynamic residual deformation of saturated soil using the mature calculation model obtained in S6.

Images