CN110442943B - Method and device for simulating underwater filling in finite element calculation of stress deformation of filling body - Google Patents
Method and device for simulating underwater filling in finite element calculation of stress deformation of filling body Download PDFInfo
- Publication number
- CN110442943B CN110442943B CN201910687239.3A CN201910687239A CN110442943B CN 110442943 B CN110442943 B CN 110442943B CN 201910687239 A CN201910687239 A CN 201910687239A CN 110442943 B CN110442943 B CN 110442943B
- Authority
- CN
- China
- Prior art keywords
- filling
- water
- boundary
- level
- calculation
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The embodiment of the invention relates to a method and a device for simulating underwater filling in finite element calculation of stress deformation of a filling body, wherein the method comprises the following steps: aiming at the water filling, dividing the water filling into a plurality of water filling calculation levels, and aiming at any water filling calculation level, executing the following steps: marking whether the lower boundary of the water in the water filling calculation level is a new boundary or an old boundary, and setting seepage boundary conditions and water pressure boundary conditions of the water boundary in the water filling calculation level; calculating water pressure increment according to the water pressure difference between the current water level and the initial water level, and calculating the water pressure generated by the current water level according to the water pressure load of the new boundary; and setting initial pore water pressure on a Gaussian point in the water filling calculation level, and calculating the dead weight load of the filling volume weight in the water filling calculation level according to a preset algorithm.
Description
Technical Field
The embodiment of the invention relates to the technical field of engineering such as dams and cofferdams in civil engineering, geotechnical engineering and hydraulic and hydroelectric engineering, in particular to a method and a device for simulating underwater filling in finite element calculation of filling body stress deformation.
Background
In the design stage of structures such as dams and cofferdams, stress deformation simulation calculation is usually required so as to demonstrate the rationality of design. The boundaries of the underwater filling bodies are constantly changed in the filling process due to the underwater filling stage of the dam, particularly the cofferdam, and the water level may also be changed in the filling process due to the difference of the water levels upstream and downstream of the filling caused by the flow of water. Moreover, the underwater filling is used as a link in the whole process of dam and cofferdam construction and operation, and the simulation of the whole process is needed.
In the related technology, generally for underwater filling, floating volume weight calculation is adopted underwater, and influence of water on stress deformation of a filling body in the filling process is considered by a wet volume weight calculation method on water. However, the simulation of the stress deformation of the filling body during the filling process is difficult to simulate such a complicated situation, and the simulation of the influence caused by the filling process in water and the water level change needs to be considered to simulate the stress deformation of the filling body and the change process of the pore water pressure in the filling body during the filling process.
Disclosure of Invention
In view of the above, to solve the above technical problems or some technical problems, embodiments of the present invention provide a method and an apparatus for simulating underwater filling in a finite element calculation of stress deformation of a filling body.
In a first aspect, an embodiment of the present invention provides a method for including a simulation of an underwater filling process in a finite element calculation of stress deformation of a filling body, where the method includes:
the method comprises the following steps of aiming at the whole construction and application process of structures such as dams, cofferdams and the like including the underwater filling process, dividing the whole process into a plurality of calculation levels including the underwater filling calculation level, and aiming at any calculation level, executing the following steps:
whether the boundary of each calculation level is a new boundary formed by filling of the calculation level or an old boundary formed by filling of the previous level is marked, and seepage boundary conditions and water pressure boundary conditions of the boundary in each calculation level are set;
during the finite element calculation of each stage, calculating the water pressure increment of the old boundary according to the water pressure difference between the current water level and the initial water level of the stage, calculating the water pressure increment of the new boundary, and calculating the water pressure generated by the current water level as the water pressure increment;
and setting the initial pore water pressure on the Gaussian point in each filling calculation level, and calculating the dead weight load of the filling volume weight in the water filling calculation level according to a preset algorithm.
In one possible embodiment, the method further comprises:
and judging whether the underwater boundary of the previous calculation level in the underwater filling calculation level disappears due to the filling of the calculation level.
In one possible embodiment, the method further comprises:
for an underwater boundary that disappears as a result of the filling of the present filling calculation level, a water pressure corresponding to the initial upstream or downstream water level of the current filling calculation level is applied in reverse on the original boundary.
In one possible embodiment, the setting of the seepage boundary condition and the water pressure boundary condition of the boundary in the filling calculation stage includes:
and setting seepage boundary conditions and water pressure boundary conditions of the boundary in the water filling calculation level according to the new boundary condition in the filling calculation level.
In one possible embodiment, the setting the initial pore water pressure at the gaussian point in the water fill calculation stage comprises:
setting the initial pore water pressure on the Gaussian point of the filling unit in the filling calculation level as the negative value of the suction corresponding to the filling saturation, wherein the saturation of the Gaussian point under water is equal to 1, and the suction is 0.
In one possible embodiment, the calculating the dead weight load of the filling volume weight in the filling calculation level in water according to a preset algorithm includes:
the dead weight load of the fill bulk weight in the fill calculation stage is calculated from the wet bulk weight (the total weight of soil particles and water per unit volume of the fill).
In a second aspect, an embodiment of the present invention provides an apparatus for simulating a filling in water in a stress deformation finite element calculation of a filling body, where the apparatus includes:
the calculation level division module is used for dividing the water filling into a plurality of water filling calculation levels;
the boundary marking module is used for marking whether the boundary in the water filling calculation level is a new boundary or an old boundary aiming at any water filling calculation level;
the condition setting module is used for setting seepage boundary conditions and water pressure boundary conditions of boundaries in the water filling calculation levels aiming at any water filling calculation level;
the first calculation module is used for calculating the water pressure increment of an old boundary and the water pressure load of a new boundary according to the water pressure difference between the current water level and the initial water level aiming at any water filling calculation level, and calculating the water pressure generated by the current water level as the water pressure increment load;
the initial pore water pressure setting module is used for setting the initial pore water pressure on a Gaussian point in a filling calculation level aiming at any filling calculation level;
and the second calculation module is used for calculating the dead weight load of the filling volume weight in the water filling calculation level according to a preset algorithm aiming at any water filling calculation level.
In one possible embodiment, the apparatus further comprises:
and the judging module is used for judging whether the underwater boundary of the upper computing level in the filling computing level disappears due to filling in the filling computing level aiming at any filling computing level.
In one possible embodiment, the apparatus further comprises:
and the boundary processing module is used for reversely applying water pressure corresponding to the initial upstream and downstream water levels of the current filling calculation level on the original boundary aiming at the underwater boundary of the previous calculation level disappeared due to the filling of the filling calculation level.
In a possible embodiment, the condition setting module is specifically configured to:
and setting seepage boundary conditions and water pressure boundary conditions of the underwater boundary in the water filling calculation level according to the new boundary condition in the water filling calculation level.
According to the technical scheme provided by the embodiment of the invention, according to the change of the water level and the boundary in the filling process, the pore water pressure and the water pressure increment of the filling body boundary are simulated and the initial pore water pressure of the filling body unit is set during the calculation of any filling calculation level, the volume weight of the filled soil body does not need to distinguish water from water, the dead weight load is uniformly calculated according to the wet volume weight, and the simulation of the seepage and the stress deformation of the complicated filling of the filling body in the water and the water level change process is realized.
Drawings
In order to more clearly illustrate the embodiments of the present specification or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments described in the embodiments of the present specification, and other drawings can be obtained by those skilled in the art according to the drawings.
FIG. 1 is a schematic flow chart of an embodiment of a method for simulating underwater filling in a finite element calculation of stress deformation of a filling body according to the present invention;
FIG. 2 is a schematic structural diagram of an apparatus for simulating filling in water in a finite element calculation of stress deformation of a filling body according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
For the convenience of understanding of the embodiments of the present invention, the following description will be further explained with reference to specific embodiments, which are not to be construed as limiting the embodiments of the present invention.
As shown in fig. 1, an implementation flow diagram of a method for simulating underwater filling in a stress deformation finite element calculation of a filling body according to an embodiment of the present invention is provided, and the method specifically includes the following steps:
s101, aiming at filling in water, dividing the filling into a plurality of filling calculation levels in water, and aiming at any filling calculation level, executing the following steps:
in the embodiment of the invention, for the underwater filling simulation, the underwater filling simulation is divided into a plurality of underwater filling calculation levels, seepage and stress deformation conditions in the filling process of a simulation filling body and the upstream and downstream water level change process are calculated, and the following steps are executed for any calculation level comprising the underwater filling calculation levels.
S102, marking whether the lower boundary of the water in the water filling calculation level is a new boundary or an old boundary, and setting seepage boundary conditions and water pressure boundary conditions of the water boundary in the water filling calculation level;
marking whether the lower boundary of the water in the water filling calculation level is a new boundary or an old boundary aiming at any water filling calculation level;
aiming at any underwater filling calculation level, judging whether an underwater boundary in the previous calculation level disappears due to filling in the current calculation level;
setting seepage boundary conditions and water pressure boundary conditions of an underwater boundary in the water filling calculation level aiming at any water filling calculation level;
the seepage boundary condition and the water pressure boundary condition of the underwater boundary in the water filling calculation level are set according to the new boundary condition in the water filling calculation level, namely, the seepage boundary condition and the water pressure boundary condition of the underwater boundary in each level of calculation, and the seepage boundary condition and the water pressure increment load are given and calculated according to the new boundary condition of each level.
S103, calculating water pressure increment according to the water pressure difference between the current water level and the initial water level by using the water pressure load of the old boundary, calculating the water pressure generated by the current water level as the incremental water pressure load by using the water pressure load of the new boundary;
for any calculation level, the water pressure increment is calculated according to the water pressure difference between the current water level and the initial water level for the water pressure load of the old boundary, and the water pressure generated by the current water level is calculated according to the water pressure load of the new boundary.
And S104, setting the initial pore water pressure on the Gaussian point in the water filling calculation level, and calculating the dead weight load of the filling volume weight in the water filling calculation level according to a preset algorithm.
Aiming at any one filling calculation level in water, setting initial pore water pressure on a Gaussian point in the filling calculation level in water in the process of forming an equation set for each calculation level in a circulation manner at the first calculation time step of each filling calculation level in water. Wherein, the initial pore water pressure on the Gaussian point in the calculation level of the filling in water is set as the negative value of the suction corresponding to the filling saturation (the filling saturation of the filling in water is 1, and the suction is 0).
And for the dead load of the filling volume weight in the water filling calculation level, calculating according to a preset algorithm. The dead weight of the filling volume weight in the water filling calculation level is calculated according to the wet volume weight, namely the dead weight is calculated by adopting the wet volume weight (the saturation of the filling body in the water is 1 and is equal to the saturated volume weight) to calculate the dead weight.
In addition, at the first calculation time step of each water filling calculation stage, for the underwater boundary of the previous calculation stage disappeared due to the filling of the current stage, the water pressure corresponding to the initial upstream and downstream water levels of the current water filling calculation stage is reversely applied on the original boundary, namely, during the first time step calculation of each calculation stage, for the underwater boundary of the previous calculation stage disappeared due to the filling of the current stage, the water pressure corresponding to the initial upstream and downstream water levels of the current stage is reversely applied on the primary boundary, so as to offset the more than actual pressure increment load transmitted through the filling interface.
Through the above description of the technical scheme provided by the embodiment of the invention, according to the change of the water level and the boundary in the filling process, the simulation of the pore water pressure and the water pressure increment of the filling body boundary and the initial pore water pressure of the filling body unit are set during the calculation of the filling calculation level in any water, the volume weight of the filling body does not need to distinguish water from water, the dead weight load is uniformly calculated, and the simulation of the seepage and the stress deformation of the complicated filling and water level change process of the filling body in the water is realized.
With respect to the method embodiment, an embodiment of the present invention further provides an embodiment of an apparatus for simulating filling in water in a finite element calculation of stress deformation of a filling body, as shown in fig. 2, the apparatus may include: a calculation level division module 210, a boundary marking module 220, a condition setting module 230, a first calculation module 240, a pressure setting module 250, and a second calculation module 260.
A calculation level division module 210, configured to divide the water filling into a plurality of water filling calculation levels;
a boundary marking module 220, configured to mark, for any submerged filling calculation level, whether a lower boundary in the submerged filling calculation level is a new boundary or an old boundary;
a condition setting module 230, configured to set, for any one of the filling calculation levels in water, a seepage boundary condition and a water pressure boundary condition of a boundary in the filling calculation level;
a first calculation module 240, configured to calculate, for any calculation level, a water pressure increment according to a water pressure difference between the current water level and the initial water level for a water pressure load of an old boundary, a water pressure load of a new boundary, and a water pressure generated by the current water level;
an initial pore water pressure setting module 250 configured to set an initial pore water pressure on a gaussian point for any filling calculation level;
and the second calculating module 260 is configured to calculate, according to a preset algorithm, a dead weight of the filling volume weight in the filling calculating stage for any filling calculating stage.
According to a specific embodiment provided by the present invention, the apparatus further comprises:
and the judging module 270 is configured to judge, for any one underwater filling calculation level, whether the underwater boundary of the previous calculation level in the underwater filling calculation level disappears due to the filling of the current level.
According to a specific embodiment provided by the present invention, the apparatus further comprises:
and a boundary processing module 280 for reversely applying a water pressure corresponding to an initial upstream or downstream water level of the current computation stage on the original boundary with respect to the underwater boundary of the previous computation stage disappeared by the filling of the current computation stage.
According to a specific embodiment of the present invention, the condition setting module 230 is specifically configured to:
and setting seepage boundary conditions and water pressure boundary conditions of the calculation stages according to the new boundary conditions in any calculation stage.
Those of skill would further appreciate that the various illustrative components and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, computer software, or combinations of both, and that the various illustrative components and steps have been described above generally in terms of their functionality in order to clearly illustrate this interchangeability of hardware and software. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the implementation. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
The steps of a method or algorithm described in connection with the embodiments disclosed herein may be embodied in hardware, a software module executed by a processor, or a combination of the two. A software module may reside in Random Access Memory (RAM), memory, Read Only Memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are merely exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (10)
1. A method for simulating an underwater filling in a stress deformation finite element calculation of a filling body, the method comprising:
aiming at the water filling, dividing the water filling into a plurality of water filling calculation levels, and aiming at any filling calculation level, executing the following steps:
marking whether the boundary in the calculation level is a new boundary or an old boundary, and setting seepage boundary conditions and water pressure boundary conditions of the underwater boundary in the filling calculation level;
calculating water pressure increment according to the water pressure difference between the current water level and the initial water level, and calculating the water pressure generated by the current water level according to the water pressure load of the new boundary;
and setting initial pore water pressure on the Gaussian points of the filling units in the calculation stage, and calculating the self-weight load of the filling unit at the current stage according to a preset algorithm.
2. The method of claim 1, further comprising:
and judging whether the underwater boundary in the previous calculation level disappears due to filling of the calculation level.
3. The method of claim 2, further comprising:
for an underwater boundary of a previous computation stage that disappears due to the filling of the current computation stage, a water pressure acting on the boundary corresponding to an initial upstream or downstream water level of the current computation stage is applied in reverse on the original boundary.
4. The method of claim 1, wherein setting seepage boundary conditions and water pressure boundary conditions for the underwater boundary in the water fill computation level comprises:
and setting seepage boundary conditions and water pressure boundary conditions of the underwater boundary in the water filling calculation level according to the new boundary condition in the water filling calculation level.
5. The method of claim 1, wherein setting an initial pore water pressure on a gaussian spot of a fill forming cell in a fill computation level in water comprises:
setting the initial pore water pressure on the Gaussian point of the filling forming unit in the calculation stage as the negative value of the suction corresponding to the filling saturation.
6. The method of claim 1, wherein calculating the deadweight load of the fill volume weight in the fill computation level in water according to a predetermined algorithm comprises:
and calculating the dead weight load of the filling volume weight in the filling calculation level in the water according to the wet volume weight.
7. A device for simulating a filling in water in a finite element calculation of stress deformation of a filling body, the device comprising:
the calculation level division module is used for dividing the water filling into a plurality of water filling calculation levels;
the boundary marking module is used for marking whether the underwater boundary is a new boundary or an old boundary aiming at any calculation level;
the condition setting module is used for setting seepage boundary conditions and water pressure boundary conditions of an underwater boundary in a calculation level aiming at any underwater filling calculation level;
the first calculation module is used for calculating the water pressure increment according to the water pressure difference between the current water level and the initial water level and the water pressure load of the new boundary and calculating the water pressure generated by the current water level aiming at the water pressure load of the old boundary at any calculation level;
the pressure setting module is used for setting the initial pore water pressure on a Gaussian point in the water filling calculation level aiming at any filling calculation level;
and the second calculation module is used for calculating the dead weight load of the filling volume weight in the calculation level according to a preset algorithm aiming at any calculation level.
8. The apparatus of claim 7, further comprising:
and the judging module is used for judging whether the underwater boundary in the previous calculation level disappears due to the filling of the current level aiming at any calculation level.
9. The apparatus of claim 8, further comprising:
and the boundary processing module is used for reversely applying water pressure corresponding to the initial upstream and downstream water levels of the current underwater filling calculation level on the original boundary aiming at the underwater boundary of the previous calculation level disappeared due to the filling of the current calculation level.
10. The apparatus according to claim 7, wherein the condition setting module is specifically configured to:
and setting seepage boundary conditions and water pressure boundary conditions of the underwater boundary in the water filling calculation level according to the new boundary condition in the water filling calculation level.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910687239.3A CN110442943B (en) | 2019-07-29 | 2019-07-29 | Method and device for simulating underwater filling in finite element calculation of stress deformation of filling body |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910687239.3A CN110442943B (en) | 2019-07-29 | 2019-07-29 | Method and device for simulating underwater filling in finite element calculation of stress deformation of filling body |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110442943A CN110442943A (en) | 2019-11-12 |
| CN110442943B true CN110442943B (en) | 2021-02-02 |
Family
ID=68431876
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910687239.3A Active CN110442943B (en) | 2019-07-29 | 2019-07-29 | Method and device for simulating underwater filling in finite element calculation of stress deformation of filling body |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110442943B (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109376429A (en) * | 2018-10-24 | 2019-02-22 | 中国水利水电第七工程局有限公司 | A kind of concrete dam template safe construction analysis method based on finite element simulation |
| CN109992900A (en) * | 2019-04-08 | 2019-07-09 | 清华大学 | A kind of more real-time online cooperative intelligent emulation modes of mass concrete and system |
-
2019
- 2019-07-29 CN CN201910687239.3A patent/CN110442943B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109376429A (en) * | 2018-10-24 | 2019-02-22 | 中国水利水电第七工程局有限公司 | A kind of concrete dam template safe construction analysis method based on finite element simulation |
| CN109992900A (en) * | 2019-04-08 | 2019-07-09 | 清华大学 | A kind of more real-time online cooperative intelligent emulation modes of mass concrete and system |
Non-Patent Citations (2)
| Title |
|---|
| Analysis of stress and deformation of rockfill and concrete face for high concrete face rockfill dam;Shun Wen Ji 等;《Applied Mechanics and Materials》;20140930;第638-640卷;第731-734页 * |
| 考虑渗流效应的大型土石围堰稳定性研究;董存军;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20130315(第3期);第C037-42页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110442943A (en) | 2019-11-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Jeuken et al. | Impact of dredging and dumping on the stability of ebb–flood channel systems | |
| Roeber et al. | Shock-capturing Boussinesq-type model for nearshore wave processes | |
| CN110362925B (en) | Method for simulating overtopping and bursting flood numerical value of earth and rockfill dam comprising reservoir area | |
| KR20230055060A (en) | Two-Dimensional Flood inundation analysis technique | |
| Altomare et al. | An experimental study on scale effects in wave reflection of low-reflective quay walls with internal rubble mound for regular and random waves | |
| CN110442943B (en) | Method and device for simulating underwater filling in finite element calculation of stress deformation of filling body | |
| CN104794753A (en) | Earth rock cofferdam breach flood risk mapping method | |
| CN105970975B (en) | A kind of base pit stand construction with basement building | |
| CN110472309B (en) | Method and device for considering three-dimensional percolation effect in two-dimensional percolation calculation model | |
| CN114861511A (en) | Method, system, equipment and medium for simulating deformation and creep of offshore seabed sediment under action of wave-induced load | |
| KR101273977B1 (en) | Dynamic calculation range type numerical simulation method of tsunami | |
| CN110795792B (en) | Method for predicting riverbed deformation of strong turbulent fluctuation area of river channel caused by engineering construction | |
| Chang et al. | Effective stress analysis of a dynamically installed anchor penetrating into rate-dependent soil | |
| Beiranvand et al. | Effect of different parameters of heterogeneous dams on safety factor using the neural network. Case study: Marvak dam | |
| CN112948915B (en) | Generalization processing method for vertical wading building in numerical simulation test | |
| Toromanovic et al. | Effects on an earth and rockfill dam undergoing dam safety measures | |
| Neves et al. | Overtopping of a porous strcuture using a smoothed particle hydrodynamics numerical model | |
| Klerk et al. | Influence of load interdependencies of flood defences on probabilities and risks at the Bovenrijn/IJssel area, The Netherlands | |
| Rasekh | Hydrodynamic pressure on culvert gates during an earthquake | |
| Androsov et al. | Tsunami modelling on unstructured grids: Verification and validation | |
| Sidorov et al. | Features of the port retaining structure calculations | |
| Wana et al. | Smoothed Particle Hydrodynamics (SPH) modelling of tsunami waves generated by a fault rupture | |
| Kim | Numerical Modeling of Solid Particle Dynamics Using the Moving Particle Semi-Implicit (MPS) Method in Seabed Penetration Scenarios | |
| CN107918722A (en) | A kind of Wall in Plain Reservoir well leaks computational methods | |
| Mitovski | Ljupcho Petkovski |
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 |