The Hyperelastic material is examined in this section.

Similar documents
Grand Canonical Ensemble

3.4 Properties of the Stress Tensor

Jones vector & matrices

From Structural Analysis to FEM. Dhiman Basu

September 27, Introduction to Ordinary Differential Equations. ME 501A Seminar in Engineering Analysis Page 1. Outline

ACOUSTIC WAVE EQUATION. Contents INTRODUCTION BULK MODULUS AND LAMÉ S PARAMETERS

8-node quadrilateral element. Numerical integration

Linear Algebra Provides a Basis for Elasticity without Stress or Strain

Economics 600: August, 2007 Dynamic Part: Problem Set 5. Problems on Differential Equations and Continuous Time Optimization

Fundamentals of Continuum Mechanics. Seoul National University Graphics & Media Lab

167 T componnt oftforc on atom B can b drvd as: F B =, E =,K (, ) (.2) wr w av usd 2 = ( ) =2 (.3) T scond drvatv: 2 E = K (, ) = K (1, ) + 3 (.4).2.2

A Note on Estimability in Linear Models

Lecture 14. Relic neutrinos Temperature at neutrino decoupling and today Effective degeneracy factor Neutrino mass limits Saha equation

COHORT MBA. Exponential function. MATH review (part2) by Lucian Mitroiu. The LOG and EXP functions. Properties: e e. lim.

Lecture 3: Phasor notation, Transfer Functions. Context

Heisenberg Model. Sayed Mohammad Mahdi Sadrnezhaad. Supervisor: Prof. Abdollah Langari

Stress-Based Finite Element Methods for Dynamics Analysis of Euler-Bernoulli Beams with Various Boundary Conditions

Basic Electrical Engineering for Welding [ ] --- Introduction ---

Electrochemical Equilibrium Electromotive Force. Relation between chemical and electric driving forces

Einstein Equations for Tetrad Fields

Optimal Ordering Policy in a Two-Level Supply Chain with Budget Constraint

From Structural Analysis to Finite Element Method

CHAPTER 33: PARTICLE PHYSICS

Quasi-Classical States of the Simple Harmonic Oscillator

Review - Probabilistic Classification

VISUALIZATION OF DIFFERENTIAL GEOMETRY UDC 514.7(045) : : Eberhard Malkowsky 1, Vesna Veličković 2

External Equivalent. EE 521 Analysis of Power Systems. Chen-Ching Liu, Boeing Distinguished Professor Washington State University

Section 11.6: Directional Derivatives and the Gradient Vector

Ch. 24 Molecular Reaction Dynamics 1. Collision Theory

Abstract Interpretation: concrete and abstract semantics

Division of Mechanics Lund University MULTIBODY DYNAMICS. Examination Name (write in block letters):.

Discrete Shells Simulation

That is, we start with a general matrix: And end with a simpler matrix:

te Finance (4th Edition), July 2017.

Green Functions, the Generating Functional and Propagators in the Canonical Quantization Approach

9.5 Complex variables

Stretching and bending deformations due to normal and shear tractions of doubly curved shells using third-order shear and normal deformable theory

COMPLEX NUMBER PAIRWISE COMPARISON AND COMPLEX NUMBER AHP

Folding of Regular CW-Complexes

Relate p and T at equilibrium between two phases. An open system where a new phase may form or a new component can be added

The van der Waals interaction 1 D. E. Soper 2 University of Oregon 20 April 2012

First looking at the scalar potential term, suppose that the displacement is given by u = φ. If one can find a scalar φ such that u = φ. u x.

14. MODELING OF THIN-WALLED SHELLS AND PLATES. INTRODUCTION TO THE THEORY OF SHELL FINITE ELEMENT MODELS

Self-Adjointness and Its Relationship to Quantum Mechanics. Ronald I. Frank 2016

CHAPTER 7d. DIFFERENTIATION AND INTEGRATION

:2;$-$(01*%<*=,-./-*=0;"%/;"-*

ME 321 Kinematics and Dynamics of Machines S. Lambert Winter 2002

A NON-LINEAR MODEL FOR STUDYING THE MOTION OF A HUMAN BODY. Piteşti, , Romania 2 Department of Automotive, University of Piteşti

Consider a system of 2 simultaneous first order linear equations

6 Finite element methods for the Euler Bernoulli beam problem

Derivation of Electron-Electron Interaction Terms in the Multi-Electron Hamiltonian

A general N-dimensional vector consists of N values. They can be arranged as a column or a row and can be real or complex.

Addition of angular momentum

Static/Dynamic Deformation with Finite Element Method. Graphics & Media Lab Seoul National University

Math 34A. Final Review

The Matrix Exponential

10/7/14. Mixture Models. Comp 135 Introduction to Machine Learning and Data Mining. Maximum likelihood estimation. Mixture of Normals in 1D

EXST Regression Techniques Page 1

5.80 Small-Molecule Spectroscopy and Dynamics

Lucas Test is based on Euler s theorem which states that if n is any integer and a is coprime to n, then a φ(n) 1modn.

Mathematical Model of Arterial Hemodynamics, Description, Computer Implementation, Results Comparison

The Matrix Exponential

1973 AP Calculus AB: Section I

Pipe flow friction, small vs. big pipes

Ερωτήσεις και ασκησεις Κεφ. 10 (για μόρια) ΠΑΡΑΔΟΣΗ 29/11/2016. (d)

The Fourier Transform

2.10 Convected Coordinates

Basic Polyhedral theory

ON THE INTEGRAL INVARIANTS OF KINEMATICALLY GENERATED RULED SURFACES *

NOTES FOR CHAPTER 17. THE BOLTZMANN FACTOR AND PARTITION FUNCTIONS. Equilibrium statistical mechanics (aka statistical thermodynamics) deals with the

Econ107 Applied Econometrics Topic 10: Dummy Dependent Variable (Studenmund, Chapter 13)

Some Useful Formulae

Lie Groups HW7. Wang Shuai. November 2015

Finite element discretization of Laplace and Poisson equations

In this section is given an overview of the common elasticity models.

Chemical Physics II. More Stat. Thermo Kinetics Protein Folding...

Phys 774: Nonlinear Spectroscopy: SHG and Raman Scattering

Physics 256: Lecture 2. Physics

ANALYSIS: The mass rate balance for the one-inlet, one-exit control volume at steady state is

10. The Discrete-Time Fourier Transform (DTFT)

Outlier-tolerant parameter estimation

INTERFACE CORNERS IN ANISOTROPIC/PIEZOELECTRIC/ VISCOELASTIC MATERIALS

Addition of angular momentum

u 3 = u 3 (x 1, x 2, x 3 )

Supplementary Materials

IV. First Law of Thermodynamics. Cooler. IV. First Law of Thermodynamics

HORIZONTAL IMPEDANCE FUNCTION OF SINGLE PILE IN SOIL LAYER WITH VARIABLE PROPERTIES

1997 AP Calculus AB: Section I, Part A

y = 2xe x + x 2 e x at (0, 3). solution: Since y is implicitly related to x we have to use implicit differentiation: 3 6y = 0 y = 1 2 x ln(b) ln(b)

Differentiation of Exponential Functions

Probability and Stochastic Processes: A Friendly Introduction for Electrical and Computer Engineers Roy D. Yates and David J.

Network Congestion Games

2F1120 Spektrala transformer för Media Solutions to Steiglitz, Chapter 1

1.9 Cartesian Tensors

Fakultät III Univ.-Prof. Dr. Jan Franke-Viebach

cycle that does not cross any edges (including its own), then it has at least

Cramér-Rao Inequality: Let f(x; θ) be a probability density function with continuous parameter

Summary: Solving a Homogeneous System of Two Linear First Order Equations in Two Unknowns

Direct Approach for Discrete Systems One-Dimensional Elements

The graph of y = x (or y = ) consists of two branches, As x 0, y + ; as x 0, y +. x = 0 is the

Transcription:

4. Hyprlastcty h Hyprlastc matral s xad n ths scton. 4..1 Consttutv Equatons h rat of chang of ntrnal nrgy W pr unt rfrnc volum s gvn by th strss powr, whch can b xprssd n a numbr of dffrnt ways (s 3.7.6): W & Jσ : d P : F& S : E& S : C& (4..1) h ntrnal nrgy s rgardd as a functon of a dformaton varabl. For xampl, h chang n nrgy s du to th dformaton whch taks plac, so tak W to b a functon of, say, th dformaton gradnt F (t), W (F). It s assumd that n th rfrnc confguraton th stran nrgy s zro, W ( I) 0, and that t grows wth dformaton, W 0 1. h chan rul gvs From 4..1, & ( F ) W W : F & (4..) F P F& ( ) : W : F& (4..3) Snc F and F & can tak on any valu ndpndnt of th othr, on must hav P (4..4) hs s a consttutv quaton rlatng th knmatc varabls to th forc varabls. From 4..1, altrnatv forms ar: ( E) S, S (4..5) h procdur usd n Eqns. 4..-4 cannot b usd for th strss powr rlaton J σ : d snc thr s no functon whos drvatv s th rat of dformaton. Instad, us th rlaton 3.5.6, σ J 1 PF, 4..4, P /, to gt 1 and that t tnds to nfnty as th matral s thr comprssd to a pont, J dt F 0 or xpandd to an nfnt rang, J dt F. 360

σ J F, σ J F m (4..6) An altrnatv rlaton can b obtand by frst drvng a rlatonshp btwn th partal drvatvs of th stran nrgy functon wth rspct to th dformaton gradnt, W /, and wth rspct to th rght Cauchy-Grn tnsor, W /. Suppos frst that th stran nrgy s a functon of th dformaton gradnt: F (4..7) h chan rul for W ( C(F) ) W gvs (4..8) Wth C F F, km F kn + F km kn δ δ F + δ δ F δ F + δ F k m kn k n km m n n (4..9) so that (and usng th fact that C s symmtrc), δ F F m n n F n nm m m + F + F + δ F n n n (4..10) or F, F k k (4..11) Now 4..6 can b r-wrttn as σ F J F (4..1) Slarly, Eqn. 4..4 cab b r-wrttn as P F (4..13) 361

h Stran Enrgy and th Rght Strtch nsor h rght strtch tnsor can b xprssd as C UU, whr U s th rght strtch, or, snc U s also symmtrc, C U U. On can s th slarty btwn ths rlaton and th rlaton C F F so, usng th sam arguts as gvn abov, on has for W W C(U), (s 1.11.36) ( ) U (4..14) U Snc U s symmtrc, U U U U U (4..15) showng that U and W / ar coaxal. 4.. Obctvty of th Consttutv Equatons * An obsrvr transformaton (s.8) rsults n F QF * * W ( C ) / /, and, so, from 4..1, and * C C, so * σ J QF F Q QσQ (4..16) whch s th obctvty rqurt for a spatal tnsor and so ths consttutv law satsfs th rqurt of matral fram-ndffrnc, whn W s a functon of C. hs vdntly holds tru also whn W s a functon of E. It dos not, howvr, hold tru n gnral whn W s a functon of F, as n th consttutv law P /. Howvr, wth W a functon of C, th consttutv quaton 4..13 can b sn to b obctv. h obctv consttutv quatons n ths scton ar ndcatd by a box around thm. 4..3 Elastcty nsors h total dffrntal of th PK strss can b wrttn as S( C) S( C) 1 1 ds : dc : dc C : dc (4..17) whr C s th fourth-ordr tnsor 36

S( C) S( E) C (4..18) and s a masur of th chang n strss whch occurs du to a chang n stran. It s calld th lastcty tnsor (n th matral dscrpton) or th tangnt modulus. Snc S and E ar symmtrc, C posssss th or symmtrs, C kl C kl Clk, and so has 36 ndpndnt componnts. Howvr, f hyprlastc condtons hold, so that S ( E) / /, thn W W ( E) C 4 (4..19) and so C posssss th maor symmtrs, componnts. C C, and only only 1 ndpndnt kl kl h rat form of th abov quatons s C C E S & d ( ) W ( ) C& W ( ) : E& dt : (4..0) or S & 1 C : C& : E& C (4..1) 4..4 A Not on th Stran Enrgy Functon Som confuson can ars n xprssons such as S, S (4..) As tond at th nd of 1.11.5, on can thr us th nrgy functon W * W *( E11, E1, E13, E1, E, E33), a functon of 6 ndpndnt varabls, n whch cas S 11 * 1 * 1 *, S1, S13, K (4..3) 11 1 13 us th nrgy functon W W ( E11, E1, E13, E1, E, E3, E31, E3, E33 ), a functon of 9 varabls, not all of thm ndpndnt, n whch cas S 11, S1, K S 1, K (4..4) 11 1 1 363

and th symmtry of W wth rspct to th strans, whch must b assumd hr, pls that th strss s also symmtrc. 4..5 Hyprlastcty wth Constrants Consdr a Hyprlastc matral whch s subct to th scalar constrant Wthout th constrant, on has ϕ ϕ 0 or & ϕ : F& 0 (4..5) & ( F ) W W : F& P : F& (4..6) h consttutv quaton P / can b drvd from 4..6 whn F & s arbtrary. Howvr, for th matral wth th constrant, th F & cannot b canclld out from ach sd; on has th condtons P : F& φ 0, : F& 0 (4..7) In gnral thn, P α (4..8) whr α s som arbtrary scalar. h strss s thrfor gvn by P and th rat of chang of ntrnal nrgy s + α (4..9) & : F + α : F& W& + αφ& (4..30) h stran nrgy s W + αφ( F) (4..31) h scond trm hr s vdntly zro and so dos not contrbut to th stran nrgy. Howvr, th strss s th drvatv of th stran nrgy wth rspct to a knmatc quantty, and th drvatv of ths last trm, th scond trm n 4..9, s not b zro. h scalar α s, or dtrs th magntud of, a worklss racton strss. It nsurs that th constrant s satsfd; t s not st a valu n th consttutv quaton, rathr, t s 364

dtrd by consdrng partcular problms, through qulbrum and boundary condtons. For N constrants φ 0, 1K N, 4..9 gnralss to P + α (4..3) N 1 Incomprssbl Matrals An ncomprssbl matral s on whos volum rmans constant throughout a moton, and so has th followng constrant: J dt F 1 or J & 0 (4..33) From 1.11.34, (dt A ) / A (dt A) A, on has J 1. hus th PK1 strsss can b wrttn as dj / df JF, whch quals F at P F pf (4..34) and th stran nrgy functon s of th form 1) W p( J (4..35) wth p α. Usng 3.5.9, S F P, on also has S pc (4..36) Usng.5.0 and.5.18b, and th strss powr s ( F EF & ) ( C E& tr ) C : E& 1 J J C C& J & Jtrd Jtr J : (4..37) 1 & & 1 p & S : C : C C : C & & P : F : F pf : F& W& pj& (4..38) consstnt wth 4..30. 365

h stran nrgy of 4..35 s xprssd n trms of F. hs can b r-xprssd n trms of C. Snc J dt F dt C IIIC, th ncomprssblty constrant can b xprssd as and th stran nrgy can b wrttn as III 0 ϕ (4..39) C 1 W p(iiic ) (4..40) h factor of 1/ has bn ncludd so that th p n 4..35 s th sam as th p n 4..40, snc, from 1.11.33, IIIC C (4..41) ladng to th sam xprsson for th PK strss as bfo, Eqn.4..36. Slarly, n trms of th Cauchy strss, on has{ Problm 1} σ F F pi (4..4) Bcaus of ts rol n ths quaton, th scalar p s calld th hydrostatc prssur. Inxtnsbl Constrant Consdr a matral whch s nxtnsbl n a drcton dfnd by a unt vctor  n th rfrnc confguraton. h constrant for ths matral s gvn by Eqn...60, hn, from 4..11, AC ˆ Aˆ 0 φ (4..43) F C FAˆ Aˆ FAˆ Aˆ (4..44) h strss s thn P F + αfaˆ Aˆ (4..45) 366

Post-contractng wth F (wth J 1 ) thn gvs σ F F + αfaˆ FAˆ (4..46) For nxtnsblty n two drctons, on has th strss Aˆ C Aˆ 0, φ Aˆ C Aˆ 0 φ (4..47) 1 1 1 σ ˆ ˆ α F F + α1fa1 FA1 + FA FA ˆ ˆ (4..48) 4..6 Problms 1. Us th quaton P pf to show that th consttutv quaton for an ncomprssbl hyprlastc matral can b wrttn n trms of th Cauchy strss as σ F F pi 367