Derivation Method for the Foundation Boundaries of Hydraulic Numerical Simulation Models Based on the Elastic Boussinesq Solution

The foundation boundaries of numerical simulation models of hydraulic structures dominated by a vertical load are investigated. The method used is based on the stress formula for fundamental solutions to semi-infinite space body elastic mechanics under a vertical concentrated force. The limit method is introduced into the original formula, which is then partitioned and analyzed according to the direction of the depth extension of the foundation. The point load will be changed to a linear load with a length of 2a. Inverse proportion function assumptions are proposed at parameter a and depth l of the calculation points to solve the singularity questions of elastic stress in a semi-infinite space near the ground. Compared with the original formula, changing the point load to a linear load with a length of 2a is more reasonable. Finally, the boundary depth criterion of a hydraulic numerical simulation model is derived and applied to determine the depth boundary formula for gravity dam numerical simulations.


Introduction
In hydraulic engineering, numerical simulation models are widely used in structural analysis [1], such as in infinite and discrete element methods.The finite element method (FEM) is the most widely used technique in hydraulic numerical simulation [2].Selecting the foundation boundaries affects the number of finite element mesh size, calculation accuracy, and computation time [3].Many scholars have proposed selection methods for foundation boundaries based on the finite element simulation model.Foundation boundaries are generally selected based on engineering experience [4] or by comparing the accuracy of stress results using FEM for different boundary conditions [5].However, these selection methods are still not applied in engineering design because of different issues in analyzed objects and calculation accuracy [6].At present, comprehensive qualitative studies remain inadequate in both China and abroad; for example, selecting calculation boundaries has not been studied based on elasticity theory [7].According to Saint-Venant's principle [8], hydraulic structure load has more influence on the stress of a point near the load and minimal influence on the stress of a point far from the load; hence, conducting a partition discussion is necessary.Research on foundation boundaries under dam structure load can analyze foundation stress under loads and search for the minimum impact areas to function as boundaries [9].Many hydraulic structure loads are produced under vertical force such as in gravity dams, water locks, and docks.The current study focuses on hydraulic structures dominated by a vertical force.The Boussinesq stress solution to a semi-infinite elastic mechanics space body under a vertical concentrated load is widely used in soil mechanics calculation [10].This study combines L'Hospital's rule and the Boussinesq equation to derive the formula for the selection principle of foundation boundaries.

Methodology
This study investigates the problem of selecting a finite element numerical model and a foundation boundary scope for hydraulic structures under a vertical load such as gravity dams and sluices.The load of hydraulic structures on the foundation consists of nonlinear vertical distributed loads resulting from a nonlinear structure shape, nonlinear materials, and construction factors [11].It can be decomposed into numerous distributed loads with different sizes (Figure 1).The uneven distribution of loads has different effects on foundation depth.The Boussinesq stress solution to semi-infinite elastic mechanics space body under a vertical concentrated load is widely used in soil mechanics calculation, particularly when the calculation point stress of the foundation is approximately zero.The coordinate of the point can determine the foundation boundaries of a hydraulic numerical simulation model.The formula, however, suffers from a singular point problem [12] near the ground.In this chapter, formula derivation is presented in three parts.Section 2.1 introduces the Boussinesq foundation stress formula for vertical concentrated loads.Section 2.2 describes the singular point problem near the ground and proposes solutions to this problem.Section 2.3 derives the stress formula under distributed loads with a parameter width of  and determines the coordinates of a point in which stress is approximately zero to obtain the foundation boundaries for an elastic numerical simulation model of a hydraulic structure.

Introduction to Common Formulas.
A hydraulic structure load on a semi-infinite elastic foundation consists of different continuous vertical distributed forces, which are decomposed into concentrated forces.Stress equations based on the half-space problem and the plane elastic mechanics problem provided by Boussinesq satisfy the objective of the aforementioned analysis problem.Suppose that a half-space has no body force and the normal concentrated force  on the plane boundary has a value of . Figure 2 shows the distribution map of a load.This problem is axisymmetric, and the coordinate origin point is the center of .The stress boundary condition is (  ) =0, ̸ =0 = 0, (  ) =0, ̸ =0 = 0.The solution is described as follows: In ( 1), (2), and (3),   is the stress in the direction of ,   is the stress in the direction of ,   is the stress in the direction of ,  is Poisson's ratio of the foundation, and  is the amount of the force.The Boussinesq solution focuses on axisymmetric problems.The stress of hydraulic structure load on a foundation is an axisymmetric problem, wherein axes  and  are replaceable when studying finite element boundaries.Equations ( 1), (2), and (3) show that 3D stresses are proportional to the square root of length factors.The vertical stress (  ) on the foundation is critical because it is the main factor that results in the settlement of buildings [13].Therefore, additional research on vertical stress (  ) to derive the foundation boundary of the finite element model is reasonable.

Singular Point of Common
Formulas.Saint-Venant's principle in elastic mechanics states that if a force system acting on a small elastic surface is replaced by another equivalent static force, then the effects on the stress and strain of this region are significant, whereas the effect on distance is negligible.The accuracy of the Boussinesq formula is high when the calculation point is far from the ground.By contrast, the vertical force increases rapidly when the calculation point is near load areas, which contradicts actual situations.Thus, the point load formula is no longer applicable near the dam and requires correction.The vertical stress concentration problem near the dam area under a point load [14,15] can be proven as follows: Given the randomness of , , and  when they satisfy the function of In the elastomeric hypothesis, the existence of stress concentration indicates the presence of a plastic expansion area.If such area exists, then the elastomeric hypothesis is contradicted.Newton solved the gravity concentration problem of two close particles by considering each particle as a sphere with a radius.To avoid stress concentration, the point load is changed to a line load with a certain length.
In Mathematical Principles and Natural Philosophy, Newton explained the law of gravitation of objects as follows.When the distance between two particles tends to be "infinitely small, " the particles are regarded as a sphere with a radius [16], and consequently, the range of the integral of the force between two objects in space can be calculated.When the particles are close to each other, the Boussinesq point stress formula presents the same problem; that is, (/ 2 ) → ∞.Thus, particle assumption is no longer true.Based on the definition of a particle provided in Newton's law of universal gravitation, the concept of a point load is relative and can be calculated based on the particle formula when two objects are far from each other.Based on this principle, the span of a load is regarded as zero when the load is far from the dam.Newton used calculus to solve problems of increasing gravity in two particles with a small distance between them.The Boussinesq vertical stress   → ∞ when the calculation point is near the load.Consequently, the point load near the dam should be changed to a line load with a certain length.

Deriving the New Equation.
To avoid the problem of stress concentration near the ground, the aforementioned method of changing the point load to a distributed load is used to derive the correction formula for the plane and space problems.The modified formula is then used to analyze the influence of the foundation scope.When the calculation point is far from the load areas, its correction width can be ignored based on Saint-Venant's principle.By contrast, when the calculation point is close to the load areas, its correction width cannot be ignored.The following sections discuss the plane and space problems.The vertical stress formula for different depths under the effect of a point load is corrected and the function of load width parameter  and depth  is proposed.The relative error between the original and the corrected formula, as well as the influence of the foundation scope, is then analyzed.0, the vertical stress of all points on line  =  0 tends to be zero.In this research, the vertical stress (  ) of point  whose distance from the original point is  is studied.The loads of hydraulic structures that are acting on the foundation are superposed by successive point loads of different sizes.

The Plane Problem.
The influence scope of an arbitrary point load can then be analyzed.Based on the preceding analysis, the point load is changed to a distributed load with a width of 2.The following section provides the hypothetical function of width parameter  and the depth of calculated point .
(1) Determining the Physical Parameters.According to Saint-Venant's principle, the farther the calculated point from the load areas, the less the error caused by the correction width parameter 2 of the distributed load and the original point load.Meanwhile, the closer the calculated point is from the load areas, the larger the width parameter 2 should be.
Assuming that the existence of critical depth  0 satisfies  ≥  0 , parameter  approximates zero and the load can be considered as a point load.By contrast, when  <  0 , the load should be considered as a distributed load.Observed through continuity, a tapering function between  and  exists.
Assuming that parameter  is inversely proportional to , then The three situations of the calculated point depth  are as follows.
A The calculation point is far from the surface with B The calculation point is close to the surface with  ̸ = 0, / 2 → ∞.
(2) Correction Formula and Comparison.The following formulas correct the vertical stress formula for the three situations of the calculation point depth , respectively, and compare precision before and after correction.
(1)  is extremely large; thus, / → ∞: with  =  tan , where  = /2, The result of the original formula is  1 = 3/2 The error is merely  2 /3 2 .Given that / → ∞, the error is approximately zero, which illustrates that the Boussinesq solution has good applicability.Hence, the error caused by the width of the distributed load can be ignored.
The integration results of the distributed load are as follows: Thus,   is bounded.In fact, as long as  ×  is bounded, no singular point will exist.
Apply the Taylor series at  = 0; then sin () =  −  3  6 +  ( 3 ) , The result of the original formula is  1 = 3/2 Given that / is extremely large when the calculation point is near the ground, the error caused by the width of the distributed load cannot be ignored.The stress concentration problem can be solved in this area, which demonstrates that the width of the distributed load cannot be ignored. ( When searching for the intersection of  1 = sinarctan  − (1/3)sin 3 arctan  and  2 =  using the software MATLAB, the result is When  < 0.106, the relative error and the width of the distributed load can be ignored.By contrast, when  > 0.106, the relative error cannot be ignored and the point load must be changed to a distributed load.Therefore, for a plane problem, the error caused by the width of the distributed load can be ignored only if / < 0.106, and the Boussinesq point stress solution can still be used.

Space Problem.
In a space problem, the point load is changed to a distributed load with a radius of .
(1)  is the radius of the distributed load.When  is extremely large, / → ∞: The result of the original formula is The preceding equations show that when  is extremely large, the error caused by the radius of the distributed load can be ignored.
(2)  is extremely small, and thus  ̸ = 0, / → ∞.Following the same analysis method, Let Following (20), when  is extremely small, the error caused by the radius of the distributed load is large.

Applications
The Three Gorges Dam located in Hubei Province, China, is the largest hydropower project in the world.The dam has two parts: the main building and the diversion.Three Gorges is a concrete gravity dam that is 3035 m long, 185 m high, 15 m wide at its crest, and 124 m wide at its bottom.Its installed capacity of 22.5 million kW ranks first in the world.In this section, the method proposed in Section 2 is first used to analyze the boundary range of the Three Gorges Dam through a numerical model.Then, the finite element numerical model of the Three Gorges Dam is established to illustrate the validity and efficiency of the proposed formulas.

Deriving Foundation Size.
Hydraulic structure loads such as in gravity dams, water locks, and docks are dominated by a continuous vertical concentration.Based on the conclusions for the plane and space problems, one issue arises, that is, the calculated value of the vertical stress (  ) at point  is small when width parameter  and calculation point depth value  satisfy the condition of / < 0.106.
The objective of the present study is to explore the boundaries of the finite element model selection problem for vertical concentration loads, which is equivalent to finding the coordinates of the points when the base stress is small.In the following section, the base boundary depth of the finite element model is analyzed using the preceding conclusions.
Because of the limit of influenced scope based on the hydraulic structure loads, so the foundation boundary belongs to the result of (11).Equation (11) shows the stress results of the definite integral as follows: Its depth boundary is  max = 4.717 = 3.16.4 shows the finite element model of the Three Gorges Dam with the following specifications: the -axis along the downstream direction, the -axis along the dam axis, and the -axis along the vertical direction.The finite element model is divided into 5393 units and 11102 nodes.The FEM material parameters are presented in Table 1.The actual dimensions of the finite element model of the Three Gorges Dam are provided in Tables 2 and 3.
Based on the aforementioned calculation principle and parameter, the specific calculation results and analyses are as follows.Given that the foundation and dam body have the same size and material along the -axis direction, the stress results of the plane and space problems along the -axis section are also the same.The following finite element model is focused on plane problems.Figure 5 illustrates the vertical stress (  ) calculation results of the finite element model of the Three Gorges Dam.All size parameters of the finite element method are chosen for the international system of units: meters.The legend of Figure 5 is the size range of   (in Pascal).This figure can indicate that the foundation vertical stress (  ) is approximately zero when depth is deeper than the depth of point .The depth of point  is the furthest depth of the load effect.Figure 5 shows the depth of point ; that is,   = 556 m.Dam height is 185 m; that is,   = 556/185 = 3.01.
The boundary range acquisition method is based on the Boussinesq formula.This formula is suitable for a semiinfinite body problem under a concentrated load stress.The finite element model of the Three Gorges Dam selected 1 to be the foundation range upstream and downstream and 4 to be the foundation depth range.According to the finite element calculation results, the foundation influence scope is limited on the loads of hydraulic structures and expanding the scope of the foundation boundary has minimal effect on calculation precision.The Boussinesq formula accurately derives the elastic mechanics formula based on a differential equation, but FEM is the solution to a partial differential equation.Moreover, numerical solution techniques are approximated through the variation method, and the error function is set to minimum values to generate a stable solution.FEM is an approximate method for solving the Boussinesq half-space problem.The two methods, which have different principles, produce slightly different results; however, they exhibit the same trend.The depth boundary defined in Section 3.1 is  = 3.16.Its relative error is less than 5%; therefore, the model exhibits high precision.Furthermore, the model proves that (23) can be used as

Figure 1 :
Figure 1: Decomposing the structure of a load.

Figure 2 :
Figure 2: Distribution map of a load.

Figure 3 :
Figure 3: Load distribution map of the plane problem.
3)  is moderate, which corresponds to depth situation C: tends to be equal to zero.Parameter  is restricted by the size of the structures.Assume that width  and height  of hydraulic structures satisfy the equation  = .Then, the deeper  is, the smaller parameter  is, which results in a smaller calculation error of the Boussinesq equation.The maximum value of load length 2 is the width  of the structure, which corresponds to the maximum depth as follows:  ) tends to be zero, and  is equivalent to the depth boundary of the finite element calculation model.Various hydraulic structures have different values for parameter .The value of  should be calculated for each project, and then the base depth boundary of the finite element numerical simulation model can be calculated.The preceding conclusions can be applied to a specific project.The Three Gorges Water Conservancy Project is the largest gravity dam project in China.In this study, the depth boundary of the finite element numerical simulation model of the aforementioned project is analyzed.The Three Gorges Dam Project is 185 m high and 124 m wide.The parameter  of Three Gorges is

Table 2 :
FEM dimensions of the dam body.

Table 3 :
FEM dimensions of the foundation.