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

Abstract

This application provides a method and system for three-dimensional stress-deformation simulation analysis of geomembrane structures, relating to the field of geomembrane structure technology. The method includes: parametrically modeling the geomembrane structure; identifying interlayer contact surfaces based on the geomembrane layer modeling data, and discretizing them according to a preset mesh to obtain a finite number of contact surface elements; dividing the finite number of contact surface elements into direct contact surface elements and indirect contact surface elements; collecting stress-deformation load data 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 fused simulation data of the geomembrane structure through data fusion. This application solves the technical problem in the prior art where the entire interlayer contact surface is treated as a homogeneous body, ignoring interlayer contact defects, resulting in insufficient accuracy of three-dimensional stress-deformation simulation results, thus improving the accuracy of three-dimensional stress-deformation simulation of geomembrane structures.

Description

Technical Field

[0001] This application relates to the field of geomembrane structure technology, specifically to a three-dimensional stress-deformation simulation analysis method and system for geomembrane structures. Background Technology

[0002] Existing simulation methods typically simplify the entire interlayer contact surface as a homogeneous body with uniform mechanical properties, ignoring the spatial non-uniform distribution of interface properties caused by various factors in actual engineering. Under this simplification, interlayer interactions are often described as a uniform contact model using uniform contact parameters. However, in practical engineering applications, due to the influence of various factors, the contact state between the geomembrane and its overlying protective layer and underlying support layer will change significantly, leading to non-uniformity of the mechanical properties between the membrane layers. This affects the stress distribution and deformation behavior of the membrane layers, making it impossible to accurately simulate the three-dimensional stress-deformation of geomembrane structures.

[0003] In summary, the existing technology suffers from a technical problem: by treating the entire interlayer contact surface as a homogeneous body and ignoring interlayer contact defects, the accuracy of the three-dimensional stress-deformation simulation results is insufficient. Summary of the Invention

[0004] The purpose of this application is to provide a three-dimensional stress-deformation simulation analysis method and system for geomembrane structures, in order to solve the technical problem in the prior art that the accuracy of the three-dimensional stress-deformation simulation results is insufficient because the entire interlayer contact surface is regarded as a homogeneous body and the interlayer contact defects are ignored.

[0005] To achieve the above objectives, this application provides a three-dimensional stress-deformation simulation analysis method and system for geomembrane structures.

[0006] Firstly, this application provides a three-dimensional stress-deformation simulation analysis method for geomembrane structures. This method is implemented using a three-dimensional stress-deformation simulation analysis system for geomembrane structures. The method includes: parametrically modeling the geomembrane structure to construct layered geomembrane modeling data; identifying interlayer contact surfaces based on the layered geomembrane modeling data, and discretizing the interlayer contact surfaces according to a preset grid to obtain a finite number of contact surface units; collecting interlayer contact defects, interlayer construction quality, and environmental pollution status corresponding to the finite number of contact surface units, and using these as hierarchical features to divide the finite number of contact surface units into direct contact surface units and indirect contact surface units; collecting stress-deformation load data to simulate the direct and indirect contact surface units 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 and indirect simulation data to obtain three-dimensional stress-deformation fused simulation data of the geomembrane structure.

[0007] Optionally, the stress contact region and non-stress contact region based on the interlayer contact surface are identified 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 units; the non-stress contact region is discretized using a second preset grid to obtain a second set of contact surface units, wherein the step size of the first preset grid is smaller than the step size of the second preset grid; a finite number of contact surface units are obtained based on the first set of contact surface units and the second set of contact surface units.

[0008] Optionally, the stress transition contact region based on the interlayer contact surface is identified under the stress deformation load data. The stress transition contact region is a 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. 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 elements are obtained based on the first set of contact surface elements, the second set of contact surface elements, and the third set of contact surface elements.

[0009] Optionally, the contribution scores of the finite number of contact surface units with respect to each graded feature of interlayer contact defects, interlayer construction quality, and environmental pollution status are calculated respectively to obtain defect contribution scores, construction contribution scores, and pollution contribution scores; the defect contribution scores, construction contribution scores, and pollution contribution scores are weighted by configuring contribution weights to obtain a comprehensive contribution score of the finite number of contact surface units; contact surface units with scores less than a preset score threshold are classified as direct contact surface units; contact surface units with scores greater than or equal to the preset score threshold are classified as indirect contact surface units.

[0010] Optionally, the equivalent thin-layer mechanical parameters of the indirect contact surface unit are constructed; the stress change influence degree based on the equivalent thin-layer mechanical parameters is simulated according to the stress deformation load data, wherein the stress change influence degree is the data change degree between the stress deformation output data obtained by the equivalent simulation and the stress deformation load data; and the indirect contact surface units with stress change influence degree less than the preset influence degree are re-marked as direct contact surface units.

[0011] Optionally, 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 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, a three-dimensional finite element method is performed on the stress-deformation input vector set according to 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.

[0012] Optionally, an indirect contact weak coupling model is constructed, which includes a stress attenuation coefficient; a material degradation coefficient is introduced to identify the equivalent material matrix of the indirect contact surface element; and a three-dimensional finite element method is used to solve the stress deformation input vector set based on 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, which is obtained by fitting and calculating the defect attenuation factor, the mass attenuation factor, and the contamination attenuation factor.

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

[0015] Secondly, this application also provides a three-dimensional stress-deformation simulation and analysis system for geomembrane structures, used to execute the three-dimensional stress-deformation simulation and analysis method for geomembrane structures as described in the first aspect. The three-dimensional stress-deformation simulation and analysis system for geomembrane structures includes: a parametric modeling module for parametrically modeling the geomembrane structure to construct geomembrane layered modeling data; a discretization module for identifying interlayer contact surfaces based on the geomembrane layered modeling data and discretizing the interlayer contact surfaces according to a preset grid to obtain a finite number of contact surface units; and a contact surface unit partitioning module for collecting data from the finite number of contact surface units. The system uses 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. A simulation analysis module is used to collect stress-deformation load data to simulate the direct and indirect contact surface units, obtaining three-dimensional stress-deformation direct simulation data and three-dimensional stress-deformation indirect simulation data. A data fusion module is used to fuse the three-dimensional stress-deformation direct and indirect simulation data to obtain three-dimensional stress-deformation fused simulation data of the geomembrane structure.

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

[0017] The geomembrane structure is parametrically modeled to construct layered geomembrane modeling data. Based on this data, interlayer contact surfaces are identified, and these surfaces are discretized according to a preset grid to obtain a finite number of contact surface units. Interlayer contact defects, construction quality, and environmental pollution status corresponding to these contact surface units are collected. These factors are used as hierarchical features to classify 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 and indirect contact surface units, obtaining three-dimensional stress-deformation direct simulation data and three-dimensional stress-deformation indirect simulation data. The three-dimensional stress-deformation direct and indirect simulation data are then fused to obtain the three-dimensional stress-deformation fused simulation data of the geomembrane structure. In other words, by constructing geomembrane layer modeling data through parametric modeling, the interlayer contact surfaces are discretized through a preset mesh. Contact units or interface units are used to simulate complex interactions such as friction, slippage, and voids between the geomembrane and the overlying protective layer and the underlying support layer. Interlayer contact defects, construction quality, and environmental pollution status are collected as grading standards. The contact surface units are divided into direct contact surface units and indirect contact surface units. When simulating contact stress deformation, a direct plus indirect fusion method is used to analyze the deformation simulation state between layers, thereby improving the accuracy of the simulation analysis data.

[0018] The above description is merely an overview of the technical solution of this application. To better understand the technical means of this application and to facilitate its implementation according to the description, and to make the above and other objects, features, and advantages of this application more apparent, specific embodiments of this application are described below. It should be understood that the content described in this section is not intended to identify key or important features of the embodiments of this application, nor is it intended to limit the scope of this application. Other features of this application will become readily apparent through the following description. Attached Figure Description

[0019] To more clearly illustrate the technical solutions in this application or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are merely exemplary. For those skilled in the art, other drawings can be obtained based on the provided drawings without creative effort.

[0020] Figure 1 This is a flowchart illustrating the three-dimensional stress-deformation simulation analysis method for geomembrane structures proposed in this application.

[0021] Figure 2 This is a schematic diagram of the three-dimensional stress-deformation simulation and analysis system for geomembrane structures according to this application.

[0022] Figure labeling: Parametric modeling module 11, Discretization module 12, Contact surface element division module 13, Simulation analysis module 14, Data fusion module 15. Detailed Implementation

[0023] This application provides a three-dimensional stress-deformation simulation analysis method and system for geomembrane structures, solving the technical problem of insufficient accuracy in three-dimensional stress-deformation simulation results caused by treating the entire interlayer contact surface as a homogeneous body and ignoring interlayer contact defects in existing technologies. The method constructs layered geomembrane modeling data through parametric modeling, discretizes the interlayer contact surface using a preset mesh, and uses contact elements or interface elements to simulate complex interactions such as friction, slippage, and voiding between the geomembrane and the overlying protective layer and underlying support layer. Interlayer contact defects, construction quality, and environmental pollution status are collected as classification standards, dividing the contact surface elements into direct contact surface elements and indirect contact surface elements. When simulating contact stress-deformation, a direct and indirect fusion approach is used to analyze the deformation simulation state between layers, improving the accuracy of the simulation analysis data.

[0024] The technical solutions of this application will now be clearly and completely described with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of this application, and not all of them. It should be understood that this application is not limited to the exemplary embodiments described herein. All other embodiments obtained by those skilled in the art based on the embodiments of this application without creative effort are within the scope of protection of this application. It should also be noted that, for ease of description, only the parts related to this application are shown in the accompanying drawings, not all of them.

[0025] Example 1, please refer to the appendix. Figure 1 This application provides a three-dimensional stress-deformation simulation and analysis method for geomembrane structures. The method is applied to a three-dimensional stress-deformation simulation and analysis system for geomembrane structures, and specifically includes the following steps:

[0026] Parametric modeling of geomembrane structures is performed to construct geomembrane layer modeling data.

[0027] Specifically, geomembrane structures are multi-layered composite systems composed of a high-performance polymer film as the core impermeable layer, along with upper and lower protective layers, drainage layers, and support layers. They are commonly found in projects such as landfills. Geomembrane structures typically refer to structural systems composed of multiple geomembrane layers, each with specific thickness, mechanical properties, and impermeability. Based on engineering design drawings and geological survey reports, the entire geomembrane structure is abstracted into several continuous layers, such as, from bottom to top, the foundation, lower protective layer, geomembrane impermeable layer, upper protective layer, drainage layer, and landfill. By analyzing the physical properties of the geomembrane, the geometric dimensions and mechanical properties of each layer are extracted. Based on the actual needs of the project, parameters are set for each layer, including the membrane thickness and the physical properties of the materials used. The parameters of each layer are organized into a structured dataset, and the attributes of each layer are stored in the form of a data table.

[0028] Parametric modeling software was used to create a 3D geometric model for each layer, driven by key parameters. During the modeling process of the geomembrane structure, the data of each membrane layer was superimposed to construct a complete layered modeling data system. The data for each membrane layer includes not only its geometric parameters but also information about the contact surfaces between that membrane and other layers, such as the friction coefficient and contact area. Finally, a 3D digital model containing all the layered structures was automatically generated, and the geomembrane layered modeling data, including the geometric information, topological relationships, and initial material properties of each layer, was output. The geomembrane layered modeling data is a collection of digital information that can completely describe the spatial distribution, geometric dimensions, material properties, and relative positions of each layer of the geomembrane structure.

[0029] By constructing parametric and layered modeling data for the geomembrane structure, the geometric characteristics and physical properties of each membrane layer are clearly defined, ensuring that the simulation accurately reflects the geomembrane structure in actual engineering projects. 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 from scratch in traditional methods.

[0030] 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.

[0031] Furthermore, this application also includes the following steps: identifying the stress contact region and non-stress contact region based on the interlayer contact surface under the stress deformation load data; discretizing the stress contact region using a first preset grid to obtain a first set of contact surface elements; discretizing the non-stress contact region 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; and obtaining a finite number of contact surface elements based on the first set of contact surface elements and the second set of contact surface elements.

[0032] Specifically, based on the geomembrane layer modeling data, all interlayer contact surfaces requiring contact analysis are automatically identified. This involves analyzing the geometric positions and physical contact states between the membrane layers to determine which areas constitute the actual contact surfaces. Interlayer contact surfaces are the contact areas between different membrane layers in a geomembrane structure, determining the mechanical interactions between layers, such as friction and shear forces.

[0033] Import stress-deformation load data, which represents the stress and deformation data of the geomembrane structure under external loads. Based on this data, the identified interlayer contact surfaces are divided into stress contact areas and non-stress contact areas. Stress contact areas are regions where there is significant contact stress or deformation between layers. These areas typically have higher stress levels and require higher calculation accuracy, such as slopes, anchorage trenches, and areas with concentrated loads. Non-stress contact areas, on the other hand, refer to regions where there is no significant contact stress or deformation between layers. These areas require relatively lower calculation accuracy, such as large flat areas at the site bottom. In other words, under stress-deformation load data, it is necessary to identify which contact surface areas are mechanically critical (stress contact areas) and which are secondary (non-stress contact areas).

[0034] In finite element analysis, mesh generation is the process of dividing the entire model space into smaller elements for numerical solution. A preset mesh refers to the mesh generation criteria set before the analysis begins, including mesh size and step size. It is typically optimized based on the importance of the region and stress distribution. In this application, a first preset mesh, a second preset mesh, and a third preset mesh are included to discretize the stress contact region, the non-stress contact region, and the stress transition contact region, respectively.

[0035] A first preset mesh, using a smaller mesh step size, is applied to partition the stress contact region. This generates a large number of small, densely packed elements within the stress contact region, forming the first set of contact surface elements. This set can depict the complex stress gradient changes and potential strain localization phenomena within this region. For non-stress contact regions identified as less important, a second preset mesh is applied, using a larger mesh step size. This generates a second set of contact surface elements. This second set of elements is smaller in number and larger in size, sufficient to describe the overall mechanical behavior of the non-stress contact region without imposing excessive computational burden. The step size of the first preset mesh is smaller than that of the second preset mesh.

[0036] For example, in a simulation of a tailings dam 100m long with a slope of 1:1.5, the stress state of the geomembrane seepage barrier layer under water pressure during the impoundment period needs to be analyzed. The stress contact area is identified as the entire dam slope, with an area of ​​approximately 850m². 2The first preset mesh was used, with a step size set to 0.15m. After discretization, approximately 37,800 fine quadrilateral contact surface elements were generated, forming the first set of contact surface elements. The non-stress contact region was identified as the flat area at the bottom of the reservoir, with an area of ​​approximately 2200m². 2 The second preset mesh was used, with a step size of 0.8m. After discretization, approximately 3400 rough quadrilateral contact surface elements were generated, which constitute the second set of contact surface elements.

[0037] Discretization yields a first set and a second set of contact surface elements, the quantity and shape of which are directly related to the stress distribution and mesh step size in the contact surface region. Based on the meshing of stress-contact and non-stress-contact regions, a finite set of contact surface elements is ultimately generated. Refined meshing for different stress regions and contact characteristics improves computational accuracy in critical areas while reducing computational complexity in non-stress-contact regions, thus minimizing computational resource waste while ensuring simulation accuracy.

[0038] Furthermore, this application also includes the following steps: identifying the stress transition contact region based on the interlayer contact surface under the stress deformation load data, wherein the stress transition contact region is a transition contact region between the stress contact region and the non-stress contact region; discretizing the stress transition contact region 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; and obtaining a finite number of contact surface elements based on the first set of contact surface elements, the second set of contact surface elements, and the third set of contact surface elements.

[0039] Specifically, after identifying the stress contact region and the non-stress contact region, the stress transition region in the interlayer contact surface is identified through analysis of stress-deformation load data. The stress transition contact region refers to the area located between the stress contact region and the non-stress contact region, in which the stress gradually decreases from a larger value. The characteristic of the stress transition contact region is that the stress state changes significantly. It does not require extremely high precision like the core stress region, nor can it be highly simplified like the non-stress region. It is a transition zone from complex to simple mechanical behavior.

[0040] The third preset mesh scheme is invoked, and an intermediate step size larger than the first preset mesh step size but smaller than the second preset mesh step size is used to discretize the stress transition contact region, thereby generating a third set of contact surface elements. The third set of contact surface elements plays a bridging role in terms of size. The fine mesh element sizes in the stress contact region are not significantly different, ensuring the smooth transmission of stress; at the same time, they are relatively close to the coarse mesh element sizes in the non-stress contact region, achieving a gradual decrease in computational scale.

[0041] The first preset mesh is used for the stress contact region with a smaller step size to calculate stress and deformation more accurately; the second preset mesh is used for the non-stress contact region with a larger step size to optimize computational efficiency; the third preset mesh is used for the stress transition contact region with a step size between the first and second preset meshes to ensure the accuracy of the transition region while controlling computational complexity.

[0042] By discretizing the stress contact region, non-stress contact region, and stress transition contact region, three sets of contact surface elements are obtained. These three sets are then seamlessly spliced ​​together in space to form a finite set of contact surface elements with continuously varying mesh density. For example, the stress transition contact region, located between the stress contact region and the non-stress contact region, is a ring approximately 2 meters wide with an area of ​​about 1500 m². 2 Using the third preset grid with a step size of 0.4m, approximately 9400 cells were generated.

[0043] By discretizing the stress transition contact region with a moderate mesh step size, computational efficiency is improved while maintaining computational accuracy. Refined modeling of the stress transition region allows the simulation results to better reflect the mechanical state in actual engineering, especially the stress distribution in the transition region. Accurately dividing the contact surface elements not only improves the accuracy of the simulation but also reduces unnecessary computational burden, thereby increasing the efficiency of the simulation process.

[0044] 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.

[0045] Furthermore, this application also includes the following steps: calculating the contribution score of each graded feature of the finite number of contact surface units with respect to interlayer contact defects, interlayer construction quality, and environmental pollution status, respectively, to obtain defect contribution score, construction contribution score, and pollution contribution score; weighting the defect contribution score, construction contribution score, and pollution contribution score by configuring contribution weights to obtain a comprehensive contribution score of the finite number of contact surface units; classifying contact surface units with scores less than a preset score threshold as direct contact surface units; and classifying contact surface units with scores greater than or equal to the preset score threshold as indirect contact surface units.

[0046] Specifically, after obtaining a limited number of contact surface units, a multi-source data acquisition and fusion process is initiated to collect data on interlayer contact defects, interlayer construction quality, and environmental pollution status. Interlayer contact defects refer to potential flaws at the contact surfaces between layers, such as uneven contact pressure, air bubbles, impurity layers, or membrane unevenness. These defects affect the friction and bearing capacity between membrane layers, thus impacting the mechanical and impermeability properties of the geomembrane structure. Interlayer construction quality refers to the degree of adherence to design requirements and standards during the laying and installation of the geomembrane and its upper and lower protective layers, such as interface smoothness, material compaction, joint quality, and interlayer adhesion. Environmental pollution status refers to the intrusion of external substances into the interlayer interface and their impact on material properties, including physical, chemical, and biological pollution. Factors include voids between membrane layers, membrane surface unevenness, air bubbles, and impurity layers.

[0047] Detailed inspection and testing of the contact surfaces between geomembrane layers are conducted to collect data on contact defects, including gaps between layers, unevenness of the membrane surface, air bubbles, and impurity layers. These defects are then located and associated with their specific contact surface units. Through on-site construction monitoring, construction reports, and quality inspections, quality data regarding the membrane construction process is collected, including construction temperature, construction speed, joint quality between membrane layers, and membrane integrity. Environmental monitoring equipment, such as gas sensors, temperature and humidity sensors, and chemical sensors, is used to collect data on pollutant concentrations and temperature and humidity changes in the environment in which the geomembrane is located, reflecting whether the geomembrane is affected by environmental factors such as chemical corrosion and microbial contamination.

[0048] The contribution scores of a finite number of contact surface units to each graded characteristic in terms of interlayer contact defects, interlayer construction quality, and environmental pollution status are calculated 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; for example, no defects are scored as 0 points, minor wrinkles as 8 points, and severe scratches or holes as 20 points, thus obtaining the defect contribution score. For interlayer construction quality, a reverse scoring is performed based on indicators such as compaction and smoothness; for example, an excellent area with a compaction degree ≥95% is scored as 0 points, a qualified area with a compaction degree of approximately 90% is scored as 6 points, and an unqualified area with a compaction degree <85% is scored as 15 points, thus obtaining the construction contribution score. For environmental pollution status, a score is given based on the nature and degree of pollutants; for example, no pollution is scored as 0 points, minor physical sedimentation as 10 points, and the presence of highly corrosive chemical pollutants as 18 points, thus obtaining the pollution contribution score.

[0049] The contribution weights are configured based on the relative importance of three tiered characteristics—interlayer contact defects, interlayer construction quality, and environmental pollution status—in determining the final contact type. The weighting relies on engineering experience and statistical data; for example, construction quality may have a higher weight in the initial stages of operation, while the weight of environmental pollution will gradually increase over long-term operation.

[0050] The configured contribution weights are used to perform a weighted calculation on the defect contribution score, construction contribution score, and pollution contribution score, ultimately obtaining a comprehensive contribution score for a finite number of contact surface units. Contact surface units with a comprehensive contribution score less than a preset threshold are classified as direct contact surface units, while those with a comprehensive contribution score greater than or equal to the preset threshold are classified as indirect contact surface units. The preset threshold is a dividing value used to classify continuous comprehensive contribution scores into direct and indirect contact categories. For example, with weights set as defect weight 0.4, construction quality weight 0.3, and pollution weight 0.3, and a preset threshold of 10 points, unit A, located in a flat area at the bottom of the reservoir, has no defects, so its defect score is 0; its compaction degree is 96%, so its construction quality score is 0; and it has no pollution, so its pollution score is 0. Therefore, the comprehensive score = (0 × 0.4) + (0 × 0.3) + (0 × 0.3) = 0 points, classifying it as a direct contact surface unit. Unit B is located on a slope with slight folds, including folds 2 cm high. Its defect score is 8, compaction is 92%, and construction quality score is 4. It is unpolluted, with a pollution score of 0. Therefore, the overall score is (8 × 0.4) + (4 × 0.3) + (0 × 0.3) = 3.2 + 1.2 + 0 = 4.4 points, classifying it as a direct contact unit. Unit C is located at the bottom of the slope in a harsh environment. It has a trace of hot welding repair, a defect score of 12, insufficient regional compaction (less than 88%), a construction quality score of 8, and leachate precipitation leading to chemical crystallization, resulting in a pollution score of 15. Therefore, the overall score is (12 × 0.4) + (8 × 0.3) + (15 × 0.3) = 4.8 + 2.4 + 4.5 = 11.7 points, classifying it as an indirect contact unit. By repeating the above operation on the remaining contact units, the final quantitative result is that 85% of the units are classified as direct contact and 15% as indirect contact.

[0051] A direct contact unit is a contact area unit in which the geomembrane and the adjacent material layer are tightly bonded without any intermediate medium. An indirect contact unit is a contact area unit in which there is a third phase medium or physical separation between the geomembrane and the adjacent material layer. The intermediate medium or separation changes the original interface contact state, causing a fundamental change in its mechanical behavior, such as impurity layers, loose layers, construction quality defects, etc.

[0052] By calculating the contribution score of each contact surface unit and combining it with weighted calculations, the impact of interlayer contact defects, construction quality, and environmental pollution on the contact surface units is comprehensively considered. This allows for the classification of contact surface units, accurately identifying which contact surface units require high-precision simulation and which can be simplified.

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

[0054] Specifically, the equivalent thin-layer mechanical parameters of the indirect contact surface unit are constructed, transforming the actual contact surface into a thin layer of a certain thickness and assigning equivalent mechanical parameters to this layer. For units with impurity layers, since the impurity layers are actual physical materials, their mechanical parameters are determined using laboratory testing methods. Through meticulous on-site sampling, samples of the impurity layers in their original state are sent to the laboratory for direct shear tests according to relevant specifications. The direct shear test measures the shear strength of the specimens under different vertical pressures, and the strength envelope is plotted to determine cohesion and internal friction angle. For units with loose zones, since the loose zones are essentially areas where the density of the original material decreases, the equivalent stiffness parameters are determined using indirect on-site testing and back analysis methods. Through ground-penetrating radar data inversion, utilizing the difference in electromagnetic wave propagation speed in materials of different densities, combined with a small amount of borehole calibration data, the equivalent elastic modulus distribution of the region is calculated. For construction quality defects, the strength reduction factor is determined through non-destructive testing or sampling tests. For welding defects, the size and type of defects are determined by ultrasonic testing, and the reduction factor is determined based on the defect rate; for material damage, tensile tests are conducted by sampling, and the strength ratio of the damaged area to the intact material is compared to determine the reduction factor; for laying defects, the geometric parameters of the folds are determined by on-site testing, and the reduction factor is determined based on the ratio of the fold height to the wavelength.

[0055] The constructed equivalent thin-layer mechanical parameters are substituted into the numerical model, and a specific simulation analysis is performed based on stress-deformation load data. During the simulation, the stress variation influence degree of each indirect contact surface element is calculated. The stress variation influence degree is defined by comparing the difference between the mechanical response calculated using the equivalent thin-layer model and the original load data, reflecting the actual importance of the defect interface in the structural system. A smaller stress variation influence degree means a smaller difference between the simulation results and the actual data, indicating higher simulation accuracy.

[0056] Stress-deformation load data refers to the stress and deformation data generated by a geomembrane structure under external loads. The accuracy of the simulation results is measured by calculating the degree of variation between the simulated stress-deformation output data and the actual stress-deformation load data. For a certain indirect contact surface element, the simulated stress-deformation data is 10 kPa and displacement is 5 mm, while the actual load data is 9.5 kPa and displacement is 4.8 mm. The influence of stress variation can be calculated by comparing the difference between these two sets of data. If the influence of stress variation is 0.05, that is, the deviation between the simulation result and the actual data is less than 5%, the simulation result is considered relatively accurate.

[0057] If the stress variation influence of an indirect contact surface element is less than a preset influence threshold, the stress transmission of that element is considered relatively stable, meeting the conditions for direct contact, and it is re-marked as a direct contact element. The preset influence threshold is a pre-defined threshold used to determine whether indirect contact surface elements need to be reclassified, typically determined based on engineering experience, safety standards, and numerical experiments. By constructing equivalent thin-layer mechanical parameters for indirect contact surface elements and calculating the stress variation influence, the simulation accuracy is effectively evaluated, and the reliability of the simulation results is ensured. If the stress variation influence of some indirect contact elements is small, it indicates that their mechanical behavior is close to that of direct contact surface elements, accurately reflecting actual working conditions. By re-marking these elements as direct contact surface elements, the accuracy and stability of the simulation model are further improved.

[0058] 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.

[0059] Furthermore, this application also includes the following steps: based on the stress-deformation load data and the state input data of the direct contact surface element, constructing a stress-deformation input vector set, including external load vector, material parameter vector, boundary condition vector, and time step parameter vector; constructing a direct contact stiffness model through normal contact stiffness, normal projection matrix, tangential contact stiffness, and tangential projection matrix; using Coulomb friction constraints, performing a three-dimensional finite element solution on the stress-deformation input vector set according to 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.

[0060] Specifically, a stress-deformation input vector set is constructed based on stress-deformation load data and the state input data of selected direct contact surface elements. For direct contact surface elements, the state input data includes external loads, material parameters, boundary conditions, and time step parameters. The stress-deformation input vector set is a collection of vectors representing all input parameters affecting stress-deformation calculations, including external loads, material properties, boundary conditions, and time step. The external load vector describes the external loads applied to the geomembrane structure, such as pressure, tension, or bending moment. External loads affect the stress state of the structure and cause deformation through these loads. The material parameter vector represents the physical and mechanical properties of each layer of material in the geomembrane structure, such as Young's modulus, Poisson's ratio, and yield strength. Material parameters directly affect the stress response of the structure. The boundary condition vector describes the boundary conditions followed by the geomembrane structure during the simulation, such as fixed constraints, free boundaries, or contact boundaries, determining how the structure is constrained by external forces. The time step parameter vector controls the magnitude of the time step in dynamic simulations and determines the accuracy and stability of the simulation in finite element analysis, especially when simulating non-static loads and time-varying processes.

[0061] A direct contact stiffness model is constructed using normal contact stiffness, normal projection matrix, tangential contact stiffness, and tangential projection matrix. Based on the contact conditions and material properties between layers, the stiffness of the contact surface in the normal direction is calculated, yielding the normal contact stiffness, which determines the contact surface's resistance to deformation under vertical pressure. High stiffness indicates that the contact surface is not easily deformed, while low stiffness means the contact surface is more prone to deformation. Based on the friction and relative sliding characteristics between the contact surfaces, the stiffness of the contact surface in the tangential direction is calculated, yielding the tangential contact stiffness. Tangential stiffness affects the frictional force between the contact surfaces, determining whether and to what extent the film layer slips. The normal and tangential projection matrices are used to decompose the forces and deformations on the contact surface into normal and tangential components, aiding in a more accurate analysis of the interactions between the contact surfaces in three-dimensional space. The direct contact stiffness model is a numerical model describing the mechanical properties of the direct contact interface. It defines the interface's resistance to deformation in the normal and tangential directions using four core parameters, accurately reflecting the stiffness characteristics of the contact surface.

[0062] Coulomb friction constraint is a numerical implementation method based on Coulomb's law of friction, used to limit the maximum shear stress between contact surfaces. Relative sliding is allowed when the interfacial shear stress reaches a certain critical value, which is determined by both the interfacial friction characteristics and the normal pressure. The constructed stress-deformation input vector set is integrated with a contact stiffness model specifically designed for direct contact characteristics. During the solution process, Coulomb friction constraint is introduced as an important boundary condition, determining 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 an adhesive state; when the shear stress reaches the critical value, the contact surface enters a sliding state.

[0063] An incremental iterative algorithm is employed, performing multiple iterations within each load increment step until the convergence criterion is met. First, the contact stiffness matrix is ​​calculated based on the current contact state. Then, the overall stiffness matrix is ​​assembled, the displacement increment is solved, and the stress state is updated based on the displacement increment. Finally, the contact state is checked for changes, and convergence is determined. Through this iterative process, three-dimensional stress-deformation direct simulation data is obtained. The stress-displacement increment records the deformation development at each calculation step, the interlayer normal stress reflects the degree of compaction at the contact surface, the interlayer shear stress shows the distribution of interfacial friction, and the stress field distribution provides a comprehensive view of the mechanical state of the entire structure. The stress-displacement increment represents the increment between stress and deformation caused by external loads or other influences; the interlayer normal stress describes the stress between layers in the normal direction, which is usually caused by external loads or compression between membrane layers; the interlayer shear stress describes the stress between layers in the tangential direction, which is usually related to friction, slip, and shear deformation; the stress field distribution represents the stress state at different locations throughout the geomembrane structure. The stress field distribution is obtained through finite element analysis and is used to analyze the stress response of the material. For example, the maximum displacement increment of stress displacement occurs in the 6th load step, with a value of 3.2 cm; the interlayer normal stress ranges from 35 to 650 kPa, reaching its 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 it is 5 to 8 kPa; in terms of stress field distribution, the maximum tensile stress of the geomembrane is 2.1 MPa, which occurs at the junction of the slope and the bank slope; 85% of the contact units are in an adhesive state, and 15% are in a sliding state.

[0064] By constructing a stress-deformation input vector set and combining it with normal and tangential contact stiffness models for accurate three-dimensional finite element analysis, the stress and deformation behavior of geomembrane structures under complex loads is simulated. Through rigorous Coulomb friction constraints and incremental iterative algorithms, the complex mechanical behavior of the contact surface, including the transition process between adhesion and slip states, is realistically reproduced, thereby improving the accuracy of the simulation results.

[0065] Furthermore, this application also includes the following steps: constructing an indirect contact weak coupling model, the indirect contact weak coupling model including a stress attenuation coefficient; introducing a material degradation coefficient to identify the equivalent material matrix of the indirect contact surface element; and performing a three-dimensional finite element solution on the stress deformation input vector set based on 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] Furthermore, this application also includes the following steps: 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.

[0067] Specifically, indirect contact surface units refer to areas where contact force transmission is weakened due to construction defects, environmental factors, etc., without direct contact. The indirect contact weak coupling model is used to simulate the contact behavior in these areas, considering the weakening effect of force transmission, such as material aging, contamination, or construction defects. Based on the specific defect types identified in the indirect contact units, the indirect contact weak coupling model is constructed. The indirect contact weak coupling model includes a stress attenuation coefficient, which describes the degree of stress reduction caused by various factors when the contact surface is subjected to force. The stress attenuation coefficient is determined by fitting three attenuation factors: defect attenuation factor, mass attenuation factor, and contamination attenuation factor. The defect attenuation factor describes the degree of stress attenuation caused by construction defects in the geomembrane structure; a high defect attenuation factor indicates significant defects in the contact surface, leading to a significant reduction in force transmission. The mass attenuation factor represents the degree of stress attenuation caused by substandard construction quality in the geomembrane structure; a high mass attenuation factor weakens the structure's load-bearing capacity. The contamination attenuation factor describes the impact of environmental pollution on the geomembrane material, leading to a decrease in its mechanical properties; a high contamination attenuation factor reduces the friction and compressive strength between the contact surfaces. For example, a slope area with continuous folds of average height 2.5 cm has an attenuation factor of 0.65 determined through calibration tests; the compaction degree of this slope area is 82%, corresponding to a mass attenuation factor of 0.75; chemical analysis shows a pH of 3.5 and excessive heavy metal ion concentration, determining a pollution attenuation factor of 0.45. Using a weighted geometric mean with weights of 0.4, 0.3, and 0.3, the stress attenuation coefficient is calculated to be 0.59.

[0068] A material degradation coefficient is introduced to identify the equivalent material matrix of indirect contact surface elements. The material degradation coefficient reflects the degree of material degradation over time or due to environmental changes, typically considering the influence of the external environment on geomembrane materials. In finite element analysis, the equivalent material matrix describes the changes in the mechanical properties of the material after degradation, contamination, defects, etc. The equivalent material matrix simplifies complex mechanical calculations, allowing it to better reflect material behavior under real-world conditions. For areas with chemical contamination or aging, the material degradation coefficient will correspondingly reduce the material's elastic modulus and strength parameters. Establishing an equivalent material matrix based on the degraded material parameters simplifies the complex multilayer system into a single-layer medium with equivalent mechanical properties.

[0069] The stress attenuation coefficient and equivalent material matrix are integrated into the finite element solver to perform three-dimensional finite element analysis on the stress-deformation input vector set. During the solution process, the stress attenuation coefficient is used to adjust the stress transfer capability of the contact surface, while the equivalent material matrix is ​​used to calculate the stiffness contribution of the elements. A joint simulation method combining the equivalent medium model and the multi-field coupled attenuation model is used to obtain three-dimensional stress-deformation indirect simulation data, accurately reflecting the special mechanical behavior of interfaces with defects. Multi-field coupling refers to the mutual influence between different physical fields, such as stress, temperature, and humidity fields. The interaction of these physical fields leads to the attenuation of stress at the contact surface. The multi-field coupled attenuation model is a mathematical model describing the stress change at the geomembrane contact surface under the combined effect of these multiple factors. Specifically, the equivalent stress of the contact surface is calculated based on the current stress state and stress attenuation coefficient; the system stiffness matrix is ​​assembled based on the equivalent material matrix; then, the equilibrium equations are solved to obtain the displacement increment; the stress-strain state is updated based on the displacement increment; and finally, the convergence of the system is checked. This process is iterated until a preset convergence criterion is met. Three-dimensional stress-deformation indirect simulation data are calculated mechanical response data, including full-field distribution information such as stress, strain, and displacement, taking into account the influence of interface defects.

[0070] By constructing an indirect contact weakly coupled model, the mechanical behavior of the membrane layer contact surface is simulated, taking into account the influence of factors such as defects, quality, and pollution on the geomembrane structure. A material degradation coefficient and an equivalent material matrix are introduced to ensure that the simulation accurately reflects the degradation and weakening phenomena that occur in the material during actual use. More accurate stress and deformation data are obtained through three-dimensional finite element analysis.

[0071] 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.

[0072] Furthermore, this application also includes the following steps: defining direct-fusion smoothing weights, indirect-fusion smoothing weights, and hybrid-fusion smoothing weights; using the direct contact surface units and indirect contact surface units as fusion nodes, fusing according to the direct-fusion smoothing weights, indirect-fusion smoothing weights, and hybrid-fusion smoothing weights to obtain three-dimensional stress-deformation fusion simulation data.

[0073] Specifically, based on element type and spatial relationship, three different fusion smoothing weights are defined. The direct-fusion smoothing weight is a weighting coefficient used in the data fusion process between directly contacting elements. During simulation, the fusion of directly contacting elements is adjusted according to this smoothing weight. A larger smoothing weight means a stronger degree of fusion between directly contacting elements, resulting in more uniform deformation and stress transfer. The indirect-fusion smoothing weight is a weighting coefficient used in the data fusion process between indirectly contacting elements, used to maintain a reasonable smooth transition when processing defect area data. Indirectly contacting elements typically involve attenuation effects, therefore, the setting of their smoothing weight affects the degree of stress transfer and deformation smoothness between contact surfaces. The hybrid-fusion smoothing weight is a weighting coefficient used in the data fusion process between directly contacting elements and indirectly contacting elements, used to handle the data connection problem between different types of elements. Hybrid fusion is usually the most complex in practical applications because these two types of elements have significant differences in mechanical behavior; the hybrid smoothing weight determines the transition method from direct contact to indirect contact.

[0074] The direct-fusion smoothing weight is set to a high value, such as 0.8, to ensure strict continuity of data within the good contact area; the indirect-fusion smoothing weight is set to a moderate value, such as 0.5, to maintain reasonable smoothness within the defect area while considering the impact of defects; the hybrid-fusion smoothing weight adopts a distance-based gradient function, such as 0.6, to form a smooth transition zone at the boundary between the direct and indirect contact areas.

[0075] A fusion network covering the entire contact surface is constructed using the boundary nodes of all directly contacting and indirectly contacting surface elements as fusion nodes. At each fusion node, an appropriate weight coefficient is automatically selected based on the element type combination in which the node is located. For nodes located inside the directly contacting region, a direct-fusion smoothing weight is used; for nodes located inside the indirectly contacting region, an indirect-fusion smoothing weight is used; and for nodes located at the boundary between the two types of regions, a hybrid-fusion smoothing weight is used.

[0076] In the specific fusion calculation, the contributions of direct and indirect contact units are weighted and summed according to their corresponding smoothing weights to ensure a smooth transition of stress and deformation results between contact surface units in the transition region. The fusion process includes stress-displacement smoothing, weighting of normal and tangential stresses, and adjustment of the continuity of shear stress and stress field distribution. The three-dimensional stress-deformation fusion simulation data is the final stress and deformation distribution data obtained by integrating direct and indirect contact simulation data after the fusion of each contact surface unit. It accurately reflects the stress distribution, deformation, and response behavior of the entire geomembrane structure under different conditions. For example, a direct-fusion smoothing weight of 0.8 indicates relatively uniform stress transmission; an indirect-fusion smoothing weight of 0.5 indicates a strong attenuation effect and poor contact stress transmission; and a mixed-fusion smoothing weight of 0.6 is used to smooth the transition of stress and deformation changes. The final three-dimensional stress-deformation fusion simulation data includes: a stress displacement increment of 1.2 mm, a normal stress of 200 kPa, and a shear stress of 120 kPa in the direct contact area; a stress displacement increment of 2.5 mm, a normal stress of 80 kPa, and a shear stress of 50 kPa in the indirect contact area; the fused stress field distribution map shows that the stress in the contact area transitions smoothly in the transition region, avoiding abrupt stress changes.

[0077] By defining and fusing smoothing weights for direct, indirect, and mixed contact elements, the stress and deformation behavior of geomembrane structures under different contact states can be effectively simulated. Weighted fusing of stress and deformation data avoids abrupt changes in stress and deformation in transition regions, improving the accuracy and reliability of the simulation results.

[0078] In summary, the three-dimensional stress-deformation simulation and analysis method for geomembrane structures provided in this application has the following technical advantages:

[0079] The geomembrane structure is parametrically modeled to construct layered geomembrane modeling data. Based on this data, interlayer contact surfaces are identified, and these surfaces are discretized according to a preset grid to obtain a finite number of contact surface units. Interlayer contact defects, construction quality, and environmental pollution status corresponding to these contact surface units are collected. These factors are used as hierarchical features to classify 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 and indirect contact surface units, obtaining three-dimensional stress-deformation direct simulation data and three-dimensional stress-deformation indirect simulation data. The three-dimensional stress-deformation direct and indirect simulation data are then fused to obtain the three-dimensional stress-deformation fused simulation data of the geomembrane structure. In other words, by constructing geomembrane layer modeling data through parametric modeling, the interlayer contact surfaces are discretized through a preset mesh. Contact units or interface units are used to simulate complex interactions such as friction, slippage, and voids between the geomembrane and the overlying protective layer and the underlying support layer. Interlayer contact defects, construction quality, and environmental pollution status are collected as grading standards. The contact surface units are divided into direct contact surface units and indirect contact surface units. When simulating contact stress deformation, a direct plus indirect fusion method is used to analyze the deformation simulation state between layers, thereby improving the accuracy of the simulation analysis data.

[0080] Example 2: Based on the same inventive concept as the three-dimensional stress-deformation simulation analysis method for geomembrane structures in Example 1, this application also provides a three-dimensional stress-deformation simulation analysis system for geomembrane structures. Please refer to the appendix. Figure 2 The three-dimensional stress-deformation simulation and analysis system for geomembrane structures includes:

[0081] The system comprises the following modules: a parametric modeling module 11 for constructing layered geomembrane modeling data; a discretization module 12 for identifying interlayer contact surfaces based on the layered geomembrane modeling data and discretizing the interlayer contact surfaces according to a preset grid to obtain a finite number of contact surface units; a contact surface unit division module 13 for collecting interlayer contact defects, interlayer construction quality, and environmental pollution status corresponding to the finite number of contact surface units, and dividing 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 hierarchical features; a simulation analysis module 14 for collecting stress-deformation load data to simulate the direct and indirect contact surface units to obtain three-dimensional stress-deformation direct simulation data and three-dimensional stress-deformation indirect simulation data; and a data fusion module 15 for fusing the three-dimensional stress-deformation direct simulation data and three-dimensional stress-deformation indirect simulation data to obtain three-dimensional stress-deformation fused simulation data of the geomembrane structure.

[0082] Furthermore, the discretization module 12 in the three-dimensional stress-deformation simulation analysis system for geomembrane structures is also used to: identify the stress contact region and non-stress contact region based on the interlayer contact surface under the stress-deformation load data; discretize the stress contact region using a first preset grid to obtain a first set of contact surface units; discretize the non-stress contact region using a second preset grid to obtain a second set of contact surface units, wherein the step size of the first preset grid is smaller than the step size of the second preset grid; and obtain a finite number of contact surface units based on the first set of contact surface units and the second set of contact surface units.

[0083] Furthermore, the discretization module 12 in the three-dimensional stress-deformation simulation analysis system for geomembrane structures is also used to: 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 a transition contact region between the stress contact region and the non-stress contact region; discretize the stress transition contact region 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; and obtain a finite number of contact surface elements based on the first set of contact surface elements, the second set of contact surface elements, and the third set of contact surface elements.

[0084] Furthermore, the contact surface unit division module 13 in the three-dimensional stress-deformation simulation analysis system for geomembrane structures is also used to: calculate the contribution score of the finite number of contact surface units for 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; perform weighted calculation on the defect contribution score, construction contribution score, and pollution contribution score by configuring contribution weights to obtain the comprehensive contribution score of the finite number of contact surface units; classify contact surface units with scores less than a preset score threshold as direct contact surface units; and classify contact surface units with scores greater than or equal to the preset score threshold as indirect contact surface units.

[0085] Furthermore, the contact surface unit partitioning module 13 in the three-dimensional stress-deformation simulation analysis system for geomembrane structures is also used to: construct the equivalent thin-layer mechanical parameters of the indirect contact surface unit; simulate the stress change influence degree based on the equivalent thin-layer mechanical parameters according to the stress-deformation load data, wherein the stress change influence degree is the data change degree between the stress-deformation output data obtained from the equivalent simulation and the stress-deformation load data; and re-mark the indirect contact surface units with a stress change influence degree less than a preset influence degree as direct contact surface units.

[0086] Furthermore, the simulation analysis module 14 in the three-dimensional stress-deformation simulation analysis system for geomembrane structures is also used to: construct a stress-deformation input vector set based on the stress-deformation load data and the state input data of the direct contact surface element, including external load vector, material parameter vector, boundary condition vector, and time step parameter vector; construct a direct contact stiffness model through normal contact stiffness, normal projection matrix, tangential contact stiffness, and tangential projection matrix; and use Coulomb friction constraints to perform three-dimensional finite element solution on the stress-deformation input vector set according to 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.

[0087] Furthermore, the simulation analysis module 14 in the three-dimensional stress-deformation simulation analysis system for geomembrane structures is also used 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 the equivalent material matrix of the indirect contact surface element; and perform three-dimensional finite element solution 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] Furthermore, the simulation analysis module 14 in the three-dimensional stress-deformation simulation analysis system for geomembrane structures is also used for: 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 pollution attenuation factor.

[0089] Furthermore, the data fusion module 15 in the three-dimensional stress-deformation simulation analysis system for geomembrane structures is also used to: define direct-fusion smoothing weights, indirect-fusion smoothing weights, and mixed-fusion smoothing weights; and use the direct contact surface unit and the indirect contact surface unit as fusion nodes to perform fusion according to the direct-fusion smoothing weights, indirect-fusion smoothing weights, and mixed-fusion smoothing weights to obtain three-dimensional stress-deformation fusion simulation data.

[0090] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on the differences from other embodiments. The three-dimensional stress-deformation simulation analysis method and specific examples for geomembrane structures in the aforementioned embodiment one are also applicable to the three-dimensional stress-deformation simulation analysis system for geomembrane structures in this embodiment. Through the foregoing detailed description of the three-dimensional stress-deformation simulation analysis method for geomembrane structures, those skilled in the art can clearly understand the three-dimensional stress-deformation simulation analysis system for geomembrane structures in this embodiment. Therefore, for the sake of brevity, it will not be described in detail here.

[0091] The above description of the disclosed embodiments enables those skilled in the art to make or use this application. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of this application. Therefore, this application is not to be limited to the embodiments shown herein, but is to be accorded 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 this application without departing from the spirit and scope of this application. Therefore, if such modifications and variations fall within the scope of this application and its equivalents, this application also intends to include such modifications and variations.

Claims

1. A method for simulating and analyzing three-dimensional stress deformation of a geomembrane structure, characterized by, 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 of three-dimensional stress and deformation and the indirect simulation data of three-dimensional stress and deformation, the three-dimensional stress and deformation fused simulation data of the geomembrane structure is obtained. 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; 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.

2. The method for three-dimensional stress deformation simulation analysis of a geomembrane structure according to claim 1, wherein 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 method for three-dimensional stress deformation simulation analysis of a geomembrane structure according to claim 2, wherein 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 method for three-dimensional stress deformation simulation analysis of geomembrane structure according to claim 1, wherein, 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.

5. The three-dimensional stress-deformation simulation and analysis method for geomembrane structures as described in claim 4, 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.

6. The three-dimensional stress-deformation simulation and analysis method for geomembrane structures as described in claim 5, 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.

7. The method for three-dimensional stress deformation simulation analysis of geomembrane structure according to claim 1, wherein, 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.

8. 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 7, 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