Research Article  Open Access
Ali Ghorbanpour Arani, Zahra Khoddami Maraghi, Mehdi Khani, Iman Alinaghian, "Free Vibration of Embedded Porous Plate Using ThirdOrder Shear Deformation and Poroelasticity Theories", Journal of Engineering, vol. 2017, Article ID 1474916, 13 pages, 2017. https://doi.org/10.1155/2017/1474916
Free Vibration of Embedded Porous Plate Using ThirdOrder Shear Deformation and Poroelasticity Theories
Abstract
This research aims at studying free vibration of rectangular plate made of porous materials in which Yfoam, Gfoam, and are used and compared with each other. To obtain the Biot formulation of the constitutive equations for a porous material, linear poroelasticity theory is used. Young modulus and density of porous plate are different in transverse direction versus porosity. In order to increase the accuracy of results in comparison with classical plate and firstorder shear deformation theories, Reddy’s theory was utilized in this research. Besides, five coupled equations of motion have been studied using Hamilton’s principle and are solved by differential quadrature method (DQM). Detailed results of this study show the significant effect of aspect ratio, thickness ratio, boundary conditions, and porosity on dimensionless frequency and deflection of porous plate. Results of this study can contribute to the design of pneumatic conveying, handling, and control systems.
1. Introduction
Porous materials are materials which can be both organic and inorganic and contain pores (cavities, channels, and interstices) which are deeper than they are. These materials are mainly used in energy absorbing systems, porous electrodes, sound absorbers, filters, insulating materials, heat exchangers, construction materials, electromagnetic shielding, membrane, and membrane support [1].
This paper focuses on the vibration analysis of porous plate regarding its application in mechanical systems. In what follows, some of the works related to this issue are reviewed. Biot and Willis [2] studied the methods of measurement to determine the elastic coefficients based on the theory of the deformation of a porous elastic solid containing a compressible fluid. They also discussed the physical interpretation of the coefficients in various alternate forms. Mansutti et al. [3] presented steady flows of nonNewtonian fluids using a porous plate with suction or injection. In this study, the flow is the counterpart of the classical “asymptotic suction” problem, within the context of the nonNewtonian fluid models. Numerical solutions were compared with the solutions from a regular perturbation approach, as well as a singular perturbation. Leclaire et al. [4] presented the vibrations of a rectangular porous plate described by two coupled equations involving the time and space derivatives of the deflection and of the relative fluidsolid motion. Using Galerkin’s variational method and classical theory of plates (nonporous), the equations were solved. The RayleighRitz decomposition on the basis of orthogonal or quasiorthogonal Eigenfunctions was used and the solutions were inserted in the vibrational equations. Results of the study showed that if the viscoelastic damping of the solid is accounted, comparison between experimental and theoretical results on the plate defection is conveniently good. A simple model of the transverse vibrations of a thin rectangular porous plate saturated by a fluid was studied by Leclaire et al. [5]. The classical theory of homogeneous plates and Biot’s stressstrain relations in an isotropic porous medium with a uniform porosity was used. They come to this conclusion that the resonance frequency is an increasing function of porosity and permeability and a decreasing function of tortuosity. A comparative study of Biot’s theory and the linear theory of porous media for wave propagation problems was presented in Schanz and Diebels’s article [6]. Regarding the results, a positive correlation between two approaches was observed in case of incompressible constituents, whereas in case of compressible constituents large differences appeared. Debowski and Magnucki [7] explored the dynamic stability of a porous rectangular plate to study an axial compressed porouscellular rectangular plate which is a generalization of sandwich or multilayer plates. The result of this study shed light on the nonlinear hypothesis of deformation of the plane cross section for porouscellular plate that includes the linear hypothesis for homogenous plate and the shear deformable effect. They proved that the mathematical model of the dynamic stability of the porouscellular rectangular plate under pulsating compression load could be reduced to Mathieu’s equation. Khorshidvand et al. [8] studied buckling analysis of radially solid circular plate made of porous material bounded with the layers of piezoelectric actuators. The results of their study showed that increasing the porosity decreases the critical load and the plate will be unstable. They also concluded that the plate with high porosity with the piezo layer has a better effect on stability of the plate in comparison to increasing the thickness of the porous plate. Reducing the porosity makes the behavior of the plate tend to incline to homogeneous/isotropic behavior. Jabbari et al. [9] investigated the buckling analysis of radially loaded solid circular plate made of porous material. Properties of their plate varied across the thickness. They derived the equilibrium and stability equations through the variational formulation based on Sander’s nonlinear straindisplacement relation. In this study, the pores of the porous plate were saturated with a fluid. They come to the conclusion that increasing porosity decreases buckling load and the plate will be unstable by then. They also found that the Monotonous porosity is more unstable than symmetric and nonsymmetric porosity and the critical buckling load will be reduced by increasing the compressibility of fluid within the pores. An analysis of buckling of thin circular FG plates made of saturated poroussoft ferromagnetic materials in transverse magnetic field was done by Jabbari et al. [10]. The energy method based on the classical plate theory with the assumption of powerlaw composition for the constituent materials was used in the research. The research indicated that increasing the porosity makes the critical magnetic field decrease, the plate will be unstable, and the critical magnetic field will increase by increasing the thickness. Finally they concluded that Monotonous porosity is more unstable than the symmetric and nonsymmetric porosity in the aforementioned research. Joubaneh et al. [11] studied the thermal buckling analysis of solid circular plate made of porous material bounded with piezoelectric sensoractuator patches in her research. The basic assumption under her research was that porous material properties vary with respect to thickness. The governing equations were obtained for clamped boundary condition based on mechanical nonlinear equilibrium and linear stability equations. The effects of porous plate’s thickness, porosity, and piezoelectric thickness on thermal stability of the plate were investigated and studied in the research. A detailed analysis of complicated and general problem of a functionally graded (FG) porous variable thickness circular plate subjected to nonaxisymmetric and nonuniform shear and normal tractions was done by Behravan Rad and Shariyat [12]. In their research, the plate was supported by a nonuniform Kerr elastic foundation based on the differential quadrature and state space vector techniques in the radial and thickness directions, respectively. The results showed that the induced magnetic quantities are mainly affected by the radial displacement component. They concluded that the material properties index, porosity, and the elastic foundation have positive effect and the compressibility parameter has an adverse and negative effect on the rigidity of the plate. Thermomechanical vibration characteristics of FG Reddy beams made of porous material subjected to various thermal loadings were investigated by Ebrahimi and Jafari [13]. They discussed the influences of several important parameters such as powerlaw exponents, porosity distributions, porosity volume fractions, and thermal effects on natural frequencies of the temperaturedependent FG beams.
Bending and vibration of FG beams were analyzed with a simple and refined trigonometric higherorder beam theory by Bourada et al. [14]. In this theory the thickness stretching effect () was also included and the transverse shear deformation effects were considered without requiring shear correction factors. In this paper the material properties varied in the thickness direction. Numerical results of the present theory were compared with other theories to show the effect of the inclusion of transverse normal strain on the deflections and stresses.
Bending and free vibration analysis of FG plates were developed using an original hyperbolic sine shear deformation theory by Neves et al. [15]. They used Carrera’s Unified Formulation to derive the equations of motion for throughthethickness deformations. They showed the meshless technique in several numerical examples.
Hebali et al. [16] presented a new quasithreedimensional (3D) hyperbolic shear deformation theory for the bending and free vibration analysis of FG plates. This theory reduced the number of unknowns by dividing the transverse displacement into bending, shear, and thickness stretching parts. Also their theory satisfied the zero traction boundary conditions on the surfaces of the plate without using shear correction factor.
None of the aforementioned papers have investigated the vibration of porous plate using linear poroelasticity theory besides changing in Young modulus and density. The significance of this paper lies in using the thirdorder shear deformation theory (TSDT) to derive the equation of motions according to thickness ratio of plate for thin and thick plates. Variation of Young modulus and density in transvers direction of plate provides computational complexity that was solved by MAPLE and MATLAB software programs. Three different materials are compared with each other and different boundary conditions are applied on porous plate. The results show the effect of aspect ratio, thickness ratio, elastic medium, and porosity on dimensionless frequency of porous plate accordingly. This paper provides an initial platform for research on free and forced vibration of porous plates according to the many uses for this material. This article will help researchers to continue the studies about forced vibration of porous plate, drain and undrain states, nonconsolidated materials, and different states of fluid in porous materials.
Porous plate can be used in aerospace industry and nuclear reactor as a vibration damper, in sea structures, vessels, and submarines due to very low density and in reformer and catalysts due to the high specific surface. The results of this research can be applied in all cases.
2. Structural Definition
A schematic diagram of a porous plate is illustrated in Figure 1 in which geometrical parameters of length , width , and thickness are shown as well.
3. Linear Poroelasticity Theory
The Biot formulation of the constitutive equations for a fluidfilled porous material assumes the linearity between the stress and the strain and reversibility of the deformation process. Detournay and Chen [17] stated the linear poroelasticity theory of Biot in two statements:(1)An increase of pore pressure induces a dilation of pore.(2)Compression of the pore causes an increase in pore pressure.According to the abovementioned assumptions, stressstrain relation for porous material is shown in [8, 9]in whichwhere and are normal and shear stresses and the terms of engineering constant are as follows [18]:in which and are Young modulus and Poisson’s ratio, respectively. (i) The constant is sometimes called the Biot modulus. The constant α is known as the ratio of the fluid volume gained (or lost) in a material element to the volume change of that element, when the pore pressure is allowed to return to its initial state and it varies between 0 and 1. These parameters are introduced as follows [4, 17]:(ii) Coefficient can be regarded as a coefficient which characterizes the consolidation state of the heterogeneous material. is between 0 for nonconsolidated materials and 1 for consolidated materials [4]; here .(iii) Porosity is defined as ratio of volume of pore space to the volume of porous material. Change in the compressibility of the fluid or porosity affects the bulk constants [17].
The variation of Young modulus and density in transverse direction of porous plate is shown as follows [10, 11]:
In this equation, is the coefficient of plate porosity, and are Young’s modulus of elasticity at and , and and are the density at and , respectively. is the relationship between the modulus of elasticity and shear modulus and is Poisson’s ratio, which is assumed to be constant across the plate thickness. Mechanical properties of the porous material vary across the thickness of the plate , in which
Simplified expressions can be extracted for the poroelastic parameters for limiting cases in particular [17].
For highly compressible fluid constituent , the approximate expressions for are
In the limit, and . In other words, the porous material behaves as an elastic material without fluid where is called porosity and is bulk modulus.
4. ThirdOrder Shear Deformation Theory (TSDT)
TSDT expands the firstorder shear deformation theory (FSDT) by assuming the following:(i)Shear strain and consequently shear stress are not constant through the plate thickness, where TSDT leads to better results for a moderately thick plate.(ii)A shear correction factor is not required for the strain equations.(iii)The displacement field accommodates a quadratic variation of the transverse shear through the thickness and the vanishing of transverse shear is mainly highlighted on the top and bottom surfaces of the plate [18, 19].The displacement field of TSDT for plates is presented here [5]: where are displacement functions along directions and are rotations about  and axes and is time. Using Hooke’s law, the linear strain field for TSDT is being calculated.
5. Energy Method in Porous Plate
In this part, the energy method is used to work out the governing equations. Strain and kinetic energy of the rectangular porous plate is calculated and presented [20, 21]:
Using (1) and (7), the strain energy of porous plate can be derived.
is the density of porous plate and is tractions. By substituting (7) into (9), the kinetic energy is obtained.
5.1. External Work
Pasternak foundation is able to consider transverse shear loads and normal loads. The effect of surrounding elastic medium on the porous plate which is simulated with Pasternak model is considered as follows [21]:where and are Winkler modulus for normal load and shear modulus for transverse shear loads, respectively. Therefore, the external work due to elastic medium is calculated as
5.2. Hamilton’s Principle
Hamilton’s principle is utilized to obtain equations of motion. The principle can be stated in an analytical form in which the first variation form of equations must be zero [20].
In this equation, , and are variations of strain energy, variation of Kinetic energy, and variation of external work. We substituted (8) and (9) and external work into (12) for TSDT and used dimensionless parameters afterwards which were introduced in
The equations of motion are worked out by setting the coefficients equal to zero:in which
6. Solution Procedure Using DQM
The first (16) is used to separate the spatial and temporal distribution:in which is the dimensionless frequency ( is the dimensional frequency). In DQM, the differential equations change into the first algebraic equations. According to what was mentioned, the partial derivations of the function () are estimated by a specific variable at discontinuous points using a set of weighting series. It is supposed that has been a function representing , and regarding variables and () when is the grid points along these variables with the following derivation [22]:
In this equation, are the weighting coefficients which are derived using Chebyshev polynomials for the positions of the grid points whose recursive formulae can be found in [23]. The standard form of the equation of motion is obtained by applying DQM and using (16) into governing (14a), (14b), (14c), (14d), and (14e) considering simply supported boundary conditions as mentioned below:
In this part, the eigenvalue problem is encountered in which the eigenvalues of are introduced as dimensionless frequencies. It is worth mentioning that is the mass matrix and is the stiffness.
In order to ensure the correct results, linear frequency for an isotropic square plate is compared in Table 1 for two modes. The results of the present work and other published papers have been matched well ([24–26]) and it shows that the correct solution has been used in the process of the research.
7. Numerical Results and Discussion
In this study, TSDT was used to work out the equations of motion of rectangular porous plate. The effects of aspect ratio, thickness ratio, boundary condition, porosity, and material properties are included in the results of the study. Detournay and Chen [17] have presented the material properties of saturated porous plate (Tennessee marble) and it is listed in Table 2. Besides, Table 2 reports on the properties of Yfoam, Gfoam, and Coustone.

Figure 2 illustrates the effect of porosity on dimensionless natural frequency of porous plate. It is crystal clear that natural frequency is reduced by increasing the porosity. In fact, according to (5), the porosity affects Young modulus of porous plate.
Based on the thickness ratio, plates are classified, where has been determined for thick plates. Classical theory of plates is applicable to very thin and moderately thin plates, while HSDT applies to thick plates. Figure 2 also shows that increasing leads to decrease in the dimensionless frequency for all values of porosity.
Figure 3 displays the effect of aspect ratio on dimensionless natural frequency of porous plate in different porosities. Aspect ratio is defined as a geometrical parameter of porous plate introducing the ratio of length to width. As it can be seen from Figure 3, the vibration frequency of porous plate does not necessarily follow the determined trend for and the results are redundant in this range. But for larger aspect ratios , the dimensionless natural frequency is increased by increasing .
Figures 4(a) and 4(b) show the impact of five different boundary conditions on frequency response of porous plate including clamped and simple support. Five boundary conditions of CCCC, CCCS, CCSS, CSSS, and SSSS have been compared in this figure in which the value of dimensionless frequency for CCCC is larger than the rest of values. Results show that CCCC boundary condition provides the maximum value for natural frequency and the minimum value is related to SSSS boundary condition. On the one hand, it can be also concluded that, in Figure 4(b) when the aspect ratio of porous plate increases, the difference among above boundary conditions decreases at a higher porosity. On the other hand, the geometrical parameters highly affect vibration behavior of porous plate.
(a)
(b)
Dimensionless frequency of porous plate versus porosity has been shown in Figure 5 for four primary modes. It is obvious from the figure that higher mode is associated with a greater frequency and also none of modes have overlap in .
Figure 6 displays the deflection of porous plate versus porosity for four primary modes in
Figure 7 illustrates the variation of dimensionless natural frequency with porosity of porous plate in different elastic mediums. Figure 7 shows the role of elastic medium in stability of system in which the lowest curve reports the values of dimensionless frequency without elastic medium. In this figure, has the minimum value to the other cases. 2th curve belongs to the Winkler foundation with spring constant