INVERSE ANALYSIS METHODS OF IDENTIFYING CRUSTAL CHARACTERISTICS USING GPS ARRYA DATA
|
|
- Gervais Hoover
- 5 years ago
- Views:
Transcription
1 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 Institute, University of Tokyo) Contents. Stress inversion method: find equilibrating stress using measured strain and partial information on stress-strain relation. Elasticity inversion method: find local elasticity using densely measured displacement
2 GPS NETWORK AND ITS DATA 4Ê 4Ê GPS stations displacement rate 5cm/yr 3Ê 3Ê 4Ê 3Ê 3Ê 4Ê
3 NEED FOR LOCAL STRESS PREDICTION Material Test of Next Generation Earthquake Prediction GPS data strain σ σ ε local ε σ relation ε regional stress 3
4 IS STRESS INVERSION POSSIBLE? 3D State Most Difficult D State Possible? 3 UNKNOWNS σ σ σ EQUATIONS equilibrium in x -direction equilibrium in x -direction CONDITION STRAIN 4
5 STRESS INVERSION σ σ σ Airy = a = a,, = a, partial information: σ + σ = f ( ε σ f ( ε ij + σ ) = κ( ε = κ( ε + ε + ε ) ) ij ) Poisson equation self-equilibrating stress G.E. B.C. B.V.P. a n, a + a, +, n = κ( ε a, = n + ε r + ) n r D bulk modulus? measured strain resultant forces 5
6 EXTENSION TO OTHER DEFORMATION STATE Dynamic State σ σ σ,, + σ + σ + σ,, = = ρ&& u = ρ&& u f ( ε) Finite Deformation State: σ j,k j,k ( X ( X + σ k k x x j j ) σ ) σ = f ( ε) j,k j,k = = ( σ = σ (X), x = x (X)) ij ij i i 6
7 NUMERICAL SIMULATION Conditions elasto-plastic material with unknown yield function prediction of stress and stress-strain relation distributed force displacement field is used as input data.. FEM computation with x48 elements 7
8 RESULTS OF INVERSION σ σ distribution of σ exact sample surface predicted principle stressplastic strain rate σ σ 3 σ σ 3 σ σ 3 σ σ 3 good agreement exact yield locus predicted 8
9 MODEL EXPERIMENT 7cm Riedel shears Torsional Shearing 8cm examples of image video recorder CCD camera Experiment Apparatus servomotor max. shear strain 9
10 OVERALL STRESS-STRAIN RELATION stress (kgf/m 3 ) time (sec) stress (kgf/m 3 ) rotation (radian) strain will be different from local relations? time (sec)
11 LOCAL STRESS-STRAIN RELATIONS σ [kpa] 4 σ [kpa] D A D A normal stress-strain ε shear stress-strain ε A: far from crack D: near crack no common relations?
12 RESULTS OF INVERSION max. shear stress (kpa) 3 yielding A B C D E max. shear strain common elasto-plastic relations?
13 APPLICATION TO GPS ARRYA DATA verification of numerical analysis method check numerical stability of solving boundary value problem check dependency of parameters application of stress inversion method geophysical interpretation of analysis results critical examination of assumption of plane state development of crust deformation monitor automatic processing of GPS array data 3
14 CONVERGENCE σ (hydrostatic stress) Δ: resolution of strain distribution (degree) 4
15 EFFECT OF REFERENCE τ (max. shear stress) 5
16 COMPARISON OF STRESS WITH STRAIN σ (hydrostatic stress) τ (max. shear stress) ε (vol. strain) γ (max. shear strain) 6
17 REGIONAL CONSTITUTIVE RELATIONS τ/γ regional stiffness (τ, γ: max. shear stress and strain) φ θ regional anisotropy (φ, γ: principle stress and strain) regional heterogeneity and anisotropy 7
18 CHANGE IN INVARIANT st invariant nd invariant 8
19 CHANGE IN REGIONAL STATE st invariant 5 ε σ nd invariant 5 5 ε σ A γ τ γ τ C 5 ε σ 5 5 B 3 5 Y A B γ τ ε σ C -5 3 γ τ 9
20 REGIONAL STRESS AND STRAIN st invariant 5 σ ε nd invariant τ γ 75 ε γ 5 γ τ A C 5 σ ε τ γ 75 B A -5-5 B σ ε t-g C
21 GPS DATA DURING GPS Data no spatial filtering to get rid of measurement noise linear interpolation between two GPS station More Sophisticated Treatment of BVP FEM with triangle element weak form regionally averaged field quantities
22 APPLICATION TO GPS NETWORK DATA BVP in Rate Form and Weak Form G.E. B.C. a&, n a& + a&,, + n a& ϕ, a&, + ϕ,a&, κ( ϕ,u& + ϕ,u& )ds = Computation of Average Quantities =κ( ε&, +ε& = n r& ) + n r& strain rate displacement rate ε& σ& = Ω = Ω Ω Ω n n u& a& dl, dl average stress rate computed by using st -order derivative
23 GPS NETWORK AND ITS DATA 4Ê 4Ê GPS stations displacement rate 5cm/yr 3Ê 3Ê 4Ê 3Ê 3Ê 4Ê 3
24 REGIONAL STRAIN RATE 6e-7 6e-7 4Ê 4e-7 4Ê e-7 4e-7 -e-7 e-7-4e-7 volumetric -6e-7 Strain max. shear Strain 3Ê 3Ê 4Ê 3Ê 3Ê 4Ê 4
25 REGIONAL STRESS RATE 6e-7 6e-7 4Ê 4e-7 4Ê e-7 4e-7 -e-7 e-7-4e-7 volumetric -6e-7 x κ [Pa] max. shear x κ [Pa] 3Ê 3Ê 4Ê 3Ê 3Ê 4Ê 5
26 COMPARISON WITH SEISMIC EVENTS? 6e-7 4Ê 4e-7 4Ê e-7 -e-7-4e-7 volumetric -6e-7 >[km], >M3 x κ [Pa] 3Ê 3Ê 4Ê 3Ê 3Ê 4Ê 6
27 REGIONAL CONSTITUTIVE RELATIONS findκ( x) s.t. τ& ( x; κ) = κ( x) γ& ( x). τ &, γ& maximum shear 4.5 κ regional stiffness. κ is originally used to relate σ & ε through σ = κ ε not too far from known geological structure
28 DRAWBACKS OF STRESS INVERSION Need to Know One Constitutive Relation bulk stress and bulk strain isotropy assumption Need to Know Boundary Traction/Resultant Force assumption of uniform stress fast decrease of non-uniform boundary traction Difficulty in Understanding Plane-Stress-State Model another analysis method needed? 8
29 DRAWBACKS OF ELASTICITY INVERSION Sensitive to Displacement Error need to make fine discretization of target body need to have some strong modes of deformation Why is it so? no mistakes in mathematics poor understanding of physics 9
30 PHYSICAL PROCESS AND MEASUREMENT GIS data data inverse analysis May not be good to pose an inverse problem from data to characteristics, because a path from characteristics to data has physical process and measurement measurement X unknown mechanical property characteristics response crust deformation physical process source unknown plate slip/ underground deformation 3
31 ESTIMATION AND INVERSION estimation data measurement inversion estimation: from data to response, i.e., estimate function for displacement from data which are measured discretely inversion: from response to characteristics, i.e., find most suitable characteristics for physical process even though source is not known characteristics response physical process source 3
32 BLOCK IN CONTINUUM nodes at which displacement is measured inner node body B x boundary node x principal of elasticity inversion block S k k displacement at all nodes given Δ Δ Δ Δ = k k no need to consider interaction of S with outside region 3
33 IDENTIFICATION OF DISPLACEMENT MODE concentrated force node block hexagonal block block block 3 material sample material sample test apply several BC s, and measure displacement at nodes of a hexagonal block.. identify displacement modes (a characteristic set of nodal displacement). identify local elasticity 33
34 IDENTIFICATION OF DISPLACEMENT MODE displacement mode elastic parameters c c.5 constraint 3 c 4 c 5 c 6 c block g. g block block unconstraint
35 ELASTICITY INVERSION METHOD body B local field variables allows Taylor series expansion block Ω nodes at which displacement is measured BASIC PROCEDURES OF INVERSION. use displacement data to determine Tayler series expansion coefficients of displacement. use equilibrium equation to estimate elastic parameters by expanding stress in Talker series uniform Poisson's ratio is assumed 35
36 DETERMINATION OF DISPLACEMENT COEFFICIENT. Taylor Expansion: { a ip } u = P i ( x) a ipf p ( x) p= {a {f ip p } = {u i } = {, x,u i,, x,u, x i,,, x u x i,,u, x i,,, L} u i,,u i,, L}. Displacement Data: u n i = P p= f pn a ip { u n i } n 3. Solution of Matrix Equation f = f ( x ) pn p a ip = A α= λ α n u n i ϕ α n ψ α p + P A β= b β i ψ β p { λ } α, ϕ α ψ α n, p : SVD of f pn fully determined undetermined 36
37 37. Elasticity Tensor Expressed in Terms of Poisson Ratio. Equation of Equilibrium and Its Taylor Expansion 3. Coefficient of Expansion: b ip = for th Order (p=) ESTIMATION OF POISSON RATIO ν ESTIMATION OF POISSON RATIO ν + ν = a a a b b M ) ( f b )) ( u (c ) ( b p p ip j, k,l ijkl i = = = x x x ijkl ijkl ijkl c c c ν + = linear equation of ν is derived
38 NUMERICAL SIMULATION (). error error..5 nd order 3rd order.5 nd order 3rd order ν ν=.5 -. ν ν=.35 ν=.5 ν=.35 frequency 8 6 measurement 4 expansion of displacement expansion of equilibrium 3 rd order th or st order ν 38
39 NUMERICAL SIMULATION (). Measurement Error: { e n i } th. P n u i = f pna ip + p= e n i e'.8.4 ν exact =.35 ν exact =.5. Find ν such that ν minimize subjected to e = (e b ( ν) = ip n i ) st e measurement expansion of displacement 3 rd order expansion of equilibrium th or st order ν exact = ν 39
40 APPLICATION OF LOCAL GPS DATA pre-wte 35.5 Velocity field Pre-Tottori Earthquake cm/year GPS array Western Tottori Earthquake (M JMA =7.3) examine change in deformation and elasticity before and after this earthquake 53 GPS observation points 8 triangular elements GPS data obtained from 997 to annual and biannual sinusoidal variations excluded post-wte Velocity field Post-Tottori Earthquake 3cm/year
41 STRAIN RATE pre-wte Dilatational strain rate post-wte Dilatational strain rate e e-7 4e-7 4e-7 dilatational 35 e-7 35 e e e-7-4e-7-4e e-7 [Strain/Year] 34-6e-7 [Strain/Year] Maximum shear strain rate Maximum shear strain rate e e-6 8e-7 8e-7 shear 35 6e e-7 4e-7 4e e-7 e-7 34 [Strain/Year] 34 [Strain/Year]
42 STRESS RATE pre-wte Dilatational stress increment post-wte Dilatational stress increment dilatational [KPa/Year] 34 - [KPa/Year] Maximum shear stress increment Maximum shear stress increment shear [KPa/Year] 34 [KPa/Year]
43 PRINICPLE AXIS pre-wte Principal axes of strain rate post-wte Principal axes of strain rate strain (extensional) (contractional).5 [Micro Strain/Year] (extensional) (contractional).5 [Micro Strain/Year] Principal axes of stress rate Principal axes of stress rate stress (extensional) (contractional) [KPa/Year] (extensional) (contractional) [KPa/Year]
44 POISSON RATIO pre-wte Poisson s ratio Pre-Tottori Earthquake Mt. Daisen post-wte Poisson s ratio Post-Tottori Earthquake Mt. Daisen Naka-Umi Naka-Umi Change of Poisson s ratio Mt. Daisen Naka-Umi difference ν has been reduced near the source faults. The comparison with other analyses are being made [%]
45 CONCLUDING REMARKS Two inverse analysis methods stress inversion find Airy s stress function by solving Poisson s equation elasticity inversion find elastic parameters by estimating displacement expansion coefficients Development of new inversion is needed for geophysics where experiments cannot be made. Application small material samples used for bio-mechanics geomaterials new image analysis with higher spatial resolution 45
Finite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More 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 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 informationUnderstand basic stress-strain response of engineering materials.
Module 3 Constitutive quations Learning Objectives Understand basic stress-strain response of engineering materials. Quantify the linear elastic stress-strain response in terms of tensorial quantities
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 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 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 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 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 informationElastoplastic Deformation in a Wedge-Shaped Plate Caused By a Subducting Seamount
Elastoplastic Deformation in a Wedge-Shaped Plate Caused By a Subducting Seamount Min Ding* 1, and Jian Lin 2 1 MIT/WHOI Joint Program, 2 Woods Hole Oceanographic Institution *Woods Hole Oceanographic
More informationMohr's Circle and Earth Stress (The Elastic Earth)
Lect. 1 - Mohr s Circle and Earth Stress 6 Mohr's Circle and Earth Stress (The Elastic Earth) In the equations that we derived for Mohr s circle, we measured the angle, θ, as the angle between σ 1 and
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 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 informationHÅLLFASTHETSLÄRA, LTH Examination in computational materials modeling
HÅLLFASTHETSLÄRA, LTH Examination in computational materials modeling TID: 2016-28-10, kl 14.00-19.00 Maximalt 60 poäng kan erhållas på denna tenta. För godkänt krävs 30 poäng. Tillåtet hjälpmedel: räknare
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 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 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 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 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 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 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 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 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 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 short review of continuum mechanics
A short review of continuum mechanics Professor Anette M. Karlsson, Department of Mechanical ngineering, UD September, 006 This is a short and arbitrary review of continuum mechanics. Most of this material
More informationMODELING OF CONCRETE MATERIALS AND STRUCTURES. Kaspar Willam. Isotropic Elastic Models: Invariant vs Principal Formulations
MODELING OF CONCRETE MATERIALS AND STRUCTURES Kaspar Willam University of Colorado at Boulder Class Meeting #2: Nonlinear Elastic Models Isotropic Elastic Models: Invariant vs Principal Formulations Elastic
More informationMicro-mechanics in Geotechnical Engineering
Micro-mechanics in Geotechnical Engineering Chung R. Song Department of Civil Engineering The University of Mississippi University, MS 38677 Fundamental Concepts Macro-behavior of a material is the average
More informationNumerical Modelling of Dynamic Earth Force Transmission to Underground Structures
Numerical Modelling of Dynamic Earth Force Transmission to Underground Structures N. Kodama Waseda Institute for Advanced Study, Waseda University, Japan K. Komiya Chiba Institute of Technology, Japan
More information1.050 Engineering Mechanics. Lecture 22: Isotropic elasticity
1.050 Engineering Mechanics Lecture 22: Isotropic elasticity 1.050 Content overview I. Dimensional analysis 1. On monsters, mice and mushrooms 2. Similarity relations: Important engineering tools II. Stresses
More informationGEO E1050 Finite Element Method Mohr-Coulomb and other constitutive models. Wojciech Sołowski
GEO E050 Finite Element Method Mohr-Coulomb and other constitutive models Wojciech Sołowski To learn today. Reminder elasticity 2. Elastic perfectly plastic theory: concept 3. Specific elastic-perfectly
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 informationSymmetry and Properties of Crystals (MSE638) Stress and Strain Tensor
Symmetry and Properties of Crystals (MSE638) Stress and Strain Tensor Somnath Bhowmick Materials Science and Engineering, IIT Kanpur April 6, 2018 Tensile test and Hooke s Law Upto certain strain (0.75),
More informationAnalysis of a portal steel frame subject to fire by use of a truss model
Analysis of a portal steel frame subject to fire by use of a truss model P. G. Papadopoulos & A. Mathiopoulou Department of Civil Engineering, Aristotle University of Thessaloniki, Greece Abstract A plane
More informationME FINITE ELEMENT ANALYSIS FORMULAS
ME 2353 - FINITE ELEMENT ANALYSIS FORMULAS UNIT I FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PROBLEMS 01. Global Equation for Force Vector, {F} = [K] {u} {F} = Global Force Vector [K] = Global Stiffness
More 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 informationSurface force on a volume element.
STRESS and STRAIN Reading: Section. of Stein and Wysession. In this section, we will see how Newton s second law and Generalized Hooke s law can be used to characterize the response of continuous medium
More informationConstitutive modelling of fabric anisotropy in sand
Geomechanics from Micro to Macro Soga et al. (Eds) 2015 Taylor & Francis Group, London, ISBN 978-1-138-02707-7 Constitutive modelling of fabric anisotropy in sand Z.W. Gao School of Engineering, University
More informationFig. 1. Circular fiber and interphase between the fiber and the matrix.
Finite element unit cell model based on ABAQUS for fiber reinforced composites Tian Tang Composites Manufacturing & Simulation Center, Purdue University West Lafayette, IN 47906 1. Problem Statement In
More informationGeology 229 Engineering Geology. Lecture 5. Engineering Properties of Rocks (West, Ch. 6)
Geology 229 Engineering Geology Lecture 5 Engineering Properties of Rocks (West, Ch. 6) Common mechanic properties: Density; Elastic properties: - elastic modulii Outline of this Lecture 1. Uniaxial rock
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 informationEshan V. Dave, Secretary of M&FGM2006 (Hawaii) Research Assistant and Ph.D. Candidate. Glaucio H. Paulino, Chairman of M&FGM2006 (Hawaii)
Asphalt Pavement Aging and Temperature Dependent Properties through a Functionally Graded Viscoelastic Model I: Development, Implementation and Verification Eshan V. Dave, Secretary of M&FGM2006 (Hawaii)
More 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 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 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 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 informationReciprocal of the initial shear stiffness of the interface K si under initial loading; reciprocal of the initial tangent modulus E i of the soil
Appendix F Notation a b B C c C k C N C s C u C wt C θ D r D 1 D 2 D 10 D 30 Reciprocal of the initial shear stiffness of the interface K si under initial loading; reciprocal of the initial tangent modulus
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 information1 Nonlinear deformation
NONLINEAR TRUSS 1 Nonlinear deformation When deformation and/or rotation of the truss are large, various strains and stresses can be defined and related by material laws. The material behavior can be expected
More informationBasic Concepts of Strain and Tilt. Evelyn Roeloffs, USGS June 2008
Basic Concepts of Strain and Tilt Evelyn Roeloffs, USGS June 2008 1 Coordinates Right-handed coordinate system, with positions along the three axes specified by x,y,z. x,y will usually be horizontal, and
More informationSummary so far. Geological structures Earthquakes and their mechanisms Continuous versus block-like behavior Link with dynamics?
Summary so far Geodetic measurements velocities velocity gradient tensor (spatial derivatives of velocity) Velocity gradient tensor = strain rate (sym.) + rotation rate (antisym.) Strain rate tensor can
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 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 informationChapter 5 Elastic Strain, Deflection, and Stability 1. Elastic Stress-Strain Relationship
Chapter 5 Elastic Strain, Deflection, and Stability Elastic Stress-Strain Relationship A stress in the x-direction causes a strain in the x-direction by σ x also causes a strain in the y-direction & z-direction
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 informationMECHANICS OF MATERIALS. EQUATIONS AND THEOREMS
1 MECHANICS OF MATERIALS. EQUATIONS AND THEOREMS Version 2011-01-14 Stress tensor Definition of traction vector (1) Cauchy theorem (2) Equilibrium (3) Invariants (4) (5) (6) or, written in terms of principal
More informationInfluences of material dilatancy and pore water pressure on stability factor of shallow tunnels
Influences of material dilatancy and pore water pressure on stability factor of shallow tunnels YANG Xiao-li( ), HUANG Fu( ) School of Civil and Architectural Engineering, Central South University, Changsha
More informationPractice Final Examination. Please initial the statement below to show that you have read it
EN175: Advanced Mechanics of Solids Practice Final Examination School of Engineering Brown University NAME: General Instructions No collaboration of any kind is permitted on this examination. You may use
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 informationGravity Tectonics Volcanism Atmosphere Water Winds Chemistry. Planetary Surfaces
Gravity Tectonics Volcanism Atmosphere Water Winds Chemistry Planetary Surfaces Gravity & Rotation Polar flattening caused by rotation is the largest deviation from a sphere for a planet sized object (as
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting Lectures & 3, 9/31 Aug 017 www.geosc.psu.edu/courses/geosc508 Discussion of Handin, JGR, 1969 and Chapter 1 Scholz, 00. Stress analysis and Mohr Circles Coulomb Failure
More information1. classic definition = study of deformed rocks in the upper crust
Structural Geology I. Introduction 1. classic definition = study of deformed rocks in the upper crust deformed includes translation, rotation, and strain (change of shape) All rocks are deformed in some
More informationMonte-Carlo Simulation of Failure Phenomena using Particle Discretization
Monte-Carlo Simulation of Failure Phenomena using Particle Discretization Kenji OGUNI Atsushi WAKAI Muneo HORI Earthquake Research Institute, University of Tokyo Motivation Deformation Failure x a PDF
More informationMathematical Background
CHAPTER ONE Mathematical Background This book assumes a background in the fundamentals of solid mechanics and the mechanical behavior of materials, including elasticity, plasticity, and friction. A previous
More informationHomework Problems. ( σ 11 + σ 22 ) 2. cos (θ /2), ( σ θθ σ rr ) 2. ( σ 22 σ 11 ) 2
Engineering Sciences 47: Fracture Mechanics J. R. Rice, 1991 Homework Problems 1) Assuming that the stress field near a crack tip in a linear elastic solid is singular in the form σ ij = rλ Σ ij (θ), it
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 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 informationGEOMETRIC NONLINEAR ANALYSIS
GEOMETRIC NONLINEAR ANALYSIS The approach for solving problems with geometric nonlinearity is presented. The ESAComp solution relies on Elmer open-source computational tool [1] for multiphysics problems.
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 informationMechanical Engineering Ph.D. Preliminary Qualifying Examination Solid Mechanics February 25, 2002
student personal identification (ID) number on each sheet. Do not write your name on any sheet. #1. A homogeneous, isotropic, linear elastic bar has rectangular cross sectional area A, modulus of elasticity
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 informationVORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS
The 4 th World Conference on Earthquake Engineering October -7, 008, Beijing, China VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS K.
More informationLecture 14: Strain Examples. GEOS 655 Tectonic Geodesy Jeff Freymueller
Lecture 14: Strain Examples GEOS 655 Tectonic Geodesy Jeff Freymueller A Worked Example Consider this case of pure shear deformation, and two vectors dx 1 and dx 2. How do they rotate? We ll look at vector
More informationStress and Strains in Soil and Rock. Hsin-yu Shan Dept. of Civil Engineering National Chiao Tung University
Stress and Strains in Soil and Rock Hsin-yu Shan Dept. of Civil Engineering National Chiao Tung University Stress and Strain ε 1 1 2 ε 2 ε Dimension 1 2 0 ε ε ε 0 1 2 ε 1 1 2 ε 2 ε Plane Strain = 0 1 2
More informationA thermo-hydro-mechanically coupled analysis of clay using a thermo-elasto-viscoplastic model
JHUWS05 A thermo-hydro-mechanically coupled analysis of clay using a thermo-elasto-viscoplastic model by F. Oka, S. Kimoto, Y.-S. Kim, N. Takada Department of Civil & Earth Resources Engineering, Kyoto
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 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 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 informationNONLINEAR 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 information3D and 2D Formulations of Incremental Stress-Strain Relations for Granular Soils
Archives of Hydro-Engineering and Environmental Mechanics Vol. 55 (2008), No. 1 2, pp. 45 5 IBW PAN, ISSN 121 726 Technical Note D and 2D Formulations of Incremental Stress-Strain Relations for Granular
More informationTowards Efficient Finite Element Model Review Dr. Richard Witasse, Plaxis bv (based on the original presentation of Dr.
Towards Efficient Finite Element Model Review Dr. Richard Witasse, Plaxis bv (based on the original presentation of Dr. Brinkgreve) Journée Technique du CFMS, 16 Mars 2011, Paris 1/32 Topics FEA in geotechnical
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 informationUsing MATLAB and. Abaqus. Finite Element Analysis. Introduction to. Amar Khennane. Taylor & Francis Croup. Taylor & Francis Croup,
Introduction to Finite Element Analysis Using MATLAB and Abaqus Amar Khennane Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Croup, an informa business
More informationReference material Reference books: Y.C. Fung, "Foundations of Solid Mechanics", Prentice Hall R. Hill, "The mathematical theory of plasticity",
Reference material Reference books: Y.C. Fung, "Foundations of Solid Mechanics", Prentice Hall R. Hill, "The mathematical theory of plasticity", Oxford University Press, Oxford. J. Lubliner, "Plasticity
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 informationUniversity of Sheffield The development of finite elements for 3D structural analysis in fire
The development of finite elements for 3D structural analysis in fire Chaoming Yu, I. W. Burgess, Z. Huang, R. J. Plank Department of Civil and Structural Engineering StiFF 05/09/2006 3D composite structures
More informationPost Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method
9210-220 Post Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method You should have the following for this examination one answer book scientific calculator No
More informationNUMERICAL ANALYSIS OF A PILE SUBJECTED TO LATERAL LOADS
IGC 009, Guntur, INDIA NUMERICAL ANALYSIS OF A PILE SUBJECTED TO LATERAL LOADS Mohammed Younus Ahmed Graduate Student, Earthquake Engineering Research Center, IIIT Hyderabad, Gachibowli, Hyderabad 3, India.
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 informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Surface and body forces Tensors, Mohr circles. Theoretical strength of materials Defects Stress concentrations Griffith failure
More information3D and Planar Constitutive Relations
3D and Planar Constitutive Relations A School on Mechanics of Fibre Reinforced Polymer Composites Knowledge Incubation for TEQIP Indian Institute of Technology Kanpur PM Mohite Department of Aerospace
More informationModule-4. Mechanical Properties of Metals
Module-4 Mechanical Properties of Metals Contents ) Elastic deformation and Plastic deformation ) Interpretation of tensile stress-strain curves 3) Yielding under multi-axial stress, Yield criteria, Macroscopic
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 informationIntroduction to Seismology Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 1.510 Introduction to Seismology Spring 008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 1.510 Introduction to
More informationWave Propagation in Fractured Poroelastic Media
Wave Propagation in Fractured Poroelastic Media WCCM, MS170: Advanced Computational Techniques in Geophysical Sciences, Barcelona, Spain, July 2014 Juan E. Santos Instituto del Gas y del Petróleo (IGPUBA),
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 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 informationComputational Inelasticity FHLN05. Assignment A non-linear elasto-plastic problem
Computational Inelasticity FHLN05 Assignment 2017 A non-linear elasto-plastic problem General instructions A written report should be submitted to the Division of Solid Mechanics no later than October
More informationStress intensity factors for an inclined and/or eccentric crack in a finite orthotropic lamina
1886 Stress intensity factors for an inclined and/or eccentric crack in a finite orthotropic lamina Abstract Stress intensity factors (SIF) are determined for an inclined and / or eccentric crack in a
More informationTheory of Plasticity. Lecture Notes
Theory of Plasticity Lecture Notes Spring 2012 Contents I Theory of Plasticity 1 1 Mechanical Theory of Plasticity 2 1.1 Field Equations for A Mechanical Theory.................... 2 1.1.1 Strain-displacement
More informationEMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading
MA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading MA 3702 Mechanics & Materials Science Zhe Cheng (2018) 2 Stress & Strain - Axial Loading Statics
More information