CN114818180B - Construction method of sea water aging immersed tube tunnel GINA water stop time-varying constitutive model - Google Patents
Construction method of sea water aging immersed tube tunnel GINA water stop time-varying constitutive model Download PDFInfo
- Publication number
- CN114818180B CN114818180B CN202210410929.6A CN202210410929A CN114818180B CN 114818180 B CN114818180 B CN 114818180B CN 202210410929 A CN202210410929 A CN 202210410929A CN 114818180 B CN114818180 B CN 114818180B
- Authority
- CN
- China
- Prior art keywords
- water stop
- aging
- gina
- gina water
- stress
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 141
- 230000032683 aging Effects 0.000 title claims abstract description 124
- 239000013535 sea water Substances 0.000 title claims abstract description 47
- 238000010276 construction Methods 0.000 title description 7
- 230000035882 stress Effects 0.000 claims abstract description 53
- 230000008859 change Effects 0.000 claims abstract description 20
- 238000000034 method Methods 0.000 claims abstract description 18
- 238000012360 testing method Methods 0.000 claims description 41
- 238000009864 tensile test Methods 0.000 claims description 17
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 10
- 239000007788 liquid Substances 0.000 claims description 10
- 230000006835 compression Effects 0.000 claims description 9
- 238000007906 compression Methods 0.000 claims description 9
- 238000003825 pressing Methods 0.000 claims description 8
- 239000004743 Polypropylene Substances 0.000 claims description 6
- 239000007789 gas Substances 0.000 claims description 6
- 239000004033 plastic Substances 0.000 claims description 6
- 229920001155 polypropylene Polymers 0.000 claims description 6
- 238000001816 cooling Methods 0.000 claims description 5
- 229910052757 nitrogen Inorganic materials 0.000 claims description 5
- -1 polypropylene Polymers 0.000 claims description 5
- 239000000314 lubricant Substances 0.000 claims description 4
- 230000036961 partial effect Effects 0.000 claims description 4
- 230000000452 restraining effect Effects 0.000 claims description 4
- 230000003679 aging effect Effects 0.000 claims description 3
- 239000011248 coating agent Substances 0.000 claims description 3
- 238000000576 coating method Methods 0.000 claims description 3
- 238000013101 initial test Methods 0.000 claims description 3
- 238000007789 sealing Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 abstract description 4
- 238000011160 research Methods 0.000 abstract description 3
- 239000000463 material Substances 0.000 description 5
- 230000015556 catabolic process Effects 0.000 description 3
- 238000006731 degradation reaction Methods 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 239000003513 alkali Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006065 biodegradation reaction Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000012669 compression test Methods 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 238000004132 cross linking Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000002427 irreversible effect Effects 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000005297 material degradation process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003647 oxidation Effects 0.000 description 1
- 238000007254 oxidation reaction Methods 0.000 description 1
- 229920000642 polymer Polymers 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000002829 reductive effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000008961 swelling Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
- G01N1/286—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q involving mechanical work, e.g. chopping, disintegrating, compacting, homogenising
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
- G01N1/42—Low-temperature sample treatment, e.g. cryofixation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N17/00—Investigating resistance of materials to the weather, to corrosion, or to light
-
- 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/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
- G01N3/10—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces generated by pneumatic or hydraulic pressure
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
- G01N1/286—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q involving mechanical work, e.g. chopping, disintegrating, compacting, homogenising
- G01N2001/2873—Cutting or cleaving
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/04—Ageing analysis or optimisation against ageing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Immunology (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Civil Engineering (AREA)
- Architecture (AREA)
- Structural Engineering (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Environmental & Geological Engineering (AREA)
- Environmental Sciences (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention relates to a method for constructing a time-varying constitutive model of a immersed tube tunnel GINA water stop with seawater ageing, which belongs to the technical field of immersed tube waterproof equipment research, and comprises the steps of obtaining stress relaxation curves at different ageing temperatures; determining an aging performance change value P of rubber used for the GINA water stop; determining a stress-strain relation curve of the GINA water stop; obtaining a stress-strain relation curve of a full aging period; and constructing a constitutive model of stress relaxation and seawater aging of the GINA water stop belt. The invention can realize the dynamic monitoring of the service state of the GINA water stop and provides a basis for the service life assessment and risk early warning of the GINA water stop.
Description
Technical Field
The invention belongs to the technical field of immersed tube waterproof equipment research, and particularly relates to a construction method of a marine aging immersed tube tunnel GINA water stop variable structure model.
Background
Along with the continuous improvement of the economic level of China, particularly the construction of immersed tube tunnels such as a bridge of a port Zhuhai and Macao, a deep and medium channel and the like, the sign that China has become a large country for constructing the immersed tube tunnels.
The pipe joint connector is a weak link for water prevention of the immersed tunnel, and the GINA water stop is used as a main water stop unit of the pipe joint connector, so that the service performance and the service life of the immersed tunnel are directly affected. The GINA water stop is used only in the 60 th century, and the GINA water stop which is established in China at present and used in established immersed tube tunnels is all from the Netherlands TRELLEBORG, VREDESTEIN or Japanese Kogyo rubber Co. At present, no engineering data show that the GINA water stop belt can be in healthy service for 100-120 years without problems. On one hand, the GINA water stop belt is directly contacted with sea water and is subjected to the effects of oxidation, swelling, saline-alkali corrosion, mechanical force, biodegradation and the like, and the high molecular polymer of the GINA water stop belt can generate irreversible physical and chemical reactions such as crosslinking, degradation, cracking and the like, so that the performance of the GINA water stop belt is reduced; on the other hand, under the load effect, the GINA water stop belt can produce stress relaxation, service performance degradation and other mechanical performances. And GINA water stops cannot be repaired and replaced during operation, and once large-scale leakage water is generated, immeasurable consequences and losses can occur. Therefore, the dynamic monitoring, evaluation and early warning of the service state of the GINA water stop and the prediction of the service life are carried out, the aging constitutive relation of the aging of the GINA water stop under the coupling effect of the environment and the load is determined in advance, and the method has great research significance.
Disclosure of Invention
The application provides a construction method of a sea water aging immersed tube tunnel GINA waterstop time-varying constitutive model based on a rubber constitutive relation of Mooney-Rivlin model representation GINA waterstop stress relaxation and sea water aging characteristics.
Firstly, carrying out a compression test on a selected GINA water stop test piece, when the compression amount reaches a designed compression amount, placing the test piece in a designed seawater aging bin for carrying out a stress relaxation test, and observing the material degradation and the mechanical property degradation of the GINA water stop under the combined action of seawater and load; cooling an unaged GINA water stop belt in liquid nitrogen, cutting the water stop belt into a dumbbell shape, carrying out a uniaxial tension test, and testing the stress-strain relation; and simulating a uniaxial tensile test process by using a Mooney-Rivlin model to construct a constitutive model of stress relaxation and seawater aging of the GINA water stop belt.
The application adopts the specific technical scheme that:
The method for constructing the time-varying constitutive model of the immersed tube tunnel GINA water stop with seawater ageing comprises the following steps:
1) Selecting a GINA water stop test piece, and performing seawater accelerated aging tests at different temperatures under the designed compression quantity on the GINA water stop test piece to obtain stress relaxation curves at different aging temperatures;
2) Determining the contact stress sigma and the initial contact stress sigma 0 of the GINA water stop after aging according to the stress relaxation curve in the step 1), and determining the aging coefficient k of the GINA water stop according to k=sigma/sigma 0; determining an aging performance change value P of rubber used for the GINA water stop according to the P=sigma-sigma 0;
3) Obtaining a lnP curve changing along with time t at normal temperature according to a time-temperature superposition principle, and obtaining an equation P=exp (f (t)), wherein f (t) is a change function of an aging performance change value P and time t; determining P= (1-k) sigma 0 according to the relation between the aging coefficient k of the GINA water stop and the aging property change value P of rubber used by the GINA water stop, and obtaining the normal temperature aging coefficient k Often times ＝1-exp(f(t))/σ0 of the GINA water stop;
4) Carrying out a uniaxial tensile test on the GINA water stop at the environment temperature of 23 ℃ and the tensile speed of 500mm/min to determine a stress-strain relation curve of the GINA water stop;
5) Correcting the stress-strain relation curve in the step (4) by using a normal temperature ageing coefficient k Often times of the GINA water stop belt to obtain a stress-strain relation curve of a full ageing period;
6) According to the stress-strain relation curve of the full aging period in the step 5), a constitutive model of the stress relaxation and seawater aging of the GINA water stop is built by combining a Mooney-Rivlin model:
σ＝2(λ2-λ-1)(f(t)+g(t)λ-1)；
Wherein t is aging time; lambda is the elongation ratio of the GINA water stop; f (t) is a function of C 10 over time t; g (t) is a function of the change of C 01 over time t, C 10 and C 01 are the Rivlin coefficients of the Mooney-Rivlin model, the values of which are determined by the stress-strain relationship curve for the full aging cycle.
Further defined, the step (6) specifically comprises:
6.1 Determining epsilon i according to the uniaxial tensile test of the GINA water stop in the step 4), determining 3 elongation ratios lambda i of the GINA water stop according to epsilon i, respectively marking as lambda 1、λ2 and lambda 3, and obtaining Green strain tensor invariants I 1 and I 2, wherein the following formula is shown in the specification:
λi＝1+εi (1)
I1＝λ1 2+λ2 2+λ3 2 (2)
I2＝λ1 2λ2 2+λ2 2λ3 2+λ1 2λ3 2 (3)
Wherein I 1 is a first strain tensor invariant and I 2 is a second strain tensor invariant; lambda 1、λ2 and lambda 3 represent the elongation ratio in the X-axis direction, the elongation ratio in the Y-axis direction, and the elongation ratio in the Z-axis direction, respectively; epsilon i is the strain, i=1, 2 or 3, wherein epsilon 1 is the strain in the X-axis direction, epsilon 2 is the strain in the Y-axis direction, and epsilon 3 is the strain in the Z-axis direction;
6.2 The strain energy density W derived based on the Mooney-Rivlin model is:
W＝C10(I1-3)+C01(I2-3) (4)
Wherein, C 10 and C 01 are Rivlin coefficients of a Mooney-Rivlin model, and the values of the Rivlin coefficients are determined by a stress-strain relation curve of a full aging period;
6.3 Through Kirchhoff stress versus Green strain): the relation between the contact stress sigma and the elongation ratio lambda after the aging of the GINA water stop is obtained as follows:
The partial derivative of I 1、I2 is calculated by the formula (4) to obtain Relationship between the contact stress sigma elongation ratio lambda after aging of GINA water stop:
σ＝2(λ2-λ-1)(C10+C01λ-1) (6)
6.4 Carrying out parameter identification on the formula (6) by a nonlinear least square method by adopting Origin software based on a stress-strain relation curve of a full aging period to obtain values C 10 and C 01 under the full life period t i, obtaining functions f (t) and g (t) of changes of C 10 and C 01 along with aging time according to curves of changes of C 10 and C 01 along with aging time, and further constructing a constitutive model of stress relaxation and seawater aging of the GINA water stop belt:
σ＝2(λ2-λ-1)(f(t)+g(t)λ-1) (7)。
Further defined, in the step (1), the conditions of the seawater accelerated aging test are:
placing a water stop with the Shore hardness of 50HS on a bottom plate, restraining the water stop by using a ballast plate and a layering, loading the water stop by using a press machine with the rated load of 3000kN according to the TB/T2843-2010 standard, controlling the ballast plate by using a middle plate after the compression amount reaches 125mm, placing a liquid gas pressure sensor on the middle plate, pressing the middle plate by using a top plate, fixing the top plate, the middle plate and the bottom plate, integrally placing the fixed top plate, the middle plate and the bottom plate in a polypropylene plastic sealing bin 1, filling 2/3 of natural seawater, and respectively collecting the contact stress of the GINA water stop in the seawater in real time within the range of 50-80 ℃.
Further defined, the uniaxial tensile test of the step (4) specifically comprises the following conditions:
Taking out the GINA water stop tape aged by the seawater in the step (1), cooling by liquid nitrogen, preparing a dumbbell-shaped test piece according to GB/T528-2009, coating a lubricant on the test piece, and carrying out a tensile test by a universal tester at a speed of 500mm/min at the room temperature of 23 ℃.
Further defined, the dumbbell test pieces had an overall length of 100mm, a thickness of 2.0mm, and an initial test length of 20.0mm for the test sections.
The application has the beneficial effects that:
1) The application is based on a Mooney-Rivlin model, comprehensively considers the influencing factors of sea water temperature and aging time, builds the constitutive model of the stress relaxation and sea water aging of the GINA water stop, can truly simulate, can realize dynamic monitoring of the service state of the GINA water stop, and provides a basis for service life assessment and risk early warning of the GINA water stop.
2) As can be seen from the investigation, at present, no material constitutive model specially reflecting the service state of the GINA water stop belt exists, and the conventional rubber constitutive model is mostly a static constitutive model, so that the time-varying characteristic of the stress-strain change of the GINA water stop belt material along with the time cannot be reflected. The technical advantages of the constitutive model are as follows: on one hand, the stress relaxation test and the uniaxial tensile test are carried out on the GINA water stop under the combined action of seawater and load, and then the Mooney-Rivlin model is applied to simulate according to the obtained test data, so that a time-varying constitutive model of the stress relaxation and the seawater aging of the GINA water stop is constructed, and the time-varying constitutive model is constructed according to the working environment and the working condition of the actual GINA water stop, so that more accurate and more practical life prediction data of the GINA water stop can be obtained; on the other hand, the construction model is simple in construction operation, clear in flow, high in accuracy and strong in feasibility, can represent the change of the constitutive relation of the GINA water stop belt along with time, can realize dynamic monitoring of the service state of the GINA water stop belt, and provides a good basis for service life assessment of the GINA water stop belt.
Drawings
FIG. 1 is a plan view of a seawater ageing test cartridge;
FIG. 2 is a perspective view of a sea water aging test chamber;
FIG. 3 is a GINA water stop test piece of different hardness;
FIG. 4 is a GINA waterstop stress relaxation curve;
FIG. 5 is a graph showing the aging coefficient k versus aging time;
FIG. 6 is a graph of stress-strain relationship;
FIG. 7 is a full aging cycle stress-strain curve;
FIG. 8 is a graph of C 10 versus aging time, and b is a graph of C 01 versus aging time;
Wherein, 1-polypropylene plastic seals the storehouse, 2-liquid gas pressure sensor, 3-GINA waterstop.
Detailed Description
The present application will be described in further detail with reference to the drawings and examples, in order to make the objects, technical solutions and advantages of the present application more apparent. It should be understood that the particular embodiments described herein are illustrative only and are not intended to limit the application, i.e., the embodiments described are merely some, but not all, of the embodiments of the application.
Accordingly, the following detailed description of the embodiments of the application is not intended to limit the scope of the application, as claimed, but is merely representative of selected embodiments of the application. All other embodiments, which can be made by a person skilled in the art without making any inventive effort, are intended to be within the scope of the present application.
It should be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a method that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such method or conventional elements.
The technical scheme of the application is further described with reference to the accompanying drawings and the embodiments.
The application provides a method for constructing a time-varying constitutive model of a immersed tube tunnel GINA water stop with seawater ageing, which comprises the following steps:
1) Selecting a GINA water stop test piece, and carrying out seawater accelerated aging tests of different temperatures under the designed compression quantity on the GINA water stop test piece, wherein the conditions of the seawater accelerated aging tests are as follows: placing a water stop with the Shore hardness of 50HS on a bottom plate, restraining the water stop by using a ballast plate and a pressing strip, loading the water stop by using a press machine with the rated load of 3000kN according to the TB/T2843-2010 standard, controlling the ballast plate by using a middle plate after the compression amount reaches 125mm, placing a liquid gas pressure sensor 2 on the middle plate, pressing the middle plate by using a top plate, fixing the top plate, the middle plate and the bottom plate, integrally placing the fixed top plate, the middle plate and the bottom plate in a polypropylene plastic sealed bin 1, filling 2/3 of natural seawater, and respectively collecting the contact stress of the GINA water stop 3 in the seawater in real time within the range of 50-80 ℃ to obtain stress relaxation curves at different aging temperatures.
2) Determining the contact stress sigma and the initial contact stress sigma 0 of the GINA water stop after aging according to the stress relaxation curve in the step 1), and determining the aging coefficient k of the GINA water stop 3 according to k=sigma/sigma 0; determining an aging performance change value P of rubber used by the GINA water stop 3 according to the P=sigma-sigma 0;
3) Obtaining a curve of lnP changing along with time t at normal temperature according to a time-temperature superposition principle, obtaining an equation P=exp (f (t)), and determining P= (1-k) sigma 0 according to the correlation between an aging coefficient k of the GINA water stop 3 and an aging property change value P of rubber used by the GINA water stop 3, thereby obtaining a normal temperature aging coefficient k Often times ＝1-exp(f(t))/σ0 of the GINA water stop 3, wherein f (t) represents a change function between P along with the change function;
4) Carrying out uniaxial tensile test on the GINA water stop 3 under the conditions that the ambient temperature is 23 ℃ and the tensile speed is 500mm/min, namely taking out the GINA water stop 3 after seawater aging, carrying out liquid nitrogen cooling, preparing a dumbbell-shaped test piece according to GB/T528-2009, coating a lubricant on the test piece, and carrying out tensile test at the speed of 500mm/min under the condition that the room temperature is 23 ℃ by using a universal tester; determining a stress-strain relation curve of the GINA water stop;
5) Correcting the stress-strain relation curve in the step (4) by using a normal temperature ageing coefficient k Often times of the GINA water stop belt to obtain a stress-strain relation curve of a full ageing period;
6) According to the stress-strain relation curve of the full aging period in the step 5), a constitutive model of the stress relaxation and seawater aging of the GINA water stop is built by combining a Mooney-Rivlin model:
σ＝2(λ2-λ-1)(f(t)+g(t)λ-1)
Wherein t is aging time; lambda is the elongation ratio of the GINA water stop; f (t) is a function of C 10 over time t; g (t) is a function of the change of C 01 with time t, C 10 and C 01 are the Rivlin coefficients of the Mooney-Rivlin model, and the values of the Rivlin coefficients are determined by a stress-strain relation curve of a full aging period;
The method comprises the following steps:
6.1 Determining epsilon i according to the uniaxial tensile test of the GINA water stop 3 in the step 4), determining elongation ratio lambda of the GINA water stop 3 directions according to epsilon i, and marking lambda 1、λ2 and lambda 3,εi as strains, i=1, 2 or 3, wherein epsilon 1 is the strain of the X axis direction, epsilon 2 is the strain of the Y axis direction and epsilon 3 is the strain of the Z axis direction; wherein λ 1 is the elongation ratio in the x direction, λ 2 is the elongation ratio in the y direction, and λ 3 is the elongation ratio in the z direction; green strain invariants I 1 and I 2 were determined as follows:
λi＝1+εi (1)
I1＝λ1 2+λ2 2+λ3 2 (2)
I2＝λ1 2λ2 2+λ2 2λ3 2+λ1 2λ3 2 (3)
Wherein I 1 is a first strain tensor invariant and I 2 is a second strain tensor invariant;
6.2 The strain energy function for superelastic rubber is obtained by the taylor formula expansion:
Assuming that the rubber material used for GINA water stop 3 is incompressible, j=1, and C ij is the rubber characteristic parameter: i+j=1, i.e., j=0 or 0 when i=1; when i=0, j=1;
For the full polynomial, when n=1, the strain energy of the remaining linear portion is preserved, i.e., as a function of strain energy density of the Mooney-Rivlin constitutive model:
W＝C10(I1-3)+C01(I2-3) (4)
Wherein, C 10 and C 01 are Rivlin coefficients of a Mooney-Rivlin model, and the values of the Rivlin coefficients are determined by a stress-strain relation curve of a full aging period;
6.3 Deriving Kirchhoff stress by biasing the elongation ratio, according to the Kirchhoff stress vs Green strain: by combining lambda 2 2＝λ3 2＝λ1 -1, the relationship between the contact stress sigma and the elongation ratio lambda after the aging of the GINA water stop is obtained:
The partial derivative of I 1、I2 is calculated by the formula (4) to obtain Relationship between the contact stress sigma elongation ratio lambda after aging of GINA water stop:
σ＝2(λ2-λ-1)(C10+C01λ-1) (6)
6.3 Carrying out parameter identification on the formula (6) by a nonlinear least square method by adopting Origin software based on a stress-strain relation curve of a full aging period to obtain values C 10 and C 01 under the full life period t i, obtaining functions f (t) and g (t) of changes of C 10 and C 01 along with aging time according to curves of changes of C 10 and C 01 along with aging time, and further constructing a constitutive model of stress relaxation and seawater aging of the GINA water stop belt:
σ＝2(λ2-λ-1)(f(t)+g(t)λ-1) (7)。
The following examples are given by way of illustration.
(1) A GINA water stop 3 with the length of 20.0cm multiplied by the width of 29.5cm multiplied by the height of 27.5cm is selected and placed on a bottom plate with the length of 60cm multiplied by the width of 35cm multiplied by the thickness of 20mm, a ballast plate with the length of 35cm multiplied by the width of 35cm multiplied by the thickness of 20mm and a pressing strip are used for restraining, a pressing machine with the rated load of 3000kN is used for loading the GINA water stop 3 with the hardness of 50HS with the shore hardness according to the standard of TB/T2843-2010, after the compression amount reaches 125mm, a middle plate with the length of 60cm multiplied by the width of 35cm multiplied by the thickness of 20mm is used for controlling the ballast plate, a Siemens QBE2103-P10 liquid gas pressure sensor 2 is placed on the middle plate, a top plate with the length of 60cm multiplied by 35cm multiplied by the thickness of 20mm is used for pressing, and a Q235 tool bolt with the size of M30 is selected for fixing the top plate, the middle plate and the bottom plate. The device was placed in a polypropylene plastic sealed silo 1 designed 80cm long by 50cm wide by 90cm high, filled with 2/3 of natural seawater, see fig. 1 and 2.
The attenuation conditions of the contact stress relaxation of the GINA water stop 3 in the seawater at 50 ℃, 60 ℃ and 70 ℃ are measured in real time through the liquid-gas pressure sensor 2 in the aging polypropylene plastic sealing bin 1, and stress relaxation curves at different aging temperatures are obtained, and the results are shown in fig. 4.
(2) Determining the contact stress sigma and the initial contact stress sigma 0 of the GINA water stop after aging according to the stress relaxation curve in the step (1), and determining the aging coefficient k of the GINA water stop according to k=sigma/sigma 0; determining an aging performance change value P of rubber used for the GINA water stop according to the P=sigma-sigma 0;
(3) Fitting according to the time Wen Diejia principle to obtain a lnP curve which changes along with t at normal temperature, obtaining P=exp (f (t)), obtaining a curve of an aging coefficient k at each temperature along with aging time, referring to fig. 5, obtaining P= (1-k) sigma 0 according to k=sigma/sigma 0＝(σ0-P)/σ0, and further determining a normal temperature aging coefficient k Often times ＝1-exp(f(t))/σ0 of the GINA water stop belt.
(4) The GINA water stop 3 was subjected to a uniaxial tensile test, the GINA water stop in the sea water aging bin was taken out, and was subjected to liquid nitrogen cooling, and a dumbbell-shaped test piece (the total length of the test piece was 100mm, the thickness of the test piece was 2.0mm, and the initial test length of the test piece test section was 20.0 mm) was prepared by using a punch mill according to GB/T528-2009, as shown in FIG. 3.
The GINA water stop 3 was subjected to uniaxial tension test using a MTSYAW6306,6306 electrohydraulic servo universal tester (3000 kN). First, a lubricant was coated on a selected test piece to reduce friction, a tensile test was performed at a speed of 500mm/min at room temperature of 23 ℃, the force measured by a sensor on a platen was divided by the original cross-sectional area of the GINA water stop 3, and the displacement measured by a displacement sensor was divided by the original height of the GINA water stop to obtain a stress-strain relationship curve of fig. 6.
(5) And correcting the stress-strain relation curve obtained by the uniaxial tensile test by adopting a normal temperature aging coefficient k Often times to obtain a stress-strain relation curve of a full aging period, as shown in figure 7.
(6) Determining the elongation ratio lambda of the GINA water stop according to the stress-strain relation curve of the full aging period in the step (5), and constructing a constitutive model of stress relaxation and seawater aging of the GINA water stop by utilizing a Mooney-Rivlin model, wherein the constitutive model specifically comprises the following steps:
6.1 First, green strain invariants I 1 and I 2 are found from the elongation ratios λ 1、λ2 and λ 3 in 3 directions, as in the formulae (1) to (3):
λi＝1+εi (1)
I1＝λ1 2+λ2 2+λ3 2 (2)
I2＝λ1 2λ2 2+λ2 2λ3 2+λ1 2λ3 2 (3)
Wherein, during uniaxial stretching, ∈ 2＝ε3 = -1, i.e. λ 2＝λ3 =0, then λ i＝λ1, i.e. λ 1＝ε1 +1;
6.2 Determining a strain energy density expression according to Green strain invariants I 1 and I 2
W＝C10(I1-3)+C01(I2-3) (4)
Solving the strain energy density;
the approximate fit using the strain energy density function formula is:
Wherein: c 10＝0.386e0.0008t,C01＝－0.203e0.0008t, t is the test time (h), j=1 assuming that the rubber material used for GINA water stop 3 is incompressible, and C ij is the rubber characteristic parameter: i+j=1, i.e., j=0 or 0 when i=1; when i=0, j=1.
6.3 Deriving Kirchhoff stress by biasing the elongation ratio, according to the Kirchhoff stress vs Green strain: by combining lambda 2 2＝λ3 2＝λ1 -1, the relationship between the contact stress sigma and the elongation ratio lambda after the aging of the GINA water stop is obtained:
Referring to Table 1, the partial derivative of I 1、I2 is calculated by equation (5) to give Relationship between the contact stress sigma elongation ratio lambda after aging of GINA water stop:
σ＝2(λ2-λ-1)(C10+C01λ-1) (6)
TABLE 1 stress-strain values and draw ratios obtained by uniaxial tensile test
6.5 Carrying out parameter identification on the formula (7) by a nonlinear least square method based on a stress-strain relation curve of a full aging period by adopting Origin software, referring to table 2, obtaining values of C 10 and C 01 at different aging times t i, and determining curves of C 10 and C 01 changing along with the aging time, referring to fig. 8;
TABLE 2 Material constants C 10 and C 01 change with aging time
According to the curves of C 10 and C 01 changing along with aging time, the functions f (t) and g (t) of C 10 and C 01 changing along with time are regressed, and then the aging time-varying GINA water stop belt stress relaxation and the seawater aging constitutive model are considered:
σ＝2(λ2-λ-1)(f(t)+g(t)λ-1) (7)；
The constructed model of the stress relaxation and seawater aging of the GINA water stop belt constructed by the method can be substituted into ANSYS finite element software through programming or commands, a Solid185 unit simulation rubber sample is adopted, the aging process of the GINA water stop belt 3 is simulated, the aging life of the GINA water stop belt 3 is predicted, and the parameters of the corrected constructed model are verified through the comparative analysis of test and numerical simulation results.
Namely, according to the change curve of the aging coefficient k with aging time at different aging temperatures (50 ℃,60 ℃, 70 ℃, 80 ℃ (in the figure, the unit of the temperature is used in the specification), the temperature of 50 ℃ is used as a reference temperature, namely alpha 323K =1, and the value b is obtained by adopting Origin software to carry out linear regression with x=1/T 0-1/T,y＝lnαT, and further, alpha T Often times ＝exp[b(1/T0 -1/T at normal temperature can be obtained through the formula lnalpha=b (1/T 0 -1/T).
According to the fitting equation lnP =f (T) of the mechanical properties at T 0 =50 ℃ with aging time, T 0,t0＝t Often times /αT Often times = 58450 days at p=σ 0-kpw is calculated, that is, the service life of the GINA water stop 3 at normal temperature is 160 years.
Through the patent of the application, the aging life of the obtained GINA water stop can reach 160 years. The life obtained by the patent of the application is more close to the real situation.
The foregoing description of the preferred embodiments of the application is not intended to be limiting, but rather is intended to cover all modifications, equivalents, and alternatives falling within the spirit and principles of the application.
Claims (5)
1. The method for constructing the time-varying constitutive model of the immersed tube tunnel GINA water stop with seawater ageing is characterized by comprising the following steps of:
1) Selecting a GINA water stop test piece, and performing seawater accelerated aging tests at different temperatures under the designed compression quantity on the GINA water stop test piece to obtain stress relaxation curves at different aging temperatures;
2) Determining the contact stress sigma and the initial contact stress sigma 0 of the GINA water stop after aging according to the stress relaxation curve in the step 1), and determining the aging coefficient k of the GINA water stop according to k=sigma/sigma 0; determining an aging performance change value P of rubber used for the GINA water stop according to the P=sigma-sigma 0;
3) Obtaining a lnP curve changing along with time t at normal temperature according to a time-temperature superposition principle, and obtaining an equation P=exp (f (t)), wherein f (t) is a change function of an aging performance change value P and time t; determining P= (1-k) sigma 0 according to the relation between the aging coefficient k of the GINA water stop and the aging property change value P of rubber used by the GINA water stop, and obtaining the normal temperature aging coefficient k Often times ＝1-exp(f(t))/σ0 of the GINA water stop;
4) Carrying out a uniaxial tensile test on the GINA water stop at the environment temperature of 23 ℃ and the tensile speed of 500mm/min to determine a stress-strain relation curve of the GINA water stop;
5) Correcting the stress-strain relation curve in the step (4) by using a normal temperature ageing coefficient k Often times of the GINA water stop belt to obtain a stress-strain relation curve of a full ageing period;
6) According to the stress-strain relation curve of the full aging period in the step 5), a constitutive model of the stress relaxation and seawater aging of the GINA water stop is built by combining a Mooney-Rivlin model:
σ＝2(λ2-λ-1)(f(t)+g(t)λ-1)；
Wherein t is aging time; lambda is the elongation ratio of the GINA water stop; f (t) is a function of C 10 over time t; g (t) is a function of C 01 over time t; c 10 and C 01 are Rivlin coefficients of the Mooney-Rivlin model, the values of which are determined by stress-strain relationship curves for the full aging period.
2. The method for constructing the time-varying constitutive model of the marine aged immersed tube tunnel GINA water stop according to claim 1, wherein the step (6) is specifically:
6.1 Determining epsilon i according to the uniaxial tensile test of the GINA water stop in the step 4), determining 3 elongation ratios lambda i of the GINA water stop according to epsilon i, respectively marking as lambda 1、λ2 and lambda 3, and obtaining Green strain tensor invariants I 1 and I 2, wherein the following formula is shown in the specification:
λi＝1+εi (1)
I1＝λ1 2+λ2 2+λ3 2 (2)
I2＝λ1 2λ2 2+λ2 2λ3 2+λ1 2λ3 2 (3)
Wherein I 1 is a first strain tensor invariant and I 2 is a second strain tensor invariant; lambda 1、λ2 and lambda 3 represent the elongation ratio in the X-axis direction, the elongation ratio in the Y-axis direction, and the elongation ratio in the Z-axis direction, respectively; epsilon i is the strain, i=1, 2 or 3, wherein epsilon 1 is the strain in the X-axis direction, epsilon 2 is the strain in the Y-axis direction, and epsilon 3 is the strain in the Z-axis direction;
6.2 According to Green strain invariants I 1 and I 2, determining that the strain energy density W is:
W＝C10(I1-3)+C01(I2-3) (4)
Wherein, C 10 and C 01 are Rivlin coefficients of a Mooney-Rivlin model, and the values of the Rivlin coefficients are determined by a stress-strain relation curve of a full aging period;
6.3 Through Kirchhoff stress versus Green strain): the relation between the contact stress sigma and the elongation ratio lambda after the aging of the GINA water stop is obtained as follows:
The partial derivative of I 1、I2 is calculated by the formula (4) to obtain Relationship between the contact stress sigma elongation ratio lambda after aging of GINA water stop:
σ＝2(λ2-λ-1)(C10+C01λ-1) (6)
6.4 Carrying out parameter identification on the formula (6) by a nonlinear least square method by adopting Origin software based on a stress-strain relation curve of a full aging period to obtain values C 10 and C 01 under the full life period t i, obtaining functions f (t) and g (t) of changes of C 10 and C 01 along with aging time according to curves of changes of C 10 and C 01 along with aging time, and further constructing a constitutive model of stress relaxation and seawater aging of the GINA water stop belt:
σ＝2(λ2-λ-1)(f(t)+g(t)λ-1) (7)。
3. The method for constructing the time-varying constitutive model of the immersed tube tunnel GINA water stop with seawater aging according to claim 2, wherein in the step (1), conditions of the seawater accelerated aging test are as follows:
Placing a water stop with the Shore hardness of 50HS on a bottom plate, restraining the water stop by using a ballast plate and a pressing strip, loading the water stop by using a press machine with the rated load of 3000kN according to the TB/T2843-2010 standard, controlling the ballast plate by using a middle plate after the compression amount reaches 125mm, placing a liquid gas pressure sensor on the middle plate, pressing the middle plate by using a top plate, fixing the top plate, the middle plate and the bottom plate, integrally placing the fixed top plate, the middle plate and the bottom plate in a polypropylene plastic sealing bin, filling 2/3 of natural seawater, and respectively collecting the contact stress of the GINA water stop in the seawater in real time within the range of 50-80 ℃.
4. The method for constructing the seawater-aged immersed tube tunnel GINA water stop time-varying constitutive model according to claim 2, wherein the uniaxial tensile test of the step (4) is specifically carried out under the following conditions:
Taking out the GINA water stop tape aged by the seawater in the step (1), cooling by liquid nitrogen, preparing a dumbbell-shaped test piece according to GB/T528-2009, coating a lubricant on the test piece, and carrying out a tensile test by a universal tester at a speed of 500mm/min at the room temperature of 23 ℃.
5. The method for constructing the marine aged immersed tube tunnel GINA water stop belt time-varying constitutive model according to claim 4, wherein the total length of the dumbbell-shaped test piece is 100mm, the thickness is 2.0mm, and the initial test length of the test section is 20.0mm.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210410929.6A CN114818180B (en) | 2022-04-19 | 2022-04-19 | Construction method of sea water aging immersed tube tunnel GINA water stop time-varying constitutive model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210410929.6A CN114818180B (en) | 2022-04-19 | 2022-04-19 | Construction method of sea water aging immersed tube tunnel GINA water stop time-varying constitutive model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114818180A CN114818180A (en) | 2022-07-29 |
| CN114818180B true CN114818180B (en) | 2024-05-03 |
Family
ID=82505165
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210410929.6A Active CN114818180B (en) | 2022-04-19 | 2022-04-19 | Construction method of sea water aging immersed tube tunnel GINA water stop time-varying constitutive model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114818180B (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110489904A (en) * | 2019-08-26 | 2019-11-22 | 石家庄铁道大学 | A kind of Studying Metal Rubber Isolator System kinetic model and parameter identification method |
| CN111089793A (en) * | 2018-10-08 | 2020-05-01 | 中国石油化工股份有限公司 | Hydrogenated nitrile butadiene rubber two-parameter constitutive model C10、C01Determination method |
-
2022
- 2022-04-19 CN CN202210410929.6A patent/CN114818180B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111089793A (en) * | 2018-10-08 | 2020-05-01 | 中国石油化工股份有限公司 | Hydrogenated nitrile butadiene rubber two-parameter constitutive model C10、C01Determination method |
| CN110489904A (en) * | 2019-08-26 | 2019-11-22 | 石家庄铁道大学 | A kind of Studying Metal Rubber Isolator System kinetic model and parameter identification method |
Non-Patent Citations (1)
| Title |
|---|
| 沉管隧道GINA橡胶止水带数值模拟分析;黄帆;;结构工程师;20100228(第01期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114818180A (en) | 2022-07-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Gillen et al. | Predicting and confirming the lifetime of o-rings | |
| Gillen et al. | Validation of improved methods for predicting long-term elastomeric seal lifetimes from compression stress–relaxation and oxygen consumption techniques | |
| Ronan et al. | Long-term stress relaxation prediction for elastomers using the time–temperature superposition method | |
| CN102865952B (en) | Nondestructive testing method for working stress of concrete | |
| Djebli et al. | Experimental analysis and damage modeling of high-density polyethylene under fatigue loading | |
| Bourchak et al. | PE-HD fatigue damage accumulation under variable loading based on various damage models. | |
| Drozdov | Swelling of thermo-responsive hydrogels | |
| CN104914041A (en) | Aging testing method of shield tunnel elastic sealing gasket finished products | |
| CN114818180B (en) | Construction method of sea water aging immersed tube tunnel GINA water stop time-varying constitutive model | |
| Gehrmann et al. | Displacement-controlled fatigue testing of rubber is not strain-controlled | |
| Uzundurukan et al. | Suction and swell characteristics of compacted clayey soils | |
| Johlitz et al. | Chemical ageing of elastomers: experiments and modelling | |
| CN105158147B (en) | A kind of device and method for being used to test sealing ring material aging | |
| Lacroix et al. | A local criterion for fatigue crack initiation on chloroprene rubber: approach in dissipation | |
| WO2023201522A1 (en) | Seawater aging immersed tunnel gina water stop time-varying constitutive model construction method | |
| Cangémi et al. | A two-phase model for the mechanical behavior of semicrystalline polymers | |
| Dewimille et al. | Behaviour of thermoplastic polymers during explosive decompressions in a petroleum environment | |
| KR100927738B1 (en) | Low cycle environmental fatigue test method using small fatigue specimen | |
| Lacasse et al. | Static and dynamic cut growth fatigue characteristics of silicone-based elastomeric sealants | |
| Zhang et al. | Characterization of ductile damage in polyethylene based on effective stress concept | |
| Jiang | A damaged evolution model for strain fatigue of ductile metals | |
| Utkan et al. | Compression behavior of typical silicone rubbers for soft robotics applications at elevated temperatures | |
| Lai et al. | Combined tension‐torsion creep experiments on polycarbonate in the nonlinear range | |
| Das et al. | Nonlinear behavior of a medium density polyethylene pipe material | |
| Sheikhi et al. | Finite element analysis of different material models for polyurethane elastomer using estimation data sets |
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 |