SOME NEW FINITE ELEMENT METHODS FOR THE ANALYSIS OF SHOCK AND ACCELERATION WAVES IN NONLINEAR MATERIALS. J. T. ODEN The University of Texas at Austin
|
|
- Herbert Lyons
- 5 years ago
- Views:
Transcription
1 ~ SOME NEW FINITE ELEMENT METHODS FOR THE ANALYSIS OF SHOCK AND ACCELERATION WAVES IN NONLINEAR MATERIALS J. T. ODEN The University of Teas at Austin L. C. WELLFORD, Jr. University of Southern California ABSTRACT In this note, we give a concise summary of a number of new discontinuous finite element schemes that we have developed for the analysis of shock and acceleration waves in nonlinearly elastic solids.. INTRODUCTION Recent theoretical and numerical results indicate that many of the traditional shortcomings of Gaerkin methods of shock wave phenomena can be overcome by using finite elements with built-in discontinuities. The theoretical foundations of discontinuous Gaerkin methods for shock waves is developed in [lj, [2J and the development of numerical schemes to implement this theory is described in [3J. The development of an independent theory for discontinuous finite element methods for acceleration waves is described in [4J. In the present paper, we give a summary account of some aspects of this theory, we describe some new techniques for both shock and acceleration waves, and we present some new numerical results. 2. VARIATIONAL PRINCIPLES FOR WAVES We consider a class of one-dimensional problems in dynamic finite elasticity. We use the following notations: X = material particle E. I, I being an open set of particles whose motion is to be studied. u(x,t) = displacement of particle X at time t, t > o. cr(,t) = first Pioa-Kirchhoff stress at X, t p = mass density of a material body B, I c B, in its reference
2 configuration, t = o. ~(I) p Cm(I) o wmo) p = Soboev space of functions with generalized partial derivatives of order ~ m in Lp(I) = space of m-times continuously differentiable functions on I with compact support on I. = closure of Coo(I) in the Soboev norm o "U II p = I \ I aku Ip dx (2.) wn(i) ~ ~ P I k<m [] Y = I(,t+)- I(X,t-)= jump of a function at particle Y. Iy = }(I(X,t+) + I(X,t-)= average value of a function at a jump. f(x,t) = Body force at X,t. u,u,u X = velocity, acceleration, and displacement gradient (at X,t), respectively. It is well-known that a shock wave in a material body B, the displacement gradient u and the velocity u eperience simple jumps. If Y denotes the particle atxwhich these jumps occur, then the following kinematical compatibility condition must hold [u] + V[u ] = 0 y X y (2.2) Here V is the intrinsic velocity of the shock: dy V = dt (2.3) From momentum balances, it can be shown that v2 = [cr]y p[ux]y (2.4) In addition, an acceleration wave represents the propagation of discontinuities in the acceleration and the second gradients, ux~. If an acceleration wave front eists at particle Z, kinematical compatiblity demands that [U]Z - W 2 [U XX ]Z = 0 (2.5) where W is the intrinsic velocity of the acceleration wave, dz W = dt (2.6) and, from momentum considerations, (2.7)
3 Now we assume that at any time t > 0 there may eist in I, N shock waves defined at the collection of surfaces (particles) {y.}n and M acceleration i=l waves at particles {Z.}~ ; of course, it is possible that one or more of the = Zi and Y. coincide. These physical surfaces define natural partitions of I. J into N+ shockless domaind J. and M+ domains devoid of acceleration waves, K. : I/{Yi} N+ =U J i i=l M+l I/{Z.} = U K. i=l (2.8) With these preliminaries, we may state the following variational principles for waves:. Theorem. Let energy be conserved in a locally integrable sense on I. Then linear momentum is weakly conserved in a body with N shock waves Y i if [f [(pu - pfj_ + (Y X N+ N i=l J. i=l ]dx + [ {pvj[u]y; _ IL + p[cr] v y } = {crv} Y. i 0 for everl velocity field v satisfying the kinematical boundary conditions, wherein I = [O,LJ. Moreover, if there are instead M fronts {Z.} at which acceleration waves occur, momentum is conserved if ' [J i= K i M.. \ 3 [(pu - pf)v + crvx]d + ~ {pwi[uxx]z.. = (2.9) (2.0) for every kinematically admissible velocity v Equation (2.9) was derived in [3] and (2.0) was proved in [4J. They suggest that in the presence of both shock and acceleration waves we should have, for every kinematically admissible velocity v, M+N+2 L +L N [(pi.i- pf)v + crv]d {pvi[u]y i k=l J K i= k k M + \ {pw 3 [u ] + [0-] }(v X ) L Zi Z.. =
4 3. FINITE ELEMENT APPROXIMATIONS We now use (2.9) - (2.l) as a basis for discontinuous finite element approimations. We iqtroduce another partition P of I: 0 = X O.< Xl ~ 000 < X G = L, with hi = X' - i- and h = m~ hi. The interval [ -, Xl] = Ii is a finite element. Net, we construct the Ifdiscrete lf spaces, < j < N: < i < G-; Q.. = I. () J. } lj J (3.) < i < M; < j < G-l; 0.. = I. n K.} lj J (3.2) ~,k(i), (~,k(z» = (V(X): V(X) E ~,k(i) (~,k(i», V(Q) = V(L) = O} (3.3) Pk(O) = {p: p is a polynomial of degree ~ 2 on Q} (3.4) These spaces also have finite-element interpolation properties: (3.5) etc.; see, for eample [5J. The semi discrete finite-element method for shock waves consists of seeking UEH~,k(I}nc (I) such that - cry [Z] } i Y. = {cr(o,t}z(o)}i L o (3.6) whereas the semi discrete finite element method for acceleration waves consists of seeking U E ~,k(i)n C (I) such that M+l f;;~. [(plj - pf)v + crvx]dx +L {PW~[UXX]Zi = M i=l [cr(v)]z. }(VX)Z. = {cr(o,t)v(o)}i L o (3.7)
5 Note that in either case the space of trial functions differs from the space of test functions. 4. APPLICATIONS We summarize briefly the results of some recent appication-s. For more details on some of these results, see [lj - [4J; some of these numerical results we describe here are new. We proceed as follows:. We use simple central and forward differences for the temporal approimations in (3.6) and (3.7); e.g. lj(n~t) = JL (U n + l _ 2U n + Un-) ~t U(n~t) = it (Un - U n - ); un = U(X,n~t) 2. For shock wave problems, we take as the constitutive equation aw aw ai a = au- = ar- au a, a =,2,3 X a X W = W(I l, 2, 3 ) = strain energy per unit initial I = (I, II, III) = principal invariants of deformation tensor. a For incompressible materials, we use aw ai am cr=-o~+pala au X du Ci =,2 where p is the hydrostatic pressure. Generally, we assume p can be eliminated via boundary conditions, such as in the case in which B is a thin cylindrical bar; I is its ais, and the stress normal to I is zero. The eamples to be cited below involve Mooney materials, for which wherein C l, C 2 = material constants. W = C l (I-3) + C 2 (II - 3) 3. For acceleration waves, we limit the present discussion to linear Hookean materials for which a = Eu X. 4. We come to the critical question of the construction of the spaces H~,k(I) and ~,k(i). We give two basic eamples: ~.l(i) Finite-Element Space. First consider a piecewise linear local approimation characterized by four parameters: U(X) = J a + bx c + dx < y > y The parameters are determined in terms of the values U, and U 2 of U at the nodal points at each end of the element and in terms of the jump S in U at the shock front Y. We find that
6 2 U(X) = \ U '4J L aa (X) + sa(i<> + Sy</>(X) a=l (4.) wherein the (X) are the usual Lagrange finite element interpolation functions Ci I(X) = - I h ' (4.2) and a(x} and </>(X) are discontinuous shape functions, l, X < Y a(x) = ; 0, X > y X > < y y (4.3) The entire collection is illustrated in Fig.. It is easily verified that these functions satisfy the kinematical compatibility conditions (2.2). '4J.5 a.5 X 2.5 ~ ~ X </>.5 i':j I: y. h ~ X. Figure. Discontinuous linear trial functions An V~,l(I) Finite-Element Space. For acceleration waves, a similar process can be used if we employ piecewise quadratic approimations of the type in Fig. 2. Here we use
7 U - Ia + bx + cx2 d + ex + fx 2 < Y > Y (4.4) together with the equations, a = U d + eh + fh 2 = U 2 d + 2eh + 4fh 2 = U 3 a + by + cy2 = d + ey + fy2 (4.5) Suppose that a shock also eists at Z which has strength S. Then where 3 U(X) = \ U L aa (X) + R8(X) + Q<t>(X) + An(X) (4.6) a=l 3 X X h 2h2 { X. 8(X) = 0, X < Y X > Y 2 3X + ~ ' X ~ Y - 2h 2h <t>(x) = ( 2 3X+L.,X>Y - "2 h 2h2 (4.7) - -+_ 2 I3(X) - - 2h 2h2 n(x) = X2, X < Y 2 0, X > y The generalized coordinates Rand Q in (4.3) are defined as follows: R = S - AZ Q = SZ - AZ 2 2 (4.8) It is easy to prove that this choice of interpoants satisfies the proper kinematical compatibility conditions at the shock front; observe that [UX]Z + (S - AY)[8 ]z + A[nz]z = S [U] = - Sy Z 2 Thus (2.2) holds. Note also that when S = 0, then R = -AZ and Q = -"2AZ. We remark that the position of a shock at n~t is conveniently obtained via the formula
8 n [a ]y psn (4.9) We also remark that special provisions must be made to handle the reflection of waves. When, for eample, a shock wave impacts a fied particle, momentum and kinematical jump conditions must be satisfied. We know that after a reflection, the particle velocity at a wall must be zero. We have, then, at a fied wall Y, p'j2[u X ] - [a] = 0 and V[U ] + 0- = 0. y y y wherein U X ' a, U are known and U and V must be calculated. In actual calculations, we use the Newton-Raphson method to solve this system of nonlinear equations n y h Figure 2. Discontinuous Quadratic Trial Functions With Discontinuity on Left Side of Internal Node 5. SOME NUMERICAL EXAMPLES. As a first eample, we mention a problem described in [3J which concerns the motion of a cylindrical bar of Mooney material,.36 in. long, an initial area of in 2, C = 24 bs/in 2, C 2 =.5 bs./in. 2, p = 0-4 b.-sec 2 /in 4 subjected to a time-varying end force F(t) at its free end which is lbs at t = 0 and is increased linearly to lbs in 2.0 sec. The other end of the bar is fied. 2h
9 A typical result is shown in Fig. 3, where our method is compared with a shock-smearing scheme. Obviously, the discontinuous finite-element scheme provides for a rather sharp definition of the structure of the shock Stress 40 (P.s. I.) , \ \,,,,, It = sec. Shock Fitting Method Parabolic Regularization Method Figure 3. Shock Structure calculated using discontinuous finite elements Particle Velocity.7 (In./Sec ) t =.25 0 sec. Eact Distance from the Free End of the Bar (In.) Solution Figure 4. Velocity profile for compressive acceleration wave.
10 Acceleration wave calculations are more delicate. In Fig. 4 we see the calculated acceleration amplitude in a illeary elastic bar of modulus 00 lbs/in 2, p = 0-4, V = 000 in./sec., h = in., and ~t = sec. The bar is fied at = 0.5, free at = 0, and is subjected to an applied force p(t) at the free end which varies linearly in time from a magnitude of 0.0 at 5 = 0 to 6.0 psi at t = 0 sec., after which it is held constant. Details of the calculation are given in [4J. Acknowledgement. We are pleased to acknowledge support of this work by the Army Research Office in Durham under Contract No. DAHC04-75-G REFERENCES. Wellford, L. C. Jr. and Oden, J. T., "A Theory of Discontinuous Finite Element Galerkin Approimations of Shock Waves in Nonlinear Elastic Solids I. Variational Methods for Nonlinear Waves," Computer Methods in Applied Mechanics and Engineering, Vol. 7, No. (to appear) -2. Wellford, L. C. Jr. and Oden, J. T., "ATheory of Discontinuous Finite Element Gaerkin Approimations of Shock Waves in Nonlinear Elastic Solids II. Accuracy and Convergence, II Computer Methods in App ied Mechani cs and Engineering, Vol. 7, No.2 to appear). 3. Wellford, L. C. Jr. and Oden, J. T" "Discontinuous Finite Element Approimations for the Analysis of Shock Waves in Nonlinearly Elastic Solids," Journal of Computational Physics ( to appear). 4. Oden, J. T. and Wellford, L. C. Jr., "Discontinuous Finite Element Approimations for the Analysis of Acceleration Waves in Elastic Solids," The Mathematics of Finite Elements with A ications, Brunel University, Ubridge, England, April, 975 to be pu ished by Academic Press, London) 5. Oden, J. T. and Reddy, J. N., Mathematical Theory of Finite Elements, Wiley Interscience, New York (to appear).
COMPOSITE PLATE THEORIES
CHAPTER2 COMPOSITE PLATE THEORIES 2.1 GENERAL Analysis of composite plates is usually done based on one of the following the ries. 1. Equivalent single-layer theories a. Classical laminate theory b. Shear
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 informationHIGHER-ORDER THEORIES
HIGHER-ORDER THEORIES THIRD-ORDER SHEAR DEFORMATION PLATE THEORY LAYERWISE LAMINATE THEORY J.N. Reddy 1 Third-Order Shear Deformation Plate Theory Assumed Displacement Field µ u(x y z t) u 0 (x y t) +
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 informationSteps in the Finite Element Method. Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen
Steps in the Finite Element Method Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen General Idea Engineers are interested in evaluating effects such as deformations, stresses,
More informationDevelopment of Truss Equations
CIVL 7/87 Chapter 3 - Truss Equations - Part /53 Chapter 3a Development of Truss Equations Learning Objectives To derive the stiffness matri for a bar element. To illustrate how to solve a bar assemblage
More informationMECh300H Introduction to Finite Element Methods. Finite Element Analysis (F.E.A.) of 1-D Problems
MECh300H Introduction to Finite Element Methods Finite Element Analysis (F.E.A.) of -D Problems Historical Background Hrenikoff, 94 frame work method Courant, 943 piecewise polynomial interpolation Turner,
More informationKINEMATICS OF CONTINUA
KINEMATICS OF CONTINUA Introduction Deformation of a continuum Configurations of a continuum Deformation mapping Descriptions of motion Material time derivative Velocity and acceleration Transformation
More informationBAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS
Journal of Computational and Applied Mechanics, Vol.., No. 1., (2005), pp. 83 94 BAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS Vladimír Kutiš and Justín Murín Department
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 informationInternational Journal of Pure and Applied Mathematics Volume 58 No ,
International Journal of Pure and Applied Mathematics Volume 58 No. 2 2010, 195-208 A NOTE ON THE LINEARIZED FINITE THEORY OF ELASTICITY Maria Luisa Tonon Department of Mathematics University of Turin
More informationLinear DifferentiaL Equation
Linear DifferentiaL Equation Massoud Malek The set F of all complex-valued functions is known to be a vector space of infinite dimension. Solutions to any linear differential equations, form a subspace
More informationLecture 7: The Beam Element Equations.
4.1 Beam Stiffness. A Beam: A long slender structural component generally subjected to transverse loading that produces significant bending effects as opposed to twisting or axial effects. MECH 40: Finite
More informationCIVL 7/8117 Chapter 4 - Development of Beam Equations - Part 2 1/34. Chapter 4b Development of Beam Equations. Learning Objectives
CIV 7/87 Chapter 4 - Development of Beam Equations - Part /4 Chapter 4b Development of Beam Equations earning Objectives To introduce the work-equivalence method for replacing distributed loading by a
More informationThe Generalized Interpolation Material Point Method
Compaction of a foam microstructure The Generalized Interpolation Material Point Method Tungsten Particle Impacting sandstone The Material Point Method (MPM) 1. Lagrangian material points carry all state
More informationIn this section, mathematical description of the motion of fluid elements moving in a flow field is
Jun. 05, 015 Chapter 6. Differential Analysis of Fluid Flow 6.1 Fluid Element Kinematics In this section, mathematical description of the motion of fluid elements moving in a flow field is given. A small
More informationChapter 3 Variational Formulation & the Galerkin Method
Institute of Structural Engineering Page 1 Chapter 3 Variational Formulation & the Galerkin Method Institute of Structural Engineering Page 2 Today s Lecture Contents: Introduction Differential formulation
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 informationInput: A set (x i -yy i ) data. Output: Function value at arbitrary point x. What for x = 1.2?
Applied Numerical Analysis Interpolation Lecturer: Emad Fatemizadeh Interpolation Input: A set (x i -yy i ) data. Output: Function value at arbitrary point x. 0 1 4 1-3 3 9 What for x = 1.? Interpolation
More informationPreliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
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 informationLARGE AMPLITUDE OSCILLATIONS OF A TUBE OF INCOMPRESSIBLE ELASTIC MATERIAL*
71 LARGE AMPLITUDE OSCILLATIONS OF A TUBE OF INCOMPRESSIBLE ELASTIC MATERIAL* BT JAMES K. KNOWLES California Institute of Technology 1. Introduction. In recent years a number of problems have been solved
More information2.1 Strain energy functions for incompressible materials
Chapter 2 Strain energy functions The aims of constitutive theories are to develop mathematical models for representing the real behavior of matter, to determine the material response and in general, to
More information13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)
13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs) A prototypical problem we will discuss in detail is the 1D diffusion equation u t = Du xx < x < l, t > finite-length rod u(x,
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 informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More information12. Stresses and Strains
12. Stresses and Strains Finite Element Method Differential Equation Weak Formulation Approximating Functions Weighted Residuals FEM - Formulation Classification of Problems Scalar Vector 1-D T(x) u(x)
More informationB.Tech. Theory Examination (Semester IV) Engineering Mathematics III
Solved Question Paper 5-6 B.Tech. Theory Eamination (Semester IV) 5-6 Engineering Mathematics III Time : hours] [Maimum Marks : Section-A. Attempt all questions of this section. Each question carry equal
More informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY Dept. of Civil and Environmental Engineering FALL SEMESTER 2014 Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
More informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY SPRING SEMESTER 207 Dept. of Civil and Environmental Engineering Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
More informationGalerkin method for the numerical solution of the RLW equation using quintic B-splines
Journal of Computational and Applied Mathematics 19 (26) 532 547 www.elsevier.com/locate/cam Galerkin method for the numerical solution of the RLW equation using quintic B-splines İdris Dağ a,, Bülent
More informationChapter 4: Fluid Kinematics
Overview Fluid kinematics deals with the motion of fluids without considering the forces and moments which create the motion. Items discussed in this Chapter. Material derivative and its relationship to
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
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 informationChapter 3. Second Order Linear PDEs
Chapter 3. Second Order Linear PDEs 3.1 Introduction The general class of second order linear PDEs are of the form: ax, y)u xx + bx, y)u xy + cx, y)u yy + dx, y)u x + ex, y)u y + f x, y)u = gx, y). 3.1)
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 informationWRT in 2D: Poisson Example
WRT in 2D: Poisson Example Consider 2 u f on [, L x [, L y with u. WRT: For all v X N, find u X N a(v, u) such that v u dv v f dv. Follows from strong form plus integration by parts: ( ) 2 u v + 2 u dx
More informationA Mixed Nonconforming Finite Element for Linear Elasticity
A Mixed Nonconforming Finite Element for Linear Elasticity Zhiqiang Cai, 1 Xiu Ye 2 1 Department of Mathematics, Purdue University, West Lafayette, Indiana 47907-1395 2 Department of Mathematics and Statistics,
More informationNumerical Analysis & Computer Programming
++++++++++ Numerical Analysis & Computer Programming Previous year Questions from 07 to 99 Ramanasri Institute W E B S I T E : M A T H E M A T I C S O P T I O N A L. C O M C O N T A C T : 8 7 5 0 7 0 6
More informationNUMERICAL SOLUTION OF DOPANT DIFFUSION EQUATIONS. M.J. Bainss, C.P. Please and P.K. Sweby. Mathematics Department University of Reading, Berkshire.
271 NUMERICAL SOLUTION OF DOPANT DIFFUSION EQUATIONS M.J. Bainss, C.P. Please and P.K. Sweby. Mathematics Department University of Reading, Berkshire. 272 SUMMARY A description is given of the application
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 informationPRE-LEAVING CERTIFICATE EXAMINATION, 2010
L.7 PRE-LEAVING CERTIFICATE EXAMINATION, 00 MATHEMATICS HIGHER LEVEL PAPER (300 marks) TIME : ½ HOURS Attempt SIX QUESTIONS (50 marks each). WARNING: Marks will be lost if all necessary work is not clearly
More informationarxiv: v1 [math.na] 1 May 2013
arxiv:3050089v [mathna] May 03 Approximation Properties of a Gradient Recovery Operator Using a Biorthogonal System Bishnu P Lamichhane and Adam McNeilly May, 03 Abstract A gradient recovery operator based
More informationLinear Hyperbolic Systems
Linear Hyperbolic Systems Professor Dr E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 8, 2014 1 / 56 We study some basic
More informationChapter 12 Plate Bending Elements. Chapter 12 Plate Bending Elements
CIVL 7/8117 Chapter 12 - Plate Bending Elements 1/34 Chapter 12 Plate Bending Elements Learning Objectives To introduce basic concepts of plate bending. To derive a common plate bending element stiffness
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationWave Propagation Through Soft Tissue Matter
Wave Propagation Through Soft Tissue Matter Marcelo Valdez and Bala Balachandran Center for Energetic Concepts Development Department of Mechanical Engineering University of Maryland College Park, MD 20742-3035
More informationA NOTE ON THE ANALYSIS OF NONLINEAR DYNAMICS OF ELASTIC MEMBRANES BY THE FINITE ELEMENT METHOD. J. T. Oden, J. E. Key, and R. B.
I A NOTE ON THE ANALYSIS OF NONLINEAR DYNAMICS OF ELASTIC MEMBRANES BY THE FINITE ELEMENT METHOD J. T. Oden, J. E. Key, and R. B. Fost Department of Engineering Mechanics The University of Alabama in Huntsville
More informationCH.3. COMPATIBILITY EQUATIONS. Multimedia Course on Continuum Mechanics
CH.3. COMPATIBILITY EQUATIONS Multimedia Course on Continuum Mechanics Overview Introduction Lecture 1 Compatibility Conditions Lecture Compatibility Equations of a Potential Vector Field Lecture 3 Compatibility
More informationTwo hours. To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER. 29 May :45 11:45
Two hours MATH20602 To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER NUMERICAL ANALYSIS 1 29 May 2015 9:45 11:45 Answer THREE of the FOUR questions. If more
More informationProceedings of the ASME th International Conference on Ocean, Offshore and Arctic Engineering OMAE2016 June 19-24, 2016, Busan, South Korea
Proceedings of the ASME 26 35th International Conference on Ocean, Offshore and Arctic Engineering OMAE26 June 9-24, 26, Busan, South Korea OMAE26-54554 LOCAL STRAIN AND STRESS CALCULATION METHODS OF IRREGULAR
More informationBEM for compressible fluid dynamics
BEM for compressible fluid dynamics L. gkerget & N. Samec Faculty of Mechanical Engineering, Institute of Power, Process and Environmental Engineering, University of Maribor, Slovenia Abstract The fully
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 information4 NON-LINEAR ANALYSIS
4 NON-INEAR ANAYSIS arge displacement elasticity theory, principle of virtual work arge displacement FEA with solid, thin slab, and bar models Virtual work density of internal forces revisited 4-1 SOURCES
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 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 informationLesson Rigid Body Dynamics
Lesson 8 Rigid Body Dynamics Lesson 8 Outline Problem definition and motivations Dynamics of rigid bodies The equation of unconstrained motion (ODE) User and time control Demos / tools / libs Rigid Body
More informationChapter 6 2D Elements Plate Elements
Institute of Structural Engineering Page 1 Chapter 6 2D Elements Plate Elements Method of Finite Elements I Institute of Structural Engineering Page 2 Continuum Elements Plane Stress Plane Strain Toda
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 informationIntroduction to Finite Element Method
Introduction to Finite Element Method Dr. Rakesh K Kapania Aerospace and Ocean Engineering Department Virginia Polytechnic Institute and State University, Blacksburg, VA AOE 524, Vehicle Structures Summer,
More informationPractice Exam 1 Solutions
Practice Exam 1 Solutions 1a. Let S be the region bounded by y = x 3, y = 1, and x. Find the area of S. What is the volume of the solid obtained by rotating S about the line y = 1? Area A = Volume 1 1
More informationUNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES
UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES A Thesis by WOORAM KIM Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the
More informationHow might we evaluate this? Suppose that, by some good luck, we knew that. x 2 5. x 2 dx 5
8.4 1 8.4 Partial Fractions Consider the following integral. 13 2x (1) x 2 x 2 dx How might we evaluate this? Suppose that, by some good luck, we knew that 13 2x (2) x 2 x 2 = 3 x 2 5 x + 1 We could then
More informationMITOCW MITRES2_002S10nonlinear_lec15_300k-mp4
MITOCW MITRES2_002S10nonlinear_lec15_300k-mp4 The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources
More 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 informationFirst order Partial Differential equations
First order Partial Differential equations 0.1 Introduction Definition 0.1.1 A Partial Deferential equation is called linear if the dependent variable and all its derivatives have degree one and not multiple
More informationDiscontinuous Galerkin and Finite Difference Methods for the Acoustic Equations with Smooth Coefficients. Mario Bencomo TRIP Review Meeting 2013
About Me Mario Bencomo Currently 2 nd year graduate student in CAAM department at Rice University. B.S. in Physics and Applied Mathematics (Dec. 2010). Undergraduate University: University of Texas at
More informationMa 221 Final Exam Solutions 5/14/13
Ma 221 Final Exam Solutions 5/14/13 1. Solve (a) (8 pts) Solution: The equation is separable. dy dx exy y 1 y0 0 y 1e y dy e x dx y 1e y dy e x dx ye y e y dy e x dx ye y e y e y e x c The last step comes
More informationMathematics 136 Calculus 2 Everything You Need Or Want To Know About Partial Fractions (and maybe more!) October 19 and 21, 2016
Mathematics 36 Calculus 2 Everything You Need Or Want To Know About Partial Fractions (and maybe more!) October 9 and 2, 206 Every rational function (quotient of polynomials) can be written as a polynomial
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 informationMath review. Math review
Math review 1 Math review 3 1 series approximations 3 Taylor s Theorem 3 Binomial approximation 3 sin(x), for x in radians and x close to zero 4 cos(x), for x in radians and x close to zero 5 2 some geometry
More informationThe Non-Linear Field Theories of Mechanics
С. Truesdell-W.Noll The Non-Linear Field Theories of Mechanics Second Edition with 28 Figures Springer-Verlag Berlin Heidelberg NewYork London Paris Tokyo Hong Kong Barcelona Budapest Contents. The Non-Linear
More informationPropagation of discontinuities in solutions of First Order Partial Differential Equations
Propagation of discontinuities in solutions of First Order Partial Differential Equations Phoolan Prasad Department of Mathematics Indian Institute of Science, Bangalore 560 012 E-mail: prasad@math.iisc.ernet.in
More informationSystems of Linear Differential Equations Chapter 7
Systems of Linear Differential Equations Chapter 7 Doreen De Leon Department of Mathematics, California State University, Fresno June 22, 25 Motivating Examples: Applications of Systems of First Order
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationMath 2250 Lab 08 Lab Section: Class ID: Name/uNID: Due Date: 3/23/2017
Math 2250 Lab 08 Lab Section: Class ID: Name/uNID: Due Date: 3/23/2017 TA: Instructions: Unless stated otherwise, please show all your work and explain your reasoning when necessary, as partial credit
More informationMathematical Tripos Part IA Lent Term Example Sheet 1. Calculate its tangent vector dr/du at each point and hence find its total length.
Mathematical Tripos Part IA Lent Term 205 ector Calculus Prof B C Allanach Example Sheet Sketch the curve in the plane given parametrically by r(u) = ( x(u), y(u) ) = ( a cos 3 u, a sin 3 u ) with 0 u
More informationENGN 2290: Plasticity Computational plasticity in Abaqus
ENGN 229: Plasticity Computational plasticity in Abaqus The purpose of these exercises is to build a familiarity with using user-material subroutines (UMATs) in Abaqus/Standard. Abaqus/Standard is a finite-element
More informationAircraft Structures Kirchhoff-Love Plates
University of Liège erospace & Mechanical Engineering ircraft Structures Kirchhoff-Love Plates Ludovic Noels Computational & Multiscale Mechanics of Materials CM3 http://www.ltas-cm3.ulg.ac.be/ Chemin
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More information7.2. Matrix Multiplication. Introduction. Prerequisites. Learning Outcomes
Matrix Multiplication 7.2 Introduction When we wish to multiply matrices together we have to ensure that the operation is possible - and this is not always so. Also, unlike number arithmetic and algebra,
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 13
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras (Refer Slide Time: 00:25) Module - 01 Lecture - 13 In the last class, we have seen how
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 informationApplications in Solid Mechanics
Companies, 4 CHAPTER 9 Applications in Solid 9. INTRODUCTION The bar and beam elements discussed in Chapters 4 are line elements, as only a single coordinate axis is required to define the element reference
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 informationQuintic B-Spline Galerkin Method for Numerical Solutions of the Burgers Equation
5 1 July 24, Antalya, Turkey Dynamical Systems and Applications, Proceedings, pp. 295 39 Quintic B-Spline Galerkin Method for Numerical Solutions of the Burgers Equation İdris Dağ 1,BülentSaka 2 and Ahmet
More informationJim Lambers MAT 419/519 Summer Session Lecture 13 Notes
Jim Lambers MAT 419/519 Summer Session 2011-12 Lecture 13 Notes These notes correspond to Section 4.1 in the text. Least Squares Fit One of the most fundamental problems in science and engineering is data
More informationCIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen
CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall
More informationMarch Algebra 2 Question 1. March Algebra 2 Question 1
March Algebra 2 Question 1 If the statement is always true for the domain, assign that part a 3. If it is sometimes true, assign it a 2. If it is never true, assign it a 1. Your answer for this question
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More informationON DIVISION ALGEBRAS*
ON DIVISION ALGEBRAS* BY J. H. M. WEDDERBURN 1. The object of this paper is to develop some of the simpler properties of division algebras, that is to say, linear associative algebras in which division
More informationBurgers equation - a first look at fluid mechanics and non-linear partial differential equations
Burgers equation - a first look at fluid mechanics and non-linear partial differential equations In this assignment you will solve Burgers equation, which is useo model for example gas dynamics anraffic
More information(0, 0), (1, ), (2, ), (3, ), (4, ), (5, ), (6, ).
1 Interpolation: The method of constructing new data points within the range of a finite set of known data points That is if (x i, y i ), i = 1, N are known, with y i the dependent variable and x i [x
More informationUniversity of Illinois at Urbana-Champaign College of Engineering
University of Illinois at Urbana-Champaign College of Engineering CEE 570 Finite Element Methods (in Solid and Structural Mechanics) Spring Semester 03 Quiz # April 8, 03 Name: SOUTION ID#: PS.: A the
More informationNonlinear Problems of Elasticity
Stuart S. Antman Nonlinear Problems of Elasticity With 105 Illustrations Springer-Verlag New York Berlin Heidelberg London Paris Tokyo Hong Kong Barcelona Budapest Contents Preface vn Chapter I. Background
More informationLaw of behavior very-rubber band: almost incompressible material
Titre : Loi de comportement hyperélastique : matériau pres[...] Date : 25/09/2013 Page : 1/9 Law of behavior very-rubber band: almost incompressible material Summary: One describes here the formulation
More informationNONLINEAR PARTIAL DIFFERENTIAL EQUATIONS IN ENGINEERING AND APPLIED SCIENCE
NONLINEAR PARTIAL DIFFERENTIAL EQUATIONS IN ENGINEERING AND APPLIED SCIENCE Proceedings of a Conference Sponsored By Office of Naval Research Held At University of Rhode Island Kingston, Rhode Island edited
More informationMechanics of Structures (CE130N) Lab 3
UNIVERSITY OF CALIFORNIA AT BERKELEY CE 130N, Spring 2009 Department of Civil and Environmental Engineering Prof. S. Govindjee and Dr. T. Koyama Structural Engineering, Mechanics and Materials Lab 3 1
More information