Effect of Columns Stiffness's and Number of Floors on the Accuracy of the Tributary Area Method
The using of finite element programs in analyzing and designing buildings are becoming very popular, but there are many engineers still using the tributary area method (TAM) in designing the structural members such as columns. This study is an attempt to investigate the accuracy of the TAM results with different load condition (gravity and lateral load), different floors numbers, and different columns stiffness's. To conduct this study, linear elastic analysis in ETABS program is used. The results from finite element method are compared to those obtained from TAM. According to the analysis of the data obtained, it can be seen that there is significance difference between the real load carried by columns and the load which is calculated by using the TAM. Thus, using 3-D models are the best choice to calculate the real load effected on columns and design these columns according to this load.
Rotor Dynamic Analysis for a Shaft Train by Using Finite Element Method
In the present paper, a large turbo-generator shaft train including a heavy-duty gas turbine engine, a coupling, and a generator is established. The method of analysis is based on finite element simplified model for lateral and torsional vibration calculation. The basic elements of rotor are the shafts and the disks which are represented as circular cross section flexible beams and rigid body elements, respectively. For more accurate results, the gyroscopic effect and bearing dynamics coefficients and function of rotation are taken into account, and for the influence of shear effect, rotor has been modeled in the form of Timoshenko beam. Lateral critical speeds, critical speed map, damped mode shapes, Campbell diagram, zones of instability, amplitudes, phase angles response due to synchronous forces of excitation and amplification factor are calculated. Also, in the present paper, the effect of imbalanced rotor and effects of changing in internal force and temperature are studied.
Active Control Improvement of Smart Cantilever Beam by Piezoelectric Materials and On-Line Differential Artificial Neural Networks
The main goal of this study is to test differential
neural network as a controller of smart structure and is to enumerate
its advantages and disadvantages in comparison with other
controllers. In this study, the smart structure has been considered as a
Euler Bernoulli cantilever beam and it has been tried that it be under
control with the use of vibration neural network resulting from
movement. Also, a linear observer has been considered as a reference
controller and has been compared its results. The considered
vibration charts and the controlled state have been recounted in the
final part of this text. The obtained result show that neural observer
has better performance in comparison to the implemented linear
Static Modeling of the Delamination of a Composite Material Laminate in Mode II
The purpose of this paper is to analyze numerically by the three-dimensional finite element method, using ABAQUS calculation code, the mechanical behavior of a unidirectional and multidirectional delaminated stratified composite under mechanical loading in Mode II. This study consists of the determination of the energy release rate G in mode II as well as the distribution of equivalent von Mises stresses along the damaged zone by varying several parameters such as the applied load and the delamination length. It allowed us to deduce that the high energy release rate favors delamination at the free edges of a stratified plate subjected to bending.
A Numerical Investigation of Lamb Wave Damage Diagnosis for Composite Delamination Using Instantaneous Phase
This paper presents a study of Lamb wave damage
diagnosis of composite delamination using instantaneous phase
data. Numerical experiments are performed using the finite element
method. Different sizes of delamination damages are modeled
using finite element package ABAQUS. Lamb wave excitation
and responses data are obtained using a pitch-catch configuration.
Empirical mode decomposition is employed to extract the intrinsic
mode functions (IMF). Hilbert–Huang Transform is applied to each
of the resulting IMFs to obtain the instantaneous phase information.
The baseline data for healthy plates are also generated using the
same procedure. The size of delamination is correlated with the
instantaneous phase change for damage diagnosis. It is observed that
the unwrapped instantaneous phase of shows a consistent behavior
with the increasing delamination size.
Structural Evaluation of Airfield Pavement Using Finite Element Analysis Based Methodology
Nondestructive deflection testing has been accepted widely as a cost-effective tool for evaluating the structural condition of airfield pavements. Backcalculation of pavement layer moduli can be used to characterize the pavement existing condition in order to compute the load bearing capacity of pavement. This paper presents an improved best-fit backcalculation methodology based on deflection predictions obtained using finite element method (FEM). The best-fit approach is based on minimizing the squared error between falling weight deflectometer (FWD) measured deflections and FEM predicted deflections. Then, concrete elastic modulus and modulus of subgrade reaction were back-calculated using Heavy Weight Deflectometer (HWD) deflections collected at the National Airport Pavement Testing Facility (NAPTF) test site. It is an alternative and more versatile method in considering concrete slab geometry and HWD testing locations compared to methods currently available.
Modeling and Simulation of Ship Structures Using Finite Element Method
The development in the construction of unconventional ships and the implementation of lightweight materials have shown a large impulse towards finite element (FE) method, making it a general tool for ship design. This paper briefly presents the modeling and analysis techniques of ship structures using FE method for complex boundary conditions which are difficult to analyze by existing Ship Classification Societies rules. During operation, all ships experience complex loading conditions. These loads are general categories into thermal loads, linear static, dynamic and non-linear loads. General strength of the ship structure is analyzed using static FE analysis. FE method is also suitable to consider the local loads generated by ballast tanks and cargo in addition to hydrostatic and hydrodynamic loads. Vibration analysis of a ship structure and its components can be performed using FE method which helps in obtaining the dynamic stability of the ship. FE method has developed better techniques for calculation of natural frequencies and different mode shapes of ship structure to avoid resonance both globally and locally. There is a lot of development towards the ideal design in ship industry over the past few years for solving complex engineering problems by employing the data stored in the FE model. This paper provides an overview of ship modeling methodology for FE analysis and its general application. Historical background, the basic concept of FE, advantages, and disadvantages of FE analysis are also reported along with examples related to hull strength and structural components.
Comparison of Electrical Parameters of Oil-Immersed and Dry-Type Transformer Using Finite Element Method
The choice evaluation between oil-immersed and dry-type transformers is often controlled by cost, location, and application. This paper compares the electrical performance of liquid- filled and dry-type transformers, which will assist the customer to choose the right and efficient ones for particular applications. An accurate assessment of the time-average flux density, electric field intensity and voltage distribution in an oil-insulated and a dry-type transformer have been computed and investigated. The detailed transformer modeling and analysis has been carried out to determine electrical parameter distributions. The models of oil-immersed and dry-type transformers are developed and solved by using the finite element method (FEM) to compare the electrical parameters. The effects of non-uniform and non-coherent voltage gradient, flux density and electric field distribution on the power losses and insulation properties of transformers are studied in detail. The results show that, for the same voltage and kilo-volt-ampere (kVA) rating, oil-immersed transformers have better insulation properties and less hysteresis losses than the dry-type.
Study on Sharp V-Notch Problem under Dynamic Loading Condition Using Symplectic Analytical Singular Element
V-notch problem under dynamic loading condition is considered in this paper. In the time domain, the precise time domain expanding algorithm is employed, in which a self-adaptive technique is carried out to improve computing accuracy. By expanding variables in each time interval, the recursive finite element formulas are derived. In the space domain, a Symplectic Analytical Singular Element (SASE) for V-notch problem is constructed addressing the stress singularity of the notch tip. Combining with the conventional finite elements, the proposed SASE can be used to solve the dynamic stress intensity factors (DSIFs) in a simple way. Numerical results show that the proposed SASE for V-notch problem subjected to dynamic loading condition is effective and efficient.
Comparison of ANN and Finite Element Model for the Prediction of Ultimate Load of Thin-Walled Steel Perforated Sections in Compression
The analysis of perforated steel members is a 3D problem in nature, therefore the traditional analytical expressions for the ultimate load of thin-walled steel sections cannot be used for the perforated steel member design. In this study, finite element method (FEM) and artificial neural network (ANN) were used to simulate the process of stub column tests based on specific codes. Results show that compared with those of the FEM model, the ultimate load predictions obtained from ANN technique were much closer to those obtained from the physical experiments. The ANN model for the solving the hard problem of complex steel perforated sections is very promising.
FEM Simulations to Study the Effects of Laser Power and Scan Speed on Molten Pool Size in Additive Manufacturing
Additive manufacturing (AM) is increasingly crucial in biomedical and aerospace industries. As a recently developed AM technique, selective laser melting (SLM) has become a commercial method for various manufacturing processes. However, the molten pool configuration during SLM of metal powders is a decisive issue for the product quality. It is very important to investigate the heat transfer characteristics during the laser heating process. In this work, the finite element method (FEM) software ANSYS® (work bench module 16.0) was used to predict the unsteady temperature distribution for resolving molten pool dimensions with consideration of temperature-dependent thermal physical properties of TiAl6V4 at different laser powers and scanning speeds. The simulated results of the temperature distributions illustrated that the ratio of laser power to scanning speed can greatly influence the size of molten pool of titanium alloy powder for SLM development.
Forced Vibration of a Planar Curved Beam on Pasternak Foundation
The objective of this study is to investigate the forced vibration analysis of a planar curved beam lying on elastic foundation by using the mixed finite element method. The finite element formulation is based on the Timoshenko beam theory. In order to solve the problems in frequency domain, the element matrices of two nodded curvilinear elements are transformed into Laplace space. The results are transformed back to the time domain by the well-known numerical Modified Durbin’s transformation algorithm. First, the presented finite element formulation is verified through the forced vibration analysis of a planar curved Timoshenko beam resting on Winkler foundation and the finite element results are compared with the results available in the literature. Then, the forced vibration analysis of a planar curved beam resting on Winkler-Pasternak foundation is conducted.
Assessment of Material Type, Diameter, Orientation and Closeness of Fibers in Vulcanized Reinforced Rubbers
In this work, the effect of material type, diameter, orientation and closeness of fibers on the general performance of reinforced vulcanized rubbers are investigated using finite element method with experimental verification. Various fiber materials such as hemp, nylon, polyester are used for different fiber diameters, orientations and closeness. 3D finite element models are developed by considering bonded contact elements between fiber and rubber sheet interfaces. The fibers are assumed as linear elastic, while vulcanized rubber is considered as hyper-elastic. After an experimental verification of finite element results, the developed models are analyzed under prescribed displacement that causes tension. The normal stresses in fibers and shear stresses between fibers and rubber sheet are investigated in all models. Large deformation of reinforced rubber sheet also represented with various fiber conditions under incremental loading. A general assessment is achieved about best fiber properties of reinforced rubber sheets for tension-load conditions.
Finite Element Modeling of Heat and Moisture Transfer in Porous Material
This paper presents a two-dimensional model to study the heat and moisture transfer through porous building materials. Dynamic and static coupled models of heat and moisture transfer in porous material under low temperature are presented and the coupled models together with variable initial and boundary conditions have been considered in an analytical way and using the finite element method. The resulting coupled model is converted to two nonlinear partial differential equations, which is then numerically solved by an implicit iterative scheme. The numerical results of temperature and moisture potential changes are compared with the experimental measurements available in the literature. Predicted results demonstrate validation of the theoretical model and effectiveness of the developed numerical algorithms. It is expected to provide useful information for the porous building material design based on heat and moisture transfer model.
A Canadian Leaf Shaped Triple Band Patch Antenna with DGS for X and C-Band Applications
A shaped single feed microstrip antenna is realized for C-Band and X-Band applications. The frequency range of C-band and X-band varies from 4 to 8 Gigahertz and 8 to 12 Gigahertz. The antenna operates under three frequency bands, one under C band and two under X-band applications. Defect on the ground called DGS (Defected Ground Structure) is made to enhance the distinctiveness of the antenna parameters. The design consists of DGS provided to improve the antenna performance. The substrate material used is of the Flame Retardant grade-4 (FR4) epoxy having high mechanical and electrical strength. The design and analysis was done using the FEM (Finite Element Method) based Ansoft HFSS (High Frequency Structural Simulator) Version 12. For the resonant frequencies of 5.21, 9.17 and 10.45, a value of reflection coefficient obtained is of -39.0, -16.0 and -30.7 dB respectively. Other constraints of antenna such as bandwidth, gain, directivity and Voltage Standing Wave Ratio (VSWR) are also conferred.
Thermal Fracture Analysis of Fibrous Composites with Variable Fiber Spacing Using Jk-Integral
In this study, fracture analysis of a fibrous composite
laminate with variable fiber spacing is carried out using Jk-integral
method. The laminate is assumed to be under thermal loading.
Jk-integral is formulated by using the constitutive relations of plane
orthotropic thermoelasticity. Developed domain independent form
of the Jk-integral is then integrated into the general purpose finite
element analysis software ANSYS. Numerical results are generated
so as to assess the influence of variable fiber spacing on mode I
and II stress intensity factors, energy release rate, and T-stress. For
verification, some of the results are compared to those obtained
using displacement correlation technique (DCT).
Electromagnetic Assessment of Submarine Power Cable Degradation Using Finite Element Method and Sensitivity Analysis
Submarine power cables used for offshore wind farms
electric energy distribution and transmission are subject to numerous
threats. Some of the risks are associated with transport, installation
and operating in harsh marine environment. This paper describes the
feasibility of an electromagnetic low frequency sensing technique for
submarine power cable failure prediction. The impact of a structural
damage shape and material variability on the induced electric field is
evaluated. The analysis is performed by modeling the cable using the
finite element method, we use sensitivity analysis in order to identify
the main damage characteristics affecting electric field variation.
Lastly, we discuss the results obtained.
Three Dimensional Finite Element Analysis of Functionally Graded Radiation Shielding Nanoengineered Sandwich Composites
In recent years, nanotechnology has played an important role in the design of an efficient radiation shielding polymeric composites. It is well known that, high loading of nanomaterials with radiation absorption properties can enhance the radiation attenuation efficiency of shielding structures. However, due to difficulties in dispersion of nanomaterials into polymer matrices, there has been a limitation in higher loading percentages of nanoparticles in the polymer matrix. Therefore, the objective of the present work is to provide a methodology to fabricate and then to characterize the functionally graded radiation shielding structures, which can provide an efficient radiation absorption property along with good structural integrity. Sandwich structures composed of Ultra High Molecular Weight Polyethylene (UHMWPE) fabric as face sheets and functionally graded epoxy nanocomposite as core material were fabricated. A method to fabricate a functionally graded core panel with controllable gradient dispersion of nanoparticles is discussed. In order to optimize the design of functionally graded sandwich composites and to analyze the stress distribution throughout the sandwich composite thickness, a finite element method was used. The sandwich panels were discretized using 3-Dimensional 8 nodded brick elements. Classical laminate analysis in conjunction with simplified micromechanics equations were used to obtain the properties of the face sheets. The presented finite element model would provide insight into deformation and damage mechanics of the functionally graded sandwich composites from the structural point of view.
High Cycle Fatigue Analysis of a Lower Hopper Knuckle Connection of a Large Bulk Carrier under Dynamic Loading
The fatigue of ship structural details is of major concern in the maritime industry as it can generate fracture issues that may compromise structural integrity. In the present study, a fatigue analysis of the lower hopper knuckle connection of a bulk carrier was conducted using the Finite Element Method by means of ABAQUS/CAE software. The fatigue life was calculated using Miner’s Rule and the long-term distribution of stress range by the use of the two-parameter Weibull distribution. The cumulative damage ratio was estimated using the fatigue damage resulting from the stress range occurring at each load condition. For this purpose, a cargo hold model was first generated, which extends over the length of two holds (the mid-hold and half of each of the adjacent holds) and transversely over the full breadth of the hull girder. Following that, a submodel of the area of interest was extracted in order to calculate the hot spot stress of the connection and to estimate the fatigue life of the structural detail. Two hot spot locations were identified; one at the top layer of the inner bottom plate and one at the top layer of the hopper plate. The IACS Common Structural Rules (CSR) require that specific dynamic load cases for each loading condition are assessed. Following this, the dynamic load case that causes the highest stress range at each loading condition should be used in the fatigue analysis for the calculation of the cumulative fatigue damage ratio. Each load case has a different effect on ship hull response. Of main concern, when assessing the fatigue strength of the lower hopper knuckle connection, was the determination of the maximum, i.e. the critical value of the stress range, which acts in a direction normal to the weld toe line. This acts in the transverse direction, that is, perpendicularly to the ship's centerline axis. The load cases were explored both theoretically and numerically in order to establish the one that causes the highest damage to the location examined. The most severe one was identified to be the load case induced by beam sea condition where the encountered wave comes from the starboard. At the level of the cargo hold model, the model was assumed to be simply supported at its ends. A coarse mesh was generated in order to represent the overall stiffness of the structure. The elements employed were quadrilateral shell elements, each having four integration points. A linear elastic analysis was performed because linear elastic material behavior can be presumed, since only localized yielding is allowed by most design codes. At the submodel level, the displacements of the analysis of the cargo hold model to the outer region nodes of the submodel acted as boundary conditions and applied loading for the submodel. In order to calculate the hot spot stress at the hot spot locations, a very fine mesh zone was generated and used. The fatigue life of the detail was found to be 16.4 years which is lower than the design fatigue life of the structure (25 years), making this location vulnerable to fatigue fracture issues. Moreover, the loading conditions that induce the most damage to the location were found to be the various ballasting conditions.
Numerical Analysis and Design of Dielectric to Plasmonic Waveguides Couplers
In this work, efficient directional coupler composed of
dielectric waveguides and metallic film has been analyzed in details
by simulations using finite element method (FEM). The structure
consists of a step-index fiber with dielectric core, silica cladding, and
a metal nanowire parallel to the core. The results show that an
efficient conversion of optical dielectric modes to long range
plasmonic is possible. Low insertion losses in conjunction with short
coupling length and a broadband operation can be achieved under
certain conditions. This kind of couplers has potential applications for
the design of photonic integrated circuits for signal routing between
dielectric/plasmonic waveguides, sensing, lithography, and optical
storage systems. A high efficient focusing of light in a very small
region can be obtained.
Implicit Eulerian Fluid-Structure Interaction Method for the Modeling of Highly Deformable Elastic Membranes
This paper is concerned with the development of a
fully implicit and purely Eulerian fluid-structure interaction method
tailored for the modeling of the large deformations of elastic
membranes in a surrounding Newtonian fluid. We consider a
simplified model for the mechanical properties of the membrane, in
which the surface strain energy depends on the membrane stretching.
The fully Eulerian description is based on the advection of a modified
surface tension tensor, and the deformations of the membrane are
tracked using a level set strategy. The resulting nonlinear problem
is solved by a Newton-Raphson method, featuring a quadratic
convergence behavior. A monolithic solver is implemented, and we
report several numerical experiments aimed at model validation and
illustrating the accuracy of the presented method. We show that
stability is maintained for significantly larger time steps.
An Implicit Methodology for the Numerical Modeling of Locally Inextensible Membranes
We present in this paper a fully implicit finite element
method tailored for the numerical modeling of inextensible fluidic
membranes in a surrounding Newtonian fluid. We consider a highly
simplified version of the Canham-Helfrich model for phospholipid
membranes, in which the bending force and spontaneous curvature
are disregarded. The coupled problem is formulated in a fully
Eulerian framework and the membrane motion is tracked using
the level set method. The resulting nonlinear problem is solved
by a Newton-Raphson strategy, featuring a quadratic convergence
behavior. A monolithic solver is implemented, and we report several
numerical experiments aimed at model validation and illustrating
the accuracy of the proposed method. We show that stability is
maintained for significantly larger time steps with respect to an
explicit decoupling method.
Experimental and Finite Element Study of Bending Fatigue Failure: A Case Study on Main Shaft of a Gyrator Crusher
This study investigates the mechanism of a Gyratory crusher-located in Golgohar mining and industrial Co. specifically with a focus on stresses distribution and fatigue failure of its main shaft. At first step, the cross section of the fractured shaft is studied, and the crack growth is analyzed. Then, the rotational motion of the shaft and the oil temperature of oil circuit of equipment are monitored. Condition monitoring is used to help finding a better modification. Based on the results of this study, the main causes of shaft failure are identified, and corrective solution is offered to increase crusher performance, especially its main shaft life. To predict the efficiency of the proposed modification, finite element simulation is performed, and its results are compared with the similar modified cases. The comparison and interpretation of simulation results confirm the efficiency of proposed corrective method.
Failure Modes and Bearing Capacity Estimation for Strip Foundations in C-ɸ Soils: A Numerical Study
In this study, typical c-ɸ soils subjected to loadings were assessed with a view to understand the general stress distribution and settlement behaviour of the soils under drained conditions. Numerical estimations of the non-dimensional bearing capacity factors, Nq and Nγ for varied angles of friction in the soil mass were obtained using PLAXIS. Ultimate bearing capacity values over a Ф range of 0-30 degrees were also computed and compared with analytical results obtained from the traditional simplified uncoupled approach of Terzaghi and Meyerhof. Results from the numerical study agree well with theoretical findings.
Dynamic Soil Structure Interaction in Buildings
Since the evolution of computational tools and simulation software, there has been considerable increase in research on Soil Structure Interaction (SSI) to decrease the computational time and increase accuracy in the results. To aid the designer with a proper understanding of the response of structure in different soil types, the presented paper compares the deformation, shear stress, acceleration and other parameters of multi-storey building for a specific input ground motion using Response-spectrum Analysis (RSA) method. The response of all the models of different heights have been compared in different soil types. Finite Element Simulation software, ANSYS, has been used for all the computational purposes. Overall, higher response is observed with SSI, while it increases with decreasing stiffness of soil.
Experimental Study on Flexural Strength of Reinforced Geopolymer Concrete Beams
This paper presents the flexural response of Reinforced Geopolymer Concrete (RGPC) beams. A commercial finite element (FE) software ABAQUS has been used to perform a structural behavior of RGPC beams. Using parameters such: stress, strain, Young’s modulus, and Poisson’s ratio obtained from experimental results, a beam model has been simulated in ABAQUS. The results from experimental tests and ABAQUS simulation were compared. Due to friction forces at the supports and loading rollers; slip occurring, the actual deflection of RGPC beam from experimental test results were slightly different from the results of ABAQUS. And there is good agreement between the crack patterns of fly ash-based geopolymer concrete generated by FE analysis using ABAQUS, and those in experimental data.
Analysis of Residual Stresses and Angular Distortion in Stiffened Cylindrical Shell Fillet Welds Using Finite Element Method
In this paper, a two-dimensional method is developed to simulate the fillet welds in a stiffened cylindrical shell, using finite element method. The stiffener material is aluminum 2519. The thermo-elasto-plastic analysis is used to analyze the thermo-mechanical behavior. Due to the high heat flux rate of the welding process, two uncouple thermal and mechanical analysis are carried out instead of performing a single couple thermo-mechanical simulation. In order to investigate the effects of the welding procedures, two different welding techniques are examined. The resulted residual stresses and distortions due to different welding procedures are obtained. Furthermore, this study employed the technique of element birth and death to simulate the weld filler variation with time in fillet welds. The obtained results are in good agreement with the published experimental and three-dimensional numerical simulation results. Therefore, the proposed 2D modeling technique can effectively give the corresponding results of 3D models. Furthermore, by inspection of the obtained residual hoop and transverse stresses and angular distortions, proper welding procedure is suggested.
Effect of Mechanical Loading on the Delamination of Stratified Composite in Mode I
The present study is based on the three-dimensional digital analysis by the finite elements method of the mechanical loading effect on the delamination of unidirectional and multidirectional stratified composites. The aim of this work is the determination of the release energy rate G in mode I and the Von Mises equivalent constraint distribution along the damaged area under the influence of several parameters such as the applied load and the delamination size. The results obtained in this study show that the unidirectional composite laminates have better mechanical resistance one the loading line than the multidirectional composite laminates.
Nonlinear Finite Element Modeling of Deep Beam Resting on Linear and Nonlinear Random Soil
An accuracy nonlinear analysis of a deep beam resting on elastic perfectly plastic soil is carried out in this study. In fact, a nonlinear finite element modeling for large deflection and moderate rotation of Euler-Bernoulli beam resting on linear and nonlinear random soil is investigated. The geometric nonlinear analysis of the beam is based on the theory of von Kàrmàn, where the Newton-Raphson incremental iteration method is implemented in a Matlab code to solve the nonlinear equation of the soil-beam interaction system. However, two analyses (deterministic and probabilistic) are proposed to verify the accuracy and the efficiency of the proposed model where the theory of the local average based on the Monte Carlo approach is used to analyze the effect of the spatial variability of the soil properties on the nonlinear beam response. The effect of six main parameters are investigated: the external load, the length of a beam, the coefficient of subgrade reaction of the soil, the Young’s modulus of the beam, the coefficient of variation and the correlation length of the soil’s coefficient of subgrade reaction. A comparison between the beam resting on linear and nonlinear soil models is presented for different beam’s length and external load. Numerical results have been obtained for the combination of the geometric nonlinearity of beam and material nonlinearity of random soil. This comparison highlighted the need of including the material nonlinearity and spatial variability of the soil in the geometric nonlinear analysis, when the beam undergoes large deflections.
Finite Element Method for Modal Analysis of FGM
Modal analysis of a FGM plate containing the ceramic phase of Al2O3 and metal phase of stainless steel 304 was performed using ABAQUS, with the assumptions that the material has an elastic mechanical behavior and its Young modulus and density are varying in thickness direction. For this purpose, a subroutine was written in FOTRAN and linked with ABAQUS. First, a simulation was performed in accordance to other researcher’s model, and then after comparing the obtained results, the accuracy of the present study was verified. The obtained results for natural frequency and mode shapes indicate good performance of user-written subroutine as well as FEM model used in present study. After verification of obtained results, the effect of clamping condition and the material type (i.e. the parameter n) was investigated. In this respect, finite element analysis was carried out in fully clamped condition for different values of n. The results indicate that the natural frequency decreases with increase of n, since with increase of n, the amount of ceramic phase in FGM plate decreases, while the amount of metal phase increases, leading to decrease of the plate stiffness and hence, natural frequency, as the Young modulus of Al2O3 is equal to 380 GPa and the Young modulus of stainless steel 304 is equal to 207 GPa.