Reservoir basin anti-seepage panel dam composite geomembrane numerical simulation method, device, medium and product

Abstract

The invention discloses a reservoir basin seepage-proofing panel dam composite geomembrane numerical simulation method, device, medium and product, and the method comprises the steps: constructing a geomembrane plane model, dividing the geomembrane plane model into a plurality of plane isoparametric units, and constructing a spatial polyhedron isoparametric unit by taking each plane isoparametric unit as a top surface to obtain a reservoir basin finite element model; constructing a stiffness matrix of a plane isoparametric unit; combining the stiffness matrixes of the plane isoparametric unit and the spatial polyhedron isoparametric unit to obtain a stiffness matrix of the reservoir basin finite element model; constructing a dam finite element model, combining the reservoir basin finite element model with the dam finite element model, and combining the stiffness matrixes of the reservoir basin finite element model and the dam finite element model to obtain a reservoir basin anti-seepage panel dam composite geomembrane model and a stiffness matrix thereof; and numerical calculation is carried out based on the reservoir basin anti-seepage panel dam composite geomembrane model and the stiffness matrix thereof. The problem that geomembrane simulation is difficult to carry out on the dam scale is solved.

Description

Technical Field

[0001] The present invention belongs to the technical field of dam safety monitoring, and in particular relates to a numerical simulation method, equipment, medium and product of a composite geomembrane for a reservoir basin anti-seepage panel dam. Background Art

[0002] Pumped-storage power stations, as power facilities that utilize hydropower for energy storage and dispatch, can help balance grid loads and improve energy efficiency. Construction efforts have been increasing in recent years. The anti-seepage performance of the reservoir area in pumped-storage power stations is directly related to the operational safety and economic benefits of the power station. Reservoir basin anti-seepage measures, a form of anti-seepage that is suitable for reservoirs with complex geological conditions and deep impermeable layers, are gradually being applied in projects.

[0003] Common reservoir basin anti-seepage measures in pumped-storage power plants include paving panels on the dam and reservoir banks, and laying a composite geomembrane on the reservoir floor. These two, combined with overlapping and waterstop features, form a complete anti-seepage system. Considering that cracking of the geomembrane could pose a threat to the overall anti-seepage effectiveness of the project, geomembrane analysis and research are often required.

[0004] Currently, research on geomembrane performance typically focuses on permeability and strength. For example, studies are conducted on geomembrane leakage and the seepage characteristics of defective geomembranes, as well as experimental investigations of the tensile properties of unidirectional and multidirectional geomembranes. Finite element methods are commonly used to calculate stress and deformation in dams. However, due to the thinness of geomembranes and the significant differences in their mechanical properties from those of the dam material, geomembrane simulations at the dam scale are difficult.

[0005] At present, the following methods are mainly used to consider the role of geomembranes:

[0006] The calculation method of embedding solid elements in geomembrane is complex, and the calculation of solid elements is large when geomembrane is established. The displacement of unit surface nodes is converted into geomembrane strain, and the calculation is not accurate. The geomembrane is simulated by thin film elements, which is large in calculation and difficult to converge. Summary of the Invention

[0007] The purpose of the present invention is to provide a numerical simulation method for composite geomembranes of reservoir basin anti-seepage panel dams, so as to solve the problem that geomembrane simulation is difficult on the dam scale due to the thin thickness of the geomembrane and the large difference in mechanical properties from the dam material.

[0008] The present invention solves the above technical problems through the following technical solutions: a numerical simulation method for composite geomembrane of reservoir basin anti-seepage panel dam, comprising:

[0009] Constructing a geomembrane plane model in a global coordinate system, dividing the geomembrane plane model into a plurality of plane isoparametric units, and constructing a corresponding spatial polyhedron isoparametric unit with each plane isoparametric unit as the top surface to obtain a reservoir basin finite element model;

[0010] A local coordinate system is constructed based on each planar isoparametric element, and the stiffness matrix of each planar isoparametric element in the global coordinate system is constructed according to the coordinate position, displacement and nodal force of each planar isoparametric element in the global coordinate system and the coordinate position of each planar isoparametric element in the local coordinate system;

[0011] The stiffness matrix of the reservoir basin finite element model is obtained by combining the stiffness matrix of each plane isoparametric element in the global coordinate system and the stiffness matrix of each spatial polyhedron isoparametric element in the global coordinate system.

[0012] Construct the finite element model of the dam under reservoir basin geological conditions in the global coordinate system;

[0013] The reservoir basin finite element model is combined with the dam finite element model to obtain a composite geomembrane model of the reservoir basin anti-seepage face dam;

[0014] Combining the stiffness matrix of the reservoir basin finite element model with the stiffness matrix of the dam finite element model to obtain the stiffness matrix of the composite geomembrane model of the reservoir basin anti-seepage face dam;

[0015] Numerical calculations are performed based on the composite geomembrane model of the reservoir basin anti-seepage panel dam and its stiffness matrix.

[0016] Furthermore, the planar isoparametric unit is a planar quadrilateral isoparametric unit or a planar triangle isoparametric unit.

[0017] Furthermore, the specific construction process of the stiffness matrix of each planar isoparametric element in the global coordinate system includes:

[0018] According to the coordinate position of each plane isoparametric element in the global coordinate system and the local coordinate system, the rotation matrix between the global coordinate system and the local coordinate system is obtained;

[0019] According to the rotation matrix and the displacement of the vertex of each plane isoparametric unit in the global coordinate system, the displacement of the vertex of each plane isoparametric unit in the local coordinate system is obtained, and then the displacement matrix of each plane isoparametric unit in the local coordinate system is obtained;

[0020] According to the rotation matrix and the nodal force of the vertex of each plane isoparametric unit in the global coordinate system, the nodal force of the vertex of each plane isoparametric unit in the local coordinate system is obtained, and then the nodal force matrix of each plane isoparametric unit in the local coordinate system is obtained;

[0021] According to the displacement matrix and nodal force matrix of each plane isoparametric element in the local coordinate system, the stiffness matrix of each plane isoparametric element in the local coordinate system is obtained;

[0022] According to the rotation matrix and the stiffness matrix of each planar isoparametric element in the local coordinate system, the stiffness matrix of each planar isoparametric element in the global coordinate system is obtained.

[0023] Furthermore, the calculation formula of the rotation matrix is:

[0024]

[0025] Among them, R represents the rotation matrix between the global coordinate system and the local coordinate system, represents the coordinate position of the vertex of the plane isoparametric unit in the local coordinate system, and (x, y, z) represents the coordinate position of the vertex of the plane isoparametric unit in the global coordinate system.

[0026] Furthermore, the specific formula of the stiffness matrix of each plane isoparametric element in the local coordinate system is:

[0027]

[0028]

[0029] in, represents the stiffness matrix of the plane isoparametric element in the local coordinate system, represents the displacement matrix of the plane isoparametric element in the local coordinate system, represents the nodal force matrix of the plane isoparametric element in the local coordinate system, D represents the displacement matrix of the plane isoparametric element in the global coordinate system, L represents the nodal force matrix of the plane isoparametric element in the global coordinate system, and C represents the transformation matrix composed of the rotation matrix R.

[0030] Furthermore, the specific formula of the stiffness matrix of each plane isoparametric element in the global coordinate system is:

[0031]

[0032] Among them, S represents the stiffness matrix of the plane isoparametric element in the global coordinate system, C represents the transformation matrix composed of the rotation matrix R, Represents the stiffness matrix of the planar isoparametric element in the local coordinate system, and the superscript T represents the transpose.

[0033] Furthermore, numerical calculations are performed based on the composite geomembrane model of the reservoir basin anti-seepage panel dam and its stiffness matrix, specifically including:

[0034] Set boundary conditions and material properties, and apply loads to the geomembrane composite model of the reservoir anti-seepage face dam.

[0035] The displacement distribution of the composite geomembrane model of the reservoir basin anti-seepage face dam is calculated based on the load and the stiffness matrix of the composite geomembrane model of the reservoir basin anti-seepage face dam;

[0036] The stress distribution of the composite geomembrane model of the reservoir basin anti-seepage face plate dam is calculated based on the displacement distribution of the composite geomembrane model of the reservoir basin anti-seepage face plate dam.

[0037] Based on the same concept, the present invention also provides an electronic device, including a memory, a processor, and a computer program / instruction stored in the memory, wherein the processor executes the computer program / instruction to implement the numerical simulation method of the composite geomembrane of the reservoir basin anti-seepage panel dam as described above.

[0038] Based on the same concept, the present invention also provides a computer-readable storage medium having a computer program / instruction stored thereon, which, when executed by a processor, implements the above-mentioned numerical simulation method for the composite geomembrane of the reservoir basin anti-seepage panel dam.

[0039] Based on the same concept, the present invention also provides a computer program product, including a computer program / instruction, which, when executed by a processor, implements the above-mentioned numerical simulation method for the composite geomembrane of the reservoir basin anti-seepage panel dam.

[0040] Compared with the prior art, the advantages of the present invention are:

[0041] The present invention constructs a reservoir basin finite element model based on geomembrane plane isoparametric units and soil space polyhedron isoparametric units, and then combines the dam finite element model to construct a composite geomembrane model of the reservoir basin anti-seepage panel dam. Through the coupling calculation of geomembrane plane isoparametric units and soil space polyhedron isoparametric units, the simulation of the geomembrane on the dam scale is realized, and the stress and displacement distribution calculation of the composite geomembrane under different loads and boundary conditions is simulated, which reflects the tensile capacity of the geomembrane and solves the problem of difficulty in simulating the geomembrane on the dam scale. The present invention converts the mechanical properties of the geomembrane into the global coordinate system through the conversion of the stiffness matrix, reflects the mechanical contribution of the geomembrane to the overall model, and improves the calculation accuracy.

[0042] The present invention provides guidance or reference for geomembrane research and calculation, is conducive to optimizing the design parameters of composite geomembranes, reduces the workload of parameter adjustment in the numerical analysis of composite geomembranes for reservoir basin anti-seepage panel dams, optimizes the design and performance analysis process of composite geomembranes, and improves the efficiency and reliability of numerical simulations. BRIEF DESCRIPTION OF THE DRAWINGS

[0043] In order to more clearly illustrate the technical solution of the present invention, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only one embodiment of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0044] Figure 1 Flowchart of the numerical simulation method of the composite geomembrane of the reservoir basin anti-seepage panel dam according to the embodiment of the present invention;

[0045] Figure 2 Schematic diagram of a coupling model of a planar quadrilateral isoparametric unit and a spatial hexahedral isoparametric unit in an embodiment of the present invention;

[0046] Figure 3 This is a plan view of the dam in an embodiment of the present invention;

[0047] Figure 4 Schematic diagram of a composite geomembrane model of a reservoir anti-seepage panel dam according to an embodiment of the present invention;

[0048] Figure 5 Schematic diagram of the bottom settlement of the reservoir basin without geomembrane at normal water level in an embodiment of the present invention; wherein the numbers represent the displacement in the Z direction, i.e., the settlement displacement, in mm;

[0049] Figure 6 Schematic diagram of the bottom settlement of the reservoir basin when a geomembrane is installed at a normal water level in an embodiment of the present invention; wherein the numbers represent the displacement in the Z direction, i.e., the settlement displacement, in mm;

[0050] Figure 7 Schematic diagram of the deformation of the reservoir bottom in the river direction without a geomembrane at a normal water level in an embodiment of the present invention; wherein the numbers represent the displacement in the X direction, i.e., the deformation displacement in the river direction, in mm;

[0051] Figure 8 Schematic diagram of the deformation of the reservoir bottom in the river direction when a geomembrane is installed at a normal water level in an embodiment of the present invention; wherein the numbers represent the displacement in the X direction, i.e., the deformation displacement in the river direction, in mm;

[0052] Figure 9 Schematic diagram of the transverse river deformation of the reservoir bottom without a geomembrane at a normal water level in an embodiment of the present invention; wherein the numbers represent the Y-direction displacement, i.e., the transverse river deformation displacement, in mm;

[0053] Figure 10 Schematic diagram of the transverse deformation of the reservoir bottom when a geomembrane is present at a normal water level in an embodiment of the present invention. The numbers represent displacement in the Y direction, i.e., transverse deformation displacement, in mm.

[0054] Explanation of reference numerals: 1- cushion layer, 2- panel, 3- main rockfill, 4- downstream rockfill, 5- pressure slope, 6- geomembrane. DETAILED DESCRIPTION

[0055] The following is a clear and complete description of the technical solutions of the present invention in conjunction with the accompanying drawings of the embodiments of the present invention. Obviously, the embodiments described are only part of the embodiments of the present invention, not all of the embodiments. All other embodiments obtained by ordinary technicians in this field based on the embodiments of the present invention without making any creative efforts are within the scope of protection of the present invention.

[0056] The following specific embodiments are used to describe the technical solution of the present application in detail. The following specific embodiments can be combined with each other, and the same or similar concepts or processes may not be described in detail in some embodiments.

[0057] Example 1

[0058] In order to solve the problem that it is difficult to simulate geomembranes on a dam scale due to the thin thickness of geomembranes and the significant difference in mechanical properties from dam materials, the present invention provides a numerical simulation method for composite geomembranes of reservoir basin anti-seepage panel dams, which can better complete the coupling calculation of geomembranes and dams, and improve calculation accuracy and efficiency. Figure 1 The flow chart of the numerical simulation method of the composite geomembrane for the anti-seepage panel dam of the reservoir basin provided by the present invention is shown. Figure 1 As shown, the numerical simulation method of the composite geomembrane of the reservoir basin anti-seepage panel dam of the present invention includes the following steps:

[0059] Step S1: Construct a geomembrane plane model in a global coordinate system, divide the geomembrane plane model into multiple plane isoparametric units, and construct a corresponding spatial polyhedron isoparametric unit with each plane isoparametric unit as the top surface to obtain a reservoir basin finite element model.

[0060] Based on the actual geomembrane geometry, a geomembrane plane model is constructed in the global coordinate system using GID software. The geomembrane plane model is then subdivided into multiple plane isoparametric elements. In this embodiment, the plane isoparametric elements can be quadrilateral or triangular isoparametric elements. Quadrilateral isoparametric elements offer higher computational accuracy than triangular isoparametric elements. The greater the number of plane isoparametric elements, the higher the numerical simulation accuracy but the more complex the calculation.

[0061] During the geomembrane simulation process, the geomembrane adopts a plane isoparametric unit. The deformation of the geomembrane in the reservoir bottom area is uniform. The displacement and stress of the geomembrane under different working conditions are calculated without considering the changes in the physical properties of the geomembrane itself and ignoring the nonlinear response of the geomembrane material under different stresses.

[0062] In order to realize the coupling calculation between geomembrane and soil, a corresponding spatial polyhedron isoparametric unit (such as Figure 2 As shown in Figure 2), spatial polyhedron isoparametric elements are used to simulate the soil, thereby obtaining a finite element model of the reservoir basin. Figure 2 As shown, when the plane isoparametric unit is a plane quadrilateral isoparametric unit (i.e., a quadrilateral corresponding to vertices 1 to 4), the spatial polyhedron isoparametric unit is a spatial hexahedron isoparametric unit (i.e., a hexahedron corresponding to vertices 1 to 8).

[0063] Step S2: Construct a local coordinate system based on each planar isoparametric element, and construct the stiffness matrix of each planar isoparametric element in the global coordinate system according to the coordinate position, displacement and nodal force of each planar isoparametric element in the global coordinate system and the coordinate position of each planar isoparametric element in the local coordinate system.

[0064] The local coordinate system of each plane isoparametric unit is constructed with the geometric center or a vertex of each plane isoparametric unit as the coordinate origin, the plane where the plane isoparametric unit is located as the XY plane, and the vertical axis to the plane isoparametric unit as the Z axis.

[0065] In a specific embodiment of the present invention, the stiffness matrix of each planar isoparametric element in the global coordinate system is constructed based on the coordinate position, displacement and nodal force of each planar isoparametric element in the global coordinate system and the coordinate position of each planar isoparametric element in the local coordinate system, specifically including:

[0066] Step S2.1: Based on the coordinate positions of each plane isoparametric element in the global coordinate system and the local coordinate system, obtain the rotation matrix between the global coordinate system and the local coordinate system. The calculation formula of the rotation matrix is:

[0067]

[0068] Among them, R represents the rotation matrix between the global coordinate system and the local coordinate system, represents the coordinate position of the vertex of the plane isoparametric unit in the local coordinate system, and (x, y, z) represents the coordinate position of the vertex of the plane isoparametric unit in the global coordinate system.

[0069] Step S2.2: Based on the rotation matrix and the displacement of the vertex of each planar isoparametric unit in the global coordinate system, the displacement of the vertex of each planar isoparametric unit in the local coordinate system is obtained, and then the displacement matrix of each planar isoparametric unit in the local coordinate system is obtained.

[0070] Taking the plane isoparametric element as a plane quadrilateral isoparametric element as an example, the calculation formula for the displacement of the vertex of each plane isoparametric element in the local coordinate system is:

[0071]

[0072] in, represents the displacement of the i-th vertex of the planar quadrilateral isoparametric element in the local coordinate system, represents the displacement of the i-th vertex of the plane quadrilateral isoparametric element in the local coordinate system, d x ,d y ,d z represents the displacement of the i-th vertex of the plane quadrilateral isoparametric unit in each coordinate axis of the global coordinate system, d i Represents the displacement of the i-th vertex of the planar quadrilateral isoparametric element in the global coordinate system.

[0073] Displacement matrix of plane quadrilateral isoparametric element in local coordinate system It can be expressed as:

[0074]

[0075] Step S2.3: Based on the rotation matrix and the nodal forces of the vertices of each planar isoparametric unit in the global coordinate system, the nodal forces of the vertices of each planar isoparametric unit in the local coordinate system are obtained, and then the nodal force matrix of each planar isoparametric unit in the local coordinate system is obtained.

[0076] Taking the plane isoparametric element as a plane quadrilateral isoparametric element as an example, the calculation formula of the nodal force of each plane isoparametric element vertex in the local coordinate system is:

[0077]

[0078] in, represents the nodal force of the i-th vertex of the planar quadrilateral isoparametric element in the local coordinate system, represents the nodal force of the i-th vertex of the plane quadrilateral isoparametric element in the local coordinate system, x ,l y ,l z represents the nodal force of the i-th vertex of the plane quadrilateral isoparametric element in each coordinate axis of the global coordinate system, l i Represents the nodal force at the i-th vertex of the planar quadrilateral isoparametric element in the global coordinate system.

[0079] Nodal force matrix of plane quadrilateral isoparametric element in local coordinate system It can be expressed as:

[0080]

[0081] Step S2.4: Based on the displacement matrix and nodal force matrix of each planar isoparametric element in the local coordinate system, the stiffness matrix of each planar isoparametric element in the local coordinate system is obtained. The specific formula is:

[0082]

[0083] in, represents the stiffness matrix of the plane isoparametric element in the local coordinate system, D represents the displacement matrix of the plane isoparametric element in the global coordinate system, L represents the nodal force matrix of the plane isoparametric element in the global coordinate system, and C represents the transformation matrix composed of the rotation matrix R. When the plane isoparametric element is a plane quadrilateral isoparametric element, the transformation matrix is:

[0084]

[0085] Step S2.5: Obtain the stiffness matrix of each planar isoparametric element in the global coordinate system based on the rotation matrix and the stiffness matrix of each planar isoparametric element in the local coordinate system.

[0086] According to formula (6) to formula (8), we can get:

[0087]

[0088] Since C is an orthogonal matrix, C -1 =C T According to formula (11) and SD = L, the stiffness matrix of each plane isoparametric element in the global coordinate system is:

[0089]

[0090] Among them, S represents the stiffness matrix of the plane isoparametric element in the global coordinate system, C represents the transformation matrix composed of the rotation matrix R, Represents the stiffness matrix of the planar isoparametric element in the local coordinate system, and the superscript T represents the transpose.

[0091] Step S3: Combining the stiffness matrix of each plane isoparametric element in the global coordinate system and the stiffness matrix of each spatial polyhedron isoparametric element in the global coordinate system to obtain the stiffness matrix of the reservoir basin finite element model.

[0092] Because the spatial polyhedron isoparametric element is based on the global coordinate system, the stiffness matrix of the spatial polyhedron isoparametric element does not require conversion between the local coordinate system and the global coordinate system. The stiffness matrix of the spatial polyhedron isoparametric element is calculated based on the material properties of the spatial polyhedron isoparametric element. The specific calculation process is based on existing technology.

[0093] According to the vertex numbering order of each plane isoparametric element and each spatial polyhedron isoparametric element, the stiffness matrix of each plane isoparametric element in the global coordinate system is combined with the stiffness matrix of each spatial polyhedron isoparametric element in the global coordinate system to obtain the stiffness matrix of the reservoir basin finite element model. By establishing the reservoir basin finite element model and its stiffness matrix, the theoretical and numerical foundations for the subsequent assembly of the composite geomembrane model of the reservoir basin anti-seepage face dam are laid.

[0094] Step S4: Construct a finite element model of the dam under reservoir basin geological conditions in a global coordinate system.

[0095] A 3D model of the dam is constructed based on the geometric shape of the dam under the geological conditions of the reservoir basin, and the 3D model of the dam is meshed to obtain a finite element model of the dam. Figure 3 and Figure 4 As shown in FIG, the finite element model of the dam includes the main structures of the dam body, such as the cushion layer 1, transition layer, main rockfill 3, downstream rockfill 4, pressure slope 5, as well as anti-seepage structures such as the face plate 2, geomembrane 6, and joint waterstop.

[0096] Step S5: combining the reservoir basin finite element model with the dam finite element model to obtain a composite geomembrane model of the reservoir basin anti-seepage face dam.

[0097] Step S6: combining the stiffness matrix of the reservoir basin finite element model with the stiffness matrix of the dam finite element model to obtain the stiffness matrix of the composite geomembrane model of the reservoir basin anti-seepage face plate dam.

[0098] The material properties of each node unit in the dam finite element model are set, and the stiffness matrix is calculated according to the material properties of each node unit, thereby obtaining the stiffness matrix of the dam finite element model.

[0099] Step S7: Numerical calculation is performed based on the composite geomembrane model of the reservoir basin anti-seepage face dam and its stiffness matrix.

[0100] In a specific embodiment of the present invention, numerical calculation is performed based on the composite geomembrane model of the reservoir anti-seepage panel dam and its stiffness matrix, specifically including:

[0101] Step S7.1: Set boundary conditions and material properties, and apply loads (such as vertical loads, water loads, and seismic loads) to the geomembrane composite model of the reservoir basin anti-seepage face dam.

[0102] Step S7.2: Calculate the displacement distribution of the composite geomembrane model of the reservoir basin anti-seepage face plate dam based on the load and the stiffness matrix of the composite geomembrane model of the reservoir basin anti-seepage face plate dam;

[0103] Step S7.3: Calculate the stress distribution of the composite geomembrane model of the reservoir basin anti-seepage panel dam according to the displacement distribution of the composite geomembrane model of the reservoir basin anti-seepage panel dam, that is, obtain the displacement and stress distribution of each node in the composite geomembrane model of the reservoir basin anti-seepage panel dam, and simulate the displacement and stress distribution of the composite geomembrane under different loads and boundary conditions, such as Figures 5 to 10 shown. Figures 5 to 10 In the global coordinate system, the direction of the water flow is the positive direction of the X axis, the direction from the right bank to the left bank along the dam axis is the positive direction of the Y axis, and the vertical direction is the positive direction of the Z axis. Figures 5 to 10 It can be seen that the geomembrane has no obvious effect on the settlement of the geomembrane at the bottom of the reservoir basin under the action of water load, which can reflect the ability of the geomembrane to adapt to the settlement deformation of the reservoir bottom; compared with the displacement of the reservoir bottom along the river and across the river without geomembrane, the calculation results using the geomembrane unit are smaller and the displacement distribution is more uniform, reflecting the tensile capacity of the geomembrane at the bottom of the reservoir; there is a large uneven deformation between the backfill and excavation areas of the reservoir bottom, which causes the geomembrane to be more obviously tensile in the boundary area between excavation and backfill; with the design allowable tensile strain of the geomembrane as the control indicator, the calculation results have a large safety margin.

[0104] The calculation results were verified during the calculation process. Specifically, the numerical simulation results were compared and analyzed with the actual engineering monitoring data. By adjusting the geomembrane parameters and verification, the design rationality and engineering applicability of the composite geomembrane were improved.

[0105] The method of the present invention can be used to optimize the design parameters of the composite geomembrane, including membrane thickness, contact interface performance and installation process, so as to improve the safety and economy of the project.

[0106] Example 2

[0107] An embodiment of the present invention also provides an electronic device, which includes: a memory, a processor and a computer program / instructions stored in the memory, and the processor executes the computer program / instructions to implement the numerical simulation method of the composite geomembrane of the reservoir basin anti-seepage panel dam in the embodiment of the present application.

[0108] Although not shown, the electronic device includes a processor that can perform various appropriate operations and processes based on programs and / or data stored in a read-only memory (ROM) or programs and / or data loaded from a storage portion into a random access memory (RAM). The processor can be a multi-core processor or can include multiple processors. In some embodiments, the processor can include a general-purpose main processor and one or more special coprocessors, such as a central processing unit, a graphics processing unit (GPU), a neural network processor (NPU), a digital signal processor (DSP), etc. Various programs and data required for device operation are also stored in RAM. The processor, ROM, and RAM are connected to each other via a bus. An input / output (I / O) interface is also connected to the bus.

[0109] The processor and memory are used together to execute the program / instructions stored in the memory. When the program / instructions are executed by the computer, the methods, steps or functions described in the above embodiments can be implemented.

[0110] Although not shown, an embodiment of the present invention further provides a computer-readable storage medium having a computer program / instruction stored thereon, which, when executed by a processor, implements the numerical simulation method for the composite geomembrane of the reservoir basin anti-seepage panel dam in the embodiment of the present application.

[0111] Computer-readable storage media include permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. The information can be computer-readable instructions, data structures, program modules, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disc read-only memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassettes, magnetic disk storage or other magnetic storage devices, or any other non-transmission media that can be used to store information that can be accessed by a computing device. As defined herein, computer-readable media does not include transitory media such as modulated data signals and carrier waves.

[0112] Although not shown, an embodiment of the present invention further provides a computer program product, including: a computer program / instruction, which, when executed by a processor, implements the numerical simulation method for the composite geomembrane of the reservoir basin anti-seepage panel dam in the embodiment of the present application.

[0113] The above disclosure is only a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any technician familiar with this technical field can easily think of changes or modifications within the technical scope disclosed in the present invention, and they should all be covered by the scope of protection of the present invention.

Claims

1. A numerical simulation method for composite geomembrane of reservoir basin anti-seepage panel dam, characterized in that: The numerical simulation method comprises: Constructing a geomembrane plane model in a global coordinate system, dividing the geomembrane plane model into a plurality of plane isoparametric units, and constructing a corresponding spatial polyhedron isoparametric unit with each plane isoparametric unit as the top surface to obtain a reservoir basin finite element model; A local coordinate system is constructed based on each planar isoparametric element, and the stiffness matrix of each planar isoparametric element in the global coordinate system is constructed according to the coordinate position, displacement and nodal force of each planar isoparametric element in the global coordinate system and the coordinate position of each planar isoparametric element in the local coordinate system; The stiffness matrix of the reservoir basin finite element model is obtained by combining the stiffness matrix of each plane isoparametric element in the global coordinate system and the stiffness matrix of each spatial polyhedron isoparametric element in the global coordinate system. Construct the finite element model of the dam under reservoir basin geological conditions in the global coordinate system; The reservoir basin finite element model is combined with the dam finite element model to obtain a composite geomembrane model of the reservoir basin anti-seepage face dam; Combining the stiffness matrix of the reservoir basin finite element model with the stiffness matrix of the dam finite element model to obtain the stiffness matrix of the composite geomembrane model of the reservoir basin anti-seepage face dam; Numerical calculations are performed based on the composite geomembrane model of the reservoir basin anti-seepage panel dam and its stiffness matrix.

2. The numerical simulation method for composite geomembrane of reservoir basin anti-seepage panel dam according to claim 1 is characterized in that: The plane isoparametric unit is a plane quadrilateral isoparametric unit or a plane triangle isoparametric unit.

3. The numerical simulation method for composite geomembrane of reservoir basin anti-seepage panel dam according to claim 1 is characterized in that: The specific construction process of the stiffness matrix of each plane isoparametric element in the global coordinate system includes: According to the coordinate position of each plane isoparametric element in the global coordinate system and the local coordinate system, the rotation matrix between the global coordinate system and the local coordinate system is obtained; According to the rotation matrix and the displacement of the vertex of each plane isoparametric unit in the global coordinate system, the displacement of the vertex of each plane isoparametric unit in the local coordinate system is obtained, and then the displacement matrix of each plane isoparametric unit in the local coordinate system is obtained; According to the rotation matrix and the nodal force of the vertex of each plane isoparametric unit in the global coordinate system, the nodal force of the vertex of each plane isoparametric unit in the local coordinate system is obtained, and then the nodal force matrix of each plane isoparametric unit in the local coordinate system is obtained; According to the displacement matrix and nodal force matrix of each plane isoparametric element in the local coordinate system, the stiffness matrix of each plane isoparametric element in the local coordinate system is obtained; According to the rotation matrix and the stiffness matrix of each planar isoparametric element in the local coordinate system, the stiffness matrix of each planar isoparametric element in the global coordinate system is obtained.

4. The numerical simulation method for composite geomembrane of reservoir basin anti-seepage panel dam according to claim 3 is characterized in that: The calculation formula of the rotation matrix is: Among them, R represents the rotation matrix between the global coordinate system and the local coordinate system, represents the coordinate position of the vertex of the plane isoparametric unit in the local coordinate system, and (x, y, z) represents the coordinate position of the vertex of the plane isoparametric unit in the global coordinate system.

5. The numerical simulation method for composite geomembrane of reservoir basin anti-seepage panel dam according to claim 3 is characterized in that: The specific formula of the stiffness matrix of each plane isoparametric element in the local coordinate system is: in, represents the stiffness matrix of the plane isoparametric element in the local coordinate system, represents the displacement matrix of the plane isoparametric element in the local coordinate system, represents the nodal force matrix of the plane isoparametric element in the local coordinate system, D represents the displacement matrix of the plane isoparametric element in the global coordinate system, L represents the nodal force matrix of the plane isoparametric element in the global coordinate system, and C represents the transformation matrix composed of the rotation matrix R.

6. The numerical simulation method for composite geomembrane of reservoir basin anti-seepage panel dam according to claim 3 is characterized in that: The specific formula of the stiffness matrix of each plane isoparametric element in the global coordinate system is: Among them, S represents the stiffness matrix of the plane isoparametric element in the global coordinate system, C represents the transformation matrix composed of the rotation matrix R, Represents the stiffness matrix of the planar isoparametric element in the local coordinate system, and the superscript T represents the transpose.

7. The numerical simulation method for composite geomembrane of reservoir basin anti-seepage face dam according to any one of claims 1 to 6, characterized in that: Numerical calculations are performed based on the composite geomembrane model of the reservoir basin anti-seepage face dam and its stiffness matrix, specifically including: Set boundary conditions and material properties, and apply loads to the geomembrane composite model of the reservoir anti-seepage face dam. The displacement distribution of the composite geomembrane model of the reservoir basin anti-seepage face dam is calculated based on the load and the stiffness matrix of the composite geomembrane model of the reservoir basin anti-seepage face dam; The stress distribution of the composite geomembrane model of the reservoir basin anti-seepage face plate dam is calculated based on the displacement distribution of the composite geomembrane model of the reservoir basin anti-seepage face plate dam.

8. An electronic device comprising a memory, a processor, and a computer program / instruction stored in the memory, characterized in that: The processor executes the computer program / instruction to implement the numerical simulation method for composite geomembrane of reservoir basin anti-seepage face dam according to any one of claims 1 to 7.

9. A computer-readable storage medium having a computer program / instruction stored thereon, characterized in that: When the computer program / instruction is executed by a processor, the numerical simulation method for composite geomembrane of reservoir basin anti-seepage face dam according to any one of claims 1 to 7 is realized.

10. A computer program product comprising a computer program / instructions, characterized in that When the computer program / instruction is executed by a processor, the numerical simulation method for composite geomembrane of reservoir basin anti-seepage face dam according to any one of claims 1 to 7 is realized.

Images