Seismic Response of Earth-Rock Dams with Innovative Antiseepage Walls on the Effect of Microscopic Fluid-Solid Coupling

: In dykes and dam projects, the microscopic ﬂuid-solid coupling effect from the interaction between soil skeleton and pore water during an earthquake is crucial to consider, as it can lead to dam safety problems. To control seepage in medium-and small-sized dams, polymer antiseepage walls have emerged as effective measures. In recent years, this method has been increasingly utilized in projects worldwide as it is essential for preventing potential dam safety issues. To address concerns related to seismic safety, this study conducts theoretical analysis, model tests, and numerical simulations to investigate the seismic response of earth-rock dams with polymer antiseepage walls, with a speciﬁc focus on the microscopic ﬂuid-solid coupling effect. The dynamic viscoelastic constitutive model used in this study incorporates Biot’s theory of dynamic consolidation and the results of dynamic mechanical analysis (DMA) of polymer materials. To validate the model, a centrifuge test is performed, and it is then utilized to study the seismic response of earth-rock dams with polymer antiseepage walls. Furthermore, the inﬂuence of factors such as ﬂuid-solid coupling, water level, polymer material density, and wall thickness on the seismic response of dams with antiseepage walls is analyzed. Finally, the seismic safety of the earth-rock dam with the polymer antiseepage wall is thoroughly examined. The results emphasize the need to consider the ﬂuid-solid coupling effect, as factors like water level and design parameters of the antiseepage wall signiﬁcantly impact the seismic response of earth-rock dams with polymer antiseepage walls.


Introduction
In the context of modern engineering and infrastructure development, sustainability principles have become paramount to ensuring a harmonious coexistence between human progress and being environmentally friendly.Recently, there has been a rise in the use of innovative non-water-reacted polymer materials for constructing antiseepage walls in dam reinforcement and emergency rescue projects [1].This material is preferred due to its lightweight nature, non-polluting characteristics, excellent anti-seepage performance, and outstanding durability [2,3].Figure 1 presents a schematic diagram that illustrates the material synthesis process, specialized grouting equipment, and construction methods.
A special hole-drilling tool with static pressing equipment is used during the construction of antiseepage walls to penetrate the dam and create a continuous ultra-thin slot [4] Following this, a two-component polymer material is injected into the slot through a grouting pipe, initiating a chemical reaction.This reaction causes the combined materials to rapidly expand, effectively filling the slot and solidifying into polymer plates [5].Eventually, the adjacent hole plates bond together, forming a continuous ultra-thin polymer wall.This process significantly enhances the dam's antiseepage performance [6].The adoption of polymer material not only addresses dam safety concerns but also aligns with the broader goals of sustainable development.By mitigating seepage in medium and small-sized dams, these antiseepage walls contribute to enhanced infrastructure performance while minimizing environmental impact [7].
Sustainability 2023, 15, x FOR PEER REVIEW 2 of 28 to rapidly expand, effectively filling the slot and solidifying into polymer plates [5].Eventually, the adjacent hole plates bond together, forming a continuous ultra-thin polymer wall.This process significantly enhances the dam's antiseepage performance [6].The adoption of polymer material not only addresses dam safety concerns but also aligns with the broader goals of sustainable development.By mitigating seepage in medium and small-sized dams, these antiseepage walls contribute to enhanced infrastructure performance while minimizing environmental impact [7].Numerous dam failures have been triggered by earthquakes.For instance, the collapse of the Tangjiashan dam occurred after the Wenchuan earthquake in 2008 [8], the Fujinuma dam was damaged following the Tōhoku earthquake in 2011 [9], the Fundão mine tailings dam was broken after the southeastern Brazil earthquake in 2015 [10], and the landslides and dam damage were caused by the Jiuzhaigou earthquake in 2017 [11].During earthquakes, earth-rock dams experience dynamic fluid-solid coupling effects, which refer to the transfer of stress (pore pressure) between the pore water and soil skeleton within the dams [12].This coupling effect results in deformation and potential damage.Notably, under strong earthquakes, the impact of this coupling effect becomes more pronounced, causing various safety concerns for dams.These concerns include the formation of cracks, landslides, permanent deformation, liquefaction, and even the risk of dam break [13], as shown in Figure 2. Numerous dam failures have been triggered by earthquakes.For instance, the collapse of the Tangjiashan dam occurred after the Wenchuan earthquake in 2008 [8], the Fujinuma dam was damaged following the T ōhoku earthquake in 2011 [9], the Fundão mine tailings dam was broken after the southeastern Brazil earthquake in 2015 [10], and the landslides and dam damage were caused by the Jiuzhaigou earthquake in 2017 [11].During earthquakes, earth-rock dams experience dynamic fluid-solid coupling effects, which refer to the transfer of stress (pore pressure) between the pore water and soil skeleton within the dams [12].This coupling effect results in deformation and potential damage.Notably, under strong earthquakes, the impact of this coupling effect becomes more pronounced, causing various safety concerns for dams.These concerns include the formation of cracks, landslides, permanent deformation, liquefaction, and even the risk of dam break [13], as shown in Figure 2.
This study investigates the seismic response of earth-rock dams with polymer antiseepage walls.Two types of dynamic fluid-solid coupling effects [14] can be identified between earth-rock dams and reservoir water under seismic conditions.The first type is the macroscopic fluid-solid coupling effect that occurs at the interface between the reservoir water and the dam.At this interface, the water and dam are clearly separated on a macroscopic scale.The second type is the microscopic fluid-solid coupling effect that occurs within the earth-rock dams.Within the dam, the water and soil cannot be distinguished on a macroscopic scale.The focus of this study is the second type of coupling effect, which involves the microscopic interaction between the soil skeleton and pore water.This study investigates the seismic response of earth-rock dams with polymer antiseepage walls.Two types of dynamic fluid-solid coupling effects [14] can be identified between earth-rock dams and reservoir water under seismic conditions.The first type is the macroscopic fluid-solid coupling effect that occurs at the interface between the reservoir water and the dam.At this interface, the water and dam are clearly separated on a macroscopic scale.The second type is the microscopic fluid-solid coupling effect that occurs within the earth-rock dams.Within the dam, the water and soil cannot be distinguished on a macroscopic scale.The focus of this study is the second type of coupling effect, which involves the microscopic interaction between the soil skeleton and pore water.
For interpreting the microscopic fluid-solid coupling effect between the soil skeleton and pore water, Terzaghi [15] first proposed the concept of effective stress in 1925 and combined it with Darcy's law to create the soil consolidation theory.This theory was suitable for solving one-dimensional soil consolidation cases but not strictly applicable in twoand three-dimensional cases.Biot then developed a more accurate three-dimensional quasi-static consolidation equation and laid the groundwork for static fluid-solid coupling theory.This theory applied Terzaghi's principle and continuum theory to fluid and solid media and explained how soil skeleton and pore water are coupled in soil mass.In 1956, Biot [16] developed the dynamic consolidation theory based on the static consolidation equation, taking into account the inertia force of the soil skeleton and pore water.This theory, which is derived from Biot's theory, could be used to compute the coupling impact of the soil skeleton and pore water under dynamic conditions.Biot's theory is more effective at describing microscopic dynamic fluid-solid coupling than early theories.As a result, numerous studies have been conducted using Biot's dynamic consolidation hypothesis.Lin et al. [17] used Biot's theory to investigate fluid wave response in porous media.They used the sandstone test to determine the link between linear porosity and modulus and then to derive the parameters for Biot's dynamic consolidation equation.Ghaboussi et al. [18,19] established a finite element method by using the variational principle based on Biot's theory.Using finite element analysis, the dynamic response and liquefaction [20] of layered sand under seismic response were investigated.Zienkiewicz et al. [21][22][23] continued to research Biot's dynamic consolidation theory and simplified the equation in For interpreting the microscopic fluid-solid coupling effect between the soil skeleton and pore water, Terzaghi [15] first proposed the concept of effective stress in 1925 and combined it with Darcy's law to create the soil consolidation theory.This theory was suitable for solving one-dimensional soil consolidation cases but not strictly applicable in two-and three-dimensional cases.Biot then developed a more accurate three-dimensional quasi-static consolidation equation and laid the groundwork for static fluid-solid coupling theory.This theory applied Terzaghi's principle and continuum theory to fluid and solid media and explained how soil skeleton and pore water are coupled in soil mass.In 1956, Biot [16] developed the dynamic consolidation theory based on the static consolidation equation, taking into account the inertia force of the soil skeleton and pore water.This theory, which is derived from Biot's theory, could be used to compute the coupling impact of the soil skeleton and pore water under dynamic conditions.Biot's theory is more effective at describing microscopic dynamic fluid-solid coupling than early theories.As a result, numerous studies have been conducted using Biot's dynamic consolidation hypothesis.Lin et al. [17] used Biot's theory to investigate fluid wave response in porous media.They used the sandstone test to determine the link between linear porosity and modulus and then to derive the parameters for Biot's dynamic consolidation equation.Ghaboussi et al. [18,19] established a finite element method by using the variational principle based on Biot's theory.Using finite element analysis, the dynamic response and liquefaction [20] of layered sand under seismic response were investigated.Zienkiewicz et al. [21][22][23] continued to research Biot's dynamic consolidation theory and simplified the equation in accordance with various unknowns.They rewrote it into four different forms: "u-p" form, "u-w" form, "u-U" form, and "u-w-p" form.Pan [24], Liu [25], and Chen [26] explored dam-foundationreservoir interactions, highlighting the impact of foundation size and dam foundation dynamics on dam acceleration.Kong [27] proposed a spectral element-based numerical method for fluid-solid seismic wave interactions.Zhang [28] bridged wavefield simulation and dam seismic analysis.Rayegani [29,30] stressed the need for advanced strategies like intelligent dampers to enhance dam seismic resistance.Jing [31] compared near-field and far-field seismic wave effects on earth-rock dams, enriching our understanding of complex interactions.These theories and methods provide new approaches for studying the dynamic coupling effect between pore water and earth-rock dams under seismic action and contribute to investigating the dynamic coupling effect of pore water and earth-rock dams under earthquake action.
Based on the aforementioned review, the polymer antiseepage wall has been widely used in the reinforcement projects of earth-rock dams as a new type of antiseepage wall.However, there are limitations in the current understanding of its dynamic interaction with the soil skeleton and pore water.Although studies have been conducted on the static fluid-solid coupling of polymer antiseepage walls in earth-rock dams, the dynamic fluid-solid coupling model based on Biot's static consolidation theory is insufficient for fully describing this dynamic interaction.Therefore, it is of high significance to adopt a coupling theory that accurately represents the fine-scale fluid-solid coupling effect in the seismic response analysis of earth-rock dams to establish a fully coupled dynamic model that incorporates the fluid-solid interaction.In this study, based on Biot's dynamic consolidation theory, a static and dynamic finite differential model of the earth-rock dam with a polymer antiseepage wall was developed to examine the effect of microscopic fluid-solid coupling.To reflect the material properties in the model, a dynamic viscoelastic constitutive model of polymer materials was developed using DMA test results.The numerical model was then validated by centrifuge test data [32] and the seismic response laws of antiseepage walls [33] were studied.Furthermore, the response of an earth-rock dam with a concrete antiseepage wall was also investigated for comparison.Finally, the seismic stability of the earth-rock dam with a polymer antiseepage wall was analyzed, and the related mechanisms were discussed.

Theoretical Research
Biot's dynamic consolidation equation is utilized to describe the microscopic dynamic fluid-solid coupling effect between the soil skeleton and pore space.Biot's dynamic consolidation equation serves as the coupling equation, and its solution is obtained using the finite difference method [34].The following assumptions are proposed to make the derivation of Biot's dynamic consolidation theory [35][36][37]: Soil skeleton is an ideal continuous porous medium, of which compressibility is ignored; 2.
Pore water is compressible and its flow obeys the generalized Darcy's law; 3.
Soils are isotropic medium and pores are interconnected; 4.
Pore size is much smaller than the seimic wavelength; 5.
The effect of temperature is ignored.
Considering the compressibility of soil particles, the revised effective stress is as follow [38]: In the formula, σ ij is the total stress, δ ij is the Kronecker function, p is the pore water pressure, α is the Biot's constant, α = 1 − Kb/Ks; Kb and Ks are soil skeleton and soil particles of the bulk modulus, and the value of α is 1, assuming the soil particles being incompressible.
The Mohr-Coulomb model is chosen as the constitutive model for the analysis.To ensure the accuracy of the software analysis, the fatigue behavior of the soil was not considered in this model.The total deformation consists of two parts: elastic deformation and plastic deformation.
The elastic deformation is represented by Hooke's law: In the formula, λ and µ are Lame constants, λ = Ev (1+v)(1−2v) , µ = E 2(1+v) , ε ij is the strain rate of soil skeleton, and ε kk is the volumetric strain, The plastic deformation is represented by the yield criterion of the Mohr-Coulomb model: In the formula, c, ϕ are the cohesion and angle of internal friction of the material.The motion of the fluid is expressed by Darcy's law: In the formula, k ij is the dynamic permeability coefficient, k ij = k/ρg, k is the permeability coefficient in soil mechanics, and k(s) is the coefficient related to fluid saturation.
The equilibrium equations include momentum balance equation and fluid mass balance equation.Momentum balance equation is listed as follows: In the formula, ..
u i is the acceleration of soil skeleton, ρ is the total effective density, ρ = nρ f + (1 − n)ρ s , ρ f is the density of water, ρ s is the density of soil skeleton, and g i is the acceleration of gravity.
The fluid mass balance equation is listed as follows: In the formula, q i,j is the velocity of fluid seepage, q i,j = ∂q i ∂x i = ∂q i ∂x + ∂q i ∂y , q v is the inflow of liquid per unit time, and ζ is the change of the fluid volume in the porous media per unit volume.
The constitutive equation of pore fluid is listed as follows: In the formula, M is Biot's modulus and ε is the volumetric strain of soil mass.Substituting the fluid mass balance equation into the pore fluid's continuity equation, the expression for the fluid continuity equation is obtained: Substituting Darcy's law into the fluid continuity equation and assuming the soil saturation being 1:

Static Fluid-Solid Coupling Model
In this paper, a static fluid-solid coupling model for earth-rock dam with polymer antiseepage wall was established, taking the earth-rock dam of Jiulong Reservoir as the prototype.

Model Building
The earth-rock dam of Jiulong Reservoir is a small dam situated in Jiulong River, Luoshan County, Xinyang City, Henan Province, China.The materials of the dam's body and foundation are mainly composed of heavy silty soil with an average permeability coefficient of 2.6 × 10 −5 cm/s.The dam is 15 m tall, 240 m long, 5.6 m wide at its top, and 90 m wide at its bottom.The earth-rock dam needed to be strengthened as there were significant leak problems after operating for 60 years.So a polymer antiseepage wall with a thickness of 0.03 m was constructed in the dam for reinforcement.The wall is 2 m deep into the dam's base, with a total length of 215 m along the dam's axis, as shown in Figure 3a.The FLAC3D 7.0 software was used for the analysis.The dimensions of the earth-rock dam model were the same as the prototype, with a width of 189.6 m and a dam base height of 15 m.The dam was evenly distributed along the dam axis direction, and the upstream water level was 8.6 m.To improve calculation efficiency, the model was set to unit thickness.
An antiseepage wall, located 0.5 m away from the dam axis line on the upstream dam slope, was added inside the earth-rock dam.The wall had a height of 17 m, extending 2 m into the dam base, and a thickness of 0.03 m.The cross-section of the earth-rock dam with a polymer antiseepage wall is depicted in Figure 3b.The grid in the antiseepage wall area was densified, while the partitioned grid met the requirements for conducting seismic analysis.For the continuity of hydrostatic calculation, the numerical model was divided into 3505 units and 7400 nodes.Hexahedral computational elements were used for calculations, as seen in Figure 3c.
significant leak problems after operating for 60 years.So a polymer antiseepage wall with a thickness of 0.03 m was constructed in the dam for reinforcement.The wall is 2 m deep into the dam's base, with a total length of 215 m along the dam's axis, as shown in Figure 3a.
The FLAC3D 7.0 software was used for the analysis.The dimensions of the earthrock dam model were the same as the prototype, with a width of 189.6 m and a dam base height of 15 m.The dam was evenly distributed along the dam axis direction, and the upstream water level was 8.6 m.To improve calculation efficiency, the model was set to unit thickness.An antiseepage wall, located 0.5 m away from the dam axis line on the upstream dam slope, was added inside the earth-rock dam.The wall had a height of 17 m, extending 2 m into the dam base, and a thickness of 0.03 m.The cross-section of the earth-rock dam with a polymer antiseepage wall is depicted in Figure 3b.The grid in the antiseepage wall area was densified, while the partitioned grid met the requirements for conducting seismic analysis.For the continuity of hydrostatic calculation, the numerical model was divided into 3505 units and 7400 nodes.Hexahedral computational elements were used for calculations, as seen in Figure 3c.

Boundary Conditions
In fluid-soil coupling analysis, the boundary conditions consist of a stress boundary, a displacement boundary, a flow potential boundary, and a flow boundary.The fixed boundaries are applied to the bottom of the dam foundation to restrict movement in three directions.The horizontal displacement and displacement along the dam axis direction are constrained by the left and right sides of the dam foundation, respectively.The upstream reservoir water exerts stress and seepage in the stress field and seepage field of the earth-rock dam.On the upstream slope and the foundation surface, it applies hydrostatic pressure and pore water pressure.Additionally, a drainage boundary with a value of 0 is set as the flow boundary on the right side of the dam foundation.In addition to considering stress and displacement boundaries, the movement of pore water relative to the soil skeleton, i.e., pore water pressure boundary and fluid flow boundary conditions, also needs to be considered.The boundary conditions of the fluid-solid coupling earth-rock dam model with a polymer antiseepage wall are as follows: Displacement boundary: Stress boundary:

Boundary Conditions
In fluid-soil coupling analysis, the boundary conditions consist of a stress boundary, a displacement boundary, a flow potential boundary, and a flow boundary.The fixed boundaries are applied to the bottom of the dam foundation to restrict movement in three directions.The horizontal displacement and displacement along the dam axis direction are constrained by the left and right sides of the dam foundation, respectively.The upstream reservoir water exerts stress and seepage in the stress field and seepage field of the earthrock dam.On the upstream slope and the foundation surface, it applies hydrostatic pressure and pore water pressure.Additionally, a drainage boundary with a value of 0 is set as the flow boundary on the right side of the dam foundation.In addition to considering stress and displacement boundaries, the movement of pore water relative to the soil skeleton, i.e., pore water pressure boundary and fluid flow boundary conditions, also needs to be considered.The boundary conditions of the fluid-solid coupling earth-rock dam model with a polymer antiseepage wall are as follows: Displacement boundary: Stress boundary: Flow potential boundary: Flow boundary: In the numerical analysis, the Mohr-Coulomb elastoplastic constitutive model is selected for soils, the linear elastic constitutive model is for concrete, and the generalized Maxwell constitutive model is for viscoelastic polymer materials.
For comparison, the dam model with a concrete antiseepage wall was established with a thickness of 0.5 m.Other geometric parameters are the same as those of the polymer antiseepage wall.The parameters of materials are shown in Table 1 [39].Darcy's law generally meets the seepage condition of a dam while an initial hydraulic gradient exists.The boundary condition of the known water head is represented by the contact surface (Γ) between the dam body and the polymer wall, listed as follows: Under seismic reaction, the contact surface (Γ) between the dam body and the polymer wall should also meet the following boundary condition.
In the formula, p is the pore water pressure, n is the normal direction, and ..
u n is the normal acceleration.
In the fluid-solid coupling calculation, the compressibility of soil particles should be ignored, and the default Biot's coefficient is set to 1. Fluid bulk modulus K f is the compressibility of fluid.The relationship of K f , fluid pressure increment ∆P and fluid volume strain ∆V f /V f is as follows: To study the dynamic viscoelastic constitutive model of polymer materials, the DMA tests [40] were carried out.In these tests, parameters of the material's dynamic viscoelastic-ity, storage modulus, and loss factor were tested.Then the dynamic viscoelastic constitutive model of polymer materials was established according to the relationship between parameters and frequency.

Derivation of Model
The generalized Maxwell model, also known as the spring-dashpot combination model, consists of multiple Maxwell units connected in parallel.Each Maxwell unit includes an elastic strain spring with an elastic modulus of E i and a viscous strain dashpot with a corresponding viscosity of η i as shown in Figure 4.This model accurately describes the stress relaxation behavior and viscoelastic properties of viscoelastic materials under dynamic loads.It is composed of multiple elements, corresponding to the complex relaxation times of high-polymer grouting materials due to their complex composition.The parameters of the model are relatively easy to obtain and can be acquired through fitting DMA test data to calculate the viscoelastic constitutive model parameters.The parameters obtained from the generalized Maxwell model can be directly applied in finite element software for efficient dynamic response analysis of materials.Thus, this model finds wide applications.Consequently, in the preliminary research and subsequent finite element analysis, the generalized Maxwell model is adopted as the dynamic viscoelastic constitutive model of polymer materials based on DMA test results.The generalized Maxwell model provides an intuitive reflection of the stress-relaxation behavior of materials.The storage modulus  and loss modulus  ″ of the generalized Maxwell model are expressed as follows: =    1 +   (18) In the formula,  is the quasi-static modulus,  is the elastic modulus of each Maxwell model, and  is the relaxation time.The relationship between relaxation time  and viscosity  is  =  / .
The expression of the stress relaxation modulus and the physical explanation for the viscoelastic behavior are provided by the Maxwell model.The number of Maxwell units and the relaxation time required vary for different viscoelastic materials.For materials with simple viscoelastic properties, a two-or three-unit Maxwell model can be used.However, more Maxwell units are required for complex models.Additionally, there is uncertainty in selecting the parameters due to the principle of equivalent combinations of different parameter sets.The high computational efficiency of the fourth-order generalized Maxwell model has been verified [41,42] In the formula, E e is the quasi-static modulus, E i is the elastic modulus of each Maxwell model, and τ i is the relaxation time.The relationship between relaxation time τ i and viscosity η i is The expression of the stress relaxation modulus and the physical explanation for the viscoelastic behavior are provided by the Maxwell model.The number of Maxwell units and the relaxation time required vary for different viscoelastic materials.For materials with simple viscoelastic properties, a two-or three-unit Maxwell model can be used.However, more Maxwell units are required for complex models.Additionally, there is uncertainty in selecting the parameters due to the principle of equivalent combinations of different parameter sets.The high computational efficiency of the fourth-order generalized Maxwell model has been verified [41,42] so it was chosen for analysis.The energy storage modulus E and loss modulus E in the fourth-order generalized Maxwell model are as follows: In analysis, the lsqcurvefit function program is written by Matlab R2019b software for fitting procedures.The relaxation time is bounded as follows for calculation efficiency: The elastic modulus of Maxwell model is also constrained as follows: 0 Three polymer densities were calculated in the numerical analysis.Take 260 kg/m 3 , for example.The fourth-orderfourth-order parameters of the dynamic viscoelastic constitutive model are shown in Table 2.In this paper, a fourth-order generalized Maxwell model is derived to determine the stress increment at the current time step.The finite difference scheme of the Maxwell model is: In the formula, .
u 1 is the velocity, F 1 is the force, k 1 is the elastic modulus, and η 1 is the viscosity coefficient.Simplifying the formula as follows: In the formula, F n 1 is the force of the previous step and F n+1 1 is the force of the current step.
Equating the above equation to the stress-strain relationship, the expression is listed as follows: In the formula, σ n+1 1d_ij is the bias stress of the current step, σ n 1d_ij is the bias stress of the previous step, G 1 is the shear modulus and ∆ε d ij is the bias strain increment.The spherical stress can be expressed as follows: In the formula, K is the bulk modulus.In the generalized Maxwell model, multiple Maxwell models are in parallel relationships with each other, and a single spring could be viewed as a simplified version of a Maxwell body.Hence, the expression of the total deflection stress of the fourth-order Maxwell is listed as follows: Then the total spherical stress expression of the fourth-order Maxwell model is: (

Dynamic Fluid-Solid Coupling Model
The dynamic fluid-solid coupling model of an earth-rock dam with a polymer antiseepage wall was established by applying the free-field boundary, adding mechanical

Dynamic Fluid-Solid Coupling Model
The dynamic fluid-solid coupling model of an earth-rock dam with a polymer antiseepage wall was established by applying the free-field boundary, adding mechanical damping, and using the dynamic viscoelastic constitutive model of polymer materials.

Boundary Conditions
In dynamic calculation, the free-field boundary condition is set to eliminate the wavereflecting effect, as shown in Figure 6, which is from the dynamic analysis manual of FLAC (Itasca).If body grids and free-field grids move in the same manner, the normal dampers are not activated; if not, the dampers take a static boundary approach to absorb the seismic energy.The body grids were assumed to be uniform.

Earthquake Wave
The El−Centro wave, the first seismic wave captured by humans, has been widely used in model tests and numerical analyses.To accurately recreate the original earthquake waveform, the El Centro wave with an acceleration amplitude of 0.4 g was selected.The acceleration time curve is shown in Figure 7.

Earthquake Wave
The El−Centro wave, the first seismic wave captured by humans, has been widely used in model tests and numerical analyses.To accurately recreate the original earthquake waveform, the El Centro wave with an acceleration amplitude of 0.4 g was selected.The acceleration time curve is shown in Figure 7.
To study the impact of different seismic wave acceleration amplitudes on the model, seismic waves with acceleration amplitudes of 0.05, 0.1, 0.2, and 0.4 g were simulated in this experiment.These seismic waves correspond to seismic intensities of VI, VII, VIII, and IX degrees, respectively.In this experiment, El-Centro waves with different peak values will be applied to the bottom of the model, ranging from smallest to largest and lasting for 20 s, as working conditions 1 to 4.

Earthquake Wave
The El−Centro wave, the first seismic wave captured by humans, has been widely used in model tests and numerical analyses.To accurately recreate the original earthquake waveform, the El Centro wave with an acceleration amplitude of 0.4 g was selected.The acceleration time curve is shown in Figure 7.To study the impact of different seismic wave acceleration amplitudes on the model, seismic waves with acceleration amplitudes of 0.05, 0.1, 0.2, and 0.4 g were simulated in this experiment.These seismic waves correspond to seismic intensities of VI, VII, VIII, and

Damping
Internal material friction and potential slippage of the contact surface cause damping.Rayleigh damping has been often used in the structures and dynamic analysis of elastomers to reduce the self-oscillation amplitude of the mode.
Assuming that the damping matrix C is related to the stiffness matrix K and the mass matrix M: In the formula, the damping constants α and β are proportional to mass and stiffness, respectively, determined by the following equation: In the formula, ζ min is the critical damping ratio, and generally set to 0.05 [43] in the calculation of geotechnical materials, and ω min is the minimum center frequency.

Model Verification
The centrifuge vibration table increases the gravity on the model by high-speed rotation, enabling the model to reach the stress state of the prototype and replicate the structural characteristics in the model, facilitating the observation of possible deformations and damages to the structure.To study the seismic response of earth-rock dams with polymer antiseepage walls and concrete antiseepage walls, centrifuge tests were carried out in the National Key Laboratory of Soft Soils and Geoenvironmental Engineering (Zhejiang University).The ZJU400 geotechnical centrifuge and vibration table for tests are shown in Figure 8.
rotation, enabling the model to reach the stress state of the prototype and replicate the structural characteristics in the model, facilitating the observation of possible deformations and damages to the structure.To study the seismic response of earth-rock dams with polymer antiseepage walls and concrete antiseepage walls, centrifuge tests were carried out in the National Key Laboratory of Soft Soils and Geoenvironmental Engineering (Zhejiang University).The ZJU400 geotechnical centrifuge and vibration table for tests are shown in Figure 8.The similar relationship between the non-uniform stress centrifuge model and the prototype can be determined based on the similarity of the stress-strain relationship, as depicted in Table 3.The geometric similarity ratio of the dynamic centrifuge test is 1/100, while the acceleration similarity ratio for the centrifuge rotation is 50 g, based on the size of the prototype and the centrifuge test model box.In terms of the dam design parameters, the dam has a height of 15 cm, a foundation thickness of 15 cm, a top width of 5 cm, and a bottom width of 73 cm.The ratio of the upstream and downstream dam slopes is 1:2.3.The similar relationship between the non-uniform stress centrifuge model and the prototype can be determined based on the similarity of the stress-strain relationship, as depicted in Table 3.The geometric similarity ratio of the dynamic centrifuge test is 1/100, while the acceleration similarity ratio for the centrifuge rotation is 50 g, based on the size of the prototype and the centrifuge test model box.In terms of the dam design parameters, the dam has a height of 15 cm, a foundation thickness of 15 cm, a top width of 5 cm, and a bottom width of 73 cm.The ratio of the upstream and downstream dam slopes is 1:2.3.Moreover, the antiseepage wall has a height of 17 cm and extends 2 cm into the dam foundation.
The dynamic fluid-solid coupling calculation was performed for working conditions 1-4.The model extracted the peak acceleration at the monitoring points corresponding to the experimental measurements to determine the simulated and experimental peak acceleration values of the impermeable wall in the earth-rock dam.The difference between these values represented the relative error, which was visualized in Figure 9. Additionally, the average relative error and variance are presented in Table 4.
From the analysis, the average relative error between the calculated results and the centrifuge test results is 10.4% and 9.9%, with the variance being 7.4% and 6.3%, respectively.The calculated peak acceleration values at each measurement point were close to the experimental values.The finite element numerical model of the earth-rock dam with a polymer antiseepage wall effectively reflects the acceleration response of the earth-rock dam under seismic action.Therefore, it is considered that the finite element numerical model is reasonable.So the dynamic simulation model was used for the following analysis.
The dynamic fluid-solid coupling calculation was performed for working conditions 1-4.The model extracted the peak acceleration at the monitoring points corresponding to the experimental measurements to determine the simulated and experimental peak acceleration values of the impermeable wall in the earth-rock dam.The difference between these values represented the relative error, which was visualized in Figure 9. Additionally, the average relative error and variance are presented in Table 4.  Table 4. Acceleration response error under different conditions.

Coupling Effect
Seismic response analysis is conducted under working condition 4 to study the effects of fluid-solid coupling of an earth-rock dam with a polymer antiseepage wall.According to Figure 10, the peak acceleration of the antiseepage wall increases progressively along its height, with the maximum acceleration observed at the top.Under the influence of the coupling effect, the peak acceleration of the wall exceeds that under non-coupling conditions.This is because, after considering the fluid-structure interaction, the upstream reservoir water has simultaneous effects on both the stress field and the seepage field of the earth-rock dam, resulting in an increased acceleration response.During an earthquake, the earth-rock dam constrains the movement of the polymer antiseepage wall.Therefore, the peak acceleration of the wall is significantly higher under the coupling effect.The compressive stress of the antiseepage wall initially decreases along the wall height and then increases, with the maximum compressive stress located at the bottom of the wall.Under coupling and non-coupling conditions, the maximum values are 202 kPa and 162 kPa, respectively, and the magnitude of the amplification is 25%.Similarly, the tensile stress initially increases along the wall height and then decreases, with the maximum tensile stress located at the bottom of the wall.Under coupling and non-coupling conditions, the maximum tensile stress values are 177 kPa and 156 kPa, respectively, showing a 13% increase in the maximum tensile stress when considering coupling effects.and 162 kPa, respectively, and the magnitude of the amplification is 25%.Similarly, the tensile stress initially increases along the wall height and then decreases, with the maximum tensile stress located at the bottom of the wall.Under coupling and non-coupling conditions, the maximum tensile stress values are 177 kPa and 156 kPa, respectively, showing a 13% increase in the maximum tensile stress when considering coupling effects.The displacement distribution pattern of the polymer antiseepage wall, as depicted in Figure 11, demonstrates the following characteristics: The displacement decreases with increasing wall height, with the maximum displacement occurring at the top of the wall.Specifically, under coupling and non-coupling conditions, the maximum displacements of the antiseepage wall are 9.94 cm and 2.99 cm, respectively.The coupling effect significantly affects the displacement of the antiseepage wall.This can be attributed to the coupling effect, which causes an increased displacement of the earth-rock dam due to the influence of the seepage field.Due to the constraint imposed by the earth-rock dam, the displacement of the polymer antiseepage wall also increases.
In conclusion, coupling effects should not be ignored in the safety evaluation and analysis of earth-rock dams [44].After considering these effects, it is found that the acceleration, stress, and displacement of the polymer antiseepage wall are relatively larger than The displacement distribution pattern of the polymer antiseepage wall, as depicted in Figure 11, demonstrates the following characteristics: The displacement decreases with increasing wall height, with the maximum displacement occurring at the top of the wall.Specifically, under coupling and non-coupling conditions, the maximum displacements of the antiseepage wall are 9.94 cm and 2.99 cm, respectively.The coupling effect significantly affects the displacement of the antiseepage wall.This can be attributed to the coupling effect, which causes an increased displacement of the earth-rock dam due to the influence of the seepage field.Due to the constraint imposed by the earth-rock dam, the displacement of the polymer antiseepage wall also increases.

Reservoir Water Level
The variation in reservoir water level not only affects the stress fields and seepage fields of the earth-rock dam but also influences the acceleration, stress, and deformation behavior of the dam under seismic response.In this study, the effect of water level on earth-rock dams with polymer antiseepage walls was investigated at four different levels: 1 m, 5 m, 9 m, and 13 m, aiming to simulate the different conditions during low and high water periods.
The peak acceleration of the polymer antiseepage wall at different reservoir water levels is shown in Figure 12a.It shows that when the water level is the same, the acceleration increases with wall height, and the greatest acceleration occurs at the top of the wall.For example, the acceleration at the top of the wall is 0.70 g, 0.72 g, 0.74 g, and 0.75 g, respectively, with the water level being 1 m, 5 m, 9 m, and 13 m.From Figure 12b,c, the compressive and tensile stresses first decrease then increase from the bottom to the top of the wall.The greatest compressive and tensile stresses occur at the bottom of the wall.The great compressive stress is 165 kPa, 175 kPa, 201 kPa, and 231 kPa, respectively, and the great tensile stress is 145 kPa, 159 kPa, 180 kPa, and 206 kPa, respectively, at the water levels of 1 m, 5 m, 9 m, and 13 m.From Figure 12, it can be observed that the peak acceleration, compressive, and tensile stresses of the dam increase with the water level.This is due to the fact that the dam experiences greater pressure when the water level is higher.As the water level of the dam increases, it experiences less effective stress in horizontal and vertical directions, thus leading to its loss of stability.
(a) In conclusion, coupling effects should not be ignored in the safety evaluation and analysis of earth-rock dams [44].After considering these effects, it is found that the acceleration, stress, and displacement of the polymer antiseepage wall are relatively larger than before.Therefore, in the subsequent calculations, the fluid-solid coupling effects between the soil skeleton and pore water will be taken into account.

Reservoir Water Level
The variation in reservoir water level not only affects the stress fields and seepage fields of the earth-rock dam but also influences the acceleration, stress, and deformation behavior of the dam under seismic response.In this study, the effect of water level on earth-rock dams with polymer antiseepage walls was investigated at four different levels: 1 m, 5 m, 9 m, and 13 m, aiming to simulate the different conditions during low and high water periods.
The peak acceleration of the polymer antiseepage wall at different reservoir water levels is shown in Figure 12a.It shows that when the water level is the same, the acceleration increases with wall height, and the greatest acceleration occurs at the top of the wall.For example, the acceleration at the top of the wall is 0.70 g, 0.72 g, 0.74 g, and 0.75 g, respectively, with the water level being 1 m, 5 m, 9 m, and 13 m.From Figure 12b,c, the compressive and tensile stresses first decrease then increase from the bottom to the top of the wall.The greatest compressive and tensile stresses occur at the bottom of the wall.The great compressive stress is 165 kPa, 175 kPa, 201 kPa, and 231 kPa, respectively, and the great tensile stress is 145 kPa, 159 kPa, 180 kPa, and 206 kPa, respectively, at the water levels of 1 m, 5 m, 9 m, and 13 m.From Figure 12, it can be observed that the peak acceleration, compressive, and tensile stresses of the dam increase with the water level.This is due to the fact that the dam experiences greater pressure when the water level is higher.As the water level of the dam increases, it experiences less effective stress in horizontal and vertical directions, thus leading to its loss of stability.
Figure 13 shows that the displacement increases from the bottom to the top of the wall, and the greatest displacement occurs at the top of the wall, with values of 3.68 cm, 7.84 cm, and 10.06 cm, respectively, at water levels of 1 m, 5 m, 9 m, and 13 m.However, when the water level is 13 m, the displacement first increases then decreases from the bottom to the top of the wall, with a maximum displacement of 12.57 cm at the middle of the wall.According to the lateral comparison, the rising water level may result in the horizontal and vertical permanent displacement of the dam and the displacement of the antiseepage wall.levels of 1 m, 5 m, 9 m, and 13 m.From Figure 12, it can be observed that the peak acceleration, compressive, and tensile stresses of the dam increase with the water level.This is due to the fact that the dam experiences greater pressure when the water level is higher.As the water level of the dam increases, it experiences less effective stress in horizontal and vertical directions, thus leading to its loss of stability.Figure 13 shows that the displacement increases from the bottom to the top of the wall, and the greatest displacement occurs at the top of the wall, with values of 3.68 cm, 7.84 cm, and 10.06 cm, respectively, at water levels of 1 m, 5 m, 9 m, and 13 m.However, when the water level is 13 m, the displacement first increases then decreases from the bottom to the top of the wall, with a maximum displacement of 12.57 cm at the middle of the wall.According to the lateral comparison, the rising water level may result in the horizontal and vertical permanent displacement of the dam and the displacement of the antiseepage wall.7.84 cm, and 10.06 cm, respectively, at water levels of 1 m, 5 m, 9 m, and 13 m.However, when the water level is 13 m, the displacement first increases then decreases from the bottom to the top of the wall, with a maximum displacement of 12.57 cm at the middle of the wall.According to the lateral comparison, the rising water level may result in the horizontal and vertical permanent displacement of the dam and the displacement of the antiseepage wall.

Wall's Density
The mechanical properties of polymer materials, including their density, play a significant role [22].This section analyzed the seismic response of the earth-rock dam with a polymer antiseepage wall at three different densities: 160 kg/m 3 , 260 kg/m 3 , and 350 kg/m 3 .This analysis takes into account the impact of density on construction costs.The input peak seismic acceleration was 0.4 g, and the water level was 8.6 m.
The peak acceleration of the polymer antiseepage wall increases slightly with density, as shown in Figure 14a.From Figure 14b, the compressive stress first decreases then increases from the bottom to the top of the wall, and the greatest compressive stress occurs at the bottom of the wall with values of 179 kPa, 199 kPa, and 265 kPa, respectively.From Figure 14c, the tensile stress has the same law as the compressive stress.The displacement of the polymer antiseepage wall is shown in Figure 15.The distribution increases from the bottom to the top of the wall, and the greatest displacement occurs at the top of the wall with values of 10.04 cm, 9.94 cm, and 9.71 cm, respectively. (a)

Wall's Density
The mechanical properties of polymer materials, including their density, play a significant role [22].This section analyzed the seismic response of the earth-rock dam with a polymer antiseepage wall at three different densities: 160 kg/m 3 , 260 kg/m 3 , and 350 kg/m 3 .This analysis takes into account the impact of density on construction costs.The input peak seismic acceleration was 0.4 g, and the water level was 8.6 m.
The peak acceleration of the polymer antiseepage wall increases slightly with density, as shown in Figure 14a.From Figure 14b, the compressive stress first decreases then increases from the bottom to the top of the wall, and the greatest compressive stress occurs at the bottom of the wall with values of 179 kPa, 199 kPa, and 265 kPa, respectively.From Figure 14c, the tensile stress has the same law as the compressive stress.The displacement of the polymer antiseepage wall is shown in Figure 15.The distribution increases from the bottom to the top of the wall, and the greatest displacement occurs at the top of the wall with values of 10.04 cm, 9.94 cm, and 9.71 cm, respectively.
This section shows that the distribution of the acceleration, stress, and displacement of the antiseepage wall at different densities is very similar.The acceleration increases with the density, while the difference is small as the wall shares little energy during seismic response due to the polymer antiseepage wall's small volume and mass compared with the dam.The compressive and tensile stresses of the wall increase with density; this may be due to the elastic modulus of the polymer material increasing, as structures with higher stiffness experience a larger seismic response than structures with lower stiffness.
at the bottom of the wall with values of 179 kPa, 199 kPa, and 265 kPa, respectively.From Figure 14c, the tensile stress has the same law as the compressive stress.The displacement of the polymer antiseepage wall is shown in Figure 15.The distribution increases from the bottom to the top of the wall, and the greatest displacement occurs at the top of the wall with values of 10.04 cm, 9.94 cm, and 9.71 cm, respectively.This section shows that the distribution of the acceleration, stress, and displacement of the antiseepage wall at different densities is very similar.The acceleration increases with the density, while the difference is small as the wall shares little energy during seismic response due to the polymer antiseepage wall's small volume and mass compared with the dam.The compressive and tensile stresses of the wall increase with density; this

Wall's Thickness
To satisfy the needs of seepage control and reinforcement in dam projects, various thicknesses of antiseepage walls can be formed using polymer grouting based on the width of the slot.In this study, wall thicknesses of 2 cm, 3 cm, 4 cm, and 5 cm were used to analyze the seismic response of earth-rock dams with polymer antiseepage walls under the fluid-solid coupling effect.
From Figure 16a, the acceleration increases from the bottom to the top of the wall, and the greatest acceleration occurs at the top of the wall.From Figure 16b,c, the compressive stress first decreases then increases with the increase in wall thickness.The greatest compressive stress occurs at the bottom of the wall; for example, the greatest compressive stress of the wall is 222 kPa, 199 kPa, 190 kPa, and 177 kPa, respectively, and the greatest tensile stress is 195 kPa, 177 kPa, 168 kPa, and 159 kPa, respectively.Figure 17 shows the displacement of polymer antiseepage walls with different wall thicknesses.The distribution of the displacement increases from the bottom to the top of the wall, and the greatest distribution occurs at the top of the wall.When the thickness is 2 cm, 4 cm, or 5 cm, the greatest displacement of the wall is 10.56 cm, 9.92 cm, and 9.87 cm, respectively.
Sustainability 2023, 15, x FOR PEER REVIEW 21 of 28 may be due to the elastic modulus of the polymer material increasing, as structures with higher stiffness experience a larger seismic response than structures with lower stiffness.

Wall's Materials
In order to compare the differences between the dam with the polymer wall and dam with the traditional concrete wall, dynamic fluid-solid coupling calculations w performed on earth-rock dams with polymer and concrete antiseepage walls under w ing conditions 1-4.The results are shown in Figure 18.The development tendency o celeration increases with the wall height and the input seismic acceleration, and the gr est acceleration occurs at the top of the wall.The peak accelerations of polymer and crete antiseepage walls are shown in Table 5.The trends of peak acceleration for the types of antiseepage walls are generally in agreement.A polymer antiseepage wall's p acceleration is lower than that of a concrete one under the same working conditions; is because polymer materials have a lower stiffness and lighter mass than concrete m These figures show that as the wall's thickness increases, the acceleration, compressive, and tensile stresses all decrease.This phenomenon may be attributed to the increased energy dissipation of the wall.Further, because the movement of the antiseepage wall is constrained by the dam during seismic response, the acceleration of the wall differs by a relatively small amount at different thicknesses.

Wall's Materials
In order to compare the differences between the dam with the polymer wall and the dam with the traditional concrete wall, dynamic fluid-solid coupling calculations were performed on earth-rock dams with polymer and concrete antiseepage walls under working conditions 1-4.The results are shown in Figure 18.The development tendency of acceleration increases with the wall height and the input seismic acceleration, and the greatest acceleration occurs at the top of the wall.The peak accelerations of polymer and concrete antiseepage walls are shown in Table 5.The trends of peak acceleration for the two types of antiseepage walls are generally in agreement.A polymer antiseepage wall's peak acceleration is lower than that of a concrete one under the same working conditions; this is because polymer materials have a lower stiffness and lighter mass than concrete materials.The greatest compressive and tensile stresses of the two types of antiseepage under working conditions 1-4 are shown in Figures 19 and 20.It can be seen that u seismic response, the polymer wall's greatest compressive and tensile stresses both at the bottom of the wall.The greatest compressive stress of the wall under is 138 kP kPa, 175 kPa, and 199 kPa, respectively, and the value of the greatest tensile stress kPa, 139 kPa, 154 kPa, and 177 kPa, respectively.By comparison, the maximum com sive stress of the concrete wall occurs at the bottom of the wall, while the maximum te stress of the concrete wall is near the bottom.The greatest compressive stresses o concrete wall are 645 kPa, 821 kPa, 1443 kPa, and 2561 kPa, respectively.The maxi tensile stress of the wall is 1004 kPa, 1069 kPa, 1242 kPa, and 1686 kPa, respectively seen from these figures that the polymer antiseepage wall experiences smaller maxi compressive and tensile stresses than the concrete wall because the modulus of elas of polymer materials is similar to the adjacent soil, and polymer materials have e meric characteristics that provide some flexibility under seismic response.The greatest compressive and tensile stresses of the two types of antiseepage walls under working conditions 1-4 are shown in Figures 19 and 20.It can be seen that under seismic response, the polymer wall's greatest compressive and tensile stresses both occur at the bottom of the wall.The greatest compressive stress of the wall under is 138 kPa, 154 kPa, 175 kPa, and 199 kPa, respectively, and the value of the greatest tensile stress is 126 kPa, 139 kPa, 154 kPa, and 177 kPa, respectively.By comparison, the maximum compressive stress of the concrete wall occurs at the bottom of the wall, while the maximum tensile stress of the concrete wall is near the bottom.The greatest compressive stresses of the concrete wall are 645 kPa, 821 kPa, 1443 kPa, and 2561 kPa, respectively.The maximum tensile stress of the wall is 1004 kPa, 1069 kPa, 1242 kPa, and 1686 kPa, respectively.It is seen from these figures that the polymer antiseepage wall experiences smaller maximum compressive and tensile stresses than the concrete wall because the modulus of elasticity of polymer materials is similar to the adjacent soil, and polymer materials have elastomeric characteristics that provide some flexibility under seismic response.The relationship between the compressive and tensile strengths of polymer materials at various densities is depicted in Figure 21.As indicated in the figure, the compressive and tensile strengths of polymer materials with a density of 260 kg/m 3 are 4540 kPa and 3350 kPa, respectively, and the maximum compressive and tensile stress experienced by the polymer antiseepage wall during an earthquake are lower than their strengths.This yields sufficient safety reserves.In comparison, the compressive and tensile strengths of the concrete material are usually 25,000 kPa and 1500 kPa, respectively.However, under working condition 4, the maximum tensile stress of the concrete antiseepage wall has exceeded its tensile strength, leading to tensile damage.The concrete damage can be predicted from Figure 20b; it can also be proved by the results of the centrifuge test shown in Figure 22.
The excavation result of the dynamic centrifuge test is shown in Figure 22.The polymer wall appears to be in good overall condition.In contrast, the concrete wall has cracks near the bottom, resulting in failure on both sides of the wall, which impacts the impermeability of the earth-rock dam.The test results are in accordance with the results of the The relationship between the compressive and tensile strengths of polymer materials at various densities is depicted in Figure 21.As indicated in the figure, the compressive and tensile strengths of polymer materials with a density of 260 kg/m 3 are 4540 kPa and 3350 kPa, respectively, and the maximum compressive and tensile stress experienced by the polymer antiseepage wall during an earthquake are lower than their strengths.This yields sufficient safety reserves.In comparison, the compressive and tensile strengths of the concrete material are usually 25,000 kPa and 1500 kPa, respectively.However, under working condition 4, the maximum tensile stress of the concrete antiseepage wall has exceeded its tensile strength, leading to tensile damage.The concrete damage can be predicted from Figure 20b; it can also be proved by the results of the centrifuge test shown in Figure 22.
The excavation result of the dynamic centrifuge test is shown in Figure 22.The polymer wall appears to be in good overall condition.In contrast, the concrete wall has cracks near the bottom, resulting in failure on both sides of the wall, which impacts the impermeability of the earth-rock dam.The test results are in accordance with the results of the The relationship between the compressive and tensile strengths of polymer materials at various densities is depicted in Figure 21.As indicated in the figure, the compressive and tensile strengths of polymer materials with a density of 260 kg/m 3 are 4540 kPa and 3350 kPa, respectively, and the maximum compressive and tensile stress experienced by the polymer antiseepage wall during an earthquake are lower than their strengths.This yields sufficient safety reserves.In comparison, the compressive and tensile strengths of the concrete material are usually 25,000 kPa and 1500 kPa, respectively.However, under working condition 4, the maximum tensile stress of the concrete antiseepage wall has exceeded its tensile strength, leading to tensile damage.The concrete damage can be predicted from Figure 20b; it can also be proved by the results of the centrifuge test shown in Figure 22.
The excavation result of the dynamic centrifuge test is shown in Figure 22.The polymer wall appears to be in good overall condition.In contrast, the concrete wall has cracks near the bottom, resulting in failure on both sides of the wall, which impacts the impermeability of the earth-rock dam.The test results are in accordance with the results of the numerical simulation, which indicates that the established numerical model is a reasonable model.

Conclusions
In this study, a dynamic microscopic fluid-solid coupling model was formul analyze the behavior of an earth-rock dam with a polymer antiseepage wall.The response of the polymer antiseepage wall was examined, and the seismic stability wall was evaluated.The conclusions of this study are as follows: (1) The peak acceleration of the antiseepage wall was observed to increase as a r the concurrent actions of reservoir water on both the stress field and seepa of the earth-rock dam.This increase in the dam's acceleration response conseq led to an elevation in the peak acceleration of the wall.Additionally, the app of microscopic fluid-solid coupling had a notable impact, resulting in a 25% a increase, respectively, in both the maximum compressive and tensile stresse rienced by the wall.Moreover, the displacement of the polymer antiseepa experienced a significant increase due to the constraints imposed by the ear dam.Consequently, it is imperative to consider the effect of fluid-solid coup the safety assessment and analysis of earth-rock dams.

Conclusions
In this study, a dynamic microscopic fluid-solid coupling model was formulated to analyze the behavior of an earth-rock dam with a polymer antiseepage wall.The seismic response of the polymer antiseepage wall was examined, and the seismic stability of the wall was evaluated.The conclusions of this study are as follows: (1) The peak acceleration of the antiseepage wall was observed to increase as a result of the concurrent actions of reservoir water on both the stress field and seepage field of the earth-rock dam.This increase in the dam's acceleration response consequently led to an elevation in the peak acceleration of the wall.Additionally, the application of microscopic fluid-solid coupling had a notable impact, resulting in a 25% and 13% increase, respectively, in both the maximum compressive and tensile stresses experienced by the wall.Moreover, the displacement of the polymer antiseepage wall experienced a significant increase due to the constraints imposed by the earth-rock dam.Consequently, it is imperative to consider the effect of fluid-solid coupling in the safety assessment and analysis of earth-rock dams.

Conclusions
In this study, a dynamic microscopic fluid-solid coupling model was formulated to analyze the behavior of an earth-rock dam with a polymer antiseepage wall.The seismic response of the polymer antiseepage wall was examined, and the seismic stability of the wall was evaluated.The conclusions of this study are as follows: (1) The peak acceleration of the antiseepage wall was observed to increase as a result of the concurrent actions of reservoir water on both the stress field and seepage field of the earth-rock dam.This increase in the dam's acceleration response consequently led to an elevation in the peak acceleration of the wall.Additionally, the application of microscopic fluid-solid coupling had a notable impact, resulting in a 25% and 13% increase, respectively, in both the maximum compressive and tensile stresses experienced by the wall.Moreover, the displacement of the polymer antiseepage wall experienced a significant increase due to the constraints imposed by the earth-rock dam.Consequently, it is imperative to consider the effect of fluid-solid coupling in the safety assessment and analysis of earth-rock dams.(2) Considering the effect of fluid-solid coupling, the peak acceleration, dynamic stress, and displacement of the polymer antiseepage wall increased progressively as the water level rose.Furthermore, the acceleration, compressive stress, and tensile stress of the polymer antiseepage wall also increased with the density of polymer materials due to the increased stiffness and seismic response of the wall.However, as the wall thickness increased, the acceleration, compressive stress, tensile stress, and displacement of the polymer antiseepage wall decreased due to the enhanced energy dissipation capacity of the wall.(3) Under the same working conditions, the peak acceleration of the polymer antiseepage wall was smaller than that of the concrete one due to the effect of microscopic fluidsolid coupling.Under the same earthquake, the compressive and tensile stresses of the polymer antiseepage wall were both smaller than their respective compressive and tensile strengths, which indicated sufficient safety reserves.while the tensile stress of the concrete antiseepage wall exceeded its tensile strength, causing the wall to be damaged.The displacement of the polymer antiseepage wall was greater than that of the concrete wall due to the flexible nature of polymer materials and their ability to self-recover from large deformations.

Figure 2 .
Figure 2. Damages of dams caused by earthquakes.

Figure 2 .
Figure 2. Damages of dams caused by earthquakes.

Figure 3 .
Figure 3. Structure and modeling of Jiulong earth-rock dam: (a) sketch of Jiulong earth-rock dam; (b) cross-section of dam model; (c) three-dimensional finite element model.

Figure 3 .
Figure 3. Structure and modeling of Jiulong earth-rock dam: (a) sketch of Jiulong earth-rock dam; (b) cross-section of dam model; (c) three-dimensional finite element model.
2) Program development The fourth-order generalized Maxwell model is developed in VC++.Development work includes modifying header files (.H) and source files (.CPP) and generating dynamic link library files (.DLL) for the main program to call.The program flow chart of the generalized Maxwell model is shown in Figure 5. Sustainability 2023, 15, x FOR PEER REVIEW 11 of 28 The fourth-order generalized Maxwell model is developed in VC++.Development work includes modifying header files (.H) and source files (.CPP) and generating dynamic link library files (.DLL) for the main program to call.The program flow chart of the generalized Maxwell model is shown in Figure 5.

Figure 5 .
Figure 5.The secondary developed generalized Maxwell model.In the diagram, the rounded rectangle represents the start or end, the rectangle represents a process, and the diamond represents a decision.

Figure 5 .
Figure 5.The secondary developed generalized Maxwell model.In the diagram, the rounded rectangle represents the start or end, the rectangle represents a process, and the diamond represents a decision.
numerical simulation, which indicates that the established numerical model is a r ble model.

Figure 21 .Figure 22 .
Figure 21.Relationship between strength and density: (a) relationship between compressiv strength and density; (b) relationship between tensile strength and density.

Figure 21 .
Figure 21.Relationship between strength and density: (a) relationship between compressive strength and density; (b) relationship between tensile strength and density.

Figure 21 .Figure 22 .
Figure 21.Relationship between strength and density: (a) relationship between compressive strength and density; (b) relationship between tensile strength and density.

Table 1 .
Mechanical parameters of materials.
so it was chosen for analysis.The energy storage modulus E′ and loss modulus E″ in the fourth-order generalized Maxwell model are as follows: The storage modulus E and loss modulus E of the generalized Maxwell model are expressed as follows:

Table 3 .
Similarity relation of earth-rock dams.