International Journal of Advanced Engineering Technology E-ISSN
|
|
- Marianna Peters
- 6 years ago
- Views:
Transcription
1 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, Faculty of Tech. & Engg., M. S. University of Baroda, Vadodara , Gujarat, India Professor,Department of Civil Engineering, Parul Institute of Engg. and Tech., Limda, Vadodara , Gujarat, India ABSTRACT Solution to fiber reinforced composite plate problem under different types of loading and boundary conditions is attempted using Integrated Force Method (IFM). Internal moments developed due to external loading are calculated by imposing simultaneously the equilibrium equations and compatibility conditions. Various closed form matrices are developed by discretizing the potential strain energy for a four noded rectangular element having 9 force and 1 displacement degrees of freedom; which is named as RECT_9F_1D element. Using symmetry, the influence of aspect ratio on internal moments and deflections is studied. Matlab is exploited for plotting the moment contours and deflected shape of transversely loaded plates. The numerical results obtained for twelve Glass Fiber Reinforced Plastic (GFRP) plate problems are compared with the available small deflection theory results to validate the proposed formulation and its computer implementation; a good agreement is found. KEY WORDS: IFM, Orthotropic plate, Matlab, Moment contours. INTRODUCTION In many practical applications of orthotropic plates, rather than following hypothetical assumptions, it is important to consider actual directional dependent bending characteristics, behavior and stiffness. Generally speaking anisotropy can be introduced in plate either by modifying material constituents or by providing ribs, stiffeners or corrugations for which an exact analytical solution in terms of deflection and internal moments exists only for some simple boundary conditions [1]. The complete solution of plates with clamped boundary conditions is possible only through numerical or approximate method such as Rayleigh Ritz method []. Other methods such as Differential Quadrature Method, Modified Spline Function Method, Superposition Technique and Modified Fourier Series Method consist of very complex and tedious mathematical formulation [3]. Thus, there is a need of a method, which can be easily applied to any type of boundary condition and can easily handle any type of loading condition. Direct Moment Calculation (DMC) approach, which is the key aspect of Integrated Force Method (IFM) is well established since last few decades [4, 5]. The present work includes basic equations of IFM for a four noded rectangular plate bending element. Different closed-form matrices such as Global Equilibrium Matrix [S], Global Flexibility Matrix [G] and Auto generated Matlab based Global Compatibility Matrix [C] are derived by discretizing the expressions of potential and complimentary strain energies. The familiar process of differentiation used in popular stiffness method to evaluate stresses from displacements is avoided in the integrated force method. These matrices are assembled using the traditional procedure. Total twelve problems of orthotropic plate bending with variation in geometrical dimensions, boundary conditions and loading are considered. Using symmetry, quarter plate is discretized in 5 x 5 mesh. Nodal displacements and internal moments are calculated. Using Matlab software [6] moment contours and deflected shape for each case are plotted with maxima/minima values. Influence of aspect ratio for fixed and simply supported edge conditions under central point load and uniform lateral pressure are studied. Results obtained by the proposed method are compared with those based on classical small deflection theory [7]. FUNDAMENTAL THEORY OF IFM The IFM equations for a continuum discretized into finite number of elements with n and m force and displacement degrees of freedom respectively, are obtained by coupling the m` number of equilibrium equation and r =n m compatibility conditions. The m equilibrium equations (EE) are written as [4, 5] [B] {F} = {P} (1) and the r compatibility conditions are written as [C] [G] {F} = {δr} () These conditions are combined to obtain the IFM governing equations as [ B] [ ][ ] { } { P} { } F = Or [S] {F} = {P}. (3) C G δr From forces {F}, displacements {X} are calculated using the following equation. {X} = [J] {[G] {F} + {β 0 } (4) where, [J] = m rows of [[S] -1 ] T, [B] is of m x n size matrix, which is sparse and unsymmetrical, [G] is the symmetrical flexibility matrix; it is a block-diagonal matrix where each block represents a flexibility matrix for an element, [C] is the compatibility matrix of size r x n, {δr} = -[C] {β} 0 is the effective deformation vector with {β} 0 being the initial deformation vector of dimension n, [S] is the IFM
2 governing unsymmetrical matrix of size n x n and [J] is the m x n size deformation coefficient matrix which is calculated from [S] matrix. IFM has following 3 matrices that are equivalent to the operators of classical elasticity theory: 1. The equilibrium matrix [B] which links internal forces to external loads,. The compatibility matrix [C] which governs the deformations and 3. The flexibility matrix [G] which relates deformations to forces. Both the equilibrium and compatibility matrices of the IFM are unsymmetrical, whereas the material constitutive matrix and the flexibility matrix are symmetrical. FORMULATION FOR A RECTANGULAR ELEMENT Element Equilibrium Matrix [Be] represents the curvatures corresponding to each internal moment. Consider a 4-noded, 1 displacement degrees of freedom (ddof) (X 1 to X 1 ) rectangular element of thickness t with dimensions as a x b along the x and y axes as shown in Figs. 1 and. The force field is chosen in terms of nine independent forces as (7) Relations between internal moments and independent forces are written as Arranging Eqs. (8) to (10) in matrix form,... (8)... (9)... (10) Fig. 1 Rectangular Element or { } [ Y]{ } (11) M = F e... (1) where {F e } is the element force vector. The variation of above moments is considered bilinear along both the directions. The displacement field satisfies the continuity condition and the selected forces also satisfy the mandatory requirement. Polynomial function for lateral displacement can be written as U Fig. Nodal Displacements The elemental equilibrium matrix written in terms of forces at grid points represents the vectorial summation of n internal forces {F} and m external loads {P}. The nodal equilibrium equation in matrix notation can be stored as rectangular matrix [B e ] of size m x n. The variational functional is evaluated as a portion of IFM functional which yields the basic elemental equilibrium matrix [B e ] in explicit form as follows: p = D M x = { } D w + M x T y w + M y xy w dxdy x y (5) M { ε}ds (6) Where, {M} T =(M x, M y, M xy ) are the in-plane internal moments and {Є} T = (13) Or w(x,y) = [A]{α}.. (14) where [A] is a function of x and y matrix and {α} is a vector of constants. Substituting coordinates, one can find constants and finally the interpolation matrix [N] following the usual finite element procedure. Here each component of [N] is associated with nodal displacements X 1, X, X 1 as shown in Fig.. Eq. 14 is now expressed in terms of nodal displacements as w(x,y) = [N]{X}... (15) By arranging all force and displacement functions properly, one can discretize the Eq. 5 to obtain the elemental equilibrium matrix as U e = {X} T [B e ]{F}... (16) where [B e ] = s [Z] T [Y] ds... (17) Here [Z] = [L][N] where [L] is the differential operator matrix, [N] is the displacement interpolation function matrix and [Y] is the force interpolation function matrix. The integration yields the following non symmetrical equilibrium matrix [B e ] for the element.
3 [B e ]= G e 77 =, G e 88 =, G e 99 = (18) Elemental Flexibility Matrix (Orthotropic) [G e ortho]: The element flexibility matrix is obtained by discretizing the complementary strain energy as [G e ortho]= s [Y] T [D ortho ][Y] dx dy (19) where [D ortho ] is the orthotropic plate material matrix of size ( 3 x 3 ) which is written as (0) where, D x = - E x h 3 /1(1- and D y = - E y h 3 /1(1- are the plate rigidities along x and y directions respectively, D xy = G xy h 3 /1 is torsional rigidity in xy plane, E x, E y and G xy are the Young s modulus of elasticity in x and y directions and shear modulus of the orthotropic material in shear plane respectively, are the Poisson s ratios along orthogonal axes and h is the thickness of plate. Substituting matrices given in Eqs. 11 and 0 in Eq. 19 and integrating within domain (a x b) yields the symmetrical flexibility matrix [G e ortho] of size (9 x 9). The nonzero components of the matrix are as follows: G e 11 =, G e 15 = G e 51 = - 11, G e =, G e 6 = G e 6 = - G e, G e 33 =, G e 37 = G e 73 = - G e 33, G e 44 =, G e 48 = G e 84 = - 44, G e 55 =, G e 66 =, Global Compatibility Matrix [C] The compatibility matrix is obtained from the deformation displacement relation ({β} = [B] T {X}). In DDR all the deformations are expressed in terms of all possible nodal displacements and the r compatibility conditions are developed in terms of internal forces i.e., F 1, F n, where n is the total number of internal forces in a given problem. The global compatibility matrix [C] can be evaluated by multiplying the coefficients of compatibility conditions developed by using Matlab and the global flexibility matrix [G]. NUMERICAL EXAMPLES Results of twelve different examples of laterally loaded glass fiber reinforced plastic plate are presented to validate the suggested formulation and computer implementation of integrated force method. For numerical solution, a plate of size m x m of thickness t = 3mm and having Young s modulus in x and y directions as 40kN/mm and 8kN/mm respectively and Poisson s ratio υ xy as 0.5 is considered. The shear modulus in xy plane is considered as 4kN/mm. Due to double symmetry, only quarter of the plate is discretized into 5 x 5 mesh. The following variations are considered in solving the rectangular plate problems. 1. Element aspect ratio R = 1.0, 1.5 and.0.. Central Point Load (CPL) of 10kN and Uniform Lateral Pressure (ULP) of 10kN/m. 3. Support conditions: Fully clamped and Simply supported. Following steps are used for solving the above problems: Step 1: A four-noded rectangular element (a x b) with 1 ddof and 9 fdof is used to discretize the problem using two- way symmetry. The elemental [B e ] matrix is obtained by substituting a = 0.1m, b = 0.1m for orthotropic case in Eq. 18 considering 5 x 5 discretization scheme. After assembling, the size of global equilibrium matrix will be 75 x 5. Step : The compatibility matrix for the twenty five elements is obtained from the displacement deformation relations (DDR) i.e. β = [B] T {X}. In the DDR, 5 deformations which correspond to 5 force variables are expressed in terms of 75 displacements (X 1, X, X 75 ). The problem requires 150 compatibility conditions [C] that are obtained by eliminating the 75 displacements from the 5 DDR using auto-generated Matlab based computer program by giving input as upper part of the global equilibrium matrix. Step 3: The global flexibility matrix [G] for the problem is obtained by diagonal concatenation of the
4 twenty five elemental flexibility matrices calculated as per Eq. 19. Step 4: By multiplying compatibility matrix [C] and global flexibility matrix [G], bottom most part of the global equilibrium matrix is obtained. Assembling both the matrices gives complete [S] matrix of size 5 x 5, which comprises of EE and CC. The forces are obtained by using MatLab s Pseudo inverse command. Step 5: The displacements are calculated by using the relation ({X} = [J][G]{F}), where [J] = m rows of matrix [[S] -1 ] T. Linear static analysis is carried out for all the plate problems and nodal displacements and internal moments are calculated at each node of 5 x 5 discretization pattern. Tables 1 and show deflection and moments Mx, My and Mxy at the centre of plate for various aspect ratios. The values are compared with those given by classical plate theory [1,7]. Two dimensional moment contours as well as three dimensional deflected shapes are drawn for plate using Matlab software. Figures 3 and 4 show deformation patterns for simply supported plate whereas Figs. 5 and 6 show the same for clamped plate under central point loading and uniform lateral pressure. Figures 7 thru 9 show moment contours for Mx, My and Mxy respectively for simply supported plate subjected to central point loading. LOAD CPL ULP LOAD CPL ULP Table 1 Results obtained for Simply Supported GFRP Plate R w c (m) Mxc ( N-m) Myc ( N-m) Mxyc( N-m) IFM Exact 1 IFM Exact 1 IFM Exact 1 IFM Exact Table Results obtained for Fully Clamped GFRP Plate R w c (m) Mxc ( N-m) Myc ( N-m) Mxyc( N-m) IFM Exact 7 IFM Exact 7 IFM Exact 7 IFM Exact Fig. 3 Deformed Shape of a SS Plate under CPL
5 Fig. 4 Deformed Shape of a SS Plate under ULP Fig. 5 Deformed Shape of a Clamped Plate under CPL Fig. 6 Deformed Shape of a Clamped Plate under ULP
6 Fig. 7 Contours of Mxx for SS Orthotropic Plate under CPL along yy Axis along yy Axis Fig. 8 Contours of Myy for SS Orthotropic Plate under CPL
7 along yy Axis Fig. 9 Contours of Mxy for SS Orthotropic Plate under CPL CONCLUSIONS A method of analysis, which imposes simultaneously equilibrium equations and compatibility conditions to calculate internal moments, is presented for the analysis of transversely loaded fiber reinforced composite plates. Results obtained for a number of examples with change in aspect ratio, boundary condition and loading condition show good agreement with theoretical results. However, with the increase in aspect ratio, the difference in results is found to increase; which may be improved by finer discretization in the longer direction. The use of Matlab software in conjunction with the software developed based on integrated force method in Visual Basic makes the whole exercise quite attractive as it provides directly graphical output in the form of deflected shape and moment contours. The method has good potential for its application to reinforced concrete slab, stiffened plates and corrugated sheets. REFERENCES 1. Chandrashekhar K., Theory of Plates, Universities Press (India) Ltd., Hyderabad, Whitney J. M., Structural Analysis of Laminated Plates, Lancaster, PA, Technomic Publishing Co., Gorman D. J., Free Vibration Analysis of Rectangular Plates using a set of Static Beam Function in Rayleigh-Ritz Method, Journal of Sound Vibration, Vol. 189, No.1, pp.81-87, Patnaik S. N., Integrated Force Method for Discrete Analysis, International Journal of Numerical Methods and Engineering, Vol. 41, pp , Doiphode G. S., Kulkarni S. M. and Patodi S. C., Improving Plate Bending Solutions using Integrated Force Method, Proceedings of the Sixth Structural Engineering Convention (SEC 008), IIT Chennai, pp. 7 35, Dec Matlab: Manual on Graphical Tool for D and 3D Plots, The Mathworks, Inc., Natick, MA , USA. 7. Bairagi N. K., Plate Analysis, Khanna Publishers, New Delhi, 1986.
CHAPTER 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 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 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 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 informationBending 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 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 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 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 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 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 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 informationJEPPIAAR ENGINEERING COLLEGE
JEPPIAAR ENGINEERING COLLEGE Jeppiaar Nagar, Rajiv Gandhi Salai 600 119 DEPARTMENT OFMECHANICAL ENGINEERING QUESTION BANK VI SEMESTER ME6603 FINITE ELEMENT ANALYSIS Regulation 013 SUBJECT YEAR /SEM: III
More informationParametric study on the transverse and longitudinal moments of trough type folded plate roofs using ANSYS
American Journal of Engineering Research (AJER) e-issn : 2320-0847 p-issn : 2320-0936 Volume-4 pp-22-28 www.ajer.org Research Paper Open Access Parametric study on the transverse and longitudinal moments
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 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 informationFINITE GRID SOLUTION FOR NON-RECTANGULAR PLATES
th International Conference on Earthquake Geotechnical Engineering June 5-8, 7 Paper No. 11 FINITE GRID SOLUTION FOR NON-RECTANGULAR PLATES A.Halim KARAŞĐN 1, Polat GÜLKAN ABSTRACT Plates on elastic foundations
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 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 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 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 informationCHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES
CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES * Governing equations in beam and plate bending ** Solution by superposition 1.1 From Beam Bending to Plate Bending 1.2 Governing Equations For Symmetric
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 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 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 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 informationELASTOPLASTIC STEEL BEAM BENDING ANALYSIS BY USING ABAQUS
11 ELASTOPLASTIC STEEL BEAM BENDING ANALYSIS BY USING ABAQUS Biswajit Jena * * Corresponding author: biswa.tech88@gmail.com Abstract: In recent years tremendous efforts have been made in development of
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 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 informationLevel 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method
9210-203 Level 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method You should have the following for this examination one answer book No additional data is attached
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 informationBilinear Quadrilateral (Q4): CQUAD4 in GENESIS
Bilinear Quadrilateral (Q4): CQUAD4 in GENESIS The Q4 element has four nodes and eight nodal dof. The shape can be any quadrilateral; we ll concentrate on a rectangle now. The displacement field in terms
More informationEffect of Specimen Dimensions on Flexural Modulus in a 3-Point Bending Test
Effect of Specimen Dimensions on Flexural Modulus in a 3-Point Bending Test M. Praveen Kumar 1 and V. Balakrishna Murthy 2* 1 Mechanical Engineering Department, P.V.P. Siddhartha Institute of Technology,
More informationCodal Provisions IS 1893 (Part 1) 2002
Abstract Codal Provisions IS 1893 (Part 1) 00 Paresh V. Patel Assistant Professor, Civil Engineering Department, Nirma Institute of Technology, Ahmedabad 38481 In this article codal provisions of IS 1893
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 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 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 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 informationOptimum Height of Plate Stiffener under Pressure Effect
The st Regional Conference of Eng. Sci. NUCEJ Spatial ISSUE vol., No.3, 8 pp 459-468 Optimum Height of Plate Stiffener under Pressure Effect Mazin Victor Yousif M.Sc Production Engineering University of
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 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 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 informationACCURATE FREE VIBRATION ANALYSIS OF POINT SUPPORTED MINDLIN PLATES BY THE SUPERPOSITION METHOD
Journal of Sound and Vibration (1999) 219(2), 265 277 Article No. jsvi.1998.1874, available online at http://www.idealibrary.com.on ACCURATE FREE VIBRATION ANALYSIS OF POINT SUPPORTED MINDLIN PLATES BY
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 informationAnalytical Strip Method for Thin Isotropic Cylindrical Shells
IOSR Journal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 2278-1684,p-ISSN: 2320-334X, Volume 14, Issue 4 Ver. III (Jul. Aug. 2017), PP 24-38 www.iosrjournals.org Analytical Strip Method for
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 informationVerification Examples. FEM-Design. version
FEM-Design 6.0 FEM-Design version. 06 FEM-Design 6.0 StruSoft AB Visit the StruSoft website for company and FEM-Design information at www.strusoft.com Copyright 06 by StruSoft, all rights reserved. Trademarks
More informationLecture 7: The Beam Element Equations.
4.1 Beam Stiffness. A Beam: A long slender structural component generally subjected to transverse loading that produces significant bending effects as opposed to twisting or axial effects. MECH 40: Finite
More informationComposite Structural Mechanics using MATLAB
Session 2520 Composite Structural Mechanics using MATLAB Oscar Barton, Jr., Jacob B. Wallace United States Naval Academy Annapolis, Md 21402 Abstract In this paper MATLAB is adopted as the programming
More informationGATE SOLUTIONS E N G I N E E R I N G
GATE SOLUTIONS C I V I L E N G I N E E R I N G From (1987-018) Office : F-16, (Lower Basement), Katwaria Sarai, New Delhi-110016 Phone : 011-65064 Mobile : 81309090, 9711853908 E-mail: info@iesmasterpublications.com,
More informationMIXED RECTANGULAR FINITE ELEMENTS FOR PLATE BENDING
144 MIXED RECTANGULAR FINITE ELEMENTS FOR PLATE BENDING J. N. Reddy* and Chen-Shyh-Tsay School of Aerospace, Mechanical and Nuclear Engineering, University of Oklahoma, Norman, Oklahoma The paper describes
More informationAnalysis of Rectangular Plate with Opening by Finite Difference Method
American Journal of Civil Engineering and Architecture, 2015, Vol. 3, No. 5, 165-173 Available online at http://pubs.sciepub.com/ajcea/3/5/3 Science and Education Publishing DOI:10.12691/ajcea-3-5-3 Analysis
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 informationAN INTEGRATED KIRCHHOFF PLATE ELEMENT BY GALERKIN METHOD FOR THE ANALYSIS OF PLATES ON ELASTIC FOUNDATION
AN INTEGRATED IRCHHOFF PLATE ELEMENT BY GALERIN METHOD FOR THE ANALYSIS OF PLATES ON ELASTIC FOUNDATION Ragesh.P.P, V.Mustafa 2, T.P.Somasundaran 3 (Research Scholar, National Institute of Technology Calicut,
More informationStatic & Dynamic. Analysis of Structures. Edward L.Wilson. University of California, Berkeley. Fourth Edition. Professor Emeritus of Civil Engineering
Static & Dynamic Analysis of Structures A Physical Approach With Emphasis on Earthquake Engineering Edward LWilson Professor Emeritus of Civil Engineering University of California, Berkeley Fourth Edition
More informationVIBRATION AND DAMPING ANALYSIS OF FIBER REINFORCED COMPOSITE MATERIAL CONICAL SHELLS
VIBRATION AND DAMPING ANALYSIS OF FIBER REINFORCED COMPOSITE MATERIAL CONICAL SHELLS Mechanical Engineering Department, Indian Institute of Technology, New Delhi 110 016, India (Received 22 January 1992,
More informationPILE SOIL INTERACTION MOMENT AREA METHOD
Pile IGC Soil 2009, Interaction Moment Guntur, INDIA Area Method PILE SOIL INTERACTION MOMENT AREA METHOD D.M. Dewaikar Professor, Department of Civil Engineering, IIT Bombay, Mumbai 400 076, India. E-mail:
More informationVibration of Thin Beams by PIM and RPIM methods. *B. Kanber¹, and O. M. Tufik 1
APCOM & ISCM -4 th December, 23, Singapore Vibration of Thin Beams by PIM and RPIM methods *B. Kanber¹, and O. M. Tufik Mechanical Engineering Department, University of Gaziantep, Turkey. *Corresponding
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 information2766. Differential quadrature method (DQM) for studying initial imperfection effects and pre- and post-buckling vibration of plates
2766. Differential quadrature method (DQM) for studying initial imperfection effects and pre- and post-buckling vibration of plates Hesam Makvandi 1, Shapour Moradi 2, Davood Poorveis 3, Kourosh Heidari
More informationELASTICITY AND FRACTURE MECHANICS. Vijay G. Ukadgaonker
THEORY OF ELASTICITY AND FRACTURE MECHANICS y x Vijay G. Ukadgaonker Theory of Elasticity and Fracture Mechanics VIJAY G. UKADGAONKER Former Professor Indian Institute of Technology Bombay Delhi-110092
More informationFIXED BEAMS CONTINUOUS BEAMS
FIXED BEAMS CONTINUOUS BEAMS INTRODUCTION A beam carried over more than two supports is known as a continuous beam. Railway bridges are common examples of continuous beams. But the beams in railway bridges
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 informationMAE 323: Chapter 6. Structural Models
Common element types for structural analyis: oplane stress/strain, Axisymmetric obeam, truss,spring oplate/shell elements o3d solid ospecial: Usually used for contact or other constraints What you need
More informationIntroduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams.
Outline of Continuous Systems. Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams. Vibrations of Flexible Strings. Torsional Vibration of Rods. Bernoulli-Euler Beams.
More informationINSTITUTE OF AERONAUTICAL ENGINEERING (Autonomous) Dundigal, Hyderabad
NSTTUTE OF AERONAUTCAL ENGNEERNG (Autonomous) Dundigal, Hyderabad - 00 043 AERONAUTCAL ENGNEERNG TUTORAL QUESTON BANK Course Name : ARCRAFT VEHCLES STRUCTURES Course Code : A2109 Class : B. Tech Semester
More informationAn Evaluation and Comparison of Models for Maximum Deflection of Stiffened Plates Using Finite Element Analysis
Marine Technology, Vol. 44, No. 4, October 2007, pp. 212 225 An Evaluation and Comparison of Models for Maximum Deflection of Stiffened Plates Using Finite Element Analysis Lior Banai 1 and Omri Pedatzur
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 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 informationME 1401 FINITE ELEMENT ANALYSIS UNIT I PART -A. 2. Why polynomial type of interpolation functions is mostly used in FEM?
SHRI ANGALAMMAN COLLEGE OF ENGINEERING AND TECHNOLOGY (An ISO 9001:2008 Certified Institution) SIRUGANOOR, TIRUCHIRAPPALLI 621 105 Department of Mechanical Engineering ME 1401 FINITE ELEMENT ANALYSIS 1.
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 information202 Index. failure, 26 field equation, 122 force, 1
Index acceleration, 12, 161 admissible function, 155 admissible stress, 32 Airy's stress function, 122, 124 d'alembert's principle, 165, 167, 177 amplitude, 171 analogy, 76 anisotropic material, 20 aperiodic
More informationToward a novel approach for damage identification and health monitoring of bridge structures
Toward a novel approach for damage identification and health monitoring of bridge structures Paolo Martino Calvi 1, Paolo Venini 1 1 Department of Structural Mechanics, University of Pavia, Italy E-mail:
More informationCHENDU COLLEGE OF ENGINEERING &TECHNOLOGY DEPARTMENT OF CIVIL ENGINEERING SUB CODE & SUB NAME : CE2351-STRUCTURAL ANALYSIS-II UNIT-1 FLEXIBILITY
CHENDU COLLEGE OF ENGINEERING &TECHNOLOGY DEPARTMENT OF CIVIL ENGINEERING SUB CODE & SUB NAME : CE2351-STRUCTURAL ANALYSIS-II UNIT-1 FLEXIBILITY METHOD FOR INDETERMINATE FRAMES PART-A(2MARKS) 1. What is
More informationMethods of Analysis. Force or Flexibility Method
INTRODUCTION: The structural analysis is a mathematical process by which the response of a structure to specified loads is determined. This response is measured by determining the internal forces or stresses
More informationTable of Contents. Preface... 13
Table of Contents Preface... 13 Chapter 1. Vibrations of Continuous Elastic Solid Media... 17 1.1. Objective of the chapter... 17 1.2. Equations of motion and boundary conditions of continuous media...
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 informationANALYSIS RADIAL FOLDED PLATE
ANALYSIS RADIAL FOLDED PLATE Ms.Pranoti Satish Bhamare 1, Prajakta S. Bramhankar 2, Pooja G. Baviskar 3 1,2,3 Assistant Professor, Department of Civil Engineering, Guru Gobind Singh College of Engineering
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 informationCivil Engineering Design (1) Analysis and Design of Slabs 2006/7
Civil Engineering Design (1) Analysis and Design of Slabs 006/7 Dr. Colin Caprani, Chartered Engineer 1 Contents 1. Elastic Methods... 3 1.1 Introduction... 3 1. Grillage Analysis... 4 1.3 Finite Element
More informationCOMPUTATIONAL ELASTICITY
COMPUTATIONAL ELASTICITY Theory of Elasticity and Finite and Boundary Element Methods Mohammed Ameen Alpha Science International Ltd. Harrow, U.K. Contents Preface Notation vii xi PART A: THEORETICAL ELASTICITY
More informationStiffness Matrices, Spring and Bar Elements
CHAPTER Stiffness Matrices, Spring and Bar Elements. INTRODUCTION The primary characteristics of a finite element are embodied in the element stiffness matrix. For a structural finite element, the stiffness
More informationAdvanced Structural Analysis EGF Section Properties and Bending
Advanced Structural Analysis EGF316 3. Section Properties and Bending 3.1 Loads in beams When we analyse beams, we need to consider various types of loads acting on them, for example, axial forces, shear
More informationStatic & Free Vibration Analysis of an Isotropic and Orthotropic Thin Plate using Finite Element Analysis (FEA)
International Journal of Current Engineering and Technology E-ISSN 2277 4106, P-ISSN 2347 5161 2015 INPRESSCO, All Rights Reserved Available at http://inpressco.com/category/ijcet Research Article Static
More informationCHAPTER OBJECTIVES Use various methods to determine the deflection and slope at specific pts on beams and shafts: 2. Discontinuity functions
1. Deflections of Beams and Shafts CHAPTER OBJECTIVES Use various methods to determine the deflection and slope at specific pts on beams and shafts: 1. Integration method. Discontinuity functions 3. Method
More informationBIAXIAL STRENGTH INVESTIGATION OF CFRP COMPOSITE LAMINATES BY USING CRUCIFORM SPECIMENS
BIAXIAL STRENGTH INVESTIGATION OF CFRP COMPOSITE LAMINATES BY USING CRUCIFORM SPECIMENS H. Kumazawa and T. Takatoya Airframes and Structures Group, Japan Aerospace Exploration Agency 6-13-1, Ohsawa, Mitaka,
More informationComputational Simulation of Dynamic Response of Vehicle Tatra T815 and the Ground
IOP Conference Series: Earth and Environmental Science PAPER OPEN ACCESS Computational Simulation of Dynamic Response of Vehicle Tatra T815 and the Ground To cite this article: Jozef Vlek and Veronika
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 informationSabah Shawkat Cabinet of Structural Engineering Walls carrying vertical loads should be designed as columns. Basically walls are designed in
Sabah Shawkat Cabinet of Structural Engineering 17 3.6 Shear walls Walls carrying vertical loads should be designed as columns. Basically walls are designed in the same manner as columns, but there are
More informationDiscrete Analysis for Plate Bending Problems by Using Hybrid-type Penalty Method
131 Bulletin of Research Center for Computing and Multimedia Studies, Hosei University, 21 (2008) Published online (http://hdl.handle.net/10114/1532) Discrete Analysis for Plate Bending Problems by Using
More informationExample 3.7 Consider the undeformed configuration of a solid as shown in Figure 3.60.
162 3. The linear 3-D elasticity mathematical model The 3-D elasticity model is of great importance, since it is our highest order hierarchical model assuming linear elastic behavior. Therefore, it provides
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationA METHOD OF LOAD INCREMENTS FOR THE DETERMINATION OF SECOND-ORDER LIMIT LOAD AND COLLAPSE SAFETY OF REINFORCED CONCRETE FRAMED STRUCTURES
A METHOD OF LOAD INCREMENTS FOR THE DETERMINATION OF SECOND-ORDER LIMIT LOAD AND COLLAPSE SAFETY OF REINFORCED CONCRETE FRAMED STRUCTURES Konuralp Girgin (Ph.D. Thesis, Institute of Science and Technology,
More informationIf the number of unknown reaction components are equal to the number of equations, the structure is known as statically determinate.
1 of 6 EQUILIBRIUM OF A RIGID BODY AND ANALYSIS OF ETRUCTURAS II 9.1 reactions in supports and joints of a two-dimensional structure and statically indeterminate reactions: Statically indeterminate structures
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 informationEffect of Mass Matrix Formulation Schemes on Dynamics of Structures
Effect of Mass Matrix Formulation Schemes on Dynamics of Structures Swapan Kumar Nandi Tata Consultancy Services GEDC, 185 LR, Chennai 600086, India Sudeep Bosu Tata Consultancy Services GEDC, 185 LR,
More informationUnit 13 Review of Simple Beam Theory
MIT - 16.0 Fall, 00 Unit 13 Review of Simple Beam Theory Readings: Review Unified Engineering notes on Beam Theory BMP 3.8, 3.9, 3.10 T & G 10-15 Paul A. Lagace, Ph.D. Professor of Aeronautics & Astronautics
More informationTHE BENDING STIFFNESSES OF CORRUGATED BOARD
AMD-Vol. 145/MD-Vol. 36, Mechanics of Cellulosic Materials ASME 1992 THE BENDING STIFFNESSES OF CORRUGATED BOARD S. Luo and J. C. Suhling Department of Mechanical Engineering Auburn University Auburn,
More informationSimplified Nonlinear Analysis of Reinforced Concrete Slabs and Beams
Journal of Civil Engineering Research 2015, 5(6): 162-171 DOI: 10.5923/j.jce.20150506.06 Simplified Nonlinear Analysis of Reinforced Concrete Slabs and Beams Graham Dean Roberts 1,*, Kuinian Li 2 1 Dams
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 informationFracture Behaviour of FRP Cross-Ply Laminate With Embedded Delamination Subjected To Transverse Load
Fracture Behaviour of FRP Cross-Ply Laminate With Embedded Delamination Subjected To Transverse Load Sriram Chintapalli 1, S.Srilakshmi 1 1 Dept. of Mech. Engg., P. V. P. Siddhartha Institute of Technology.
More information