Abstract
The general solution of the axial symmetry elastic space problem in cement concrete pavement is an elementary question. Combining the Love method and the Southwell method, the Southwell operator was used in the variable selection process, the displacement function was introduced to express the displacement component by the Love method, the expression of stress indicated by the displacement function was obtained by combining the displacement components, geometric equations, and physical equations, then the stress was substituted into the equilibrium equation, the biharmonic equation of displacement function was obtained by mathematical operation, and finally, a new general solution of the axisymmetric elastic space problem in cement concrete pavement was obtained. According to the proposed general solution, mechanical calculation for a cement concrete pavement slab on Winkler foundation was carried out, the displacement curve at the top surface of the slab and the stress curve at the top, middle, and bottom surface of the slab were derived, and the results showed that the method is feasible in the calculation of cement concrete pavement. The general solution provided a new method for solving the axisymmetric elastic space problem of cement concrete pavement.
1. Introduction
The problem of axisymmetric elastic space refers to the elastic body that the geometry, constraint condition, and load are symmetric to a certain axis, and the stress, strain, and displacement of the axisymmetric elastic space problem are also symmetrical to the axis [1–3]. In fact, an axisymmetric problem is a special form of spatial problem where all variables are independent of θ and only related to r and z. Therefore, the shear stress, shear strain, and displacement along the θ direction are all zero [4, 5].
The common method for solving the general solution of axisymmetric elastic space in the pavement is usually the Love method or the Southwell method [6, 7]. Both the Love method [8] and the Southwell method [9] use the introduction of the displacement function to represent the displacement component, and by combining the displacement component represented by the displacement function, and the geometric and physical equation, then the expression for stress represented by the displacement function is obtained. Substituting the obtained stress into the equilibrium equation, a harmonic equation about the displacement function is obtained, and solving this equation gives the expression of the displacement function; substituting the obtained expression into the displacement function to represent the displacement and stress, the general solution of the axisymmetric elastic space problem is obtained [10–12].
The solution of displacement and stress of the cement concrete pavement belongs to the axisymmetric elastic space problem [13]. Scholars have achieved many results in this area [14, 15]; however, it is difficult to find the stress function for pavement [16, 17]. The author adopted the Southwell operator without variable substitution but directly adopted the idea of the Love method to obtain a general solution of the axisymmetric elastic space problem, and this general solution undergoes certain form changes to derive the general solutions of the Southwell and Love methods. According to the proposed general solution, the mechanical calculation was carried out for a cement concrete pavement slab on Winkler foundation, and the displacement curve at the top surface of the slab and the stress curve at the top, middle, and bottom surfaces of the slab were obtained to verify the feasibility of the method in the cement concrete pavement. The general solution provided a new method for solving the axisymmetric elastic space problem.
2. The Derivation of General Solution
The cylindrical coordinate system is adopted as solving the axisymmetric elastic space problem in pavement engineering. The basic equations are shown as follows.
Equilibrium equations are as follows:
Geometric equations are as follows:
Constitutive equations are as follows:where μ is Poisson’s ratio, E is elasticity modulus, u and are displacement of pavement, and is shear modulus.
The compatibility equation is as follows:where .
The stress can be obtained as follows:
Substituting equation (5) into equations (1) and (4) yields
According to equations (2), (3), and (5), the displacement component can be expressed as follows:where , is the Southwell operator.
From equations (2), (3), (7), and (8), the expression of the displacement component and stress component is as follows:
Substituting equations (9)–(12) into equation (1), displacement must be satisfied by the following equation:
According to Hankel transform, we can obtain the following equation:
According to Hankel inverse transform, equation (14) can be rewritten as the following equation:
The differential operation is performed on both sides of equation (14) r and z, and we can get the following equations:
According to equations (16)–(18), the displacement of the cement concrete pavement can be expressed as follows:
Combining equations (13) and (20), we can draw the conclusion as the following equation:
According to Hankel transform, equation (21) can be rewritten as the following equation:
General solution is constructed to solve differential equation (22).
Substituting equation (23) into equation (14) yields:
Substituting equation (23) into equations (18) and (19) yields
Substituting equations (25) and (26) into equations (7)–(11), denoted by A = ξ2Aξ, B = ξBξ, C = ξ2Cξ, D = ξDξ, yieldwhere .
3. The Relationship between General Solutions
Parameters are identified as AL, BL, CL, and DL. Denoted by A = AL + (4μ − 1)BL, B = −BL, , and , we will just get the Love solution as follows by the equation of (27):where .
Parameters are identified as AS, BS, CS, and DS. Denoted by A = As + 2μBs, B = −Bs, C = Cs + 2μDs, and D = Ds, we will just get the Southwell solution as follows by the equation of (27):where .
Denoted by As = −AL − (2μ − 1)BL, Bs = BL, Cs = CL + (2μ − 1)DL, and Ds = DL, we can convert the Southwell solution into Love solution.
Combining these results leads to the general transformation relationship between solutions as shown in Table 1.
4. Application of General Solution on Winkler Foundation
As shown in Figure 1, δ is radius of the circle, and is vertical circular uniform distributed load.
Boundary conditions are as follows:
Combining equations (26) and (31), we can get
Denoted by z = 0, from equations (30) and (31), we can get
Combining equations (30) and (33), we can get the following equation:
Substituting equations (33) and (34) into equation (26) yieldwhere x = ξδ.
5. The Proposed General Solution Apply to Winkler Foundation
5.1. Model
The arrangement of the circular uniform load on the Winkler foundation is shown in Figure 2.
Boundary conditions are as follows:
Substituting equation (26) into equations (36) and (37) yields
According to Hankel transform, we can obtain the following equation:where h is the thickness of the cement concrete slab; k is Winkler foundation model.
The solution of equation (39) can obtain the expression for A, B, C, and D about ξ, μ, E, h, and k. Substituting A, B, C, and D into equation (27), the displacement and stress of cement concrete pavement can be obtained.
5.2. Examples
The calculation diagram of cement concrete pavement is shown in Figure 3, the specific calculation parameters are E = 11500 MPa, h = 0.18 m, μ = 0.15, k = 1.45 × 107 N/m3, p = 700 kN/m2, r = 0.151 m, I choose integral interval 0–10, 0–20, 0–30, 0–40, 0–50, 0–500 and calculate displacement at point a, stress at point a, b, c, and the results are summarized in Tables 1–4.
Based on Table 2, displacement at point a is about 0.86 mm. According to Tables 3−5, stress at point a is −1.1 MPa, stress at point b is −0.028 MPa, and stress at point c is 1.09 MPa.
5.3. Variation of Displacement and Stress with r
The relative stiffness radius l of cement concrete pavement is as follows:
l ∈ (0.8, 1.5), stresses and displacement vary with cement concrete pavement board size, and the range of finite size slab instead of infinite size slab is referred to Zheng [6]:
Calculation range of r on the Winkler foundation is 0 to 4.4 m. According to the former calculation models, the displacement at point a on the top surface of the slab is calculated to vary with r as shown in Figure 4, and the stresses at points a, b, and c vary with r as shown in Figures 5−7, respectively.
From Figure 4, the displacement at point a decreases rapidly with the increase of r, and the displacement at point a gradually tends to a constant value when the r is greater than 3 m. From Figures 5–7, it can be seen that the stresses at each point occurred great changing amplitude with the increase of r. In particular, the stress of the plate changes from compressive stress to tensile stress (Figure 5) or changes from tensile to compressive stresses (Figures 6 and 7) with the increase of r.
6. Conclusions
Using the compatibility equation of the Southwell method and the calculating idea of the Love method, a new general solution of the axisymmetric elastic space problem was derived in pavement engineering. By variable transform, the general solution can be transformed into the Southwell solution or Love solution.
According to the boundary conditions of the pavement and the characteristics of the general solution, the general solution is divided into the general solution for solving the composite boundary, the general solution for the stress problem, and the general solution for the displacement problem, and these three types of problem are represented in the form of the Love general solution, the Southwell general solution, and the proposed general solution in the paper. The general solution provides a broader idea for analyzing the axisymmetric elastic space problem of cement concrete pavement.
Using the general solution proposed by the author, the mechanical calculation of a cement concrete pavement slab on Winkler foundation as an illustration is carried out, the displacement curve at the top surface of the slab and the stress curves at the top, middle, and bottom surface of the slab were obtained, and the result showed that the method is feasible in the calculation of cement concrete pavement.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest.
Acknowledgments
This work was supported by the Science and Technology Project of Gansu Provincial Department of Transportation in 2022 (GSSJTT202203).