Inverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros

Size: px
Start display at page:

Download "Inverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros"

Transcription

1 Inverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros

2 Computational Design Forward design: direct manipulation of design parameters Level of abstraction Exploration of design spaces Inverse design: automatically infer design parameters from functional specifications Optimization-based

3 Computational Design of Actuated Deformable Characters, Skouras et al, 2013 Inverse Design Example Design an object that can deform like this: Need model that predicts the way in which objects deform Design parameters: material parameters, points of application for forces or constraints

4 A Lightweight Introduction to the Finite Element Method

5 The plan Start with a quick review of mass-spring systems Modeling deformable objects using continuum mechanics and the finite element method

6 Mass Spring Systems One way of modeling deformable objects is as a network of masses and springs Spring Mass

7 Mass Spring Systems 1 mass point, 1 spring l 0 l y i X y i x Spring length before/after l = x - y i l 0 = X - y i Deformation Measure Elastic Energy Forces Stiffness Direction Vector f int (x) = Deformation Measure

8 Mass Spring Systems Things to note: In the absence of external forces, elastic objects return to reference (undeformed) configuration Hyperelastic assumption: deformation energy depends only on deformed configuration: elastic forces are conservative

9 Mass Spring Systems The acceleration of a point mass is given by: ma = f int (x) + f ext Use standard methods to integrate this system forward in time to compute motion Or compute static equilibrium configuration Principle of minimum potential energy A mechanical system in static equilibrium will assume a state of minimum potential energy: find x to minimize W(x), or equivalently, such that f int (x) + f ext = 0

10 Mass Spring Systems f int (x) =

11 Mass Spring Systems: Simulation Pseudocode a0 //Array of accelerations v0 //Array of velocities p0 //Array of positions For each particle, p mmass of the particle fsum of all spring forces acting on the particle a[p]f/m End vintegrate(a) pintegrate(v)

12 Mass Spring Systems Can be used to model arbitrary deformable objects, but Behavior depends on tessellation Find good spring layout Find good spring constants Different types of springs interfere Limited accuracy No explicit volume or area preservation

13 Alternative Start from continuum mechanics Discretize with Finite Elements Decompose model into simple elements Setup & solve system of algebraic equations Advantages Accurate and controllable material behavior Largely independent of mesh structure

14 Mass Spring vs Continuum Mechanics Mass spring systems require: 1. Measure of Deformation 2. Material Model 3. Deformation Energy 4. Internal Forces Continuum Mechanics relies on conceptually identical pieces of information

15 Continuum Mechanics: 3D Deformations For a deformable body, identify: undeformed state 3 R described by positions deformed state 3 ' R described by positions Displacement field u describes ' u : ' x X u(x ) in terms of X x z u(x) ' X x u(x ) y x

16 Continuum Mechanics: 3D Deformations Consider material points 1 and X 2 and d X such that 2 X 1 is infinitesimal Now consider deformed vector x d' x x X u X ) X u( ) z 2 X 1 d 1 X 2 d d y X d 2 ( 2 1 X1 u( X1 d) u( X1) u( x1) ud u( x1) d' x 1 d' ' x 2 ( I u)d u Deformation x gradient F= X xu xv xw u y y y v w u z v z w z

17 So Displacement field transforms points Jacobian transforms directions x X u(x ) d' ( I u) d Fd Deformation x gradient F= X

18 Displacement Field and Deformation Gradient In general, displacement field is not explicitly described. Nevertheless, toy examples:

19 Displacement Field and Deformation Gradient In general, displacement field is not explicitly described. Nevertheless, toy examples:

20 Displacement Field and Deformation Gradient In general, displacement field is not explicitly described. Nevertheless, toy examples:

21 Displacement Field and Deformation Gradient In general, displacement field is not explicitly described. Nevertheless, toy examples:

22 Displacement Field and Deformation Gradient In general, displacement field is not explicitly described. Nevertheless, toy examples:

23 Measure of deformations Displacement field transforms points x X u(x ) Jacobian transforms directions d' ( I u) d Fd How can we describe deformation at given point?

24 Back to spring deformation Deformation measure (strain): Equivalent to F! Undeformed spring: = 1 Undeformed (infinitesimal) continuum volume: = I

25 Strain (description of deformation in terms of relative displacement) Deformation measure (strain): Desirable property: if spring is undeformed, strain is 0 Can we find a similar measure that would work for an arbitrary volume?

26 3D Nonlinear Strain Recall: Displacement field transforms points Jacobian (i.e. deformation gradient) transforms vectors x X u(x ) d' ( I u) d Fd Idea: to quantify change in shape, measure change in squared length for any vector d' 2 d 2 d' T d' T d d d T ( F T F I) d Green strain 1 T E ( F F I) 2

27 Strain (description of deformation in terms of relative displacement) Deformation measure (strain): Property: if spring is undeformed, strain is 0 Can we find a similar measure that would work for an arbitrary volume? Green strain 1 T E ( F F I) 2

28 3D Linear Strain Green strain is quadratic in displacements 2 1 w v w u w w v v u v w u v u u z z y z x y z y y x x z x y x ) ( ) ( u u u u I F F E T T T Neglecting quadratic terms leads to the linear Cauchy strain (small strain) Written out: Notation ε = 1 2 u + ut = 1 2 F + Ft I

29 3D Linear Strain Linear Cauchy strain Geometric interpretation y,v x,u i : normal strains i : shear strains 2 1 w v w u w w v v u v w u v u u z z y z x y z y y x x z x y x x y xy

30 Cauchy vs. Green strain Nonlinear Green strain is rotation-invariant Apply incremental rotation R to given deformation F to obtain F = RF Then E' 1 2 ( F' T F' I) 1 T T T ( F R RF I) ( F F I) E Linear Cauchy strain is not rotation-invariant ε = 1 2 F + F t 1 2 F + Ft = ε artifacts for larger rotations

31 Mass Spring vs Continuum Mechanics Mass spring systems: 1. Measure of Deformation 2. Material Model 3. Deformation Energy 4. Internal Forces Continuum Mechanics: 1. Measure of Deformation: Green or Cauchy strain 2. Material Model

32 Material Model: linear isotropic material Material model links strain to energy (and stress) Linear isotropic material (generalized Hooke s law) Energy density Ψ = 1 2 λtr(ε)2 +μtr(ε 2 ) tr ε = ε ii Lame parameters λ and μ are material constants related to Poisson Ratio and Young s modulus Interpretation tr ε 2 = ε F 2 penalizes all strain components equally tr ε 2 penalizes dilatations, i.e., volume changes

33 Volumetric Strain (dilatation) Consider a cube with side length a For a given deformation ε, the volumetric strain is ΔV/V0 = (a 1 + ε 11 a 1 + ε 22 a 1 + ε 33 a 3 )/a 3 = ε 11 + ε 22 + ε 33 + O ε 2 tr ε a

34 Young s Modulus and Poisson Ratio Energy density: Ψ = 1 2 λtr(ε)2 +μtr(ε 2 ) Lame parameters λ and μ are material constants related to the fundamental physical parameters: Poisson s Ratio and Young s modulus ( Young s modulus (E): measure of stiffness Poisson s ratio (ν): measure of volume preservation 10 Kg OBJECT E= GPa 10 Kg OBJECT E=60 KPa

35 Linear isotropic material Energy density: Ψ = 1 2 λtr(ε)2 +μtr(ε 2 ) Problem: Cauchy strain is not invariant under rotations artifacts for rotations Solutions: Corotational elasticity Nonlinear elasticity

36 Material Model: non-linear isotropic model Replace Cauchy strain with Green strain St. Venant-Kirchhoff material (StVK) Energy density: Ψ StVK = 1 2 λtr(e)2 +μtr(e 2 ) Rotation invariant!

37 Problems with StVK StVK softens under compression Energy Deformation Ψ StVK = 1 2 λtr(e)2 +μtr E 2

38 Advanced nonlinear materials Green Strain E = 1 2 Ft F I = 1 2 (C I) Split into deviatoric (i.e. shape changing) and volumetric deformations Volumetric: J = det F Deviatoric: C = det F 2/3 C Neo-Hookean material: Ψ NH = μ 2 tr C 3 μln (J) + λ 2 ln J 2 Energy

39 Different Models St. Venant-Kirchoff Neo-Hookean Linear

40 Mass Spring vs Continuum Mechanics Mass spring systems: 1. Measure of Deformation 2. Material Model 3. Deformation Energy 4. Internal Forces Continuum Mechanics: 1. Measure of Deformation: Green or Cauchy strain 2. Material Model: linear, StVK, Neo-Hookean 3. Energy Density (!)

41 Finite Element Discretization Divide domain into discrete elements, e.g., tetrahedra Explicitly store displacement values at nodes (x i ). Displacement field everywhere else obtained through interpolation: x X = N i X x i Compute Deformation Gradient: F = x X = N x i X i i X t

42 Basis Functions 4 nodes, undeformed positions X i, deformed positions x i basis functions N i X j : R 3 R Interpolate at vertices (delta-property) N i X j = δ ij Common choice: linear basis functions N i x, y, z = a i x + b i y + c i z + d i Compute N i (and N i X ) through X 4 V e X 3 X 1 Ω e X 2

43 From energy density to deformation energy Integrate energy density over the entire element: W e = Ψ F Ω e Deformation gradient F = If basis functions are linear: x X X = x i i N i X t X 4 F is linear in x i (nodal displacements) F is constant throughout element: V e X 3 W e = Ψ F = Ψ F V e Ω e X 1 Ω e X 2

44 Mass Spring vs Finite Element Method Mass spring systems: 1. Measure of Deformation 2. Material Model 3. Deformation Energy 4. Internal Forces Continuum Mechanics: 1. Measure of Deformation: Green or Cauchy strain 2. Material Model: linear, StVK, Neo-Hookean 3. Deformation Energy: integrate energy density 4. Internal Forces

45 An example Goal: find equilibrium configuration i.e., f i = 0 i Given x with f x 0, find Δx such that f x + Δx = 0 f x + Δx = f x + KΔx + O(Δx 2 ) Solve KΔx = f x for Δx Stiffness matrix K = f x

46 So, now you to can do this:

47 Computational Design of Actuated Deformable Characters, Skouras et al, 2013 Inverse Design Example Design an object that can deform like this: Need model that predicts the way in which objects deform Design parameters: material parameters, points of application for forces or constraints

48 Computational Design of Actuated Deformable Characters, Skouras et al, 2013 Inverse Design Example Design an object that can deform like this: FEM Model Design parameters: material parameters, points of application for forces or constraints

49 Pipeline + Input and Target Shapes Actuator Location Optimization Material Optimization Fabricated Deformable Model

50 Mathematical Formulation Design Parameters E x, p = E match x, x target subject to f total x, p = 0

51 Actuator Location Optimization

52 Material Distribution Optimization

53 Results Input Animation Fabricated Prototype

54 Results Rest Pose Target Pose Stiff Soft

55 Questions? Read paper for next class!

Continuum Mechanics and the Finite Element Method

Continuum Mechanics and the Finite Element Method Continuum Mechanics and the Finite Element Method 1 Assignment 2 Due on March 2 nd @ midnight 2 Suppose you want to simulate this The familiar mass-spring system l 0 l y i X y i x Spring length before/after

More information

Constitutive models. Constitutive model: determines P in terms of deformation

Constitutive models. Constitutive model: determines P in terms of deformation Constitutive models Constitutive model: determines P in terms of deformation Elastic material: P depends only on current F Hyperelastic material: work is independent of path strain energy density function

More information

Measurement of deformation. Measurement of elastic force. Constitutive law. Finite element method

Measurement of deformation. Measurement of elastic force. Constitutive law. Finite element method Deformable Bodies Deformation x p(x) Given a rest shape x and its deformed configuration p(x), how large is the internal restoring force f(p)? To answer this question, we need a way to measure deformation

More information

Finite Element Method in Geotechnical Engineering

Finite 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 information

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost

Game 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 information

Performance Evaluation of Various Smoothed Finite Element Methods with Tetrahedral Elements in Large Deformation Dynamic Analysis

Performance Evaluation of Various Smoothed Finite Element Methods with Tetrahedral Elements in Large Deformation Dynamic Analysis Performance Evaluation of Various Smoothed Finite Element Methods with Tetrahedral Elements in Large Deformation Dynamic Analysis Ryoya IIDA, Yuki ONISHI, Kenji AMAYA Tokyo Institute of Technology, Japan

More information

Fundamentals of Linear Elasticity

Fundamentals 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 information

Soft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies

Soft 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 information

3D Elasticity Theory

3D Elasticity Theory 3D lasticity Theory Many structural analysis problems are analysed using the theory of elasticity in which Hooke s law is used to enforce proportionality between stress and strain at any deformation level.

More information

Elements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004

Elements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004 Elements of Continuum Elasticity David M. Parks Mechanics and Materials II 2.002 February 25, 2004 Solid Mechanics in 3 Dimensions: stress/equilibrium, strain/displacement, and intro to linear elastic

More information

Surface force on a volume element.

Surface force on a volume element. STRESS and STRAIN Reading: Section. of Stein and Wysession. In this section, we will see how Newton s second law and Generalized Hooke s law can be used to characterize the response of continuous medium

More information

Mechanics PhD Preliminary Spring 2017

Mechanics 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 information

Macroscopic theory Rock as 'elastic continuum'

Macroscopic theory Rock as 'elastic continuum' Elasticity and Seismic Waves Macroscopic theory Rock as 'elastic continuum' Elastic body is deformed in response to stress Two types of deformation: Change in volume and shape Equations of motion Wave

More information

A Numerical Study of Finite Element Calculations for Incompressible Materials under Applied Boundary Displacements

A Numerical Study of Finite Element Calculations for Incompressible Materials under Applied Boundary Displacements A Numerical Study of Finite Element Calculations for Incompressible Materials under Applied Boundary Displacements A Thesis Submitted to the College of Graduate Studies and Research in Partial Fulfillment

More information

Mechanics of materials Lecture 4 Strain and deformation

Mechanics of materials Lecture 4 Strain and deformation Mechanics of materials Lecture 4 Strain and deformation Reijo Kouhia Tampere University of Technology Department of Mechanical Engineering and Industrial Design Wednesday 3 rd February, 206 of a continuum

More information

Mathematical Background

Mathematical Background CHAPTER ONE Mathematical Background This book assumes a background in the fundamentals of solid mechanics and the mechanical behavior of materials, including elasticity, plasticity, and friction. A previous

More information

1.050 Engineering Mechanics. Lecture 22: Isotropic elasticity

1.050 Engineering Mechanics. Lecture 22: Isotropic elasticity 1.050 Engineering Mechanics Lecture 22: Isotropic elasticity 1.050 Content overview I. Dimensional analysis 1. On monsters, mice and mushrooms 2. Similarity relations: Important engineering tools II. Stresses

More information

INVERSE ANALYSIS METHODS OF IDENTIFYING CRUSTAL CHARACTERISTICS USING GPS ARRYA DATA

INVERSE ANALYSIS METHODS OF IDENTIFYING CRUSTAL CHARACTERISTICS USING GPS ARRYA DATA Problems in Solid Mechanics A Symposium in Honor of H.D. Bui Symi, Greece, July 3-8, 6 INVERSE ANALYSIS METHODS OF IDENTIFYING CRUSTAL CHARACTERISTICS USING GPS ARRYA DATA M. HORI (Earthquake Research

More information

1 Nonlinear deformation

1 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 information

NDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16.

NDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16. CAVITY INSPECTION NDT&E Methods: UT VJ Technologies NDT&E Methods: UT 6. NDT&E: Introduction to Methods 6.1. Ultrasonic Testing: Basics of Elasto-Dynamics 6.2. Principles of Measurement 6.3. The Pulse-Echo

More information

MITOCW MITRES2_002S10nonlinear_lec15_300k-mp4

MITOCW 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 information

Linear Elasticity ( ) Objectives. Equipment. Introduction. ε is then

Linear Elasticity ( ) Objectives. Equipment. Introduction. ε is then Linear Elasticity Objectives In this lab you will measure the Young s Modulus of a steel wire. In the process, you will gain an understanding of the concepts of stress and strain. Equipment Young s Modulus

More information

You may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.

You may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator. MATHEMATICAL TRIPOS Part III Thursday 1 June 2006 1.30 to 4.30 PAPER 76 NONLINEAR CONTINUUM MECHANICS Attempt FOUR questions. There are SIX questions in total. The questions carry equal weight. STATIONERY

More information

Continuum mechanics V. Constitutive equations. 1. Constitutive equation: definition and basic axioms

Continuum mechanics V. Constitutive equations. 1. Constitutive equation: definition and basic axioms Continuum mechanics office Math 0.107 ales.janka@unifr.ch http://perso.unifr.ch/ales.janka/mechanics Mars 16, 2011, Université de Fribourg 1. Constitutive equation: definition and basic axioms Constitutive

More information

Comparison of Models for Finite Plasticity

Comparison of Models for Finite Plasticity Comparison of Models for Finite Plasticity A numerical study Patrizio Neff and Christian Wieners California Institute of Technology (Universität Darmstadt) Universität Augsburg (Universität Heidelberg)

More information

Variational principles in mechanics

Variational principles in mechanics CHAPTER Variational principles in mechanics.1 Linear Elasticity n D Figure.1: A domain and its boundary = D [. Consider a domain Ω R 3 with its boundary = D [ of normal n (see Figure.1). The problem of

More information

PEAT SEISMOLOGY Lecture 2: Continuum mechanics

PEAT SEISMOLOGY Lecture 2: Continuum mechanics PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a

More information

Fig. 1. Circular fiber and interphase between the fiber and the matrix.

Fig. 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 information

International Journal of Pure and Applied Mathematics Volume 58 No ,

International Journal of Pure and Applied Mathematics Volume 58 No , International Journal of Pure and Applied Mathematics Volume 58 No. 2 2010, 195-208 A NOTE ON THE LINEARIZED FINITE THEORY OF ELASTICITY Maria Luisa Tonon Department of Mathematics University of Turin

More information

ENGN 2340 Final Project Report. Optimization of Mechanical Isotropy of Soft Network Material

ENGN 2340 Final Project Report. Optimization of Mechanical Isotropy of Soft Network Material ENGN 2340 Final Project Report Optimization of Mechanical Isotropy of Soft Network Material Enrui Zhang 12/15/2017 1. Introduction of the Problem This project deals with the stress-strain response of a

More information

Basic Equations of Elasticity

Basic Equations of Elasticity A Basic Equations of Elasticity A.1 STRESS The state of stress at any point in a loaded bo is defined completely in terms of the nine components of stress: σ xx,σ yy,σ zz,σ xy,σ yx,σ yz,σ zy,σ zx,andσ

More information

Engineering Sciences 241 Advanced Elasticity, Spring Distributed Thursday 8 February.

Engineering Sciences 241 Advanced Elasticity, Spring Distributed Thursday 8 February. Engineering Sciences 241 Advanced Elasticity, Spring 2001 J. R. Rice Homework Problems / Class Notes Mechanics of finite deformation (list of references at end) Distributed Thursday 8 February. Problems

More information

2.1 Strain energy functions for incompressible materials

2.1 Strain energy functions for incompressible materials Chapter 2 Strain energy functions The aims of constitutive theories are to develop mathematical models for representing the real behavior of matter, to determine the material response and in general, to

More information

Introduction to Seismology Spring 2008

Introduction to Seismology Spring 2008 MIT OpenCourseWare http://ocw.mit.edu 12.510 Introduction to Seismology Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Stress and Strain

More information

Classification of Prostate Cancer Grades and T-Stages based on Tissue Elasticity Using Medical Image Analysis. Supplementary Document

Classification of Prostate Cancer Grades and T-Stages based on Tissue Elasticity Using Medical Image Analysis. Supplementary Document Classification of Prostate Cancer Grades and T-Stages based on Tissue Elasticity Using Medical Image Analysis Supplementary Document Shan Yang, Vladimir Jojic, Jun Lian, Ronald Chen, Hongtu Zhu, Ming C.

More information

4 NON-LINEAR ANALYSIS

4 NON-LINEAR ANALYSIS 4 NON-INEAR ANAYSIS arge displacement elasticity theory, principle of virtual work arge displacement FEA with solid, thin slab, and bar models Virtual work density of internal forces revisited 4-1 SOURCES

More information

CIVL4332 L1 Introduction to Finite Element Method

CIVL4332 L1 Introduction to Finite Element Method CIVL L Introduction to Finite Element Method CIVL L Introduction to Finite Element Method by Joe Gattas, Faris Albermani Introduction The FEM is a numerical technique for solving physical problems such

More information

Constitutive Equations

Constitutive Equations Constitutive quations David Roylance Department of Materials Science and ngineering Massachusetts Institute of Technology Cambridge, MA 0239 October 4, 2000 Introduction The modules on kinematics (Module

More information

Hervé DELINGETTE. INRIA Sophia-Antipolis

Hervé DELINGETTE. INRIA Sophia-Antipolis Hervé DELINGETTE INRIA Sophia-Antipolis French National Research Institute in Computer Science and Automatic Control Asclepios : 3D Segmentation, Simulation Platform, Soft Tissue Modeling http://www.inria.fr/asclepios

More information

Load Cell Design Using COMSOL Multiphysics

Load Cell Design Using COMSOL Multiphysics Load Cell Design Using COMSOL Multiphysics Andrei Marchidan, Tarah N. Sullivan and Joseph L. Palladino Department of Engineering, Trinity College, Hartford, CT 06106, USA joseph.palladino@trincoll.edu

More information

Analytical Mechanics: Elastic Deformation

Analytical Mechanics: Elastic Deformation Analytical Mechanics: Elastic Deformation Shinichi Hirai Dept. Robotics, Ritsumeikan Univ. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 1 / 60 Agenda Agenda

More information

Nonlinear Theory of Elasticity. Dr.-Ing. Martin Ruess

Nonlinear Theory of Elasticity. Dr.-Ing. Martin Ruess Nonlinear Theory of Elasticity Dr.-Ing. Martin Ruess geometry description Cartesian global coordinate system with base vectors of the Euclidian space orthonormal basis origin O point P domain of a deformable

More information

The Finite Element Method II

The Finite Element Method II [ 1 The Finite Element Method II Non-Linear finite element Use of Constitutive Relations Xinghong LIU Phd student 02.11.2007 [ 2 Finite element equilibrium equations: kinematic variables Displacement Strain-displacement

More information

SEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by

SEMM 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 information

Toward 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 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 information

Simulation of elasticity, biomechanics and virtual surgery. Joseph M. Teran

Simulation of elasticity, biomechanics and virtual surgery. Joseph M. Teran Contents Simulation of elasticity, biomechanics and virtual surgery Joseph M. Teran 1 Simulation of elasticity, biomechanics and virtual surgery 3 Introduction 3 Real-time computing 4 Continuum mechanics

More information

Mechanics of Earthquakes and Faulting

Mechanics of Earthquakes and Faulting Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Standard Solids and Fracture Fluids: Mechanical, Chemical Effects Effective Stress Dilatancy Hardening and Stability Mead, 1925

More information

Linear Cosserat elasticity, conformal curvature and bounded stiffness

Linear Cosserat elasticity, conformal curvature and bounded stiffness 1 Linear Cosserat elasticity, conformal curvature and bounded stiffness Patrizio Neff, Jena Jeong Chair of Nonlinear Analysis & Modelling, Uni Dui.-Essen Ecole Speciale des Travaux Publics, Cachan, Paris

More information

BAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS

BAR 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 information

Conservation of mass. Continuum Mechanics. Conservation of Momentum. Cauchy s Fundamental Postulate. # f body

Conservation of mass. Continuum Mechanics. Conservation of Momentum. Cauchy s Fundamental Postulate. # f body Continuum Mechanics We ll stick with the Lagrangian viewpoint for now Let s look at a deformable object World space: points x in the object as we see it Object space (or rest pose): points p in some reference

More information

**********************************************************************

********************************************************************** Department of Civil and Environmental Engineering School of Mining and Petroleum Engineering 3-33 Markin/CNRL Natural Resources Engineering Facility www.engineering.ualberta.ca/civil Tel: 780.492.4235

More information

FEM for elastic-plastic problems

FEM for elastic-plastic problems FEM for elastic-plastic problems Jerzy Pamin e-mail: JPamin@L5.pk.edu.pl With thanks to: P. Mika, A. Winnicki, A. Wosatko TNO DIANA http://www.tnodiana.com FEAP http://www.ce.berkeley.edu/feap Lecture

More information

Non-linear and time-dependent material models in Mentat & MARC. Tutorial with Background and Exercises

Non-linear and time-dependent material models in Mentat & MARC. Tutorial with Background and Exercises Non-linear and time-dependent material models in Mentat & MARC Tutorial with Background and Exercises Eindhoven University of Technology Department of Mechanical Engineering Piet Schreurs July 7, 2009

More information

Chapter 2. Rubber Elasticity:

Chapter 2. Rubber Elasticity: Chapter. Rubber Elasticity: The mechanical behavior of a rubber band, at first glance, might appear to be Hookean in that strain is close to 100% recoverable. However, the stress strain curve for a rubber

More information

Elements of Rock Mechanics

Elements of Rock Mechanics Elements of Rock Mechanics Stress and strain Creep Constitutive equation Hooke's law Empirical relations Effects of porosity and fluids Anelasticity and viscoelasticity Reading: Shearer, 3 Stress Consider

More information

EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading

EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading MA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading MA 3702 Mechanics & Materials Science Zhe Cheng (2018) 2 Stress & Strain - Axial Loading Statics

More information

Sparse monolithic compliant mechanisms using continuum structural topology optimization

Sparse monolithic compliant mechanisms using continuum structural topology optimization INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2005; 62:1579 1605 Published online 3 December 2004 in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/nme.1224

More information

1 Useful Definitions or Concepts

1 Useful Definitions or Concepts 1 Useful Definitions or Concepts 1.1 Elastic constitutive laws One general type of elastic material model is the one called Cauchy elastic material, which depend on only the current local deformation of

More information

A short review of continuum mechanics

A short review of continuum mechanics A short review of continuum mechanics Professor Anette M. Karlsson, Department of Mechanical ngineering, UD September, 006 This is a short and arbitrary review of continuum mechanics. Most of this material

More information

UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES

UNCONVENTIONAL 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 information

Constitutive models: Incremental plasticity Drücker s postulate

Constitutive models: Incremental plasticity Drücker s postulate Constitutive models: Incremental plasticity Drücker s postulate if consistency condition associated plastic law, associated plasticity - plastic flow law associated with the limit (loading) surface Prager

More information

Chapter 5 Structural Elements: The truss & beam elements

Chapter 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 information

Chapter 3: Stress and Equilibrium of Deformable Bodies

Chapter 3: Stress and Equilibrium of Deformable Bodies Ch3-Stress-Equilibrium Page 1 Chapter 3: Stress and Equilibrium of Deformable Bodies When structures / deformable bodies are acted upon by loads, they build up internal forces (stresses) within them to

More information

NONLINEAR CONTINUUM FORMULATIONS CONTENTS

NONLINEAR CONTINUUM FORMULATIONS CONTENTS NONLINEAR CONTINUUM FORMULATIONS CONTENTS Introduction to nonlinear continuum mechanics Descriptions of motion Measures of stresses and strains Updated and Total Lagrangian formulations Continuum shell

More information

Bending of Simply Supported Isotropic and Composite Laminate Plates

Bending of Simply Supported Isotropic and Composite Laminate Plates Bending of Simply Supported Isotropic and Composite Laminate Plates Ernesto Gutierrez-Miravete 1 Isotropic Plates Consider simply a supported rectangular plate of isotropic material (length a, width b,

More information

Mathematical Modeling. of Large Elastic-Plastic Deformations

Mathematical Modeling. of Large Elastic-Plastic Deformations Applied Mathematical Sciences, Vol. 8, 04, no. 60, 99 996 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/ams.04.447 Mathematical Modeling of Large Elastic-Plastic Deformations L. U. Sultanov Kazan

More information

University of Sheffield The development of finite elements for 3D structural analysis in fire

University of Sheffield The development of finite elements for 3D structural analysis in fire The development of finite elements for 3D structural analysis in fire Chaoming Yu, I. W. Burgess, Z. Huang, R. J. Plank Department of Civil and Structural Engineering StiFF 05/09/2006 3D composite structures

More information

Law of behavior very-rubber band: almost incompressible material

Law of behavior very-rubber band: almost incompressible material Titre : Loi de comportement hyperélastique : matériau pres[...] Date : 25/09/2013 Page : 1/9 Law of behavior very-rubber band: almost incompressible material Summary: One describes here the formulation

More information

ME FINITE ELEMENT ANALYSIS FORMULAS

ME 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 information

MITOCW MITRES2_002S10nonlinear_lec05_300k-mp4

MITOCW MITRES2_002S10nonlinear_lec05_300k-mp4 MITOCW MITRES2_002S10nonlinear_lec05_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 information

Constitutive Relations

Constitutive Relations Constitutive Relations Dr. Andri Andriyana Centre de Mise en Forme des Matériaux, CEMEF UMR CNRS 7635 École des Mines de Paris, 06904 Sophia Antipolis, France Spring, 2008 Outline Outline 1 Review of field

More information

Using MATLAB and. Abaqus. Finite Element Analysis. Introduction to. Amar Khennane. Taylor & Francis Croup. Taylor & Francis Croup,

Using 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 information

Stress analysis of a stepped bar

Stress 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 information

Mechanics of Earthquakes and Faulting

Mechanics of Earthquakes and Faulting Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Overview Milestones in continuum mechanics Concepts of modulus and stiffness. Stress-strain relations Elasticity Surface and body

More information

A new closed-form model for isotropic elastic sphere including new solutions for the free vibrations problem

A new closed-form model for isotropic elastic sphere including new solutions for the free vibrations problem A new closed-form model for isotropic elastic sphere including new solutions for the free vibrations problem E Hanukah Faculty of Mechanical Engineering, Technion Israel Institute of Technology, Haifa

More information

MECH 5312 Solid Mechanics II. Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso

MECH 5312 Solid Mechanics II. Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso MECH 5312 Solid Mechanics II Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso Table of Contents Thermodynamics Derivation Hooke s Law: Anisotropic Elasticity

More information

FVM for Fluid-Structure Interaction with Large Structural Displacements

FVM for Fluid-Structure Interaction with Large Structural Displacements FVM for Fluid-Structure Interaction with Large Structural Displacements Željko Tuković and Hrvoje Jasak Zeljko.Tukovic@fsb.hr, h.jasak@wikki.co.uk Faculty of Mechanical Engineering and Naval Architecture

More information

Simulation in Computer Graphics Elastic Solids. Matthias Teschner

Simulation in Computer Graphics Elastic Solids. Matthias Teschner Simulation in Computer Graphics Elastic Solids Matthias Teschner Outline Introduction Elastic forces Miscellaneous Collision handling Visualization University of Freiburg Computer Science Department 2

More information

Rock Rheology GEOL 5700 Physics and Chemistry of the Solid Earth

Rock Rheology GEOL 5700 Physics and Chemistry of the Solid Earth Rock Rheology GEOL 5700 Physics and Chemistry of the Solid Earth References: Turcotte and Schubert, Geodynamics, Sections 2.1,-2.4, 2.7, 3.1-3.8, 6.1, 6.2, 6.8, 7.1-7.4. Jaeger and Cook, Fundamentals of

More information

Lecture 15 Strain and stress in beams

Lecture 15 Strain and stress in beams Spring, 2019 ME 323 Mechanics of Materials Lecture 15 Strain and stress in beams Reading assignment: 6.1 6.2 News: Instructor: Prof. Marcial Gonzalez Last modified: 1/6/19 9:42:38 PM Beam theory (@ ME

More information

1 Hooke s law, stiffness, and compliance

1 Hooke s law, stiffness, and compliance Non-quilibrium Continuum Physics TA session #5 TA: Yohai Bar Sinai 3.04.206 Linear elasticity I This TA session is the first of three at least, maybe more) in which we ll dive deep deep into linear elasticity

More information

in this web service Cambridge University Press

in this web service Cambridge University Press CONTINUUM MECHANICS This is a modern textbook for courses in continuum mechanics. It provides both the theoretical framework and the numerical methods required to model the behavior of continuous materials.

More information

Glossary. Glossary of Symbols. Glossary of Roman Symbols Glossary of Greek Symbols. Contents:

Glossary. Glossary of Symbols. Glossary of Roman Symbols Glossary of Greek Symbols. Contents: Glossary Glossary of Symbols Contents: Glossary of Roman Symbols Glossary of Greek Symbols Glossary G-l Glossary of Roman Symbols The Euclidean norm or "two-norm." For a vector a The Mooney-Rivlin material

More information

Transactions on Engineering Sciences vol 6, 1994 WIT Press, ISSN

Transactions on Engineering Sciences vol 6, 1994 WIT Press,   ISSN Large strain FE-analyses of localized failure in snow C.H. Liu, G. Meschke, H.A. Mang Institute for Strength of Materials, Technical University of Vienna, A-1040 Karlsplatz 13/202, Vienna, Austria ABSTRACT

More information

Exercise: concepts from chapter 8

Exercise: concepts from chapter 8 Reading: Fundamentals of Structural Geology, Ch 8 1) The following exercises explore elementary concepts associated with a linear elastic material that is isotropic and homogeneous with respect to elastic

More information

Lecture 8. Stress Strain in Multi-dimension

Lecture 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 information

Mechanical Properties of Materials

Mechanical Properties of Materials Mechanical Properties of Materials Strains Material Model Stresses Learning objectives Understand the qualitative and quantitative description of mechanical properties of materials. Learn the logic of

More information

CRITERIA FOR SELECTION OF FEM MODELS.

CRITERIA 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 information

Lecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity

Lecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity Lecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity by Borja Erice and Dirk Mohr ETH Zurich, Department of Mechanical and Process Engineering, Chair of Computational Modeling

More information

MMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM

MMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM B Course Content: A INTRODUCTION AND OVERVIEW Numerical method and Computer-Aided Engineering; Phsical problems; Mathematical models; Finite element method;. B Elements and nodes, natural coordinates,

More information

Finite element analysis of FRP pipelines dynamic stability

Finite element analysis of FRP pipelines dynamic stability Boundary Elements and Other esh Reduction ethods XXXVIII 39 Finite element analysis of FRP pipelines dynamic stability D. G. Pavlou & A. I. Nergaard Department of echanical and Structural Engineering and

More information

On Verification and Validation of Spring Fabric Model

On Verification and Validation of Spring Fabric Model On Verification and Validation of Spring Fabric Model Zheng Gao, Qiangqiang Shi, Yiyang Yang, Bernard Moore, and Xiaolin Li Department of Applied Mathematics and Statistics Stony Brook University Stony

More information

Static Equilibrium; Elasticity & Fracture

Static Equilibrium; Elasticity & Fracture Static Equilibrium; Elasticity & Fracture The Conditions for Equilibrium Statics is concerned with the calculation of the forces acting on and within structures that are in equilibrium. An object with

More information

An Assessment of the LS-DYNA Hourglass Formulations via the 3D Patch Test

An Assessment of the LS-DYNA Hourglass Formulations via the 3D Patch Test 5 th European LS-DYNA Users Conference Code Developments An Assessment of the LS-DYNA Hourglass Formulations via the 3D Patch Test Authors Leonard E. Schwer, Schwer, Engineering & Consulting Services Samuel

More information

Two problems in finite elasticity

Two problems in finite elasticity University of Wollongong Research Online University of Wollongong Thesis Collection 1954-2016 University of Wollongong Thesis Collections 2009 Two problems in finite elasticity Himanshuki Nilmini Padukka

More information

Vaango Theory Guide. Version Version October 1, The Utah Uintah team. and. Biswajit Banerjee

Vaango Theory Guide. Version Version October 1, The Utah Uintah team. and. Biswajit Banerjee Vaango Theory Guide Version Version 17.10 October 1, 2017 The Utah Uintah team and Biswajit Banerjee Copyright 2015-2017 Parresia Research Limited The contents of this manual can and will change significantly

More information

Discrete Analysis for Plate Bending Problems by Using Hybrid-type Penalty Method

Discrete 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 information

Mechanics of Biomaterials

Mechanics of Biomaterials Mechanics of Biomaterials Lecture 7 Presented by Andrian Sue AMME498/998 Semester, 206 The University of Sydney Slide Mechanics Models The University of Sydney Slide 2 Last Week Using motion to find forces

More information

COMPUTATIONAL MODELING OF SHAPE MEMORY MATERIALS

COMPUTATIONAL MODELING OF SHAPE MEMORY MATERIALS COMPUTATIONAL MODELING OF SHAPE MEMORY MATERIALS Jan Valdman Institute of Information Theory and Automation, Czech Academy of Sciences (Prague) based on joint works with Martin Kružík and Miroslav Frost

More information

MODELING OF CONCRETE MATERIALS AND STRUCTURES. Kaspar Willam. Isotropic Elastic Models: Invariant vs Principal Formulations

MODELING OF CONCRETE MATERIALS AND STRUCTURES. Kaspar Willam. Isotropic Elastic Models: Invariant vs Principal Formulations MODELING OF CONCRETE MATERIALS AND STRUCTURES Kaspar Willam University of Colorado at Boulder Class Meeting #2: Nonlinear Elastic Models Isotropic Elastic Models: Invariant vs Principal Formulations Elastic

More information