Numerical Stability Analysis of Composite Plates Thermally Stiffened by Finite Element Method
|
|
- Beryl Lambert
- 6 years ago
- Views:
Transcription
1 Proceedings of COBE 005 Copyright 005 by ABC 8th International Congress of echanical Engineering ovember 6-, 005, Ouro Preto, G umerical Stability Analysis of Composite Plates hermally Stiffened by Finite Element ethod Rodolfo azutti eleiro Instituto ecnológico de Aeronáutica IA São José dos Campos SP, CEP 8-900, Brazil rodomaz@yahoo.com, rodomaz@ita.br José Antônio Hernandes Instituto ecnológico de Aeronáutica IA São José dos Campos SP, CEP 8-900, Brazil hernandes@ita.br Abstract. In the study of laminated plates it is interesting to consider the thermal stiffening effects caused by the difference between cure and operation temperatures. his paper presents a finite element formulation, according to the Reissner-indlin theory, of the buckling problem of laminated composite plates under thermal residual stresses. A geometric stiffness matrix is obtained and applied to several critical load buckling problems of composite plates. he numerical results are compared to results available in the literature. Keywords: finite elements, buckling, laminated, thermal analysis, plate stability. Introduction In the prediction of plate structures behavior the critical buckling load is of vital importance. For anisotropic laminates residual stresses caused by difference of cure and service temperatures may have an effect in the critical buckling load. he thermal residual stresses appear because each lamina has different expansion coefficients in fiber direction and normal to fiber direction. When a plate undergoes a difference of temperature there are two possible behaviors, it can deform freely or be under thermal residual stresses if it is restrained. In symmetric plates thermal residual stress appears only when the plate is constrained or reinforced. he subject of composite plates under residual thermal stress has potential for applications in aeronautical structures. Contrary to the first feeling, laminates can be tailored such that there is an increase of the buckling load for some range of service temperature. his study is performed with a quadratic triangular finite element based on Reissner-indlin theory, proposed by Sze et al. (997) and which has its formulation described in detail by Goto (00), where also its capability was enhanced to deal with membrane behavior.. Formulation and Results Figure. Global coordinate system considered his work is subdivided in three parts: the buckling problem description, finite element geometric stiffness formulation and numerical results... Buckling Problem Description and Formulation he theory of buckling deals principally with conditions under which equilibrium ceases to be stable. A structure is said to be buckled when it loses stability. his can be seen when two or more possible stability positions appear for a certain load, and probably none of them are stable.
2 Buckling of laminated composite plates can happen caused basically by three types of loads: compression in fiber direction, compression normal to fiber direction (in plate plane) and by shear (also in plate plane). Any of this loads and any combination of them will have a critical load, which is the load when the plate loses stability and become buckled. o calculate the critical load we will use the total potential energy criterion (El aschie, 990), which is especially interesting for finite element analysis. his criterion states that the second derivate of total potential energy of the plate with respect to each displacement is zero when the applied load approaches the critical load. π = {}{ ε σ} d + { γ} { τ} d { u} { f } da () A In Eq.(), {ε} is the strain vector, {σ} is the stress vector, {γ} is the out of plane shear strain vector, {τ} is the out of plane shear stress vector, {u} is the displacement vector, {f} is the external forces vector and π is the total potential energy. denotes the volume and A the area of the plate. So, to evaluate the critical load we apply a load case of interest. Appling a multiplier on the applied stresses, with the second derivate of π we can calculate the multiplier factor of the loads, which is the critical load... Finite element formulation for buckling analysis Sze et al. (997) proposed the AS6, a six node triangular element (Figure ) based on Reissner-indlin theory for bending of plates. his element uses quadratic interpolation for bending degrees of freedom (DOF) and linear interpolation for out of plane shear strain. Lucena eto et al. (00) formulated the element explicitly, i.e., using analytical integration and added membrane DOF to it. Figure. AS6 in global coordinates and local coordinates o solve the buckling problem we will consider the pre-buckling stress state. Immediately before the buckling of the plate acting stresses are due to a loading identical to the critical load. In order to determine this stress state two linear independent static stress analyses must be done: one for thermal loading and another for mechanical loading. Knowing the stress state we can compute the geometric stiffness matrices associated to each loading and then finally find the solution of the resulting eigenvalue problem. he linear static finite element analysis development described in Goto s (00) master thesis is omitted here. he following quantities defined in the local plate coordinate system (Figure ) are needed: z k+ 3 ( Aij, Bij, Dij, Eij ) = Qij (, z, z, z ) dz,,,,6 z k = z k+ ( Gij, Fij ) = kij Qij (, z) dz,, 4,5 k = z k k i j = () i j = (3) he matrices A, B, D and E are the stiffness matrices for the in plane behavior, G and F are the out of plane shear stiffness matrices, Q is the layer constitutive matrix, is the number of layers, z k+ is the height at top of layer k and z k is the height at bottom of layer. With these quantities, we can define the generalized loads:
3 Proceedings of COBE 005 Copyright 005 by ABC 8th International Congress of echanical Engineering ovember 6-, 005, Ouro Preto, G { } { } { } [ ] [ ] [ ] [ ] [ ] [ ] A B { ε } B D { m = κ} L D E { Q} { } [ G] [ ] { γ F } = (4) (5) In Eq.(4) and Eq.(5), {}, {}, {L}, {Q} and {} are the stress resultants vectors, {ε m } is the strain vector at the mid-plane of the plate (z=0) and {κ} is the flexural strain vector. For the buckling problem, we must use the nonlinear parts of the strain vectors. So, defining the strains vectors with {ε L } and {γ L } as the linear parts, and {ε } and {γ } as the nonlinear parts (erri, 000): L { ε} = { ε } + { ε } L { γ } = { γ } + { γ } (6) u ' x θ y' x ε = εm + z κ = v' y + z θx' y u + v θ θ L { } { } { } { } ' y ' x y' y x' x ( u ( ) ( )) ' x+ v ' x+ w ' x+ z u ' xθy ' x v ' xθx ' x + z θy ' x+ θx ' x u ( ) ( ) ' y v ' y w ' y z u ' y y ' y v ' y x ' y z y ' y x ' y ( θ θ θ θ ) ( θ θ θ θ ) ε = ( θ θ + θ + θ ) u u + v v + w v + z + + z u + u v v L { γ } { γ } ' x ' y ' x ' y ' x ' y y' y y' x x' y x' x ' x y' y ' y y' x ' x x' y ' y x' x w' x + θ y = w' y θx ( ) ( ) u' xθy v' xθx z θy' xθy θx' xθ + + x = u' yθy v' yθx + z θy' yθy + θx' yθx (7) (8) (9) (0) With: i ' j i = j () ore details on these definitions and plate theory can be found in Reddy (945)... Linear initial pre-buckling thermal stress analysis he thermal load stress resultants vector are given by
4 k+ ({ },{ } ) = [ ]{ α} (, ) z () Q z dz k= z k he thermal load vectors are defined by Eq.(3), where the subscript m and b refer to membrane and bending displacement degrees of freedom. [B m ] and [B b ] are matrices relating the strain vectors with the nodal DOF vectors by Eq.(4), according to Goto (00), {α} is the thermal expansion coefficients vector and is the difference between work and cure temperature. { f } = [ B ] dxdy{ },{ f } = [ B ] dxdy{ } (3) m m b b A A { ε } = [ B ]{ d }, { } [ B ]{ d } m m m κ = (4) b Assembling the system equations of equilibrium, we solve the linear system: b [ K]{ d} = { f } (5) atrix [K] is the linear stiffness of the plate, defined by Eq.(6) from Eq.(), while {d} contains the structure nodal displacements. { d} [ K]{ d} = { ε} { σ} d + { γ} { τ} d (6) With the displacement vector we can calculate the deformation at any point inside an element by Eq.(7) and Eq.(4), and from the deformation with Eq.(4) and Eq.(5) we get the stress resultants and the difference between these resultants (no subscript) and the thermal resultants (with subscript ), are the following thermal residual stress resultants (subscript R). { LR} { QR} { } { L} { Q} { } { R} { } { } =, = (7) { R} { } { }.. Linear initial pre-buckling applied loading analysis R he pre-buckling stress analysis for the applied loading vector, {f P }, is carried out with the solution of the linear system given by Eq.(5). his analysis is distinct from the thermal initial stress analysis because we need separate stress states. Again, with the displacement vector {d}, from Eq.(7), Eq.(4) and Eq.(4) we calculate the stress resultants for the pre-buckling state, denoted by the subscript P, which are { P },{ P },{L P },{Q P } and { P }...3 Buckling eigenvalue problem With all the pre-buckling information available we can now formulate the linear buckling equation. he total potential energy for buckling is given by π = ε σ + σ + γ τ + τ {} ({ P} { } ) d {}({ P} { } ) d (8) his equation has only initial pre-buckling stresses. From Eq.(6) into Eq.(8), and with (9) we have:
5 Proceedings of COBE 005 Copyright 005 by ABC 8th International Congress of echanical Engineering ovember 6-, 005, Ouro Preto, G { d}[ KG ]{ d} = { ε } { σ} d + { γ } {} τ d (9) ({ d} [ K]{ d} { d} [ KG ]{ d} { d} [ KGP ]{ d} ) π = + + (0) Where [K G ] is the geometric stiffness matrix, which is composed of two parts, associated to thermal and mechanical loading, according to the formulation presented by Almeida and Hansen (997), and also used by erri (000). A geometric stiffness matrix is said to be inconsistent if the interpolation functions used to generate it are different from the ones used for the linear stiffness. For AS6, to avoid mesh locking effect (Reddy, 993), the shear deformations use linear interpolation functions for the linear stiffness (Sze et al., 997), which results from a quadratic interpolation where the mid side nodes displacements are average of edge displacements. Here we use quadratic interpolation functions for all the terms in the geometric stiffness matrix. his is done because the mesh locking effect does not occur with the non-linear deformation terms that appear in the geometric stiffness matrix. Defining the vector {e g }, for the AS6: { eg} { u' x u' y v' x v' y w' x w' y θ y' x θy' y θx' x θx' y θy θx} = () { eg} = [ ] [ ] [ 3] [ 4] [ 5] [ 6] { d} = [ ]{ d} () [ ] i nix ' niy ' nix ' n ' iy n n = ix ' iy ' nix ' niy ' 0 ni nix ' niy ' 0 0 ni 0 (3) In Eq.(3) n i are the element displacement interpolation functions. Recalling the definition of the stress resultants, (Eq.() to Eq.(4)), we define [Ψ]: Qx Q 0 x Q y Qy Ψ = L L 0 0 x L L 0 0 y L L 0 x L L 0 y Qx Qx x y Qy Qy x y 0 0 [ ] (4) Finally we obtain
6 ( ) [ ] = [ ] [ Ψ][ ] K dxdy G (5) A Eq.(5) can be formulated for any initial stress state which is explicit in the stress resultants appearing in [Ψ], Eq.(4). In our case [K G ] is due to thermal stress resultants and [K GP ] is due to stress resultants from mechanical loading. A numerical integration is used for Eq.(5), similarly to what was done for the linear stiffness by Sze et al. (997). he variation of π, Eq.(0), leads to π = det ([ K ] + λ[ KG] + λp[ KGP] ) = 0 { d} (6) In the eigenproblem of Eq.(6), two buckling load factors appear, λ P, associated with the mechanical applied loading and λ, associated with the thermal loading. hree buckling problem can be defined depending on which of these factors are free and which at fixed values. When only λ is free (λ P fixed) we have the thermal buckling problem with initial mechanical stress. When only λ P is free (λ fixed) we have the mechanical buckling problem with thermal initial stress. When none of them is fixed, λ P nor λ, they should be merged in only one eigenvalue and the geometric stiffness matrices should be added. here are as much critical loads and buckling modes as nodal displacement degrees of freedom, but we need the first (minimum module) critical load, corresponding to the minimum module eigenvalue. Finally solving the eigenproblem for the first eigenvalue we find the critical load by multiplying the applied load by the corresponding eigenvalue..3 Analysis and results wo study cases are presented, the first case being for the purpose of finite element geometric stiffness validation. he second case considers a reinforced composite plate under thermal residual stresses and uniform compression load. he material properties are defined in able (), and the plate lamination sequences are illustrated in Figure (3) for case two. able. aterial properties and dimensions considered Property alue for case alue for case Longitudinal modulus of elasticity, E 3.4 GPa 54.5 GPa ransversal modulus of elasticity, E GPa.3 GPa In-plane Poisson's ratio, υ ransversal shear modulus, G 3.9 GPa 3.36 GPa ransversal shear modulus, G GPa 6.98 GPa In-plane shear modulus, G 5.8 GPa 6.98 GPa Longitudinal thermal expansion coefficient, α X 0-6 /ºC ransversal thermal expansion coefficient, α - 3. X 0-6 /ºC Lamina side, a 000 mm 360 mm Base plate lamina thickness, t b mm 0.5 mm Reinforced plate lamina thickness, t r mm Figure 3. Plate geometry and lamination sequences
7 Proceedings of COBE 005 Copyright 005 by ABC 8th International Congress of echanical Engineering ovember 6-, 005, Ouro Preto, G.3. Unidirectional compression buckling case A uniform compression is applied in the X-direction (Figure ), without thermal loading, therefore λ is zero. Once Eq.(6) is solved for the buckling load factor, λ P, its product by the applied distributed load gives the buckling load. he plate is simply supported. he results in able, obtained with a 5x5 finite element mesh (900 elements), show a very good agreement when compared to converged finite element solutions obtained with SC-astran, with 00x00 (0000 elements) finite element mesh (Isoni, 004). able. Buckling load factor Laminate AS6 astran (Isoni) Error [0] % [90] % [45] % [0/90/0] % [90/0/90] % [0/45/0] % [90/45/90] % [45/0/45] % [45/90/45] % [45/-45/45] % [0/90]s % [90/0]s % [0/45]s % [90/45]s % [45/0]s % [45/90]s % [45/-45]s %.3. emperature and compression buckling case his case considers the buckling of reinforced plate due to an uniform distributed load applied to its horizontal edges (Figure 3), which corresponds to the eigenproblem solution for λ P from Eq.(6), with λ fixed and equal to unity. he plate is simply supported on all four edges, and the in plane displacements normal to each edge are constrained, except at the superior edge where the load is applied. o calculate the thermal residual stresses the plate is considered unconstrained, but to prevent rigid body motion some displacement degrees of freedom must be properly constrained.,5 ormalized Buckling Loads.,0,5,0 Reinforcement 5mm 0,5 Reinforcement 30mm Reinforcement 45mm 0, emperature Difference (ºC) Figure 4. ormalized critical load by temperature difference
8 hree different reinforcement widths, b, are used, 5, 30 and 45mm, and θ is 90º. he results shown in Figure (4) are the buckling load factors normalized with respect to buckling load factor of the uniform thickness plate (without reinforcement). he thickness t b is scaled for each value of b to maintain the plate volume the same as the uniform thickness plate. herefore, values greater than one indicate buckling load factors greater than that of the uniform plate. he temperature difference in X-axis of Figure (4) is in fact negative (- ). 3. Acknowledgements he authors would like to thank CAPES for the financial support for the first author s master degree scholarship, which generates this work. 4. References Almeida, S. F.., Hansen, J. S.,997, Enhanced Elastic Buckling Loads of Composite Plates With ailored hermal Residual Stresses, Journal of Applied echanics, ol. 64, pp El aschie,. S.,990, Stress, Stability, and Chaos in Structural Engineering: an Energy Approach, Ed. cgraw- Hill Goto, S..,00, Um elemento triangular plano para placas e cascas laminadas, aster hesis, São José dos Campos:IPE Lucena eto, E., Goto, S.., e Kataoka Filho,., 00, "Um Elemento riangular Eficiente para Placas Laminadas", Campinas, Brazil, nd Iberian Latin-American Congress on Computational ethods in Engineering, pp. 8 Paiva, G. I. S.,004, Flambagem de Placas Laminadas Simplesmente Apoiadas Utilizando a eoria Clássica e o étodo de Rayleigh Ritz, aster Degree hesis, São José dos Campos:IA Reddy, J..,945, echanics of Laminated Composite Plates: heory and Analysis, Ed. CRC Press Reddy, J..,993, An Introduction to the Finite Element ethod Second Edition, Ed. cgraw-hill Sze, K. Y., Zhu, D., Chen, D.,997, Quadratic riangular C 0 Plate Bending Element, International Journal for umerical ethods in Engineering, ol. 40, pp erri, A. L.,000, Analise da Flambagem por Cisalhamento de Placas Laminadas na Presença de ensões érmicas, aster Degree hesis, São José dos Campos:IA 5. Responsibility notice he author(s) is (are) the only responsible for the printed material included in this paper.
Bending of Simply Supported Isotropic and Composite Laminate Plates
Bending of Simply Supported Isotropic and Composite Laminate Plates Ernesto Gutierrez-Miravete 1 Isotropic Plates Consider simply a supported rectangular plate of isotropic material (length a, width b,
More informationGEOMETRIC NONLINEAR ANALYSIS
GEOMETRIC NONLINEAR ANALYSIS The approach for solving problems with geometric nonlinearity is presented. The ESAComp solution relies on Elmer open-source computational tool [1] for multiphysics problems.
More informationThermal buckling and post-buckling of laminated composite plates with. temperature dependent properties by an asymptotic numerical method
hermal buckling and post-buckling of laminated composite plates with temperature dependent properties by an asymptotic numerical method F. Abdoun a,*, L. Azrar a,b, E.M. Daya c a LAMA, Higher School of
More informationAdaptation of the Lanczos Algorithm for the Solution of Buckling Eigenvalue Problems
Original Article Adaptation of the Lanczos Algorithm for the Solution of Buckling Eigenvalue Problems Abstract An adaptation of the conventional Lanczos algorithm is proposed to solve the general symmetric
More informationUNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES
UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES A Thesis by WOORAM KIM Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationNonlinear bending analysis of laminated composite stiffened plates
Nonlinear bending analysis of laminated composite stiffened plates * S.N.Patel 1) 1) Dept. of Civi Engineering, BITS Pilani, Pilani Campus, Pilani-333031, (Raj), India 1) shuvendu@pilani.bits-pilani.ac.in
More informationQuintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationStress analysis of a stepped bar
Stress analysis of a stepped bar Problem Find the stresses induced in the axially loaded stepped bar shown in Figure. The bar has cross-sectional areas of A ) and A ) over the lengths l ) and l ), respectively.
More informationUsing MATLAB and. Abaqus. Finite Element Analysis. Introduction to. Amar Khennane. Taylor & Francis Croup. Taylor & Francis Croup,
Introduction to Finite Element Analysis Using MATLAB and Abaqus Amar Khennane Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Croup, an informa business
More informationENGN 2340 Final Project Report. Optimization of Mechanical Isotropy of Soft Network Material
ENGN 2340 Final Project Report Optimization of Mechanical Isotropy of Soft Network Material Enrui Zhang 12/15/2017 1. Introduction of the Problem This project deals with the stress-strain response of a
More informationGeneral elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationOPTIMIZATION OF LAMINATED PLATES SUBJECTED TO IMPULSE LOADS
OPTIMIZATION OF LAMINATED PLATES SUBJECTED TO IMPULSE LOADS Rafael Thiago Luiz Ferreira rthiago@ita.br José Antônio Hernandes hernandes@ita.br Eliseu Lucena Neto eliseu@ita.br Paulo Ivo Braga de Queiroz
More informationInternational Journal of Advanced Engineering Technology E-ISSN
Research Article INTEGRATED FORCE METHOD FOR FIBER REINFORCED COMPOSITE PLATE BENDING PROBLEMS Doiphode G. S., Patodi S. C.* Address for Correspondence Assistant Professor, Applied Mechanics Department,
More informationFig. 1. Circular fiber and interphase between the fiber and the matrix.
Finite element unit cell model based on ABAQUS for fiber reinforced composites Tian Tang Composites Manufacturing & Simulation Center, Purdue University West Lafayette, IN 47906 1. Problem Statement In
More informationThe stiffness of plates
The stiffness of plates 1. Introduction The word plate is a collective term for elements in which forces can be transferred in two directions. Floors, walls, bridge slabs and laminates are all plates.
More informationEML4507 Finite Element Analysis and Design EXAM 1
2-17-15 Name (underline last name): EML4507 Finite Element Analysis and Design EXAM 1 In this exam you may not use any materials except a pencil or a pen, an 8.5x11 formula sheet, and a calculator. Whenever
More informationSANDWICH COMPOSITE BEAMS for STRUCTURAL APPLICATIONS
SANDWICH COMPOSITE BEAMS for STRUCTURAL APPLICATIONS de Aguiar, José M., josemaguiar@gmail.com Faculdade de Tecnologia de São Paulo, FATEC-SP Centro Estadual de Educação Tecnológica Paula Souza. CEETEPS
More informationFinite element analysis of FRP pipelines dynamic stability
Boundary Elements and Other esh Reduction ethods XXXVIII 39 Finite element analysis of FRP pipelines dynamic stability D. G. Pavlou & A. I. Nergaard Department of echanical and Structural Engineering and
More informationUNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich
UNIVERSITY OF SASKATCHEWAN ME 313.3 MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich A CLOSED BOOK EXAMINATION TIME: 3 HOURS For Marker s Use Only LAST NAME (printed): FIRST
More informationModule-6: Laminated Composites-II. Learning Unit-1: M6.1. M 6.1 Structural Mechanics of Laminates
Module-6: Laminated Composites-II Learning Unit-1: M6.1 M 6.1 Structural Mechanics of Laminates Classical Lamination Theory: Laminate Stiffness Matrix To this point in the development of classical lamination
More informationA HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS
A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A. Kroker, W. Becker TU Darmstadt, Department of Mechanical Engineering, Chair of Structural Mechanics Hochschulstr. 1, D-64289 Darmstadt, Germany kroker@mechanik.tu-darmstadt.de,
More informationTABLE OF CONTENTS. Mechanics of Composite Materials, Second Edition Autar K Kaw University of South Florida, Tampa, USA
Mechanics of Composite Materials, Second Edition Autar K Kaw University of South Florida, Tampa, USA TABLE OF CONTENTS 1. INTRODUCTION TO COMPOSITE MATERIALS 1.1 Introduction... 1.2 Classification... 1.2.1
More informationFLEXURAL RESPONSE OF FIBER RENFORCED PLASTIC DECKS USING HIGHER-ORDER SHEAR DEFORMABLE PLATE THEORY
Asia-Pacific Conference on FRP in Structures (APFIS 2007) S.T. Smith (ed) 2007 International Institute for FRP in Construction FLEXURAL RESPONSE OF FIBER RENFORCED PLASTIC DECKS USING HIGHER-ORDER SHEAR
More informationEsben Byskov. Elementary Continuum. Mechanics for Everyone. With Applications to Structural Mechanics. Springer
Esben Byskov Elementary Continuum Mechanics for Everyone With Applications to Structural Mechanics Springer Contents Preface v Contents ix Introduction What Is Continuum Mechanics? "I Need Continuum Mechanics
More informationGeneric Strategies to Implement Material Grading in Finite Element Methods for Isotropic and Anisotropic Materials
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Engineering Mechanics Dissertations & Theses Mechanical & Materials Engineering, Department of Winter 12-9-2011 Generic
More informationAERSYS KNOWLEDGE UNIT
-7016 1. INTRODUCTION The scope of this document is to provide a clarification and a deeper understanding of the two different ways to move the mid plane of the element out of the nodal plane. Although
More informationModule III - Macro-mechanics of Lamina. Lecture 23. Macro-Mechanics of Lamina
Module III - Macro-mechanics of Lamina Lecture 23 Macro-Mechanics of Lamina For better understanding of the macromechanics of lamina, the knowledge of the material properties in essential. Therefore, the
More informationDynamic Response Of Laminated Composite Shells Subjected To Impulsive Loads
IOSR Journal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 2278-1684,p-ISSN: 2320-334X, Volume 14, Issue 3 Ver. I (May. - June. 2017), PP 108-123 www.iosrjournals.org Dynamic Response Of Laminated
More informationDYNAMIC FAILURE ANALYSIS OF LAMINATED COMPOSITE PLATES
Association of Metallurgical Engineers of Serbia AMES Scientific paper UDC:669.1-419:628.183=20 DYNAMIC FAILURE ANALYSIS OF LAMINATED COMPOSITE PLATES J. ESKANDARI JAM 1 and N. GARSHASBI NIA 2 1- Aerospace
More informationDynamic Analysis of Laminated Composite Plate Structure with Square Cut-Out under Hygrothermal Load
Dynamic Analysis of Laminated Composite Plate Structure with Square Cut-Out under Hygrothermal Load Arun Mukherjee 1, Dr. Sreyashi Das (nee Pal) 2 and Dr. A. Guha Niyogi 3 1 PG student, 2 Asst. Professor,
More informationCOPYRIGHTED MATERIAL. Index
Index A Admissible function, 163 Amplification factor, 36 Amplitude, 1, 22 Amplitude-modulated carrier, 630 Amplitude ratio, 36 Antinodes, 612 Approximate analytical methods, 647 Assumed modes method,
More informationPresented By: EAS 6939 Aerospace Structural Composites
A Beam Theory for Laminated Composites and Application to Torsion Problems Dr. BhavaniV. Sankar Presented By: Sameer Luthra EAS 6939 Aerospace Structural Composites 1 Introduction Composite beams have
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 13
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras (Refer Slide Time: 00:25) Module - 01 Lecture - 13 In the last class, we have seen how
More informationHIGHER-ORDER THEORIES
HIGHER-ORDER THEORIES Third-order Shear Deformation Plate Theory Displacement and strain fields Equations of motion Navier s solution for bending Layerwise Laminate Theory Interlaminar stress and strain
More informationINTERNATIONAL JOURNAL OF APPLIED ENGINEERING RESEARCH, DINDIGUL Volume 2, No 1, 2011
Interlaminar failure analysis of FRP cross ply laminate with elliptical cutout Venkateswara Rao.S 1, Sd. Abdul Kalam 1, Srilakshmi.S 1, Bala Krishna Murthy.V 2 1 Mechanical Engineering Department, P. V.
More informationMaterials and Structures. Indian Institute of Technology Kanpur
Introduction to Composite Materials and Structures Nachiketa Tiwari Indian Institute of Technology Kanpur Lecture 34 Thermal Stresses in Plates Lecture Overview Introduction Mechanical and Thermal Strains
More informationMechanical Behavior of Composite Tapered Lamina
International Journal of Engineering Research and Development e-issn: 2278-067X, p-issn: 2278-800X, www.ijerd.com Volume 10, Issue 8 (August 2014), PP.19-27 Mechanical Behavior of Composite Tapered Lamina
More informationPassive Damping Characteristics of Carbon Epoxy Composite Plates
Journal of aterials Science and Engineering A 6 (1-2) (2016) 35-42 doi: 10.17265/2161-6213/2016.1-2.005 D DAVID PUBLISHIG Passive Damping Characteristics of Carbon Epoxy Composite Plates Dileep Kumar K
More informationEffect of Thermal Stresses on the Failure Criteria of Fiber Composites
Effect of Thermal Stresses on the Failure Criteria of Fiber Composites Martin Leong * Institute of Mechanical Engineering Aalborg University, Aalborg, Denmark Bhavani V. Sankar Department of Mechanical
More informationPrediction of Elastic Constants on 3D Four-directional Braided
Prediction of Elastic Constants on 3D Four-directional Braided Composites Prediction of Elastic Constants on 3D Four-directional Braided Composites Liang Dao Zhou 1,2,* and Zhuo Zhuang 1 1 School of Aerospace,
More informationInstitute of Structural Engineering Page 1. Method of Finite Elements I. Chapter 2. The Direct Stiffness Method. Method of Finite Elements I
Institute of Structural Engineering Page 1 Chapter 2 The Direct Stiffness Method Institute of Structural Engineering Page 2 Direct Stiffness Method (DSM) Computational method for structural analysis Matrix
More informationAEROELASTIC ANALYSIS OF SPHERICAL SHELLS
11th World Congress on Computational Mechanics (WCCM XI) 5th European Conference on Computational Mechanics (ECCM V) 6th European Conference on Computational Fluid Dynamics (ECFD VI) E. Oñate, J. Oliver
More informationChapter 12 Plate Bending Elements. Chapter 12 Plate Bending Elements
CIVL 7/8117 Chapter 12 - Plate Bending Elements 1/34 Chapter 12 Plate Bending Elements Learning Objectives To introduce basic concepts of plate bending. To derive a common plate bending element stiffness
More informationCRITERIA FOR SELECTION OF FEM MODELS.
CRITERIA FOR SELECTION OF FEM MODELS. Prof. P. C.Vasani,Applied Mechanics Department, L. D. College of Engineering,Ahmedabad- 380015 Ph.(079) 7486320 [R] E-mail:pcv-im@eth.net 1. Criteria for Convergence.
More informationModeling of the Bending Stiffness of a Bimaterial Beam by the Approximation of One-Dimensional of Laminated Theory
. Flores-Domínguez Int. Journal of Engineering Research and Applications RESEARCH ARTICLE OPEN ACCESS odeling of the Bending Stiffness of a Bimaterial Beam by the Approimation of One-Dimensional of Laminated
More informationFree Vibration Response of a Multilayer Smart Hybrid Composite Plate with Embedded SMA Wires
11(2014) 279-298 Free Vibration Response of a Multilayer Smart Hybrid Composite Plate with Embedded SMA Wires Abstract In this paper, free vibration response of a hybrid composite plate was studied. Effects
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 11
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Module - 01 Lecture - 11 Last class, what we did is, we looked at a method called superposition
More informationComputational Analysis for Composites
Computational Analysis for Composites Professor Johann Sienz and Dr. Tony Murmu Swansea University July, 011 The topics covered include: OUTLINE Overview of composites and their applications Micromechanics
More informationFREE VIBRATION OF AXIALLY LOADED FUNCTIONALLY GRADED SANDWICH BEAMS USING REFINED SHEAR DEFORMATION THEORY
FREE VIBRATION OF AXIALLY LOADED FUNCTIONALLY GRADED SANDWICH BEAMS USING REFINED SHEAR DEFORMATION THEORY Thuc P. Vo 1, Adelaja Israel Osofero 1, Marco Corradi 1, Fawad Inam 1 1 Faculty of Engineering
More informationInfluence of residual stresses in the structural behavior of. tubular columns and arches. Nuno Rocha Cima Gomes
October 2014 Influence of residual stresses in the structural behavior of Abstract tubular columns and arches Nuno Rocha Cima Gomes Instituto Superior Técnico, Universidade de Lisboa, Portugal Contact:
More informationHIGHER-ORDER THEORIES
HIGHER-ORDER THEORIES THIRD-ORDER SHEAR DEFORMATION PLATE THEORY LAYERWISE LAMINATE THEORY J.N. Reddy 1 Third-Order Shear Deformation Plate Theory Assumed Displacement Field µ u(x y z t) u 0 (x y t) +
More informationVORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS
The 4 th World Conference on Earthquake Engineering October -7, 008, Beijing, China VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS K.
More informationDepartment of Structural, Faculty of Civil Engineering, Architecture and Urban Design, State University of Campinas, Brazil
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm A SIMPLIFIED FORMULATION FOR STRESS AND TRACTION BOUNDARY IN- TEGRAL EQUATIONS USING THE
More informationD : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown.
D : SOLID MECHANICS Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown. Q.2 Consider the forces of magnitude F acting on the sides of the regular hexagon having
More informationTransactions of the VŠB Technical University of Ostrava, Mechanical Series. article No Roland JANČO *
Transactions of the VŠB Technical University of Ostrava, Mechanical Series No. 1, 013, vol. LIX article No. 1930 Roland JANČO * NUMERICAL AND EXACT SOLUTION OF BUCKLING LOAD FOR BEAM ON ELASTIC FOUNDATION
More informationUnit 18 Other Issues In Buckling/Structural Instability
Unit 18 Other Issues In Buckling/Structural Instability Readings: Rivello Timoshenko Jones 14.3, 14.5, 14.6, 14.7 (read these at least, others at your leisure ) Ch. 15, Ch. 16 Theory of Elastic Stability
More informationBHAR AT HID AS AN ENGIN E ERI N G C O L L E G E NATTR A MPA LL I
BHAR AT HID AS AN ENGIN E ERI N G C O L L E G E NATTR A MPA LL I 635 8 54. Third Year M E C H A NICAL VI S E M ES TER QUE S T I ON B ANK Subject: ME 6 603 FIN I T E E LE ME N T A N A L YSIS UNI T - I INTRODUCTION
More informationRODS: THERMAL STRESS AND STRESS CONCENTRATION
RODS: HERML SRESS ND SRESS CONCENRION Example 5 rod of length L, cross-sectional area, and modulus of elasticity E, has been placed inside a tube of the same length L, but of cross-sectional area and modulus
More informationTHERMAL STRESS ANALYSIS OF A FIBER-EPOXY COMPOSITE MATERIAL
HERMAL SCIENCE, Year 2011, Vol. 15, No. 2, pp. 559-563 559 HERMAL SRESS ANALYSIS OF A FIBER-EPOXY COMPOSIE MAERIAL by Radoljub P. OMI] a, Aleksandar S. SEDMAK b*, Dobrivoje M. ]AI] c, Marko V. MILOŠ b,
More informationFinite Element Analysis of Graphite/Epoxy Composite Pressure Vessel
Journal of Materials Science and Chemical Engineering, 2017, 5, 19-28 http://www.scirp.org/journal/msce ISSN Online: 2327-6053 ISSN Print: 2327-6045 Finite Element Analysis of Graphite/Epoxy Composite
More informationCOROTATIONAL NONLINEAR DYNAMIC ANALYSIS OF LAMINATED COMPOSITE SHELL STRUCTURES
Copyright c 009 by ABCM January 04-08, 010, Foz do Iguaçu, PR, Brazil COROTATIONAL NONLINEAR DYNAMIC ANALYSIS OF LAMINATED COMPOSITE SHELL STRUCTURES Felipe Schaedler de Almeida, schaedleralmeida@gmail.com
More informationQUESTION BANK Composite Materials
QUESTION BANK Composite Materials 1. Define composite material. 2. What is the need for composite material? 3. Mention important characterits of composite material 4. Give examples for fiber material 5.
More informationKINK BAND FORMATION OF FIBER REINFORCED POLYMER (FRP)
KINK BAND FORMATION OF FIBER REINFORCED POLYMER (FRP) 1 University of Science & Technology Beijing, China, niukm@ustb.edu.cn 2 Tsinghua University, Department of Engineering Mechanics, Beijing, China,
More informationVibration Behaviour of Laminated Composite Flat Panel Under Hygrothermal Environment
International Review of Applied Engineering Research. ISSN 2248-9967 Volume 4, Number 5 (2014), pp. 455-464 Research India Publications http://www.ripublication.com/iraer.htm Vibration Behaviour of Laminated
More informationStructural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian
Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian ahmadian@iust.ac.ir Dynamic Response of MDOF Systems: Mode-Superposition Method Mode-Superposition Method:
More informationFLEXIBILITY METHOD FOR INDETERMINATE FRAMES
UNIT - I FLEXIBILITY METHOD FOR INDETERMINATE FRAMES 1. What is meant by indeterminate structures? Structures that do not satisfy the conditions of equilibrium are called indeterminate structure. These
More informationCHAPTER 14 BUCKLING ANALYSIS OF 1D AND 2D STRUCTURES
CHAPTER 14 BUCKLING ANALYSIS OF 1D AND 2D STRUCTURES 14.1 GENERAL REMARKS In structures where dominant loading is usually static, the most common cause of the collapse is a buckling failure. Buckling may
More informationLecture 8. Stress Strain in Multi-dimension
Lecture 8. Stress Strain in Multi-dimension Module. General Field Equations General Field Equations [] Equilibrium Equations in Elastic bodies xx x y z yx zx f x 0, etc [2] Kinematics xx u x x,etc. [3]
More informationChapter 6 2D Elements Plate Elements
Institute of Structural Engineering Page 1 Chapter 6 2D Elements Plate Elements Method of Finite Elements I Institute of Structural Engineering Page 2 Continuum Elements Plane Stress Plane Strain Toda
More informationShishir Kumar Sahu, Ph. D. Professor Civil Engineering Department National Institute of Technology, Rourkela Orissa, India
Shishir Kumar Sahu, Ph. D. Professor Civil Engineering Department National Institute of Technology, Rourkela Orissa, India Introduction Review Of Literature Aim & Scope Of Present Investigation Finite
More informationNUMERICAL ANALYSIS OF A PILE SUBJECTED TO LATERAL LOADS
IGC 009, Guntur, INDIA NUMERICAL ANALYSIS OF A PILE SUBJECTED TO LATERAL LOADS Mohammed Younus Ahmed Graduate Student, Earthquake Engineering Research Center, IIIT Hyderabad, Gachibowli, Hyderabad 3, India.
More informationINTERNATIONAL JOURNAL OF PURE AND APPLIED RESEARCH IN ENGINEERING AND TECHNOLOGY
INTERNATIONAL JOURNAL OF PURE AND APPLIED RESEARCH IN ENGINEERING AND TECHNOLOGY A PATH FOR HORIZING YOUR INNOVATIVE WORK SPECIAL ISSUE FOR INTERNATIONAL LEVEL CONFERENCE "ADVANCES IN SCIENCE, TECHNOLOGY
More informationMODIFIED HYPERBOLIC SHEAR DEFORMATION THEORY FOR STATIC FLEXURE ANALYSIS OF THICK ISOTROPIC BEAM
MODIFIED HYPERBOLIC SHEAR DEFORMATION THEORY FOR STATIC FLEXURE ANALYSIS OF THICK ISOTROPIC BEAM S. Jasotharan * and I.R.A. Weerasekera University of Moratuwa, Moratuwa, Sri Lanka * E-Mail: jasos91@hotmail.com,
More informationChapter 5 CENTRIC TENSION OR COMPRESSION ( AXIAL LOADING )
Chapter 5 CENTRIC TENSION OR COMPRESSION ( AXIAL LOADING ) 5.1 DEFINITION A construction member is subjected to centric (axial) tension or compression if in any cross section the single distinct stress
More informationPassive Damping Characteristics of Carbon Epoxy Composite Plates
Journal of Materials Science and Engineering A 6 (-) 35-4 doi:.765/6-63/6.-.5 D DAVID PUBLISHING Passive Damping Characteristics of Carbon Epoxy Composite Plates Dileep Kumar K * and V V Subba Rao Faculty
More informationFINITE ELEMENT ANALYSIS OF COMPOSITE MATERIALS
FINITE ELEMENT ANALYSIS OF COMPOSITE MATERIALS Ever J. Barbero Department of Mechanical and Aerospace Engineering West Virginia University USA CRC Press Taylor &.Francis Group Boca Raton London New York
More informationBUCKLING MODE CLASSIFICATION OF MEMBERS WITH OPEN THIN-WALLED CROSS-SECTIONS
CIMS 4 Fourth International Conference on Coupled Instabilities in Metal Structures Rome, Italy, 27-29 September, 24 BUCKLING MODE CLASSIFICATION OF MEMBERS WITH OPEN THIN-WALLED CROSS-SECTIONS S. ÁDÁNY,
More informationPractice Final Examination. Please initial the statement below to show that you have read it
EN175: Advanced Mechanics of Solids Practice Final Examination School of Engineering Brown University NAME: General Instructions No collaboration of any kind is permitted on this examination. You may use
More informationCOMPARISON OF PLATE MODELS FOR ANALYSIS OF LAMINATED COMPOSITES
COMPARISON OF PLATE MODELS FOR ANALYSIS OF LAMINATED COMPOSITES P. M. Mohite and C. S. Upadhyay** Department of Aerospace Engineering, IIT Kanpur 0806, INDIA, e-mail: mohite@iitk.ac.in Assistant Professor,
More informationMechanics of Inflatable Fabric Beams
Copyright c 2008 ICCES ICCES, vol.5, no.2, pp.93-98 Mechanics of Inflatable Fabric Beams C. Wielgosz 1,J.C.Thomas 1,A.LeVan 1 Summary In this paper we present a summary of the behaviour of inflatable fabric
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationNONLINEAR STRUCTURAL DYNAMICS USING FE METHODS
NONLINEAR STRUCTURAL DYNAMICS USING FE METHODS Nonlinear Structural Dynamics Using FE Methods emphasizes fundamental mechanics principles and outlines a modern approach to understanding structural dynamics.
More informationURL: < >
Citation: Vo, Thuc, Thai, Huu-Tai and Inam, Fawad (213) Axial-flexural coupled vibration and buckling of composite beams using sinusoidal shear deformation theory. Archive of Applied Mechanics, 83 (4).
More informationComputational Stiffness Method
Computational Stiffness Method Hand calculations are central in the classical stiffness method. In that approach, the stiffness matrix is established column-by-column by setting the degrees of freedom
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 06
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 06 In the last lecture, we have seen a boundary value problem, using the formal
More informationInverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros
Inverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros Computational Design Forward design: direct manipulation of design parameters Level of abstraction Exploration
More informationComposite Structures. Indian Institute of Technology Kanpur
Mechanics of Laminated Composite Structures Nachiketa Tiwari Indian Institute of Technology Kanpur Lecture 23 Analysis of an Orthotropic Ply Lecture Overview Introduction Engineering constants for an 2
More informationAdvanced Vibrations. Distributed-Parameter Systems: Approximate Methods Lecture 20. By: H. Ahmadian
Advanced Vibrations Distributed-Parameter Systems: Approximate Methods Lecture 20 By: H. Ahmadian ahmadian@iust.ac.ir Distributed-Parameter Systems: Approximate Methods Rayleigh's Principle The Rayleigh-Ritz
More informationSimulation of Stress Concentration Problems in Laminated Plates by Quasi-Trefftz Finite Element Models
1677 Simulation of Stress Concentration Problems in Laminated Plates by Quasi-Trefftz Finite Element Models Abstract Hybrid quasi-trefftz finite elements have been applied with success to the analysis
More informationTheoretical Manual Theoretical background to the Strand7 finite element analysis system
Theoretical Manual Theoretical background to the Strand7 finite element analysis system Edition 1 January 2005 Strand7 Release 2.3 2004-2005 Strand7 Pty Limited All rights reserved Contents Preface Chapter
More informationON THE NUMERICAL ANALYSIS OF COMPOSITE MATERIAL
The 4th International Conference Advanced Composite Materials Engineering COMAT 8- October, Brasov, Romania O THE UMERICAL AALYSIS OF COMPOSITE MATERIAL D.D. icoara Transilvania University, City Brasov,
More informationCOMPUTATIONAL MODELING FOR THE COUPLING ROD-ANODE IN ALUMINUM INDUSTRY
COMPUTATIONAL MODELING FOR THE COUPLING ROD-ANODE IN ALUMINUM INDUSTRY (1) Adriano G. Batista; (1) Fábio A. Cunha; (2) M. Sena (1) Departamento de Engenharia Mecânica, Centro Tecnológico, Universidade
More informationAccelerated Testing Methodology for Long Term Durability of CFRP
IFREMER-ONR Workshop on Durability of Composites in a Marine Environment August 23 24, 22 IFREMER Centre, Nantes, France Accelerated esting Methodology for Long erm Durability of CFRP Masayuki Nakada*,
More informationUNIVERSITY OF HAWAII COLLEGE OF ENGINEERING DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING
UNIVERSITY OF HAWAII COLLEGE OF ENGINEERING DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING ACKNOWLEDGMENTS This report consists of the dissertation by Ms. Yan Jane Liu, submitted in partial fulfillment
More informationThermal buckling analysis of shear deformable laminated orthotropic plates by differential quadrature
Steel and Composite Structures, Vol. 12, No. 2 (2012) 129-147 129 hermal buckling analysis of shear deformable laminated orthotropic plates by differential quadrature S. Moradi 1 and Mohammad Hassan Mansouri
More informationD : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each.
GTE 2016 Q. 1 Q. 9 carry one mark each. D : SOLID MECHNICS Q.1 single degree of freedom vibrating system has mass of 5 kg, stiffness of 500 N/m and damping coefficient of 100 N-s/m. To make the system
More informationReview of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis
uke University epartment of Civil and Environmental Engineering CEE 42L. Matrix Structural Analysis Henri P. Gavin Fall, 22 Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods
More informationUncorrected Proof Available online at
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering00 (2017) 000 000 www.elsevier.com/locate/procedia Sustainable Civil Engineering Structures and Construction Materials, SCESCM
More informationDynamic and buckling analysis of FRP portal frames using a locking-free finite element
Fourth International Conference on FRP Composites in Civil Engineering (CICE8) 22-24July 8, Zurich, Switzerland Dynamic and buckling analysis of FRP portal frames using a locking-free finite element F.
More information