CN113125263B - Forecasting method for stress deformation of silica sol cured non-breakable sandy soil - Google Patents
Forecasting method for stress deformation of silica sol cured non-breakable sandy soil Download PDFInfo
- Publication number
- CN113125263B CN113125263B CN202110413274.3A CN202110413274A CN113125263B CN 113125263 B CN113125263 B CN 113125263B CN 202110413274 A CN202110413274 A CN 202110413274A CN 113125263 B CN113125263 B CN 113125263B
- Authority
- CN
- China
- Prior art keywords
- strain
- stress
- epsilon
- shear strain
- breakable
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/24—Investigating strength properties of solid materials by application of mechanical stress by applying steady shearing forces
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Analytical Chemistry (AREA)
- Pure & Applied Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention relates to a forecasting method for stress deformation of silica sol solidified non-breakable sand, which is characterized in that a model is established for a silica gel-sand complex formed after silica sol seeps non-breakable sand from the angle of particle friction and cementation failure, the volume deformation and the volume deformation rate in the stress process are more accurately forecasted, and the typical characteristics of the volume deformation rate can be forecasted.
Description
Technical Field
The invention belongs to the field of geotechnical engineering research, and particularly relates to a prediction method for stress deformation of silica sol-cured non-breakable sandy soil.
Background
The silica sol is formed by suspending nano-scale silica particles in water, has viscosity similar to that of water, and is quickly infiltrated through sandy soil, and then the silica particles are gathered into silica gel with a three-dimensional network structure, so that a solidified soil body is cemented. The sand on land is generally non-breakable quartz sand, and after the quartz sand is solidified by silica sol, cementing damage and failure can occur in the stress process. The conventional method for forecasting the deformation of the silica sol solidified quartz sand has a large error, so that a method is lacked, and the volume deformation rate of a silica gel-quartz sand complex formed after the silica sol seeps into the quartz sand can be accurately forecasted in the stress process.
Disclosure of Invention
The invention provides a forecasting method for the stress deformation of silica sol cured non-breakable sandy soil, which aims to accurately forecast the volume deformation and the volume deformation rate of a silica gel-quartz sand complex formed after silica sol seeps quartz sand in the stress process.
The invention relates to some abbreviations and symbols, the following are notes:
σ 1 : vertical stress to which the aggregate of particles is subjected;
σ 2 and σ 3 : horizontal stress, σ, to which the assembly of particles is subjected 2 And σ 3 The direction of (A) is vertical;
ε 1 、ε 2 and ε 3 : strain and respectively stress sigma 1 、σ 2 And σ 3 The directions are the same;
ε v : bulk strain, epsilon v =ε 1 +ε 2 +ε 3 ;
dε v : a bulk strain increment;
dε s : a shear strain increment;
d: for the solidified soil sample, the work generated by friction and silica gel cementation failure when the solidified soil sample is stressed;
E B : energy given off by failure of the silica gel bond;
m: material parameters and equal to critical stress ratio;
beta: is equal to the average effective stress p 0 Parameter of interest, β ═ k β1 p 0 +k β2 ;
k β1 ,k β2 : is a material parameter;
ρ 2 :ρ 2 (ε s ) Is a variable related to shear strain;
τ 1 ,τ 2 : d epsilon in Legendre transform v And d ε s The corresponding conjugate variable is τ 1 ,τ 2 And p ═ τ 1 ,q-ρ 2 (ε s )=τ 2 ;
k B1 and k B2 : parameter, and B 4 =-k B1 p 0 +k B2 。
The technical scheme of the invention is as follows: a forecasting method for stress deformation of silica sol cured non-breakable sandy soil comprises the following steps:
step 1: setting the vertical stress of the silica sol cured non-breakable sandy soil to be sigma 1 Stress on the horizontal plane is respectively sigma 2 And σ 3 Where σ is 2 And σ 3 Is perpendicular to the direction of the particle assembly, the strain of the particle assembly is epsilon 1 、ε 2 And ε 3 Wherein strain ε 1 、ε 2 And ε 3 Respectively in the direction of the stress sigma 1 、σ 2 And σ 3 The directions are the same; defining average effective stress p, shear stress q, stress ratio eta, and bulk strain epsilon v And shear strain epsilon s :
ε v =ε 1 +ε 2 +ε 3
And 2, step: let the work D generated by friction and silica gel bond failure be:
m is a material parameter and is equal to a critical stress ratio, beta is the initial average effective stress p 0 Parameter of interest, p 2 (ε s ) D ε is a variable related to shear strain v D ε being the increase in bulk strain s For shear strain increments, E B The energy given off for the failure of the silicone gel bond,is E B For epsilon s Partial derivatives of (d);
taking a parameter alpha:
let d ε in Legendre transform v And d ε s The corresponding conjugate variable is τ 1 ,τ 2 ,τ 1 ,τ 2 And p, q are related as follows:
legendre transformation is carried out on the formula (1), and the formulas (2) and (3) are substituted to obtain the following equation:
and step 3: taking each shear strain increment, the ratio of the bulk strain to the shear strain increment is:
and 4, step 4: calculating the current shear strain ε s To the next shear strain epsilon s +dε s Delta d epsilon of body strain v :
And 5: calculating shear strain ε s +dε s Body strain at time: d epsilon obtained by calculation in the step 4 v With current shear strain epsilon s Adding the corresponding bulk strains to obtain epsilon s +dε s Time bulk strain epsilon v Namely: epsilon v ←ε v +dε v ;
Step 6: calculating the shear strain ε by equation (4) s +dε s Corresponding toWill epsilon s +dε s Calculated by substituting formula (4)
And 7: repeating the step 4 to the step 6 to obtain the body strain sum corresponding to each shear strain in the loading process
in the above formula B y1 ,B y2 And B y3 Are parameters. Get B y1 When the number is equal to 1, the alloy is put into a container,in relation to shear strainThe actual independent parameter is B y2 ,B y3 。
Preferably, in the formula (5) (. rho) 2 The relationship to shear strain is:
in the above formula B 1 ,B 2 ,B 3 ,B 4 And B 5 Are parameters.
in the above formula, the independent parameter is B 1 ,B 4 ,B 5 。
Preferably, β is the sum of the mean effective stress p 0 Parameter of interest, β ═ k β1 p 0 +k β2
k β1 And k β2 Is a material parameter.
Preferably, B 4 Is and initial mean stress p 0 A parameter of interest, and B 4 =-k B1 p 0 +k B2 And k is B1 And k B2 Is a material parameter.
When the horizontal effective stress sigma 2 And σ 3 When the stress is constant and equal, the vertical effective stress sigma can be obtained 1 And average effective stress and shearStress:
preferably, in step 7, if the stress-strain relationship is obtained for the triaxial test, the stress-strain relationship is obtained byCalculating the shear stress q as:
the invention has the beneficial effects that: the volume deformation and the volume deformation rate of a silica gel-quartz sand complex formed after silica sol seeps quartz sand are more accurately forecasted in the stress process.
Drawings
FIG. 1 is a schematic diagram of a silica gel-silica sand composite formed after silica sol has percolated through the silica sand, and subjected to vertical and horizontal stresses;
FIG. 2 is a graph showing the variation of shear stress with shear strain and the variation of bulk strain with shear strain;
FIG. 3 is a graph of the ratio of bulk strain increase to shear strain increase versus stress.
1 is a silica gel-quartz sand composite.
Detailed Description
In order to make the technical means, innovative features, objectives and effects of the present invention apparent, the present invention will be further described with reference to the following detailed drawings.
The invention relates to some abbreviations and symbols, the following are notes:
σ 1 : vertical stress to which the assembly of particles is subjected
σ 2 And σ 3 : particle assembly receiverTo horizontal stress, σ 2 And σ 3 Is directed perpendicularly to
ε 1 、ε 2 And epsilon 3 : strain and respectively stress sigma 1 、σ 2 And σ 3 Same direction
p 0 : initial mean effective stress
ε v : bulk strain, epsilon v =ε 1 +ε 2 +ε 3
dε v : increase in bulk strain
dε s : increment of shear strain
D: work produced by friction and silica gel cementation failure when stressed on a solidified soil sample
E B : energy released by failure of silica gel cementation
M: material parameter equal to critical stress ratio
Beta: is equal to the initial mean effective stress p 0 Parameter of interest, β ═ k β1 p 0 +k β2
k β1 ,k β2 : as a parameter of the material
ρ 2 :ρ 2 (ε s ) As a variable related to shear strain
τ 1 ,τ 2 : d epsilon in Legendre transform v And d ε s The corresponding conjugate variable is τ 1 ,τ 2 ,
And p ═ τ 1 ,q-ρ 2 (ε s )=τ 2
B 1 ,B 2 ,B 3 ,B 4 And B 5 :ρ 2 A parameter related to the shear strain of the steel,
k B1 And k B2 : parameter, and B 4 =-k B1 p 0 +k B2
Example 1
Performing seepage curing on Fujian quantan standard sand by using 30% silica sol, wherein the sand is quartz sand, and performing drainage triaxial tests with confining pressure of 200kPa and 300kPa, wherein the stress on a horizontal plane is sigma respectively 2 And σ 3 And σ is 2 And σ 3 Equal to the confining pressure, volume deformation and volume deformation rate are predicted based on the method proposed herein.
The technical scheme of the invention is as follows: a forecasting method for stress deformation of silica sol cured non-breakable sandy soil comprises the following steps:
step 1: setting the vertical stress of the silica sol cured non-breakable sandy soil as sigma shown in figure 1 1 Stress on the horizontal plane is respectively σ 2 And σ 3 Where σ is 2 And σ 3 Is perpendicular to the direction of the particle assembly and the strain of the particle assembly is epsilon 1 、ε 2 And ε 3 Wherein strain ε 1 、ε 2 And ε 3 Respectively in the direction of the stress sigma 1 、σ 2 And σ 3 The directions are the same; defining average effective stress p, shear stress q, stress ratio eta, and bulk strain epsilon v And shear strain epsilon s :
ε v =ε 1 +ε 2 +ε 3
And 2, step: let the work D generated by friction and silica gel bond failure be:
m is a material parameter and is equal to a critical stress ratio, beta is the initial average effective stress p 0 Parameter of interest, p 2 (ε s ) D ε is a variable related to shear strain v D ε being the increase in bulk strain s For shear strain increments, E B The energy given off for the failure of the silicone gel bond,is E B For epsilon s Partial derivatives of (d);
taking a parameter alpha:
let d ε in Legendre transform v And d ε s The corresponding conjugate variable is τ 1 ,τ 2 ,τ 1 ,τ 2 And p, q are related as follows:
legendre transformation is carried out on the formula (1), and the formulas (2) and (3) are substituted to obtain the following equation:
and 3, step 3: taking the ratio of the body strain to the shear strain increment (i.e. the body strain rate or volume deformation rate) in each shear strain increment as:
and 4, step 4: calculating the current shear strain ∈ s To the next shear strain epsilon s +dε s Delta d epsilon of body strain v :
Here, theIn relation to shear strain ofIn the formula B y1 ,B y2 And B y3 Taking B as a parameter y1 When the number is equal to 1, the alloy is put into a container,in relation to shear strainActual independent parameter is B y2 ,B y3 ;
Here ρ 2 Is composed ofIn the formula B 1 ,B 2 ,B 3 ,B 4 And B 5 Is a parameter; get B 2 =2,The actual independent parameter is B 1 ,B 4 ,B 5 (ii) a And B 4 Is the sum of the initial mean stress p 0 Parameter of interest, B 4 =-k B1 p 0 +k B2 And k is B1 And k B2 Is a material parameter;
and 5: calculating shear strain epsilon s +dε s Body strain at time:
d epsilon obtained by calculation in the step 4 v With current shear strain epsilon s Adding the corresponding bulk strains to obtain epsilon s +dε s Time bulk strain epsilon v Namely: epsilon v ←ε v +dε v ;
Will epsilon s +dε s Corresponding to formula (4) calculationIn the formula (4), beta is the sum of the mean effective stress p 0 Parameter of interest, β ═ k β1 p 0 +k β2 ,k β1 And k β2 Is a material parameter; and calculating the shear stress from the triaxial test conditions
And 7: repeating the step 4 to the step 6 to obtain the body strain corresponding to each shear strain in the loading process,And shear stress q.
The following table lists the values of the above calculated parameters.
Table 1 independent parameter values in the calculation process
M | k β1 | k β2 | y B2 | y B3 | B 1 | k B1 | k B2 | B 5 |
1.4 | -0.251×10 -3 | 175.8 | 9 | 1.1 | 15 | -0.006×10 -3 | 3 | 0.023 |
FIG. 2 shows the curve of shear stress with shear strain and the curve of body strain with shear strain, which shows that the body strain forecast is more in line with the test condition; FIG. 3 shows the ratio of the increase in bulk strain to the increase in shear strain (i.e., the bulk strain rate or referred to as the bulk strain rate) versus stress curve whereThe bending section obtained by the test is forecasted before the peak valueAfter the peak value, a downward hook returning curve obtained by the test is well forecasted.
Claims (7)
1. A forecasting method for stress deformation of silica sol solidified non-breakable sandy soil is characterized in that: the method specifically comprises the following steps:
step 1: setting the vertical stress of the silica sol solidified non-breakable sandy soil to be sigma 1 Stress on the horizontal plane is respectively sigma 2 And σ 3 Where σ is 2 And σ 3 Is perpendicular to the direction of the particle assembly and the strain of the particle assembly is epsilon 1 、ε 2 And ε 3 In which strain epsilon 1 、ε 2 And ε 3 Respectively in the direction of the stress sigma 1 、σ 2 And σ 3 The directions are the same; defining average effective stress p, shear stress q, stress ratio eta, and bulk strain epsilon v And shear strain epsilon s :
ε v =ε 1 +ε 2 +ε 3
Step 2: let the work D generated by friction and silica gel bond failure be:
m is a material parameter and is equal to a critical stress ratio, beta is the initial average effective stress p 0 Parameter of interest, p 2 (ε s ) D ε is a variable related to shear strain v D ε is the increase in bulk strain s For shear strain increments, E B The energy given off for the failure of the silicone gel bond,is E B For epsilon s Partial derivatives of (d);
taking a parameter alpha:
let d ε in Legendre transform v And d ε s The corresponding conjugate variable is τ 1 ,τ 2 ,τ 1 ,τ 2 And p, q are related as follows:
legendre transformation is carried out on the formula (1), and the formulas (2) and (3) are substituted to obtain the following equation:
and step 3: taking each shear strain increment, the ratio of the bulk strain to the shear strain increment is:
and 4, step 4: calculating the current shear strain ∈ s To the next shear strain epsilon s +dε s Delta d epsilon of body strain v :
And 5: calculating shear strain ε s +dε s Body strain at time: d epsilon obtained by calculation in the step 4 v With current shear strain epsilon s Adding the corresponding bulk strains to obtain epsilon s +dε s Time bulk strain epsilon v Namely: epsilon v ←ε v +dε v ;
Step 6: calculating the shear strain ε by equation (4) s +dε s Corresponding toWill epsilon s +dε s Calculated by substituting formula (4)
4. The method for forecasting the stress deformation of the silica sol-cured non-breakable sandy soil according to claim 1, wherein the method comprises the following steps: in the formula (5) (. rho) 2 The relationship to shear strain is:
in the above formula B 1 ,B 2 ,B 3 ,B 4 And B 5 Are parameters.
6. The method for forecasting the stress deformation of the silica sol-cured non-breakable sandy soil according to claim 1, wherein the method comprises the following steps: beta is the sum of the initial mean effective stress p 0 Parameter of interest, β ═ k β1 p 0 +k β2 ,k β1 And k β2 Is a material parameter.
7. The method for forecasting the stress deformation of the silica sol-cured non-breakable sandy soil according to claim 1, wherein the method comprises the following steps: b is 4 Is the sum of the initial mean effective stress p 0 A parameter of interest, and B 4 =-k B1 p 0 +k B2 And k is B1 And k B2 Is a material parameter.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110413274.3A CN113125263B (en) | 2021-04-16 | 2021-04-16 | Forecasting method for stress deformation of silica sol cured non-breakable sandy soil |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110413274.3A CN113125263B (en) | 2021-04-16 | 2021-04-16 | Forecasting method for stress deformation of silica sol cured non-breakable sandy soil |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113125263A CN113125263A (en) | 2021-07-16 |
CN113125263B true CN113125263B (en) | 2022-09-09 |
Family
ID=76776925
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110413274.3A Active CN113125263B (en) | 2021-04-16 | 2021-04-16 | Forecasting method for stress deformation of silica sol cured non-breakable sandy soil |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113125263B (en) |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4854175A (en) * | 1988-02-29 | 1989-08-08 | The Research Foundation Of State University Of New York | Simple shear device for testing earthen materials and powders |
US8265915B2 (en) * | 2007-08-24 | 2012-09-11 | Exxonmobil Upstream Research Company | Method for predicting well reliability by computer simulation |
CN110411804B (en) * | 2019-09-02 | 2024-04-16 | 上海交通大学 | Test sample for mechanical properties of contact surface of soil body and structure, preparation method and test method |
CN111024486B (en) * | 2019-12-19 | 2020-07-24 | 南京航空航天大学 | Creep behavior prediction method for unidirectional ceramic matrix composite |
-
2021
- 2021-04-16 CN CN202110413274.3A patent/CN113125263B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN113125263A (en) | 2021-07-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110261573B (en) | Dynamic evaluation method for stability of high-position rocky landslide | |
Damgaard et al. | Dynamic response sensitivity of an offshore wind turbine for varying subsoil conditions | |
CN106568660A (en) | Method for predicting residual fatigue life of composite material adhesive bonding repair structure | |
CN106021709A (en) | Early concrete cracking risk assessment and control method | |
CN112461657B (en) | Rapid prediction method for critical damage stress of roadbed soil | |
CN113125263B (en) | Forecasting method for stress deformation of silica sol cured non-breakable sandy soil | |
CN110610041B (en) | Method for judging limit strain of instability and damage of shaft | |
CN107066743A (en) | A kind of extracting method of the Shearer Helical Drum load based on distinct element method | |
CN111563343B (en) | Method for determining elastic modulus of rock-fill concrete | |
CN104156590B (en) | A kind of method for building up of Mg-based nanocomposite thixotroping Plastic Forming constitutive model | |
CN114722681A (en) | Simulation prediction method for ground settlement caused by shield construction | |
CN112651153A (en) | Method for determining material parameters of crystal plastic finite element model | |
CN106951661A (en) | The separation computational methods of cemented gravel dam strain gauge actual measurement strain | |
CN113155612B (en) | Deformation prediction method for microfiber mixed silica sol solidified calcareous sand | |
CN111859570A (en) | Dynamic reliability assessment method for bridge crane structure | |
CN107862142B (en) | Method for analyzing mechanical strength of slotted casing | |
CN113125262B (en) | Method for quickly forecasting deformation of breakable calcareous sand in loading process | |
Xu et al. | Time-varying reliability analysis based on improved toughness exhaustion model and probability density evolution method to predict fatigue damage life | |
CN111398019B (en) | Method for rapidly judging relative size of rock damage under loads with different strain rates | |
Zhang et al. | Stability analysis of channel slope based on FEM strength reduction | |
CN117556644B (en) | Forward slope quantitative discrimination method for steep forward layer rock high slope | |
CN113392558B (en) | Rock damage determination method based on maximum shear strain and main tensile strain | |
CN115455685B (en) | Method for calculating buffeting response of long cable structure under two-dimensional turbulent wind excitation | |
Chen et al. | Feasibility study of a steel-UHPFRC hybrid tower for offshore wind turbines | |
CN117195507B (en) | Simulation prediction method and system for landslide river-plugging impact dam area building process |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |