Research on the identification of modulus parameters for rigid pavement panel foundations

Abstract: Based on the Pasternak two-parameter foundation model, this study establishes a dynamic model for analyzing the deflection of slabs under moving loads. The variational method and reciprocal equal work theorem are employed to obtain the deflection solution, while the Matlab program is utilized with practical examples to calculate the foundation response modulus considering transfer shear force, bending moment, and their combined effects. By fitting measured dynamic deflections with theoretical predictions using the least square method, the response modulus K of the foundation is determined, demonstrating that our proposed model effectively captures the realistic behavior of rigid pavements.

Keywords: Pasternak model; Least squares method; Foundation reaction modulus

1 Introduction

With the rapid development of China's civil aviation industry, higher performance and safety requirements for airport pavement are inevitable. Currently, the commonly used foundation models include the E. Winkler foundation model and the elastic half-space foundation model. Although the Winkler foundation model theory is simple and intuitive, requiring only one elastic parameter to express ideal elastic foundation characteristics and relating displacement of any point solely to applied stress at that point, it fails to consider soil-ground continuity. In contrast, the elastic half-space foundation model overemphasizes stress diffusion between soil and ground, making calculations more cumbersome and impractical for engineering applications. Consequently, mechanics experts have developed several more reasonable foundation calculation models through extensive experimentation. Zhang Weigang and Chen Xiangming et al utilized the Kelvin viscoelastic model under moving load conditions to obtain a relatively accurate foundation reaction modulus [1,2]. The viscoelastic Winkler dynamic model verified by Xing Yaozhong better reflects actual airport operations [3]. V. Patil employed Pasternak's two-parameter soil medium as a model to investigate material nonlinearity effects on pavement response in supporting soil medium [4]. Yakshansh Kumar applied a novel finite element-based cyclic response model for rigid pavement design while establishing the period range of moving load within rigid pavement parameters [5]. Airport pavements consist of multiple concrete panels working collectively. This paper selects the Pasternak two-parameter foundation model based on independent parameters: foundation reaction modulus K and shear modulus G; thereby establishing a dynamic rigid pavement model considering joint load transfer capacity.

2.1 The formulation of the equation of motion

Figure 1 The schematic of a two-parameter foundation model

Figure 2 The forces on the shear layer

Based on the two-parameter foundation, a dynamic model of rigid thin plate pavement on the two-parameter foundation is established, considering boundary shear and bending under moving load [6]. The dynamic equation of the two-parameter foundation incorporates two independent parameters: reaction modulus K and shear modulus G. Figure 1 illustrates the schematic diagram of the two-parameter foundation, while Figure 2 depicts the forces acting on the shear layer [7].

According to the Kirchhoff thin plate theory [8] and the theory of elastic mechanics, we establish the motion equation for plates subjected to a two-parameter foundation under dynamic loading as follows:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): D{\nabla }^4w+M\frac{{\partial }^2w}{\partial t^2}+ Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): Kw-G{\nabla }^2w=P

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \left(1\right)

D — the flexural stiffness of the plate, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): D=\frac{Eh^3}{12(1-{\mu }^2)}

µ —Poisson's ratio

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): M=\rho h

—The areal density

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\nabla }^2=\frac{{\partial }^2}{\partial x^2}+\frac{{\partial }^2}{\partial y^2}

—The Laplace operator in two dimensions

K—Foundation reaction modulus

G—Foundation shear modulus

The aforementioned equation represents the differential equation of motion for the track panel system subjected to a general moving load, denoted as Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): P(x,y,t)

. In the context of the investigated form of moving load in this study, the equality  Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): P(x,y,t)
can be expressed as follows

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): P(x,y,t)=P_0\delta (x-v_xt)\delta (y-\eta )

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \left(2\right)

The function Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \delta

is represented by  Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \delta (.)
, and the uniform gliding speed of the aircraft is denoted by  Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): v_x
.

2.2 Boundary conditions (shear and bending)

The plate dimensions are defined as length a and width b, with simultaneous transfer of shear force and bending moment between the plates in the actual working state. The boundary shear force is generated by the spring support force, where the elastic support rigidity coefficient Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): k_1

represents the boundary condition parameter. Similarly, the bending moment is generated by the spring-supporting moment, with the elastic supporting rigidity coefficient  Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): k_2
representing another boundary condition parameter. Assuming that  Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): Y=0,b
corresponds to the side where the seam is located, we can express its general formula for boundary conditions as follows:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \begin{array}{c} V_{y=0,b}=D{\left[\frac{{\partial }^3w}{\partial y^3}+\frac{(2-\mu ){\partial }^3w}{\partial x^2\partial y}\right]}_{y=0,b}=-k_1{\left[w\right]}_{y=0,b}\\ M_{y\vert y=0,b}=D\left(\frac{{\partial }^2w}{\partial y^2}+\mu \frac{{\partial }^2w}{\partial x^2}\right)\vert {}_{y=0,b}=-k_2{\vert \frac{\partial w}{\partial y}\vert }_{y=0,b} \end{array}

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \left(3\right)

3 The solution to the equation of motion

3.1 Solving time factors

The variational formulation of equation (1) for plate motion is as follows

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\int }_0^a{\int }_0^b\left\{D{\nabla }^4w+M\frac{{\partial }^2w}{\partial t^2}+\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \left. Kw-G{\nabla }^2w-P\right\}\delta wdxdy=0

 Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \left(4\right)

The lateral load concentration stands out among them. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): P=P(x,y,t)

， Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): w=w(x,y,t)

Select a function in the specified format as the solution for the variational equation mentioned above.

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \lbrace \begin{array}{c} w(x,y,t)=T(t)W(x,y)\\ P(x,y,t)=B(t)W(x,y) \end{array}

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \left(5\right)

The mode function of the plate is represented by Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): W(x,y)

, and the coefficients of the deflection function  Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): w(x,y,t)
and load  Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): P(x,y,t)
are denoted as  Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): T(t)
and  Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): B(t)
respectively. The variation corresponding to the deflection function  Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): w(x,y,t)
is as follows:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \delta w=W(x,y)T(t)

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \left(6\right)

Substituting equations (5) and (6) into (4) gives:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \overset{..}{T}(t)+w_{mn}{}^2T(t)=\frac{1}{M}B(t)

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \left(7\right)

Where:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \begin{array}{c} w_{mn}{}^2=\frac{1}{M}\frac{I_1}{I_2}\\ I_1={\int }_0^a{\int }_0^b\left[D{\nabla }^4W(x,y)+KW(x,y)-G{\nabla }^2W(x,y)\right]W(x,y)dxdy\\ I_2={\int }_0^a{\int }_0^bW^2(x,y)dxdy \end{array}

The computable expression for T(t) is as follows:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): T(t)=(c_1cosw_{mn}t+c_2sinw_{mn}t)+\frac{1}{Mw_{mn}}{\int }_0^tB(\tau )sinw_{mn}(t- Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \tau )d\tau

 Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \left(8\right)

The natural frequency of undamped free vibration of the thin plate, denoted as Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): w_{mn}

is considered. In equation (8), the first term represents the intrinsic vibration characteristics of the sheet, while the second term accounts for its response under external excitation.

Where c₁ and c₂ are constants, B(t) is dependent on the initial conditions of motion and the characteristics of the load. In this study, we consider a concentrated load P₀, irrespective of its mass. Consequently, B(t) can be expressed as follows:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): B(t)=\frac{P_0W(v_xt,\eta )}{{\int }_0^a{\int }_0^bW^2(x,y)dxdy}

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \left(9\right)

According to formulas (8) and (9), the dynamic deflection of the plate under the influence of a moving concentrated load P, with zero initial conditions (c₁ =c₂), can be determined as follows:

Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): w\left(x,y,t\right)=\sum_{m=1}^{\infty }\sum_{n=1}^{\infty }\frac{W\left(x,y\right)}{Mw_{mn}{\int }_0^a{\int }_0^bW^2\left(x,y\right)dxdy}{\int }_0^tP_0W\left(v_x\tau ,\eta \right)sinw_{mn}\left(t-\right. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \left. \tau \right)d\tau

 Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \left(10\right)

3.2 Solving coordinate factors