Dear Label_Name,

Greetings! I hope you are well.

I am writing on behalf of an international publishing house, Science and Education Publishing, USA (http://www.sciepub.com).

Science and Education Publishing (SciEP) is a multidisciplinary open access publisher of Journals covering the fields of Science, Technology and Medicine. All SciEP publications are subject to high-quality peer review, editorial and production processes and are freely available online.

Welcome to the April, 2013 edition of SciEP newsletter, providing monthly summaries of articles on science & education related topics of interest to scientist and professionals. The mission of our journals is to provide an open access forum for publishing new findings on issues pertinent to scientific and academic field. We also invite you to submit original research papers, review articles, short reports, or case studies are considered for publication in our Multidisciplinary Journals.

American Journal of Public Health Research
DOI: 10.12691/ajphr-1-1-3

The Implementation of a Tele-Homecare System with Service Oriented Architecture and HL7 Message Transmission Standard
Ching-Sung Wang

American Journal of Medical and Biological Research
DOI: 10.12691/ajmbr-1-1-2

The Influence of Ethanol on Pyruvate Kinases Activity in Vivo, in Vitro, in Silico
Lelevich Sergey Vladimirovich, Khrustalev Vladislav Victorovich, Barkovsky Eugene Victorovich, Shedogubova Tatyana Aleksandrovna

American Journal of Medical Sciences and Medicine
DOI: 10.12691/ajmsm-1-1-1

Clinical and Histological Features in Wegener Granulomatosis
Li Cavoli Gioacchino, Passantino Rita, Bono Luisa, Tortorici Calogera, Ferrantelli Angelo, Giammarresi Carlo, Rotolo Ugo

American Journal of Infectious Diseases and Microbiology
DOI: 10.12691/ajidm-1-1-2

The Curse of Lymphatic Filariasis: Would the Continual Use of Diethylcarbamazine Eliminate this Scourge in Papua New Guinea?
Anne Lanham, Lillian Mwanri

We have released issues of some journals on February 2013.
American Journal of Public Health Research
American Journal of Medical Sciences and Medicine
American Journal of Medical and Biological Research
American Journal of Pharmacological Sciences
American Journal of Clinical Medicine Research
American Journal of Nursing Research
American Journal of Cancer Prevention
American Journal of Medicine Studies
American Journal of Infectious Diseases and Microbiology
American Journal of Cardiovascular Disease Research
American Journal of Biomedical Research
American Journal of Mechanical Engineering
American Journal of Microbiological Research
Journal of Preventive Medicine
American Journal of Food and Nutrition
American Journal of Hypertension Research
American Journal of Sports Science and Medicine
Journal of Cancer Research and Treatment
American Journal of Zoological Research
Journal of Biomedical Engineering and Technology
American Journal of Microbiological Research
Biomedical Science and Engineering

We are pleased to announce that submissions are now being accepted for peer review for all the journals across all our subject areas in SciEP.

Please submit your manuscript via our online submission system:
http://www.sciepub.com/author/signin

To unsubscribe from our mailing list, please click on the following link: unsubscribe

I am looking forward to hearing from you.

Kind regards,

David Rain
Email:join@sciepub.com
Science and Education Publishing (http://www.sciepub.com)

Free Vibration Analysis of Clamped Skew Composite Plate Using Finite Element Method

Mohannad Aziz Kadhom*, Luay S. Al-Ansari, Faeq Hamid Gburi

Technical Institute of AL-Diwaniyah, AL-Furat Al-Awsat Technical University (ATU), Diwaniyah, Iraq

Abstract: This study is presented of two models were built to simulate the skew composite plates and these models are SHELL281 and SOLID186 element. Three types of laminates arrangement were used to study the effect of skew angle, fiber orientation and length-to-width ratio on the first five natural frequencies of the clamped composite plates. The skew angle was changed from (0°) to (75°) with (15°) steps while the fiber orientation was changed from (0°) to (90°) with (15°) steps and the length-to-width ratio was (0.5, 1, 1.25, 1.5 and 2). The results showed that the first five frequency parameters increase for three layouts when the skew angle increases at constant fiber orientation and constant length-to-width ratio. Also, the fiber orientation has a slight effect on the first five frequency parameters in (0/β/0) layout. The effect of length-to-width ratio on the first five frequency parameters depends on the skew angle (smaller or larger than (60o)), fiber orientation and the laminates layout. Finally, the values of the first five frequency parameters of (β/0/β) and (β/β/β) layouts are close to each other.

Keywords: free vibration, skew plate, frequency parameter, length-to-width ratio, ANSYS software, composite plate, finite element method

Received September 03, 2021; Revised October 19, 2021; Accepted October 28, 2021

Citation: Mohannad Aziz Kadhom, Luay S. Al-Ansari, Faeq Hamid Gburi, Free Vibration Analysis of Clamped Skew Composite Plate Using Finite Element Method. American Journal of Mechanical Engineering, 9 (1). 41-56. January 2021. doi:10.12691/ajme-9-1-5

1. Introduction

Composite materials can be defined as a combination of two or more materials in a macroscopic scale by mechanical and chemical bonding in order to produce a new material that has desired mechanical properties. Because of increasing strength-to-weight ratio and specific modulus, composite plates are widely used in the manufacturing of structural members. The laminated composite plate is a kind of composite materials applications and it is used commonly in automobiles, robotics, day to day appliances, building industry, aerospace structures, nuclear structures, marine vehicles and ships. The laminated composite skew plates are also used in these applications and the analysis of this type of composite plate is very difficult for different reasons such as increased number of material constants, increased number of layers and different structural couplings. Due to these difficulties, the exact mathematical solutions of the free vibration problems of skew composite plates are not available, therefore, the approximate methods are used to analyze the free vibration problem.

In the last few decades, many researchers investigated extensively the free vibration behavior of skew plates. The finite element method is one of the methods used to analyze the free vibration problem of the skew composite plates. [1] used method to estimate the natural frequencies of simply-supported rhombic plate. For clamped rhombic plates, [2] used FEM to obtain natural frequencies. [3] studied mode shapes and natural frequencies of skew plates with clamped boundary condition by adopting Galerkin FEM. [4] calculated the modes of free-free CFRP plates considering the effect of transverse shear deformation using the finite element method. [5] used finite element method to study the effect of skew angles, fiber orientation and boundary conditions on the vibration behavior of rhombic orthotropic plates. [6] used a Co finite element method to investigate the free vibration problem of sandwich and laminated beam and plate. [7] calculated mode shapes and natural frequencies of isotropic skew plates by adopting the hierarchical finite element method. [8] determined (a) the effects of skew angle, (b) angle of fiber orientation, (c) aspect ratio (length-to-width ratio) and (d) length-to-thickness ratio on the natural frequency of one lamina glass/epoxy and graphite/epoxy skew plates. [9] combined the Mindlin plate theory and the finite element method to study the effect of partially supported boundaries on the free vibration behavior of right triangular orthotropic plates. [10] used C0 continuous isoparametric quadrilateral element to study the free vibration of thin and thick cantilever skew plates using. [11,12] used finite element method to study the free vibration problem of the deep cylindrical panel and thin isotropic conical panels. [13] studied the free vibration of laminated skew plates using a general high precision triangular bending element. [14] used a nine-node degenerated shell element to analyze the free vibration of isotropic plates based on the first-order shear deformation theory. [15] used FEM to study the free vibrations problem of skew composite plates with large amplitude. [16] developed a new finite element model to simulate anisotropic and laminated doubly curved moderately thick composite shell panels using the first-order shear deformation theory. [17] Used isoparametric finite element model for studying the free vibration behavior of different types of skew plates. [18] studied the static and free vibration of general shell structures using the element-free Galerkin method. [19] used a higher-order finite element to analyze the free vibration problem of FG curved panels. [20] used finite element method to study flexural natural frequencies of simply-supported (SSSS) and clamped (CCCC) skew plates assuming isotropic and laminated composite materials. They studied the effects of aspect ratio, length-to-thickness-ratio and the skew angle the natural frequency of isotropic skew plates while they studied the additional factors such as fiber orientation angle and numbers of layers. They used these factors to compare between CQUAD4 element and CQUAD8 element and studied the accuracy of their results. [21] investigated experimentally the free vibration of skew laminated plates and compared between experimental and finite element results. [22] used ANSYS APDL software to study free vibration behaviour of rectangular composite laminate plate for different boundary conditions. He studied the effects of modulus ratios side-to-thickness ratios and aspect ratios on the first five natural frequencies of plate for different boundary conditions. [23] used a commercial finite element program (ANSYS 14.5) to compute the natural frequency of simply supported and clamped skew plates with circular cutout. He used 8-Nodes isoparametric shell element to study the effects of skew angles, side to thickness ratios, number of plies and type of layer arrangements on the frequency of skew plates. [24] used a strong formulation finite element method (SFEM) and linear theory for moderately thick structures in order to study the first three natural frequencies of laminated plates. [25] used the finite element method to study the free vibration analysis of finite square and skew laminated plates using variable stiffness plies for implementing the curvilinear fiber reinforcement. It is assumed that due to employing tow placement technology in each ply of variable stiffness composite laminated plate, the fiber reinforcement orientation angle is changed linearly with respect to longitudinal geometry coordinate. They compared their results with similar items available in the literature to show the accuracy of new formulation at different curvilinear fiber angles, different geometries and various boundary conditions. [26] analyzed the free vibration problem of simply supported stiffened plates with different geometries using finite element technique. They studied the skew angle effects and the stiffener properties on the free vibration of circular, square and rectangular (with different aspect ratio) plates. [27] Analyzed delamination size on the vibration behavior experimentally and theoretically. They adopted finite element method to simulate the composite plate with different delamination sizes, boundary conditions and layers arrangement using two software (ANSYS and ABAQUS). Also, they used Rayleigh-Ritz method to calculate the natural frequencies of the composite plate. Experimentally, they measured the natural frequency values for carbon fiber-reinforced polymer. They found that the natural frequency decreases with increasing delamination size for all boundary conditions.

In this work, two models based on FEM were used to simulate three layers laminated composite skew plates with clamped boundary condition. These two models were built by ANSYS APDL software version 17.2 using SHELL281 and SOLID186 element. These models were used to investigate the influences of skew angle, length-to-width ratio (aspect ratio), angle of fiber orientation and layers layout of the natural frequency of skew plate. This use results in finding the relation between these parameters especially the layers layout and fiber orientation for the suggestion layers layout. Also, the calculation accuracy of the two elements is studied by comparing the results of the two models.

This paper consists of five sections in addition to “introduction”, the next section (i.e. section 2) deals with the problem distribution and cases studied in this work. The finite element models and the two elements used in this work are described in Section (3). In Section (4), the comparison between the natural frequency of the new finite element models with those available in the literature to check the calculation accuracy of the new models. The results are displayed and discussed in Section (5). Finally, the conclusions and future works are summarized in Section (6).

2. Theory and Formula Problem Description and Studied Cases:

In this work, three layers laminated composite plate was considered to study the influences of the skew angle, length-to-width ratio (aspect ratio), angle of fiber orientation and layers layout of the natural frequency of clamped skew composite plate. The geometry, dimensions and layers layout of the skew composite plate are illustrated in Figure 1 [28]. The width of skew plate (a) is 0.3 m and the aspect ratio (a/b) is 0.5, 1.0, 1.25, 1.5 and 2. The skew angle (θ) is 0°, 15°, 30°, 45°, 60° and 75°.

The material properties of the layer of composites are Longitudinal Modulus of Elasticity (El) = 38.07 GPa; Transverse Modulus of Elasticity (Et)=8.1 GPa; Modulus of Rigidity G₁₂= G₁₃= G₂₃= 3.05 GPa; Poisson Ratio ν₁₂= ν₁₃= ν₂₃=0.22; Density (ρ) =2200 kg/m³.

Figure 1. Dimension of Skew Plate

3. Finite Element Models

In this work, ANSYS APDL version (17.2) was used to simulate the three layers laminated composite skew plates. Two different models were used the first one used SOLID186 element and the second model used SHELL281 (see Figure 2). “SOLID 186 is a higher order 3-D 20-node solid element that exhibits quadratic displacement behavior. The element is defined by 20 nodes having three degrees of freedom per node: translations in the nodal x, y, and z directions. The element supports plasticity, hyper elasticity, creep, stress stiffening, large deflection, and large strain capabilities. It also has mixed formulation capability for simulating deformations of nearly incompressible elastoplastic materials, and fully incompressible hyperplastic materials”. The number of elements in first model is about (275000 - 300000) elements and the number of nodes is about (550000 – 600000) nodes. “SHELL 281 is well-suited for linear, large rotation, and/or large strain nonlinear applications. Change in shell thickness is accounted for in nonlinear analyses. The element accounts for follower (load stiffness) effects of distributed pressures” [29,30,31,32]. The number of elements in first model is about (3700- 4000) elements and the number of nodes is about (11500 – 14000) nodes.

Figure 2. Dimensions of CompositeSkew Plate

4. Validation

In order to check the validation for the present elements, the non-dimensional fundamental frequency coefficient (Kf) is compared with those available in the literature. The non-dimensional fundamental frequency coefficient (Kf) is calculated as:

Where:

a – The width of plate (m),

t – The thickness of plate (m),

– The frequency (rad/sec),

Et- Transverse Modulus of Elasticity (N/m²).

-Mass density per unit area.

Table 1 shows a good agreement between the non-dimensional fundamental frequency coefficient results of the elements used in this work and those of available in the literature for the simply-supported square skew plate with anti-symmetric angle-ply arrangement.

Table 1. Fundamental Frequency of Simply-Supported Square Skew Plate with Anti-symmetric Angle-ply Arrangement

5. Result Discussions

In this work, three parameters affecting the natural frequencies of the skew laminated composite plate were studied, namely skew angle, fiber orientation and length-to-width ratio (aspect ratio). Three types of laminates layout were considered and these layouts were (0/β/0), (β/0/β) and (β/β/β).

Figure 3 shows the variation of the first frequency parameter of the three layouts due to skew angle and fiber orientation for different length-to-width ratio. The fiber orientation does not affect the first frequency parameters in (0/β/0) layout but the skew angle has a significant effect on the first frequency parameter. When the skew angle increases, the first frequency parameter increases (see Figure 4). In (β/0/β) and (β/β/β) layouts, both skew angle and fiber orientation effect on first frequency parameter. When the skew angle increases, the first frequency parameter increases at the same fiber orientation (see Figure 4). While the first frequency parameter decreases when the fiber orientation increases for the skew angle smaller than or equals (62°) (see Figure 4). But for the skew angle larger than (62°), the first frequency parameter increases when the fiber orientation increases at length-to-width ratio smaller than (1) and decreases at length-to-width ratio larger than or equals (1). This behavior can be analyzed as:

1. In general, the Young Modulus of laminated composite plate decreases when the fiber orientation (β) increases. This variation in Young Modulus leads to decreasing the first frequency parameter.

2. The increase in the skew angle leads to an increase in the first frequency parameter because of increasing the effect of boundary condition.

3. In this work, the length of plate is considered constant and equals (0.3 m). The first frequency parameter increases by increasing the length-to-width ratio due to the increase in the area of the plate.

4. The combination of these three factors leads to varying the first frequency parameter.

Figure 5 and Figure 6 show the variation of the second frequency parameter of the three layouts due to skew angle and fiber orientation for different length-to-width ratio. In (0/β/0) layout, the fiber orientation does not affect the second frequency parameters while the first frequency parameter increases when the skew angle increases (see Figure 6). In (β/0/β) and (β/β/β) layouts, the second frequency parameters decrease when the fiber orientation increases for the skew angle larger than (62°) at a length-to-width ratio larger than or equals (1) (see Figure 6). While the second frequency parameters increase when the fiber orientation increases for the skew angle smaller than (62°).

The same variation in third, fourth and fifth frequency parameters are shown in Figure 7 - Figure 12. The same effects of the skew angle, fiber orientation and length-to-width ratio on the first five frequency parameters are shown in spite of the increase in mode number.

Figure 3. The First Prequency Parameter Variation with Skew Angle and Fiber Orientation for Different Length-to Width Ratio (3D-Graph)

Figure 4. The First Frequency Parameter Variation with Skew Angle and Fiber Orientation for Different Length-to-Width Ratio (2D-Graph)

Figure 5. The Second Frequency Parameter Variation with Skew Angle and Fiber Orientation for Different Length-to-Width Ratio (3D-Graph)

Figure 6. The Second Frequency Parameter Variation with Skew Angle and Fiber Orientation for Different Length-to-Width Ratio (2D-Graph)

Figure 7. The Third Frequency Parameter Variation with Skew Angle and Fiber Orientation for Different Length-to-Width Ratio (3D-Graph)

Figure 8. The Third Frequency Parameter Variation with Skew Angle and Fiber Orientation for Different Length-to-Width Ratio (2D-Graph)

Figure 9. The Fourth Frequency Parameter Variation with Skew Angle and Fiber Orientation for Different Length-to-Width Ratio (3D-Graph)

Figure 10. The Fourth Frequency Parameter Variation with Skew Angle and Fiber Orientation for Different Length-to-Width Ratio (2D-Graph)

Figure 11. The Fifth Frequency Parameter Variation with Skew Angle and Fiber Orientation for Different Length-to-Width Ratio (3D-Graph)

Figure 12. The Fifth Frequency Parameter Variation with Skew Angle and Fiber Orientation for Different Length-to-Width Ratio (2D-Graph)

6. Conclusions

In this work, the first five frequency parameters of three layout laminated skew composite plates were studied. From the results, the following points can be concluded:

The skew angle affects the first five frequency parameters in the same way. The first five frequency parameters increase for three layouts when the skew angle increases at constant fiber orientation and constant length-to-width ratio.

1. The fiber orientation has a slight effect on the first five frequency parameters in (0/β/0) layout. But in (β/0/β) and (β/β/β) layouts, the fiber orientation effects on the first five frequency parameters and its effect depends on the skew angle (smaller or larger than (60°)) and length-to-width ratio (smaller or larger than (1)).

2. The length-to-width ratio affects the first five frequency parameters and this effect depends on the skew angle (smaller or larger than (60°)), fiber orientation and the laminates layout.

3. The laminates layouts affect the values of the first five frequency parameters. The values of the first five frequency parameters of (β/0/β) and (β/β/β) layouts are close to each other.

4. The shell and solid models are very closed to each other and the shell model is simpler and faster than the solid model.

In future work, the four and five laminates layouts can be investigated to study the effects of skew angle, fiber orientation, length-to-width ratio, types of supporting and laminate arrangement on the first five frequency parameters.

References

[1]   | Laura. P.A., and J. Grosson., 1968. Fundamental Frequency of Vibration of Rhombic Plates, Journal of Acoustical Society of America, vol. 44, pp. 823-824.
[2]   | Monforton. G.R., 1968. Finite Element Displacement and Vibration Analysis of Skew Plates, Report18, Division of Solid Mechanics, Structural and Mechanical Design Case, Western Reserve University, Cleveland, Ohio.
[3]   | Durvasula. S., 1969. Natural Frequencies and Modes of Clamped Skew Plates, American Institute of Aeronautics and Astronautics Journal, vol.7, pp.1164-1166.
[4]   | Cawley P, Adams RD (1978). The predicted and experimental natural modes of free-free CFRP plates. Journal of Composite Materials 12: 336-347.
[5]   | Malhotra. S.K., N. Ganesan., and M.A. Veluswami., 1988. Effect of Fiber Orientation and Boundary Conditions on The Vibration Behavior of Orthotropic Rhombic Plates, Journal of Composites, vol.19, No.2, pp.127-132.
[6]   | Mallikarjuna, Kant T. Free vibration of symmetrical laminated plates using a higher order theory with finite element technique. International Journal for Numerical Methods in Engineering, vol.28, pp. 1875-1889, 1989.
[7]   | Bardell. N.S., 1992. The Free Vibration of Skew Plates Using the Hierarchical Finite Element Method, Computer and Structures, vol. 45, pp 841-874.
[8]   | Krishnan. A., and Deshpande J.V., 1992. Vibration of Skew Laminates, Journal of Sound and Vibrations, vol.153, No.2, pp. 351-358.
[9]   | Mirza, S., and Y. Alizadeh. 1994. “Dynamics of Partially Supported Right Triangular Orthotropic Plates.” Journal of Sound and Vibration 178 (4): 567-575.
[10]   | McGee OG, Butalia TS. (1994). Natural frequencies of shear deformable cantilevered skew thick plates. Journal of Sound and Vibration. 176: 351-376.
[11]   | Bardell NS, Dunsdon JM, Langley RS. (1997). On the free vibration of completely free, open, cylindrically curved, isotropic shell panels. Journal of Sound and Vibration 207: 647-669.
[12]   | Bardell NS, Dunsdon JM, Langley RS. (1998). Free vibration of thin, isotropic, open conical panels. Journal of Sound and Vibration. 217: 297-320.
[13]   | Krishna Reedy. A.R., and R. Palaninathan., 1999. Free Vibration of Skew Laminates, Computers and Structures, vol.70, pp. 415- 423.
[14]   | Lee SJ, Han SE. (2001). Free vibration analysis of plates and shells with a nine-node assumed natural degenerated shell element. Journal of Sound and Vibration. 241(4): 605-633.
[15]   | Moloy. K. Singha., and M. Ganapathi., 2004. Large Amplitude Free Flexural Vibrations of Laminated Skew Plates, Journal of Non-Linear Mechanics, vol.39, pp.1709-1720.
[16]   | Hossain SJ, Sinha PK, Sheikh AH (2004).. A finite element formulation for the analysis of laminated composite shells. Computers and Structures. 82:1623-1638.
[17]   | Garg. A.K., R.K. Khare., and T. Kant., 2006. Free Vibration of Skew Fiber-Reinforced Composite and Sandwich Laminates Using a Shear Deformable Finite Element Model, Journal of Sandwich Structures and Materials, vol.8, pp. 33-53.
[18]   | Liu L, Chua LP, Ghista DN. (2006). Element free Galerkin method for static and dynamic analysis of spatial shell structures. Journal of Sound and Vibration. 295:388-406.
[19]   | Pradyumna S, Bandyopadhyay JN. (2008). Free vibration analysis of functionally graded curved panels using a higher-order finite element formulation. Journal of Sound and Vibration 318: 176-192.
[20]   | C.V. Srinivasa, Y.J. Suresh and W.P. Prema Kumar,” Free flexural vibration studies on laminated composite skew plates”; International Journal of Engineering, Science and Technology Vol. 4, No. 4, 2012, pp. 13-24.
[21]   | Srinivasa. C. V., Suresh.Y.J., Prema Kumar W.P., Experimental and finite element studies on free vibration of skew plates, International Journal of Advanced Structural Engineering, Vol.6, pp.1-11, 2014.
[22]   | Deepesh Bhavsar;” Free Vibration Analysis of Composite Laminate for Different Simply Supported Boundary Condition”; International Journal for Scientific Research & Development (IJSRD) Vol. 3, Issue 01, 2015 | ISSN (online): 2321-0613.
[23]   | K. Sai Vivek; “Free Vibration of Skew Laminated Composite Plates with Circular Cutout by Finite Element Method”; International Journal Of Modern Engineering Research (IJMER); Vol. 6; Iss. 6; June 2016 ;16
[24]   | Nicholas Fantuzzi, Francesco Tornabene, Michele Bacciocchi and Antonio J.M. Ferreira;” On the Convergence of Laminated Composite Plates of Arbitrary Shape through Finite Element Models”; Journal of Composites Science 2018, 2, 16.
[25]   | Vahid Khalafia and Jamshid Fazilatib. “Free vibration analysis of variable stiffness composite laminated thin skew plates using IGA”; Journal of Theoretical and Applied Vibration and Acoustics 4(2) 171-188 (2018).
[26]   | Bhargav Reddy Isanaka, M. Abdul Akbar, Biraja Prasad Mishra and Vinod Kushvaha. “Free vibration analysis of thin plates: Bar e versus Stiffened”; Engineering Research Express 2 (2020) 015014.
[27]   | M. Imran, R. Khan, and S. Badshah;” Experimental, analytical, and finite element vibration analyses of delaminated composite plates”; Scientia Iranica, Transactions B: Mechanical Engineering 28 (2021) 231-240.
[28]   | Muhammad Imran, Rafiullah Khan and Saeed Badshah,” A Review on the Vibration Analysis of Laminated Composite Plate”; Pakistan Journal of Scientific and Industrial Research Series A: Physical Sciences 2019 62A(2) 124-134.
[29]   | ANSYS Mechanical APDL Element Reference, ANSYS, Inc., 2016.
[30]   | T. Stolarski, Y.Nakasone and S. Yoshimoto, “Engineering Analysis with ANSYS Software”; Elsevier Butterworth-Heinemann; 2006.
[31]   | Jones. R.M., 1975. Mechanics of Composite Materials, McGraw-Hill, New York.
[32]   | F.H. Gburi, L.S. Al-Ansari, M.A. Khadom, A.A. Al-Saffar, “Free Vibration of Skew Isotropic Plate Using ANSYS”, Journal of Mechanical Engineering Research and Developments, Vol. 43, No.5, pp. 472-486, 2020.