Solutions of problems in computational physics. Waldemar Rachowicz Cracow University of Technology ul. Warszawska 24, Cracow, Poland
|
|
- Tyrone Benson
- 5 years ago
- Views:
Transcription
1 Solutions of problems in computational physics Waldemar Rachowicz Cracow University of Technology ul. Warszawska 24, Cracow, Poland
2 Crack problem. Rvatshev s functions and BEM a) b) Figure: Zadanie szczeliny: a) przemieszczenia u z, b) naprȩżenia σ zz 19 March, 2018 Solutions of problems in computational physics 2 / 44
3 Crack problem a) b) Figure: Zadanie szczeliny: a) przemieszczenie u z, b) naprȩżenie σ zz wzd luż osi szczeliny 19 March, 2018 Solutions of problems in computational physics 3 / 44
4 Linear elasticity as a model elliptic bvp Formulation of linear elasticity: strong σ = f in Ω ɛ = 1 2 ( u + T u) σ = 2µɛ + λtrɛi u = û on Γ D σn = ˆt on Γ N 9 8 >= >< >; >: variational Find u V + û : a(u, v) = l(v) v V Z a(u, v) = [µɛ(u) : ɛ(v) + λtrɛ(u)trɛ(v)]dx Z Ω Z l(v) = f vdx + Ω ˆt n ds Γ N V = {v [H 1 (Ω)] 1 : v = 0 on Γ D }. Bilinear form a(u, v) is continuous and V -coercive, linear functional l(v) is continuous: M > 0 : a(u, v) M u 1,Ω v 1,Ω u, v H 1 (Ω) α > 0 : a(v, v) α v 2 1,Ω v V H1 (Ω) c > 0 : l(v) c v 1,Ω v H 1 (Ω) 19 March, 2018 Solutions of problems in computational physics 4 / 44
5 Aspects of adaptive methods: adaptivity Types: h-, p- and hp-adaptivity. Typical feed-back adaptive algorithm: 1. Solve the problem on the current mesh (99% of the cost). 2. Find error indicators η K 3. Stop if ( K η2 K )1/2 T OL 4. Refine elements K for which: η K > α max L η L, 0 < α < Go to 1. Goal oriented adaptivity: η K η u K ηg K 19 March, 2018 Solutions of problems in computational physics 5 / 44
6 Compressible Navier-Stokes equations Principles of conservation of mass, momentum and total energy: ϱ t + m i t [m j ] = 0 (0.1) x j + [ ] mi m j + δ ij p = [τ ij ] x j ϱ x j [ ] mj (e + p) = [ ] m i τ ij x j ϱ x j ϱ q i e t + Constitutive equations: p ϱθ = R, ι = C vϱθ, p = (γ 1)ι, γ = 1.4 τ ij = 2µε ij + λε kk, µ = µ 0 ( θ θ 0 ) 1.5 c 1 + θ 0 c 1 + θ, q i = κθ,i, κ = C p Pr µ, 19 March, 2018 Solutions of problems in computational physics 6 / 44
7 Flow around a plate, M = 3, Re = 1000 a Figure: Flow around a plate, M = 3, Re = (a) Distribution of density, (b) elevation of pressure b 19 March, 2018 Solutions of problems in computational physics 7 / 44
8 Flow around a cylinder, M = 8 a Figure: Flow around a cyliner, M = 8. (a) Distribution of density, (b) elevation of pressure b 19 March, 2018 Solutions of problems in computational physics 8 / 44
9 Flow around a cylinder with an impinging shock, M = 8 a Figure: Flow around a cyliner with an impinging shock, M = 8. (a) Distribution of density, (b) contour map of pressure 19 March, 2018 Solutions of problems in computational physics 9 / 44 b
10 Incompressible Navier-Stokes equations { u t + (u )u + p = ν u, u = 0. a b Figure: Incompressible flow around a NACA-81 profile. (a) Distribution of pressure, (b) elevation of velocity u y 19 March, 2018 Solutions of problems in computational physics 10 / 44
11 Automatic h- or hp-adaptivity, L. Demkowicz, J. Kurtz, D. Pardo, W. Rachowicz Automatic h- or hp-adaptivity: We control two meshes: the coarse mesh h, p the fine mesh h/2, p + 1 Algorithm: 1. Investigate interpolation error of u h/2,p+1 on the elements of the coarse mesh with various trial refinements. 2. Select refinements resulting in the largest reduction of interpolation error per 1 new dof. 3. Perform selected refinements on the coarse mesh. 4. Obtain the fine mesh. 5. Solve on fine mesh (99% of the cost). 6. Go to 1. Goal oriented version of the algorithm exists. 19 March, 2018 Solutions of problems in computational physics 11 / 44
12 Cantilever beam Figure: Cantilever beam. Coarse and fine hp meshes, distribution of σ yz, σ 0. Ω = [0, 1] [0, 4] [0, 1], load: t z = 1 for 3.5 < y < March, 2018 Solutions of problems in computational physics 12 / 44
13 Hyperboloidal shell Figure: Hyperboloidal shell. Coarse and fine hp meshes; σ φφ, σ zz, σ 0. Ω : r2 z2 = 1, z [ 3, 1], t = 0.05; load: t x = 1 19 March, 2018 Solutions of problems in computational physics 13 / 44
14 Streamer - cable for aquisition of seismic data a b Figure: (a) structure of a segment (gel, skin, spacer, rope, hydrophone) (b) h-adaptive mesh 19 March, 2018 Solutions of problems in computational physics 14 / 44
15 Streamer Time-harmonic elasto-dynamics: σ ρω 2 u = 0 in Ω u n n = iωρ wc w u n on Γ N u = û on Γ D hydrophone E = 4.0 GPa, ν = 0.4, ρ = 1500 kg/m 3 gel E =.24 GPa, ν = 0.43, ρ = 1040 kg/m 3 skin E = 0.02 GPa, ν = 0.45, ρ = 1200 kg/m 3, spacer E = 1.8 GPa, ν = 0.3, ρ = 1200 kg/m 3, rope E = 41 GPa, ν = 0.3, ρ = 1400 kg/m 3 19 March, 2018 Solutions of problems in computational physics 15 / 44
16 Streamer a b c Figure: Distribution of pressure [db] (a) over a plane of symmetry, (b) over a hydrophone, (c) along a center line (h and hp simulations) 19 March, 2018 Solutions of problems in computational physics 16 / 44
17 Goal-oriented adaptivity Hyperboloidal shell: r 2 a 2 z2 c 2 = 1, a in = 1, a out = 1.05, c in = c out = 2, z [ 3, 1] Quantities of interest: F (u) = wσ ij ζ k dv, D( S) i, j - shell local directions S - small part of middle surface, meas(w) = 1 i, j = 1, k = 1: F = (m 11 ) av - bending moment i, j = 1, k = 0: F = (n 11 ) av - membrane force Tangential orders p 4 to avoid locking! 19 March, 2018 Solutions of problems in computational physics 17 / 44
18 Hyperboloidal shell, goal-oriented adaptivity a b c d e f Figure: (a,b,c) coarse and fine mesh, u r of Green s function for m 11 ; (d,e,f) for n 11. Results: (m 11 ) av = ± 0.7%, (n 11 ) av = 5.03 ± 0.06% 19 March, 2018 Solutions of problems in computational physics 18 / 44
19 Cantilever beam - goal-oriented adaptivity a b c d Figure: (a,b) an h- and hp-adaptive meshes, (c,d) u x and σ xx for Green s function. Q.o.i. F = ω σ 11dx, ω = [0.9375, 1.0] 3. Result: F = ± 0.03% (h-adapt.), F = ± 0.02% (hp-adapt.) 19 March, 2018 Solutions of problems in computational physics 19 / 44
20 BVP in electromagnetics Time-harmonic Maxwell equations: E = jωµh, H = 1 jωµ E, H = jωɛe + σe + J imp, ˆɛ = ɛ jσ/ω Boundary-value problem: 1 / µ E ω2ˆɛe = jωj imp in Ω, F, n E = M imp S on Γ D, n E = jωj imp S on Γ N, Variational formulation: Find E such that n E = M imp S on Γ D and: Z Z 1 Ω µ E F dx ω 2ˆɛE Z Z F dx = jω J imp F ds + jω Ω Ω {z } a(e,f ) Ω J imp S Γ N F ds, {z } l(f ) F : F t = 0, on Γ D. 19 March, 2018 Solutions of problems in computational physics 20 / 44
21 Time-harmonic Maxwell equations, FEM a b Figure: Rozpraszanie fal p laskich na wnȩce, λ = 0.5. Pole elektryczne E y : a) czȩść rzeczywista, b) czȩść urojona 19 March, 2018 Solutions of problems in computational physics 21 / 44
22 Time-harmonic Maxwell equations, FEM Figure: Rozpraszanie fal p laskich na wnȩce, λ = 0.5. Przekroj na rozpraszanie radarowe (RCS) 19 March, 2018 Solutions of problems in computational physics 22 / 44
23 Time-harmonic Maxwell equations, BEM + FEM We solve the Stratton-Chu integral equation: Z Z Z «E(y) = y [γ + D k (E) G ds γ + N (E) G ds + 1 y Γ [γ + N Γ Γ k (E)]G ds. Γ γ + D E := E n, γ+ N := k 1 ( E) n, G(x y) = e jkr, r = x y. 4πr a b Figure: Rozpraszanie na kuli przewodz acej. Sk ladowa Re(E y ) rozwi azania: a) czȩść rzeczywista b) czȩść urojona 19 March, 2018 Solutions of problems in computational physics 23 / 44
24 Time-harmonic linear acoustics, BEM + FEM Burton-Miller formulation, α [0, 1], α = (1 α)j/k: Z» α α p p(y) + (y) = α p G p «G + α 2 G p p «G ds, 2 2 n y Γ n x n x n y n x n x n y G(x y) = ejkr, r = x y. 4πr a b Figure: Rozpraszanie na walcowej pow loce sprȩżystej - ciśnienie: a) czȩść rzeczywista, b) czȩść urojona 19 March, 2018 Solutions of problems in computational physics 24 / 44
25 Inverse medium scattering problems in electromagnetics Figure: Transmitters and receivers surrounding the scatterer 19 March, 2018 Solutions of problems in computational physics 25 / 44
26 Inverse medium scattering in electromagnetics Inverse medium scattering: Consider incident waves E in m, m = 1,..., M For each of waves E in m we measure the scattered waves E mn at N observation points, x n, n = 1,..., N. Based on the measurements E mn we wish to evaluate the distribution of ˆɛ (x) by matching the measurements with the numerical simulations for trial ˆɛ(x) { } ˆɛ (x) = arg min E mn (ˆɛ) E mn 2 + M(ˆɛ) m,n where M is a (possible) stabilizing functional. 19 March, 2018 Solutions of problems in computational physics 26 / 44
27 Ex. kidney-like shape, complex, complete data a) b) Figure: Reconstruction of complex-valued ˆɛ(x) with complete measurement data, kidney-like scatterer with a tumor: a) exact Re(ˆɛ), b) reconstructed Re(ˆɛ), data E mn perturbed by random noise 19 March, 2018 Solutions of problems in computational physics 27 / 44
28 Ex. kidney-like shape, complex, complete data c) d) Figure: Reconstruction of complex-valued ˆɛ(x) with complete measurement data, kidney-like scatterer with a tumor: c) exact Im(ˆɛ), d) reconstructed Im(ˆɛ); data E mn perturbed by random noise 19 March, 2018 Solutions of problems in computational physics 28 / 44
29 Ex. scatterer in two half-spaces, complex-valued ˆɛ Figure: Reconstruction of complex-valued ˆɛ(x) with incomplete measurement data, kidney-like scatterer with a tumor in two subspaces: a)meshes, b) exact Re(ˆɛ), c) reconstructed Re(ˆɛ), data E mn perturbed by random noise 19 March, 2018 Solutions of problems in computational physics 29 / 44 a) b) c)
30 Ex. scatterer in two half-spaces, complex-valued ˆɛ a) b) Figure: Reconstruction of complex-valued ˆɛ(x) with incomplete measurement data, kidney-like scatterer with a tumor in two subspaces: a) exact Im(ˆɛ), b) reconstructed Im(ˆɛ), E mn perturbed by random noise 19 March, 2018 Solutions of problems in computational physics 30 / 44
31 Limb-like scatterer in a cylinder, complex-valued ˆɛ a) b) c) Figure: Reconstruction of complex-valued ˆɛ(x) with incomplete measurement data, limb-like scatterer with tumors: a) meshes, b) exact Re(ˆɛ), c) reconstructed Re(ˆɛ), data E mn perturbed by random noise 19 March, 2018 Solutions of problems in computational physics 31 / 44
32 Limb-like scatterer in a cylinder, complex-valued ˆɛ a) b) Figure: Reconstruction of complex-valued ˆɛ(x) with incomplete measurement data, limb-like scatterer with tumors: a) exact Im(ˆɛ), b) reconstructed Im(ˆɛ), data E mn perturbed by random noise 19 March, 2018 Solutions of problems in computational physics 32 / 44
33 Ex. 2 PEC spheres, r = 0.555λ, Linear Sampling Method (Peter Monk) Figure: Inverse scattering recovery: two PEC spheres 19 March, 2018 Solutions of problems in computational physics 33 / 44
34 Example: two PEC hexahedra, a = 1.12λ, Linear Sampling Method Figure: Inverse scattering recovery: two PEC hexahedra 19 March, 2018 Solutions of problems in computational physics 34 / 44
35 Finite anisotropic elasticity + contact + plasticity. The notions of finite elasticity: x(x, t) - Lagrangian description F = X x - deformation gradient C = F T F - right Cauchy Green tensor J = det(c) 1/2 = det(f ) > 0 - volume ratio Multiplicative decomposition of deformation gradient: F = J 1/3 F, C = J 2/3 C Fields of fibre directions: M i (X), i = 1, 2, M i = 1. Structural tensors: A i = M i M i, i = 1, 2 Invariants: I 1 = tr(c), I 2 = 1 2 [(tr(c)2 tr(c 2 ), I 4i = C : A i, Ī 1 = J 2/3 I 1, Ī 2 = J 1/3 I 2, Ī 4i = J 2/3 I 4i. 19 March, 2018 Solutions of problems in computational physics 35 / 44
36 Formulation of finite elasticity bvp Split of strain energy function (R. Flory): 2 Ψ = Ψ vol (J) + Ψ iso (Ī 1, Ī 2 ) + W fi (Ī 4i ) i=1 The second Piola-Kirchhoff stress tensor: S = 2 Ψ C = pjc 1 + 2J 2/3 {Ψ iso,1 Dev[I] + Ψ iso,2 Dev[ C 2 ] + 2 i=1 W fi,4dev[a i ]} Strain energies: Ψ vol = K 4 (J 1)2, Ψ iso = µ 2 (Ī 1 3), Ψ fi = { k 1i 2k 2i exp [k2i (Ī 4i 1) 2 ] 1 }. G. Holzapfel, T.Gasser, R. Ogden HGO 19 March, 2018 Solutions of problems in computational physics 36 / 44
37 Figure: Contact of two elastic cylinders solved on meshes of orders p = 1, p = 2 and p = 4 (in subsequent rows). The columns show: a) undeformed meshes, b) deformed meshes, c) the mean pressure, and d) effective stress σ 0 19 March, 2018 Solutions of problems in computational physics 37 / 44 p=1 p=2 p=4 a b c d
38 a b c d e f Figure: Half-cylindrical die intruding into elastic block and a thin-walled tube pressed by a flexible plate: a,d) an h-adaptive mesh on the reference configuration, b,e) displacements u z. c,f) effective stress σ 0 19 March, 2018 Solutions of problems in computational physics 38 / 44
39 The idea of angioplasty a b c Figure: Angioplasty: a) artery, stent and catheter, and displacements u r on (b) the reference and (c) the deformed configuration 19 March, 2018 Solutions of problems in computational physics 39 / 44
40 Angioplasty a Figure: Angioplasty. Distribution of: (a) pressure, and (b) stresses σ rφ in the vicinity of the contact of the tube b 19 March, 2018 Solutions of problems in computational physics 40 / 44
41 Angioplasty r r Location of artery surface Location of artery surface 2 2 1,8 1,8 1,6 1,6 1,4 1,4 1,2 1 0,8 before loading after loading after unloading 1,2 1 0,8 before loading after loading after unloading 0,6 0,6 0,4 0, a z 0,4 0, b z Figure: Location of internal surface of the tube at different stages of loading process for tubes types (a) and (b) 19 March, 2018 Solutions of problems in computational physics 41 / 44
42 Angioplasty, a quadratic mesh a b c Figure: Angioplasty. Approximation on a quadratic mesh: (a) the mesh, (b) radial displacements u r on the artery, and on the balloon (c) 19 March, 2018 Solutions of problems in computational physics 42 / 44
43 Mixed FEM formulation for highly anisotropic elastic materials (with Adam Zdunek) a Figure: c b d Pressurization od a 2-layer toroidal tube (the data correspond to a rabbit carotid artery). (a) An initial p = 2 mesh, (b) detail of an h-adaptive mesh, (c) effective stress devσ 0 on a uniform p = 4 mesh, (d) effective stress on a detail of an h-adaptive mesh 19 March, 2018 Solutions of problems in computational physics 43 / 44
44 Dziȩkujȩ! 19 March, 2018 Solutions of problems in computational physics 44 / 44
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 informationFinal Ph.D. Progress Report. Integration of hp-adaptivity with a Two Grid Solver: Applications to Electromagnetics. David Pardo
Final Ph.D. Progress Report Integration of hp-adaptivity with a Two Grid Solver: Applications to Electromagnetics. David Pardo Dissertation Committee: I. Babuska, L. Demkowicz, C. Torres-Verdin, R. Van
More informationMultiscale Analysis of Vibrations of Streamers
Multiscale Analysis of Vibrations of Streamers Leszek Demkowicz Joint work with S. Prudhomme, W. Rachowicz, W. Qiu and L. Chamoin Institute for Computational Engineering and Sciences (ICES) The University
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 informationBasic concepts to start Mechanics of Materials
Basic concepts to start Mechanics of Materials Georges Cailletaud Centre des Matériaux Ecole des Mines de Paris/CNRS Notations Notations (maths) (1/2) A vector v (element of a vectorial space) can be seen
More informationFinite Elements for Large Strains - A double mixed (M 2 ) Formulation
Finite Elements for Large Strains - A double mixed (M 2 ) Formulation Motivation Development of user friendly elements robustness simple treatment of incompressible materials complex geometries geometrical
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 informationThe Finite Element Method
The Finite Element Method 3D Problems Heat Transfer and Elasticity Read: Chapter 14 CONTENTS Finite element models of 3-D Heat Transfer Finite element model of 3-D Elasticity Typical 3-D Finite Elements
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 informationReview of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis
uke University epartment of Civil and Environmental Engineering CEE 42L. Matrix Structural Analysis Henri P. Gavin Fall, 22 Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods
More 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 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 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 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 informationThe Finite Element Method for Computational Structural Mechanics
The Finite Element Method for Computational Structural Mechanics Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) January 29, 2010 Martin Kronbichler (TDB) FEM for CSM January
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 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 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 posteriori error estimation in the FEM
A posteriori error estimation in the FEM Plan 1. Introduction 2. Goal-oriented error estimates 3. Residual error estimates 3.1 Explicit 3.2 Subdomain error estimate 3.3 Self-equilibrated residuals 3.4
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 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 informationTheories of Straight Beams
EVPM3ed02 2016/6/10 7:20 page 71 #25 This is a part of the revised chapter in the new edition of the tetbook Energy Principles and Variational Methods in pplied Mechanics, which will appear in 2017. These
More informationa x Questions on Classical Solutions 1. Consider an infinite linear elastic plate with a hole as shown. Uniform shear stress
Questions on Classical Solutions. Consider an infinite linear elastic plate with a hole as shown. Uniform shear stress σ xy = T is applied at infinity. Determine the value of the stress σ θθ on the edge
More informationOn the Numerical Modelling of Orthotropic Large Strain Elastoplasticity
63 Advances in 63 On the Numerical Modelling of Orthotropic Large Strain Elastoplasticity I. Karsaj, C. Sansour and J. Soric Summary A constitutive model for orthotropic yield function at large strain
More informationA HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS
A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A. Kroker, W. Becker TU Darmstadt, Department of Mechanical Engineering, Chair of Structural Mechanics Hochschulstr. 1, D-64289 Darmstadt, Germany kroker@mechanik.tu-darmstadt.de,
More informationMATH45061: SOLUTION SHEET 1 V
1 MATH4561: SOLUTION SHEET 1 V 1.) a.) The faces of the cube remain aligned with the same coordinate planes. We assign Cartesian coordinates aligned with the original cube (x, y, z), where x, y, z 1. The
More informationNonlinear 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 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 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 informationFinite Elements for Elastic Shell Models in
Elastic s in Advisor: Matthias Heinkenschloss Computational and Applied Mathematics Rice University 13 April 2007 Outline Elasticity in Differential Geometry of Shell Geometry and Equations The Plate Model
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 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 informationCHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES
CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES * Governing equations in beam and plate bending ** Solution by superposition 1.1 From Beam Bending to Plate Bending 1.2 Governing Equations For Symmetric
More informationDuality method in limit analysis problem of non-linear elasticity
Duality method in limit analysis problem of non-linear elasticity Igor A. Brigadnov Department of Computer Science North-Western State Technical University Millionnaya 5, St. Petersburg, 191186, Russia
More informationMHA042 - Material mechanics: Duggafrågor
MHA042 - Material mechanics: Duggafrågor 1) For a static uniaxial bar problem at isothermal (Θ const.) conditions, state principle of energy conservation (first law of thermodynamics). On the basis of
More informationMéthodes de Galerkine Discontinues pour l imagerie sismique en domaine fréquentiel
Méthodes de Galerkine Discontinues pour l imagerie sismique en domaine fréquentiel Julien Diaz, Équipe Inria Magique 3D Séminaire de Mathématiques Appliquées, Laboratoire de Mathématiques Nicolas Oresme
More informationSEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by
SEMM Mechanics PhD Preliminary Exam Spring 2014 1. Consider a two-dimensional rigid motion, whose displacement field is given by u(x) = [cos(β)x 1 + sin(β)x 2 X 1 ]e 1 + [ sin(β)x 1 + cos(β)x 2 X 2 ]e
More 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 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 informationThe Plane Stress Problem
The Plane Stress Problem Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) February 2, 2010 Martin Kronbichler (TDB) The Plane Stress Problem February 2, 2010 1 / 24 Outline
More informationMixed Finite Element Methods. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016
Mixed Finite Element Methods Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016 Linear elasticity Given the load f : Ω R n, find the displacement u : Ω R n and the
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
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 informationCommon pitfalls while using FEM
Common pitfalls while using FEM J. Pamin Instytut Technologii Informatycznych w Inżynierii Lądowej Wydział Inżynierii Lądowej, Politechnika Krakowska e-mail: JPamin@L5.pk.edu.pl With thanks to: R. de Borst
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
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 informationComputational models of diamond anvil cell compression
UDC 519.6 Computational models of diamond anvil cell compression A. I. Kondrat yev Independent Researcher, 5944 St. Alban Road, Pensacola, Florida 32503, USA Abstract. Diamond anvil cells (DAC) are extensively
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 informationIntrinsic finite element modeling of a linear membrane shell problem
arxiv:3.39v [math.na] 5 Mar Intrinsic finite element modeling of a linear membrane shell problem Peter Hansbo Mats G. Larson Abstract A Galerkin finite element method for the membrane elasticity problem
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 informationPerformance comparison between hybridizable DG and classical DG methods for elastic waves simulation in harmonic domain
March 4-5, 2015 Performance comparison between hybridizable DG and classical DG methods for elastic waves simulation in harmonic domain M. Bonnasse-Gahot 1,2, H. Calandra 3, J. Diaz 1 and S. Lanteri 2
More informationA note on finite elastic deformations of fibre-reinforced non-linearly elastic tubes
Arch. Mech., 67, 1, pp. 95 109, Warszawa 2015 Brief Note A note on finite elastic deformations of fibre-reinforced non-linearly elastic tubes M. EL HAMDAOUI, J. MERODIO Department of Continuum Mechanics
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 informationA new approach for Kirchhoff-Love plates and shells
A new approach for Kirchhoff-Love plates and shells Walter Zulehner Institute of Computational Mathematics JKU Linz, Austria AANMPDE 10 October 2-6, 2017, Paleochora, Crete, Greece Walter Zulehner (JKU
More informationStructural Analysis of Large Caliber Hybrid Ceramic/Steel Gun Barrels
Structural Analysis of Large Caliber Hybrid Ceramic/Steel Gun Barrels MS Thesis Jon DeLong Department of Mechanical Engineering Clemson University OUTLINE Merger of ceramics into the conventional steel
More informationCourse Syllabus: Continuum Mechanics - ME 212A
Course Syllabus: Continuum Mechanics - ME 212A Division Course Number Course Title Academic Semester Physical Science and Engineering Division ME 212A Continuum Mechanics Fall Academic Year 2017/2018 Semester
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More information2008 by authors and 2008 Springer Science+Business Media
Antti H. Niemi, Harri Hakula, and Juhani Pitkäranta. 28. Point load on a shell. In: Karl Kunisch, Günther Of, and Olaf Steinbach (editors). Numerical Mathematics and Advanced Applications. Proceedings
More information7. Hierarchical modeling examples
7. Hierarchical modeling examples The objective of this chapter is to apply the hierarchical modeling approach discussed in Chapter 1 to three selected problems using the mathematical models studied in
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
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 informationA consistent dynamic finite element formulation for a pipe using Euler parameters
111 A consistent dynamic finite element formulation for a pipe using Euler parameters Ara Arabyan and Yaqun Jiang Department of Aerospace and Mechanical Engineering, University of Arizona, Tucson, AZ 85721,
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More information1 Exercise: Linear, incompressible Stokes flow with FE
Figure 1: Pressure and velocity solution for a sinking, fluid slab impinging on viscosity contrast problem. 1 Exercise: Linear, incompressible Stokes flow with FE Reading Hughes (2000), sec. 4.2-4.4 Dabrowski
More informationFile name: Supplementary Information Description: Supplementary Figures, Supplementary Tables and Supplementary References
File name: Supplementary Information Description: Supplementary Figures, Supplementary Tables and Supplementary References File name: Supplementary Movie 1 Description: The movie shows compression behaviour
More informationApproximation of fluid-structure interaction problems with Lagrange multiplier
Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline
More informationFREE VIBRATION OF AXIALLY LOADED FUNCTIONALLY GRADED SANDWICH BEAMS USING REFINED SHEAR DEFORMATION THEORY
FREE VIBRATION OF AXIALLY LOADED FUNCTIONALLY GRADED SANDWICH BEAMS USING REFINED SHEAR DEFORMATION THEORY Thuc P. Vo 1, Adelaja Israel Osofero 1, Marco Corradi 1, Fawad Inam 1 1 Faculty of Engineering
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 informationCOMPARISON OF PLATE MODELS FOR ANALYSIS OF LAMINATED COMPOSITES
COMPARISON OF PLATE MODELS FOR ANALYSIS OF LAMINATED COMPOSITES P. M. Mohite and C. S. Upadhyay** Department of Aerospace Engineering, IIT Kanpur 0806, INDIA, e-mail: mohite@iitk.ac.in Assistant Professor,
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 informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More informationA primer on Numerical methods for elasticity
A primer on Numerical methods for elasticity Douglas N. Arnold, University of Minnesota Complex materials: Mathematical models and numerical methods Oslo, June 10 12, 2015 One has to resort to the indignity
More informationEnergy and Momentum Conserving Algorithms in Continuum Mechanics. O. Gonzalez
Energy and Momentum Conserving Algorithms in Continuum Mechanics O. Gonzalez Département de Mathématiques École Polytechnique Fédérale de Lausanne Lausanne, Switzerland INTRODUCTION MAIN CONCERN Time integration
More informationMechanics 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 informationNumber of pages in the question paper : 05 Number of questions in the question paper : 48 Modeling Transport Phenomena of Micro-particles Note: Follow the notations used in the lectures. Symbols have their
More informationRing-shaped crack propagation in a cylinder under nonsteady cooling
High Performance Structures and Materials III 5 Ring-shaped crack propagation in a cylinder under nonsteady cooling V. A. Zhornik, Yu. A. Prokopenko, A. A. Rybinskaya & P. A. Savochka Department of Theoretical
More informationDynamics of the Mantle and Lithosphere ETH Zürich Continuum Mechanics in Geodynamics: Equation cheat sheet
Dynamics of the Mantle and Lithosphere ETH Zürich Continuum Mechanics in Geodynamics: Equation cheat sheet or all equations you will probably ever need Definitions 1. Coordinate system. x,y,z or x 1,x
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationEITN90 Radar and Remote Sensing Lecture 5: Target Reflectivity
EITN90 Radar and Remote Sensing Lecture 5: Target Reflectivity Daniel Sjöberg Department of Electrical and Information Technology Spring 2018 Outline 1 Basic reflection physics 2 Radar cross section definition
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 informationCracks Jacques Besson
Jacques Besson Centre des Matériaux UMR 7633 Mines ParisTech PSL Research University Institut Mines Télécom Aγνωστ oς Θεoς Outline 1 Some definitions 2 in a linear elastic material 3 in a plastic material
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 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 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 informationElectromagnetic Theorems
Electromagnetic Theorems Daniel S. Weile Department of Electrical and Computer Engineering University of Delaware ELEG 648 Electromagnetic Theorems Outline Outline Duality The Main Idea Electric Sources
More informationSHAPE SENSITIVITY ANALYSIS AND OPTIMIZATION FOR A CONTACT PROBLEM IN THE MANUFACTURING PROCESS DESIGN
SHAPE SENSITIVITY ANALYSIS AND OPTIMIZATION FOR A CONTACT PROBLEM IN THE MANUFACTURING PROCESS DESIGN 4 th World Congress of Structural and Multidisciplinary Optimization June 4-8, Dalian, China Nam H.
More informationHIGHER-ORDER THEORIES
HIGHER-ORDER THEORIES Third-order Shear Deformation Plate Theory Displacement and strain fields Equations of motion Navier s solution for bending Layerwise Laminate Theory Interlaminar stress and strain
More informationChemnitz Scientific Computing Preprints
Arnd Meyer The linear Naghdi shell equation in a coordinate free description CSC/13-03 Chemnitz Scientific Computing Preprints Impressum: Chemnitz Scientific Computing Preprints ISSN 1864-0087 (1995 2005:
More information16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations
6.2 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive quations Constitutive quations For elastic materials: If the relation is linear: Û σ ij = σ ij (ɛ) = ρ () ɛ ij σ ij =
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 informationRegularity Theory a Fourth Order PDE with Delta Right Hand Side
Regularity Theory a Fourth Order PDE with Delta Right Hand Side Graham Hobbs Applied PDEs Seminar, 29th October 2013 Contents Problem and Weak Formulation Example - The Biharmonic Problem Regularity Theory
More informationCIVL4332 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 informationNon-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, Politecnico di Milano, February 17, 2017, Lesson 5
Non-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, attilio.frangi@polimi.it Politecnico di Milano, February 17, 2017, Lesson 5 1 Politecnico di Milano, February 17, 2017, Lesson 5 2 Outline
More informationImplementing a Partitioned Algorithm for Fluid-Structure Interaction of Flexible Flapping Wings within Overture
10 th Symposimum on Overset Composite Grids and Solution Technology, NASA Ames Research Center Moffett Field, California, USA 1 Implementing a Partitioned Algorithm for Fluid-Structure Interaction of Flexible
More informationDeformation of bovine eye fluid structure interaction between viscoelastic vitreous, non-linear elastic lens and sclera
Karel October Tůma 24, Simulation 2018 of a bovine eye 1/19 Deformation of bovine eye fluid structure interaction between viscoelastic vitreous, non-linear elastic lens and sclera Karel Tůma 1 joint work
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 informationSymmetric Bending of Beams
Symmetric Bending of Beams beam is any long structural member on which loads act perpendicular to the longitudinal axis. Learning objectives Understand the theory, its limitations and its applications
More informationPlane and axisymmetric models in Mentat & MARC. Tutorial with some Background
Plane and axisymmetric models in Mentat & MARC Tutorial with some Background Eindhoven University of Technology Department of Mechanical Engineering Piet J.G. Schreurs Lambèrt C.A. van Breemen March 6,
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 information(2) 2. (3) 2 Using equation (3), the material time derivative of the Green-Lagrange strain tensor can be obtained as: 1 = + + +
LAGRANGIAN FORMULAION OF CONINUA Review of Continuum Kinematics he reader is referred to Belytscho et al. () for a concise review of the continuum mechanics concepts used here. he notation followed here
More information