Free Vibration of Embedded Porous Plate Using Third-Order Shear Deformation and Poroelasticity Theories

Copyright © 2017 Ali Ghorbanpour Arani et al.This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. This research aims at studying free vibration of rectangular platemade of porousmaterials in which Y-foam,G-foam, andCoustone 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 comparisonwith classical plate andfirst-order 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.


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 non-Newtonian 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 non-Newtonian 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 fluid-solid motion.Using Galerkin's variational method and classical theory of plates (nonporous), the equations were solved.The Rayleigh-Ritz decomposition on the basis of orthogonal or quasiorthogonal Eigen-functions 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 stress-strain 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 porous-cellular 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 porous-cellular 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 porous-cellular 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 strain-displacement 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 porous-soft 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 power-law 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 sensor-actuator 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 power-law exponents, porosity distributions, porosity volume fractions, and thermal effects on natural frequencies of the temperature-dependent FG beams.
Bending and vibration of FG beams were analyzed with a simple and refined trigonometric higher-order beam theory by Bourada et al. [14].In this theory the thickness stretching effect ( ̸ = 0) 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 throughthe-thickness deformations.They showed the meshless technique in several numerical examples.
Hebali et al. [16] presented a new quasi-three-dimensional (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 third-order 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.

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.

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 above-mentioned assumptions, stress-strain relation for porous material is shown in [8,9] in which where   and   are normal and shear stresses and the terms of engineering constant are as follows [18]: (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  = 1.
(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,  1 is the coefficient of plate porosity,  1 and  0 are Young's modulus of elasticity at  = ℎ/2 and  = −ℎ/2, and  1 and  0 are the density at  = ℎ/2 and  = −ℎ/2, respectively. = /(1 + 2) 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 Simplified expressions can be extracted for the poroelastic parameters for limiting cases in particular [17].
For highly compressible fluid constituent (  / ≪ 1, : fluid), the approximate expressions for  are In the limit,   → 0 and  → 0. In other words, the porous material behaves as an elastic material without fluid where  is called porosity and  is bulk modulus.

Third-Order Shear Deformation Theory (TSDT)
TSDT expands the first-order 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].

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.

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

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].

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 , , ,  1 , and  2 regarding variables  and  (0 <  < 1, 0 <  < 1) 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 ( Ẍ +  = 0) 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 eigen-values 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  = Ω 2 √ ⋅ ℎ/ √ ,  = /(1 −  2 ) 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][25][26]) and it shows that the correct solution has been used in the process of the research.

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 Y-foam, G-foam, 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 0.05 < ℎ/ < 0.3 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  < 4 and the results are redundant in this range.But for larger aspect ratios ( ≥ 4), the dimensionless natural frequency is increased by increasing .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.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 0.1 ≤  1 ≤ 0.5.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   = 10 8 .Pasternak modules create the larger values for natural frequency of porous plate due to normal and shear modules (  = 10 8 ,   = 10 8 ).
According to Figure 8, dimensionless natural frequency of porous plate changes three materials of Y-foam, G-foam, and .Y-foam, G-foam, and  are known as porous materials with different Young modulus and porosities.According to Table 1, comparisons have been made: Young modulus of  is larger than foams and the porosity of foams is higher than .Among foams, Yfoam has a larger Young modulus and a lower porosity, so the most dimensionless natural frequency belongs to  plate and the lowest one is for G-foam.

Conclusion
The significance and novelty of this research lie in studying free vibration of embedded porous which is a relatively new and prosperous topic in this field.The linear poroelasticity theory was used to work out the stress-strain relations in porous material while the Young modulus and density vary in transverse direction of porous plate.To evaluate the effect of normal and shear modulus on the stability of porous plate, Pasternak foundation was developed.TSDT with five unknown functions was used to derive the equation of  frequency of porous plate.On the other hand, porosity reduces the density and Young modulus of matter.Thus, porous materials have the lower natural frequency due to porosity.The porosity in porous materials also creates the damping that it can be used in the aerospace industry.
(ii) Pasternak with normal (  ) and shear (  ) modulus can be regarded as important factors in stability of porous plate in which the dimensionless frequency increases with increasing the stiffness.So the vibration frequency of porous plate can be increased by providing a hard bed.
(iii) Aspect ratio and thickness ratio as geometrical parameters have an efficient and vital role in vibration behavior of porous plate.Therefore, choosing the right dimensions and thickness can be achieved to a certain natural frequency.
(iv) Result shed light on the ground that CCCC boundary condition provides maximum value for natural frequency of porous plate and minimum value relates to SSSS boundary condition.It is evident that CCCC boundary condition increases the stiffness of porous plate.
(v)  as a porous material has larger natural frequency in comparison to Y-foam and G-foam due to high Young modulus and density.
The result of this study can benefit the designers and manufacturers of aircrafts, cars, automobiles, and pneumatic conveying and material handling systems including aerators, airlocks, bins, blenders, blowers, boosters, control systems, air and gas compressors, ship, and barge loaders and unloaders.

Figure 1 :
Figure 1: A schematic diagram of a porous plate.

Figure 2 :
Figure 2: Effect of porosity and thickness ratio on dimensionless natural frequency.

Figure 3 :
Figure 3: Effect of porosity and aspect ratio on dimensionless natural frequency.

Figures 4 (Figure 4 : 5 e 1 (
Figure3displays 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 Figure3, the vibration frequency of porous plate does not necessarily follow the determined trend for  < 4 and the results are redundant in this range.But for larger aspect ratios ( ≥ 4), 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

Figure 5 :
Figure 5: Variation of dimensionless natural frequency versus porosity for four primary modes.

Figure 6
Figure 6  displays the deflection of porous plate versus porosity for four primary modes in  1 = 0.2,  = 0.2,  = 1, SSSS.Figure7illustrates the variation of dimensionless natural frequency with porosity of porous plate in different elastic mediums.Figure7shows 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   = 10 8 .Pasternak modules create the larger values for natural frequency of porous plate due to normal and shear modules (  = 10 8 ,   = 10 8 ).According to Figure8, dimensionless natural frequency of porous plate changes three materials of Y-foam, G-foam, and .Y-foam, G-foam, and  are known as porous materials with different Young modulus and porosities.According to Table1, comparisons have been made: Young modulus of  is larger than foams and the porosity of foams is higher than .Among foams, Yfoam has a larger Young modulus and a lower porosity, so the most dimensionless natural frequency belongs to  plate and the lowest one is for G-foam.

Figure 7 :Figure 8 :
Figure 7: Effect of elastic medium on dimensionless natural frequency.

Table 1 :
Comparison of the results.

Table 2 :
Mechanical properties of porous materials.