Steady-State Response Analysis of the Incompressible Nonlocal Saturated Poroelastic Beam under a Vertical Harmonic Load

Based on the nonlocal theory and the theory of saturated porous media, the mathematical and physical model and governing equations of the steady-state response of the incompressible nonlocal saturated poroelastic beam under vertical harmonic loading are established with assumption of the movement of the liquid-phase fluid only in the axial direction of the beam and considering the nonlocal effects such as particle size, pore size, and pore dynamic stress. (e dynamic response of a saturated poroelastic cantilever beam with permeability at both ends under a vertical harmonic concentrated force at the free end is studied. In the frequency domain, the analytical expressions of deflection amplification factor and equivalent couple amplification factor of liquid fluid pressure are given.(e effects of nonlocal coefficient τ, mechanical parameter α, and geometric parameter β on the deflection amplification factor and equivalent couple amplification factor at the midpoint of the nonlocal saturated poroelastic cantilever beam are studied. (e results show that the steady-state vibration of the incompressible nonlocal saturated poroelastic cantilever beam has resonance. When the nonlocal effect is considered, the deflection amplification factor and the equivalent couple amplification factor are larger, so the influence of the nonlocal effect on the steady-state response of the beam should not be ignored. (e geometric parameter β has significant effect on the peak positions of the curves of the deflection amplification factor and the equivalent couple amplification factor varying with frequency.


Introduction
Saturated porous structures are widely used in civil engineering, aviation, transportation, and other engineering fields because of their good sound absorption and energy consumption. e research on the mechanical behavior of saturated porous structures has significant engineering application value and academic value. erefore, since Biot put forward the theory of saturated porous media [1,2], the research on vibration and wave propagation in saturated porous media has attracted the attention of many scholars [3][4][5][6][7][8][9]. Since Biot theory assumes that the wave length of saturated porous media is larger than the pore size, the influence of pore size effect on wave propagation is not considered and the influence of pore size is very significant at high frequency [10]. As a method of engineering description, Biot theory has some defects in the model [11] and it is difficult to determine the elastic constant and inertial coupling coefficient in the Biot dynamic equation. e porous media theory based on the mixture axiom and the concept of volume fraction is more rigorous and reasonable. e concept of volume fraction in porous media theory can avoid the incompatibility of the model and can also consider the effects of dynamic, material, and geometric nonlinearity. Aiming at the static and dynamic problems of saturated porous structures, Busse et al. [12] established the dynamic mathematical model of the saturated poroelastic Mindlin plate based on Biot theory and gave the virtual work principle with rotation angle, deflection, and pore stress as basic unknowns. Iesan and Nappa [13] studied the thermal stress of poroelastic cylinders and circular tubes. Li et al. [14] studied the transverse vibration of fluid saturated poroelastic beams and obtained the solutions of free vibration and forced vibration. Cederbaum et al. [15] studied the mechanical behavior of saturated poroelastic beams and rods and pointed out that Mandel-Cryer phenomenon also exists in the deformation response of poroelastic beams and rods. Xiao and Chen [16] established the nonlinear governing equation of large deflection of the incompressible saturated poroelastic beam under the constraint that the pore fluid only diffuses along the axial direction of the beam and studied the nonlinear bending problem of the saturated poroelastic cantilever beam with fixed impermeability at one end and free permeability at another end under a sudden constant transverse concentrated load at the free end. However, in the analysis of mechanical behavior of saturated porous structures, the scale effects such as solid skeleton particle size, pore size, and pore stress are not considered. As the engineering structure is developing towards miniaturization and intelligence, it is more and more important to consider the influence of structural scale effect on mechanical properties. e nonlocal theory in [17,18] believes that the stress at a point in the continuum is related not only to the strain at that point but also to the strain and deformation history of all points in the area near the point. e nonlocal theory considering the scale effect is mainly used to study the mechanical problems of nanostructures. In recent years, scholars began to combine the nonlocal theory with the Biot saturated soil theory [19,20] and apply it to the field of geotechnical engineering. In this paper, the dynamic control equation of the incompressible nonlocal saturated poroelastic beam under a vertical simple harmonic load is established by combining the porous medium theory and nonlocal theory and the steady-state response of the saturated poroelastic cantilever beam with permeable ends under simple harmonic concentrated force at the free end is studied.

Mathematical Model and Basic Equation
e incompressible two-phase saturated poroelastic cantilever beam under the action of vertical simple harmonic force shown in Figure 1 is studied, where i is an imaginary unit, ω is the frequency of the simple harmonic force, Q 0 is the amplitude, and the surface of the cantilever beam is impermeable. If the mechanical properties of saturated poroelastic beams are described by using the porous medium theory without considering the influence of volume force, the governing equation is Here, σ SE is the effective stress of the solid skeleton, p is the pressure of the liquid phase, ρ S is the macroscopic density of the solid phase, ρ F is the macroscopic density of the liquid phase, κ is the solid-liquid phase coupling coefficient, u S is the displacement tensor of the solid phase, v FS represents the velocity of the liquid phase relative to the solid phase, and n F is the volume fraction of the liquid phase, from which the volume fraction of the solid phase is n S � 1 − n F . e nonlocal theory proposed by Eringen holds that the stress at a point depends not only on the strain at that point but also on the weighted average of all point strains in the whole region. Here, considering the influence of nonlocal effects such as particle size, pore size, and pore dynamic stress of the solid skeleton of the saturated poroelastic beam, it is assumed that the solid skeleton of the saturated poroelastic beam is linear elastic, uniform, and isotropic material. Without the influence of physical force, the basic equation of Eringen's nonlocal theory is [17] where σ ij is the nonlocal stress tensor, ε ij and u i are the nonlocal strain and displacement tensor, respectively, C ijkl is the fourth-order elastic tensor, and V χ(|r − r ′ |, τ) is the nonlocal modulus, which determines the weight coefficient of the corresponding point strain, where τ is the nonlocal kernel function. e influence of the pore dynamic stress, pore size, and particle size of the saturated poroelastic beam is reflected by the nonlocal kernel function; |r − r ′ | is the Euclidean distance. At the same time, Eringen also gave the differential formula of Equation (4), that is, where σ SE ij and σ LSE ij represent the nonlocal stress tensor and classical stress tensor, respectively. e 0 represents the nonlocal material constants, a is the internal characteristic length which is related to solid skeleton particle size and pore size, and ∇ 2 is the Laplace operators. e constitutive equation of solid skeleton under small deformation is where σ LSE is the classical effective stress of the solid skeleton, ε S is the strain tensor of phase skeleton, I is the identity matrix, and λ S , μ S are the Lame constants.

Dynamic Control Equation of the Incompressible Nonlocal Saturated
Poroelastic Beam e dynamic problem of the incompressible nonlocal saturated poroelastic beam under a vertical simple harmonic Advances in Civil Engineering load is one-dimensional, and the displacement of the saturated poroelastic beam is where w S (x, t) is the vertical displacement of the solid phase of the saturated poroelastic beam. Assuming that the liquid phase moves only along the axial direction, the velocities in x, y, and z directions satisfy Assuming a small amount of stress σ SE zz � σ SE yy ≈ 0, it can be obtained from equations (6) and (7) that where ε S kk is the volume strain of the solid skeleton, σ SE xx and σ SE xz are the normal stress and shear stress on the cross section of the saturated poroelastic beam, respectively, and E S and G are the elastic modulus and shear modulus of the solid phase of the saturated poroelastic beam, respectively. λ S and μ S are the lame constants and satisfy Let Here, M SE x and Q SE x are the solid bending moment and shear force on the cross section of saturated poroelastic beam, respectively, and M p is the equivalent couple of the pore fluid pressure.
Ignoring the rotational inertia of the pore liquid to the solid skeleton, multiplying the equation in the ox direction of equation (1) by z and integrating it on the cross-sectional area A of the beam, one obtains If the equation in the oz direction of equation (1) is integrated on the cross-sectional area A of the beam, we have For the rotational inertia of the solid skeleton, multiplying the equation in the ox direction of equation (2) by z and integrating it on the cross-sectional area A of the beam, we obtain It can be obtained from equations (3) and (10) that Both ends of equation (18) are multiplied by z and integrated on the cross-sectional area A of the beam, and we have It can be obtained from equations (17) and (19) that Multiplying both ends of equation (11) by z and integrating it on area A, we have where I � J A z 2 dy dz is the moment of inertia of the cross section. Integrating both ends of equation (12) on area A, we can obtain Calculating 1 − (e 0 a) 2 ∇ 2 on both ends of equations (15) and (16), respectively, and considering equations (21) and (22), we obtain Advances in Civil Engineering 3 It can be obtained from equations (20) and (21) that Sorting out equations (24)-(26), one obtains Equations (25) and (26) are the vertical dynamic control equations of the incompressible nonlocal saturated poroelastic beams.

Steady-State Response of the Incompressible Nonlocal Saturated Poroelastic Cantilever Beam under a Vertical Simple Harmonic Load
Here, taking the steady-state response of the incompressible nonlocal saturated poroelastic cantilever beam shown in Figure 2 as an example, the same method can be used to analyze the steady-state response of the incompressible nonlocal saturated poroelastic cantilever beam with other boundaries. A simple harmonic concentrated force Q(t) � q 0 e iωt is acting on the free end of an incompressible nonlocal saturated poroelastic cantilever beam. e corresponding boundary conditions are

Advances in Civil Engineering
Under the action of the simple harmonic load, the incompressible nonlocal saturated poroelastic cantilever beam vibrates in a steady state. Considering the harmony of the problem, the parameters meet w S � w S e iωt , We introduce the dimensionless variables and parameters as follows: e parameter α mainly reflects the influence of mechanical parameters such as liquid-solid coupling coefficient and elastic modulus on the steady-state response of the beam, the parameter β mainly reflects the influence of geometric parameters such as the size of the saturated porous cantilever beam. Substituting the above dimensionless quantity into dimensionless equations (30)-(32), we obtain the following initial value problem: where α 1 � iωα + ω 2 τ 2 β/1 − iωατ 2 and β 1 � − βω 2 /1 − iωατ 2 . By solving equations (32) and (33), the general solutions are where 5 , and A 6 are undetermined coefficients. e undetermined coefficients can be determined by boundary condition equations (34) and (35). e deflection amplification factor g, w and the equivalent couple amplification factor of liquid pore pressure are introduced as follows: (37)

Numerical Examples
Based on the obtained deflection amplification factor w and equivalent couple amplification factor g of the incompressible nonlocal saturated poroelastic cantilever beam under vertical simple harmonic concentrated force, the steady-state dynamic response characteristics of the incompressible nonlocal saturated poroelastic cantilever beam are analyzed through numerical examples. If there is no special description, the dimensionless quantities are τ � 0.05, α � 20, and β � 10. Figures  3-8 show the curves of mid-span deflection amplification factor w and equivalent couple amplification factor g of the incompressible nonlocal saturated poroelastic cantilever beam varying with dimensionless frequency ωL/v for different values of nonlocal coefficient τ, mechanical parameter α, and geometric parameter β. Evidently, there are obvious peaks and troughs in the curves, which shows that there is resonance in the vertical vibration of the incompressible nonlocal saturated poroelastic cantilever beam under a simple harmonic load. When the frequency is large, because the structure has no time to reflect, the deflection amplification factor of the cantilever beam gradually approaches zero, while the equivalent couple amplification factor tends to a stable value. Figures 9-14 show the variation curves of the deflection amplification factor and Advances in Civil Engineering equivalent couple amplification factor at different positions of the beam varying with x/L when the dimensionless frequency ωL/v � 1.0 for different values of nonlocal coefficient τ, mechanical parameter α, and geometric parameter β. Since the external concentrated force acts at x/L � 1, the deflection amplification factor increases gradually with the increase in x/L and the equivalent couple amplification factor in the middle of the span is the largest, which has a certain relationship with its boundary conditions. e influence of nonlocal coefficient τ on the deflection amplification factor and equivalent couple amplification factor is shown in Figures 3, 4, 9, and 10; τ � 0 is the result of ignoring the influence of nonlocal effects such as solid particle size and pore size. e deflection amplification factor and equivalent couple amplification factor are larger when nonlocal effects are considered than when they are not considered. erefore, the influence of nonlocal effects on the steady-state response of the beam should not be ignored, and the deflection amplification factor and equivalent couple amplification factor increase with the increase of nonlocal factors. In addition, when τ is small, because the left end of the beam is a fixed end, the influence of the nonlocal coefficient τ on the deflection amplification factor and the equivalent couple amplification factor is relatively small. e influence of the mechanical parameter α on the deflection amplification factor and equivalent couple amplification factor is shown in Figures 5, 6, 11, and 12. When the dimensionless frequency is small (ωL/v ＜ 1.0), the larger the parameter α is, the smaller the deflection amplification factor and the larger the valley value will be. When the dimensionless frequency is large (ωL/v ＜ 1.0), the larger the parameter α is, the larger the deflection amplification factor and the greater the peak value will be. e equivalent couple amplification factor increases with the increase in parameter α. But parameter α has little effect on the peak position of the curve of the deflection amplification factor and equivalent couple amplification factor varying with frequency. e influence of the geometric parameter β on the deflection amplification factor and equivalent couple amplification factor is shown in Figures 7, 8, 13 and 14.
e larger the parameter β is, the larger the deflection amplification factor and the smaller the valley value will be, when the dimensionless frequency is small. e larger the parameter β is, the smaller the deflection amplification factor and the smaller the peak value will be, when the dimensionless frequency is large. When the dimensionless frequency is small, the equivalent couple amplification factor increases with the parameter β; when the frequency is large, it shows the opposite, but the difference is relatively small. Different from the parameter α, the parameter β has a significant influence on the peak position of the curve of the deflection amplification factor and equivalent couple amplification factor varying      with frequency. With the increase in the parameter β, the frequency corresponding to the peak value of the curve of the deflection amplification factor and equivalent couple amplification factor will increase.

Conclusions
Combining the nonlocal elastic theory with the porous medium theory and considering the nonlocal effects such as solid skeleton particle size, pore size, and pore dynamic stress, the dynamic control equation of the incompressible nonlocal saturated poroelastic cantilever beam is established. Taking the saturated porous cantilever beam as an example, the response of the beam under the vertical harmonic concentrated force at the free end is studied. e main conclusions are as follows: (1) there is resonance in the vertical vibration of the incompressible nonlocal saturated poroelastic cantilever beam under a simple harmonic load; (2) the deflection amplification factor and equivalent couple amplification factor with nonlocal effect are larger than those without nonlocal effect, so the influence of the nonlocal effect on the steady-state response of the beam should not be ignored; (3) the mechanical parameter α and geometric parameter β have a significant influence on the deflection amplification factor and the equivalent couple amplification factor, and the influence law has a certain relationship with the frequency; (4) the parameter β has a significant influence on the peak position of the curve of the deflection amplification factor and equivalent couple amplification factor varying with frequency, and α has little effect on the peak position.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e author declares no conflicts of interest.