Hervé DELINGETTE. INRIA Sophia-Antipolis
|
|
- Arnold Bell
- 5 years ago
- Views:
Transcription
1 Hervé DELINGETTE INRIA Sophia-Antipolis French National Research Institute in Computer Science and Automatic Control Asclepios : 3D Segmentation, Simulation Platform, Soft Tissue Modeling CardioSense3D : Cardiac Simulation, SOFA : Open Software Plateform for Medical Simulation
2 Towards the virtual physiological human Physically-based simulation of biological tissues
3 Soft Tissue Characterization Biomechanical behavior of biological tissue is very complex Most biological tissue is composed of several components : Fluids : water or blood Fibrous materials : muscle fiber, neuronal fibers, Membranes : interstitial tissue, Glisson capsule Parenchyma : liver or brain
4 Soft Tissue Characterization To characterize a tissue, its stressstrain relationship is studied Height h Radius r Stress σ = Strain ε = F πr h*-h h Height h* Force F Radius r* Rest Position Cylindrical Piece of Tissue Deformed Position
5 Linear Elastic Material Simplest Material behaviour Only valid for small deformations (less than 5%) 800 Stress Strain
6 Biological Tissue Many complex phenomena arises σ σ V 0 V 1 σ loading V V3 Hysteresis Unloading ε Visco-elasticity ε Linear Domain Non-Linearity Slope =Young Modulus ε Anisotropy
7 Continuum Mechanics Continuum Mechanics Fluid Mechanics Structural Mechanics Linear Elasticity (Anisotropic, heterogeneous) Hyperelasticity Elasticity Plasticity
8 Basics of Continuum Mechanics Deformation Function X Ωa ( X ) R Displacement Function U φ ( X ) = φ( X ) X 3 Rest Position Ω Φ Deformed Position X U(X)
9 Basics of Continuum Mechanics The local deformation is captured by the deformation gradient : = = X X X X X X X X X X F j i ij φ φ φ φ φ φ φ φ φ φ X F = φ Ω X φ(x) Rest Position Deformed Position F(X) is the local affine transformation that maps the neighborhood of X into the neighborhood of φ(x)
10 Basics of Continuum Mechanics Distance between point may not be preserved Rest Position X+dX X Ω Deformed Position φ(x) φ(x+dx) Distance between deformed points ( ) ( ) ( ) T ds X dx X dx ( T = φ + φ φ φ)dx Right Cauchy-Green Deformation tensor C T = φ φ Measures the change of metric in the deformed body
11 Basics of Continuum Mechanics ( ) T Example : Rigid Body φ X motion = RX entails + no deformation T F X ) = φ ( X ) = R ( C = R R = Id Strain tensor captures the amount of deformation It is defined as the distance between C and the Identity matrix E = 1 1 ( T φ φ Id ) = ( C Id )
12 Strain Tensor Diagonal Terms : ε i Capture the length variation along the 3 axis Off-Diagonal Terms :γ i Capture the shear effect along the 3 axis E = 1 1 ε γ γ x xy xz 1 1 γ ε γ y xy yz 1 1 γ γ ε z xz yz
13 Linearized Strain Tensor Use displacement φ ( X ) = rather Id + than U ( Xdeformation ) 1 T T E = U + U + U U ( ) Assume small displacements 1 T E = U + U Lin ( )
14 Hyperelastic Energy The energy required to deform a body is a function of the invariants of strain tensor E : Trace E = I 1 Rest Position Ω Deformed Position Trace E*E= I Determinant of E = I 3 ( φ ) (,, ) 1 W 3 = w I I I dx Ω Total Elastic Energy
15 Linear Elasticity Isotropic Energy ( λ, µ ) : Lamé coefficients λ w X ) = Lin µ ( ) tr E + tr E ( Lin Hooke s Law w (X ) : density of elastic energy Advantage : Quadratic function of displacement λ µ w = ( div U ) + µ U rot U Drawback : Not invariant with respect to global rotation Extension for anisotropic materials
16 Shortcomings of linear elasticity Non valid for «large rotations and displacements»
17 St-Venant Kirchoff Elasticity Isotropic Energy ( λ, µ ) : Lamé coefficients Advantage : Generalize linear elasticity Invariant to global rotations Drawback : Poor behavior in compression Quartic function of displacement Extension for anisotropic materials λ w ( X ) = µ E ( ) tr E + tr
18 St Venant Kirchoff vs Linear Elasticity Linear St Venant Kirchoff Rest Position Linear St Venant Kirchoff
19 Other Hyperelastic Material w ( X ) = µ tre + f I Neo-Hookean Model w( X ) = µ e tre + f ( I 3 ) Fung Isotropic Model µ tre k k w X ) = e + e k Fung Anisotropic Model Veronda-Westman Mooney-Rivlin : ( ) ( ( I 1) ( ) + f ( I ) 1 1 ( γ tre ) e + c tre + f ( ) w( X ) = c I w + 3 ( ) ( X ) = c10tre + c01tre f I3 3
20 Estimating material parameters Complex for biological tissue : Heterogeneous and anisotropic materials Tissue behavior changes between in-vivo and in-vitro Ethics clearance for performing experimental studies Effect of preconditioning Potential large variability across population
21 Soft Tissue Characterization Different possible methods In vitro rheology In vivo rheology Elastometry Solving Inverse problems
22 Soft Tissue Characterization In vitro rheology can be performed in a laboratory. Technique is mature Not realistic for soft tissue (perfusion, )
23 Soft Tissue Characterization In vivo rheology can provide stress/strain relationships at several locations Influence of boundary conditions not well understood Source : Cimit, Boston USA
24 Soft Tissue Characterization Elastometry (MR, Ultrasound) mesure property inside any organ non invasively validation? Only for linear elastic materials Source Echosens, Paris
25 Soft Tissue Characterization Inverse Problems well-suited for surgery simulation (computational approach) require the geometry before and after deformation
26 Soft Tissue Characterization Still difficult to find reliable soft tissue material parameters Example : Liver soft tissue characterization
27 Discretisation techniques Four main approaches : Volumetric Mesh Based Surface Mesh Based Meshless Particles
28 Different types of meshes Surface Elements : Triangle 3, 1 nodes and more Quad 4, 8 nodes and more Volume Elements Tetrahedron 4, 10 nodes Prismatic 6, 15 nodes and more Hexahedron 8, 0 nodes and more
29 Structured vs Unstructured meshes Example 1 : Liver meshed with hexahedra 3 months work (courtesy of ESI) Example : Liver meshed with tetrahedra Automatically generated (1s)
30 Volumetric Mesh Discretization Classical Approaches : Finite Element Method (weak form) Rayleigh Ritz Method (variational form) Finite Volume Method (conservation eq.) Finite Differences Method (strong form) FEM, RRM, FVM are equivalent when using linear elements
31 Rayleigh-Ritz Method Step1 : Choose Finite Element (e.g. linear tetrahedron) Mesh discrediting the domain of computation Hyperelastic Material with its parameters Boundary Conditions ( tr E ) + tr λ w ( X ) = µ E Tetrahedron 4 nodes Fixed nodes Young Modulus Poisson Coefficient
32 Rayleigh-Ritz Method Step Write the elastic energy required to deform a single P 1 element Q 1 u ( P ) = Q P = U i i i i P 4 K P 3 P = Q 4 Q 3 [ ] T K jk U k t i WT i U j jk 1 = µ Q T T ( λ M M + µ M M + ( M M )[ Id ]) Ti [ jk ] i k j i j k i j k 3x3 36.V( Ti ) u( X ) = tre 4 i= 1 λ ( X ) u( i P i M i λi ( X ) = 6V ( T ) = M i Ui 6V ( T ) i )
33 W Rayleigh-Ritz Method Step3 Sum to get the total elastic energy ( ) ( ) T U w I, I I dx = W = U K U = 1, 3 Ω h Write the conservation of energy W ( U ) F T U + ρ ( X ) ( X g )dx Internal Energy = Ω Nodal Forces T i T Gravity Potential Energy i
34 Rayleigh-Ritz Method Step3 Linear Elasticity Write first variation of the energy : KU = R Static case M U&& + C U& + KU = R (t ) Dynamic case HyperElasticity=NonLinear Elasticity K ( U ) = R M U&& + C U& + K = ( U ) R (t ) Static case Dynamic case
35 Surface-Based Methods Only consider the mesh surface under some hypothesis : Linear Elastic Material (sometimes homogeneous) Only interact with organ surface Pros : No need to produce volumetric meshes Much faster than volumetric computation Cons : Only linear material No cutting
36 Surface-Based Methods Possible approaches : Boundary Element Models (BEM) Based on the Green Function of the linear elastic operator Requires homogeneous material Matrix Condensation Full Matrix inversion Iterative Precomputed Generation Solve 3*Ns equations F=KU
37 Other Methods Meshless Methods Use only points inside and specific shape functions Can better optimize location of DOFs Can cope with large deformations Deformation accuracy unknown Particles Smooth Particles Hydrodynamics that interact based on a state equation
38 Dynamic evolution and numerical integration Dynamic evolution Discrete models = lumped mass particles submitted to forces Newtonian evolution (1 st order differential system): { δp= V.dt δv= M -1 F(P,V).dt Explicit schemes: { Euler: δp= V t.dt δv= M -1 F(P t, V t ).dt Runge-Kutta: several evaluations to better extrapolate the new state [press9] Semi-Implicit schemes: Euler: δp= V t+dt.dt Verlet [teschner04] Unstable for large time-step!! { { δv= M -1 F(P t, V t ).dt P t+dt = P t P t-dt.+m -1 F(P t, V t ).dt V t+dt = (P t+dt P t.)dt-1
39 Evolution Implicit schemes [terzopoulos87], [baraff98], [desbrun99], [volino01], [hauth01] First-order expansion of the force: F(P t+dt, V t+dt ) F(P t, V t ) + F/ PδP + F/ VδV Euler implicit { δp= V t+dt.dt H = I - M -1 F/ V dt - M -1 F/ P dt with δv= H -1 Y Y = M -1 F(P t, V t ) + M -1 F/ P V t dt Backward differential formulas (BDF) : Use of previous states Unconditionally stable for any time-step But requires the inversion of a large sparse system Choleski decomposition + relaxation Conjugate gradient Speed and accuracy can be improve through preconditioning (alteration of H)
40 Towards Realistic Interactive Simulation Surgery Simulation must cope with several difficult technical issues : Soft Tissue Deformation Collision Detection Collision Response Haptics Rendering Real-time Constraints : 5Hz for visual rendering Hz for haptic rendering
41 Example of Soft Tissue Models Pre-computed Elastic Model Tensor-Mass and Relaxation-based Model Non-Linear Tensor-Mass Model Computational Efficiency Cutting Simulation Large Displacements - - +
42 Precomputed linear elastic model 9517 Tetrahedra
43 Tensor-Mass Models (low resolution) N = 1394 (634 Tétraèdres)
44 Non-linear Tensor-Mass Models
45 Simulation of surgical gestures Gliding Gripping Cutting (pliers) Cutting (US)
46 Hepatic Surgery Simulation
47 Cardiac Simulation Volumetric Conditions: 4 Cardiac Phases: Filling Isovolumetric Contraction Ejection Isovolumetric Relaxation Pressure Field in the endocardium Isovolumetric Constraint of myocardium Slowed 6 times Color: Action potential 47
48 More Information in Asclepios : 3D Segmentation, Simulation Platform, Soft Tissue Modeling CardioSense3D : Cardiac Simulation, SOFA : Open Software Plateform for Medical Simulation
Surgery Simulation. Olivier Clatz. Hervé Delingette. Asclepios Research project INRIA Sophia Antipolis, France
1 Surgery Simulation Olivier Clatz Hervé Delingette Asclepios Research project INRIA Sophia Antipolis, France Motivations of surgery simulation 2 Increasing complexity of therapy and especially surgery
More informationMultiplicative Jacobian Energy Decomposition Method for Fast Porous Visco-Hyperelastic Soft Tissue Model
Multiplicative Jacobian Energy Decomposition Method for Fast Porous Visco-Hyperelastic Soft Tissue Model Stéphanie Marchesseau 1, Tobias Heimann 1, Simon Chatelin 2, Rémy Willinger 2, and Hervé Delingette
More informationConstitutive 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 informationInverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros
Inverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros Computational Design Forward design: direct manipulation of design parameters Level of abstraction Exploration
More informationContinuum Mechanics and the Finite Element Method
Continuum Mechanics and the Finite Element Method 1 Assignment 2 Due on March 2 nd @ midnight 2 Suppose you want to simulate this The familiar mass-spring system l 0 l y i X y i x Spring length before/after
More 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 informationMeasurement 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 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 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 information2.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 informationDEFORMATION AND FRACTURE ANALYSIS OF ELASTIC SOLIDS BASED ON A PARTICLE METHOD
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm DEFORMATION AND FRACTURE ANALYSIS OF ELASTIC SOLIDS BASED ON A PARTICLE METHOD R. A. Amaro
More informationNDT&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 informationChapter 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 informationModule 4 : Nonlinear elasticity Lecture 25 : Inflation of a baloon. The Lecture Contains. Inflation of a baloon
Lecture 25 : Inflation of a baloon The Lecture Contains Inflation of a baloon 1. Topics in finite elasticity: Hyperelasticity of rubber, elastomers, and biological tissues with examples, M. F Beatty, App.
More informationPerformance 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 informationElements 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 informationContinuum Mechanics. Continuum Mechanics and Constitutive Equations
Continuum Mechanics Continuum Mechanics and Constitutive Equations Continuum mechanics pertains to the description of mechanical behavior of materials under the assumption that the material is a uniform
More informationContinuum 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 informationSimulation 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 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 informationBasic 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 informationDiscontinuous Galerkin methods for nonlinear elasticity
Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity
More information3D 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 informationConstitutive 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 informationComparison 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 informationMacroscopic theory Rock as 'elastic continuum'
Elasticity and Seismic Waves Macroscopic theory Rock as 'elastic continuum' Elastic body is deformed in response to stress Two types of deformation: Change in volume and shape Equations of motion Wave
More informationEngineering 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 informationA Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials
Dublin, October 2010 A Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials FracMan Technology Group Dr Mark Cottrell Presentation Outline Some Physical
More informationBiomechanics. Soft Tissue Biomechanics
Biomechanics cross-bridges 3-D myocardium ventricles circulation Image Research Machines plc R* off k n k b Ca 2+ 0 R off Ca 2+ * k on R* on g f Ca 2+ R0 on Ca 2+ g Ca 2+ A* 1 A0 1 Ca 2+ Myofilament kinetic
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationFVM 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 informationLectures on. Constitutive Modelling of Arteries. Ray Ogden
Lectures on Constitutive Modelling of Arteries Ray Ogden University of Aberdeen Xi an Jiaotong University April 2011 Overview of the Ingredients of Continuum Mechanics needed in Soft Tissue Biomechanics
More informationEXPERIMENTAL IDENTIFICATION OF HYPERELASTIC MATERIAL PARAMETERS FOR CALCULATIONS BY THE FINITE ELEMENT METHOD
Journal of KONES Powertrain and Transport, Vol. 7, No. EXPERIMENTAL IDENTIFICATION OF HYPERELASTIC MATERIAL PARAMETERS FOR CALCULATIONS BY THE FINITE ELEMENT METHOD Robert Czabanowski Wroclaw University
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 informationClassification 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 informationVariational 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 informationMechanics 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 informationA monolithic fluid structure interaction solver Verification and Validation Application: venous valve motion
1 / 41 A monolithic fluid structure interaction solver Verification and Validation Application: venous valve motion Chen-Yu CHIANG O. Pironneau, T.W.H. Sheu, M. Thiriet Laboratoire Jacques-Louis Lions
More informationTwo 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 informationFETI Methods for the Simulation of Biological Tissues
SpezialForschungsBereich F 32 Karl Franzens Universita t Graz Technische Universita t Graz Medizinische Universita t Graz FETI Methods for the Simulation of Biological Tissues Ch. Augustin O. Steinbach
More informationNon-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 informationReal-time simulation of soft tissue deformation for surgical simulation
Real-time simulation of soft tissue deformation for surgical simulation A thesis submitted in fulfilment of the requirements for the degree of Doctor of Philosophy Jinao Zhang BEng (MechEng) (Hons), RMIT
More information3.2 Hooke s law anisotropic elasticity Robert Hooke ( ) Most general relationship
3.2 Hooke s law anisotropic elasticity Robert Hooke (1635-1703) Most general relationship σ = C ε + C ε + C ε + C γ + C γ + C γ 11 12 yy 13 zz 14 xy 15 xz 16 yz σ = C ε + C ε + C ε + C γ + C γ + C γ yy
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationIntegrating Tensile Parameters in Mass-Spring System for Deformable Object Simulation
RAPPORT RECHERCHE LIRIS 1 Integrating Tensile Parameters in Mass-Spring System for Deformable Object Simulation Vincent Baudet, Michael Beuve, Fabrice Jaillet, Behzad Shariat, and Florence Zara Abstract
More informationComparative Study of Variation of Mooney- Rivlin Hyperelastic Material Models under Uniaxial Tensile Loading
Comparative Study of Variation of Mooney- Rivlin Hyperelastic Material Models under Uniaxial Tensile Loading A. N. Jadhav 1, Dr. S.R. Bahulikar, N.H. Sapate 3 1 M Tech Design Engg, Mechanical Engineering,
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 informationLecture 8: Tissue Mechanics
Computational Biology Group (CoBi), D-BSSE, ETHZ Lecture 8: Tissue Mechanics Prof Dagmar Iber, PhD DPhil MSc Computational Biology 2015/16 7. Mai 2016 2 / 57 Contents 1 Introduction to Elastic Materials
More informationConstitutive model of brain tissue suitable for finite element analysis of surgical procedures
Journal of Biomechanics 32 (1999 531 537 Technical Note Constitutive model of brain tissue suitable for finite element analysis of surgical procedures Karol Miller* Department of Mechanical and Materials
More informationConservation 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 informationLecture 8. Stress Strain in Multi-dimension
Lecture 8. Stress Strain in Multi-dimension Module. General Field Equations General Field Equations [] Equilibrium Equations in Elastic bodies xx x y z yx zx f x 0, etc [2] Kinematics xx u x x,etc. [3]
More informationMechanics 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 informationConstitutive Equations
Constitutive quations David Roylance Department of Materials Science and ngineering Massachusetts Institute of Technology Cambridge, MA 0239 October 4, 2000 Introduction The modules on kinematics (Module
More informationA 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 informationFast computation of soft tissue deformations in real-time simulation of surgery with HEML: a comparison of computation time with conventional FEM
Noname manuscript No. (will be inserted by the editor) Fast computation of soft tissue deformations in real-time simulation of surgery with HEML: a comparison of computation time with conventional FEM
More informationAdvanced Simulation of Sealings CADFEM GmbH Rainer Rauch
Advanced Simulation of Sealings CADFEM GmbH Rainer Rauch Recent developments in ANSYS V12 for the simulation of sealings Element technology Material models Contact Robust design and optimization -1- New
More informationLinear 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 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 informationBy drawing Mohr s circle, the stress transformation in 2-D can be done graphically. + σ x σ y. cos 2θ + τ xy sin 2θ, (1) sin 2θ + τ xy cos 2θ.
Mohr s Circle By drawing Mohr s circle, the stress transformation in -D can be done graphically. σ = σ x + σ y τ = σ x σ y + σ x σ y cos θ + τ xy sin θ, 1 sin θ + τ xy cos θ. Note that the angle of rotation,
More informationABHELSINKI UNIVERSITY OF TECHNOLOGY
ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for the Morley plate element Jarkko Niiranen Department of Structural
More informationElements 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 informationSIMULATION OF MECHANICAL TESTS OF COMPOSITE MATERIAL USING ANISOTROPIC HYPERELASTIC CONSTITUTIVE MODELS
Engineering MECHANICS, Vol. 18, 2011, No. 1, p. 23 32 23 SIMULATION OF MECHANICAL TESTS OF COMPOSITE MATERIAL USING ANISOTROPIC HYPERELASTIC CONSTITUTIVE MODELS Tomáš Lasota*, JiříBurša* This paper deals
More informationConstitutive 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 information202 Index. failure, 26 field equation, 122 force, 1
Index acceleration, 12, 161 admissible function, 155 admissible stress, 32 Airy's stress function, 122, 124 d'alembert's principle, 165, 167, 177 amplitude, 171 analogy, 76 anisotropic material, 20 aperiodic
More informationDeformable Materials 2 Adrien Treuille. source: Müller, Stam, James, Thürey. Real-Time Physics Class Notes.
Deformable Materials 2 Adrien Treuille source: Müller, Stam, James, Thürey. Real-Time Physics Class Notes. Goal Overview Strain (Recap) Stress From Strain to Stress Discretization Simulation Overview Strain
More informationTensor Visualization. CSC 7443: Scientific Information Visualization
Tensor Visualization Tensor data A tensor is a multivariate quantity Scalar is a tensor of rank zero s = s(x,y,z) Vector is a tensor of rank one v = (v x,v y,v z ) For a symmetric tensor of rank 2, its
More informationFull-field measurements and identification for biological soft tissues: application to arteries in vitro
Centre for Health Engineering CNRS UMR 5146 INSERM IFR 143 Prof. Stéphane Avril Full-field measurements and identification for biological soft tissues: application to arteries in vitro using single-gage
More informationA unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation
A unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation S. Bordère a and J.-P. Caltagirone b a. CNRS, Univ. Bordeaux, ICMCB,
More informationLecture 1 NONLINEAR ELASTICITY
Lecture 1 NONLINEAR ELASTICITY Soft-Matter Engineering: Mechanics, Design, & Modeling Mechanics Rigid/Hard Systems bending or torsion (but not stretching) small strains (~0.1% metals; ~1% plastics) linearized
More informationContinuum Mechanics and Theory of Materials
Peter Haupt Continuum Mechanics and Theory of Materials Translated from German by Joan A. Kurth Second Edition With 91 Figures, Springer Contents Introduction 1 1 Kinematics 7 1. 1 Material Bodies / 7
More informationInternational 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 informationAN ANISOTROPIC PSEUDO-ELASTIC MODEL FOR THE MULLINS EFFECT IN ARTERIAL TISSUE
XI International Conference on Computational Plasticity. Fundamentals and Applications COMPLAS XI E. Oñate, D.R.J. Owen, D. Peric and B. Suárez (Eds) AN ANISOTROPIC PSEUDO-ELASTIC MODEL FOR THE MULLINS
More informationProceedings of the ASME th International Conference on Ocean, Offshore and Arctic Engineering OMAE2016 June 19-24, 2016, Busan, South Korea
Proceedings of the ASME 26 35th International Conference on Ocean, Offshore and Arctic Engineering OMAE26 June 9-24, 26, Busan, South Korea OMAE26-54554 LOCAL STRAIN AND STRESS CALCULATION METHODS OF IRREGULAR
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 informationExample 3.7 Consider the undeformed configuration of a solid as shown in Figure 3.60.
162 3. The linear 3-D elasticity mathematical model The 3-D elasticity model is of great importance, since it is our highest order hierarchical model assuming linear elastic behavior. Therefore, it provides
More informationDEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS
DEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS Mohsen Safaei, Wim De Waele Ghent University, Laboratory Soete, Belgium Abstract The present work relates to the
More informationExistence of minimizers for the pure displacement problem in nonlinear elasticity
Existence of minimizers for the pure displacement problem in nonlinear elasticity Cristinel Mardare Université Pierre et Marie Curie - Paris 6, Laboratoire Jacques-Louis Lions, Paris, F-75005 France Abstract
More informationBenchmarkingfiniteelement simulation of rigid indenters in elastomers S.J. Jerrams, N. Reece-Pinchin
Benchmarkingfiniteelement simulation of rigid indenters in elastomers S.J. Jerrams, N. Reece-Pinchin Abstract Verifications of finite element techniques applied to elastomers are difficult to achieve since
More informationSoft Tissue Simulation Based on Measured Data
Soft Tissue Simulation Based on Measured Data M. Hauth 1, J. Gross 1,2, W. Straßer 1, and G.F. Buess 2 1 WSI/GRIS, University of Tübingen 2 MIC, Department of General Surgery, University Hospital Tübingen
More informationMECH 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 informationCAEFEM v9.5 Information
CAEFEM v9.5 Information Concurrent Analysis Corporation, 50 Via Ricardo, Thousand Oaks, CA 91320 USA Tel. (805) 375 1060, Fax (805) 375 1061 email: info@caefem.com or support@caefem.com Web: http://www.caefem.com
More informationNon linear Biomechanical Model of the Liver
Non linear Biomechanical Model of the Liver Stéphanie Marchesseau, Simon Chatelin, Hervé Delingette To cite this version: Stéphanie Marchesseau, Simon Chatelin, Hervé Delingette. Non linear Biomechanical
More informationCRITERIA FOR SELECTION OF FEM MODELS.
CRITERIA FOR SELECTION OF FEM MODELS. Prof. P. C.Vasani,Applied Mechanics Department, L. D. College of Engineering,Ahmedabad- 380015 Ph.(079) 7486320 [R] E-mail:pcv-im@eth.net 1. Criteria for Convergence.
More informationRock 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 informationINVERSE 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 informationNatural States and Symmetry Properties of. Two-Dimensional Ciarlet-Mooney-Rivlin. Nonlinear Constitutive Models
Natural States and Symmetry Properties of Two-Dimensional Ciarlet-Mooney-Rivlin Nonlinear Constitutive Models Alexei Cheviakov, Department of Mathematics and Statistics, Univ. Saskatchewan, Canada Jean-François
More informationYou 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 informationChapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature
Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte
More informationTesting Elastomers and Plastics for Marc Material Models
Testing Elastomers and Plastics for Marc Material Models Presented by: Kurt Miller Axel Products, Inc. axelproducts.com We Measure Structural Properties Stress Strain Time-Temperature Test Combinations
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 informationConstitutive Relations
Constitutive Relations Andri Andriyana, Ph.D. 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
More informationENGN 2340 Final Project Report. Optimization of Mechanical Isotropy of Soft Network Material
ENGN 2340 Final Project Report Optimization of Mechanical Isotropy of Soft Network Material Enrui Zhang 12/15/2017 1. Introduction of the Problem This project deals with the stress-strain response of a
More informationAnalytical 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 informationLinearized theory of elasticity
Linearized theory of elasticity Arie Verhoeven averhoev@win.tue.nl CASA Seminar, May 24, 2006 Seminar: Continuum mechanics 1 Stress and stress principles Bart Nowak March 8 2 Strain and deformation Mark
More informationPhysics of Continuous media
Physics of Continuous media Sourendu Gupta TIFR, Mumbai, India Classical Mechanics 2012 October 26, 2012 Deformations of continuous media If a body is deformed, we say that the point which originally had
More informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More informationElasticity in two dimensions 1
Elasticity in two dimensions 1 Elasticity in two dimensions Chapters 3 and 4 of Mechanics of the Cell, as well as its Appendix D, contain selected results for the elastic behavior of materials in two and
More informationDevelopment and application of time-lapse ultrasonic tomography for laboratory characterisation of localised deformation in hard soils / soft rocks
Development and application of time-lapse ultrasonic tomography for laboratory characterisation of localised deformation in hard soils / soft rocks Erika Tudisco Research Group: Stephen A. Hall Philippe
More informationA Review On Methodology Of Material Characterization And Finite Element Modelling Of Rubber-Like Materials
IOSR Journal of Engineering (IOSRJEN) ISSN (e): 50-301, ISSN (p): 78-8719 PP 06-10 www.iosrjen.org A Review On Methodology Of Material Characterization And Finite Element Modelling Of Rubber-Like Materials
More informationModelling Rubber Bushings Using the Parallel Rheological Framework
Modelling Rubber Bushings Using the Parallel Rheological Framework Javier Rodríguez 1, Francisco Riera 1, and Jon Plaza 2 1 Principia, Spain; 2 Cikatek, Spain Abstract: Bushings are anti vibration components
More information