Finite element analysis of FRP pipelines dynamic stability
|
|
- Arron Hicks
- 5 years ago
- Views:
Transcription
1 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 aterials Science, University of Stavanger, Norway Abstract A finite element model is derived for flow-induced instability analysis of multilayered, filament wound FRP pipelines resting on equispaced elastic supports. The critical flow velocity causing instability is obtained by taking into account the anisotropic properties of pipes material, the elastic properties of pipeline supports and flow parameters. A representative example is solved and discussed. Keywords: FRP pipeline, flow induced instability, critical flow velocity, FE. Introduction Even though fiber reinforced polymeric (FRP) materials, apart from high strength and low density provide superior behavior than steel in corrosion, reducing thus the maintenance cost, their use for pipeline installations is not as wide as should be expected. However, in the last decade the increasing rate of adopting FRP s for pipelines means that the industry is not far from extended applicability of composite materials for pipeline projects []. Since steel is still the main material used for long pipeline infrastructures, the flow-induced instability has systematically been analyzed only for such isotropic material. Application of the accumulated dynamic analysis know-how developed for steel pipelines on FRP ones is not possible due to their anisotropic behavior. In the present work, finite element (FE) analysis is carried out for critical flow velocities estimation of filament wound, multi-layered FRP pipelines resting on elastic supports. doi:0.495/be380
2 40 Boundary Elements and Other esh Reduction ethods XXXVIII Classical lamination theory o o o Based on Kirchhoff assumption, the correlation of strains,, and x y xy o o o curvatures k, k, k of the middle plane of a laminate (Figure ) with the x y xy corresponding force (per unit length) and moment (per unit length) resultants N x, N y, N, xy x, y, is given by the following equation (e.g. []): xy a a a b b b N 0 x 6 6 x 0 y a a a6 b b b N 6 y 0 xy a6 a6 66 b6 b6 b66 N xy 0 k x b b b6 d d d6 x 0 k y b b b6 d d d 6 y 0 k xy b6 b6 b66 d6 d6 d66 xy () Figure : Layout of an FRP laminate. The components a ij, b ij, d ij of the above matrix can be derived by the following equation: a a a b b b B B B a a a b b b B B B a a b b b B B B b b b d d d B B B D D D b b b6 d d d 6 B B B D D D b6 b6 b66 d6 d6 d66 B B B D D D ()
3 Boundary Elements and Other esh Reduction ethods XXXVIII 4 In equation () the following definitions are valid: N (3) A Q ( z z ) ij ijk k k k B Q ( z z ) N ij ijk k k k (4) D Q ( z z ) N 3 3 ij ijk k k 3 k (5) In the above equations, N is the number of layers composing the laminate, z k represents the distances of any layer k from the end surface of the laminate, and the parameters Q are given for each lamina by the following equations: ij Q Q m Q S n m Q n (6) Q Q Q Q n m Q n m (7) 66 Q Q Q Q nm Q Q Q n m (8) Q Q n Q Q n m Q m (9) Q Q Q Q n m Q Q Q nm (0) Q Q Q Q S n m Q n m () The parameters m, n depends on the fiber orientation angle ϑ of each lamina (Figure ): m cos () n sin (3) The parameters Qij are correlated to material properties in the principal coordinate system x x (Figure ) of each individual lamina composing the laminate: Q E (4) vv
4 4 Boundary Elements and Other esh Reduction ethods XXXVIII Q ve ve v v v v Q E (5) (6) vv Q G (7) 66 E (8) E E, E are the elasticity moduli, G is the shear modulus and ν, ν are the Poisson s ratios in the principal coordinate system x -x. 3 otion equation of multi-layered filament wound FRP pipes According to [3], the motion equation of filament wound FRP pipe is given by: 4 wxt (, ) wxt (, ) wxt (, ) wxt (, ) Seq U U ( m) 0 4 x x xt t (9), m is the mass per unit length for the liquid and the pipe s material respectively, U is the mean velocity of the liquid, w(x,t) is the elastic deflection of the pipe and S eq is the equivalent bending stiffness of the pipe for bending moment calculation: The parameter dw ( x) S (0) eq dx S can be obtained from [3]: eq S eq D D a d () In the above equations (x) is the bending moment, and D is the diameter of the pipe. Substituting (e.g. [4]) into equation (9) the solution () the following equation can be obtained: wxt (, ) Re uxe ( ) it, () 4 dux ( ) dux ( ) dux ( ) Seq U Ui ( m) u( x) 0 (3) 4 dx dx dx
5 Boundary Elements and Other esh Reduction ethods XXXVIII Element equation of a pipe segment Following a standard mathematical procedure the differential equation (3) can be transformed to the following matrix differential equation of first order: and In equation (4) the following is valid: d Y. Y (4) dx Y u x u x u x u x (5) ( ) ( ) ( ) ( ) T ( m) Ui U 0 Seq Seq S eq (6) Therefore, the solution of the matrix equation (4) can be written: L is the length of the pipe segment. [ ] Y e L. Y (7) Y u L u L u L u L ( ) ( ) ( ) ( ) T Y u u u u (0) (0) (0) (0) T (8) (9) Taking into account the following relations of solid mechanics: ( x) u( x) (30) m( x) S u ( x) (3) eq f( x) S u ( x) (3) where ( x), mx ( ), f( x) is the slope, bending moment and shear force distribution, the matrix equation (7) can now be written in the following form: eq
6 44 Boundary Elements and Other esh Reduction ethods XXXVIII u u L e m / S m / S eq eq f / S eq f / S eq (33) Let s assume that the matrix exponential matrix: [ ] L e is expressed by the following e L c c c c c c c c c c c c c c c c (34) Performing some rearrangements in eq. (33), the following expression can be obtained: f u m f u F U m (35) In equation (35) the following is valid: F U c4 / Seq c3 / Seq 0 0 c / S c / S eq 3 eq c34 / Seq c33 / Seq 0 c44 / Seq c43 / Seq 0 c c 0 c c 0 Seqc3 Seqc3 0 0 Seqc4 Seqc4 0 0 Therefore, the following equation for a segment - can be derived: f u m k f u m (36) (37) (38)
7 Boundary Elements and Other esh Reduction ethods XXXVIII 45 In equation (38) we have: k F U (39) 4 Pipeline resting on elastic supports Let s assume a pipeline resting on equispaced elastic supports (Figure ). The elastic supports can be linear springs with elastic constant k. Since the length and the material properties of all pipe segments are equal, the equation (38) represents the element equation of each pipe element. Let s express the equation (38) for two successive nodes i,j in the following form: Figure : Nodes of a periodically multi-supported pipeline resting on elastic supports. f k k k k u i 3 4 i m i k k k3 k 4 i f j k u 3 k3 k33 k 34 j m j k4 k4 k43 k 44 j (40) Expansion of the above equation in global coordinates for any individual element i-j and superposition of the expanded element equations of all pipeline segments, yields: F k k k3 k u k k k3 k F k3 k3 k33 k k34 k k3 k u k4 k4 k43k k44 k k3 k F k3 k3 k33k k34 k k3 k u k4 k4 k43 k k44 k k3 k FN k3 k3 k33 k k34 k k3 k4 un N k4 k4 k43 k k44 k k3 k 4 N F N k3 k3 k33 k34 u N N k4 k4 k43 k 44 N (4)
8 46 Boundary Elements and Other esh Reduction ethods XXXVIII or in an abbreviated form: F G d (4) Nx NxN Nx The above equation can also be written as d G I 0 (43) Nx NxN NxN Nx FNx Equation (43) is a NxN algebraic system with respect to the following 4N unknowns d u u u u u T Nx 3 3 N N N N F F F F F F T Nx 3 3 N N N N (44) In order for the algebraic system provided by equation (43) to be solvable, it should be completed by N more equations. Let s consider the multi-supported pipeline demonstrated in the Figure. The boundary conditions for this example are u 0 0 F ku 0 F3 ku3 3 0 FN kun N 0 un 0 N 0 (45) The above boundary conditions can be written in the following matrix form: d A B 0 Nx (46) Nx NxN NxN FNx In equation (46) the following is valid:
9 Boundary Elements and Other esh Reduction ethods XXXVIII 47 A NxN k k k (47) B NxN (48) Combination of equations (43), (46) yields the following 4Nx4N algebraic system: NxN NxN Nx 4Nx4N (49) ANxN BNxN FNx G I d The condition for non-trivial solution of the equation (49) is: 0 GNxN INxN A B det 0 NxN NxN (50) The above equation can be solved numerically, providing the values of ω versus the values of U. Critical are the values of U yielding the transition of Im{ω} from a positive value (stability) to a negative one (instability).
10 48 Boundary Elements and Other esh Reduction ethods XXXVIII 5 Implementation in a representative example An S-Glass/Epoxy pipeline is resting on one linear spring located at the middle of its length and on its fixed ends (Figure 3). The interior diameter of the pipeline is Figure 3: S-Glass/epoxy filament wound pipeline resting on a linear spring. D=0.0 m and the laminated wall is composed of N=50 layers with thickness 0.5 mm. The material properties in the principal directions of each layer are E=39 GPa, E=8.6 GPa, G=3.8 GPa, ν=0., and the masses per unit length m and have values m=4.94 Kg/m and =7.85 Kg/m. With the aid of the procedure described above the effect on the value of critical flow velocity of the span between supports, fiber orientation and value of springs elastic constant is going to be estimated. Taking into account the material data, the [abd] matrix given in equation () can be initially obtained for ϑ=±30, 45, 60 deg. Then, the parameter Seq can be calculated by the equation (), yielding Seq(30 )=3437 Nm, Seq(45 )=964 Nm, and Seq(60 )=465 Nm. The vectors {d} and {F} contain the following nodal displacements and forces corresponding to the example of Figure 3: d u u u (5) 3 3 F F F F (5) 3 3 The boundary conditions are incorporated in the matrices [A] and [B]: A k (53)
11 Boundary Elements and Other esh Reduction ethods XXXVIII 49 B (54) The matrix [G] can be obtained by the implementation of equation (4) to the considered example. Therefore, its members contain functions of the variables U and ω. Using the matrices [A], [B], [G], the equation (50) can be derived yielding a correlation of U versus ω. For four values of spring elastic constant, namely k=35, 45, 55, 65 KN/m, the equation (50) is solved numerically in order to calculate the eigen-frequencies ω for three values of fiber orientation, ϑ=±30, 45, 60 deg, and for three values of pipe span, L=0, 5, 0 m for each value of ϑ. Starting from a small initial value U=.0 m/s for the flow velocity which is changed incrementally with step.0 m/s, the eigen-frequency ω is calculated using the standard commercial program athematica. The critical flow velocity values corresponding to the transition of Im{ω} from a positive value (stability) to a negative one (instability) are demonstrated in the Figures 4(a c). It should be noted that the numerical results indicated many changes of Im{ω} from positive to negative values. Therefore, the pipeline has more than one critical flow velocity. In the Figures 4(a c) only the first critical flow velocity is adopted. According to the above results the following main conclusions can be drawn. For each value of fiber orientation angle, short pipe segments increases the value of critical flow velocity, improving thus the stability of pipeline. The percentage of U cr shift is higher for the smallest value of ϑ and decreases progressively for the higher values of ϑ. (a) For each value of span length L, higher values of critical fluid velocities appeared for the smaller values of ϑ. This result was expected since small fiber orientation angles improve the bending stiffness of the pipeline, reducing thus the bending deformation. (b) For the larger values of ϑ and L, the critical flow velocity seems to be independent of the elastic modulus of the spring support. Therefore, stiff supports do not affect the stability of flexible pipelines. However, for pipelines with higher bending stiffness (e.g. ϑ=30 deg), stiff elastic springs (k=35 45 KN/m) seems to increase the critical flow velocity. 6 Conclusions. A finite element methodology for critical flow velocity yielding flow-induced instability in FRP filament wound pipelines has been proposed.. The exhibited procedure takes into account the fiber orientation, wall thickness, elastic modulus of spring supports, and flow parameters.
12 50 Boundary Elements and Other esh Reduction ethods XXXVIII Figure 4: Critical flow velocity versus spring elastic constant for L=0 m, L=5 m, and L=0 m. (a) Fiber orientation angle ϑ=±60 deg, (b) Fiber orientation angle ϑ=±45 deg, (c) Fiber orientation angle ϑ=±30 deg 3. The condition for non-trivial solution of the derived pipeline equation yields the eigen-frequencies ω versus flow velocity U. The value of U causing transition of Im{ω} from a positive to a negative value, defines the boundary between pipes stability and instability. 4. An implementation of the methodology to a simple pipeline resting on an elastic spring, indicated that flow-induced instability can be avoided by short pipeline spans and small fiber orientation angles. Very flexible pipelines cannot be affected by very stiff supports, while stiff pipelines can improve their stability for certain values of stiff supports.
13 Boundary Elements and Other esh Reduction ethods XXXVIII 5 References [] Tiratsoo, J., Editorial, Journal of Pipeline engineering, (6), 03. [] Kollár, L.P., Springer, G.S., echanics of composite structures, Cambridge University Press: Cambridge, 009. [3] Pavlou, D.G., Composite materials in piping applications, Destech Publications: Lancaster, 03. [4] Zaitsev, V.F., Polyanin, A.D., Handbook of exact solutions for ordinary differential equations, CRC Press; Florida 003.
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 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 informationME FINITE ELEMENT ANALYSIS FORMULAS
ME 2353 - FINITE ELEMENT ANALYSIS FORMULAS UNIT I FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PROBLEMS 01. Global Equation for Force Vector, {F} = [K] {u} {F} = Global Force Vector [K] = Global Stiffness
More informationExternal Pressure... Thermal Expansion in un-restrained pipeline... The critical (buckling) pressure is calculated as follows:
External Pressure... The critical (buckling) pressure is calculated as follows: P C = E. t s ³ / 4 (1 - ν ha.ν ah ) R E ³ P C = Critical buckling pressure, kn/m² E = Hoop modulus in flexure, kn/m² t s
More informationAssumptions: beam is initially straight, is elastically deformed by the loads, such that the slope and deflection of the elastic curve are
*12.4 SLOPE & DISPLACEMENT BY THE MOMENT-AREA METHOD Assumptions: beam is initially straight, is elastically deformed by the loads, such that the slope and deflection of the elastic curve are very small,
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 informationComposites Design and Analysis. Stress Strain Relationship
Composites Design and Analysis Stress Strain Relationship Composite design and analysis Laminate Theory Manufacturing Methods Materials Composite Materials Design / Analysis Engineer Design Guidelines
More informationPost Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method
9210-220 Post Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method You should have the following for this examination one answer book scientific calculator No
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 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 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 informationComposite Plate Documentation. John Robinson
Composite Plate Documentation Release v ersion John Robinson February 18, 2016 Contents 1 Indices and tables 3 2 How to Use Composite Plate 5 3 Modules 7 Python Module Index 15 i ii Contents: Contents
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 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 informationDETERMINATION OF EI FOR PULTRUDED GFRP SHEET PILE PANELS. Abstract
DETERMINATION OF EI FOR PULTRUDED GFRP SHEET PILE PANELS Yixin Shao, Cynthia Giroux and Zeid Bdeir McGill University Montreal, Quebec, Canada Abstract The flexural rigidity, EI, plays an especially important
More informationVIBRATION CONTROL OF RECTANGULAR CROSS-PLY FRP PLATES USING PZT MATERIALS
Journal of Engineering Science and Technology Vol. 12, No. 12 (217) 3398-3411 School of Engineering, Taylor s University VIBRATION CONTROL OF RECTANGULAR CROSS-PLY FRP PLATES USING PZT MATERIALS DILEEP
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 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 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 informationLinear elastic analysis of thin laminated beams with uniform and symmetric cross-section
Applied and Computational Mechanics 2 (2008) 397 408 Linear elastic analysis of thin laminated beams with uniform and symmetric cross-section M. Zajíček a, a Faculty of Applied Sciences, UWB in Pilsen,
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 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 informationThe Effect of Distribution for a Moving Force
Paper Number 66, Proceedings of ACOUSTICS 2011 2-4 November 2011, Gold Coast, Australia The Effect of Distribution for a Moving Force Ahmed M. Reda (1,2), Gareth L. Forbes (2) (1) Atkins, Perth, Australia
More informationMATERIAL PROPERTIES. Material Properties Must Be Evaluated By Laboratory or Field Tests 1.1 INTRODUCTION 1.2 ANISOTROPIC MATERIALS
. MARIAL PROPRIS Material Properties Must Be valuated By Laboratory or Field ests. INRODUCION he fundamental equations of structural mechanics can be placed in three categories[]. First, the stress-strain
More informationSoftware Verification
EXAMPLE 16 racked Slab Analysis RAKED ANALYSIS METHOD The moment curvature diagram shown in Figure 16-1 depicts a plot of the uncracked and cracked conditions, Ψ 1 State 1, and, Ψ State, for a reinforced
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 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 informationSoftware Verification
PROGRAM NAME: SAFE 014 EXAMPLE 16 racked Slab Analysis RAKED ANALYSIS METHOD The moment curvature diagram shown in Figure 16-1 depicts a plot of the uncracked and cracked conditions, 1 State 1, and, State,
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 informationEFFECT OF LAMINATION ANGLE AND THICKNESS ON ANALYSIS OF COMPOSITE PLATE UNDER THERMO MECHANICAL LOADING
Journal of MECHANICAL ENGINEERING Strojnícky časopis, VOL 67 (217), NO 1, 5-22 EFFECT OF LAMINATION ANGLE AND THICKNESS ON ANALYSIS OF COMPOSITE PLATE UNDER THERMO MECHANICAL LOADING Arnab Choudhury 1,
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 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 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 informationA short review of continuum mechanics
A short review of continuum mechanics Professor Anette M. Karlsson, Department of Mechanical ngineering, UD September, 006 This is a short and arbitrary review of continuum mechanics. Most of this material
More informationPart D: Frames and Plates
Part D: Frames and Plates Plane Frames and Thin Plates A Beam with General Boundary Conditions The Stiffness Method Thin Plates Initial Imperfections The Ritz and Finite Element Approaches A Beam with
More informationModule 2. Analysis of Statically Indeterminate Structures by the Matrix Force Method
Module 2 Analysis of Statically Indeterminate Structures by the Matrix Force Method Lesson 8 The Force Method of Analysis: Beams Instructional Objectives After reading this chapter the student will be
More informationComposite angle ply laminates and netting analysis
10.1098/rspa.2002.1066 FirstCite e-publishing Composite angle ply laminates and netting analysis By J. T. Evans and A. G. Gibson School of Mechanical and Systems Engineering, University of Newcastle upon
More informationChapter 5 Structural Elements: The truss & beam elements
Institute of Structural Engineering Page 1 Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 2 Chapter Goals Learn how to formulate the Finite Element Equations
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
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 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 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 informationPES Institute of Technology
PES Institute of Technology Bangalore south campus, Bangalore-5460100 Department of Mechanical Engineering Faculty name : Madhu M Date: 29/06/2012 SEM : 3 rd A SEC Subject : MECHANICS OF MATERIALS Subject
More informationCOMPOSITE PLATE THEORIES
CHAPTER2 COMPOSITE PLATE THEORIES 2.1 GENERAL Analysis of composite plates is usually done based on one of the following the ries. 1. Equivalent single-layer theories a. Classical laminate theory b. Shear
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 informationContinuum Mechanics and the Finite Element Method
Continuum Mechanics and the Finite Element Method 1 Assignment 2 Due on March 2 nd @ midnight 2 Suppose you want to simulate this The familiar mass-spring system l 0 l y i X y i x Spring length before/after
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 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 informationFinite Difference Dynamic Analysis of Railway Bridges Supported by Pasternak Foundation under Uniform Partially Distributed Moving Railway Vehicle
, October 21-23, 2015, San Francisco, USA Finite Difference Dynamic Analysis of Railway Bridges Supported by Pasternak Foundation under Uniform Partially Distributed Moving Railway Vehicle M. C. Agarana
More informationChapter 4-b Axially Loaded Members
CIVL 222 STRENGTH OF MATERIALS Chapter 4-b Axially Loaded Members AXIAL LOADED MEMBERS Today s Objectives: Students will be able to: a) Determine the elastic deformation of axially loaded member b) Apply
More informationComparison of Ply-wise Stress-Strain results for graphite/epoxy laminated plate subjected to in-plane normal loads using CLT and ANSYS ACP PrepPost
Comparison of Ply-wise Stress-Strain results for graphite/epoxy laminated plate subjected to in-plane normal loads using CLT and ANSYS ACP PrepPost 1 Mihir A. Mehta, 2 Satyen D. Ramani 1 PG Student, Department
More informationConstitutive Equations
Constitutive quations David Roylance Department of Materials Science and ngineering Massachusetts Institute of Technology Cambridge, MA 0239 October 4, 2000 Introduction The modules on kinematics (Module
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 informationMechanics of Materials II. Chapter III. A review of the fundamental formulation of stress, strain, and deflection
Mechanics of Materials II Chapter III A review of the fundamental formulation of stress, strain, and deflection Outline Introduction Assumtions and limitations Axial loading Torsion of circular shafts
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 informationTopic 5: Finite Element Method
Topic 5: Finite Element Method 1 Finite Element Method (1) Main problem of classical variational methods (Ritz method etc.) difficult (op impossible) definition of approximation function ϕ for non-trivial
More informationUsing the finite element method of structural analysis, determine displacements at nodes 1 and 2.
Question 1 A pin-jointed plane frame, shown in Figure Q1, is fixed to rigid supports at nodes and 4 to prevent their nodal displacements. The frame is loaded at nodes 1 and by a horizontal and a vertical
More informationNumerical Stability Analysis of Composite Plates Thermally Stiffened by Finite Element Method
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
More informationChapter 11. Displacement Method of Analysis Slope Deflection Method
Chapter 11 Displacement ethod of Analysis Slope Deflection ethod Displacement ethod of Analysis Two main methods of analyzing indeterminate structure Force method The method of consistent deformations
More informationWORKBOOK MECHANICS OF MATERIALS AND ELEMENTS OF ENGINEERING STRUCTURES
WORKBOOK MECHANICS OF MATERIALS AND ELEMENTS OF ENGINEERING STRUCTURES LUBLIN 014 Authors: Sylwester Samborski, Andrzej Teter and Marcin Bocheński Desktop publishing: Sylwester Samborski, Andrzej Teter
More informationIDE 110 Mechanics of Materials Spring 2006 Final Examination FOR GRADING ONLY
Spring 2006 Final Examination STUDENT S NAME (please print) STUDENT S SIGNATURE STUDENT NUMBER IDE 110 CLASS SECTION INSTRUCTOR S NAME Do not turn this page until instructed to start. Write your name on
More informationStructural Dynamics. Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma).
Structural Dynamics Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma). We will now look at free vibrations. Considering the free
More informationME 7502 Lecture 2 Effective Properties of Particulate and Unidirectional Composites
ME 75 Lecture Effective Properties of Particulate and Unidirectional Composites Concepts from Elasticit Theor Statistical Homogeneit, Representative Volume Element, Composite Material Effective Stress-
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationSEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by
SEMM Mechanics PhD Preliminary Exam Spring 2014 1. Consider a two-dimensional rigid motion, whose displacement field is given by u(x) = [cos(β)x 1 + sin(β)x 2 X 1 ]e 1 + [ sin(β)x 1 + cos(β)x 2 X 2 ]e
More informationDYNAMIC RESPONSE OF SYNTACTIC FOAM CORE SANDWICH USING A MULTIPLE SCALES BASED ASYMPTOTIC METHOD
ECCM6-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Seville, Spain, -6 June 4 DYNAMIC RESPONSE OF SYNTACTIC FOAM CORE SANDWICH USING A MULTIPLE SCALES BASED ASYMPTOTIC METHOD K. V. Nagendra Gopal a*,
More informationInfluence of the filament winding process variables on the mechanical behavior of a composite pressure vessel
Influence of the filament winding process variables on the mechanical behavior of a composite pressure vessel G. Vargas 1 & A. Miravete 2 1 Universidad Pontificia Bolivariana, Facultad de Ingeniería Mecánica,
More informationBEAM DEFLECTION THE ELASTIC CURVE
BEAM DEFLECTION Samantha Ramirez THE ELASTIC CURVE The deflection diagram of the longitudinal axis that passes through the centroid of each cross-sectional area of a beam. Supports that apply a moment
More informationA new closed-form model for isotropic elastic sphere including new solutions for the free vibrations problem
A new closed-form model for isotropic elastic sphere including new solutions for the free vibrations problem E Hanukah Faculty of Mechanical Engineering, Technion Israel Institute of Technology, Haifa
More information16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations
6.2 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive quations Constitutive quations For elastic materials: If the relation is linear: Û σ ij = σ ij (ɛ) = ρ () ɛ ij σ ij =
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 informationChapter 4 Deflection and Stiffness
Chapter 4 Deflection and Stiffness Asst. Prof. Dr. Supakit Rooppakhun Chapter Outline Deflection and Stiffness 4-1 Spring Rates 4-2 Tension, Compression, and Torsion 4-3 Deflection Due to Bending 4-4 Beam
More informationLoss of Accuracy Using Smeared Properties in Composite Beam Modeling
Purdue University Purdue e-pubs Open Access Theses Theses and Dissertations Spring 2015 Loss of Accuracy Using Smeared Properties in Composite Beam Modeling Ning Liu Purdue University Follow this and additional
More informationDESIGN OF LAMINATES FOR IN-PLANE LOADING
DESIGN OF LAMINATES FOR IN-PLANOADING G. VERCHERY ISMANS 44 avenue F.A. Bartholdi, 72000 Le Mans, France Georges.Verchery@m4x.org SUMMARY This work relates to the design of laminated structures primarily
More informationPreliminaries: Beam Deflections Virtual Work
Preliminaries: Beam eflections Virtual Work There are several methods available to calculate deformations (displacements and rotations) in beams. They include: Formulating moment equations and then integrating
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 information6. Bending CHAPTER OBJECTIVES
CHAPTER OBJECTIVES Determine stress in members caused by bending Discuss how to establish shear and moment diagrams for a beam or shaft Determine largest shear and moment in a member, and specify where
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 informationMechanics 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 What programs are in PROMAL? Master Menu The master menu screen with five separate applications from
More informationBasic Energy Principles in Stiffness Analysis
Basic Energy Principles in Stiffness Analysis Stress-Strain Relations The application of any theory requires knowledge of the physical properties of the material(s) comprising the structure. We are limiting
More informationFIXED BEAMS IN BENDING
FIXED BEAMS IN BENDING INTRODUCTION Fixed or built-in beams are commonly used in building construction because they possess high rigidity in comparison to simply supported beams. When a simply supported
More informationAnalysis of the geometrical dependence of auxetic behavior in reentrant structures by finite elements
Acta Mech. Sin. (2016) 32(2):295 300 DOI 10.1007/s10409-015-0534-2 RESEARCH PAPER Analysis of the geometrical dependence of auxetic behavior in reentrant structures by finite elements V. H. Carneiro 1
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 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 informationANALYSIS OF COMPOSITE RECTANGULAR PLATES BASED ON THE CLASSICAL LAMINATED PLATE THEORY
The 4th International Conference Advanced Composite Materials Engineering COMAT 2012 18-20 October 2012, Brasov, Romania ANALYSIS OF COMPOSITE RECTANGULAR PLATES BASED ON THE CLASSICAL LAMINATED PLATE
More informationMacroscopic theory Rock as 'elastic continuum'
Elasticity and Seismic Waves Macroscopic theory Rock as 'elastic continuum' Elastic body is deformed in response to stress Two types of deformation: Change in volume and shape Equations of motion Wave
More informationApplication of Laplace Iteration method to Study of Nonlinear Vibration of laminated composite plates
(3) 78 795 Application of Laplace Iteration method to Study of Nonlinear Vibration of laminated composite plates Abstract In this paper, free vibration characteristics of laminated composite plates are
More informationParameter Design of High Speed Coupling for 6 MW Wind Turbine Considering Torsional Vibration
Parameter Design of High Speed Coupling for 6 MW Wind Turbine Considering Torsional Vibration JongHun Kang 1, Junwoo Bae 2, Seungkeun Jeong 3, SooKeun Park 4 and Hyoung Woo Lee 1 # 1 Department of Mechatronics
More informationneeded to buckle an ideal column. Analyze the buckling with bending of a column. Discuss methods used to design concentric and eccentric columns.
CHAPTER OBJECTIVES Discuss the behavior of columns. Discuss the buckling of columns. Determine the axial load needed to buckle an ideal column. Analyze the buckling with bending of a column. Discuss methods
More informationMECHANICS OF COMPOSITE STRUCTURES
MECHANICS OF COMPOSITE STRUCTURES LÁSZLÓ P. KOLLÁR Budapest University of Technology and Economics GEORGE S. SPRINGER Stanford University PUBLISHED BY THE PRESS SYNDICATE OF THE UNIVERSITY OF CAMBRIDGE
More informationEDEM DISCRETIZATION (Phase II) Normal Direction Structure Idealization Tangential Direction Pore spring Contact spring SPRING TYPES Inner edge Inner d
Institute of Industrial Science, University of Tokyo Bulletin of ERS, No. 48 (5) A TWO-PHASE SIMPLIFIED COLLAPSE ANALYSIS OF RC BUILDINGS PHASE : SPRING NETWORK PHASE Shanthanu RAJASEKHARAN, Muneyoshi
More informationBI-STABLE COMPOSITE SHELLS
AIAA 2-1385 BI-STABLE COMPOSITE SHELLS K. Iqbal and S. Pellegrino Department of Engineering, University of Cambridge, Trumpington Street, Cambridge, CB2 1PZ, U.K. Abstract This paper is concerned with
More informationCHAPTER 4: BENDING OF BEAMS
(74) CHAPTER 4: BENDING OF BEAMS This chapter will be devoted to the analysis of prismatic members subjected to equal and opposite couples M and M' acting in the same longitudinal plane. Such members are
More informationAn integrated approach to the design of high performance carbon fibre reinforced risers - from micro to macro - scale
An integrated approach to the design of high performance carbon fibre reinforced risers - from micro to macro - scale Angelos Mintzas 1, Steve Hatton 1, Sarinova Simandjuntak 2, Andrew Little 2, Zhongyi
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 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 informationNumerical simulation of the coil spring and investigation the impact of tension and compression to the spring natural frequencies
Numerical simulation of the coil spring and investigation the impact of tension and compression to the spring natural frequencies F. D. Sorokin 1, Zhou Su 2 Bauman Moscow State Technical University, Moscow,
More informationISHIK UNIVERSITY DEPARTMENT OF MECHATRONICS ENGINEERING
ISHIK UNIVERSITY DEPARTMENT OF MECHATRONICS ENGINEERING QUESTION BANK FOR THE MECHANICS OF MATERIALS-I 1. A rod 150 cm long and of diameter 2.0 cm is subjected to an axial pull of 20 kn. If the modulus
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 informationBENDING VIBRATIONS OF ROTATING NON-UNIFORM COMPOSITE TIMOSHENKO BEAMS WITH AN ELASTICALLY RESTRAINED ROOT
BENDING VIBRATIONS OF ROTATING NON-UNIFORM COMPOSITE TIMOSHENKO BEAMS WITH AN EASTICAY RESTRAINED ROOT Sen Yung ee and Jin Tang Yang Department of Mechanical Engineering, National Cheng Kung University,
More information