Non-Linear Vibrations of Tow Placed Variable Stiffness Composite Laminates

Total Page:16

File Type:pdf, Size:1020Kb

Non-Linear Vibrations of Tow Placed Variable Stiffness Composite Laminates PROGRAMA DE DOUTORAMENTO EM ENGENHARIA MECÂNICA Non-linear Vibrations of Tow Placed Variable Stiffness Composite Laminates Hamed Akhavan Doctoral Program in Mechanical Engineering Supervisor: Pedro Manuel Leal Ribeiro October 2015 c Hamed Akhavan, 2015 Abstract Vibrations of structures and mechanical systems have plenty of advantages, but also disadvantages. Notably, vibrations can lead to failures or damage of equipment. Therefore, knowing more about the vibratory behaviour helps engineers to control damage and to take preventive measures in order to avoid damages or malfunctions. For a long time, optimising material properties was a goal in material design. An important aspect in the optimisation of fibrous composite materials is the fibre angles in different layers of the laminate. In traditional composite laminates, the fibre orientation is constant in a ply, but some researchers paid attention to the possibility of changing the fibre orientation angle inside a ply. A laminate with plies which have curvilinear fibre paths is considered as a variable stiffness composite laminate (VSCL). Due to the variation of in-plane stiffness in a VSCL, these laminates find special applications in naval, aerospace, or other structures (e.g. in solar array wings, aircraft wings, fuselages and tail sections). Vibrations of VSCL plates are analysed in this thesis. To define a curvilinear fibre path, it is assumed that the fibre angle varies linearly along one of the Cartesian axes, say x, from a value T0 at the centre of the panel to the value T1 at the panel edges, x = ±a=2, where the panel has length a. The fibre path defined in this fashion is called the reference fibre path; other fibre paths in the ply can be produced by shifting - or copying - this path with a constant distance through direction y. During the manufacture process, there is a limitation on the fibre curvature radius; if the fibre path curvature is less than a determined amount, there will be kinking or buckling of the tows in the ply. In this thesis, a p-version finite element based on a third-order shear deformation theory is used to model the displacements. Non-linear terms based on von Kármán strain–displacement relations are included. Using the principle of virtual work, the equations of motion are obtained in matrix form. The order of the equations of motion (regarded as the order of the full model) can be reduced using static condensation and modal coordinates (reduced-order models). The transient responses of the equations of motion of VSCL models are obtained using Newmark method; the same method is employed to obtain solutions in static problems, using time as an auxiliary parameter. The program is implemented in an in-house Visual Fortran code. To find the periodic oscillations of VSCL plates in free and forced vibrations, the shooting method is used, together with a Runge-Kutta-Fehlberg method for integration, again in an in-house Visual Fortran computational code. In a traditional composite laminate the stiffness does not change - at the macro scale - in the laminate, while, in a VSCL laminate stiffnesses change spatially with respect to at least one coordinate. It is found that this change in local stiffness changes the natural frequencies and mode shapes of vibration of the laminate, both in the linear and the non-linear regimes. Using VSCL, the stress field and the strengths in the laminate may change in a different way from what would occur in a laminate with rectilinear fibres. Consequently, the load magnitudes at which damage onset appears, as well as the position of damage onset, can be changed. i ii The effect of geometrical imperfections (out-of-planarity) is included in the p-version finite element too. The free and forced, non-linear, periodic vibrations of VSCL plates are analysed using displacement time-histories, phase plane plots, frequency spectra and shapes assumed by the plates. Bifurcations from main branches of the backbone of frequency response curves are detected and analysed. The effect of imperfection of the plate on the hardening and softening of the free vibration backbone curves are detailed. Resumo Materiais compósitos laminados são frequentemente utilizados em elementos estruturais sujeitos a vibrações. Um aspeto importante na otimização de materiais compósitos reforçados com fibras é a orientação destas últimas. Em laminados compósitos tradicionais, a orientação das fibras dis- postas na mesma lâmina é constante, mas alguns investigadores têm explorado a possibilidade de alterar a orientação das fibras numa lâmina, utilizando fibras curvilíneas. Um laminado com fibras curvilíneas é um laminado compósito com rigidez variável (VSCL do inglês Variable Stiffness Composite Laminates). Vibrações e a deformação em regime estático de placas VSCL são analisadas nesta tese. Considera-se, em particular, o caso de fibras cujo ângulo é uma função bilinear de um dos eixos cartesianos. Limitações impostas à curvatura das fibras pelo processo de fabrico e pela ocorrência de encurvadura são tomadas em consideração. Um aspeto fundamental desta tese, é o desenvolvimento de um elemento finito da versão p para placas VSCL, baseado numa teoria de deformação de corte de terceira ordem. Os termos não-lineares das relações entre deformações e deslocamentos de von Kármán são considerados, para se conseguir estudar grandes deslocamentos (regime geometricamente não-linear). Usando o princípio dos trabalhos virtuais, as equações de movimento são obtidas na forma matricial. O número de equações de movimento do modelo completo foi reduzido utilizando condensação estática e coordenadas modais (modelos de ordem reduzida). No que concerne à resolução das equações do movimento, foram essencialmente utilizados dois métodos. Para obter soluções em regime estático e respostas dinâmicas em regime transiente, as equações do movimento foram resolvidas pelo método de Newmark, no primeiro caso utilizando o tempo como um parâmetro artificial auxiliar. Para encontrar as oscilações periódicas de placas, tanto em vibrações livres como forçadas, foi utilizado o método do disparo (método shooting), em conjunto com um método Runge-Kutta-Fehlberg para a integração no tempo das equações do movimento. Em todos os casos, recorreu-se a códigos computacionais escritos em Fortran 90. Enquanto numa placa laminada compósita com fibras retilíneas a rigidez não muda - à escala macro - de um ponto para o outro do plano médio, numa placa VSCL a rigidez muda espacial- mente. Verificou-se que esta alteração na rigidez local provoca alterações nas frequências naturais e nos modos de vibração do laminado, tanto em regime linear como no regime geometricamente não linear. Usando VSCL, o campo de tensões e a resistência do laminado muda de uma forma diferente do que ocorreria num laminado com fibras retilíneas. Por conseguinte, as magnitudes de carga para a qual dano se inicia, bem como o local em que se dá esse início, são alterados. O efeito de imperfeições geométricas (alteração da planaridade) foi incluído no elemento finito da versão p desenvolvido. Vibrações periódicas nos regimes livre e forçado, não lineares, de pla- cas VSCL foram analisadas utilizando gráficos dos deslocamentos em função do tempo, planos de fase, espectros de frequência e as formas assumidas pelas placas. Bifurcações de ramos prin- cipais da curvas dorsais (backbone curves) e de curvas de resposta de frequência foram detetadas e analisadas. O efeito de imperfeições geométricas da placa sobre a inclinação das curvas dorsais iii iv foi descrito. Acknowledgments I am grateful to Almighty Allah, whose blessings have always been enormous, who gave me the ability to complete this task. I express my heartiest thanks to my parents and family members who encouraged me at every stage of life and saw in me a successful person. It is a great honor and privilege for me to record this deep sense of gratitude and indebtedness to my supervisor Professor Pedro Ribeiro1 who guided me during all the course of this work. He was always there to help me find my way out of any problems. His critical observations helped me improve the analytical skills required for the research process as well as the writing skills required for disseminating the results. I am also thankful to Professor Marcelo Moura1 and to Professor José Rodrigues1 who ex- tended every help to me to have a better understanding of the subject. Here, I am going to appreci- ate the support from Professor Bogdan Epureanu2 who received me very well during a short visit at the University of Michigan, Ann Arbor. I do not forget all the supports from Professor Antonio Ferreira1 during my study in Porto. I appreciate his great helps. Joining the program in University of Porto proved to be a door to make friends from around the world, know new cultures and understand new perspectives of life. I would also like to mention my colleagues and friends whose silent support led me to achieve the goal of completing this work. They are Dr. Julija Vasiljevska, Dr. Rosa Crespo, Dr. Stanislav Stoykov, Dr. Hadi Fanaee, Dr. Marzieh Moeenfard, Dr. Arif Ur Rahman, Dr. Gloria Lema, Dr. Jakov Krstulovic, Dr. David Miranda, Dr. Mariana Seabra, Dr. Mohsen Mirkhalaf, Eng. Tiago Lima, Eng. Romina Fernandes, Eng. Raul Moreira, Eng. Gonçalo Carneiro, and Eng. João Reis. Also, during my stay in Ann Arbor, I had very good time with Dr. Hossein Rokni and his wife. It was a wonderful experience with all of them. During part of this work, I benefited from a research grant of project "Vibrations of variable stiffness composite laminated panels", reference PTDC / EME-PME / 098967 / 2008, a project funded by the Portuguese Science and Technology Foundation (Fundação para a Ciência e a Tec- nologia - FCT3), with national funds (PIDDAC). Also, I would like to mention the support by another grant from FCT with reference number SFRH / BD / 81707 / 2011. These supports are gratefully acknowledged. Finally, but most importantly, I wish to thank my wife Dr.
Recommended publications
  • VIBRATIONS of THIN PLATE with PIEZOELECTRIC ACTUATOR: THEORY and EXPERIMENTS Parikshit Mehta Clemson University, [email protected]
    Clemson University TigerPrints All Theses Theses 12-2009 VIBRATIONS OF THIN PLATE WITH PIEZOELECTRIC ACTUATOR: THEORY AND EXPERIMENTS Parikshit Mehta Clemson University, [email protected] Follow this and additional works at: https://tigerprints.clemson.edu/all_theses Part of the Engineering Mechanics Commons Recommended Citation Mehta, Parikshit, "VIBRATIONS OF THIN PLATE WITH PIEZOELECTRIC ACTUATOR: THEORY AND EXPERIMENTS" (2009). All Theses. 707. https://tigerprints.clemson.edu/all_theses/707 This Thesis is brought to you for free and open access by the Theses at TigerPrints. It has been accepted for inclusion in All Theses by an authorized administrator of TigerPrints. For more information, please contact [email protected]. VIBRATIONS OF THIN PLATE WITH PIEZOELECTRIC ACTUATOR: THEORY AND EXPERIMENTS A Thesis Presented to the Graduate School of Clemson University In Partial Fulfillment of the Requirements for the Degree Master of Science Mechanical Engineering by Parikshit Mehta December 2009 Accepted by: Dr. Nader Jalili, Committee Chair Dr. Mohammed Daqaq Dr. Marian Kennedy ABSTRACT Vibrations of flexible structures have been an important engineering study owing to its both deprecating and complimentary traits. These flexible structures are generally modeled as strings, bars, shafts and beams (one dimensional), membranes and plate (two dimensional) or shell (three dimensional). Structures in many engineering applications, such as building floors, aircraft wings, automobile hoods or pressure vessel end-caps, can be modeled as plates. Undesirable vibrations of any of these engineering structures can lead to catastrophic results. It is important to know the fundamental frequencies of these structures in response to simple or complex excitations or boundary conditions. After their discovery in 1880, piezoelectric materials have made their mark in various engineering applications.
    [Show full text]
  • Buckling and Free Vibrations of Nanoplates – Comparison of Nonlocal Strain and Stress Approaches
    Article Buckling and Free Vibrations of Nanoplates – Comparison of Nonlocal Strain and Stress Approaches Małgorzata Chwał and Aleksander Muc * Institute of Machine Design, Cracow University of Technology, ul. Warszawska 24, 31-155 Kraków, Poland; [email protected] * Correspondence: [email protected] Received: 27 February 2019; Accepted: 1 April 2019; Published: 3 April 2019 Abstract: The buckling and free vibrations of rectangular nanoplates are considered in this paper. The refined continuum transverse shear deformation theory (third and first order) is introduced to formulate the fundamental equations of the nanoplate. In addition, the analysis involves the nonlocal strain and stress theories of elasticity to take into account the small-scale effects encountered in nanostructures/nanocomposites. Hamilton’s principle is used to establish the governing equations of the nanoplate. The Rayleigh-Ritz method is proposed to solve eigenvalue problems dealing with the buckling and free vibration analysis of the nanoplates considered. Some examples are presented to investigate and illustrate the effects of various formulations. Keywords: nonlocal stress theory; nonlocal strain theory; buckling; free vibrations; nanoplates; higher order shear deformation theories 1. Introduction Concerning nano-objects, it is observed that classical (local) continuum theories presenting scale-free relations are inadequate for the mechanical analysis because of the absence of any material length scale parameters, see e.g. works by Peddieson et al. [1], Lu et al. [2], and Barretta et al. [3]. At the nanometer scale, the size effects often become prominent. Some non-classical continuum approaches have been applied to describe the mechanical behavior of small-sized structures, such as couple-stress theory, e.g.
    [Show full text]
  • (TH) Institut Für Baustatik a Stabilized One–Point Integrated Quadrilateral
    Universit¨at Karlsruhe (TH) Institut fur¨ Baustatik A stabilized one–point integrated quadrilateral Reissner–Mindlin plate element F. Gruttmann, W. Wagner Mitteilung 7(2003) BAUSTATIK Universit¨at Karlsruhe (TH) Institut fur¨ Baustatik A stabilized one–point integrated quadrilateral Reissner–Mindlin plate element F. Gruttmann, W. Wagner Mitteilung 7(2003) BAUSTATIK Prof. Dr.–Ing. W. Wagner Telefon: (0721) 608–2280 Institut fur¨ Baustatik Telefax: (0721) 608–6015 c Universit¨at Karlsruhe E–mail: [email protected] Postfach 6980 Internet: http://www.bs.uni-karlsruhe.de 76128 Karlsruhe A stabilized one–point integrated quadrilateral Reissner–Mindlin plate element F. Gruttmann W. Wagner Institut f¨ur Werkstoffe und Mechanik im Bauwesen Institut f¨ur Baustatik Technische Universit¨at Darmstadt Universit¨at Karlsruhe (TH) Alexanderstraße 7 Kaiserstraße 12 64283 Darmstadt 76131 Karlsruhe Germany Germany Contents 1 Introduction 2 2 Basic Equations 3 2.1 Variational Formulation ............................. 3 2.2 Finite Element Equations ............................ 4 3 Examples 8 3.1 Constant bending patch test ........................... 9 3.2 Square plate, test of mesh distorsion ...................... 11 3.2.1 Clamped square plate subjected to a concentrated load ........ 11 3.2.2 Simply supported square plate subjected to uniform load ....... 13 3.3 Square plate subjected to uniform load, test of convergence behaviour .... 14 3.4 Corner supported square plate .......................... 14 3.4.1 Load case 1: uniform load ........................ 14 3.4.2 Load case 2: frequency analysis ..................... 16 3.5 Clamped circular thick plate subjected to a concentrated load ........ 17 4 Conclusions 18 A Appendix, The bending patch test 19 B Appendix, Explicit representation of the one-point integrated stiffness matrix 21 1 Abstract A new quadrilateral Reissner–Mindlin plate element with 12 element degrees of freedom is presented.
    [Show full text]
  • Scherbatiuk Vibration Of.Pdf
    VIBRATION OF SIMPLY SUPPORTED ORTHOTROPIC RECTANGULAR LAMINATED COMPOSITE CROSS-PLY PLATES BY KEVIN DANIEL SCHERBATIUK A thesis Submitted to the Faculty of Graduate Studies In Partial Fulfillment of the Requirements For the Deeree of MASTER OF SCIENCE Department of Civil Engineering University of Manitoba Winnipeg, Manitoba, Canada @ December.200l NationalLibrarY nationale l*l 3få".tg.:" Acquisitions and Acquisitions et Bibiiographicservices servicesbibliographiques 395 Wetlington Str€€t 395, rue Wellington OttawaON KlA0l{4 OttawaON K1AON4 Canacta Canada Your ñl€ Vote úlêrênca Ou frle Noüe rëlércnca The author has granted a non- L'auteur a accordé une licence non exclusive licence allowing the exclusive permettånt à la National Library of Canada to Bibliothèque nationale du Canada de reproduce, loan, distribute or sell reproduire, prêter, distribuer ou copies of this thesis in microform, vendre des copies de cette thèse sous paper or electronic formats. la forme de microfiche/filn, de reproduction sur papier ou sur format électronique. The author retains oumership of the L'auteur conserve la propriété du copynght in this thesis. Neither the droit d'autew qui protège cette thèse. thesis nor substantial exhacts from it Ni la thèse ni des extraits substantiels may be printed or otherwise de celle-ci ne doivent être imprimés reproduced without the author's ou autrement reproduits sâns son permission. aatorisation. 0-612-80016-4 Canad'ä THE UNIVERSITY OF MANITOBA FACULTY OF GRADUATE STUDIES COPYRIG"""*'SSION PAGE Vibration of simply supported
    [Show full text]
  • 6.1 Plate Theory
    Section 6.1 6.1 Plate Theory 6.1.1 Plates A plate is a flat structural element for which the thickness is small compared with the surface dimensions. The thickness is usually constant but may be variable and is measured normal to the middle surface of the plate, Fig. 6.1.1 lateral load middle surface of plate Fig. 6.1.1: A plate 6.1.2 Plate Theory Plates subjected only to in-plane loading can be solved using two-dimensional plane stress theory1 (see Book I, §3.5). On the other hand, plate theory is concerned mainly with lateral loading. One of the differences between plane stress and plate theory is that in the plate theory the stress components are allowed to vary through the thickness of the plate, so that there can be bending moments, Fig. 6.1.2. M Fig. 6.1.2: Stress distribution through the thickness of a plate and resultant bending moment Plate Theory and Beam Theory Plate theory is an approximate theory; assumptions are made and the general three dimensional equations of elasticity are reduced. It is very like the beam theory (see Book 1 although if the in-plane loads are compressive and sufficiently large, they can buckle (see §6.7) Solid Mechanics Part II 120 Kelly Section 6.1 I, §7.4) – only with an extra dimension. It turns out to be an accurate theory provided the plate is relatively thin (as in the beam theory) but also that the deflections are small relative to the thickness. This last point will be discussed further in §6.10.
    [Show full text]
  • Physically and Geometrically Non-Linear Vibrations of Thin Rectangular Plates Ivan Breslavsky, Marco Amabili, Mathias Legrand
    Physically and geometrically non-linear vibrations of thin rectangular plates Ivan Breslavsky, Marco Amabili, Mathias Legrand To cite this version: Ivan Breslavsky, Marco Amabili, Mathias Legrand. Physically and geometrically non-linear vibrations of thin rectangular plates. International Journal of Non-Linear Mechanics, Elsevier, 2013, 58, pp.30-40. 10.1016/j.ijnonlinmec.2013.08.009. hal-00864370 HAL Id: hal-00864370 https://hal.archives-ouvertes.fr/hal-00864370 Submitted on 26 Feb 2015 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Distributed under a Creative Commons Attribution - ShareAlike| 4.0 International License Physically and Geometrically Non-linear Vibrations of Thin Rectangular Plates Ivan Breslavsky, Marco Amabili*, Mathias Legrand Abstract Static deflection as well as free and forced large-amplitude vibrations of thin rectangular rubber plates under uniformly distributed pressure are investigated. Both physical, through a neo-Hookean constitutive law to describe the non-linear elastic deformation of the material, and geometrical non-linearities are accounted for. The deflections of a thin initially flat plate are described by the von Karm´ an´ non-linear plate theory. A method for building a local model, which approximates the plate behavior around a deformed configuration, is proposed.
    [Show full text]
  • Note on Mathematical Development of Plate Theories
    Advanced Studies in Theoretical Physics Vol. 9, 2015, no. 1, 47 - 55 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/astp.2015.411150 Note on Mathematical Development of Plate Theories Patiphan Chantarawichit1, Parichat Kongtong2 Yos Sompornjaroensuk1,* and Jakarin Vibooljak3 1 Department of Civil Engineering Mahanakorn University of Technology Bangkok, 10530, Thailand * Correspondence author 2 Department of Automotive Engineering Thai-Nichi Institute of Technology Bangkok, 10250, Thailand 3 Vibool Choke Ltd., Part Bangkok, 10110, Thailand Copyright © 2014 P. Chantarawichit et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Abstract In this article, the history of mathematical development for plate analysis is reviewed with emphasis on providing the plate’s governing differential equations in rectangular coordinates. Attention is mainly devoted to theories based on the Kirchhoff, Reissner, Mindlin, and Levinson theories for isotropic plates with uniform thickness. Therefore, it is expected that this review paper should be useful for scientists and researchers in assisting them to understand and classify relevant existing plate’s theories quickly. Mathematics Subject Classification: 35Q74, 74K20 Keywords: Partial differential equations, Thin plate, Thick plate, Kirchhoff’s plate, Reissner’s plate, Mindlin’s plate, Levinson’s plate. 48 Patiphan Chantarawichit et
    [Show full text]
  • Information to Users
    INFORMATION TO USERS The most advanced technology has been used to photograph and reproduce this manuscript from the microfilm master. UMI films the text directly from the original or copy submitted. Thus, some thesis and dissertation copies are in typewriter face, while others may be from any type of computer printer. The quality of this reproduction is dependent upon the quality of the copy submitted. Broken or indistinct print, colored or poor quality illustrations and photographs, print bleedthrough, substandard margins, and improper alignment can adversely affect reproduction. In the unlikely event that the author did not send UMI a complete manuscript and there are missing pages, these will be noted. Also, if unauthorized copyright material had to be removed, a note will indicate the deletion. Oversize materials (e.g., maps, drawings, charts) are reproduced by sectioning the original, beginning at the upper left-hand corner and continuing from left to right in equal sections with small overlaps. Each original is also photographed in one exposure and is included in reduced form at the back of the book. Photographs included in the original manuscript have been reproduced xerographically in this copy. Higher quality 6" x 9" black and white photographic prints are available for any photographs or illustrations appearing in this copy for an additional charge. Contact UMI directly to order. University Microfilms International A Bell & Howell Information Company 300 North Zeeb Road, Ann Arbor, Ml 48106-1346 USA 313/761-4700 800/521-0600 Order Number 9105122 Thick plate theories, with applications to vibration Hanna, Nagy Fouad, Ph.D. The Ohio State University, 1990 Copyright ©1991 by Hanna, Nagy Fouad.
    [Show full text]
  • Free and Forced Vibrations of Monolithic and Composite Rectangular Plates with Interior Constrained Points
    Free and Forced Vibrations of Monolithic and Composite Rectangular Plates With Interior Arka P. Chattopadhyay Department of Biomedical Engineering and Mechanics, Constrained Points M/C 0219 Virginia Polytechnic Institute and State University, The restriction of deformations to a subregion of a system undergoing either free or 495 Old Turner Street, forced vibration due to an irregularity or discontinuity in it is called mode localization. Blacksburg, VA 24061 Here, we study mode localization in free and forced vibration of monolithic and unidirec- e-mail: [email protected] tional fiber-reinforced rectangular linearly elastic plates with edges either simply sup- ported (SS) or clamped by using a third-order shear and normal deformable plate theory Romesh C. Batra1 (TSNDT) with points on either one or two normals to the plate midsurface constrained Honorary Fellow from translating in all three directions. The plates studied are symmetric about their mid- Department of Biomedical Engineering surfaces. The in-house developed software based on the finite element method (FEM) is and Mechanics, first verified by comparing predictions from it with either the literature results or those M/C 0219 Virginia Polytechnic computed by using the linear theory of elasticity and the commercial FE software ABAQUS. Institute and State University, New results include: (i) the localization of both in-plane and out-of-plane modes of vibra- 495 Old Turner Street, tion, (ii) increase in the mode localization intensity with an increase in the length/width Blacksburg, VA 24061 ratio of a rectangular plate, (iii) change in the mode localization characteristics with the e-mail: [email protected] fiber orientation angle in unidirectional fiber reinforced laminae, (iv) mode localization due to points on two normals constrained, and (iv) the exchange of energy during forced harmonic vibrations between two regions separated by the line of nearly stationary points that results in a beats-like phenomenon in a subregion of the plate.
    [Show full text]