The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems. Prof. Dr. Eleni Chatzi Lecture 1-20 September, 2017
|
|
- Jocelyn Robinson
- 6 years ago
- Views:
Transcription
1 The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Prof. Dr. Eleni Chatzi Lecture 1-20 September, 2017 Institute of Structural Engineering Method of Finite Elements II 1
2 Course Information Instructors Prof. Dr. Eleni Chatzi, Dr. Giuseppe Abbiati, Dr. Konstantinos Agathos, Assistant Adrian Egger, HIL E19.5, Course Website Lecture Notes and Demos will be posted in: Performance Evaluation - Final Project (100%) Topic Categories announced on You are free to shape your own project! Computer Labs are already be pre-announced and it is advised to bring a laptop along for those sessions Help us Structure the Course! Participate in our online survey: Institute of Structural Engineering Method of Finite Elements II 2
3 Course Information Suggested Reading Nonlinear Finite Elements for Continua and Structures, by T. Belytschko, W. K. Liu, and B. Moran, John Wiley and Sons, 2000 The Finite Element Method: Linear Static and Dynamic Finite Element Analysis, by T. J. R. Hughes, Dover Publications, 2000 Nonlinear finite element analysis of solids and structures,crisfield, M.A., Remmers, J.J. and Verhoosel, C.V., John Wiley Sons, 2012 Computational methods for plasticity: theory and applications, De Souza Neto, E.A., Peric, D. and Owen, D.R., 2011 Lecture Notes by Carlos A. Felippa Nonlinear Finite Element Methods (ASEN 6107): colorado.edu/engineering/cas/courses.d/nfem.d/home.html Institute of Structural Engineering Method of Finite Elements II 3
4 Learning Goals By the end of the lecture, the goal is to learn: When is a problem nonlinear & what are types of nonlinearity? What types of nonlinearity exist? What are the differences, but also similarities, between linear and nonlinear problems? How to formulate & solve nonlinear problems What algorithms to use & what are typical problems? Institute of Structural Engineering Method of Finite Elements II 4
5 Learning Goals What we know so far We can recognize linearity As the process following the fundamental rule output = constant input F M=FL When is this linear? M w L θ when M = EI = EI φ EI φ Institute of Structural Engineering Method of Finite Elements II 5
6 Learning Goals Beyond Linearity F M=FL What if? M plasticity w L θ We have excessive stresses in a given section? Material Nonlinearity EI φ F M=FL What if? w θ We have excessive stresses in a given section? L Fracture Phenomena Institute of Structural Engineering Method of Finite Elements II 6
7 Learning Goals Beyond Linearity F w F θ M=FL & how does the structure actually deform? M =? L Large Displacements ΔL w F θ M=FL & how does the structure actually deform when heated? L Thermomechanics Problems khanacademy.org Institute of Structural Engineering Method of Finite Elements II 7
8 Significance of the Lecture The importance of Modeling Models are abstract representations of the real world Simulation models are extensively used in civil engineering practice. Such models allow the user to understand structural system performance predict structural behavior diagnose damage optimize design, etc Institute of Structural Engineering Method of Finite Elements II 8
9 Significance of the Lecture The importance of Modeling Models are abstract representations of the real world Structures may fail because The models used for their design do not properly represent the complexity of the physics Overly simplified models may not predict response under extreme loads The values of the input parameters have not been selected properly Institute of Structural Engineering Method of Finite Elements II 9
10 Course Outline FEM II - Lecture Schedule Introduction - From Linear to Nonlinear Considerations Block #1 Material Nonlinearity I - Constitutive Modeling Material Nonlinearity II - FE Implementation Computer Lab I Block #2 Large Deformations I Block #3 Dynamics I Large Deformations II Nonlinear Dynamics II Institute of Structural Engineering Method of Finite Elements II 10
11 Course Outline FEM II - Lecture Schedule Computer Lab II Block #4 Fracture I - General principles Fracture II - The extended Finite Element Method Block #5 Thermomechanics I Special Topics - Coupling/Penalty Methods (Lagrange Multipliers), Contact, Imposition of Essential BCs Computer Lab III Institute of Structural Engineering Method of Finite Elements II 11
12 Review of the Finite Element Method (FEM) FEM: Numerical Technique for approximating the solution of continuous systems. The MFE is based on the physical discretization of the observed domain, thus reducing the number of the degrees of freedom; and shifting from an analytical to a numerical formulation. Continuous Discrete q y+dy h1 q x dy q x+dx dx q y h2 L Flow of water Permeable Soil ( ) k 2 φ 2 x + 2 φ 2 y Impermeable Rock = 0 Laplace Equation F = KX Direct Stiffness Method Institute of Structural Engineering Method of Finite Elements II 12
13 Review of the Linear Finite Element Method (FEM) How is the Physical Problem formulated? A linear elastic system under specific loads and constraints at its boundaries can be described by the following fundamental equations: Linear balance (Equilibrium) equation. Linear constitutive model. Kinematic Equation, or Static balance equation or, if we consider modal analysis as an equivalent static analysis, at least stationary. Example of a bridge FE model. Institute of Structural Engineering Method of Finite Elements II 13
14 Review of the Linear Finite Element Method (FEM) How is the Physical Problem formulated? A linear elastic system under specific loads and constraints at its boundaries can be described by the following fundamental equations: Linear balance (Equilibrium) equation. Linear constitutive model. Kinematic Equation, or Static balance equation or, if we consider modal analysis as an equivalent static analysis, at least stationary. Example of a bridge FE model. Institute of Structural Engineering Method of Finite Elements II 13
15 Review of the Linear Finite Element Method (FEM) How is the Physical Problem formulated? A linear elastic system under specific loads and constraints at its boundaries can be described by the following fundamental equations: Linear balance (Equilibrium) equation. Linear constitutive model. Kinematic Equation, or Static balance equation or, if we consider modal analysis as an equivalent static analysis, at least stationary. Example of a bridge FE model. Institute of Structural Engineering Method of Finite Elements II 13
16 Review of the Linear Finite Element Method (FEM) The Strong Form of a Linear Problem On the basis of the previous fundamental euations, a governing differential equation, also known as the strong form of the problem may be extracted: 1 Equilibrium Equations ex. f (x) = σa = q(l x) dσ dx + q/a = 0 2 Constitutive Model ex. σ = Eɛ 3 Kinematics Relationship ex. ɛ = du dx Axial bar Example q x q f(x) Institute of Structural Engineering Method of Finite Elements II 14 L-x
17 Examples of Linear Balance Equations Mechanical Axial bar Example q Thermomechanical θ1 θ(x) x φ(x)>0 θ2 Force balance: dσ dx + P = 0 u : displacement [m] σ : stress [Pa] P : body load [ ] N m 3 Heat flux balance: dφ dx + R = 0 θ : temperature [K] φ : heat flux [ W m 2 ] R : heat source [ W m 3 ] Institute of Structural Engineering Method of Finite Elements II 15
18 Examples of Linear Constitutive Models Mechanical Thermomechanical Hooke s law: σ = E du dx u : displacement [m] x : space coordinate [m] du dx : strain [ ] m m σ : stress [Pa] E : Young modulus [Pa] Fourier s law: φ = k dθ dx θ : temperature [K] x : space coordinate [m] dθ dx : temperature grad. [ ] K m φ : heat flux [ ] W m 2 k : thermal cond. [ ] W mk Institute of Structural Engineering Method of Finite Elements II 16
19 Examples of Problems Equations Mechanical Thermomechanical Force balance: E d 2 u dx 2 + P = 0 u : displacement [m] E : Young modulus [Pa] P : body load [ N m 3 ] Heat flux balance: k d 2 θ dx 2 + R = 0 θ : temperature [K] k : thermal cond. [ ] W mk R : heat source [ ] W m 3 Institute of Structural Engineering Method of Finite Elements II 17
20 Boundary Value Problems (BVPs) The Strong Form Mechanical BVP Thermal BVP E d 2 u dx 2 + P = 0 { u 0,L = u 0,L du dx 0,L = du 0,L dx u : displacement [m] E : Young modulus [Pa] P : body load [ N m 3 ] k d 2 θ dx 2 + R = 0 { θ 0,L = θ 0,L dθ dx 0,L = dθ 0,L dx θ : temperature [K] k : thermal cond. [ ] W mk R : heat source [ ] W m 3 Institute of Structural Engineering Method of Finite Elements II 18
21 Boundary Value Problems (BVPs) The Strong Form Mechanical BVP Thermal BVP E d 2 u dx 2 + P = 0 { u 0,L = u 0,L du dx 0,L = du 0,L dx FEM Formulation: K u n u n u u = F u n u 1 n u 1 k d 2 θ dx 2 + R = 0 { θ 0,L = θ 0,L dθ dx 0,L = dθ 0,L dx FEM Formulation: K θ n θ n θ θ = F θ n θ 1 n θ 1 Institute of Structural Engineering Method of Finite Elements II 18
22 Boundary Value Problems (BVPs) The Strong Form Mechanical BVP Thermal BVP E d 2 u dx 2 + P = 0 { u 0,L = u 0,L du dx 0,L = du 0,L dx FEM Formulation: K u n u n u u = F u n u 1 n u 1 k d 2 θ dx 2 + R = 0 { θ 0,L = θ 0,L dθ dx 0,L = dθ 0,L dx FEM Formulation: K θ n θ n θ θ = F θ n θ 1 n θ 1 Institute of Structural Engineering Method of Finite Elements II 18
23 Weak Form - 1D FEM So how to solve the Strong Form The strong form requires strong continuity on the dependent field variables (usually displacements). Whatever functions define these variables have to be differentiable up to the order of the PDE that exist in the strong form of the system equations. Obtaining the exact solution for a strong form of the system equation is a quite difficult task for practical engineering problems. The finite difference method can be used to solve the system equations of the strong form and obtain an approximate solution. However, this method usually works well for problems with simple and regular geometry and boundary conditions. Alternatively we can use the finite element method on a weak form of the system. This weak form is usually obtained through energy principles which is why it is also known as variational form. Institute of Structural Engineering Method of Finite Elements II 19
24 Weak Form - 1D FEM From Strong Form to Weak form Three are the approaches commonly used to go from strong to weak form: Principle of Virtual Work Principle of Minimum Potential Energy Methods of weighted residuals (Galerkin, Collocation, Least Squares methods, etc) Visit our Course page on Linear FEM for further details: method-of-finite-elements-i.html Institute of Structural Engineering Method of Finite Elements II 20
25 Weak Form - 1D FEM From Strong Form to Weak form - Approach #1 Principle of Virtual Work For any set of compatible small virtual displacements imposed on the body in its state of equilibrium, the total internal virtual work is equal to the total external virtual work. The principle of virtual displacements facilitates the derivation of discretized equations: where Ω R u (u) = F u δɛ T σdω = Ω δu T P vol u dω + δu T P sur u dγ Γ σ : stress state generated by volume P vol u surface P sur u loads. δu and δɛ : compatible variations of displacement u and strain ɛ fields. and Institute of Structural Engineering Method of Finite Elements II 21
26 Weak Form - 1D FEM From Strong Form to Weak form - Approach #1 Principle of Virtual Work For any set of compatible small virtual displacements imposed on the body in its state of equilibrium, the total internal virtual work is equal to the total external virtual work. The principle of virtual displacements facilitates the derivation of discretized equations: where Ω R u (u) = F u δɛ T σdω = Ω δu T P vol u dω + δu T P sur u dγ Γ σ : stress state generated by volume P vol u surface P sur u loads. δu and δɛ : compatible variations of displacement u and strain ɛ fields. and Institute of Structural Engineering Method of Finite Elements II 21
27 Weak Form - 1D FEM From Strong Form to Weak form - Approach #1 Principle of Virtual Work For any set of compatible small virtual displacements imposed on the body in its state of equilibrium, the total internal virtual work is equal to the total external virtual work. The principle of virtual displacements facilitates the derivation of discretized equations: where Ω R u (u) = F u δɛ T σdω = Ω δu T P vol u dω + δu T P sur u dγ Γ σ : stress state generated by volume P vol u surface P sur u loads. δu and δɛ : compatible variations of displacement u and strain ɛ fields. and Institute of Structural Engineering Method of Finite Elements II 21
28 Weak Form - 1D FEM From Strong Form to Weak form - Approach #3 Consider the following (1D) generic strong form: ( d c(x) du ) + f(x) = 0 dx dx c(0) d dx u(0) = C 1 (Neumann BC) & u(l) = 0 (Dirichlet BC) Galerkin s Method Given an arbitrary weight function w, where S = {u u C 0, u(l) = 0}, S 0 = {w w C 0, w(l) = 0} C 0 is the collection of all continuous functions. Multiplying by w and l integrating over Ω w(x)[(c(x)u (x)) + f (x)]dx = 0 0 [w(0)(c(0)u (0) + C 1 ] = 0 Institute of Structural Engineering Method of Finite Elements II 22
29 Weak Form - 1D FEM Using the divergence theorem (integration by parts) we reduce the order of the differential: l 0 l wg dx = [wg] l 0 gw dx 0 The weak form is then reduced to the following problem. Also, in what follows we assume constant properties c(x) = c = const. Find u(x) S such that: l 0 w cu dx = l 0 wfdx + w(0)c 1 S = {u u C 0, u(l) = 0} S 0 = {w w C 0, w(l) = 0} Institute of Structural Engineering Method of Finite Elements II 23
30 Weak Form Notes: 1 Natural (Neumann) boundary conditions, are imposed on the secondary variables like forces and tractions. For example, u y (x 0, y 0 ) = u 0. 2 Essential (Dirichlet) or geometric boundary conditions, are imposed on the primary variables like displacements. For example, u(x 0, y 0 ) = u 0. 3 A solution to the strong form will also satisfy the weak form, but not vice versa. Since the weak form uses a lower order of derivatives it can be satisfied by a larger set of functions. 4 For the derivation of the weak form we can choose any weighting function w, since it is arbitrary, so we usually choose one that satisfies homogeneous boundary conditions wherever the actual solution satisfies essential boundary conditions. Note that this does not hold for natural boundary conditions! Institute of Structural Engineering Method of Finite Elements II 24
31 FE formulation: Reminder Let s recall the steps to the Galerking approach 1 Step #1: - Discretization 2 Step #2: - Aprroximate u(x) 3 Step #3: - Aprroximate w(x) 4 Step #4: - Substitute in the weak form w(x) 5 Step #5: - Obtain the discrete system matrices Institute of Structural Engineering Method of Finite Elements II 25
32 FE formulation: Discretization How to derive a solution to the weak form? Step #1:Follow the FE approach: Divide the body into finite elements, e, connected to each other through nodes. e x 1 e x 2 e Then break the overall integral into a summation over the finite elements: [ x e ] 2 x e w cu 2 dx wfdx w(0)c 1 = 0 x1 e x1 e e Institute of Structural Engineering Method of Finite Elements II 26
33 1D FE formulation: Galerkin Method Step #2: Approximate the continuous displacement field u(x) using a discrete equivalent: Galerkin s method assumes that the approximate (or trial) solution, u, can be expressed as a linear combination of the nodal point displacements u i, where i refers to the corresponding node number. u(x) u h (x) = i N i (x)u i = N(x)u where bold notation signifies a vector and N i (x) are the shape functions. In fact, the shape function can be any mathematical formula that helps us infer (interpolate) what happens at points that lie within the nodes of the mesh. In the 1-D case that we are using as a reference, N i (x) are defined as 1st degree polynomials indicating a linear interpolation. Institute of Structural Engineering Method of Finite Elements II 27
34 1D FE formulation: Galerkin s Method Shape function Properties: Bounded and Continuous One for each node Ni e(x j e) = δ ij, where { 1 if i = j δ ij = 0 if i j The shape functions can be written as piecewise functions of the x coordinate: This is not a convenient notation. x x i 1 Instead of using the global coordinate, x i 1 x < x i x i x i 1 x, things become simplified when N i (x) = xi + 1 x, x i x < x using coordinate ξ referring to the i+1 x i+1 x i local system of the element (see page 0, otherwise 25). Institute of Structural Engineering Method of Finite Elements II 28
35 1D FE formulation: Galerkin s Method Step #3: Approximate w(x) using a discrete equivalent: The weighting function, w is usually (although not necessarily) chosen to be of the same form as u w(x) w h (x) = i N i (x)w i = N(x)w i.e. for 2 nodes: N = [N 1 N 2 ] u = [u 1 u 2 ] T w = [w 1 w 2 ] T Alternatively we could have a Petrov-Galerkin formulation, where w(x) is obtained through the following relationships: w(x) = i (N i + δ he σ δ = coth( Pee 2 ) 2 Pe e dn i dx )w i coth = ex + e x e x e x Institute of Structural Engineering Method of Finite Elements II 29
36 1D FE formulation: Galerkin s Method Step #4: Substituting into the weak formulation and rearranging terms we obtain the following in matrix notation: l 0 l 0 w cu dx l (w T N T ) c(nu) dx 0 wfdx w(0)c 1 = 0 l 0 w T N T fdx w T N(0) T C 1 = 0 Since w, u are vectors, each one containing a set of discrete values corresponding at the nodes i, it follows that the above set of equations can be rewritten in the following form, i.e. as a summation over the w i, u i components (Einstein notation): l 0 ( l f 0 j i ) dn i (x) u i c dx j w j N j (x)dx j dn j (x) w j dx dx w j N j (x)c 1 = 0 x=0 Institute of Structural Engineering Method of Finite Elements II 30
37 1D FE formulation: Galerkin s Method This is rewritten as, [ ( l w j 0 j i 0 i cu i dn i (x) dx ) ] dn j (x) fn j (x)dx + (N j (x)c 1 ) x=0 = 0 dx The above equation has to hold w j since the weighting function w(x) is an arbitrary one. Therefore the following system of equations has to hold: ( ) l dn i (x) dn j (x) cu i fn j (x)dx + (N j (x)c 1 ) dx dx x=0 = 0 j = 1,..., n After reorganizing and moving the summation outside the integral, this becomes: [ ] l c dn i(x) dn j (x) l u i = fn j (x)dx + (N j (x)c 1 ) dx dx x=0 = 0 j = 1,..., n 0 i 0 Institute of Structural Engineering Method of Finite Elements II 31
38 1D FE formulation: Galerkin s Method We finally obtain the following discrete system in matrix notation: Ku = f where writing the integral from 0 to l as a summation over the subelements we obtain: K = A e K e K e = x e 2 x1 e [ ] dn(x) T c dn(x) dx = dx dx x e 2 x e 1 B T cbdx f = A e f e f e = x e 2 x e 1 N T fdx + N T h x=0 where A is not a sum but an assembly, and, x denotes differentiation with respect to x. In addition, B = dn(x) is known as the strain-displacement matrix. dx Institute of Structural Engineering Method of Finite Elements II 32
39 Transitioning to the Nonlinear Case How to treat the Nonlinear Case? Linear Problems K d= F or Pd ( ) = F Stiffness matrix K is constant Pd ( + d) = Pd ( ) + Pd ( ) P( α d) = α P( d) = αf If the load is doubled, displacement is doubled, too Superposition is possible 2F F d K 2d Nonlinear Problems F Pd ( ) = F, P(2 d) 2F K T How to find d for a given F? d i source: Nam-Ho Kim, EGM6352, Lecture Notes Institute of Structural Engineering Method of Finite Elements II 33 d
Chapter 3 Variational Formulation & the Galerkin Method
Institute of Structural Engineering Page 1 Chapter 3 Variational Formulation & the Galerkin Method Institute of Structural Engineering Page 2 Today s Lecture Contents: Introduction Differential formulation
More informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Prof. Dr. Eleni Chatzi Dr. Giuseppe Abbiati, Dr. Konstantinos Agathos Lecture 1-21 September, 2017 Institute of Structural Engineering
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 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 informationIntroduction to Finite Element Method
Introduction to Finite Element Method Dr. Rakesh K Kapania Aerospace and Ocean Engineering Department Virginia Polytechnic Institute and State University, Blacksburg, VA AOE 524, Vehicle Structures Summer,
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 informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Thermomechanics
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Thermomechanics Prof. Dr. Eleni Chatzi Dr. Giuseppe Abbiati, Dr. Konstantinos Agathos Lecture 13-14 December, 2017 1 / 30 Forewords
More informationPlates and Shells: Theory and Computation. Dr. Mostafa Ranjbar
Plates and Shells: Theory and Computation Dr. Mostafa Ranjbar Outline -1-! This part of the module consists of seven lectures and will focus on finite elements for beams, plates and shells. More specifically,
More informationCIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen
CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall
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 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 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 informationLecture 15: Revisiting bars and beams.
3.10 Potential Energy Approach to Derive Bar Element Equations. We place assumptions on the stresses developed inside the bar. The spatial stress and strain matrices become very sparse. We add (ad-hoc)
More informationMathematical Methods - Lecture 9
Mathematical Methods - Lecture 9 Yuliya Tarabalka Inria Sophia-Antipolis Méditerranée, Titane team, http://www-sop.inria.fr/members/yuliya.tarabalka/ Tel.: +33 (0)4 92 38 77 09 email: yuliya.tarabalka@inria.fr
More informationFinite Element Method-Part II Isoparametric FE Formulation and some numerical examples Lecture 29 Smart and Micro Systems
Finite Element Method-Part II Isoparametric FE Formulation and some numerical examples Lecture 29 Smart and Micro Systems Introduction Till now we dealt only with finite elements having straight edges.
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationMulti-Point Constraints
Multi-Point Constraints Multi-Point Constraints Multi-Point Constraints Single point constraint examples Multi-Point constraint examples linear, homogeneous linear, non-homogeneous linear, homogeneous
More informationSteps in the Finite Element Method. Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen
Steps in the Finite Element Method Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen General Idea Engineers are interested in evaluating effects such as deformations, stresses,
More informationMECh300H Introduction to Finite Element Methods. Finite Element Analysis (F.E.A.) of 1-D Problems
MECh300H Introduction to Finite Element Methods Finite Element Analysis (F.E.A.) of -D Problems Historical Background Hrenikoff, 94 frame work method Courant, 943 piecewise polynomial interpolation Turner,
More informationBAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS
Journal of Computational and Applied Mechanics, Vol.., No. 1., (2005), pp. 83 94 BAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS Vladimír Kutiš and Justín Murín Department
More informationVirtual Work and Variational Principles
Virtual Work and Principles Mathematically, the structural analysis problem is a boundary value problem (BVP). Forces, displacements, stresses, and strains are connected and computed within the framework
More information1 Nonlinear deformation
NONLINEAR TRUSS 1 Nonlinear deformation When deformation and/or rotation of the truss are large, various strains and stresses can be defined and related by material laws. The material behavior can be expected
More informationContents. Prologue Introduction. Classical Approximation... 19
Contents Prologue........................................................................ 15 1 Introduction. Classical Approximation.................................. 19 1.1 Introduction................................................................
More informationEshan V. Dave, Secretary of M&FGM2006 (Hawaii) Research Assistant and Ph.D. Candidate. Glaucio H. Paulino, Chairman of M&FGM2006 (Hawaii)
Asphalt Pavement Aging and Temperature Dependent Properties through a Functionally Graded Viscoelastic Model I: Development, Implementation and Verification Eshan V. Dave, Secretary of M&FGM2006 (Hawaii)
More informationPLAXIS. Scientific Manual
PLAXIS Scientific Manual 2016 Build 8122 TABLE OF CONTENTS TABLE OF CONTENTS 1 Introduction 5 2 Deformation theory 7 2.1 Basic equations of continuum deformation 7 2.2 Finite element discretisation 8 2.3
More informationGalerkin Finite Element Model for Heat Transfer
Galerkin Finite Element Model for Heat Transfer Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Table of Contents 1 Notation remarks 1 2 Local differential
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 informationSlender Structures Load carrying principles
Slender Structures Load carrying principles Basic cases: Extension, Shear, Torsion, Cable Bending (Euler) v017-1 Hans Welleman 1 Content (preliminary schedule) Basic cases Extension, shear, torsion, cable
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 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 informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 3: Finite Elements in 2-D Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods 1 / 18 Outline 1 Boundary
More informationIntroduction to Finite Element Method. Dr. Aamer Haque
Introduction to Finite Element Method 4 th Order Beam Equation Dr. Aamer Haque http://math.iit.edu/~ahaque6 ahaque7@iit.edu Illinois Institute of Technology July 1, 009 Outline Euler-Bernoulli Beams Assumptions
More informationContent. Department of Mathematics University of Oslo
Chapter: 1 MEK4560 The Finite Element Method in Solid Mechanics II (January 25, 2008) (E-post:torgeiru@math.uio.no) Page 1 of 14 Content 1 Introduction to MEK4560 3 1.1 Minimum Potential energy..............................
More informationIntroduction to fracture mechanics
Introduction to fracture mechanics Prof. Dr. Eleni Chatzi Dr. Giuseppe Abbiati, Dr. Konstantinos Agathos Lecture 6-9 November, 2017 Institute of Structural Engineering, ETH Zu rich November 9, 2017 Institute
More informationa x Questions on Classical Solutions 1. Consider an infinite linear elastic plate with a hole as shown. Uniform shear stress
Questions on Classical Solutions. Consider an infinite linear elastic plate with a hole as shown. Uniform shear stress σ xy = T is applied at infinity. Determine the value of the stress σ θθ on the edge
More informationAlternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE. University of Liège Aerospace & Mechanical Engineering
University of Liège Aerospace & Mechanical Engineering Alternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE Van Dung NGUYEN Innocent NIYONZIMA Aerospace & Mechanical engineering
More informationMechanics of Structures (CE130N) Lab 3
UNIVERSITY OF CALIFORNIA AT BERKELEY CE 130N, Spring 2009 Department of Civil and Environmental Engineering Prof. S. Govindjee and Dr. T. Koyama Structural Engineering, Mechanics and Materials Lab 3 1
More informationMechanical Design in Optical Engineering
OPTI Buckling Buckling and Stability: As we learned in the previous lectures, structures may fail in a variety of ways, depending on the materials, load and support conditions. We had two primary concerns:
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
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 informationLecture 27 Introduction to finite elements methods
Fall, 2017 ME 323 Mechanics of Materials Lecture 27 Introduction to finite elements methods Reading assignment: News: Instructor: Prof. Marcial Gonzalez Last modified: 10/24/17 7:02:00 PM Finite element
More informationStructural Analysis of Truss Structures using Stiffness Matrix. Dr. Nasrellah Hassan Ahmed
Structural Analysis of Truss Structures using Stiffness Matrix Dr. Nasrellah Hassan Ahmed FUNDAMENTAL RELATIONSHIPS FOR STRUCTURAL ANALYSIS In general, there are three types of relationships: Equilibrium
More informationC:\Users\whit\Desktop\Active\304_2012_ver_2\_Notes\4_Torsion\1_torsion.docx 6
C:\Users\whit\Desktop\Active\304_2012_ver_2\_Notes\4_Torsion\1_torsion.doc 6 p. 1 of Torsion of circular bar The cross-sections rotate without deformation. The deformation that does occur results from
More informationGeometric Misfitting in Structures An Interval-Based Approach
Geometric Misfitting in Structures An Interval-Based Approach M. V. Rama Rao Vasavi College of Engineering, Hyderabad - 500 031 INDIA Rafi Muhanna School of Civil and Environmental Engineering Georgia
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationCHAPTER 5. Beam Theory
CHPTER 5. Beam Theory SangJoon Shin School of Mechanical and erospace Engineering Seoul National University ctive eroelasticity and Rotorcraft Lab. 5. The Euler-Bernoulli assumptions One of its dimensions
More informationConcepts. 3.1 Numerical Analysis. Chapter Numerical Analysis Scheme
Chapter 3 Concepts The objective of this work is to create a framework to implement multi-disciplinary finite element applications. Before starting, it is necessary to explain some basic concepts of the
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 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 information2C9 Design for seismic and climate changes. Jiří Máca
2C9 Design for seismic and climate changes Jiří Máca List of lectures 1. Elements of seismology and seismicity I 2. Elements of seismology and seismicity II 3. Dynamic analysis of single-degree-of-freedom
More informationWeighted Residual Methods
Weighted Residual Methods Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy of Sciences
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 informationUnit M1.5 Statically Indeterminate Systems
Unit M1.5 Statically Indeterminate Systems Readings: CDL 2.1, 2.3, 2.4, 2.7 16.001/002 -- Unified Engineering Department of Aeronautics and Astronautics Massachusetts Institute of Technology LEARNING OBJECTIVES
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 information2.20 Fall 2018 Math Review
2.20 Fall 2018 Math Review September 10, 2018 These notes are to help you through the math used in this class. This is just a refresher, so if you never learned one of these topics you should look more
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 informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems. Prof. Dr. Eleni Chatzi Lecture ST1-19 November, 2015
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Prof. Dr. Eleni Chatzi Lecture ST1-19 November, 2015 Institute of Structural Engineering Method of Finite Elements II 1 Constitutive
More informationFINITE VOLUME METHOD: BASIC PRINCIPLES AND EXAMPLES
FINITE VOLUME METHOD: BASIC PRINCIPLES AND EXAMPLES SHRUTI JAIN B.Tech III Year, Electronics and Communication IIT Roorkee Tutors: Professor G. Biswas Professor S. Chakraborty ACKNOWLEDGMENTS I would like
More informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
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 informationUniversity of Groningen
University of Groningen Nature-inspired microfluidic propulsion using magnetic actuation Khaderi, S. N.; Baltussen, M. G. H. M.; Anderson, P. D.; Ioan, D.; den Toonder, J.M.J.; Onck, Patrick Published
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 informationFundamental Solutions and Green s functions. Simulation Methods in Acoustics
Fundamental Solutions and Green s functions Simulation Methods in Acoustics Definitions Fundamental solution The solution F (x, x 0 ) of the linear PDE L {F (x, x 0 )} = δ(x x 0 ) x R d Is called the fundamental
More informationLecture 12: Finite Elements
Materials Science & Metallurgy Part III Course M6 Computation of Phase Diagrams H. K. D. H. Bhadeshia Lecture 2: Finite Elements In finite element analysis, functions of continuous quantities such as temperature
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 information1. Let a(x) > 0, and assume that u and u h are the solutions of the Dirichlet problem:
Mathematics Chalmers & GU TMA37/MMG800: Partial Differential Equations, 011 08 4; kl 8.30-13.30. Telephone: Ida Säfström: 0703-088304 Calculators, formula notes and other subject related material are not
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 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 informationResidual Force Equations
3 Residual Force Equations NFEM Ch 3 Slide 1 Total Force Residual Equation Vector form r(u,λ) = 0 r = total force residual vector u = state vector with displacement DOF Λ = array of control parameters
More informationDevelopment of Truss Equations
CIVL 7/87 Chapter 3 - Truss Equations - Part /53 Chapter 3a Development of Truss Equations Learning Objectives To derive the stiffness matri for a bar element. To illustrate how to solve a bar assemblage
More informationDue Monday, September 14 th, 12:00 midnight
Due Monday, September 14 th, 1: midnight This homework is considering the analysis of plane and space (3D) trusses as discussed in class. A list of MatLab programs that were discussed in class is provided
More informationRitz Method. Introductory Course on Multiphysics Modelling
Ritz Method Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy of Sciences Warsaw Poland
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 informationENGN2340 Final Project: Implementation of a Euler-Bernuolli Beam Element Michael Monn
ENGN234 Final Project: Implementation of a Euler-Bernuolli Beam Element Michael Monn 12/11/13 Problem Definition and Shape Functions Although there exist many analytical solutions to the Euler-Bernuolli
More informationIndeterminate Analysis Force Method 1
Indeterminate Analysis Force Method 1 The force (flexibility) method expresses the relationships between displacements and forces that exist in a structure. Primary objective of the force method is to
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationOther state variables include the temperature, θ, and the entropy, S, which are defined below.
Chapter 3 Thermodynamics In order to complete the formulation we need to express the stress tensor T and the heat-flux vector q in terms of other variables. These expressions are known as constitutive
More informationFinite-Elements Method 2
Finite-Elements Method 2 January 29, 2014 2 From Applied Numerical Analysis Gerald-Wheatley (2004), Chapter 9. Finite-Elements Method 3 Introduction Finite-element methods (FEM) are based on some mathematical
More informationDevelopment of discontinuous Galerkin method for linear strain gradient elasticity
Development of discontinuous Galerkin method for linear strain gradient elasticity R Bala Chandran Computation for Design and Optimizaton Massachusetts Institute of Technology Cambridge, MA L. Noels* Aerospace
More informationFINITE ELEMENT METHOD: APPROXIMATE SOLUTIONS
FINITE ELEMENT METHOD: APPROXIMATE SOLUTIONS I Introduction: Most engineering problems which are expressed in a differential form can only be solved in an approximate manner due to their complexity. The
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 16
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Module - 01 Lecture - 16 In the last lectures, we have seen one-dimensional boundary value
More informationMITOCW MITRES2_002S10nonlinear_lec15_300k-mp4
MITOCW MITRES2_002S10nonlinear_lec15_300k-mp4 The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources
More informationElasto-plastic concrete beam analysis by 1-dimensional Finite Element Method
AALBORG UNIVERSITY MASTER S THESIS Elasto-plastic concrete beam analysis by 1-dimensional Finite Element Method Authors: Niels F. Overgaard Martin B. Andreasen Supervisors: Johan Clausen Lars V. Andersen
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 informationFinite Element Method
Finite Element Method Finite Element Method (ENGC 6321) Syllabus Objectives Understand the basic theory of the FEM Know the behaviour and usage of each type of elements covered in this course one dimensional
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More informationModelling mechanical systems with nite elements
Modelling mechanical systems with nite elements Bastian Kanning February 26, 2010 Outline 1 Introduction 2 Virtual displacements 3 Getting started The ansatz functions 4 Mass and stiness elementwise 5
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 informationChapter 12 Elastic Stability of Columns
Chapter 12 Elastic Stability of Columns Axial compressive loads can cause a sudden lateral deflection (Buckling) For columns made of elastic-perfectly plastic materials, P cr Depends primarily on E and
More information12. Stresses and Strains
12. Stresses and Strains Finite Element Method Differential Equation Weak Formulation Approximating Functions Weighted Residuals FEM - Formulation Classification of Problems Scalar Vector 1-D T(x) u(x)
More informationHigh order material point method
Literature study High order material point method February 1, 216 R.P.W.M. Tielen Technische Universiteit Delft HIGH ORDER MATERIAL POINT METHOD FEBRUARY 1, 216 by R.P.W.M. Tielen Supervisors: Ir. E.D.
More informationWeighted Residual Methods
Weighted Residual Methods Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Table of Contents Problem definition. oundary-value Problem..................
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 informationPSEUDO ELASTIC ANALYSIS OF MATERIAL NON-LINEAR PROBLEMS USING ELEMENT FREE GALERKIN METHOD
Journal of the Chinese Institute of Engineers, Vol. 27, No. 4, pp. 505-516 (2004) 505 PSEUDO ELASTIC ANALYSIS OF MATERIAL NON-LINEAR PROBLEMS USING ELEMENT FREE GALERKIN METHOD Raju Sethuraman* and Cherku
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Variational Problems of the Dirichlet BVP of the Poisson Equation 1 For the homogeneous
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationDHANALAKSHMI COLLEGE OF ENGINEERING, CHENNAI DEPARTMENT OF MECHANICAL ENGINEERING ME 6603 FINITE ELEMENT ANALYSIS PART A (2 MARKS)
DHANALAKSHMI COLLEGE OF ENGINEERING, CHENNAI DEPARTMENT OF MECHANICAL ENGINEERING ME 6603 FINITE ELEMENT ANALYSIS UNIT I : FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PART A (2 MARKS) 1. Write the types
More information