Three-dimensional stress deformation simulation analysis method and system for geomembrane structure

Abstract

The invention provides a geomembrane structure-oriented three-dimensional stress deformation simulation analysis method and system, and relates to the technical field of geomembrane structures, and the method comprises the following steps: carrying out parametric modeling on a geomembrane structure; identifying an interlayer contact surface according to the geomembrane layered modeling data, and performing discretization processing according to a preset grid to obtain a limited number of contact surface units; dividing a limited number of contact surface units into direct contact surface units and indirect contact surface units; stress deformation load data are collected for simulation to obtain three-dimensional stress deformation direct simulation data and three-dimensional stress deformation indirect simulation data; and obtaining three-dimensional stress deformation fusion simulation data of the geomembrane structure through data fusion. The technical problem that the precision of a three-dimensional stress deformation simulation result is insufficient due to the fact that the whole interlayer contact surface is regarded as a homogeneous body and the interlayer contact defect is ignored in the prior art is solved, and the three-dimensional stress deformation simulation precision of the geomembrane structure is improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of geomembrane structure, and particularly relates to a three-dimensional stress deformation simulation analysis method and system for a geomembrane structure. BACKGROUND

[0002] The existing simulation method usually simplifies the entire interlayer contact surface as a homogeneous body with uniform mechanical properties, ignoring the non-uniform distribution of interface properties in space caused by various factors in actual engineering. Under this simplification, the interaction between layers is often described as a uniform contact model using uniform contact parameters. However, in actual engineering applications, the contact state between the geomembrane and the overlying protective layer and the underlying supporting layer will change significantly due to the action of various factors, resulting in non-uniformity of the mechanical properties between the membrane layers, affecting the stress distribution and deformation behavior of the membrane layers, and failing to accurately simulate the three-dimensional stress deformation of the geomembrane structure.

[0003] In summary, the prior art has the technical problem of insufficient accuracy of three-dimensional stress deformation simulation results due to the entire interlayer contact surface being regarded as a homogeneous body and ignoring the interlayer contact defects. SUMMARY

[0004] The purpose of the present application is to provide a three-dimensional stress deformation simulation analysis method and system for a geomembrane structure, to solve the technical problem of insufficient accuracy of three-dimensional stress deformation simulation results in the prior art due to the entire interlayer contact surface being regarded as a homogeneous body and ignoring the interlayer contact defects.

[0005] In order to achieve the above purpose, the present application provides a three-dimensional stress deformation simulation analysis method and system for a geomembrane structure.

[0006] In a first aspect, the application provides a three-dimensional stress deformation simulation analysis method for a geomembrane structure, which is implemented by a three-dimensional stress deformation simulation analysis system for a geomembrane structure. The three-dimensional stress deformation simulation analysis method for a geomembrane structure includes: parameterized modeling of a geomembrane structure to construct geomembrane layered modeling data; identifying interlayer contact surfaces according to the geomembrane layered modeling data, and discretizing the interlayer contact surfaces according to a preset grid to obtain a finite number of contact surface elements; collecting interlayer contact defects, interlayer construction quality, and environmental pollution states corresponding to the finite number of contact surface elements, and dividing the finite number of contact surface elements into direct contact surface elements and indirect contact surface elements by taking the interlayer contact defects, the interlayer construction quality, and the environmental pollution states as hierarchical features; collecting stress deformation load data to simulate the direct contact surface elements and the indirect contact surface elements to obtain three-dimensional stress deformation direct simulation data and three-dimensional stress deformation indirect simulation data; and fusing the three-dimensional stress deformation direct simulation data and the three-dimensional stress deformation indirect simulation data to obtain three-dimensional stress deformation fusion simulation data of the geomembrane structure.

[0007] Optionally, stress contact areas and non-stress contact areas based on the interlayer contact surfaces under the stress deformation load data are identified; the stress contact areas are discretized using a first preset grid to obtain a first group of contact surface elements, and the non-stress contact areas are discretized using a second preset grid to obtain a second group of contact surface elements, wherein the step length of the first preset grid is smaller than that of the second preset grid; and a finite number of contact surface elements are obtained according to the first group of contact surface elements and the second group of contact surface elements.

[0008] Optionally, a stress transition contact area based on the interlayer contact surfaces under the stress deformation load data is identified, and the stress transition contact area is a transition contact area between the stress contact areas and the non-stress contact areas; the stress transition contact area is discretized using a third preset grid to obtain a third group of contact surface elements, wherein the step length of the third preset grid is greater than that of the first preset grid and smaller than that of the second preset grid; and a finite number of contact surface elements are obtained according to the first group of contact surface elements, the second group of contact surface elements, and the third group of contact surface elements.

[0009] Optionally, the contribution degree score of each grading feature of the limited number of contact surface units with respect to the interlayer contact defect, the interlayer construction quality, and the environmental pollution state is calculated respectively to obtain a defect contribution degree score, a construction contribution degree score, and a pollution contribution degree score; the defect contribution degree score, the construction contribution degree score, and the pollution contribution degree score are weighted and calculated by configuring a contribution degree weight to obtain a contribution degree comprehensive score of the limited number of contact surface units; the contact surface units less than a preset score threshold are divided into direct contact surface units; the contact surface units greater than or equal to the preset score threshold are divided into indirect contact surface units.

[0010] Optionally, equivalent thin layer mechanical parameters of the indirect contact surface units are constructed; stress variation influence degrees based on the equivalent thin layer mechanical parameters are simulated according to the stress deformation load data, the stress variation influence degrees being data variation degrees between equivalent simulation stress deformation output data and stress deformation load data; the indirect contact surface units with a stress variation influence degree less than a preset influence degree are re-labeled as direct contact surface units.

[0011] Optionally, a stress deformation input vector set is constructed based on the stress deformation load data and state input data of the direct contact surface units, including an external load vector, a material parameter vector, a boundary condition vector, and a time step parameter vector; a direct contact stiffness model is constructed through normal contact stiffness, normal projection matrix, tangential contact stiffness, and tangential projection matrix; three-dimensional finite element solving is performed on the stress deformation input vector set according to the direct contact stiffness model by using Coulomb friction constraint to obtain three-dimensional stress deformation direct simulation data, including stress displacement increment, interlayer normal stress, interlayer shear stress, and stress field distribution.

[0012] Optionally, an indirect contact weak coupling model is constructed, the indirect contact weak coupling model including a stress attenuation coefficient; a material degradation coefficient is introduced to identify an equivalent material matrix of the indirect contact surface units; three-dimensional finite element solving is performed on the stress deformation input vector set according to the stress attenuation coefficient of the indirect contact weak coupling model and the equivalent material matrix to obtain three-dimensional stress deformation indirect simulation data.

[0013] Optionally, the indirect contact weak coupling model includes a stress attenuation coefficient, the stress attenuation coefficient being obtained by fitting and calculating a defect attenuation factor, a quality attenuation factor, and a pollution attenuation factor.

[0014] Optionally, a direct-fusion smoothing weight, an indirect-fusion smoothing weight, and a hybrid-fusion smoothing weight are defined; the direct contact surface units and the indirect contact surface units are taken as fusion nodes, and fusion is performed according to the direct-fusion smoothing weight, the indirect-fusion smoothing weight, and the hybrid-fusion smoothing weight to obtain three-dimensional stress deformation fusion simulation data.

[0015] In a second aspect, the application further provides a three-dimensional stress deformation simulation analysis system for geomembrane structure, used for performing the three-dimensional stress deformation simulation analysis method for geomembrane structure as described in the first aspect, wherein the three-dimensional stress deformation simulation analysis system for geomembrane structure comprises: a parameterized modeling module, used for performing parameterized modeling on the geomembrane structure to construct geomembrane layered modeling data; a discretization processing module, used for identifying interlayer contact surfaces according to the geomembrane layered modeling data, and performing discretization processing on the interlayer contact surfaces according to a preset grid to obtain a finite number of contact surface elements; a contact surface element division module, used for collecting interlayer contact defects, interlayer construction quality and environmental pollution states corresponding to the finite number of contact surface elements, and dividing the finite number of contact surface elements into direct contact surface elements and indirect contact surface elements by taking the interlayer contact defects, the interlayer construction quality and the environmental pollution states as hierarchical features; a simulation analysis module, used for collecting stress deformation load data to simulate the direct contact surface elements and the indirect contact surface elements to obtain three-dimensional stress deformation direct simulation data and three-dimensional stress deformation indirect simulation data; and a data fusion module, used for fusing the three-dimensional stress deformation direct simulation data and the three-dimensional stress deformation indirect simulation data to obtain three-dimensional stress deformation fusion simulation data of the geomembrane structure.

[0016] The one or more technical solutions provided in the application have at least the following technical effects or advantages:

[0017] The geotechnical membrane layered modeling data is constructed by parameterized modeling of the geotechnical membrane structure; the interlayer contact surface is identified according to the geotechnical membrane layered modeling data, and the interlayer contact surface is discretized by a preset grid to obtain a limited number of contact surface units; the interlayer contact defects, the interlayer construction quality and the environmental pollution state corresponding to the limited number of contact surface units are collected, and the interlayer contact defects, the interlayer construction quality and the environmental pollution state are taken as grading features to divide the limited number of contact surface units into direct contact surface units and indirect contact surface units; stress deformation load data is collected to simulate the direct contact surface units and the indirect contact surface units to obtain three-dimensional stress deformation direct simulation data and three-dimensional stress deformation indirect simulation data; and the three-dimensional stress deformation direct simulation data and the three-dimensional stress deformation indirect simulation data are fused to obtain three-dimensional stress deformation fusion simulation data of the geotechnical membrane structure. That is, the geotechnical membrane layered modeling data is constructed by parameterized modeling, the interlayer contact surface is discretized by a preset grid, the contact units or interface units are used to simulate the complex interaction between the geotechnical membrane and the overlying protective layer and the underlying supporting layer, such as friction, sliding and voiding, the interlayer contact defects, the construction quality and the environmental pollution state are collected as grading standards, the contact surface units are divided into direct contact surface units and indirect contact surface units, and when the contact stress deformation simulation is performed, the direct and indirect fusion method is used to analyze the deformation simulation state between the layers, thereby improving the data accuracy of the simulation analysis.

[0018] The above description is only a summary of the technical solutions of the present application. In order to more clearly understand the technical means of the present application, the specific embodiments of the present application can be implemented according to the content of the specification, and in order to make the above and other purposes, features and advantages of the present application more obvious and easy to understand, the following specific embodiments of the present application are described. It should be understood that the content described in this part is not intended to identify the key or important features of the embodiments of the present application, nor is it intended to limit the scope of the present application. Other features of the present application will become apparent from the following description. BRIEF DESCRIPTION OF DRAWINGS

[0019] In order to more clearly illustrate the technical solutions in the present application or the prior art, the following will briefly introduce the drawings needed to be used in the embodiments or the prior art description. Obviously, the drawings in the following description are only exemplary, and other drawings can be obtained by the provided drawings without creative labor for those skilled in the art.

[0020] Figure 1 The flowchart of the three-dimensional stress deformation simulation analysis method for the geotechnical membrane structure of the present application.

[0021] Figure 2 The structural diagram of the three-dimensional stress deformation simulation analysis system for the geotechnical membrane structure of the present application.

[0022] Reference signs: parameterized modeling module 11, discretization processing module 12, contact surface unit division module 13, simulation analysis module 14, data fusion module 15. DETAILED DESCRIPTION

[0023] The present application provides a three-dimensional stress deformation simulation analysis method and system for geomembrane structure, which solves the technical problem of insufficient precision of three-dimensional stress deformation simulation results in the prior art due to the entire interlayer contact surface being regarded as a homogeneous body and the interlayer contact defects being ignored. The present application constructs geomembrane layered modeling data through parameterized modeling, discretizes the interlayer contact surface through a preset grid, simulates the complex interaction such as friction, sliding and voiding between the geomembrane and the overlying protective layer and the underlying supporting layer by using contact elements or interface elements, collects the interlayer contact defects, construction quality and environmental pollution state as grading standards, divides the contact surface units into direct contact surface units and indirect contact surface units, and adopts a direct and indirect fusion method to analyze the deformation simulation state between layers when simulating the contact stress deformation, thereby improving the data accuracy of simulation analysis.

[0024] The technical solutions in the present application will be described clearly and completely below with reference to the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present application, rather than all the embodiments of the present application. It should be understood that the present application is not limited by the example embodiments described herein. Based on the embodiments of the present application, all other embodiments obtained by those of ordinary skill in the art without creative work fall within the scope of protection of the present application. In addition, it should be noted that, for convenience of description, only parts related to the present application are shown in the accompanying drawings, rather than all parts.

[0025] Embodiment one, please refer to the accompanying Figure 1 The present application provides a three-dimensional stress deformation simulation analysis method for geomembrane structure, wherein the three-dimensional stress deformation simulation analysis method for geomembrane structure is applied to a three-dimensional stress deformation simulation analysis system for geomembrane structure, and the three-dimensional stress deformation simulation analysis method for geomembrane structure specifically includes the following steps:

[0026] Parameterized modeling is performed on the geomembrane structure to construct geomembrane layered modeling data.

[0027] Specifically, the geomembrane structure is a multi-layer composite system composed of a high-performance polymer film as the core impermeable layer, and upper and lower protective layers, drainage layers, support layers, etc., commonly used in landfill projects. The geomembrane structure generally refers to a structural system composed of multiple layers of geomembranes, each layer of membrane having a certain thickness, mechanical properties and impermeability. According to the engineering design drawings and geological survey reports, the entire geomembrane structure is abstracted into several continuous layers, such as the foundation, the lower protective layer, the geomembrane impermeable layer, the upper protective layer, the drainage layer and the landfill, etc. By analyzing the physical properties of the geomembrane, the geometric dimensions and mechanical properties of each layer of membrane are extracted, and parameters are set for each layer of membrane according to the actual needs of the project, including the thickness of the membrane and the physical properties of the materials used. The parameters of each layer of membrane are organized into a structured data set, and the properties of each layer of membrane are saved in the form of a data table.

[0028] A parameterized modeling software is used to create a three-dimensional geometric model for each layer driven by key parameters. During the modeling process of the geomembrane structure, the data of each layer of membrane is superimposed to build a complete layered modeling data system. The data of each layer of membrane not only includes the geometric parameters of the membrane, but also includes the contact surface information between the membrane and other layers of membrane, such as the friction coefficient of the membrane and the contact area. Finally, a three-dimensional digital model containing all the layered structures is automatically generated, and the construction of the geomembrane layered modeling data containing the geometric information, topological relationship and initial material properties of each layer is output. The geomembrane layered modeling data is a set of digital information that can completely describe the spatial distribution, geometric dimensions, material properties and mutual position relationship of each layer of the geomembrane structure.

[0029] By parameterizing the modeling of the geomembrane structure and constructing the layered modeling data, the geometric characteristics and physical properties of each layer of membrane are clearly defined, ensuring that the simulation process accurately reflects the geomembrane structure in the actual project. When the design scheme changes, only a few parameters need to be modified, and the model can be automatically updated within seconds, avoiding the huge workload of starting over in traditional methods.

[0030] According to the geomembrane layered modeling data, the interlayer contact surface is identified, and the interlayer contact surface is discretized to obtain a finite number of contact surface elements according to a preset grid.

[0031] Further, the application further includes the following steps: identifying the stress contact area and the non-stress contact area based on the interlayer contact surface under the stress deformation load data; discretizing the stress contact area to obtain a first group of contact surface elements using a first preset grid, and discretizing the non-stress contact area to obtain a second group of contact surface elements using a second preset grid, wherein the step size of the first preset grid is smaller than the step size of the second preset grid; and obtaining a finite number of contact surface elements according to the first group of contact surface elements and the second group of contact surface elements.

[0032] Specifically, according to the geomembrane layer modeling data, all the interlayer contact surfaces that need to be analyzed are automatically identified, that is, through the analysis of the geometric position and physical contact state between the membrane layers, it is clear which areas are the actual contact surfaces. The interlayer contact surface is the contact area between different membrane layers in the geomembrane structure, which determines the mechanical interaction between the layers, such as friction, shear force, etc.

[0033] The stress deformation load data is imported, that is, the stress and deformation data of the geomembrane structure under the action of external load. According to the stress deformation load data, the identified interlayer contact surface is divided into stress contact area and non-stress contact area. The stress contact area is the area between the layers where there is obvious contact stress or deformation, and the stress in these areas is usually large, requiring high calculation accuracy, such as slopes, anchoring trenches, and load concentration parts; while the non-stress contact area refers to the area where there is no significant contact stress or deformation between the layers, and the calculation accuracy requirement of these areas is relatively low, such as large flat areas at the bottom of the field. That is, it is identified that under the stress deformation load data, which contact surface area is the key area in mechanics, that is, the stress contact area, and which is the secondary area, that is, the non-stress contact area.

[0034] In finite element analysis, meshing is the process of dividing the entire model space into small elements for numerical solution. The preset mesh refers to the meshing standard set before the analysis, including the size of the mesh, the step length, etc., which is usually optimized according to the importance of the area and the stress distribution, including the first preset mesh, the second preset mesh, and the third preset mesh, which correspond to the stress contact area, the non-stress contact area, and the stress transition contact area for discretization processing, respectively.

[0035] The first preset mesh is applied, that is, a smaller mesh step is used to divide the stress contact area, that is, a large number of small and closely arranged elements are generated in the stress contact area, forming the first group of contact surface elements, which can depict the complex stress gradient changes and potential strain localization phenomena in this area. For the non-stress contact area identified as secondary, the second preset mesh is applied, that is, a larger mesh step is used for division, generating the second group of contact surface elements, which are small in number and large in size, sufficient to describe the overall mechanical behavior of the non-stress contact area, but will not cause excessive burden in calculation. The step length of the first preset mesh is smaller than that of the second preset mesh.

[0036] Exemplarily, assuming that in a simulation of a geomembrane impermeable layer of a tailings dam body with a length of 100m and a slope gradient of 1:1.5, the stress state under water pressure during impoundment needs to be analyzed. The stress contact area is identified as the entire dam slope, with an area of about 850m 2, the first preset grid is adopted, and the step length is set to 0.15 m. After discretization, about 37800 fine quadrilateral contact surface units, i.e., the first group of contact surface units, are generated. The non-stress contact area is identified as the flat area at the bottom of the library, with an area of about 2200 m 2 , the second preset grid is adopted, and the step length is set to 0.8 m. After discretization, about 3400 rough quadrilateral contact surface units, i.e., the second group of contact surface units, are generated.

[0037] Through the discretization process, the first group of contact surface units and the second group of contact surface units are obtained, the number and shape of which are directly related to the stress distribution of the contact surface area and the grid step length. According to the grid division of the stress contact area and the non-stress contact area, a limited set of contact surface units is finally generated. Fine grid division is performed for different stress areas and contact characteristics, which improves the calculation accuracy in key areas and reduces the calculation complexity in non-stress areas, thereby reducing the waste of computing resources while ensuring simulation accuracy.

[0038] Further, the application further includes the following steps: identifying a stress transition contact area based on the interlayer contact surface under the stress deformation load data, the stress transition contact area being a transition contact area between the stress contact area and the non-stress contact area; discretizing the stress transition contact area using a third preset grid to obtain a third group of contact surface units, wherein the step length of the third preset grid is greater than the step length of the first preset grid and less than the step length of the second preset grid; and obtaining a limited number of contact surface units based on the first group of contact surface units, the second group of contact surface units and the third group of contact surface units.

[0039] Specifically, after identifying the stress contact area and the non-stress contact area, the stress transition area in the interlayer contact surface is identified through analysis of the stress deformation load data. The stress transition contact area refers to the area between the stress contact area and the non-stress contact area, in which the stress gradually decreases from a larger value. The stress transition contact area is characterized by a significant change in stress state, neither as high-precision as the core stress area nor as highly simplified as the non-stress area, but a transition zone from complex to simple mechanical behavior.

[0040] The third preset grid scheme is called, an intermediate step length greater than the first preset grid step length but smaller than the second preset grid step length is adopted to discretize the stress transition contact area, thereby generating the third group of contact surface units. The third group of contact surface units plays a role of bridging the gap in size, the fine grid unit size of the stress contact area is not greatly different, ensuring smooth transmission of stress; at the same time, the rough grid unit size of the non-stress contact area is relatively close, achieving a smooth decrease in calculation scale.

[0041] The first preset grid is used for the stress contact area, has a small step length, and is used for more accurate calculation of stress and deformation; the second preset grid is used for the non-stress contact area, has a large step length, and is used for optimization of calculation efficiency; and the third preset grid is used for the stress transition contact area, has a step length between the first preset grid and the second preset grid, ensures the accuracy of the transition area, and controls the calculation complexity.

[0042] Through the discretization processing of the stress contact area, the non-stress contact area, and the stress transition contact area, a first group of contact surface units, a second group of contact surface units, and a third group of contact surface units are obtained, the first group of contact surface units, the second group of contact surface units, and the third group of contact surface units are seamlessly spliced in space, and collectively constitute a finite set of contact surface units with continuously changing grid density. For example, the stress transition contact area is located between the stress contact area and the non-stress contact area, is a ring belt with a width of about 2 meters, and has an area of about 1500 m 2 . The third preset grid is used, the step length is 0.4 m, and about 9400 units are generated.

[0043] Through the discretization processing of the stress transition contact area with a moderate grid step length, the calculation efficiency is improved while the calculation accuracy is ensured. The refined modeling of the stress transition area enables the simulation results to better reflect the mechanical state in actual engineering, especially the stress distribution in the transition area. Through accurate division of the contact surface units, not only the accuracy of the simulation is improved, but also unnecessary calculation burden is reduced, and the efficiency of the simulation process is improved.

[0044] The interlayer contact defects, the interlayer construction quality, and the environmental pollution state corresponding to the finite contact surface units are collected, and the interlayer contact defects, the interlayer construction quality, and the environmental pollution state are used as grading features to divide the finite contact surface units into direct contact surface units and indirect contact surface units.

[0045] Further, the application further includes the following steps: calculating the contribution degree score of each grading feature in the interlayer contact defects, the interlayer construction quality, and the environmental pollution state of the finite contact surface units respectively, obtaining the defect contribution degree score, the construction contribution degree score, and the pollution contribution degree score; performing weighted calculation on the defect contribution degree score, the construction contribution degree score, and the pollution contribution degree score by configuring the contribution degree weight, obtaining the contribution degree comprehensive score of the finite contact surface units; dividing the contact surface units less than the preset score threshold into direct contact surface units; and dividing the contact surface units greater than or equal to the preset score threshold into indirect contact surface units.

[0046] Specifically, after obtaining a limited number of contact surface units, a multi-source data collection and fusion process is initiated to collect interlayer contact defects, interlayer construction quality, and environmental pollution status. Interlayer contact defects are defects that can exist between layers, such as uneven contact pressure, bubbles, impurity layers, or unevenness of the film layer, which can affect the friction and bearing capacity between the film layers, thereby affecting the mechanical properties and impermeability of the geomembrane structure. Interlayer construction quality is the degree of execution of design requirements and standard specifications during the laying and installation of geomembranes and their upper and lower protective layers, such as interface flatness, material compaction, joint quality, and interlayer adhesion. Environmental pollution status is the condition in which external substances invade the interlayer interface and affect the material properties, including physical, chemical, and biological pollution. This includes factors such as voids between film layers, uneven film surfaces, bubbles, and impurity layers.

[0047] Through detailed inspection and detection of the interlayer contact surface of the geomembrane, data on contact defects are collected, including factors such as voids between film layers, uneven film surfaces, bubbles, and impurity layers, and are located and associated with the specific contact surface unit they are in. Through field construction monitoring, construction reports, and quality inspections, quality data during the construction process of the film layers are collected, including construction temperature, construction speed, joint quality between film layers, and film layer integrity. Through environmental monitoring equipment such as gas sensors, temperature and humidity sensors, and chemical sensors, data on pollutant concentrations and temperature and humidity changes in the environment where the geomembrane is located are collected, reflecting whether the geomembrane is affected by environmental factors such as chemical corrosion and microbial pollution.

[0048] The contribution scores of each graded feature in the interlayer contact defects, interlayer construction quality, and environmental pollution status of the limited number of contact surface units are calculated respectively to obtain defect contribution scores, construction contribution scores, and pollution contribution scores. For interlayer contact defects, they are quantified according to their type and severity, such as no defects scoring 0 points, slight wrinkles scoring 8 points, and severe scratches or holes scoring 20 points, thereby obtaining the defect contribution scores. For interlayer construction quality, reverse scoring is performed according to indicators such as compaction and flatness, such as excellent areas with compaction ≥ 95% scoring 0 points, qualified areas with compaction about 90% scoring 6 points, and unqualified areas with compaction < 85% scoring 15 points, thereby obtaining the construction contribution scores. For environmental pollution status, scoring is performed according to the nature and degree of pollution, such as no pollution scoring 0 points, slight physical accumulation scoring 10 points, and strong chemical corrosive pollutants scoring 18 points, thereby obtaining the pollution contribution scores.

[0049] The contribution degree weight is configured based on the relative importance of the three grading characteristics, i.e., interlayer contact defects, interlayer construction quality, and environmental pollution state, in determining the final contact type. The configuration of the weight relies on engineering experience and statistical data. For example, at the initial stage of operation, the weight of construction quality may be relatively high; and after long-term operation, the weight of environmental pollution gradually increases.

[0050] The contribution degree weight is called to perform weighted calculation on the defect contribution score, the construction contribution score, and the pollution contribution score, and finally obtain the contribution degree comprehensive score of a limited number of contact surface units. The contact surface units with a contribution degree comprehensive score less than a preset score threshold are divided into direct contact surface units, and the contact surface units with a contribution degree comprehensive score greater than or equal to the preset score threshold are divided into indirect contact surface units. The preset score threshold is a demarcation value for dividing the continuous contribution degree comprehensive score into two categories of direct contact and indirect contact. For example, assuming that the weights are defect weight 0.4, construction quality weight 0.3, and pollution weight 0.3, and the preset score threshold is 10 points. Unit A is located in a flat area at the bottom of the reservoir, has no defects, and thus has a defect score of 0; has a compaction degree of 96% and a construction quality score of 0; has no pollution and a pollution score of 0; and thus has a comprehensive score of (0×0.4) + (0×0.3) + (0×0.3) = 0 points, and is divided into a direct contact surface unit. Unit B is located on a slope, has a slight wrinkle, and has a wrinkle with a height of 2 cm, a defect score of 8, a compaction degree of 92%, a construction quality score of 4, no pollution, and a pollution score of 0; and thus has a comprehensive score of (8×0.4) + (4×0.3) + (0×0.3) = 3.2 + 1.2 + 0 = 4.4 points, and is divided into a direct contact surface unit. Unit C is located at the bottom of the slope, has a poor environment, has a hot welding repair mark, a defect score of 12, a compaction degree of less than 88%, a construction quality score of 8, and a pollution score of 15 due to chemical crystallization caused by leachate precipitation; and thus has a comprehensive score of (12×0.4) + (8×0.3) + (15×0.3) = 4.8 + 2.4 + 4.5 = 11.7 points, and is divided into an indirect contact surface unit. By repeating the above operation on the remaining contact surface units, finally, 85% of the units are quantified as direct contact, and 15% of the units are quantified as indirect contact.

[0051] The direct contact surface unit is a contact area unit between the geomembrane and the adjacent material layer, which is tightly attached without intermediate medium. The indirect contact surface unit is a contact area unit between the geomembrane and the adjacent material layer, which has a third phase medium or physical spacing between the geomembrane and the adjacent material layer. The intermediate medium or spacing changes the original interface contact state, resulting in a fundamental change in mechanical behavior, such as impurity layer, loose layer, construction quality defect, etc.

[0052] By calculating the contribution score of each contact surface unit and combining the weight for weighted calculation, the influence of interlayer contact defects, construction quality and environmental pollution state on the contact surface unit is comprehensively considered, so as to classify the contact surface unit and accurately divide which contact surface unit needs high-precision simulation and which contact surface unit can be simplified.

[0053] Further, the application further includes the following steps: constructing equivalent thin layer mechanical parameters of the indirect contact surface unit; simulating stress change influence degree based on the equivalent thin layer mechanical parameters according to the stress deformation load data, the stress change influence degree being a data change degree between equivalent simulation obtained stress deformation output data and stress deformation load data; and re-labeling the indirect contact surface unit with a stress change influence degree less than a preset influence degree as a direct contact surface unit.

[0054] Specifically, the equivalent thin layer mechanical parameters of the indirect contact surface unit are constructed, the actual contact surface is converted into a thin layer with a certain thickness, and the equivalent mechanical parameters of the layer are assigned. For the unit with impurity layer, since the impurity layer is a physically existing material, the mechanical parameters are determined by indoor test method. Through fine sampling on site, the impurity layer sample in original state is sent to the laboratory for direct shear test according to relevant specifications. The direct shear test determines the shear strength of the sample under different vertical pressures, and draws the strength envelope to determine the cohesion and internal friction angle. For the unit with loose area, since the loose area is essentially a region with reduced density of the original material, the equivalent stiffness parameters are determined by indirect field test and back analysis method. Through geological radar data inversion, the difference in propagation speed of electromagnetic waves in different density materials is used, combined with a small amount of drilling calibration data, to back calculate the equivalent elastic modulus distribution of the region. For construction quality defects, the strength reduction coefficient is determined by non-destructive testing or sampling test. For welding defects, the defect size and type are determined by ultrasonic testing, and the reduction coefficient is determined according to the defect rate; for material damage, the tensile test is carried out by sampling, and the strength ratio of the damaged part to the intact material is compared to determine the reduction coefficient; for laying defects, the wrinkle geometric parameters are determined by field detection, and the reduction coefficient is determined according to the wrinkle height to wavelength ratio.

[0055] The constructed equivalent thin layer mechanical parameters are substituted into the numerical model, and special simulation analysis is carried out according to the stress deformation load data. In the simulation process, the stress change influence degree of each indirect contact surface unit is calculated, the stress change influence degree is defined by comparing the difference between the mechanical response calculated by the equivalent thin layer model and the original load data, and reflects the actual importance of the defect interface in the structure system. Smaller stress change influence degree means smaller difference between the simulation result and the actual data, and higher simulation accuracy.

[0056] The stress deformation load data is stress and deformation data generated by the geomembrane structure under the action of external load. The data variation between the stress deformation output data obtained by calculation simulation and the actual stress deformation load data is used to measure the accuracy of the simulation result. For a certain indirect contact surface element, the stress deformation data obtained by simulation is stress 10 kPa and displacement 5 mm, while the actual load data is stress 9.5 kPa and displacement 4.8 mm. Then the stress variation influence degree can be calculated by comparing the difference between the two sets of data. If the stress variation influence degree is 0.05, that is, the deviation between the simulation result and the actual data is less than 5%, it means that the simulation result is relatively accurate.

[0057] If the stress variation influence degree of a certain indirect contact surface element is less than the preset influence degree, it is considered that the stress transmission of the element is relatively stable and meets the condition of direct contact, and it is re-labeled as a direct contact element. The preset influence degree is a threshold value preset for judging whether the indirect contact surface element needs to be reclassified, which is usually determined based on engineering experience, safety standards and numerical tests. By constructing the equivalent thin layer mechanics parameters of the indirect contact surface element and calculating the stress variation influence degree, the simulation accuracy is effectively evaluated, and the reliability of the simulation result is ensured. If the stress variation influence degree of some indirect contact elements is small, it means that their mechanical behavior is close to that of direct contact surface elements, which accurately reflects the actual working condition. By re-labeling these elements as direct contact surface elements, the accuracy and stability of the simulation model are further improved.

[0058] The stress deformation load data is used to simulate the direct contact surface elements and the indirect contact surface elements to obtain three-dimensional stress deformation direct simulation data and three-dimensional stress deformation indirect simulation data.

[0059] Further, the application further includes the following steps: based on the stress deformation load data and the state input data of the direct contact surface element, a stress deformation input vector set is constructed, including an external load vector, a material parameter vector, a boundary condition vector and a time step parameter vector; a direct contact stiffness model is constructed through normal contact stiffness, normal projection matrix, tangential contact stiffness and tangential projection matrix; three-dimensional finite element solving is performed on the stress deformation input vector set according to the direct contact stiffness model by using Coulomb friction constraint, to obtain three-dimensional stress deformation direct simulation data, including stress displacement increment, interlayer normal stress, interlayer shear stress and stress field distribution.

[0060] Specifically, based on the stress deformation load data and the state input data of the screened direct contact surface elements, a stress deformation input vector set is constructed. For direct contact surface elements, state input data includes external load, material parameters, boundary conditions and time step parameters. The stress deformation input vector set is a vector set used to represent all input parameters affecting stress deformation calculation, including external load, material properties, boundary conditions and time step, etc. The external load vector is used to describe the external load applied to the geomembrane structure, such as pressure, tension or bending moment, which affects the stress state of the structure and causes deformation through these loads; the material parameter vector is used to represent the physical and mechanical properties of each layer of material in the geomembrane structure, such as Young's modulus, Poisson's ratio, yield strength, etc., which directly affect the stress response of the structure; the boundary condition vector is used to describe the boundary conditions followed by the geomembrane structure during simulation, such as fixed constraint, free boundary or contact boundary, etc., which determines how the structure is limited by the outside world; the time step parameter vector is used to control the size of the time step in dynamic simulation, which determines the accuracy and stability of the simulation in finite element analysis, especially when simulating non-static load and time-varying process.

[0061] A direct contact stiffness model is constructed by normal contact stiffness, normal projection matrix, tangential contact stiffness and tangential projection matrix. According to the contact condition and material properties between layers, the stiffness of the contact surface in the normal direction is calculated to obtain the normal contact stiffness, which determines the deformation resistance of the contact surface when subjected to vertical pressure. High stiffness means that the contact surface is not easy to deform, and low stiffness means that the contact surface is more likely to deform. According to the friction condition and relative sliding characteristics between contact surfaces, the stiffness of the contact surface in the tangential direction is calculated to obtain the tangential contact stiffness, which affects the friction between contact surfaces and determines whether the membrane layer slips and the degree of slip. The normal projection matrix and the tangential projection matrix are used to decompose the force and deformation on the contact surface into normal and tangential components, which helps to more accurately analyze the interaction between contact surfaces in three-dimensional space. The direct contact stiffness model is a numerical model for describing the mechanical properties of the direct contact interface, which defines the resistance to deformation of the interface in the normal and tangential directions through four core parameters, and can accurately reflect the stiffness characteristics of the contact surface.

[0062] Coulomb friction constraint is a numerical implementation method based on Coulomb friction law, which is used to limit the maximum shear stress between contact surfaces. When the interface shear stress reaches a certain critical value, relative sliding is allowed, which is determined by the interface friction characteristics and the normal pressure. The constructed stress deformation input vector set is integrated with the contact stiffness model established specifically for direct contact characteristics. During the solution process, Coulomb friction constraint is introduced as an important boundary condition, which determines the contact behavior by judging the stress state at each contact point: when the shear stress is less than the maximum static friction, the contact surface is in the sticking state; when the shear stress reaches the critical value, the contact surface enters the sliding state.

[0063] Incremental iterative algorithm is adopted, and multiple iterations are performed within each load increment step until the convergence criterion is met. First, the contact stiffness matrix is calculated according to the current contact state, then the global stiffness matrix is assembled, the displacement increment is solved, and the stress state is updated according to the displacement increment, finally the contact state is checked for changes and convergence is judged. Through the iteration process, three-dimensional stress deformation direct simulation data is obtained, among which the stress displacement increment records the deformation development of each calculation step, the interlayer normal stress reflects the compression degree of the contact surface, the interlayer shear stress shows the distribution of interface friction, and the stress field distribution provides the overall mechanical state of the entire structure. Stress displacement increment represents the increment between stress and deformation caused by external load or other influences; interlayer normal stress is used to describe the stress between layers in the normal direction, which is usually caused by external load or compression between film layers; interlayer shear stress is used to describe the stress between layers in the tangential direction, which is usually related to friction, slip and shear deformation; stress field distribution is used to represent the stress state at different positions in the entire geomembrane structure, which is obtained by finite element solution and used to analyze the stress response of materials. For example, the maximum displacement increment of stress displacement increment occurs at the 6th load step, with a value of 3.2 cm; the interlayer normal stress ranges from 35 to 650 kPa, reaching the maximum value at the bottom of the embankment; the interlayer shear stress in the slope area reaches 12 to 18 kPa, and in the bottom area reaches 5 to 8 kPa; the stress field distribution shows that the maximum tensile stress of the geomembrane is 2.1 MPa, occurring at the junction of the slope and the bank slope; 85% of the contact elements are in the sticking state, and 15% are in the sliding state.

[0064] By constructing the stress deformation input vector set and combining the normal and tangential contact stiffness models, accurate three-dimensional finite element solution is carried out to simulate the stress and deformation behavior of geomembrane structure under complex load. Through strict Coulomb friction constraint and incremental iterative algorithm, the complex mechanical behavior of the contact surface is truly reproduced, including the transition process of sticking-sliding state, which improves the accuracy of the simulation results.

[0065] Further, the application further comprises the steps of: constructing an indirect contact weak coupling model, the indirect contact weak coupling model comprising a stress attenuation coefficient; introducing a material degradation coefficient to identify an equivalent material matrix of the indirect contact surface element; and performing three-dimensional finite element solving on the stress deformation input vector set according to the stress attenuation coefficient of the indirect contact weak coupling model and the equivalent material matrix to obtain three-dimensional stress deformation indirect simulation data.

[0066] Further, the application further comprises the steps of: the indirect contact weak coupling model comprising a stress attenuation coefficient, the stress attenuation coefficient being obtained by fitting calculation of a defect attenuation factor, a quality attenuation factor and a pollution attenuation factor.

[0067] Specifically, the indirect contact surface element refers to those areas that do not directly contact and the transmission of contact force is weakened due to construction defects, environmental factors, etc. The indirect contact weak coupling model is used to simulate the contact behavior of these areas, taking into account that the transmission of force will be affected by weakening factors such as material aging, pollution or construction defects. Based on the specific defect type of the indirect contact element identified in the early stage, the indirect contact weak coupling model is constructed. The indirect contact weak coupling model comprises a stress attenuation coefficient for describing the degree of stress attenuation of the contact surface due to various factors when subjected to force. The stress attenuation coefficient is determined by fitting by comprehensively considering three attenuation factors, i.e. a defect attenuation factor, a quality attenuation factor and a pollution attenuation factor. The defect attenuation factor describes the degree of stress attenuation in the geomembrane structure due to construction defects, and a high defect attenuation factor indicates that the contact surface has a large defect, resulting in a significant weakening of force transmission; the quality attenuation factor represents the degree of stress attenuation in the geomembrane structure due to substandard construction quality, and a high quality attenuation factor indicates that the bearing capacity of the structure is weakened; the pollution attenuation factor describes the influence of environmental pollution on the geomembrane material, resulting in a decrease in its mechanical properties, and a high pollution attenuation factor indicates a decrease in the friction force and compression resistance between the contact surfaces. For example, there are continuous wrinkles in the slope area with an average height of 2.5 cm, and the attenuation factor is determined to be 0.65 through calibration test; the compaction degree detection value of this slope area is 82%, and the corresponding quality attenuation factor is 0.75; chemical analysis shows that the pH value is 3.5 and the heavy metal ion concentration exceeds the standard, and the pollution attenuation factor is determined to be 0.45. The stress attenuation coefficient is calculated to be 0.59 by using weighted geometric mean with weights of 0.4, 0.3 and 0.3, respectively.

[0068] A material degradation coefficient is introduced to identify the equivalent material matrix of the indirect contact surface element, which reflects the degree of material degradation over time or environmental changes, and generally considers the impact of the external environment on the geomembrane material. In finite element analysis, the equivalent material matrix is used to describe the changes in the mechanical properties of the material after being affected by degradation, pollution, defects, and other factors. Through the equivalent material matrix, complex mechanical calculations are simplified, making it better reflect the material behavior under real working conditions. For areas with chemical pollution or aging, the material degradation coefficient will reduce the elastic modulus and strength parameters of the material. According to the material parameters after degradation, the equivalent material matrix is established, simplifying the complex multi-layer system into a single-layer medium with equivalent mechanical properties.

[0069] The stress attenuation coefficient and the equivalent material matrix are integrated into the finite element solver to perform three-dimensional finite element solving on the stress deformation input vector set. During the solving process, the stress attenuation coefficient is used to adjust the stress transfer ability of the contact surface, and the equivalent material matrix is used to calculate the stiffness contribution of the element. Through the joint simulation method of the equivalent medium model + multi-field coupling attenuation model, three-dimensional stress deformation indirect simulation data are obtained, accurately reflecting the special mechanical behavior of the defect interface. Multi-field coupling refers to the mutual influence between different physical fields, such as stress field, temperature field, and humidity field. The interaction of these physical fields leads to the attenuation of the stress of the contact surface. The multi-field coupling attenuation model is a mathematical model that describes the stress change of the geomembrane contact surface under the joint action of multiple factors. Specifically, the equivalent stress of the contact surface is calculated based on the current stress state and the stress attenuation coefficient, the system stiffness matrix is assembled based on the equivalent material matrix, then the displacement increment is obtained by solving the balance equation, and the stress and strain state is updated according to the displacement increment, finally the convergence of the system is checked. Iteration is continued until the preset convergence criterion is met. The three-dimensional stress deformation indirect simulation data are the calculated mechanical response data, including the full-field distribution information of stress, strain, displacement, and other fields considering the influence of the interface defects.

[0070] By constructing an indirect contact weak coupling model, the influence of defects, mass, and pollution on the geomembrane structure is considered, and the mechanical behavior of the membrane layer contact surface is simulated. The introduction of material degradation coefficient and equivalent material matrix enables the simulation to accurately reflect the degradation and weakening phenomena of the material during actual use. Through three-dimensional finite element solving, more accurate stress and deformation data are obtained.

[0071] Fusion of the three-dimensional stress deformation direct simulation data and the three-dimensional stress deformation indirect simulation data obtains the three-dimensional stress deformation fusion simulation data of the geomembrane structure.

[0072] Further, the application also comprises the following steps: defining direct-fusion smoothing weight, indirect-fusion smoothing weight and hybrid-fusion smoothing weight; fusing the direct contact surface unit and indirect contact surface unit as fusion nodes according to the direct-fusion smoothing weight, indirect-fusion smoothing weight and hybrid-fusion smoothing weight to obtain three-dimensional stress deformation fusion simulation data.

[0073] Specifically, based on the cell type and spatial position relationship, three different fusion smoothing weights are defined. The direct-fusion smoothing weight is a weight coefficient used in the data fusion process of direct contact surface units and direct contact surface units. In the simulation process, the fusion of direct contact surface units will be adjusted according to this smoothing weight. A larger smoothing weight means that the fusion between direct contact surface units is stronger, and the deformation and stress transfer will be more uniform. The indirect-fusion smoothing weight is a weight coefficient used in the data fusion process of indirect contact surface units and indirect contact surface units, which is used to maintain a reasonable smooth transition when processing defect area data. Indirect contact surface units usually involve attenuation effects, so the setting of their smoothing weight affects the stress transfer and deformation smoothing degree between contact surfaces. The hybrid-fusion smoothing weight is a weight coefficient used in the data fusion process of direct contact surface units and indirect contact surface units, which is used to handle the data connection problem between different types of units. Hybrid fusion is usually the most complex in practical applications, because the two types of units have significant differences in mechanical behavior, and the hybrid smoothing weight determines the transition mode from direct contact to indirect contact.

[0074] The direct-fusion smoothing weight is high, such as 0.8, which ensures the strict continuity of the data inside the good contact area; the indirect-fusion smoothing weight is moderate, such as 0.5, which considers the influence of defects while maintaining reasonable smoothing inside the defect area; the hybrid-fusion smoothing weight adopts a distance-based gradient function, such as 0.6, which forms a smooth transition zone at the junction of the direct-indirect contact area.

[0075] All the boundary nodes of the direct contact surface units and the indirect contact surface units are taken as fusion nodes to construct a fusion network covering the entire contact surface. At each fusion node, the corresponding weight coefficient is automatically selected according to the node type combination where the node is located. For nodes located inside the direct contact area, the direct-fusion smoothing weight is used; for nodes located inside the indirect contact area, the indirect-fusion smoothing weight is used; for nodes located at the junction of the two types of areas, the hybrid-fusion smoothing weight is used.

[0076] In the specific fusion calculation, according to the contribution degree of the direct contact unit and the indirect contact unit, the corresponding smoothing weight is used for weighted summation, so as to ensure that the stress and deformation results between the contact surface units are smoothly transitioned at the transition area. The fusion process includes smoothing of stress displacement, weighting of normal stress and tangential stress, and continuity adjustment of shear stress and stress field distribution. The three-dimensional stress deformation fusion simulation data is obtained by synthesizing the direct contact and indirect contact simulation data after completing the fusion of each contact surface unit, and is the final stress and deformation distribution data, which accurately reflects the stress distribution, deformation condition and response behavior of the entire geomembrane structure under different conditions. For example, the direct-fusion smoothing weight is set to 0.8, indicating that the stress transmission is relatively uniform; the indirect-fusion smoothing weight is set to 0.5, indicating that the attenuation effect is strong and the contact stress transmission is poor; the mixed-fusion smoothing weight is set to 0.6, which is used to smoothly transition the change of stress and deformation. The final three-dimensional stress deformation fusion simulation data includes: the stress displacement increment in the direct contact area is 1.2mm, the normal stress is 200kPa, and the shear stress is 120kPa; the stress displacement increment in the indirect contact area is 2.5mm, the normal stress is 80kPa, and the shear stress is 50kPa; the stress field distribution diagram after fusion shows that the stress in the contact area is smoothly transitioned in the transition area, avoiding stress mutation.

[0077] By defining the smoothing weight of the direct contact, indirect contact and mixed contact unit and performing fusion, the stress and deformation behavior of the geomembrane structure under different contact states can be effectively simulated. By weighted fusion of stress and deformation data, the stress and deformation mutation in the transition area is avoided, and the accuracy and reliability of the simulation results are improved.

[0078] In summary, the three-dimensional stress deformation simulation analysis method for the geomembrane structure provided in the present application has the following technical effects:

[0079] The geotechnical membrane layered modeling data is constructed by parameterized modeling of the geotechnical membrane structure; the interlayer contact surface is identified according to the geotechnical membrane layered modeling data, and the interlayer contact surface is discretized to obtain a limited number of contact surface units according to a preset grid; the interlayer contact defects, the interlayer construction quality and the environmental pollution state corresponding to the limited number of contact surface units are collected, and the interlayer contact defects, the interlayer construction quality and the environmental pollution state are taken as grading features to divide the limited number of contact surface units into direct contact surface units and indirect contact surface units; stress deformation load data is collected to simulate the direct contact surface units and the indirect contact surface units to obtain three-dimensional stress deformation direct simulation data and three-dimensional stress deformation indirect simulation data; and the three-dimensional stress deformation direct simulation data and the three-dimensional stress deformation indirect simulation data are fused to obtain three-dimensional stress deformation fusion simulation data of the geotechnical membrane structure. That is, the geotechnical membrane layered modeling data is constructed by parameterized modeling, the interlayer contact surface is discretized by a preset grid, the complex interactions such as friction, sliding and voiding between the geotechnical membrane and the overlying protective layer and the underlying supporting layer are simulated by using contact units or interface units, the interlayer contact defects, the construction quality and the environmental pollution state are collected as grading standards, the contact surface units are divided into direct contact surface units and indirect contact surface units, and when the contact stress deformation simulation is performed, a direct and indirect fusion method is used to analyze the deformation simulation state between layers, thereby improving the data accuracy of simulation analysis.

[0080] In the second embodiment, based on the same inventive concept as the three-dimensional stress deformation simulation analysis method for the geotechnical membrane structure in the first embodiment, the present application also provides a three-dimensional stress deformation simulation analysis system for the geotechnical membrane structure. Please refer to the accompanying drawings Figure 2 The three-dimensional stress deformation simulation analysis system for the geotechnical membrane structure comprises:

[0081] The parameterized modeling module 11 is configured to perform parameterized modeling on the geomembrane structure to construct geomembrane layered modeling data; the discretization processing module 12 is configured to identify an interlayer contact surface according to the geomembrane layered modeling data, and perform discretization processing on the interlayer contact surface according to a preset grid to obtain a finite number of contact surface elements; the contact surface element division module 13 is configured to collect interlayer contact defects, interlayer construction quality and environmental pollution states corresponding to the finite number of contact surface elements, divide the finite number of contact surface elements into direct contact surface elements and indirect contact surface elements by taking the interlayer contact defects, the interlayer construction quality and the environmental pollution states as grading features; the simulation analysis module 14 is configured to collect stress deformation load data to simulate the direct contact surface elements and the indirect contact surface elements to obtain three-dimensional stress deformation direct simulation data and three-dimensional stress deformation indirect simulation data; and the data fusion module 15 is configured to fuse the three-dimensional stress deformation direct simulation data and the three-dimensional stress deformation indirect simulation data to obtain three-dimensional stress deformation fusion simulation data of the geomembrane structure.

[0082] Further, the discretization processing module 12 in the three-dimensional stress deformation simulation analysis system for the geomembrane structure is further configured to: identify stress contact areas and non-stress contact areas based on the interlayer contact surface under the stress deformation load data; perform discretization processing on the stress contact areas by using a first preset grid to obtain a first group of contact surface elements, and perform discretization processing on the non-stress contact areas by using a second preset grid to obtain a second group of contact surface elements, wherein the step length of the first preset grid is smaller than the step length of the second preset grid; and obtain a finite number of contact surface elements according to the first group of contact surface elements and the second group of contact surface elements.

[0083] Further, the discretization processing module 12 in the three-dimensional stress deformation simulation analysis system for the geomembrane structure is further configured to: identify a stress transition contact area based on the interlayer contact surface under the stress deformation load data, the stress transition contact area being a transition contact area between the stress contact areas and the non-stress contact areas; perform discretization processing on the stress transition contact area by using a third preset grid to obtain a third group of contact surface elements, wherein the step length of the third preset grid is greater than the step length of the first preset grid and smaller than the step length of the second preset grid; and obtain a finite number of contact surface elements according to the first group of contact surface elements, the second group of contact surface elements and the third group of contact surface elements.

[0084] Further, the contact surface element division module 13 in the geotechnical membrane structure-oriented three-dimensional stress deformation simulation analysis system is further configured to: calculate the contribution degree score of each grading feature of the finite contact surface elements with respect to the interlayer contact defects, the interlayer construction quality, and the environmental pollution state respectively, to obtain a defect contribution degree score, a construction contribution degree score, and a pollution contribution degree score; perform weighted calculation on the defect contribution degree score, the construction contribution degree score, and the pollution contribution degree score by configuring a contribution degree weight, to obtain a contribution degree comprehensive score of the finite contact surface elements; divide the contact surface elements less than a preset score threshold into direct contact surface elements; and divide the contact surface elements greater than or equal to the preset score threshold into indirect contact surface elements.

[0085] Further, the contact surface element division module 13 in the geotechnical membrane structure-oriented three-dimensional stress deformation simulation analysis system is further configured to: construct equivalent thin-layer mechanical parameters of the indirect contact surface elements; simulate a stress change influence degree based on the equivalent thin-layer mechanical parameters according to the stress deformation load data, the stress change influence degree being a data change degree between equivalent simulation obtained stress deformation output data and the stress deformation load data; and re-mark the indirect contact surface elements with a stress change influence degree less than a preset influence degree as direct contact surface elements.

[0086] Further, the simulation analysis module 14 in the geotechnical membrane structure-oriented three-dimensional stress deformation simulation analysis system is further configured to: based on the stress deformation load data and the state input data of the direct contact surface elements, construct a stress deformation input vector set based on stress deformation, including an external load vector, a material parameter vector, a boundary condition vector, and a time step parameter vector; construct a direct contact stiffness model through normal contact stiffness, normal projection matrix, tangential contact stiffness, and tangential projection matrix; perform three-dimensional finite element solving on the stress deformation input vector set according to the direct contact stiffness model by using Coulomb friction constraint, to obtain three-dimensional stress deformation direct simulation data, including stress displacement increment, interlayer normal stress, interlayer shear stress, and stress field distribution.

[0087] Further, the simulation analysis module 14 in the geotechnical membrane structure-oriented three-dimensional stress deformation simulation analysis system is further configured to: construct an indirect contact weak coupling model, the indirect contact weak coupling model including a stress attenuation coefficient; introduce a material degradation coefficient to identify an equivalent material matrix of the indirect contact surface elements; perform three-dimensional finite element solving on the stress deformation input vector set according to the stress attenuation coefficient of the indirect contact weak coupling model and the equivalent material matrix, to obtain three-dimensional stress deformation indirect simulation data.

[0088] Further, the simulation analysis module 14 in the three-dimensional stress deformation simulation analysis system for geomembrane structure is further used for: the indirect contact weak coupling model comprises a stress attenuation coefficient, and the stress attenuation coefficient is obtained by fitting calculation of a defect attenuation factor, a mass attenuation factor and a pollution attenuation factor.

[0089] Further, the data fusion module 15 in the three-dimensional stress deformation simulation analysis system for geomembrane structure is further used for: defining a direct-fusion smoothing weight, an indirect-fusion smoothing weight and a hybrid-fusion smoothing weight; taking the direct contact surface element and the indirect contact surface element as fusion nodes, and fusing according to the direct-fusion smoothing weight, the indirect-fusion smoothing weight and the hybrid-fusion smoothing weight to obtain three-dimensional stress deformation fusion simulation data.

[0090] The embodiments in the specification are described in a progressive manner, and each embodiment focuses on the difference from other embodiments. The three-dimensional stress deformation simulation analysis method for geomembrane structure in the first embodiment and the specific examples are also applicable to the three-dimensional stress deformation simulation analysis system for geomembrane structure in the present embodiment. Through the foregoing detailed description of the three-dimensional stress deformation simulation analysis method for geomembrane structure, those skilled in the art can clearly know the three-dimensional stress deformation simulation analysis system for geomembrane structure in the present embodiment. Therefore, for the sake of brevity of the specification, the three-dimensional stress deformation simulation analysis system for geomembrane structure in the present embodiment will not be described in detail.

[0091] The above description of disclosed embodiments enables a person skilled in the art to implement or use the present application. Various modifications to these embodiments will be apparent to those skilled in the art, and the general principles defined herein can be implemented in other embodiments without departing from the spirit or scope of the present application. Therefore, the present application will not be limited to the embodiments shown herein, but will conform to the widest scope consistent with the principles and novel features disclosed herein.

[0092] Obviously, those skilled in the art can make various modifications and variations to the present application without departing from the spirit and scope of the present application. Thus, if these modifications and variations of the present application fall within the scope of the present application and its equivalent technologies, the present application also intends to include these modifications and variations.

Claims

1. A three-dimensional stress-deformation simulation and analysis method for geomembrane structures, characterized in that, include: Parametric modeling of geomembrane structures is performed to construct layered modeling data for geomembrane structures; Based on the geomembrane layer modeling data, the interlayer contact surfaces are identified, and the interlayer contact surfaces are discretized according to a preset grid to obtain a finite number of contact surface units; The interlayer contact defects, interlayer construction quality, and environmental pollution status corresponding to the finite number of contact surface units are collected. The interlayer contact defects, interlayer construction quality, and environmental pollution status are used as classification features to divide the finite number of contact surface units into direct contact surface units and indirect contact surface units. Stress-deformation load data are collected to simulate the direct contact surface unit and the indirect contact surface unit, thereby obtaining three-dimensional stress-deformation direct simulation data and three-dimensional stress-deformation indirect simulation data. By integrating the direct simulation data and indirect simulation data of three-dimensional stress and deformation, the three-dimensional stress and deformation fusion simulation data of the geomembrane structure is obtained.

2. The three-dimensional stress-deformation simulation and analysis method for geomembrane structures as described in claim 1, characterized in that, The method further includes discretizing the interlayer contact surface to obtain a finite number of contact surface elements: Identify the stress contact area and non-stress contact area based on the interlayer contact surface under the stress deformation load data; The stress contact region is discretized using a first preset grid to obtain a first set of contact surface elements, and the non-stress contact region is discretized using a second preset grid to obtain a second set of contact surface elements, wherein the step size of the first preset grid is smaller than the step size of the second preset grid. Based on the first group of contact surface units and the second group of contact surface units, a finite number of contact surface units are obtained.

3. The three-dimensional stress-deformation simulation and analysis method for geomembrane structures as described in claim 2, characterized in that, The method further includes discretizing the interlayer contact surface to obtain a finite number of contact surface elements: Identify the stress transition contact region based on the interlayer contact surface under the stress deformation load data, wherein the stress transition contact region is the transition contact region between the stress contact region and the non-stress contact region; The stress transition contact region is discretized using a third preset grid to obtain a third set of contact surface elements, wherein the step size of the third preset grid is greater than the step size of the first preset grid and less than the step size of the second preset grid. A finite number of contact surface units are obtained based on the first group of contact surface units, the second group of contact surface units, and the third group of contact surface units.

4. The three-dimensional stress-deformation simulation and analysis method for geomembrane structures as described in claim 1, characterized in that, The method involves using interlayer contact defects, interlayer construction quality, and environmental pollution status as classification characteristics to divide the finite number of contact surface units into direct contact surface units and indirect contact surface units. Calculate the contribution score of each of the finite number of contact surface units to each graded feature in terms of interlayer contact defects, interlayer construction quality, and environmental pollution status, and obtain the defect contribution score, construction contribution score, and pollution contribution score. By configuring contribution weights, the defect contribution score, construction contribution score, and pollution contribution score are weighted and calculated to obtain the comprehensive contribution score of the finite number of contact surface units. Contact surface units that are below a preset scoring threshold are classified as direct contact surface units; Contact surface units that are greater than or equal to a preset scoring threshold are classified as indirect contact surface units.

5. The three-dimensional stress-deformation simulation and analysis method for geomembrane structures as described in claim 1, characterized in that, The method for acquiring stress-deformation load data and simulating the direct contact surface element to obtain three-dimensional stress-deformation direct simulation data includes: Based on the stress-deformation load data and the state input data of the direct contact surface unit, a stress-deformation input vector set is constructed, including external load vector, material parameter vector, boundary condition vector, and time step parameter vector; A direct contact stiffness model is constructed using normal contact stiffness, normal projection matrix, tangential contact stiffness, and tangential projection matrix. Using Coulomb friction constraints, the stress-deformation input vector set is solved in three dimensions using the direct contact stiffness model to obtain three-dimensional stress-deformation direct simulation data, including stress-displacement increment, interlayer normal stress, interlayer shear stress, and stress field distribution.

6. The three-dimensional stress-deformation simulation and analysis method for geomembrane structures as described in claim 5, characterized in that, The method for acquiring stress-deformation load data and simulating the indirect contact surface element to obtain three-dimensional stress-deformation indirect simulation data includes: Construct an indirect contact weak coupling model, which includes a stress attenuation coefficient; A material degradation coefficient is introduced to identify the equivalent material matrix of the indirect contact surface unit; Based on the stress attenuation coefficient of the indirect contact weak coupling model and the equivalent material matrix, the stress deformation input vector set is solved by three-dimensional finite element method to obtain three-dimensional stress deformation indirect simulation data.

7. The three-dimensional stress-deformation simulation and analysis method for geomembrane structures as described in claim 6, characterized in that, The indirect contact weak coupling model includes a stress attenuation coefficient, which is obtained by fitting and calculating the defect attenuation factor, the mass attenuation factor, and the contamination attenuation factor.

8. The three-dimensional stress-deformation simulation and analysis method for geomembrane structures as described in claim 1, characterized in that, The method for obtaining fused three-dimensional stress-deformation simulation data of the geomembrane structure by integrating the direct and indirect three-dimensional stress-deformation simulation data includes: Define direct-fusion smoothing weights, indirect-fusion smoothing weights, and hybrid-fusion smoothing weights; Using the direct contact surface unit and the indirect contact surface unit as fusion nodes, fusion is performed according to the direct-fusion smoothing weight, the indirect-fusion smoothing weight and the hybrid-fusion smoothing weight to obtain three-dimensional stress-deformation fusion simulation data.

9. The three-dimensional stress-deformation simulation and analysis method for geomembrane structures as described in claim 4, characterized in that, After dividing the elements into direct contact surface elements and indirect contact surface elements, the method also includes: Construct the equivalent thin-layer mechanical parameters of the indirect contact surface unit; The stress variation influence degree is simulated based on the stress deformation load data and the equivalent thin-layer mechanical parameters. The stress variation influence degree is the data variation degree between the stress deformation output data obtained from the equivalent simulation and the stress deformation load data. Indirect contact surface elements whose stress change influence is less than the preset influence value are re-marked as direct contact surface elements.

10. A three-dimensional stress-deformation simulation and analysis system for geomembrane structures, characterized in that, The steps for implementing the three-dimensional stress-deformation simulation analysis method for geomembrane structures according to any one of claims 1 to 9, wherein the three-dimensional stress-deformation simulation analysis system for geomembrane structures comprises: The parametric modeling module is used to perform parametric modeling of geomembrane structures and construct geomembrane layer modeling data. The discretization module is used to identify the interlayer contact surface based on the geomembrane layer modeling data, and to discretize the interlayer contact surface according to a preset grid to obtain a finite number of contact surface units. The contact surface unit division module is used to collect the interlayer contact defects, interlayer construction quality and environmental pollution status corresponding to the finite number of contact surface units, and to divide the finite number of contact surface units into direct contact surface units and indirect contact surface units using the interlayer contact defects, interlayer construction quality and environmental pollution status as classification features. The simulation analysis module is used to collect stress and deformation load data to simulate the direct contact surface unit and the indirect contact surface unit, and obtain three-dimensional stress and deformation direct simulation data and three-dimensional stress and deformation indirect simulation data. The data fusion module is used to fuse the direct simulation data of three-dimensional stress and deformation and the indirect simulation data of three-dimensional stress and deformation to obtain the three-dimensional stress and deformation fused simulation data of the geomembrane structure.

Images