Surgery Simulation. Olivier Clatz. Hervé Delingette. Asclepios Research project INRIA Sophia Antipolis, France
|
|
- Merilyn Johns
- 5 years ago
- Views:
Transcription
1 1 Surgery Simulation Olivier Clatz Hervé Delingette Asclepios Research project INRIA Sophia Antipolis, France
2 Motivations of surgery simulation 2 Increasing complexity of therapy and especially surgery Increasing need for training surgeons and residents Medical malpractice has become socially and economically unacceptable Increasing need for objective evaluation of surgeons (see Cordis Nitanol endovascular carotid stent) Natural extension of surgery planning
3 3 Need for Training Hand-eye Synchronisation Camera being manipulated by an assistant Long instruments going through a fixed point in the abdomen
4 4 Current Training Techniques Mechanical Simulators Animals Procedure Duration (min) Average Duration of an hysterectomy with laparoscopy Cadavers Patients (source, Department of Obstetrics and Gynecology, Helsinki University Central Hospital) Total number of patients Resident Training (source, Ayudamos) Example of a mechanical simulator (source, US Surgical Corporation) Surgery Training on a pig at EITS
5 Goal for simulation : Training versus Rehearsal 5 Training: Modelling a standard patient for teaching classical or rare situations Rehearsal: Modelling a specific patient to plan and rehearse a delicate intervention, and evaluate consequences beforehand
6 Simulator Workflow 6 Position Collision Contact Deformation Force Force
7 7 Different Technical Issues Mesh Reconstruction from Images Soft Tissue Modeling Tissue Cutting Collision Detection Contact Modeling Surface Rendering Haptic Feedback With realtime constraints
8 8 Different Technical Issues Mesh Reconstruction from Images Soft Tissue Modeling Tissue Cutting Collision Detection Contact Modeling Surface Rendering Haptic Feedback
9 Liver Reconstruction 9 Deformation from a reference model reconstructed from the «Visible Human Project»
10 10 Different Technical Issues Mesh Reconstruction from Images Soft Tissue Modeling Tissue Cutting Collision Detection Contact Modeling Surface Rendering Haptic Feedback
11 11 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
12 Estimating material parameters 12 Complex for biological tissue : Heterogeneous and anisotropic materials Tissue behavior changes between in-vivo and in-vitro Effect of preconditioning Potential large variability across population Ethics clearance for performing experimental studies
13 13 Soft Tissue Characterization Different possible methods In vitro rheology In vivo rheology Elastometry Solving Inverse problems
14 Soft Tissue Characterization 14 In vitro rheology can be performed in a laboratory. Technique is mature Not realistic for soft tissue (perfusion, )
15 Soft Tissue Characterization 15 In vivo rheology can provide stress/strain relationships at several locations Influence of boundary conditions not well understood Source : Cimit, Boston USA
16 Soft Tissue Characterization 16 Elastometry (MR, Ultrasound) mesure property inside any organ non invasively validation? Only for linear elastic materials Source Echosens, Paris
17 Soft Tissue Characterization 17 Inverse Problems well-suited for surgery simulation (computational approach) requires geometry & BC before and after deformation
18 Soft Tissue Characterization 18 To characterize a tissue, its stressstrain relationship is studied eight h Radius r Stress σ = Strain ε = F πr 2 h*-h h Height h* Force F Radius r* Rest Position Cylindrical Piece of Tissue
19 Linear Elastic Material 19 Simplest Material behaviour Only valid for small deformations (less than 5%) Contrainte Déplacement
20 Biological Tissue 20 Far more complex phenomena arises σ σ σ loading V 0 V 1 V 2 V 3 Hysteresis Unloading ε Visco-elasticity ε Linear Domain Non-Linearity Slope =Young Modulus ε Anisotropy
21 Continuum Mechanics 21 Continuum Mechanics Fluid Mechanics Structural Mechanics Linear Elasticity (Anisotropic, heterogeneous) Hyperelasticity Elasticity Plasticity
22 Material Modeling: Basics of Continuum Mechanics 22 Deformation Function Displacement Function Rest Position Φ X U Ωa φ ( X ) R ( X ) = φ( X ) X 3 Ω Deformed Position X U(X)
23 23 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)
24 Basics of Continuum Mechanics 24 Distance between point may not be preserved Rest Position Ω Deformed Position X+dX X φ(x) φ(x+dx) Distance between deformed points ( ) 2 ( ) ( ) T ds X dx X dx ( T = φ + φ φ φ)dx Cauchy-Green Deformation tensor C 2 T = φ φ Measures the change of metric in the deformed body
25 Basics of Continuum Mechanics 25 Example : Rigid Body motion entails no deformation φ ( X ) = RX + T F ( X ) = φ ( X ) = R T Strain tensor captures the amount of deformation It is defined as the distance between C and the Identity matrix E = 1 2 C = R 1 2 R = Id ( T φ φ Id ) = ( C Id )
26 Strain Tensor 26 Diagonal Terms : ε i Capture the length variation along the 3 axis Off-Diagonal Terms :γ i Capture the shear effect along the 3 axis E = ε γ γ x xy xz γ ε γ y xy yz γ γ ε z xz yz
27 Linearized Strain Tensor 27 Use displacement rather than deformation E φ 1 = U 2 ( X ) = Id + U ( X ) ( T T + U + U U ) Assume small displacements E Lin = 2 1 ( T U + U )
28 Hyperelastic Energy 28 The energy required to deform a body is a function of the invariants of strain tensor E : Trace E = = I 1 Trace E*E= I 2 Determinant of E = I 3 Rest Position Ω Deformed Position ( φ ) (,, ) 1 2 W 3 = w I I I dx Ω Total Elastic Energy
29 Isotropic Energy ( λ, µ ) w (X ) : Lamé coefficients Linear Elasticity : density of elastic energy λ w X ) = Lin µ 2 ( ) 2 2 tr E + tr E ( Lin Hooke s Law 29 Advantage : Quadratic function of displacement λ 2 2 µ w = ( div U ) + µ U rot U 2 2 Drawback : Not invariant with respect to global rotation Extension for anisotropic materials 2
30 Shortcomings of linear elasticity Non valid for «large rotations and displacements» 30
31 St-Venant Kirchoff Elasticity 31 Isotropic Energy λ w ( X ) = 2 µ E ( λ, µ ) : Lamé coefficients Advantage : Generalize linear elasticity Invariant to global rotations Drawback : Poor behavior in compression Quartic function of displacement Extension for anisotropic materials ( ) 2 2 tr E + tr
32 St Venant Kirchoff vs Linear Elasticity 32 Linear St Venant Kirchoff Rest Position Linear St Venant Kirchoff
33 Other Hyperelastic Material Neo-Hookean Model Fung Isotropic Model w ( X ) = µ + I tre f ( ) 3 2 = µ ( ) w( X ) e tre + f I Fung Anisotropic Model Veronda-Westman Mooney-Rivlin : µ tre k 4 w( X ) = e + 2 k 1 2 k2 ( I 1) ( e ) + f ( I ) 1 1 ( γ tre ) 2 e + c tre + f ( ) w( X ) = c I w ( ) 2 ( X ) = c10tre + c01tre f I3
34 Discretisation techniques 34 Four main approaches : Volumetric Mesh Based Surface Mesh Based Meshless Particles
35 Different types of meshes Surface Elements : 35 Triangle 3, 12 nodes and more Quad 4, 8 nodes and more Volume Elements Tetrahedron 4, 10 nodes Prismatic 6, 15 nodes and more Hexahedron 8, 20 nodes and more
36 Structured vs Unstructured meshes Example 1 : Liver meshed with hexahedra 36 3 months work (courtesy of ESI) Example 2: Liver meshed with tetrahedra Automatically generated (10s)
37 Volumetric Mesh Discretization 37 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
38 Rayleigh-Ritz Method Step1 : Choose Finite Element (e.g. linear tetrahedron) Mesh discrediting the domain of computation Hyperelastic Material with its parameters Boundary Conditions 38 Tetrahedron 4 nodes Fixed nodes 2 2 ( tr E ) + tr λ w ( X ) = µ E 2 Young Modulus Poisson Coefficient
39 Rayleigh-Ritz Method Step2 Write the elastic energy required to deform a single element 39 P 1 Q 1 u ( P ) = Q P = U i i i i P 4 WT i 1 P 3 P 2 = jk U t j Q 4 Q 3 [ ] T K jk U k i Q 2 T T ( λ M M + µ M M + ( M M )[ Id ]) Ti [ K 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 )
40 Rayleigh-Ritz Method 40 Step3 Sum to get the total elastic energy W ( ) ( ) T U w I, I I dx = W = U K U = 1 2, 3 Ω h Write the conservation of energy T i T i W ( U ) F T U + ρ ( X ) ( X g )dx Internal Energy = Ω Nodal Forces Gravity Potential Energy
41 Step3 Rayleigh-Ritz Method Write first variation of the energy : Linear Elasticity KU = R Static case 41 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
42 Surface-Based Methods 42 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
43 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 Evolution 43 Runge-Kutta: several evaluations to better extrapolate the new state [press92] Unstable for large time-step!! Semi-Implicit schemes: { { Euler: δp= V t+dt.dt δv= M -1 F(P t, V t ).dt P t+dt = 2P t P t-dt.+m -1 F(P t, V t ).dt2 V t+dt = (P t+dt P t.)dt-1
44 Evolution 44 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 2 with δv= H -1 Y Y = M -1 F(P t, V t ) + M -1 F/ P V t dt 2 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)
45 Example of Soft Tissue Models 45 Pre-computed Elastic Model Tensor-Mass and Relaxation-based Model Non-Linear Tensor-Mass Model Computational Efficiency Cutting Simulation Large Displacements - - +
46 Precomputed linear elastic model Tetrahedra
47 47 Different Technical Issues Mesh Reconstruction from Images Soft Tissue Modeling Tissue Cutting Collision Detection Contact Modeling Surface Rendering Haptic Feedback
48 Different algorithms for cutting tetrahedral meshes 48 Split of tetrahedra [Bielser, 2000] [Mohr, 2000] [Nienhuys, 2001] + Accurate, realistic - Decrease of Mesh Quality Removing Tetrahedra [Forest, 2002] + Keeps a good mesh quality - Gross cut
49 49 Proposed Technique Remove Tetrahedra Refine Mesh before removing material
50 50 Dynamic Refinement Example of refine-before-cutting
51 Refinement by Edge Split 51 Decomposition of tetrahedra by edge split
52 Topological Singularities 52 Removing a tetrahedron may create a singularity (zero thickness at edge and vertices) (see [forest]) Singularity
53 Non-linear Tensor-Mass Models 53
54 54 Different Technical Issues Mesh Reconstruction from Images Soft Tissue Modeling Tissue Cutting Collision Detection Contact Modeling Surface Rendering Haptic Feedback
55 55 Previous Work A lot of research on Collision Detection Hierarchy of oriented bounding boxes: Gottshalk & al. - Obb-tree: A hierarchical structure for rapid interference detection - SIGGRAPH 96 public domain package RAPID Very efficient, but needs pre-computation
56 56 The Rendering Process Camera = viewing volume + projection Two steps: geometry & rasterization
57 Collision Detection and Rendering analogy 57 a tool collides the organ a part of the organ is inside the tool if we define a camera with a viewing volume that matches the tool geometry, the organ will be in the picture.
58 58 Different Technical Issues Mesh Reconstruction from Images Soft Tissue Modeling Tissue Cutting Collision Detection Contact Modeling Surface Rendering Haptic Feedback
59 Tool-Soft Tissue Interaction 59 Prevent penetration of tool inside the soft tissue Detect intersections Push explicitly mesh vertices outside the tool
60 First Approach [Picinbono, 2001] 60 2 different tools : tip and handle Compute average normal in the neighborhood of the contact Projection of vertices in this plane Handle Tip
61 Collision Processing 61 Contact with the tip of the instrument Projection on the plane defined by the tip of the instrument and the average normal of intersected triangles
62 Collision Processing 62 Contact with the handle of the tool Projection on the plane defined by the tool direction and the averaged normal direction
63 Collision Processing 63 Perform 2 detections simultaneously
64 Possible interactions 64 Slip on the surface
65 Limitations of this approach 65 Same normal vector for all triangles in the same neighborhood Leads to instabilities when handling a complex geometry
66 66 New approach Three steps Prevent vertices to collide with the tool axis Move vertices near the tip of the tool Move vertices outside the volume of the tool
67 Example 67
68 68 Different Technical Issues Mesh Reconstruction from Images Soft Tissue Modeling Tissue Cutting Collision Detection Contact Modeling Surface Rendering Haptic Feedback
69 69 Haptic Feedback Principle Give a realistic sense of contact with the soft tissue Motivation Increase realism Naturally limit the amplitude of hand motion Pitfalls Frequency update of haptics > 500 Hz Mouvement non contraint par le retour d effort Frequency update of deformable models 30 Hz
70 70 First approach [Picinbono, 2001] Force F Computed Forces Estimated Forces Deformable Model Frequency Update 20 Hz Force Computation Position Force Feedback Frequency Update 300 Hz Force Extrapolation Example of extrapolation t Unstable if complex geometry Difficult extrapolation for large hand motions P' P F p t F ( F ) n ( ) = Fn + n n 1 Pn Pn 1 t n t < t n+1
71 Local Model [Mendoza, 2001] [Balaniuk, 1999] [Mark, 1996] 71 Modèle local Deformable Model Position Retour d effort Update Frequency 20 Hz Computation of a local model Update Frequency 300 Hz Force Computation from a local model Smooth Transition from one local model to the next
72 72 Computing the local model Described as a set of planes One model for the tip One model for the handle
73 73 Force Computation Proportional to the penetration of the tool tip in the planes described by the local model F = k.( EndP O P n r ). P
74 Tensor-Mass Models (low resolution) 74 N = 1394 (6342 Tétraèdres)
75 Simulation of surgical gestures 75 Gliding Gripping Cutting (pliers) Cutting (US)
76 SOFA : 76
77 SOFA : 77
78 Thank you 78
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationAn Overview of Fluid Animation. Christopher Batty March 11, 2014
An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationTraction on the Retina Induced by Saccadic Eye Movements in the Presence of Posterior Vitreous Detachment
Traction on the Retina Induced by Saccadic Eye Movements in the Presence of Posterior Vitreous Detachment Colangeli E., Repetto R., Tatone A. and Testa A. Grenoble, 24 th October 2007 Table of contents
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 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 informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
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 informationExam paper: Biomechanics
Exam paper: Biomechanics Tuesday August 10th 2010; 9.00-12.00 AM Code: 8W020 Biomedical Engineering Department, Eindhoven University of Technology The exam comprises 10 problems. Every problem has a maximum
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 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 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 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 informationNUMERICAL MODELING OF INSTABILITIES IN SAND
NUMERICAL MODELING OF INSTABILITIES IN SAND KIRK ELLISON March 14, 2008 Advisor: Jose Andrade Masters Defense Outline of Presentation Randomized porosity in FEM simulations Liquefaction in FEM simulations
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 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 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 informationInstabilities and Dynamic Rupture in a Frictional Interface
Instabilities and Dynamic Rupture in a Frictional Interface Laurent BAILLET LGIT (Laboratoire de Géophysique Interne et Tectonophysique) Grenoble France laurent.baillet@ujf-grenoble.fr http://www-lgit.obs.ujf-grenoble.fr/users/lbaillet/
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 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 informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
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 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 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 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 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 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 informationNUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH DYNAMIC FRACTURE
NUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH DYNAMIC FRACTURE Kevin G. Wang (1), Patrick Lea (2), and Charbel Farhat (3) (1) Department of Aerospace, California Institute of Technology
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 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 informationA FINITE-VOLUME DISCRETIZATION FOR DEFORMATION OF FRACTURED MEDIA
A FINITE-VOLUME DISCRETIZATION FOR DEFORMATION OF FRACTURED MEDIA Eren Ucar a, Eirik Keilegavlen a, Inga Berre a,b, Jan Martin Nordbotten a,c a Department of Mathematics, University of Bergen, Bergen,
More informationEntropy generation and transport
Chapter 7 Entropy generation and transport 7.1 Convective form of the Gibbs equation In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
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 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 informationA Finite Element Model for Numerical Analysis of Sintering
A Finite Element Model for Numerical Analysis of Sintering DANIELA CÂRSTEA High-School Group of Railways, Craiova ION CÂRSTEA Department of Computer Engineering and Communication University of Craiova
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 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 informationMechanics 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 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 informationNonlinear Equations for Finite-Amplitude Wave Propagation in Fiber-Reinforced Hyperelastic Media
Nonlinear Equations for Finite-Amplitude Wave Propagation in Fiber-Reinforced Hyperelastic Media Alexei F. Cheviakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon, Canada
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationBACKGROUNDS. Two Models of Deformable Body. Distinct Element Method (DEM)
BACKGROUNDS Two Models of Deformable Body continuum rigid-body spring deformation expressed in terms of field variables assembly of rigid-bodies connected by spring Distinct Element Method (DEM) simple
More informationModeling, Simulating and Rendering Fluids. Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan
Modeling, Simulating and Rendering Fluids Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan Applications Mostly Hollywood Shrek Antz Terminator 3 Many others Games Engineering Animating Fluids is
More informationFINAL EXAMINATION. (CE130-2 Mechanics of Materials)
UNIVERSITY OF CLIFORNI, ERKELEY FLL SEMESTER 001 FINL EXMINTION (CE130- Mechanics of Materials) Problem 1: (15 points) pinned -bar structure is shown in Figure 1. There is an external force, W = 5000N,
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 informationReal Time Physics Class Notes. Matthias Müller, NVIDIA Jos Stam, Autodesk Doug James, Cornell University Nils Thürey, ETH Zurich
Real Time Physics Class Notes Matthias Müller, NVIDIA Jos Stam, Autodesk Doug James, Cornell University Nils Thürey, ETH Zurich Contents 1 Introduction 5 1.1 Real-time vs. Off-line Physics........................
More information