Abstract
In this work, we have investigated the nonlinear instability behavior of a thin laminated composite circular plate subjected to a tensile periodic load. A nonlinear system with quadratic and cubic nonlinear terms, referred to the effect of large deflection, is obtained. The method of multiple time scales is used to obtain periodic solutions for the second- and third-order systems, which directly leads to the solvability conditions of the two orders including the nonlinearity terms. Different cases of resonance with the instability criteria are taken into account. For the numerical application, the effect of the tensile periodic load was studied on a plate of graphite fabric-carbon material.
Similar content being viewed by others
References
Bert, C.W., Birman, V.: Dynamic instability of shear deformable antisymmetric angle-ply plates. Int. J. Solids Struct. 23(7), 1053–1061 (1987)
Darabi, M., Darvizeh, M., Darvizeh, A.: Non-linear analysis of dynamic stability for functionally graded cylindrical shells under periodic axial loading. Compos. Struct. 83, 201–211 (2008)
Bolotin, V.V.: The Dynamic Stability of Elastic Systems. Holden-Day, San Francisco (1964)
Sahu, S.K., Datta, P.K.: Research advances in the dynamic stability behavior of plates and shells: 1987–2005 part I: Conservative systems. Appl. Mech. Rev. 60(2), 65–75 (2007)
Srinivasan, R.S., Chellapandi, P.: Dynamic stability of rectangular laminated composite plates. Comput. Struct. 24(2), 233–238 (1986)
Bert, C.W., Birman, V.: Dynamic instability of shear deformable antisymmetric angle-ply plates. Int. J. Solids Struct. 23(7), 1053–1061 (1987)
Reddy, J.N.: Mechanics of laminated composite plates and shells: theory and analysis,2nd ed., CRC Press (2003)
Moorthy, J., Reddy, J.N., Plaut, R.H.: Parametric instability of laminated composite plates with transverse shear deformation. Int. J. Solids Struct. 26(7), 801–811 (1990)
Sridhar, S., Mook, D.T., Nayfeh, A.H.: Nonlinear resonances in the forced responses of plates. Part I: Symmetric responses of circular plates. J. Sound Vib. 41, 359–373 (1975)
Sridhar, S., Mook, D.T., Nayfeh, A.H.: Nonlinear resonances in the forced responses of plates. Part II: Asymmetric responses of circular plates. J. Sound Vib. 59, 159–170 (1978)
Hadian, J., Nayfeh, A.H.: Modal interaction in circular plates. J. Sound Vib. 142, 279–292 (1990)
Walker, M.: Multiobjective design of laminated plates for maximum stability using the finite element method. Comput. Struct. 54(2), 389–393 (2001)
Hu, H.T., Tzeng, W.L.: Buckling analysis of skew laminate plates subjected to uniaxial inplane loads. Thin-Walled Struct. 38(1), 53–77 (2000)
Baltaci, A., Sarikanat, M., Yildiz, H.: Static stability of laminated composite circular plates with holes using shear deformation theory. Finite Elem. Anal. Des. 43(11–12), 839–846 (2007)
Powmya, A., Narasimhan, M.C.: Free vibration analysis of axisymmetric laminated composite circular and annular plates using Chebyshev collocation. Int. J. Adv. Struct. Eng. (IJASE) 7(2), 129–141 (2015)
Coman, C.D., Haughton, D.M.: Localized wrinkling instabilities in radially stretched annular thin films. Acta Mech. 185(3–4), 179–200 (2006)
Coman, C.D., Haughton, D.M.: On some approximate methods for the tensile instabilities of thin annular plates. J. Eng. Math. 56(1), 79–99 (2006)
Popov, A.A.: Parametric resonance in cylindrical shells: a case study in the nonlinear vibration of structural shells. Eng. Struct. 25, 789–799 (2003)
Amabili, M.: Nonlinear Vibrations and Stability of Shells and Plates. Cambridge University Press, New York (2008)
Fung, Y.C.: Foundation of Solid Mechanics. Prentice-Hall, Englewood Cliffs (1965)
Nayfeh, A.H.: Introduction to perturbation techniques. Wiley, Hoboken (2011)
Zakaria, Kadry, Sirwah, Magdy A., Alkharashi, Sameh: Instability through porous media of three layers superposed conducting fluids. Eur. J. Mech. B Fluids 28, 259–270 (2009)
Allaire, G., Delgado, G.: Stacking sequence and shape optimization of laminated composite plates via a level-set method. J. Mech. Phys. Solids 97, 168–196 (2016)
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendices
Appendix A1
The elements of the transformed reduced stiffness matrix \([{\bar{Q}}_{ij}]^{(k)}\) are
where [T] is the transformation matrix between the principal material coordinates and the plate’s coordinates given by
and \(\vartheta \) is the orientation angle.
Appendix A2
In general, we considered unidirectional high-strength fiber-reinforced laminae of graphite/carbon materials, so the composite laminate corresponds to a collection of laminae which are orientated in different directions. This sequence of different orientations is termed the lamination scheme or stacking sequence. The stacking sequence and material properties of individual lamina provide additional flexibility to the designers to tailor the stiffness and strength of the laminate. It should be noticed that the material parameters of a composite beam or plate are varying in three cases:
-
In a circular plate with bi-dimensional varying thickness, i.e., \(h=h(r)\),
-
In a plate made of a functionally graded polar orthotropic material,
-
In a plate made of anisotropic layers.
To clarify why the physical properties of an orthotropic composite plate are invariant and moreover independent of the radial and circumferential coordinates, we will give an example of how a plate of anisotropic layers is fabricated
-
Let \(\Omega \) be a symmetric laminated composite structure composed of the superposition of 2N anisotropic layers, each one of constant thickness \(\varepsilon >0\), but this time each one with a variable shape \({\mathcal {M}}_i \subset {\mathcal {M}}\), where \({\mathcal {M}}\) is a regular sub-domain of \(R^2\). By abuse of notation, we denote \(\Omega \) as:
$$\begin{aligned} \Omega = \{{\mathcal {M}}_{i}\},~~~i= -N,\ldots ,-1,1,\ldots .N. \end{aligned}$$Since we suppose \(\Omega \) symmetric, i.e., \({\mathcal {M}}_{-i} = {\mathcal {M}}_{i}\), we consider only N layers, so we rather write \(\Omega = \{{\mathcal {M}}_{i}\},~~~i=1,2,\ldots ,N\). The index i grows from the inside to the outside of the laminated composite structure, see Fig. 2.
We considered each laminate to be an orthotropic material and unidirectionally reinforced. For the purpose of this study, each layer is a non-homogeneous two-phase material, where each “hole” is filled with another “weak” material with different physical properties (weight, electric or heat conductivity, etc.). We will denote this weak material as \(A_{0}\) (in brown color). Let \(\tau _{i}\) be the characteristic function of the i-layer. According to the classical theory of plates, the composite structure \(\Omega \) is characterized by the superposition of the elastic properties of each layer. Thus, the extensional stiffness tensor A reads [23]:
where \(A_i\) is the extensional stiffness of the i-layer, meanwhile the bending stiffness tensor D reads
Now, to fabricate a laminated composite plate with invariant material properties, we assumed each ply of an orthotropic material and unidirectionally reinforced. Each hole will inject with the same fiber material, which leads to the characteristic function \(\tau _{i}\) for each i-layer will be invariant and will not be changing at any point of the injected ply and so for all-composite layers. Further, we use the following parameterized expression of the extensional stiffness tensor \(A_i\) for a laminate with two orientations [23]:
where \(\xi = (\xi _{ij})\in \{0,1\},~i=1,2,\ldots ,N,~j=1,2\) is a binary matrix of orientations and satisfies the following identity:
Besides the definition of the proportion of plies of each fiber orientation within the composite laminate [23]:
By inserting the parameterized expression of \(A_i\) into Eq. (2), with \(A_0 = A_i\), we get
where h is the total thickness of the composite plate. In our model, the proportion of plies \(p_j\) is constant, since we have only two plies. Then we can put \(p_j=1/2,~~j=1,2\). Thus,
Similarly, the bending stiffness tensor D can be represented as
The above results prove that the physical properties of an orthotropic composite plate are invariant and independent of the space coordinates.
Appendix A3
The constants \(a_{ij}, b_{i,j}\), and \(d_{ij},~~i,j=1,2,6,\) are
\([d_{ij}]=-[B_{ij}][b_{ij}]+[D_{ij}]\), \(i,j=1,2\) where \(a_{66}=1/A_{66}\),
\(b_{66}=B_{66}/A_{66}\),
\(d_{66}=(A_{66}D_{66}-B_{66}^{2})/A_{66}\),
\(\Delta = A_{11}A_{22}-A_{12}^{2}\).
Appendix A4
The first-order approximation solutions of Eq. (55) take the following forms:
Similarly, the reader can deduce the factor \( q_{22}\).
Appendix A5
Rights and permissions
About this article
Cite this article
Elmorabie, K.M., Zakaria, K. & Sirwah, M.A. Nonlinear instability of a thin laminated composite circular plate subjected to a tensile periodic load. Acta Mech 231, 5213–5238 (2020). https://doi.org/10.1007/s00707-020-02783-8
Received:
Revised:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00707-020-02783-8