Free Vibration Analysis of A Sandwich Cylindrical Shell with A Functionally Graded Auxetic Honeycomb Core via The zig-zag Theory
Subject Areas : Structural Mechanics
1 - Fars province, Abadeh, Islamic Azad University
Keywords: Free vibration, Rotating shell, Auxetic honeycomb structure, Functionally graded material,
Abstract :
In the present paper, a semi-analytical solution is presented for the free vibration analysis of a sandwich cylindrical shell with a re-entrant auxetic honeycomb (AH) core fabricated from metal-ceramic functionally graded materials (FGM). It is assumed that the volume fraction of the ceramic phase in the functionally graded auxetic honeycomb (FGAH) core increases from zero at the inner surface of the core to one at the outer one according to various patterns including power-law function (P-FGM), sigmoid function (S-FGM), and exponential function (E-FGM). The FGAH core is covered with an isotropic homogenous inner face layer made of metal and an isotropic homogenous outer one made of ceramic. The sandwich shell is modeled via Murakami’s zig-zag theory, and the governing equations are derived using Hamilton’s principle. An exact solution is presented for a simply supported shell via the Navier method to find the natural frequencies of the shell. The effects of various parameters on the natural frequencies are studied such as material gradation, the thickness-to-radius ratio, the core-to-face layers thickness ratio, and geometric factors of the auxetic cells. It is found that for each vibrational mode, an optimal ratio can be found between the thickness of the FGAH core and the thickness of the shell which leads to the highest natural frequency.
[1] A. Ghorbanpour Arani, M. Jamali, A.H. Ghorbanpour Arani, R. Kolahchi, M. Mosayyebi, Electro-magneto wave propagation analysis of viscoelastic sandwich nanoplates considering surface effects, Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science 231(2) (2017) 387-403.
[2] E. Haghparast, A.A. Ghorbanpour Arani, A. Ghorbanpour Arani, Effect of fluid–structure interaction on vibration of moving sandwich plate with Balsa wood core and nanocomposite face sheets, International Journal of Applied Mechanics 12(07) (2020) 2050078.
[3] A.H. Ghorbanpour Arani, M. Abdollahian, A. Ghorbanpour Arani, Nonlinear dynamic analysis of temperature-dependent functionally graded magnetostrictive sandwich nanobeams using different beam theories, Journal of the Brazilian Society of Mechanical Sciences and Engineering 42 (2020) 1-20.
[4] A.A. Ghorbanpour Arani, Z. Khoddami Maraghi, A. Ghorbanpour Arani, The Frequency Response of Intelligent Composite Sandwich Plate Under Biaxial In-Plane Forces, Journal of Solid Mechanics 15(1) (2023).
[5] A.H. Ghorbanpour Arani, A. Rastgoo, M. Sharafi, R. Kolahchi, A. Ghorbanpour Arani, Nonlocal viscoelasticity based vibration of double viscoelastic piezoelectric nanobeam systems, Meccanica 51 (2016) 25-40.
[6] A.H. Ghorbanpour Arani, A. Rastgoo, A. Hafizi Bidgoli, R. Kolahchi, A. Ghorbanpour Arani, Wave propagation of coupled double-DWBNNTs conveying fluid-systems using different nonlocal surface piezoelasticity theories, Mechanics of Advanced Materials and Structures 24(14) (2017) 1159-1179.
[7] P. Sourani, A. Ghorbanpour Arani, M. Hashemian, S. Niknejad, Nonlinear dynamic stability analysis of CNTs reinforced piezoelectric viscoelastic composite nano/micro plate under multiple physical fields resting on smart foundation, Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science 238(10) (2024) 4307-4342.
[8] A. Ghorbanpour Arani, N. Miralaei, A. Farazin, M. Mohammadimehr, An extensive review of the repair behavior of smart self-healing polymer matrix composites, Journal of Materials Research 38(3) (2023) 617-632.
[9] K. Torabi, H. Afshari, F. Hajiaboutalebi, Vibration and flutter analyses of cantilever trapezoidal honeycomb sandwich plates, Journal of Sandwich Structures & Materials 21(8) (2019) 2887-2920.
[10] S. Galehdari, H. Khodarahmi, A. Atrian, Design and analysis of graded honeycomb shock absorber for increasing the safety of passengers in armored vehicles exposed to mine explosion, Journal of Solid Mechanics 9(2) (2017) 370-383.
[11] H. Alikhani, S. Derakhshan, H. Khoramishad, Experimental and Numerical Investigations of the Effect of Impact Angle and Impactor Geometry on the High-Velocity Impact Response of Aluminum Honeycomb Structures, Journal of Solid Mechanics 14(4) (2022).
[12] T.D. Qing, C.Y. Zhi, Wave propagation in sandwich panel with auxetic core, Journal of Solid Mechanics 2(4) (2010) 393-402.
[13] N.D. Duc, K. Seung-Eock, P.H. Cong, N.T. Anh, N.D. Khoa, Dynamic response and vibration of composite double curved shallow shells with negative Poisson's ratio in auxetic honeycombs core layer on elastic foundations subjected to blast and damping loads, International Journal of Mechanical Sciences 133 (2017) 504-512.
[14] H. Eipakchi, F.M. Nasrekani, Vibrational behavior of composite cylindrical shells with auxetic honeycombs core layer subjected to a moving pressure, Composite Structures 254 (2020) 112847.
[15] N.V. Nguyen, H. Nguyen-Xuan, T.N. Nguyen, J. Kang, J. Lee, A comprehensive analysis of auxetic honeycomb sandwich plates with graphene nanoplatelets reinforcement, Composite Structures 259 (2021) 113213.
[16] P. Xiao, Z. Yifeng, S. Jie, S. Zheng, Global buckling analysis of composite honeycomb sandwich plate with negative Poisson’s ratio (CHSP-NPR) using variational asymptotic equivalent model, Composite Structures 264 (2021) 113721.
[17] F. Xu, K. Yu, L. Hua, In-plane dynamic response and multi-objective optimization of negative Poisson's ratio (NPR) honeycomb structures with sinusoidal curve, Composite Structures 269 (2021) 114018.
[18] Q.-H. Pham, P.-C. Nguyen, T.T. Tran, Free vibration response of auxetic honeycomb sandwich plates using an improved higher-order ES-MITC3 element and artificial neural network, Thin-Walled Structures 175 (2022) 109203.
[19] T.Q. Quan, V.M. Anh, V. Mahesh, N.D. Duc, Vibration and nonlinear dynamic response of imperfect sandwich piezoelectric auxetic plate, Mechanics of Advanced Materials and Structures 29(1) (2022) 127-137.
[20] N.D. Dat, T.Q. Quan, N.D. Duc, Vibration analysis of auxetic laminated plate with magneto-electro-elastic face sheets subjected to blast loading, Composite Structures 280 (2022) 114925.
[21] Y. Li, B. Liu, Thermal buckling and free vibration of viscoelastic functionally graded sandwich shells with tunable auxetic honeycomb core, Applied Mathematical Modelling 108 (2022) 685-700.
[22] A. Pakrooyan, P. Yousefi, K. Khorshidi, M.M. Najafizadeh, A. Nezamabadi, Free vibration analysis of an auxetic honeycomb sandwich plate placed at the wall of a fluid tank, Ocean Engineering 263 (2022) 112353.
[23] P.H. Cong, V.D. Trung, N.D. Khoa, N.D. Duc, Vibration and nonlinear dynamic response of temperature-dependent FG-CNTRC laminated double curved shallow shell with positive and negative Poisson’s ratio, Thin-Walled Structures 171 (2022) 108713.
[24] B. Liu, S. Li, Y. Li, Bending of FGM sandwich plates with tunable auxetic core using DQM, European Journal of Mechanics-A/Solids 97 (2023) 104838.
[25] S. Sadikbasha, V. Pandurangan, Crashworthiness of auxetic hexachiral core under in-plane loading, Materials Today: Proceedings (2023).
[26] B. Nečemer, Ž. Božić, S. Glodež, Fatigue resistance of the auxetic honeycombs, Procedia Structural Integrity 46 (2023) 68-73.
[27] M. Rai, A. Chawla, S. Mukherjee, Parametric study of re-entrant honeycomb cored auxetic sandwich panel exposed to blast loading, Materials Today: Proceedings (2023).
[28] M. Sarafraz, H. Seidi, F. Kakavand, N.S. Viliani, Free vibration and buckling analyses of a rectangular sandwich plate with an auxetic honeycomb core and laminated three-phase polymer/GNP/fiber face sheets, Thin-Walled Structures 183 (2023) 110331.
[29] M. Sarafraz, H. Seidi, F. Kakavand, N.S. Viliani, Flutter characteristics of a rectangular sandwich plate with laminated three-phase polymer/GNP/fiber face sheets and an auxetic honeycomb core in yawed supersonic fluid flow, Journal of the Brazilian Society of Mechanical Sciences and Engineering 45(4) (2023) 197.
[30] M. Koizumi, The concept of FGM, Ceramic transactions 34 (1993) 3-10.
[31] J.-h. Zhang, B.-j. Dong, B. He, Y. Sun, Free vibrations and impact resistance of a functionally graded honeycomb sandwich plate, Shock and Vibration 2021 (2021) 1-15.
[32] F. Li, W. Yuan, C. Zhang, Free vibration and sound insulation of functionally graded honeycomb sandwich plates, Journal of Sandwich Structures & Materials 24(1) (2022) 565-600.
[33] H. Amirabadi, A. Darakhsh, M. Sarafraz, H. Afshari, Vibration characteristics of a cylindrical sandwich shell with an FG auxetic honeycomb core and nanocomposite face layers subjected to supersonic fluid flow, Noise & Vibration Worldwide (2024) 09574565241252994.
[34] H. Amirabadi, A. Mottaghi, M. Sarafraz, H. Afshari, Free vibrational behavior of a conical sandwich shell with a functionally graded auxetic honeycomb core, Journal of Vibration and Control (2024) 10775463241240215.
[35] X. Zhu, J. Zhang, W. Zhang, J. Chen, Vibration frequencies and energies of an auxetic honeycomb sandwich plate, Mechanics of Advanced Materials and Structures 26(23) (2019) 1951-1957.
[36] H. Murakami, Laminated composite plate theory with improved in-plane responses, (1986).
[37] A. Ghorbanpour Arani, F. Kiani, H. Afshari, Free and forced vibration analysis of laminated functionally graded CNT-reinforced composite cylindrical panels, Journal of Sandwich Structures & Materials 23(1) (2021) 255-278.
[38] A. Ghorbanpour Arani, F. Kiani, H. Afshari, Aeroelastic analysis of laminated FG-CNTRC cylindrical panels under yawed supersonic flow, International Journal of Applied Mechanics 11(06) (2019) 1950052.
[39] J.N. Reddy, Mechanics of laminated composite plates and shells: theory and analysis, CRC press 2003.
[40] V. Mohammadlou, Z. Khoddami Maraghi, A. Ghorbanpour Arani, Thermoelastic analysis of axisymmetric conical shells: Investigating stress–strain response under uniform heat flow with semi-coupled approach, Numerical Heat Transfer, Part A: Applications (2024) 1-22.
[41] J.N. Reddy, Energy principles and variational methods in applied mechanics, John Wiley & Sons 2017.
[42] A. Ghorbanpour Arani, E. Haghparast, A.H. Ghorbanpour Arani, Size‐dependent vibration of double‐bonded carbon nanotube‐reinforced composite microtubes conveying fluid under longitudinal magnetic field, Polymer Composites 37(5) (2016) 1375-1383.
[43] H. Li, K.-Y. Lam, T.-Y. Ng, Rotating shell dynamics, Elsevier 2005.
[44] C. Loy, K. Lam, J. Reddy, Vibration of functionally graded cylindrical shells, International Journal of Mechanical Sciences 41(3) (1999) 309-324.
Journal of Solid Mechanics Vol. 16, No. 1 (2024) pp. 48-64 DOI: 10.60664/jsm.2024.3111794 |
Research Paper Free Vibration Analysis of A Sandwich Cylindrical Shell with A Functionally Graded Auxetic Honeycomb Core via The zig-zag Theory |
H. Amirabadi 1 | |
Department of Mechanical Engineering, Abadeh Branch, Islamic Azad University, Abadeh, Iran | |
Received 16 November 2023; Received in revised form 29 June 2024; Accepted 2 July 2024 | |
| ABSTRACT |
| In the present paper, a semi-analytical solution is presented for the free vibration analysis of a sandwich cylindrical shell with a re-entrant auxetic honeycomb (AH) core fabricated from metal-ceramic functionally graded materials (FGM). It is assumed that the volume fraction of the ceramic phase in the functionally graded auxetic honeycomb (FGAH) core increases from zero at the inner surface of the core to one at the outer one according to various patterns including power-law function (P-FGM), sigmoid function (S-FGM), and exponential function (E-FGM). The FGAH core is covered with an isotropic homogenous inner face layer made of metal and an isotropic homogenous outer one made of ceramic. The sandwich shell is modeled via Murakami’s zig-zag theory, and the governing equations are derived using Hamilton’s principle. An exact solution is presented for a simply supported shell via the Navier method to find the natural frequencies of the shell. The effects of various parameters on the natural frequencies are studied such as material gradation, the thickness-to-radius ratio, the core-to-face layers thickness ratio, and geometric factors of the auxetic cells. It is found that for each vibrational mode, an optimal ratio can be found between the thickness of the FGAH core and the thickness of the shell which leads to the highest natural frequency.
|
| Keywords: Free vibration; Rotating shell; Auxetic honeycomb structure; Functionally graded material. |
1 INTRODUCTION
T
HE main aims of utilizing sandwich structures rather than single-layer ones are to achieve high stiffness-to-mass and strength-to-mass ratios [1, 2] or to attain a smart structure whose mechanical properties can be affected by an external factor [3-8]. Structures that benefit from low density and tunable elastic constants can be selected as good choices to be used as a core in a sandwich structure. Since the core is not subjected to intense bending loads, it does not necessarily need to be of high rigidity and strength. Honeycombs are fragile but low-density structures that can be utilized as a core to attain a sandwich structure with an excellent value of stiffness-to-weight ratio. Thus, honeycomb structures have been used in various engineering fields such as transportation, civil, mechanical, and aerospace engineering. A fair amount of work has been presented in recent years to investigate the mechanical behavior of three-layered sandwich structures with honeycomb cores. As the most well-known type of honeycomb structure, the hexagonal one has been extensively used in various engineering fields [9-11]. In a hexagonal honeycomb structure, the Poisson’s ratios are always positive values. However, other kinds of honeycomb structures have been proposed in recent years whose Poisson’s ratios are negative values which are called auxetic honeycomb (AH) structures.
A fair amount of work has been presented in recent years associated with the mechanical analysis of sandwich structures with AH cores. Qing Tian and Zhi Chun [12] examined the wave propagation characteristics of a sandwich rectangular plate with an AH core and isotropic homogenous face layers. They examined the effects of the geometric parameters of the AH core on the wave propagation characteristics of the plate. The free and forced vibration analyses of a doubly-curved sandwich panel with an AH core subjected to blast load were studied by Duc et al. [13]. They focused on the impacts of the geometric parameters of the AH core on the dynamic deflection and natural frequencies of the shell. The dynamic response of a sandwich cylindrical shell with an AH core subjected to moving internal pressure was investigated by Eipakchi and Naserkani [14]. The effects of the geometric parameters of the AH core on the critical values of the frequency and velocity of the moving pressure were examined by them. The dynamic buckling and free vibration analyses of sandwich rectangular plates with an AH core and polymeric face sheets reinforced with graphene nanoplatelets (GNP) were studied by Nguyen et al. [15]. They studied the dependency of the natural frequencies and the stability regions on the geometric factors of the AH core. Xiao et al. [16] inspected the mechanical buckling analysis of a sandwich rectangular plate with an AH core and two laminated composite face sheets. They examined the effects of geometric factors of the AH core on the critical buckling load of the plate. Xu et al. [17] examined the free in-plane vibration behavior of an AH structure with curved sinusoidal walls. They focused on the optimization of such a structure and inspected the impacts of utilizing it on the energy absorption capacity of the structure and the amplitude of its oscillations. The free vibration analysis of a sandwich rectangular plate with an AH core and two FGM face sheets was studied by Pham et al. [18]. They focused on the effects of the geometric factors of the AH core on the natural frequencies of the plate. The nonlinear free and forced vibration analyses of an imperfect sandwich rectangular plate with an AH core and two piezoelectric face sheets were examined by Quan et al. [19]. They studied the dependency of the natural frequencies of the plate and the dynamic deflection on the geometric factors of the AH core.
Dat et al. [20] examined the nonlinear free and forced vibration analyses of a sandwich rectangular plate with an AH core and magneto-electro-elastic face sheets exposed to an explosive load. They inspected the influences of geometric factors of the AH core on the natural frequencies of the plate and the dynamic deflection. The thermal buckling and free vibration analyses of a viscoelastic doubly-curved sandwich shell with a tunable AH core and FGM face layers were examined by Li and Liu [21]. They focused on the effects of the geometric parameters of the AH core on the critical temperature and the natural frequencies of the shell. As a model of the wall of a fluid-filled tank, Pakrooyan et al. [22] investigated the free vibration analysis of a sandwich rectangular plate with an AH core and isotropic homogenous face sheets in contact with quiescent fluid. They studied the dependency of the natural frequencies on the geometric parameters of the AH core. Cong et al. [23] presented a parametric investigation to study the nonlinear free and forced vibration analyses of a doubly-curved sandwich panel with an AH core and two laminated polymeric face layers reinforced with carbon nanotubes (CNTs). They considered the temperature-dependency of the thermo-mechanical properties and examined the effect of geometric factors of the AH core on the natural frequencies and dynamic response of the shell. Liu et al. [24] studied the static bending analysis of a sandwich rectangular plate with a tunable AH core and FGM face sheets. The dependencies of the static deformation and stress distribution on the geometric parameters of the AH core were studied by them. The crashworthiness of hexachiral AH structures exposed to an in-plane loading was studied by Sadikbasha and Pandurangan [25]. They discussed the dependency of the energy absorption capacity of sandwich structures with hexachiral AH core on the geometric parameters of the cells. Necemer et al. [26] used the ANSYS software and analyzed the fatigue resistance of several AH structures. They compared five types of AH structures including re-entrant, S-shaped, star-shaped, chiral, and double arrowhead to check which one has the best fatigue resistance. Rai et al. [27] studied the dynamic response of a sandwich panel with an AH core subjected to explosive loading. They tried to optimize the energy absorption capacity of such a structure. Sarafraz et al. [28, 29] presented numerical solutions to analyze the mechanical buckling, the free vibration, and the flutter (aeroelastic stability) analyses of a sandwich rectangular plate with an AH core and laminated three-phase polymer/GNP/fiber face sheets. They investigated the effects of the geometrical parameters of the AH core and the mass fraction of the fibers and GNPs on the critical buckling load, the natural frequencies, and the critical aerodynamic pressure of the plate.
The concept of FGM was introduced in the 1980s by some materials scientists in Japan [30]. This type of non-homogenous material is usually produced by the composition of two different materials including a ceramic phase and a metal phase whose volume fractions vary continuously in one, two, or even three directions. The main advantage of FGMs is their high resistance against simultaneous thermal and mechanical loads. Therefore, these materials have been used in various engineering fields such as civil, aerospace, and mechanical engineering. As a new idea to achieve low mass and high endurance against thermo-mechanical loading, recently, some researchers have proposed to produce AH structures from FGMs. To the best knowledge of the authors, there are few numbers of papers associated with the influences of FGAH cores on the mechanical behavior of sandwich plates and shells [31-34]. For the first time, an analytical solution is presented on the free vibration analysis of a sandwich cylindrical shell with a re-entrant FGAH core and two homogenous face layers. The sandwich shell is modeled based on the zig-zag theory. Three patterns are considered to describe the variations of the volume fractions of the metal and ceramic including power-law, sigmoid, and exponential functions. The effects of various factors on the natural frequencies of the shell are investigated including material gradation, the thickness-to-radius ratio, the core-to-face layers thickness ratio, and geometric factors of the auxetic cells. The results of this work can be useful in the design, analysis, and optimization of future aerospace structures.
2 MATHEMATICAL MODELING
2.1. Material Properties
As shown in Fig. 1, a three-layered sandwich cylindrical shell of mean radius R, length L, and total thickness h is considered. The sandwich shell consists of a re-entrant FGAH core of thickness hc and isotropic homogeneous face layers of the same thickness h1=h3=hf=0.5(h-hc). The inner and outer layers of the shell are fabricated from metal and ceramic, respectively.
The material properties of the non-homogeneous material utilized to fabricate the AH core vary from a metal-rich surface at the inner side of the core to a ceramic-rich one at the outer surface of the core according to a power-law function as [33]
| (1) |
where P stands for either the density (ρs), the elastic modulus (Es), or the Poisson’s ratio (νs), and the subscripts m and c indicate the material properties of the metal and ceramic phases, respectively. In Eq. (1), the dimensionless parameter p is called the material index which determines how the volume fractions of metal and ceramic vary in the thickness direction.
Fig. 1
A sandwich cylindrical shell with an FGAH core and isotropic homogeneous face layers.
The variations of material properties of the non-homogeneous material utilized to fabricate the AH core through the thickness direction are depicted in Fig. 2 for Pc/Pm=2 and several various values of the material index.
Since an FGM is an isotropic material, the shear modulus (Gs) of the material used to produce the AH core can be calculated through the relation below:
| (2) |
Referring to the geometric parameters of an individual cell of a re-entrant FGAH structure illustrated in Fig. 1, the elastic constants and the density of the FGAH core can be attained using the following relation [35]:
| (3) |
|
Fig. 2 Variations of material properties of the material utilized to fabricate the AH core through the thickness direction. |
where
| (4) |
Eq. (3) indicates that the FGAH core is an orthotropic non-homogenous structure fabricated from an isotropic non-homogenous material.
The density and the elastic constants of the isotropic homogeneous inner and outer face layers can be presented as follows:
| (5) |
where Gm and Gc are the shear moduli of the metal and ceramic phases, respectively.
2.2. Deformation, Strain, and Stress
Based on Murakami’s zig-zag theory, the deformation field can be described as follows [36]:
| (6) |
where u1, u2, and u3 show displacement along x, θ, and z directions, respectively, and φx and φθ are the rotations about θ and x axes, respectively. The variables sx and sθ are used to incorporate the zig-zag effect, and f(z) is defined as follows [36]:
| (7) |
where hk represents the thickness of the kth layer, zk and zk+1 are the transverse position of the bottom and top of the kth face layer, respectively, NL stands for the number of layers which is 3 in this paper, and H is the well-known Heaviside function.
It is noteworthy that by removing sx, sθ, and f(z) the well-known first-order shear deformation theory (FSDT) is attained. A schematical comparison between the displacement fields in the zig-zag and the FSDT theories is depicted in Fig. 3.
The normal (εij) and shear (γij=2εij) components of the strain tensor can be calculated through the following relations [37, 38]:
| (8) |
where
| (9) |
It should be noted that Eq. (7) is attained by considering the assumption of shallow shells:
| (10) |
The non-zero components of the stress tensor (σij) in the kth layer of the shell can be calculated as [39, 40]
| (11) |
where ks=5/6 is the well-known shear correction factor and
| (12) |
Fig. 3
The deformation in the zig-zag and the FSDT theories.
2.3. Hamilton’s Principle
According to Hamilton’s principle, the governing equations and boundary conditions regarding the vibration analysis of a structure can be obtained using the relation below [41, 42]:
| (13) |
in which δ is the well-known variational operator, t shows time, [t1,t2] represents an arbitrary time interval, and. U, T, and Wn.c. indicate the strain energy of the shell, the kinetic energy of the shell, and the work done by non-conservative external loads, respectively.
The variation of the strain energy of the shell can be described as follows [41]:
| (14) |
in which V shows the volume of the shell and is defined as
| (15) |
where dS=Rdxdθ is the surface of the shell at the middle surface (z=0).
Utilizing Eqs. (8), (14), and (15), the following relation can be presented:
| (16) |
where
| (17) |
By substituting Eqs. (8) and (11) into Eq. (17), one can write the relation below:
(18) |
|
in which the stiffness coefficients are defined as follows:
| (19) |
The equation below describes the variation of the kinetic energy of the shell [41]:
| (20) |
Utilizing Eqs. (6) and (15), Eq. (20) can be represented as follows:
| (21) |
where
| (22) |
Substituting Eqs. (16) and (21) into Eq. (13) and considering Wn.c.=0 for the free vibration analysis, the following relations can be attained as the governing equations:
| (23) |
Also, the following equations describe the boundary conditions for simply supported edges at x=0 & L:
|
| ||
[1] Corresponding author. Tel.: +987144351093.
E-mail address: h.amirabadi@iauabadeh.ac.ir (Hossein Amirabadi)