Cement mortar loss simulation method and system based on tidal action

Abstract

The invention discloses a cement mortar loss simulation method and system based on tidal action, and relates to the technical field of engineering material durability research, and the method comprises the following steps: S1, obtaining tide simulation parameters and cement mortar structure parameters input by a user, the tide simulation parameters being used for defining a simulated ocean hydrodynamic environment, and the cement mortar structure parameters being used for simulating the ocean hydrodynamic environment; the cement mortar structure parameters are used for defining to-be-analyzed engineering structure physical attributes; and S2, based on the parameters of the cement mortar structure, generating a three-dimensional geometric model of the cement mortar structure and a surrounding soil body in a preset coordinate system. The method comprises the following steps: acquiring tide and cement mortar structure parameters, constructing a three-dimensional geometric model and a finite element grid, establishing dynamic boundary conditions in combination with tide and wave characteristics, acquiring pore water pressure distribution by using transient hydrodynamic analysis, quantifying mass loss through a loss calculation model, and dynamically updating material attributes to simulate coupling influence.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of durability research of engineering materials, in particular to a cement mortar loss simulation method and system based on tidal action. BACKGROUND

[0002] In coastal engineering, cement mortar is often used for cementation of structures such as foundation pit support and embankment reinforcement, and its stability is directly related to engineering safety. The periodic rise and fall of seawater, wave impact and pore water pressure change caused by tidal action can cause cement mortar to lose due to penetration, corrosion and other reasons, resulting in a decrease in structural strength and even accidents such as collapse, which is a problem that needs to be focused on in coastal engineering.

[0003] However, current research on cement mortar loss under tidal action relies on field sampling or simple simulation, and lacks a special simulation method and system that can accurately reproduce the dynamic characteristics of tides, such as water level cycles, wave patterns and pore water pressure fluctuations, making it difficult to quantify the impact of different tidal parameters on loss rate and range, resulting in a lack of scientific basis for predicting and preventing loss in engineering, and a high safety risk. SUMMARY

[0004] To overcome the deficiencies of the prior art, the present application provides a cement mortar loss simulation method and system based on tidal action, which solves the problem that the prior art usually relies on field sampling or simple simulation when used relatively, and lacks a method that can accurately reproduce the dynamic characteristics of tides.

[0005] To achieve the above purpose, the present application is implemented by the following technical scheme: a cement mortar loss simulation method based on tidal action, comprising: S1, obtaining tidal simulation parameters and cement mortar structure parameters input by a user, the tidal simulation parameters being used to define the simulated marine hydrodynamic environment, and the cement mortar structure parameters being used to define the physical properties of the engineering structure to be analyzed; S2, generating a three-dimensional geometric model of the cement mortar structure and its surrounding soil in a preset coordinate system based on the cement mortar structure parameters, and performing spatial discretization processing on the three-dimensional geometric model to obtain a finite element calculation grid containing multiple elements and nodes; S3, based on the tidal simulation parameters, establishing boundary conditions for describing hydrodynamic changes, and applying the boundary conditions to the fluid domain boundary of the finite element calculation grid; S4, calling a hydrodynamic calculation module to perform transient hydrodynamic analysis on the finite element calculation grid to solve and obtain the water level, wave-induced hydrodynamic pressure at each node and the pore water pressure distribution inside the soil body within the entire simulation time; S5, a cement mortar loss calculation model is established, based on the obtained pore water pressure distribution, hydraulic gradient and seepage velocity in each unit at each time step are calculated, and according to a preset loss criterion function, the mass loss rate of the cement mortar is determined; S6, the total mass loss amount of each unit in the entire simulation duration is cumulatively calculated, and the simulation results are processed to generate a visual output file representing the loss range and degree of the cement mortar.

[0006] Further, the tidal simulation parameters include tidal type, tidal period, tidal amplitude, initial water level, wave type, wave height, wave period and wave direction angle; the cement mortar structure parameters include elastic modulus, Poisson's ratio, density, permeability coefficient, porosity of the cement mortar and contact properties of the cement mortar and soil contact surface.

[0007] Further, the boundary condition for describing hydrodynamic changes is established, and the boundary condition is applied to the fluid domain boundary of the finite element calculation grid, including: Based on the tidal type, tidal period and tidal amplitude, a tidal function for representing periodic water level changes is generated, and is applied to the far-field water area boundary of the finite element calculation grid as a time-varying pressure boundary or water level boundary; Based on the wave type, wave height and wave period, a linear wave theory or a nonlinear wave theory is selected, a wave velocity field and a pressure field function are generated, and are superimposed on the tidal function to serve as input conditions for the far-field water area boundary.

[0008] Further, the hydrodynamic calculation module is called to obtain the pore water pressure distribution inside the soil, including: On the finite element calculation grid, control equations based on unsaturated porous medium seepage theory are established; The hydrodynamic pressure acting on the structure surface by the tide and wave is applied to the solid-liquid coupling interface of the finite element calculation grid as a transient load; Using a time stepping algorithm, the control equations are iteratively solved to obtain the pore water pressure value of each node in the finite element calculation grid at each time step, forming a spatiotemporal distribution dataset of the pore water pressure.

[0009] Further, the loss criterion function is a function of the seepage velocity and the critical starting flow velocity of the cement mortar, and the mass loss rate of the cement mortar is determined, including: For any unit in the finite element calculation grid that belongs to the cement mortar, when the calculated seepage velocity inside the unit exceeds the preset critical starting flow velocity threshold, it is determined that the unit has lost; Based on the difference of permeation velocity exceeding the threshold of critical start-up flow rate, the mass loss rate components caused by chemical dissolution and physical erosion are calculated respectively through preset dissolution model and erosion model; The chemical dissolution rate component and the physical erosion rate component are linearly superimposed to obtain the total mass loss rate of the unit at the current time step.

[0010] Further, the simulation method further comprises: After the calculation of each time step is completed, the material properties of the mass loss unit are dynamically updated according to the calculated cement mortar mass loss amount; The updated material properties include reducing the density and elastic modulus of the unit, and increasing the porosity and permeability coefficient of the unit; The updated material properties are used as the initial conditions for the next time step calculation to realize the coupled influence simulation of the loss process on the physical properties of the structure and the hydraulic response.

[0011] Further, when the cumulative mass loss amount of any unit reaches a preset structure failure percentage, the unit is removed from the set of effective bearing units and is treated as a fluid unit in subsequent calculations to simulate the change of water flow path after local collapse of the structure.

[0012] Further, the generation of a visual output file representing the loss range and degree of cement mortar includes: A three-dimensional cloud chart of the cumulative loss amount of cement mortar is generated, and the three-dimensional cloud chart uses different colors to indicate the spatial distribution and numerical size of the loss amount; A curve graph of the loss amount of cement mortar at key positions over time is generated for analyzing the loss development history of a specific area; An analysis report of the overall strength, stability safety factor of the structure before and after simulation is output, and the report contains a quantitative safety factor reduction value.

[0013] Further, the method further comprises: A plurality of different sets of tidal simulation parameters are set as comparative analysis working conditions; For each set of comparative analysis working conditions, the above-mentioned steps are repeatedly executed to obtain corresponding simulation results; The simulation results under different working conditions are compared to quantitatively analyze the influence law of single parameter changes such as tidal period, tidal amplitude, wave height on the loss rate and the total loss amount of cement mortar, and a parameter sensitivity analysis atlas is generated.

[0014] The application also provides a cement mortar loss simulation system based on tidal action, characterized by comprising: a parameter and model definition module configured to obtain user-inputted tidal simulation parameters and cement mortar structure parameters, and generate a three-dimensional geometric model and a finite element calculation grid based on the cement mortar structure parameters; a hydrodynamic boundary application module configured to establish boundary conditions for describing hydrodynamic changes based on the tidal simulation parameters, and apply the boundary conditions to a fluid domain boundary of the finite element calculation grid; a hydrodynamic analysis module configured to call a hydrodynamic calculation module to perform transient hydrodynamic analysis on the finite element calculation grid to obtain water level, hydrodynamic pressure at each node, and pore water pressure distribution inside the soil body; a loss coupling calculation module configured to establish a cement mortar loss calculation model, calculate hydraulic gradient and seepage velocity based on the obtained pore water pressure distribution, and determine and accumulate the mass loss rate and total mass loss amount of the cement mortar according to a preset loss criterion function; a result generation and display module configured to process simulation results to generate visual output files and analysis reports representing the loss range, degree and structure safety factor changes of the cement mortar.

[0015] Advantages The present application can accurately reproduce the cement mortar loss process under the action of tides, quantify the influence of different tidal parameters on loss, and provide scientific basis for engineering loss prediction and prevention, and reduce safety risks. BRIEF DESCRIPTION OF DRAWINGS

[0016] Figure 1 The method flowchart of the present application; Figure 2 The system structure diagram of the present application. DETAILED DESCRIPTION

[0017] The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor fall within the scope of protection of the present application.

[0018] Please refer to Figure 1The application provides a cement mortar loss simulation method based on tidal action, comprising: S1, acquiring tidal simulation parameters input by a user and cement mortar structure parameters, the tidal simulation parameters being used to define a simulated marine hydrodynamic environment, and the cement mortar structure parameters being used to define physical properties of an engineering structure to be analyzed; S2, generating a three-dimensional geometric model of the cement mortar structure and its surrounding soil body in a preset coordinate system based on the cement mortar structure parameters, and performing spatial discretization processing on the three-dimensional geometric model to obtain a finite element calculation grid containing multiple units and nodes; S3, establishing boundary conditions for describing hydrodynamic changes based on the tidal simulation parameters, and applying the boundary conditions to a fluid domain boundary of the finite element calculation grid; S4, calling a hydrodynamic calculation module to perform transient hydrodynamic analysis on the finite element calculation grid to obtain water level, wave-induced hydrodynamic pressure at each node and pore water pressure distribution in the soil body within the entire simulation time length; S5, establishing a cement mortar loss calculation model, calculating hydraulic gradient and seepage velocity in each unit at each time step based on the obtained pore water pressure distribution, and determining the mass loss rate of the cement mortar according to a preset loss criterion function; S6, cumulatively calculating the total mass loss amount of each unit within the entire simulation time length, and processing the simulation results to generate a visual output file representing the loss range and degree of the cement mortar.

[0019] Specifically, in actual operation, the user can input the tidal simulation parameters and the cement mortar structure parameters through a dedicated input interface. For example, the tidal simulation parameters can be obtained from historical data of a marine observation station and input, and the cement mortar structure parameters can be determined according to engineering design drawings and material test reports and then input.

[0020] When generating the three-dimensional geometric model, a commonly used finite element modeling software such as ANSYS, ABAQUS, etc. can be used to construct the three-dimensional geometric model of the structure and its surrounding soil body in a preset Cartesian coordinate system according to the cement mortar structure parameters. When performing spatial discretization processing on the three-dimensional geometric model, tetrahedral units can be used for division, and the unit size is determined according to the complexity of the structure and the calculation accuracy requirement, for example, the unit size can be set to 0.1 m for a complex structure region, and the unit size can be set to 1 m for a simple surrounding soil body region, so as to obtain the finite element calculation grid containing multiple units and nodes.

[0021] When the boundary conditions are established, according to the tidal period, tidal amplitude and the like in the tidal simulation parameters, corresponding functions are generated by programming to describe the water power change, and the boundary conditions represented by the functions are applied to the fluid domain boundary of the finite element calculation grid through the boundary condition application function of the finite element software.

[0022] When the water power calculation module is called, the relevant algorithm of transient water power analysis is built in the module, and according to the applied boundary conditions, iterative calculation is performed on the finite element calculation grid, so that the water level, hydrodynamic pressure at each node and the pore water pressure distribution in the soil body within the entire simulation time are solved.

[0023] When the cement mortar loss calculation model is established, based on the pore water pressure distribution, the hydraulic gradient is calculated through the hydraulic gradient calculation formula The hydraulic gradient is calculated, wherein The hydraulic gradient is calculated, wherein The hydraulic gradient is calculated, wherein The permeation velocity is calculated, wherein The permeation velocity is calculated, wherein The permeation velocity is calculated, wherein

[0024] When the total mass loss amount is accumulated and calculated, the mass loss rate at each time step is multiplied by the time step, and then the results within the entire simulation time are accumulated. When the simulation results are processed, the three-dimensional cloud chart, curve chart and the like can be generated by using the visualization software to generate the visualization output file.

[0025] Through the embodiment, the loss process of the cement mortar under the action of the tide can be accurately reproduced, and a scientific basis for the prediction and prevention of the loss in the engineering is provided.

[0026] The application further provides tidal simulation parameters, the tidal simulation parameters including a tidal type, a tidal period, a tidal amplitude, an initial water level, a wave type, a wave height, a wave period and a wave direction angle; and cement mortar structure parameters including an elastic modulus, a Poisson's ratio, a density, a permeation coefficient, a porosity of the cement mortar and a contact property of a contact surface between the cement mortar and the soil body.

[0027] Specifically, the tide type can be divided into semi-diurnal tide, diurnal tide, etc., for example, a certain coastal area is semi-diurnal tide, then input this type. The tide period is the time of one tide rise and fall, for example, the period of semi-diurnal tide is about 12 hours. The tide amplitude is the maximum water level difference of tide rise and fall, for example, the tide amplitude of a certain area is 2m. The initial water level is the water level at the beginning of simulation, which can be determined according to the actual water level at a certain time, for example, the initial water level is 1m. The wave type includes wind wave, swell, etc., the wave height is the height of the wave, for example, the wave height is 1.5m. The wave period is the time of one wave fluctuation, for example, the wave period is 6s. The wave direction angle is the angle between the wave propagation direction and the north direction, for example, the wave direction angle is 30 degrees.

[0028] The elastic modulus of cement mortar is generally between 20-40GPa, for example, the elastic modulus of a certain type of cement mortar is 30GPa. The Poisson's ratio is usually between 0.2-0.3, for example, 0.25. The density is about 2400-2500kg / m 3 , for example, 2450kg / m 3 . The permeability coefficient is determined according to the density of the cement mortar, for example, 1×10 -8 m / s. The porosity is generally between 10%-20%, for example, 15%. The contact properties of the cement mortar and the soil contact surface include friction coefficient, bond strength, etc., for example, the friction coefficient is 0.3, and the bond strength is 1MPa.

[0029] By specifying the specific values of these parameters, the simulated marine hydrodynamic environment and the physical properties of the engineering structure can be more accurately defined, and the simulation accuracy can be improved.

[0030] The application further proposes to establish boundary conditions for describing hydrodynamic changes, and apply the boundary conditions to the fluid domain boundary of the finite element calculation grid, including: Based on the tide type, tide period and tide amplitude, a tide function is generated to represent the periodic water level change, and is applied as a time-varying pressure boundary or water level boundary to the far-field water area boundary of the finite element calculation grid; Based on the wave type, wave height and wave period, a linear wave theory or a nonlinear wave theory is selected to generate a wave velocity field and a pressure field function, and the wave velocity field and the pressure field function are superimposed on the tide function to serve as input conditions for the far-field water area boundary.

[0031] Specifically, when generating the tide function, for semi-diurnal tide, according to the tide period and the tide amplitude , it can be expressed as , wherein is the water level at time , h is the initial water level, Time. The tidal function is applied to the far-field water area boundary of the finite element calculation grid as the water level boundary varying with time, for example, in the finite element software, by setting the water level of the boundary varying with time as the function.

[0032] When the wave type is wind wave and the wave height is small, the linear wave theory can be selected to generate the wave velocity field and the pressure field function. The wave velocity field function can be expressed as: wherein is the horizontal velocity along the wave direction, is the gravity acceleration, taking 9.8 m / s 2 , is the wave height, is the wave circular frequency, is the wave period, is the wave number, is the wave length, is the vertical coordinate, is the water depth, is the horizontal coordinate along the wave direction. The wave pressure field function can be expressed as wherein is the wave pressure, is the seawater density, taking 1025 kg / m 3 .

[0033] The wave velocity field and the pressure field function are superimposed on the tidal function, for example, in the input condition of the far-field water area boundary, while the water level change caused by the tide and the velocity and pressure changes caused by the wave are taken into account together as the boundary condition.

[0034] Through the embodiment, the change of the hydrodynamic force can be more accurately described, and the simulated marine hydrodynamic environment is closer to the actual situation.

[0035] The application further proposes to call the hydrodynamic calculation module to solve and obtain the pore water pressure distribution inside the soil body, comprising: establishing a control equation based on the unsaturated porous medium seepage theory on the finite element calculation grid; the hydrodynamic pressure acting on the structure surface by the tide and the wave is taken as a transient load and applied to the solid-liquid coupling interface of the finite element calculation grid; a time stepping algorithm is adopted to iteratively solve the control equation, so as to obtain the pore water pressure value of each node in the finite element calculation grid at each time step, and form a time-space distribution data set of the pore water pressure.

[0036] Specifically, the control equation based on the unsaturated porous medium seepage theory can be expressed as: wherein is the porosity, is the saturation, is the density of water, is the seepage velocity, is the time. The equation describes the flow rule of water in the porous medium.

[0037] The hydrodynamic pressure generated by the joint action of tides and waves is applied as a transient load to the solid-liquid coupling interface, for example, by the load application function of the finite element software, the calculated hydrodynamic pressure values at different times are applied to the corresponding interface nodes.

[0038] The Newmark-β method in the time stepping algorithm is used to iteratively solve the control equation, and the time step can be determined according to the tidal period and the wave period, for example, 0.1s. In each time step, the pore water pressure value of each node in the finite element calculation grid at this time is obtained by iterative calculation, and the pore water pressure values at different time steps are integrated to form a time-space distribution data set of the pore water pressure.

[0039] In this way, the pore water pressure distribution inside the soil body can be accurately solved, which lays a foundation for subsequent calculation of cement mortar loss.

[0040] The application further proposes that the loss criterion function is a function of the seepage velocity and the critical starting flow velocity of the cement mortar, and the mass loss rate of the cement mortar is determined, comprising: For any unit in the finite element calculation grid that belongs to the cement mortar, when the calculated seepage velocity inside the unit exceeds the preset critical starting flow velocity threshold, it is determined that the unit has loss; Based on the seepage velocity difference exceeding the critical starting flow velocity threshold, the mass loss rate components caused by chemical dissolution and physical erosion are calculated through the preset dissolution model and the erosion model, respectively; The chemical dissolution rate component and the physical erosion rate component are linearly superimposed to obtain the total mass loss rate of the unit at the current time step.

[0041] Specifically, the preset critical starting flow velocity threshold can be determined according to the material properties and experimental data of the cement mortar, for example, 0.01m / s. When the calculated seepage velocity of a unit is 0.015m / s, which exceeds the threshold, it is determined that the unit has loss.

[0042] The dissolution model can be represented as wherein is the mass loss rate component caused by chemical dissolution, is the chemical dissolution coefficient, which can be determined by experiment, for example, 0.002kg / (m•s), is the seepage velocity, is the critical starting flow velocity threshold. The erosion model can be represented as wherein is the mass loss rate component caused by physical erosion, is the physical erosion coefficient, for example, 0.003 kg / (m•s).

[0043] Total mass loss rate As in the above example, then ,

[0044] Total mass loss rate .

[0045] The present application further proposes that the simulation method further comprises: After the calculation of each time step is completed, the material properties of the unit where mass loss occurs are dynamically updated according to the calculated mass loss of the cement mortar; The updated material properties include reducing the density and elastic modulus of the unit, and increasing the porosity and permeability coefficient of the unit; The updated material properties are used as the initial conditions for the calculation of the next time step, realizing the simulation of the coupled influence of the mass loss process on the physical properties of the structure and its hydraulic response.

[0046] Specifically, after the calculation of each time step is completed, for the unit where mass loss occurs, the material properties are dynamically updated according to the mass loss thereof. For example, when the mass loss of a unit reaches 10% of the initial mass of the unit, the density thereof is reduced by 10%, from the original 2450 kg / m 3 to 2205 kg / m 3 ; the elastic modulus is reduced by 20%, from 30 GPa to 24 GPa; the porosity is increased by 20%, from 15% to 18%; and the permeability coefficient is increased by 50%, from 1×10 -8 m / s to 1.5×10 -8 m / s.

[0047] These updated material properties are input into the calculation of the next time step as initial conditions. In this way, during the simulation process, the physical properties of the structure will change continuously with the mass loss process, thereby affecting the hydraulic response, realizing the simulation of the coupled influence of the two.

[0048] The present application further proposes that when the cumulative mass loss of any unit reaches a preset percentage of the failure of the structure, the unit is removed from the set of effective load-bearing units, and in subsequent calculations, it is treated as a fluid unit, to simulate the change in flow path after the local collapse of the structure.

[0049] Specifically, the preset structural failure percentage can be determined according to engineering safety requirements and structural characteristics of the cement mortar, for example, 50%. When the cumulative mass loss of a unit reaches 50% of its initial mass, the unit is removed from the set of effective load-bearing units.

[0050] In subsequent calculations, the unit is regarded as a fluid unit, and its material properties adopt the properties of a fluid, such as a density of 1025 kg / m 3 , and has no load-bearing capacity. In this way, the flow path of the water flow in the area changes, simulating the change in the water flow path after the local collapse of the structure.

[0051] In this way, the influence of the local collapse of the structure on the entire water flow field can be more accurately simulated, improving the authenticity of the simulation.

[0052] The application further proposes to generate visual output files representing the loss range and degree of the cement mortar, including: generate a three-dimensional cloud chart of the cumulative loss of the cement mortar, and the three-dimensional cloud chart uses different colors to indicate the spatial distribution and numerical size of the loss; generate a curve graph of the loss of the cement mortar at key positions over time, for analyzing the development history of the loss in a specific area; output a comparative analysis report of the overall strength and stability safety factor of the structure before and after simulation, and the report contains quantitative safety factor reduction values.

[0053] Specifically, when generating a three-dimensional cloud chart, a visualization software such as Tecplot can be used to assign different colors to different loss ranges according to the cumulative loss data of each unit, for example, a loss of 0-10 kg is represented by blue, a loss of 10-20 kg is represented by green, and a loss of 20 kg or more is represented by red, directly displaying the spatial distribution and numerical size of the loss.

[0054] The key positions can be selected from weak parts of the structure, contact parts with the soil, etc., such as corners of the structure. The curve graph of the loss at these positions over time is generated, with the abscissa representing time and the ordinate representing loss, and the trend and development history of the loss in a specific area over time can be clearly analyzed through the curve graph.

[0055] When outputting the comparative analysis report, the overall strength of the structure, such as compressive strength and tensile strength, and the stability safety factor before and after simulation are calculated. For example, the stability safety factor of the structure before simulation is 3.0, and the stability safety factor after simulation is 2.0, the safety factor reduction value is 1.0, and the changes of these quantitative indicators are described in detail in the report.

[0056] By generating these visual output files and analysis reports, the simulation results can be more intuitively displayed, providing convenience for engineers to understand the loss and make decisions.

[0057] The application further proposes that the method further comprises: Setting multiple groups of different tidal simulation parameters as comparative analysis conditions; For each group of comparative analysis conditions, repeating the above steps to obtain corresponding simulation results; Comparing the simulation results under different conditions, quantitatively analyzing the influence of single parameters such as tidal period, tidal amplitude, wave height on the cement mortar loss rate and the final loss amount, and generating a parameter sensitivity analysis atlas.

[0058] Specifically, when setting the comparative analysis conditions, for example, three groups of different tidal periods are set, which are 12h, 18h and 24h respectively, and other tidal simulation parameters remain unchanged; three groups of different tidal amplitudes are set, which are 1m, 2m and 3m respectively, and other parameters remain unchanged; and three groups of different wave heights are set, which are 0.5m, 1.0m and 1.5m respectively, and other parameters remain unchanged.

[0059] For each group of conditions, simulation calculation is performed according to the above steps to obtain corresponding results such as cement mortar loss rate and final loss amount.

[0060] Comparing the results of different conditions, for example, analyzing the change of loss rate and final loss amount when the tidal period changes from 12h to 24h, calculating the change rate and other quantitative indicators. According to these quantitative indicators, parameter sensitivity analysis atlas such as column chart and line chart is generated to intuitively show the influence degree of each parameter on the loss.

[0061] Through this comparative analysis, the influence of each tidal parameter on the cement mortar loss can be determined, which provides a basis for taking prevention and control measures in engineering according to different tidal conditions.

[0062] Please refer to Figure 2 The application further proposes a cement mortar loss simulation system based on tidal action, comprising: A parameter and model definition module is used to obtain the tidal simulation parameters and cement mortar structure parameters input by the user, and generate a three-dimensional geometric model and a finite element calculation grid based on the cement mortar structure parameters; A hydrodynamic boundary application module is used to establish boundary conditions for describing the change of hydrodynamic force based on the tidal simulation parameters, and apply the boundary conditions to the fluid domain boundary of the finite element calculation grid; A hydrodynamic analysis module is used to call a hydrodynamic calculation module to perform transient hydrodynamic analysis on the finite element calculation grid to obtain the water level, hydrodynamic pressure at each node and the pore water pressure distribution inside the soil body; a seepage coupling calculation module for establishing a cement mortar seepage calculation model, calculating hydraulic gradient and seepage velocity based on the obtained pore water pressure distribution, and determining and accumulating the mass seepage rate and total mass seepage amount of the cement mortar according to a preset seepage criterion function; a result generation and display module for processing the simulation results to generate visual output files and analysis reports representing the seepage range, degree and structural safety coefficient variation of the cement mortar.

[0063] Specifically, the parameter and model definition module can be realized by developing a dedicated user interface. After the user inputs the relevant parameters on the interface, the module calls the modeling program to automatically generate a three-dimensional geometric model according to the cement mortar structure parameters, and generates a finite element calculation grid using a meshing algorithm, such as an adaptive meshing algorithm that automatically adjusts the grid size according to the complexity of the structure.

[0064] The hydrodynamic boundary application module contains a boundary condition generation program that generates hydrodynamic boundary conditions according to the input tidal simulation parameters, preset formulas and algorithms, and then applies the boundary conditions to the corresponding fluid domain boundaries through the interface with the finite element calculation grid.

[0065] The hydrodynamic analysis module integrates a hydrodynamic calculation program that can receive finite element calculation grid and boundary condition data, perform calculations using a transient hydrodynamic analysis algorithm, and output water level, dynamic water pressure and pore water pressure distribution data for each node.

[0066] The seepage coupling calculation module contains a seepage calculation model and related algorithms. After receiving the pore water pressure distribution data, it calculates the hydraulic gradient and seepage velocity, determines the mass seepage rate according to the seepage criterion function, and accumulates the total mass seepage amount, while considering the impact of the seepage process on the structure parameters to achieve coupled calculation.

[0067] The result generation and display module receives the result data output by the seepage coupling calculation module, calls visualization processing programs and report generation programs, generates three-dimensional cloud maps, curve graphs and other visual output files, and analysis reports containing structural safety coefficient variations, and displays them to the user through the display interface.

[0068] In summary, the application obtains the tidal and cement mortar structure parameters, constructs a three-dimensional geometric model and a finite element grid, establishes a dynamic boundary condition combined with the characteristics of tides and waves, obtains the pore water pressure distribution by using transient hydrodynamic analysis, quantifies the mass loss by using a loss calculation model, and dynamically updates the material properties to simulate the coupling effect. At the same time, multiple groups of working conditions are supported to compare and analyze the influence of tidal parameters, and visual results and safety factor reports are generated. The method can accurately reproduce the cement mortar loss process under the action of tides, quantify the influence of different tidal parameters on the loss, provide a scientific basis for engineering loss prediction and prevention, and reduce safety risks.

[0069] It should be noted that, in this text, the relationship terms such as first and second are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between the entities or operations. Moreover, the term "include", "contain" or any other variant thereof is intended to cover non-exclusive inclusion, so that the process, method, article or equipment including a series of elements not only includes those elements, but also includes other elements not explicitly listed or inherent to such process, method, article or equipment.

[0070] Although the embodiments of the present application have been shown and described, it will be understood by those of ordinary skill in the art that various changes, modifications, substitutions and variations can be made to the embodiments without departing from the principles and spirit of the present application, and the scope of the present application is defined by the appended claims and their equivalents.

Claims

1. A method for simulating cement mortar loss based on tidal forces, characterized in that, include: S1. Obtain the tidal simulation parameters and cement mortar structure parameters input by the user. The tidal simulation parameters are used to define the simulated marine hydrodynamic environment, and the cement mortar structure parameters are used to define the physical properties of the engineering structure to be analyzed. S2. Based on the cement mortar structure parameters, generate a three-dimensional geometric model of the cement mortar structure and the surrounding soil in a preset coordinate system, and perform spatial discretization processing on the three-dimensional geometric model to obtain a finite element calculation mesh containing multiple elements and nodes. S3. Based on the tidal simulation parameters, establish boundary conditions to describe hydrodynamic changes, and apply the boundary conditions to the fluid domain boundary of the finite element calculation grid. S4. Call the hydrodynamic calculation module to perform transient hydrodynamic analysis on the finite element calculation grid, and solve for the distribution of water level, wave-induced hydrodynamic pressure and pore water pressure inside the soil at each node throughout the entire simulation period. S5. Establish a cement mortar loss calculation model. Based on the obtained pore water pressure distribution, calculate the hydraulic gradient and seepage velocity in each unit at each time step, and determine the cement mortar mass loss rate according to the preset loss criterion function. S6. Accumulate and calculate the total mass loss of each unit over the entire simulation period, process the simulation results, and generate a visual output file characterizing the range and extent of cement mortar loss.

2. The method for simulating cement mortar loss based on tidal action according to claim 1, characterized in that, The tidal simulation parameters include tidal type, tidal period, tidal amplitude, initial water level, wave type, wave height, wave period, and wave direction angle; the cement mortar structural parameters include the elastic modulus, Poisson's ratio, density, permeability coefficient, porosity, and contact properties of the cement mortar-soil interface.

3. The method for simulating cement mortar loss based on tidal action according to claim 1, characterized in that, The process of establishing boundary conditions to describe hydrodynamic changes and applying these boundary conditions to the fluid domain boundary of the finite element computation mesh includes: Based on the tidal type, tidal period, and tidal amplitude, a tidal function is generated to characterize the periodic water level changes, and this function is applied as a time-varying pressure boundary or water level boundary to the far-field water boundary of the finite element calculation grid. Based on the wave type, wave height, and wave period, either linear wave theory or nonlinear wave theory is selected to generate wave velocity field and pressure field functions, which are then superimposed on the tidal function and used together as input conditions for the far-field water boundary.

4. The method for simulating cement mortar loss based on tidal action according to claim 1, characterized in that, The process of calling the hydrodynamic calculation module to obtain the pore water pressure distribution inside the soil includes: On the finite element calculation mesh, a governing equation based on the seepage theory of unsaturated porous media is established; The hydrodynamic pressure exerted on the structural surface by the combined action of tides and waves is applied as a transient load to the solid-liquid coupling interface of the finite element calculation mesh. A time-stepping algorithm is used to iteratively solve the governing equations, obtaining the pore water pressure value of each node in the finite element computational grid at each time step, thus forming a spatiotemporal distribution dataset of pore water pressure.

5. The method for simulating cement mortar loss based on tidal action according to claim 1, characterized in that, The loss criterion function is a function of the infiltration rate and the critical initiation flow rate of the cement mortar. Determining the mass loss rate of the cement mortar includes: For any cell belonging to cement mortar in the finite element calculation mesh, if the calculated infiltration velocity inside the cell exceeds the preset critical initiation flow velocity threshold, the cell is determined to have lost flow. Based on the permeation velocity difference exceeding the critical initiation flow rate threshold, the mass loss rate components caused by chemical dissolution and physical scouring are calculated respectively using the preset dissolution model and scouring model. The total mass loss rate of the unit at the current time step is obtained by linearly superimposing the chemical dissolution rate component and the physical scouring rate component.

6. The method for simulating cement mortar loss based on tidal action according to claim 1, characterized in that, The simulation method further includes: After the calculation is completed at each time step, the material properties of the unit where the loss occurred are dynamically updated based on the calculated amount of cement mortar loss. The updated material properties include reducing the unit's density and elastic modulus, and increasing the unit's porosity and permeability coefficient. By using the updated material properties as the initial conditions for the calculation of the next time step, the coupled influence of the loss process on the structural physical properties and its hydraulic response can be simulated.

7. The method for simulating cement mortar loss based on tidal action according to claim 6, characterized in that, When the cumulative mass loss of any unit reaches the preset percentage of structural failure, the unit is removed from the set of effective load-bearing units and treated as a fluid unit in subsequent calculations to simulate the change in water flow path after partial structural collapse.

8. The method for simulating cement mortar loss based on tidal action according to claim 1, characterized in that, The generated visualization output file, characterizing the extent and degree of cement mortar loss, includes: A three-dimensional cloud map of the cumulative loss of cement mortar is generated, and the three-dimensional cloud map uses different colors to indicate the spatial distribution and numerical value of the loss. Generate a graph showing the change in cement mortar loss at key locations over time, which can be used to analyze the loss development process in a specific area. The output structure provides a comparative analysis report on its overall strength, stability, and safety factor before and after the simulation. The report includes a quantified safety factor reduction value.

9. The method for simulating cement mortar loss based on tidal action according to claim 1, characterized in that, The method further includes: Multiple sets of different tidal simulation parameters were set as comparative analysis conditions; For each set of comparative analysis conditions, repeat the steps described in claims 1 to 8 to obtain the corresponding simulation results; By comparing simulation results under different working conditions, the influence of changes in single parameters such as tidal period, tidal amplitude, and wave height on the cement mortar loss rate and final total loss is quantitatively analyzed, and parameter sensitivity analysis charts are generated.

10. A cement mortar loss simulation system based on tidal forces, characterized in that, include: The parameter and model definition module is used to obtain the tidal simulation parameters and cement mortar structure parameters input by the user, and generate a three-dimensional geometric model and finite element calculation mesh based on the cement mortar structure parameters. A hydrodynamic boundary application module is used to establish boundary conditions describing hydrodynamic changes based on the tidal simulation parameters, and to apply the boundary conditions to the fluid domain boundary of the finite element calculation mesh. The hydrodynamic analysis module is used to call the hydrodynamic calculation module to perform transient hydrodynamic analysis on the finite element calculation mesh, and solve for the water level, dynamic water pressure and pore water pressure distribution inside the soil at each node. The loss coupling calculation module is used to establish a cement mortar loss calculation model. Based on the obtained pore water pressure distribution, it calculates the hydraulic gradient and seepage velocity, and determines and accumulates the cement mortar mass loss rate and total mass loss according to the preset loss criterion function. The results generation and display module is used to process the simulation results and generate visual output files and analysis reports that characterize the range, degree and structural safety factor changes of cement mortar loss.

Images