DYNAMIC FINITE DEFORMATION THERMO-VISCOELASTICITY USING ENERGY-CONSISTENT TIME-INTEGRATION
|
|
- Roger Perkins
- 6 years ago
- Views:
Transcription
1 European Congress on Computational Methods in Applied Sciences and Engineering (ECCOMAS ) J. Eberhardsteiner et.al. (eds.) Vienna, Austria, September -4, DYNAMIC FINITE DEFORMATION THERMO-VISCOELASTICITY USING ENERGY-CONSISTENT TIME-INTEGRATION Melanie Krüger, Michael Groß & Peter Betsch University of Siegen Paul-Bonatz-Straße 9-, 5768 Siegen melanie.krueger@uni-siegen.de, peter.betsch@uni-siegen.de Chemnitz University of Technology Straße der Nationen 6, 9 Chemnitz michael.gross@mb.tu-chemnitz.de Keywords: enhanced GENERIC, enhanced TC algorithm, partitioned discrete gradients Abstract. The main goal of the present work is the description of a dynamic finite deformation thermo-viscoelastic continuum in the enhanced GENERIC (General equations for nonequilibrium reversible irreversible coupling) format. Therefore the integration is done with partitioned discrete derivatives for the thermodynamically consistent system. The system of partial differential equations is described in an enhancement of the so called GENERIC format. This GENERIC format was introduced for thermo-elastodynamic systems. The considered variables of the system are the Poissonian variables, which are the linear momentum, the configuration, the entropy and the internal variable. There are two constitutive equations for the thermo-viscoelastic continuum necessary. The thermal evolution equation is described with Fourier s law of isotropic heat conduction and the viscous part is given by the fourth order compliance tensor. The enhanced numerical stability of the newly developed structure-preserving integrators in comparison to standard integrators is demonstrated by means of numerical examples.
2 INTRODUCTION Structure preserving integrators are meanwhile well-known as integrators, which lead to enhanced robustness and long stability. An often used approach is the discrete derivative of Gonzalez []. Furthermore, the Hamiltonian formulation leads for elastodynamics to an energy-momentum scheme (see Gonzalez []). Öttinger [9] introduced for thermodynamics a GENERIC (General Equations for Non-Equilibrium Reversible Irreversible Coupling) format which includes a dissipative term. The state variables are the Poissonian variables (configuration, linear momenta and entropy). This framework is given for closed systems. In Romero [,, 3] the GENERIC format is applied with partitioned discrete gradients (see Gonzalez []) for a thermoelastic double pendulum and a thermoelastic continuum. This integrator is called the TC (Thermodynamically Consistent) integrator and preserves the underlying structural properties. In Krüger et al.[5] a comparison of the TC integrator and two other structure preserving integrators is presented. In the present work, the GENERIC framework and the TC integrator will be enhanced to thermoviscoelastic systems. For a double pendulum, the enhanced GENERIC format and the enhanced TC algorithm is considered in Krüger et al.[6] THERMOVISCOELASTIC CONTINUUM The motion of a continuum B with a particle P and the boundary B over a t is described with two configurations (see Figure ). Here, the reference configuration is given ϕ X u x B t x B t X e i φ B P B φ t Figure : Configurations of a continuum at t = and the current configuration B t is defined at a t >. and B t are the corresponding boundaries. The particle P is transformed with the nonlinear mappings φ and φ t to the reference configuration and the current configuration B t (see Holzapfel [4]), respectively: X = φ(p,t) x = φ t (P,t) () X andxdenote the position vectors of the pointsx andx, respectively. The nonlinear mapping ϕ maps the position vectorxto the current configurationxby: x = ϕ(x,t) ()
3 The displacement field u is defined by the difference of the position vectors: u(x,t) = x(x,t) X (3) The velocity and acceleration field can be derived by the first and second partial derivative of the mappingϕ: v(x,t) = ϕ(x,t) t =: ϕ a(x,t) = ϕ(x,t) t =: ϕ (4) These derivatives can also be written v = ẋ and a = ẍ. Hence, a dot denotes a partial derivative at a fixed position vector X. The deformation of the continuum is given by the deformation gradient F = ϕ(x,t) = Gradϕ(X,t) (5) X With the definition of the Jacobian determinant J = detf >, the volume element of the reference configuration Ṿ can be related to the current volume ṿ = J Ṿ. The deformation gradient leads to the strain of continuum, which is denoted by the right Cauchy-Green strain tensor C = F T F (6) In order to describe thermal behavior, the entropy is chosen as state variable: s = s(x,t) (7) This variable was also selected in the works of Romero [, 3]. The underlying constitutive law of the internal energy e yields the temperature θ = e s The viscous (history dependend) material behavior is described by an isotropic internal variable (see Reese et al. [] & Groß [3]): (8) C i = F T i F i (9) B z X F i F e ϕ X F x B t Figure : Intermediate configuration of the continuum 3
4 F e describes the elastic deformation gradient. The internal variablec i is symmetric and characterizes the internal structure of the continuum with irreversible (dissipative) effects. Note that the internal deformation gradient F i is given on the intermediate configuration B z (see Figure ), which is a transformation of two inelastic states. In Noll [7, 8] the material isomorphism is introduced, which establishes a connection between two particles. The isotropic material behavior leads to the elastic strain measure 3 PHYSICAL STRUCTURE C e = C C i () The total energy H of the continuum is given by the total kinetic energy T and the total internal energy E: H = T (p)+e(c,s,c i ) () The total kinetic energy T and the total internal energy E are the integrals over the domain of the kinetic energiet and the internal energy e: T (p) = T (p) E(C,s,C i ) = e(c,s,c i ) () Furthermore, the total kinetic energy T is given in terms of the linear momenta p = ρ v and the constant densityρ by: T (p) = p p ρ (3) The total entropy of the systems is defined by the local entropy s functions: S = s (4) As stability criterion, the Lyapunov functionlis introduced: The temperaturθ is the reference temperature. 4 STRONG EVOLUTION EQUATIONS L = H θ S (5) The strong evolution equations are the two equations of motion, which are well known from elastodynamics, the thermal evolution equation and the viscous evolution equation: ϕ = ρ p ṗ = DivP ṡ = [ ] DivQ D int θ Ċ i = C i V : Σ vis The first Piola-Kirchhoff stress tensor P depends on the second Piola-Kirchhoff stress tensor S: P = F S 4 S = e C (6) (7)
5 As constitutive equation for the first Piola-Kirchhoff heat flux Q the Fouriers law of heat conduction is used with the assumption of isotropy: The isotropic heat conduction tensork is defined by Q = K Gradθ (8) K = κj C (9) where κ > denotes the heat conduction parameter and the Jacobian determinantj = detc is evaluated with the right Cauchy-Green strain tensor C. The viscous evolution equation leads to the viscous Mandel stressσ vis, which depends on the inelastic stress tensorγ : Σ vis = C i Γ This inelastic stressγ can be recovered in the internal dissipation Γ = e C i () D int = Γ : Ċ i () The fourth order compliance tensorv is split into a deviatoric and volumetric part: with V = V dev IdevT + V vol n dim I vol () I devt = I T I vol I vol = n dim I I I T = I I The deviatoric and volumetric parameter V dev and V vol as well as the dimension n dim define this compliance tensor with the restrictions 5 ENHANCED GENERIC V dev > V vol > (3) V dev n dim (4) The GENERIC format was introduced by Öttinger [9]. This format is given for thermal and mechanical isolated systems. For the thermoviscoelastic continuum, the GENERIC system will be enhanced. The initial value problem is stated by ż = [ L(z)+L vis (z) ] δh(z)+ [ M(z)+M vis (z) ] δs(z) z(t = ) = z (5) The state vectorz R 3 is defined by the mappingϕ, the linear momentap, the entropy s and the internal variable in Voigt notationc i : ϕ z = p s (6) [C i ] vn 5
6 The matrices L and M denote the skew-symmetric Poisson matrix and the symmetric positivesemidefinite friction matrix, respectively. The viscous evolution equation leads to the enhanced symmetric matricesl vis and M vis. These matrices are as follows: 3 3 I L = I M = 6 DivQ 6 θ L vis = [V C i C i +V C i C i ] vn M vis = 6 θ Dint Here, the parameterv and V are shorthand notations for the expressions (7) V = V dev V = V vol n dim V dev n dim (8) The functional derivatives of the total energy δh and total entropy δs with respect to the state vectorzare DivP (3 ) p δh = ρ δs = (3 ) = const. (9) θ [Γ] vn (6 ) The related degeneracy conditions for the enhanced GENERIC format are similar to the degeneracy conditions of the GENERIC format, which means: [ δh MδS = δs ] L+L vis δh = (3) Furthermore, the following enhanced degeneracy condition is necessary: δh M vis δs = δh L vis δh (3) 5. Structural properties The enhanced GENERIC format is linked to certain structural properties, which are fullfilled neglegting external forces. The first special property is the total energy consistency: Ḣ = δh ż (3) 6
7 Including the skew-symmetry of L and the degeneracy condition of Eq. (3), the rate of the total energy yields Ḣ = δh L vis δh +δh M vis δs (33) = due to the enhanced degeneracy condition Eq. (3). Analogously, the rate of the total entropy Ṡ is given as follows: Ṡ = δs ż B [ = δs ] (34) M+M vis δs Here, the degeneracy condition of Eq. (3) is inserted. The evaluation of the Matrices M and M vis, as well as the vectorδs leads to Ṡ = [κj θ 3 Grad θ θ C Grad θ ] +Dint D tot (35) = θ Generally, the rate of the Lyapunov functionv implies a stable equilibrium state, iff Here, the rate of the Lyapunov function is once more given by: V (36) V = Ḣ θ Ṡ (37) Inserting the two properties of Eqs. (33) and (35), the stable equilibrium state is guaranteed by: D V tot = θ (38) θ 5. Weak evolution equations The weak evolution equations of the enhanced GENERIC format can be derived by the total energy balance of Eq. (3). The functional dervative of the total energy δh yields the test functionsw z. The first and second entry will be replaced by the strong evolution equations: w z = w ϕ w p w s [w Ci ] vn = ṗ ϕ θ [Γ] vn (39) 7
8 The four weak equations of the enhanced GENERIC are w z ż = w ϕ p+ w p DivP ρ B + w Ci : C i V : Γ vis w s θ [ DivQ D int ] (4) 6 DISCRETIZATION IN TIME The weak evolution equations are now discretized in with the enhanced TC (Thermodynamically consistent) integrator. The TC integrator for thermoelastic systems was introduced by Romero [] for a thermoelastic double pendulum. This integrator is now enhanced for the thermoviscoelastic continuum. The intervali = [, T] is split into finite elements of the numbern = [,...,n tp ] with the intervals I n = [t n,t n+ ]. The step size is denoted by h n = t n+ t n. The enhanced TC integrator is based upon the G-equivariant functional derivative of Gonzalez []. The discrete weak evolution equations are given by H n+ H n z n+ z n = w z h n h n B ([ = w z ] L+L vis G H + [ M+M vis] G S ) (4) where 3 3 I L = I , Lvis = ] [V C C +V C C in+ in+ in+ in+ vn M = DivQ 6 θ, Mvis = 6 θ Dint The specific evaluation of the internal variablesc in+, the temperaturesθ is done in the following way: Q and the internal DissipationD int ( ) n+ θ = [ ] ( )n +( ) n+ Q = D s e D int = κ detc n+ C n+ = C in+ Γ (4), the heat flux vector Grad θ : V : C in+ Γ Here, the partitioned discrete gradient operator D( ) is used (see Gonzalez []) for the derivative of the internal energy e with respect to the entropy s. The G-equivariant functional derivatives (43) 8
9 of the total energy and entropy are given by ( ) Div F n+ S G D H = T p p n+ θ [ ] Γ vn3 The second Piola-Kirchhoff stress tensor S with partitioned discrete gradients: S 7 DISCRETIZATION IN SPACE (3 ) (6 ) G (3 ) S = D s s (44) and the inelastic stress tensor Γ are also evaluated = D C e Γ = D Ci e (45) The spatial discretization of the weak evolution equations is performed with the Finite- Element-Method. Therefore, the continuum body B is approximated with n e elements of Ω e B h : B B h = e= Ω e (46) The discrete boundary B h is analogously given by B h = n e e= Ω e. The isoparametric concept leads to ansatz functions, which are used for the geometry as well as for the field variables. These ansatz functionsn A (ξ) are defined on a reference elementω with normalized coordinatesξ. This leads to the following approximated test and trial functions: X e = N A X ea A= w e ϕ = A= x e = N A x ea w e p = A= p e = N A p ea A= s e = N A s ea A= A= ws e = A= N A w ea ϕ N A w ea p N A w ea s The viscous test and trial functions are evaluated on the element level. Hence, it is not necessary to approximate these functions in space. We denote the test function by w e C i and the internal variable by C e i. The deformation gradientf e is approximated as follows: The gradientsj e and J e are defined as J e = Xe ξ (47) F e = j e J e (48) j e = xe ξ The ordinary finite element approximation of the testfunction w e s in Eq. (47) pose a problem. An admissible testfunction is the temperature θ e, but this testfunction is not approximated in 9 (49)
10 Eq. (47) and is not consistent with Eq. (3). Therefore, a projection of the testfunction θ Pe leading to stability (see Romero []) with the node vectorθ pea = w ea s : θ pe = N A θ pea (5) A= based on an additional discrete equation: w θp θ p = w θp θ (5) is necessary. The testfunction of this equation can be approximated like the other testfunctions and is given by: w θ p = s n+ s n h n (5) 8 DISCRETE EVOLUTION EQUATIONS The discrete weak evolution equations are now given for the enhanced GENERIC system. The global system includes the equations of motion, the thermal evolution equation and the projection: e=a,b= e=a,b= w ea ϕ e= A,B= xeb n+ xeb n HAB = h n w ea AB peb n+ p eb n p H = h n w ea AB seb n+ s H seb n = h n e=a,b= w ea θ p HAB θ peb = e= A,B= e= A= e= A= e= A= The local system is given by the viscous evolution equation: e= Ω w e C i : The internal forces F intea F intea = F e n+ Ω C e i n+ C e i n h n = e= and external forces F extea GradN A : S e detj e w ea ϕ MAB ρ w ea p w ea s ( w ea θ p F extea p eb n+ ) F intea (T extea T intea) Ω N A θ e detj e w e C i : Ci e V e : C n+ i e Γ e n+ Ω are given below: F extea = B= H AB r T eb (53) (54) (55) The vectort eb is the node vector of the first Piola-Kirchhoff stress vectort e = A= NA T ea. Analogously, the thermal evolution equation exhibits internal and external parts with the exter-
11 nal heat fluxq e = [[ GradN A T intea = T extea = Ω n af B= Ω A= NA Q ea θ pe : NA θ pe Gradθ pe θ pe NA N B detjr e QeB ] K e Gradθ pe NA D inte θ pe ] detj e (56) The scalars H AB, M AB and H ρ r AB denote the boundary and volume integrals over the ansatz functions with the determinant over the boundariesdetjr e = X e, ξ X e, η : H AB = Ω N A N B detj e 9 NUMERICAL EXAMPLES M AB ρ = ρ H AB H AB r = Ω N A N B detj e r (57) The internal energy e is of Simo-Pister type and can be split into a compressible, a viscous and a thermal part: e(c, s, C i ) = ψ com (C)+ψ vis (C, C i )+e the (C, s) (58) with ψ com = µ (trc n dim ln)+ψ vol (J) ψ vis = µ e (trc e n dim lnj e )+ψ vol (J e ) e the = k [θ(c,s) θ ]+θ n dim β ψvol (J) J µ and λ are the Lamé parameters, µ e and λ e are the viscous Lamé parameters, k is the heat coefficient and β is the coupling parameter. The volumetric free energy ψ vol (J ( ) ) is given by: (59) ψ vol (J ( ) ) = λ ( ) [ln J ( ) +(J ( ) ) ] (6) The parameters for the following examples are given below: λ = 3 µ = 75 ρ = 8.93 λ e = 3 µ e = 75 V dev = V vol = 5 κ = k = 5 θ = 3 β =. (6) The continuum body is a disk (see Figure 3). This disk has an inner radiusr =.8 and an outer radius r =.. The thickness of the disk is t =.4. The initial temperature is logarithmic distributed from the temperature at the outer ring of 38 K to a temperature at the inner ring of 3 K. For the first example we prescribe mechanical Dirichlet boundaries for the inner ring, which means a fixed position, and an initial angular velocity ωz = 6 elsewhere (see Figure 4). In the second example we state mechanical Neumann-Boundaries given by the traction vector T h n+ = B NB t B sin ( π 8 t n+ ) with t B = [ ] T, see Figure 5 and in the
12 r r Figure 3: Continuum disk and the initial temperature contour third example we prescribe thermal Neumann-boundaries with the boundary heat flux Q h = n+ ( ) naf B NB π sin t 8 n+, see Figure 6. These boundary conditions are applied for 8 s. After that the traction vector and the heat flux vanish and the system is totally closed. The step size is t =. s. ω z H S 5 4 V Figure 4: Example with Dirichlet boundaries and initial angular velocity, total energy H, total entropy S and Lyapunov function V In Figure 4 the structural properties of the enhanced GENERIC system are shown for the first example. Here, external forces are neglected, this means that the disk is thermally and mechanically closed. The temperature balances to one temperature for the whole continuum. The total energy H is conserved for the total simulation of 3 s. The total entropy S is increasing and the Lyapunov functionv is decreasing, which indicates a stable system. Figure 5 shows the result for the disk with Neumann boundaries related with the traction
13 T H 3.4 T S V Figure 5: Example with Neumann boundaries for the traction vector, total energy H, total entropy S and Lyapunov function V vector T h. This Neumann boundaries affect the total energy H, as well as the total entropy n+ S and the Lyapunov function V. After 8 s the loads are removed and the structural properties of the enhanced GENERIC format are visible. Q H S 8 6 V Figure 6: Example with Neumann boundaries for the heat flux, total energy H, total entropy S and Lyapunov function V The third example is the disk with Neumann boundaries related to the boundary heat flux Q h. The results are shown in Figure 6. Once again these boundaries affect the structural n+ properties (H, S and V ). After the 8 s of loads, the total energy H is conserved, the total 3
14 entropys increases and the Lyapunov functionv decreases. CONCLUSIONS Structure preserving integrators have become more and more important in the last decades. The preservation of the structural properties of the underlying system leads to long stability and computational robustness. In this work, we have newly proposed the enhanced GENERIC framework in connection with an enhanced TC integrator for thermoviscoelastic systems. This framework is based upon the Poissonian formulation of thermodynamics, which uses the configuration, the linear momenta, the entropy and the internal variables as state variables. The numerical examples show that the structural properties are preserved for closed systems. This means a constant total energy, an increasing total entropy and a decreasing Lyapunov function. REFERENCES [] O. Gonzalez: Design and analysis of conserving integrators for nonlinear Hamiltonian systems with symmetry. Ph.D. Thesis, Department of Mechanical Engineering, Stanford University, March 996. [] O. Gonzalez: Exact Energy and Momentum Conserving Algorithms for General Models in Nonlinear Elasticity. Comput. Methods Appl. Mech. Engrg., 9 (), [3] M. Groß: Higher-order accurate and energy-momentum consistent discretisation of dynamic finite deformation thermo-viscoelasticity. Habilitationsschrift, urn:nbn:de:hbz: , Institut für Numerische Mechanik, Universität Siegen, Oktober 8. [4] G.A. Holzapfel: Nonlinear solid mechanics. A continuum approach for engineering. John Wiley & Sons Ltd.,. [5] M. Krüger, M. Groß and P. Betsch: A comparison of structure-preserving integrators for discrete thermoelastic systems. Comput. Mech., 47 (), 7 7. [6] M. Krüger, M. Groß and P. Betsch: Energy-consistent -integration for dynamic finite deformation thermo-viscoelasticity. Proc. Appl. Math. Mech., (), [7] W. Noll: A mathematical theory of the mechanical behavior of continuous media. Arch. Rational Mech. Anal., (958), [8] W. Noll: A new mathematical theory of simple materials. Arch. Rational Mech. Anal., 48 (97), 5. [9] H.C. Öttinger: Beyond equilibrium thermodynamics. John Wiley & Sons Inc., 5. [] S. Reese and S. Govindjee: A theory of finite viscoelasticity and numerical aspects. Int. J. Solids Structures, 35 (998), [] I. Romero: Thermodynamically consistent -stepping algorithms for non-linear thermomechanical systems. Int. J. Numer. Meth. Engng, 79 (9),
15 [] I. Romero: Algorithms for coupled problems that preserve symmetries and the laws of thermodynamics. Part I: Monolithic integrators and their application to finite strain thermoelasticity. Comput. Methods Appl. Mech. Engrg., (), [3] I. Romero: Algorithms for coupled problems that preserve symmetries and the laws of thermodynamics. Part II: Fractional step methods. Comput. Methods Appl. Mech. Engrg., (),
thermo-viscoelasticity using
thermo-viscoelasticity Chair of Computational Mechanics University of Siegen Chair of Applied Mechanics and Dynamics Chemnitz University of Technology ECCOMAS, 13.9.1 1 Thermoviscoelastic continuum framework
More informationEnergy-consistent time-integration for dynamic finite deformation. thermoviscoelasticity. GAMM Annual Meeting,
for dynamic finite thermo-viscoelasticity Chair of Computational Mechanics University of Siegen GAMM Annual Meeting, 9.04.0 Thermo-viscoelastic double pendulum (closed system) Thermodynamically Consistent
More informationDYNAMIC FINITE DEFORMATION VISCOELASTICITY IN PRINCIPAL STRETCHES: ENERGY-CONSISTENT TIME INTEGRATION USING MIXED FINITE ELEMENTS
COMPDYN 9 ECCOMAS Thematic Conference on Computational Methods in Structural Dynamics and Earthquake Engineering M. Papadrakakis, N.D. Lagaros, M. Fragiadakis (eds.) Rhodes, Greece, -4 June 9 DYNAMIC FINITE
More informationContents of the presentation
Peter Betsch An energy-entropy-consistent finite thermo-viscoelasticity 1 Peter Betsch Chair of Applied Chemnitz University of Technology March 1, 014 1 Chair of Computational Mechanics Institute of Mechanics
More informationOn the consistent discretization in time of nonlinear thermo-elastodynamics
of nonlinear thermo-elastodynamics Chair of Computational Mechanics University of Siegen GAMM Annual Meeting, 3.03.00 0000 0000 Thermodynamic double pendulum Structure preserving integrators Great interest
More information4 Constitutive Theory
ME338A CONTINUUM MECHANICS lecture notes 13 Tuesday, May 13, 2008 4.1 Closure Problem In the preceding chapter, we derived the fundamental balance equations: Balance of Spatial Material Mass ρ t + ρ t
More informationSimulation of Thermomechanical Couplings of Viscoelastic Materials
Simulation of Thermomechanical Couplings of Viscoelastic Materials Frank Neff 1, Thomas Miquel 2, Michael Johlitz 1, Alexander Lion 1 1 Institute of Mechanics Faculty for Aerospace Engineering Universität
More informationA Characterization of Conserved Quantities in Non-Equilibrium Thermodynamics
Entropy 2013, 15, 5580-5596; doi:10.3390/e15125580 OPEN ACCESS entropy ISSN 1099-4300 www.mdpi.com/journal/entropy Article A Characterization of Conserved Quantities in Non-Equilibrium Thermodynamics Ignacio
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 informationDRIVING FORCE IN SIMULATION OF PHASE TRANSITION FRONT PROPAGATION
Chapter 1 DRIVING FORCE IN SIMULATION OF PHASE TRANSITION FRONT PROPAGATION A. Berezovski Institute of Cybernetics at Tallinn Technical University, Centre for Nonlinear Studies, Akadeemia tee 21, 12618
More informationComparison of Models for Finite Plasticity
Comparison of Models for Finite Plasticity A numerical study Patrizio Neff and Christian Wieners California Institute of Technology (Universität Darmstadt) Universität Augsburg (Universität Heidelberg)
More informationContinuum mechanics V. Constitutive equations. 1. Constitutive equation: definition and basic axioms
Continuum mechanics office Math 0.107 ales.janka@unifr.ch http://perso.unifr.ch/ales.janka/mechanics Mars 16, 2011, Université de Fribourg 1. Constitutive equation: definition and basic axioms Constitutive
More informationA hierarchy of higher order and higher grade continua Application to the plasticity and fracture of metallic foams
A hierarchy of higher order and higher grade continua Application to the plasticity and fracture of metallic foams Samuel Forest Centre des Matériaux/UMR 7633 Mines Paris ParisTech /CNRS BP 87, 91003 Evry,
More informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part III Thursday 1 June 2006 1.30 to 4.30 PAPER 76 NONLINEAR CONTINUUM MECHANICS Attempt FOUR questions. There are SIX questions in total. The questions carry equal weight. STATIONERY
More informationOn the Numerical Modelling of Orthotropic Large Strain Elastoplasticity
63 Advances in 63 On the Numerical Modelling of Orthotropic Large Strain Elastoplasticity I. Karsaj, C. Sansour and J. Soric Summary A constitutive model for orthotropic yield function at large strain
More informationThe Effect of Evolving Damage on the Finite Strain Response of Inelastic and Viscoelastic Composites
Materials 2009, 2, 858-894; doi:0.3390/ma204858 Article OPEN ACCESS materials ISSN 996-944 www.mdpi.com/journal/materials The Effect of Evolving Damage on the Finite Strain Response of Inelastic and Viscoelastic
More informationON THE CONSTITUTIVE MODELING OF THERMOPLASTIC PHASE-CHANGE PROBLEMS C. Agelet de Saracibar, M. Cervera & M. Chiumenti ETS Ingenieros de Caminos, Canal
On the Constitutive Modeling of Thermoplastic Phase-change Problems C. AGELET DE SARACIBAR y & M. CERVERA z ETS Ingenieros de Caminos, Canales y Puertos Edificio C1, Campus Norte, UPC, Jordi Girona 1-3,
More informationComputational non-linear structural dynamics and energy-momentum integration schemes
icccbe 2010 Nottingham University Press Proceedings of the International Conference on Computing in Civil and Building Engineering W Tizani (Editor) Computational non-linear structural dynamics and energy-momentum
More informationOther state variables include the temperature, θ, and the entropy, S, which are defined below.
Chapter 3 Thermodynamics In order to complete the formulation we need to express the stress tensor T and the heat-flux vector q in terms of other variables. These expressions are known as constitutive
More informationMIXED FINITE ELEMENT FORMULATIONS FOR THE GALERKIN-BASED TIME INTEGRATION OF FINITE ANISOTROPIC ELASTODYNAMICS
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 18, Glasgow, UK MIXED FINITE ELEMENT FORMULATIONS FOR THE GALERKIN-BASED
More informationStability of Thick Spherical Shells
Continuum Mech. Thermodyn. (1995) 7: 249-258 Stability of Thick Spherical Shells I-Shih Liu 1 Instituto de Matemática, Universidade Federal do Rio de Janeiro Caixa Postal 68530, Rio de Janeiro 21945-970,
More informationA COMPUTATIONAL FRAMEWORK FOR LARGE DEFORMATION PROBLEMS IN FLEXIBLE MULTIBODY DYNAMICS
MULTIBODY DYNAMICS 0, ECCOMAS Thematic Conference J.C. Samin, P. Fisette (eds.) Brussels, Belgium, 4-7 July 0 A COMPUTATIONAL FRAMEWORK FOR LARGE DEFORMATION PROBLEMS IN FLEXIBLE MULTIBODY DYNAMICS Nicolas
More informationCH.9. CONSTITUTIVE EQUATIONS IN FLUIDS. Multimedia Course on Continuum Mechanics
CH.9. CONSTITUTIVE EQUATIONS IN FLUIDS Multimedia Course on Continuum Mechanics Overview Introduction Fluid Mechanics What is a Fluid? Pressure and Pascal s Law Constitutive Equations in Fluids Fluid Models
More informationNonlinear stability of steady flow of Giesekus viscoelastic fluid
Nonlinear stability of steady flow of Giesekus viscoelastic fluid Mark Dostalík, V. Průša, K. Tůma August 9, 2018 Faculty of Mathematics and Physics, Charles University Table of contents 1. Introduction
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 informationSome Aspects of a Discontinuous Galerkin Formulation for Gradient Plasticity at Finite Strains
Some Aspects of a Discontinuous Galerkin Formulation for Gradient Plasticity at Finite Strains Andrew McBride 1 and B. Daya Reddy 1,2 1 Centre for Research in Computational and Applied Mechanics, University
More informationIn this section, thermoelasticity is considered. By definition, the constitutive relations for Gradθ. This general case
Section.. Thermoelasticity In this section, thermoelasticity is considered. By definition, the constitutive relations for F, θ, Gradθ. This general case such a material depend only on the set of field
More informationThis introductory chapter presents some basic concepts of continuum mechanics, symbols and notations for future reference.
Chapter 1 Introduction to Elasticity This introductory chapter presents some basic concepts of continuum mechanics, symbols and notations for future reference. 1.1 Kinematics of finite deformations We
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 informationNon-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3
Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Tommaso Ruggeri Department of Mathematics and Research Center of Applied Mathematics University of Bologna January 21, 2017 ommaso
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 informationVULCAN 2000: A FINITE ELEMENT SYSTEM FOR THE SIMULATION OF CASTING PROCESSES C. Agelet de Saracibar, M. Cervera & M. Chiumenti ETS Ingenieros de Camin
VULCAN 2000: A FINITE ELEMENT SYSTEM FOR THE SIMULATION OF CASTING PROCESSES C. AGELET DE SARACIBAR y & M. CERVERA z ETS Ingenieros de Caminos, Canales y Puertos International Center for Numerical Methods
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 informationContinuum Mechanics Fundamentals
Continuum Mechanics Fundamentals James R. Rice, notes for ES 220, 12 November 2009; corrections 9 December 2009 Conserved Quantities Let a conseved quantity have amount F per unit volume. Examples are
More information02 - introduction to vectors and tensors. me338 - syllabus. introduction tensor calculus. continuum mechanics. introduction. continuum mechanics
02 - introduction to vectors and tensors me338 - syllabus holzapfel nonlinear solid mechanics [2000], chapter 1.1-1.5, pages 1-32 02-1 introduction 2 the potato equations is the branch of mechanics concerned
More informationContinuum Mechanics and Theory of Materials
Peter Haupt Continuum Mechanics and Theory of Materials Translated from German by Joan A. Kurth Second Edition With 91 Figures, Springer Contents Introduction 1 1 Kinematics 7 1. 1 Material Bodies / 7
More informationDynamics of the Mantle and Lithosphere ETH Zürich Continuum Mechanics in Geodynamics: Equation cheat sheet
Dynamics of the Mantle and Lithosphere ETH Zürich Continuum Mechanics in Geodynamics: Equation cheat sheet or all equations you will probably ever need Definitions 1. Coordinate system. x,y,z or x 1,x
More informationINDEX 363. Cartesian coordinates 19,20,42, 67, 83 Cartesian tensors 84, 87, 226
INDEX 363 A Absolute differentiation 120 Absolute scalar field 43 Absolute tensor 45,46,47,48 Acceleration 121, 190, 192 Action integral 198 Addition of systems 6, 51 Addition of tensors 6, 51 Adherence
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 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 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 informationAn Unconditionally Stable, Energy-Momentum Consistent Implementation of the Material-Point Method
Disclaimer The Center for High Performance Computing HPC@UM provides a focus for high performance computing and communication at the University of ew Mexico UM. HPC@UM is committed to innovative research
More informationCVEN 5161 Advanced Mechanics of Materials I
CVEN 5161 Advanced Mechanics of Materials I Instructor: Kaspar J. Willam Revised Version of Class Notes Fall 2003 Chapter 1 Preliminaries The mathematical tools behind stress and strain are housed in Linear
More informationSecond-gradient theory : application to Cahn-Hilliard fluids
Second-gradient theory : application to Cahn-Hilliard fluids P. Seppecher Laboratoire d Analyse Non Linéaire Appliquée Université de Toulon et du Var BP 132-83957 La Garde Cedex seppecher@univ-tln.fr Abstract.
More informationOn Entropy Flux Heat Flux Relation in Thermodynamics with Lagrange Multipliers
Continuum Mech. Thermodyn. (1996) 8: 247 256 On Entropy Flux Heat Flux Relation in Thermodynamics with Lagrange Multipliers I-Shih Liu Instituto de Matemática Universidade do rio de Janeiro, Caixa Postal
More information07 - balance principles. me338 - syllabus balance principles balance principles. cauchy s postulate. cauchy s lemma -t
me338 - syllabus holzapfel nonlinear solid mechanics [2000], chapter 4, pages 131-161 1 2 cauchy s postulate cauchy s lemma cauchy s postulate stress vector t to a plane with normal n at position x only
More informationIranian Journal of Mathematical Sciences and Informatics Vol.2, No.2 (2007), pp 1-16
Iranian Journal of Mathematical Sciences and Informatics Vol.2, No.2 (2007), pp 1-16 THE EFFECT OF PURE SHEAR ON THE REFLECTION OF PLANE WAVES AT THE BOUNDARY OF AN ELASTIC HALF-SPACE W. HUSSAIN DEPARTMENT
More informationA review of Continuum Thermodynamics
A review of Continuum Thermodynamics Amabile Tatone 1 1 Disim, University of L Aquila, Italy November 2017 Summary Thermodynamics of continua is not a simple subject. It deals with the interplay between
More informationRHEOLOGICAL COEFFICIENTS FOR MEDIA WITH MECHANICAL RELAXATION PHENOMENA
Communications to SIMAI Congress, ISSN 187-915, Vol. (7) DOI: 1.1685/CSC6157 RHEOLOGICAL COEFFICIENTS FOR MEDIA WITH MECHANICAL RELAXATION PHENOMENA A. CIANCIO, V. CIANCIO Department of Mathematics, University
More informationComputational Modelling of Isotropic Multiplicative Growth
Copyright c 2005 Tech Science Press CMES, vol.8, no.2, pp.119-134, 2005 Computational Modelling of Isotropic Multiplicative Growth G.Himpel,E.Kuhl,A.Menzel,P.Steinmann 1 Abstract: The changing mass of
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 informationComputational modelling of isotropic multiplicative growth
Copyright c 0000 Tech Science Press CMES, vol.0, no.0, pp.1-14, 0000 Computational modelling of isotropic multiplicative growth G. Himpel, E. Kuhl, A. Menzel, P. Steinmann 1 Abstract: The changing mass
More informationNonlinear Waves in Solid Continua with Finite Deformation
Nonlinear Waves in Solid Continua with Finite Deformation by Jason Knight Submitted to the graduate degree program in Department of Mechanical Engineering and the Graduate Faculty of the University of
More informationInternal Variables and Generalized Continuum Theories
Internal Variables and Generalized Continuum Theories Arkadi Berezovski, Jüri Engelbrecht and Gérard A. Maugin Abstract The canonical thermomechanics on the material manifold is enriched by the introduction
More informationSchur decomposition in the scaled boundary finite element method in elastostatics
IOP Conference Series: Materials Science and Engineering Schur decomposition in the scaled boundary finite element method in elastostatics o cite this article: M Li et al 010 IOP Conf. Ser.: Mater. Sci.
More informationThe Rotating Inhomogeneous Elastic Cylinders of. Variable-Thickness and Density
Applied Mathematics & Information Sciences 23 2008, 237 257 An International Journal c 2008 Dixie W Publishing Corporation, U. S. A. The Rotating Inhomogeneous Elastic Cylinders of Variable-Thickness and
More informationCONSTITUTIVE THEORIES FOR THERMOELASTIC SOLIDS IN LAGRANGIAN DESCRIPTION USING GIBBS POTENTIAL YUSSHY MENDOZA
CONSTITUTIVE THEORIES FOR THERMOELASTIC SOLIDS IN LAGRANGIAN DESCRIPTION USING GIBBS POTENTIAL BY YUSSHY MENDOZA B.S. in Mechanical Engineering, Saint Louis University, 2008 Submitted to the graduate degree
More informationThe Class of Simo & Pister-type Hyperelasticity Relations
The Class of Simo & Pister-type Hyperelasticity Relations Stefan Hartmann Technical Report Series Fac3-0-0 Impressum Publisher: Fakultät für Mathematik/Informatik und Maschinenbau, Technische Universität
More informationConstitutive models: Incremental (Hypoelastic) Stress- Strain relations. and
Constitutive models: Incremental (Hypoelastic) Stress- Strain relations Example 5: an incremental relation based on hyperelasticity strain energy density function and 14.11.2007 1 Constitutive models:
More informationCourse Syllabus: Continuum Mechanics - ME 212A
Course Syllabus: Continuum Mechanics - ME 212A Division Course Number Course Title Academic Semester Physical Science and Engineering Division ME 212A Continuum Mechanics Fall Academic Year 2017/2018 Semester
More informationHamed Babaei a, Valery I. Levitas b. Abstract
Phase-field approach for stress- and temperature-induced phase transformations that satisfies lattice instability conditions. Part 2. Simulations of phase transformations Si I Si II Hamed Babaei a, Valery
More informationThe Finite Element Method II
[ 1 The Finite Element Method II Non-Linear finite element Use of Constitutive Relations Xinghong LIU Phd student 02.11.2007 [ 2 Finite element equilibrium equations: kinematic variables Displacement Strain-displacement
More informationin this web service Cambridge University Press
CONTINUUM MECHANICS This is a modern textbook for courses in continuum mechanics. It provides both the theoretical framework and the numerical methods required to model the behavior of continuous materials.
More information( ) Notes. Fluid mechanics. Inviscid Euler model. Lagrangian viewpoint. " = " x,t,#, #
Notes Assignment 4 due today (when I check email tomorrow morning) Don t be afraid to make assumptions, approximate quantities, In particular, method for computing time step bound (look at max eigenvalue
More information1 Useful Definitions or Concepts
1 Useful Definitions or Concepts 1.1 Elastic constitutive laws One general type of elastic material model is the one called Cauchy elastic material, which depend on only the current local deformation of
More informationA NOTE ON RELATIONSHIP BETWEEN FIXED-POLE AND MOVING-POLE APPROACHES IN STATIC AND DYNAMIC ANALYSIS OF NON-LINEAR SPATIAL BEAM STRUCTURES
European Congress on Computational Methods in Applied Sciences and Engineering (ECCOMAS 212) J. Eberhardsteiner et.al. (eds.) Vienna, Austria, September 1-14, 212 A NOTE ON RELATIONSHIP BETWEEN FIXED-POLE
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 informationFormulation of a nonlinear porosity law for fully saturated porous media at finite strains
Author manuscript, published in "Journal of the Mechanics and Physics of Solids 61, 2 (2013) 537-556" DOI : 10.1016/j.jmps.2012.09.010 Formulation of a nonlinear porosity law for fully saturated porous
More informationME185 Introduction to Continuum Mechanics
Fall, 0 ME85 Introduction to Continuum Mechanics The attached pages contain four previous midterm exams for this course. Each midterm consists of two pages. As you may notice, many of the problems are
More informationAn introduction to implicit constitutive theory to describe the response of bodies
An introduction to implicit constitutive theory to describe the response of bodies Vít Průša prusv@karlin.mff.cuni.cz Mathematical Institute, Charles University in Prague 3 July 2012 Balance laws, Navier
More informationChemnitz Scientific Computing Preprints
Arnd Meyer The Koiter shell equation in a coordinate free description CSC/1-0 Chemnitz Scientific Computing Preprints ISSN 1864-0087 Chemnitz Scientific Computing Preprints Impressum: Chemnitz Scientific
More informationComparison of the multiplicative decompositions F = F θ F M and F = F M F θ in finite strain thermoelasticity
Comparison of the multiplicative decompositions F = F θ F and F = F F θ in finite strain thermoelasticity Stefan Hartmann Technical Report Series Fac3-12-01 Impressum Publisher: Fakultät für athematik/informatik
More information2.073 Notes on finite-deformation theory of isotropic elastic-viscoplastic solids
2.073 Notes on finite-deformation theory of isotropic elastic-viscoplastic solids Lallit Anand Department of Mechanical Engineering Massachusetts Institute of Technology Cambridge, MA 02139, USA April
More informationFETI Methods for the Simulation of Biological Tissues
SpezialForschungsBereich F 32 Karl Franzens Universita t Graz Technische Universita t Graz Medizinische Universita t Graz FETI Methods for the Simulation of Biological Tissues Ch. Augustin O. Steinbach
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 informationVISCOELASTIC PROPERTIES OF FILLED RUBBER. EXPERIMENTAL OBSERVATIONS AND MATERIAL MODELLING
Engineering MECHANICS, Vol. 14, 2007, No. 1/2, p. 81 89 81 VISCOELASTIC PROPERTIES OF FILLED RUBBER. EXPERIMENTAL OBSERVATIONS AND MATERIAL MODELLING Bohdana Marvalova* The paper presents an application
More informationCONSTITUTIVE THEORIES: BASIC PRINCIPLES
1 6.161.6 CONSTITUTIVE THEORIES: BASIC PRINCIPLES I-Shih Liu Instituto de Matemática, Universidade Federal do Rio de Janeiro, Brasil Keywords: Frame of reference, change of observer, Euclidean objectivity,
More informationRole of thermodynamics in modeling the behavior of complex solids
IWNET Summer School 2015 Role of thermodynamics in modeling the behavior of complex solids Bob Svendsen Material Mechanics RWTH Aachen University Microstructure Physics and Alloy Design Max-Planck-Institut
More informationMathematical Model for Thermal Shock Problem of a Generalized Thermoelastic Layered Composite Material with Variable Thermal Conductivity
COMPUTATIONAL METHODS IN SCIENCE AND TECHNOLOGY ( 65-7 (6 Mathematical Model for Thermal Shock Problem of a Generalized Thermoelastic Layered Composite Material with Variable Thermal Conductivity H. M.
More informationConstitutive Relations
Constitutive Relations Andri Andriyana, Ph.D. Centre de Mise en Forme des Matériaux, CEMEF UMR CNRS 7635 École des Mines de Paris, 06904 Sophia Antipolis, France Spring, 2008 Outline Outline 1 Review of
More informationSome improvements of Xfem for cracked domains
Some improvements of Xfem for cracked domains E. Chahine 1, P. Laborde 2, J. Pommier 1, Y. Renard 3 and M. Salaün 4 (1) INSA Toulouse, laboratoire MIP, CNRS UMR 5640, Complexe scientifique de Rangueil,
More informationCourse No: (1 st version: for graduate students) Course Name: Continuum Mechanics Offered by: Chyanbin Hwu
Course No: (1 st version: for graduate students) Course Name: Continuum Mechanics Offered by: Chyanbin Hwu 2011. 11. 25 Contents: 1. Introduction 1.1 Basic Concepts of Continuum Mechanics 1.2 The Need
More informationReceiver. Johana Brokešová Charles University in Prague
Propagation of seismic waves - theoretical background Receiver Johana Brokešová Charles University in Prague Seismic waves = waves in elastic continuum a model of the medium through which the waves propagate
More informationNUMERICAL SIMULATION OF THE INELASTIC SEISMIC RESPONSE OF RC STRUCTURES WITH ENERGY DISSIPATORS
NUMERICAL SIMULATION OF THE INELASTIC SEISMIC RESPONSE OF RC STRUCTURES WITH ENERGY DISSIPATORS ABSTRACT : P Mata1, AH Barbat1, S Oller1, R Boroschek2 1 Technical University of Catalonia, Civil Engineering
More information4.3 Momentum Balance Principles
4.3 Momentum Balance Principles 4.3.1 Balance of linear angular momentum in spatial material description Consider a continuum body B with a set of particles occupying an arbitrary region Ω with boundary
More informationA New Approach in Material Modeling Towards Shear Cutting of Carbon Fiber Reinforced Plastics
DOI: 1.24352/UB.OVGU-218-8 TECHNISCHE MECHANIK, 38, 1, (218), 91 96 submitted: July 31, 217 A New Approach in Material Modeling Towards Shear Cutting of Carbon Fiber Reinforced Plastics L. Poggenpohl,
More informationConfigurational Forces as Basic Concepts of Continuum Physics
Morton E. Gurtin Configurational Forces as Basic Concepts of Continuum Physics Springer Contents 1. Introduction 1 a. Background 1 b. Variational definition of configurational forces 2 с Interfacial energy.
More informationThermodynamics of non-simple elastic materials
Journal of Elasticity, Vol. 6, No.4, October 1976 NoordhotT International Publishing - Leyden Printed in The Netherlands Thermodynamics of non-simple elastic materials R. C. BATRA ME Department, The University
More informationBiomechanics. Soft Tissue Biomechanics
Biomechanics cross-bridges 3-D myocardium ventricles circulation Image Research Machines plc R* off k n k b Ca 2+ 0 R off Ca 2+ * k on R* on g f Ca 2+ R0 on Ca 2+ g Ca 2+ A* 1 A0 1 Ca 2+ Myofilament kinetic
More informationAn object oriented framework: From flexible multibody dynamics to fluid-structure interaction.
The 2 nd Joint International Conference on Multibody System Dynamics May 29 June 1, 2012, Stuttgart, Germany An object oriented framework: From flexible multibody dynamics to fluid-structure interaction.
More informationEnergy and Momentum Conserving Algorithms in Continuum Mechanics. O. Gonzalez
Energy and Momentum Conserving Algorithms in Continuum Mechanics O. Gonzalez Département de Mathématiques École Polytechnique Fédérale de Lausanne Lausanne, Switzerland INTRODUCTION MAIN CONCERN Time integration
More informationTHE BREZIS-EKELAND-NAYROLES VARIATIONAL PRINCIPALE FOR ELASTOPLASTIC PROBLEMS
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 2018, Glasgow, UK THE BREZIS-EKELAND-NAYROLES VARIATIONAL PRINCIPALE
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 informationCONTINUUM MECHANICS. lecture notes 2003 jp dr.-ing. habil. ellen kuhl technical university of kaiserslautern
CONTINUUM MECHANICS lecture notes 2003 jp dr.-ing. habil. ellen kuhl technical university of kaiserslautern Contents Tensor calculus. Tensor algebra.................................... Vector algebra.................................
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 informationFEAP - - A Finite Element Analysis Program
FEAP - - A Finite Element Analysis Program Version 8.5 Theory Manual Robert L. Taylor Department of Civil and Environmental Engineering University of California at Berkeley Berkeley, California 94720-1710
More informationCVEN 7511 Computational Mechanics of Solids and Structures
CVEN 7511 Computational Mechanics of Solids and Structures Instructor: Kaspar J. Willam Original Version of Class Notes Chishen T. Lin Fall 1990 Chapter 1 Fundamentals of Continuum Mechanics Abstract In
More informationHydrodynamics, Thermodynamics, and Mathematics
Hydrodynamics, Thermodynamics, and Mathematics Hans Christian Öttinger Department of Mat..., ETH Zürich, Switzerland Thermodynamic admissibility and mathematical well-posedness 1. structure of equations
More informationSimulating Thin Shells with MPM
Simulating Thin Shells with MPM Biswajit Banerjee Center for the Simulation of Accidental Fires and xplosions University of Utah March 14, 2005 Outline The Problem. Shell Theory. MPM Formulation. Results.
More informationarxiv: v1 [cs.na] 21 Jul 2010
DEVELOPMENT OF THREE DIMENSIONAL CONSTITUTIVE THEORIES BASED ON LOWER DIMENSIONAL EXPERIMENTAL DATA Satish Karra, K. R. Rajagopal, College Station arxiv:1007.3760v1 [cs.na] 1 Jul 010 Abstract. Most three
More information