NONLINEAR CONTINUUM FORMULATIONS CONTENTS
|
|
- Cody Tate
- 5 years ago
- Views:
Transcription
1 NONLINEAR CONTINUUM FORMULATIONS CONTENTS Introduction to nonlinear continuum mechanics Descriptions of motion Measures of stresses and strains Updated and Total Lagrangian formulations Continuum shell element Post-buckling of composite panels
2 DESCRIPTIONS OF MOTION t = t > X 3 x 3 X x X Reference configuration Material description x Current configuration x ( X, t), ( X,) X f= f( X, t) Spatial description X X( x, t) f= f( x, t) JN Reddy Continuum Shell Element
3 MATERIAL TIME DERIVATIVE Material description d d f [ f( X, t)] = [ f ( X, t)] = = X fixed dt dt t Spatial description d dx j [ f( x, t)] = [ f( x, t)] + [ f( x, t)] dt t dt x = [ f( x, t)] + v [ f( x, t)] j t x = [ f( x,)] t + v f( x,) t t j j JN Reddy Continuum Shell Element 3
4 LAGRANGIAN DESCRIPTION OF MOTION Incremental descriptions of deformed configurations Total Lagrangian description Updated Lagrangian description JN Reddy Continuum Shell Element 4
5 Kinematics of Deformation x3, X3 Reference configuration X= ( X,) Deformed configuration Particle X, occupying position x at time t : x = ( X, t ) Deformed configuration x, X Particle X, occupying position x at time t x : = ( X, t ) x, X JN Reddy Continuum Elements 5
6 Kinematics of Deformation (continued) x3 X 3 (time t = ) ( X) T dx,, dx = F dx dx = F dx F dx XQ XP Q dx P x Q u P x P uq x X (time t) _ Q dx _ P x X u x X dxdx-dxdx= dxedx - T EF F I T E= u u u u E ij u u j u u X X X X i m m j i i j T E Green strain tensor JN Reddy Continuum Elements 6
7 Definitions MEASURES OF STRAIN u= x x or u i = x i x i u= x x or u i = x i x i u= u u or u i = u i u i E ij = ( u i E ij = Green Incremental Strain Tensor ε ij d x i d x j =( ds) ( ds) = ( E ij E ij ) d x i d x j ( e ij + η ij )d x i d x j + u j x j + u k ) u k x i x i x j ( u i + u j x j + u k ) u k x i x i x j e ij = ( ui x j + u j + u k u k x i x i + x j Linear Green Incremental Strain Tensor e ij = ( ui x j + u k ) u k x i x j η ij = u k u k x i x j Nonlinear Incremental Strain Tensor Infinitesmal Green-Lagrange Incremental Strain Tensor JN Reddy u ) j x i Continuum Elements 7
8 Point forces Distributed force F F f STRESS VECTOR F Plane 3 F 4 Plane Continuous rigid support Roller support F f F f a nˆ t (nˆ ) nˆ Plane Stress vector t ( nˆ ) f lim a a ( nˆ ) ( nˆ ) t t JN Reddy Continuum Formulations 8
9 CAUCHY STRESS TENSOR x 3 x A tn ˆ B nˆ s t eˆ x 3 ê 3 s t ê ê tn ˆ eˆ nˆ s x t x x t 3 x 3 t x 3 xz 3 zx t σ nˆ or t n ; t eˆ, σ t eˆ eˆ eˆ ( nˆ ) ( nˆ ) i ij j i ji j j j ij i j x 33 zz xx eˆ 3 t 3 s 3 3 zy 3 yz yy yx xy
10 FIRST AND SECOND PIOLA-KIRCHHOFF STRESS TENSORS df TdAPNˆ dapda df tda nˆ da da Nˆ ( X) nˆ df PdAda da da danˆ da, danˆ da Reference configuration Current configuration JN Reddy d J d a F A df df d SdA F f F P A S A d d d d F PS, SJ F σf T P First Piola-Kirchhoff stress tensor S Second Piola-Kirchhoff stress tensor Continuum Formulations
11 CONSERVATION PRICIPLES Conservation of Mass rdv = r dv rjdv = r dv ò ò ò ò o r = rj Conservation of Linear Momentum d r vdv = r fdv + tda ò ò ò dt dv dt T r = s +r f Continuum Formulations
12 Stress Measures (continued) Cauchy Stress Tensor Configuration in which it occurs σ ij = σ ij, σ ij = σ ij Configuration in which it is measured First Piola-Kirchhoff Stress Tensor t= σ ˆn, T= P ˆn d f= td A= T d A ( σ ˆn)d A= ( P ˆn)d A
13 Stress Measures (continued) Second Piola-Kirchhoff Stress Tensor ( ˆn S)d A=d f d f= J d f J = ( x x F, ) x J= det x S = F ( P) S J F s F T F ( x) x, Deformation gradient tensor ( x) X 3
14 DESCRIPTION OF MOTION C C x 3 x 3 x 3 Q C d s P u u Q d s P u Q d s P x x x x x x Motion: x= x(x,x,x 3,t) Notation: coordinates ofa point: x, x, x volum es: V, V, V areas: A, A, A density: ρ, ρ, totaldisplacem ents ofa point: u, ρ u 4
15 Various Strain Components Displacements: u= x x or u i = x i x i u= x x or u i = x i x i u= u u or u i = u i u i Green-Lagrange Strain Tensor ( ds) =d x d x= d x i d x j i ( ds) =d x d x= d x i d x i ( ds) =d x d x= d x i d x i E ij d x i d x j =( ds) ( ds) E ij d x i d x j =( ds) ( ds) d x i = x i x j d x j, d x i = x i x j d x j JN Reddy Continuum Formulations 5
16 Strain Measures (continued) Green-Lagrange Strain Tensor E ij = E ij = ( u i + u j x j + u k ) u k x i x i x j ( u i + u j x j + u k ) u k x i x i x j Green-Lagrange Incremental Strain Tensor Linear Green-Lagrange Incremental Strain Tensor e ij = ε ij d x i d x j =( ds) ( ds) ( ui x j + = ( E ij E ij ) d x i d x j ( e ij + η ij )d x i d x j u j x i + u k x i u k x j + u k ) u k x i x j JN Reddy Continuum Formulations 6
17 Nonlinear Green-Lagrange Incremental StrainTensor η ij = u k x i u k x j Infinitesimal Green-Lagrange Incremental Strain Tensor e ij = ( ui x j + u ) j x i Updated Green-Lagrange Strain (increment) Tensor ( ε ij )d x i d x j = ( ds ) ( ds ) ) x k η ij = x i δ ij x j = ( ui + u j x j + u ) k u k x i x i x j ε ij = ( x k Strain Measures (continued) e ij + η ij e ij = ( uk u k x i x j ( ui + x j ) u ) j x i 7
18 CONSERVATION OF MASS AND CONSTITUTIVE RELATIONS Conservation of Mass ρ= ρ J, ρ= ρ J Constitutive Relations S ij = C ijk l ε kl, S ij = C ijk l ε kl C ijk l = ρ ρ x i x p x j x q x k x r x l x s C pqrs C ijk l = ρ ρ x i x ρ x j x q x k x r x l x s C pqrs JN Reddy Continuum Formulations 8
19 Virtual Work Statements General Virtual Work Statement (in configuration to be determined) δw σ:δ( e) d V δ R= V = σ ij δ( e ij )d V δ R= V δ R= f δu d V+ t δu d S V S = f i δu i d V+ t i δu i d S V Infinitesimal Green-Lagrange incremental strain tensor components e S u u i ij xj j x i Continuum Shell Element 9
20 JN Reddy Virtual Work Statement for the Total Lagrangian Formulation Write various integrals over configuration as equivalent integrals over the reference configuration σ ij δ( e ij )d V= V S ij δ( E ij )d V V f i δu i d V= V f i δu i d V V t i δu i d S= S t i δu i d S S S ij δ( E ij )d V δ( R)= V δ( R)= f i δu i d V+ t i δu i d S V V S S ij δ( E ij )d V δ( R)= Continuum Formulations
21 Total Lagrangian Formulation. Virtual Work Statement: E ij = δ( R)= V f i δu i d V+. Incremental Decompositions: e ij = ( u i x j + u j x i + u k S ij = S t i δu i d S S ij = S ij + S ij S ij = C ijkl E kl, S ij = C ijkl ε kl ( ) ( ui x j + E ij = E ij + ε ij, ε ij = e ij + η ij E ij = V ) u k x i, x j, u j + u k x i x i u k S ij δ( E ij )d V δ( R)= u k x i u k x j + x j ( u i + u j x j + u k u k x i x i x j ( ) ( ) ( ) x i J x m ), η ij = ) σ mn ( x j x n ) u k x i u k x j
22 Total Lagrangian Formulation (continued) 3. Virtual Work Statement with Incremental Decompositions: d( Eij ) d( Eij ) d( eij ) d( eij ), d( eij ) d( eij ) d( hij ) S ij δ( ε ij )d V+ S ij δ( η ij )d V V V =δ( R) S ij δ( e ij )d V V 4. Linearized Virtual Work Statement with Incremental Decompositions: S ij = C ijk l ε kl C ijk l e kl, δ( ε ij ) δ( e ij ) C ijkl e kl δ( e ij )d V+ V S ij δ( η ij )d V V =δ( R) S ij δ( e ij )d V V
23 x =x, x =y, u = u, u = v, u = ū, u = v FINITE ELEMENT MODEL OF THE TOTAL LAGRANGIAN FORMULATION C ijkl e kl δ( e ij )d V= V e ij = ( ui + x j e xx { e} = e yy e xy = ū v ū + v V u j x i + u k x i + u {δ e} T [ C]{ e} d V u k x j + ū + v u k ) u k x i x j u ū + v v u ū + v v v + ū u + v v Total Lagrangian Element 3
24 FINITE ELEMENT MODEL (continued) = ū v ū + v + u ū + v u ū + v v u ū + v v v + ū u + v v [ D] { ū v } u u u + u [ ] D u v v v + v { ū v } e xx { e} = e yy = ([D]+[D u ]){ū} e xy Total Lagrangian Element 4
25 FINITE ELEMENT MODEL (continued) V [D ]= [ D ]= {δ e} = {δ e} T [ C]{ e} d V V V = {δū} T ([D]+ [D u ]) T [ C]([D]+ [D u ]){ū} d V V, [D u ]=, δ e xx δ e yy δ e xy [ C]= u u u + u v C C C C C 66 =([D]+ [D u ]){δū} v v + v Total Lagrangian Element 5
26 V = V = = FINITE ELEMENT MODEL (continued) S ij δ( η ij )d V= {δ η} T { S} d V [ V ( δū ū S xx + δ v ( δū + ū S xy + δ v T V V δū δū δ v δ v V {δ η} T { S} d V ) v + S yy ( δū v + ū δū + v S xx S xy S xy S yy S xx S xy S xy S yy η ij = ū + δ v )] δ v d V ū ū v v ) v d V {δū} T [ D] T [ S][ D]{ū} d V (9.6.) u k x i u k x j Total Lagrangian Element 6
27 FINITE ELEMENT MODEL (continued) [ S]= {u} = {ū} = [Ψ]= S xx S xy S xy S yy S xx S xy S xy S yy { u v} = { n j= u } jψ j (x) n j= v jψ j (x) { n j= ūjψ } j (x) n j= v jψ j (x) =[Ψ]{ } { } ū = =[Ψ]{ } v [ ] ψ ψ ψ n ψ ψ ψ n { } T ={u,v,u,v,,u n,v n } { } T = {ū, v,ū, v,,ū n, v n } Total Lagrangian Element 7
28 FINITE ELEMENT MODEL (continued) V V {δ e} T [ C]{ e} d V = {δū} T ([D]+ [D u ]) T [ C]([D]+ [D u ]){ū} d V V = {δ } T {([D]+ [D u ])[Ψ]} T [ C]([D]+ [D u ])[Ψ]{ }d V V = {δ } T [B L ] T [ C][B L ]{ }d V (9.6.5) V {δ η} T { S} d V = {δū} T [ D ] T [ S][ D ]{ū} d V V = {δ } T [B NL ] T [ S][B NL]{ }d V V Total Lagrangian Element 8
29 FINITE ELEMENT MODEL (continued) δ( R)= = = δ( R)= = V V V V V S ij δ( e ij )d V {δ e} T { S} d V {δ } T [B L ] T { S} d V f i δu i d V+ t i δu i d S S {δ } T [Ψ] T { f} d V+ {δ } T [Ψ] T { t} d ( S Total Lagrangian Element 9
30 FINITE ELEMENT MODEL (continued) [K L ]= [K NL ]= { F}= { f} = [D ]= V V V { f x f y ([K L ]+ [K NL ]){ }= { F} { F} [B L ] T [ C][B L ]d V, { F}= [B L ] T { S} d V V [B NL ] T [ S][B NL ]d V [Ψ] T { f} d V+ [Ψ] T { t} d S }, { t} =, [D u ]= S { t x t y } u u u + u v v v + v Total Lagrangian Element 3
31 FINITE ELEMENT MODEL (continued) [ D ]= [ S]=, [ C]= S xx S xy S xy S yy S xx S xy S xy S yy C C C C C 66, { S}= S xx S yy S xy [B L ]= [B L ]+ [B u L ]+ [B v L ] [B L ]= ψ ψ ψ ψ ψ ψ ψ n ψ ψ ψ n ψ n ψ n Total Lagrangian Element 3
32 FINITE ELEMENT MODEL (continued) [B u L ]= [B v L ]= [B NL ]= u u ψ u ψ u ψ n u ψ u ψ u ψ n ψ + u ψ u ψ + u ψ u ψ n + u ψ n v v v ψ (9.6.7b) ψ v ψ v ψ n ψ v ψ v ψ n + v ψ v ψ + v ψ v ψ n + v ψ ψ ψ ψ ψ ψ ψ n ψ n ψ ψ ψ n ψ n (9.6.7c) (9.6.8) ψ n Total Lagrangian Element 3
33 FINITE ELEMENT MODEL (continued) [ ]{ } [K L ]+ [K N ] [K L ] {ū} [K L ] [K L ]+ [K N ] { v} { ( K ij L =h e C + u ) ( ) ψ i ψ j u Ω + ψ i ψ j C e ( + C + u ) ( u ψi ψ j + ψ ) i ψ j [( + C 66 + u ) ψi + u ] ψ i [( + u ) ψj + u ]} ψ j dx dy { ( K ij L =h e C + u ) ( v ψ i ψ j Ω + C + v ) u ψ i ψ j e [( + C + u )( + v ) ψi ψ j + u ] v ψ i ψ j [( + C 66 + u ) ψi + u ] ψ i [( + v ) ψj + v ]} ψ j dx dy = K ji L 33 = { } { F } { F } { F } { F }
34 FINITE ELEMENT MODEL (continued) K L ij K N ij =h e =h e Ω e Ω e { ( ) v ψ i C + C ( + v + C 66 [( + v [( + v ) ψj + v [ ψ i ψ j S xx + S xy + S yy ψ i ψ j Fi =h e f x ψ i dx dy + h e Ω e Fi =h e f y ψ i dx dy + h e Ω e {( Fi =h e + u ) ψi Ω e + [( + u ( ψ j + C + v ) ( v ψi ψ j + ψ i ) ψi + v ] ψ i ]} ψ j dx dy ( ψi ψ j ] dx dy = K N ij Γ e Γ e t x ψ i ds t y ψ i ds S xx + u ) ψi + u ψ i + ψ i ψ i ] S yy S xy } dx dy ) ψ i ) ψ j ) ψ j ψ j (9.6.3a) 34
35 6. NUMERICAL RESULTS Transverse deflection versus load (cantilever beam under uniformly distributed load).. Without iteration q Load, q O With iteration a = in. b = in. h = in. E = ksi. =. q =.5 b h v y a x, u Transverse deflection, v Total Lagrangian Element 35
36 NUMERICAL RESULTS Load, q Normal stresses vs. load for a cantilever beam under uniform load (updated). Cauchy stress Cauch stress (five Q8 elements) a = in. b = in. h = in. E = ksi. =. q =.5 Second Piola-Kirchhoff stress (five Q8 elements) (updated Lagrangian formulation) Axial stress, sigma-xx (Sxx) 36
37 SHELL FINITE ELEMENTS Degenerate Continuum Shell Elements It is obtained by a degeneration process from the 3D continuum equations. Less expensive than the solid element but more expensive than the shell-theory element. Also exhibit locking problems. Shell Theory Elements Based on a curvilinear description of the continuum using an specific shell theory. Analytical integration of energy terms over the thickness. Also exhibit locking problems. Continuum Shell Element 37
38 CONTINUUM SHELL ELEMENTS The degenerate continuum shell element approach was first developed by Ahmad et al. from a threedimensional solid element by a process which the 3D elasticity equations are expressed in terms of midsurface nodal variables. It uses displacements and rotational DOFs. Continuum Shell Element 38
39 CONTINUUM SHELL ELEMENTS Assumptions are:. Fibers remain straight. The stress normal to the midsurface vanishes (the plane stress condition) Discretization of the reference surface by using Cartesian coordinates where n k k Xi k(, ) ( Xi ) TOP ( Xi ) k BOTTOM (, ): k (,, ): Lagrangian interpolation functions Natural coordinates Continuum Shell Element 39
40 CONTINUUM SHELL ELEMENTS n k k xi ψk(, ) ( xi ) TOP ( xi ) k BOTTOM Discretization of the final configuration of the reference surface Equivalent to the discretization of the displacement of the reference surface as u u h n k k i k (, ) i k e3i k Continuum Shell Element 4
41 Continuum Shell Element (continued) Virtual Work Statement for the Dynamic Case ρ t+ t ü i δ( t+ t u i ) dv + C ijrs l rs δ( l ij ) dv V V + t S ij δ( η ij ) dv = t+ t R t S ij δ( l ij ) dv V Finite Element Model x i = n k= ψ k x k i, t x i = t u i = n k= n k= ψ k t u k i, u i = V ψ k t x k i, t+ t x i = n k= n k= ψ k u k i i=,,3 ψ k t+ t k k, t [M ] t+ t { l }+ ( t [K L ]+ t [K NL ]){ l }= t+ t {R} t {F} Continuum Shell Element 4
42 Numerical Examples Continuum Shell Element 4
43 Numerical Examples Continuum Shell Element 43
44 Numerical Examples 44
45 Numerical Examples Continuum Shell Element 45
46 POSTBUCKLING AND FAILURE ANALYSIS Nonlinear FE analysis Comparison with experimental results of Starnes and Rouse Progressive failure analysis a = 5.8 cm ( in.), b = 7.8 cm (7 in.), h k =.4 mm (.55 in.) 4-ply laminate: [45/-45/ /45/- 45/ /45/-45//9] s E = 3 Gpa (9, ksi), E = 3 Gpa (,89 ksi), u constant w w x u w w x u x P b w a y G = 6.4 Gpa (93 ksi), ν =.38 (graphite-epoxy) Continuum Shell Element46
47 Experimental Setup and Failure Region (Starnes & Rouse, NASA Langley) (a) Typical panel with test fixture (b) A transverse shear failure mode Continuum Shell Element 47
48 POST-BUCKLING OF A COMPOSITE PANEL UNDER IN-PLANE LOADING Finite Element Models Mesh has six elements per buckling mode half wave in each direction 7. in Elements used: () Four-node C -based flat shell element (STAGS) () Nine-node continuum shell element (Chao & Reddy) (3) For-node and nine-node ANS shell elements (Park & Stanley). in All meshes have 7 elements (9 nodes for the meshes with four-node elements and 35 nodes for meshes with nine-node elements). Continuum Shell Element48
49 Comparison of the Experimental (Moire) and Analytical Out-of-Plane Deflection Patterns Experimental Theoretical (FEM) Continuum Shell Element49
50 Comparison of the Experimental and Analytical Solutions for End Shortening Continuum Shell Element 5
51 Postbuckling of a Composite Panel Analysis Test End shortening, u /u cr.5..5 Analysis Test Maximum deflection, w /h Continuum Shell Element 5
52 Postbuckling of a Composite Panel (Stress Distributions) P/P cr =. (MPa) Axial Stress, σ xx (ksi) P/P cr =.5 P/P cr =. P a a y x y/b b -5 Continuum Shell Element 5
53 Postbuckling of a Composite Panel (Stress Distributions) Shear Stress, σ zx (ksi) from constitutive relations from equilibrium relations P/P cr =. P/P cr =.5 P/P cr = y/b (MPa) Continuum Shell Element 53
54 Summary Introduction to nonlinear continuum mechanics Descriptions of motion Measures of stresses and strains Updated and Total Lagrangian formulations Continuum shell element Post-buckling of composite panels Continuum Shell Element 54
The Finite Element Method
The Finite Element Method 3D Problems Heat Transfer and Elasticity Read: Chapter 14 CONTENTS Finite element models of 3-D Heat Transfer Finite element model of 3-D Elasticity Typical 3-D Finite Elements
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More 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 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 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 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 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 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 informationABHELSINKI UNIVERSITY OF TECHNOLOGY
ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for the Morley plate element Jarkko Niiranen Department of Structural
More informationMechanics of materials Lecture 4 Strain and deformation
Mechanics of materials Lecture 4 Strain and deformation Reijo Kouhia Tampere University of Technology Department of Mechanical Engineering and Industrial Design Wednesday 3 rd February, 206 of a continuum
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 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 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 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 informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Overview Milestones in continuum mechanics Concepts of modulus and stiffness. Stress-strain relations Elasticity Surface and body
More informationChapter 5 Elastic Strain, Deflection, and Stability 1. Elastic Stress-Strain Relationship
Chapter 5 Elastic Strain, Deflection, and Stability Elastic Stress-Strain Relationship A stress in the x-direction causes a strain in the x-direction by σ x also causes a strain in the y-direction & z-direction
More informationFINAL EXAMINATION. (CE130-2 Mechanics of Materials)
UNIVERSITY OF CLIFORNI, ERKELEY FLL SEMESTER 001 FINL EXMINTION (CE130- Mechanics of Materials) Problem 1: (15 points) pinned -bar structure is shown in Figure 1. There is an external force, W = 5000N,
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 informationNonlinear Theory of Elasticity. Dr.-Ing. Martin Ruess
Nonlinear Theory of Elasticity Dr.-Ing. Martin Ruess geometry description Cartesian global coordinate system with base vectors of the Euclidian space orthonormal basis origin O point P domain of a deformable
More informationNonlinear bending analysis of laminated composite stiffened plates
Nonlinear bending analysis of laminated composite stiffened plates * S.N.Patel 1) 1) Dept. of Civi Engineering, BITS Pilani, Pilani Campus, Pilani-333031, (Raj), India 1) shuvendu@pilani.bits-pilani.ac.in
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 information2 Introduction to mechanics
21 Motivation Thermodynamic bodies are being characterized by two competing opposite phenomena, energy and entropy which some researchers in thermodynamics would classify as cause and chance or determinancy
More informationStress, Strain, Mohr s Circle
Stress, Strain, Mohr s Circle The fundamental quantities in solid mechanics are stresses and strains. In accordance with the continuum mechanics assumption, the molecular structure of materials is neglected
More informationLinearized theory of elasticity
Linearized theory of elasticity Arie Verhoeven averhoev@win.tue.nl CASA Seminar, May 24, 2006 Seminar: Continuum mechanics 1 Stress and stress principles Bart Nowak March 8 2 Strain and deformation Mark
More 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 informationBasic Equations of Elasticity
A Basic Equations of Elasticity A.1 STRESS The state of stress at any point in a loaded bo is defined completely in terms of the nine components of stress: σ xx,σ yy,σ zz,σ xy,σ yx,σ yz,σ zy,σ zx,andσ
More informationA short review of continuum mechanics
A short review of continuum mechanics Professor Anette M. Karlsson, Department of Mechanical ngineering, UD September, 006 This is a short and arbitrary review of continuum mechanics. Most of this material
More 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 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 informationChapter 3 Stress, Strain, Virtual Power and Conservation Principles
Chapter 3 Stress, Strain, irtual Power and Conservation Principles 1 Introduction Stress and strain are key concepts in the analytical characterization of the mechanical state of a solid body. While stress
More informationVariational principles in mechanics
CHAPTER Variational principles in mechanics.1 Linear Elasticity n D Figure.1: A domain and its boundary = D [. Consider a domain Ω R 3 with its boundary = D [ of normal n (see Figure.1). The problem of
More informationGEOMETRIC NONLINEAR ANALYSIS
GEOMETRIC NONLINEAR ANALYSIS The approach for solving problems with geometric nonlinearity is presented. The ESAComp solution relies on Elmer open-source computational tool [1] for multiphysics problems.
More informationMEC-E8001 FINITE ELEMENT ANALYSIS
MEC-E800 FINIE EEMEN ANAYSIS 07 - WHY FINIE EEMENS AND IS HEORY? Design of machines and structures: Solution to stress or displacement by analytical method is often impossible due to complex geometry,
More informationGeneric Strategies to Implement Material Grading in Finite Element Methods for Isotropic and Anisotropic Materials
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Engineering Mechanics Dissertations & Theses Mechanical & Materials Engineering, Department of Winter 12-9-2011 Generic
More informationIntroduction to Seismology Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 12.510 Introduction to Seismology Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Stress and Strain
More informationMECH 5312 Solid Mechanics II. Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso
MECH 5312 Solid Mechanics II Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso Table of Contents Preliminary Math Concept of Stress Stress Components Equilibrium
More informationBy drawing Mohr s circle, the stress transformation in 2-D can be done graphically. + σ x σ y. cos 2θ + τ xy sin 2θ, (1) sin 2θ + τ xy cos 2θ.
Mohr s Circle By drawing Mohr s circle, the stress transformation in -D can be done graphically. σ = σ x + σ y τ = σ x σ y + σ x σ y cos θ + τ xy sin θ, 1 sin θ + τ xy cos θ. Note that the angle of rotation,
More informationMeasurement of deformation. Measurement of elastic force. Constitutive law. Finite element method
Deformable Bodies Deformation x p(x) Given a rest shape x and its deformed configuration p(x), how large is the internal restoring force f(p)? To answer this question, we need a way to measure deformation
More 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 informationConstitutive Relations
Constitutive Relations Dr. Andri Andriyana Centre de Mise en Forme des Matériaux, CEMEF UMR CNRS 7635 École des Mines de Paris, 06904 Sophia Antipolis, France Spring, 2008 Outline Outline 1 Review of field
More informationConsider an elastic spring as shown in the Fig.2.4. When the spring is slowly
.3 Strain Energy Consider an elastic spring as shown in the Fig..4. When the spring is slowly pulled, it deflects by a small amount u 1. When the load is removed from the spring, it goes back to the original
More informationBACKGROUNDS. Two Models of Deformable Body. Distinct Element Method (DEM)
BACKGROUNDS Two Models of Deformable Body continuum rigid-body spring deformation expressed in terms of field variables assembly of rigid-bodies connected by spring Distinct Element Method (DEM) simple
More informationLecture M1 Slender (one dimensional) Structures Reading: Crandall, Dahl and Lardner 3.1, 7.2
Lecture M1 Slender (one dimensional) Structures Reading: Crandall, Dahl and Lardner 3.1, 7.2 This semester we are going to utilize the principles we learnt last semester (i.e the 3 great principles and
More informationA HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS
A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A. Kroker, W. Becker TU Darmstadt, Department of Mechanical Engineering, Chair of Structural Mechanics Hochschulstr. 1, D-64289 Darmstadt, Germany kroker@mechanik.tu-darmstadt.de,
More 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 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 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 informationCH.2. DEFORMATION AND STRAIN. Continuum Mechanics Course (MMC) - ETSECCPB - UPC
CH.. DEFORMATION AND STRAIN Continuum Mechanics Course (MMC) - ETSECCPB - UPC Overview Introduction Deformation Gradient Tensor Material Deformation Gradient Tensor Inverse (Spatial) Deformation Gradient
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 informationFVM for Fluid-Structure Interaction with Large Structural Displacements
FVM for Fluid-Structure Interaction with Large Structural Displacements Željko Tuković and Hrvoje Jasak Zeljko.Tukovic@fsb.hr, h.jasak@wikki.co.uk Faculty of Mechanical Engineering and Naval Architecture
More informationPREDICTION OF BUCKLING AND POSTBUCKLING BEHAVIOUR OF COMPOSITE SHIP PANELS
FONDATĂ 1976 THE ANNALS OF DUNAREA DE JOS UNIVERSITY OF GALATI. FASCICLE IX. METALLURGY AND MATERIALS SCIENCE N 0. 007, ISSN 15 08X PREDICTION OF BUCKLING AND POSTBUCKLING BEHAVIOUR OF COMPOSITE SHIP PANELS
More informationCONSERVATION OF MASS AND BALANCE OF LINEAR MOMENTUM
CONSERVATION OF MASS AND BALANCE OF LINEAR MOMENTUM Summary of integral theorems Material time derivative Reynolds transport theorem Principle of conservation of mass Principle of balance of linear momentum
More informationApplications of Eigenvalues & Eigenvectors
Applications of Eigenvalues & Eigenvectors Louie L. Yaw Walla Walla University Engineering Department For Linear Algebra Class November 17, 214 Outline 1 The eigenvalue/eigenvector problem 2 Principal
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More 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 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 informationChapter 5 Structural Elements: The truss & beam elements
Institute of Structural Engineering Page 1 Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 2 Chapter Goals Learn how to formulate the Finite Element Equations
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 informationMECH 5312 Solid Mechanics II. Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso
MECH 5312 Solid Mechanics II Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso Table of Contents Thermodynamics Derivation Hooke s Law: Anisotropic Elasticity
More informationLarge Thermal Deflections of a Simple Supported Beam with Temperature-Dependent Physical Properties
Large Thermal Deflections of a Simple Supported Beam with Temperature-Dependent Physical Properties DR. ŞEREF DOĞUŞCAN AKBAŞ Civil Engineer, Şehit Muhtar Mah. Öğüt Sok. No:2/37, 34435 Beyoğlu- Istanbul,
More informationMechanical Properties of Materials
Mechanical Properties of Materials Strains Material Model Stresses Learning objectives Understand the qualitative and quantitative description of mechanical properties of materials. Learn the logic of
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 information3. The linear 3-D elasticity mathematical model
3. The linear 3-D elasticity mathematical model In Chapter we examined some fundamental conditions that should be satisfied in the modeling of all deformable solids and structures. The study of truss structures
More informationConstitutive Equations
Constitutive quations David Roylance Department of Materials Science and ngineering Massachusetts Institute of Technology Cambridge, MA 0239 October 4, 2000 Introduction The modules on kinematics (Module
More informationLecture 3: The Concept of Stress, Generalized Stresses and Equilibrium
Lecture 3: The Concept of Stress, Generalized Stresses and Equilibrium 3.1 Stress Tensor We start with the presentation of simple concepts in one and two dimensions before introducing a general concept
More informationM E 320 Professor John M. Cimbala Lecture 10
M E 320 Professor John M. Cimbala Lecture 10 Today, we will: Finish our example problem rates of motion and deformation of fluid particles Discuss the Reynolds Transport Theorem (RTT) Show how the RTT
More informationStability Analysis of Laminated Composite Thin-Walled Beam Structures
Paper 224 Civil-Comp Press, 2012 Proceedings of the Eleventh International Conference on Computational Structures echnolog, B.H.V. opping, (Editor), Civil-Comp Press, Stirlingshire, Scotland Stabilit nalsis
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 informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationGeometry-dependent MITC method for a 2-node iso-beam element
Structural Engineering and Mechanics, Vol. 9, No. (8) 3-3 Geometry-dependent MITC method for a -node iso-beam element Phill-Seung Lee Samsung Heavy Industries, Seocho, Seoul 37-857, Korea Hyu-Chun Noh
More informationAnalytical Mechanics: Elastic Deformation
Analytical Mechanics: Elastic Deformation Shinichi Hirai Dept. Robotics, Ritsumeikan Univ. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 1 / 60 Agenda Agenda
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 informationAA242B: MECHANICAL VIBRATIONS
AA242B: MECHANICAL VIBRATIONS 1 / 57 AA242B: MECHANICAL VIBRATIONS Dynamics of Continuous Systems These slides are based on the recommended textbook: M. Géradin and D. Rixen, Mechanical Vibrations: Theory
More informationContinuum Mechanics. Continuum Mechanics and Constitutive Equations
Continuum Mechanics Continuum Mechanics and Constitutive Equations Continuum mechanics pertains to the description of mechanical behavior of materials under the assumption that the material is a uniform
More informationA FINITE ELEMENT MODEL FOR THE ANALYSIS OF DELAMINATIONS IN FRP SHELLS
TRENDS IN COMPUTATIONAL STRUCTURAL MECHANICS W.A. Wall, K.-U. Bletzinger and K. Schweizerhof (Eds.) c CIMNE, Barcelona, Spain 2001 A FINITE ELEMENT MODEL FOR THE ANALYSIS OF DELAMINATIONS IN FRP SHELLS
More informationCONSTITUTIVE RELATIONS FOR LINEAR ELASTIC SOLIDS
Chapter 9 CONSTITUTIV RLATIONS FOR LINAR LASTIC SOLIDS Figure 9.1: Hooke memorial window, St. Helen s, Bishopsgate, City of London 211 212 CHAPTR 9. CONSTITUTIV RLATIONS FOR LINAR LASTIC SOLIDS 9.1 Mechanical
More informationPLAXIS. Scientific Manual
PLAXIS Scientific Manual 2016 Build 8122 TABLE OF CONTENTS TABLE OF CONTENTS 1 Introduction 5 2 Deformation theory 7 2.1 Basic equations of continuum deformation 7 2.2 Finite element discretisation 8 2.3
More informationCOMPARISON OF PLATE MODELS FOR ANALYSIS OF LAMINATED COMPOSITES
COMPARISON OF PLATE MODELS FOR ANALYSIS OF LAMINATED COMPOSITES P. M. Mohite and C. S. Upadhyay** Department of Aerospace Engineering, IIT Kanpur 0806, INDIA, e-mail: mohite@iitk.ac.in Assistant Professor,
More informationAccepted Manuscript. R.C. Batra, J. Xiao S (12) Reference: COST Composite Structures. To appear in:
Accepted Manuscript Finite deformations of curved laminated St. Venant-Kirchhoff beam using layerwise third order shear and normal deformable beam theory (TSNDT) R.C. Batra, J. Xiao PII: S0263-8223(12)00486-2
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 informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationNDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16.
CAVITY INSPECTION NDT&E Methods: UT VJ Technologies NDT&E Methods: UT 6. NDT&E: Introduction to Methods 6.1. Ultrasonic Testing: Basics of Elasto-Dynamics 6.2. Principles of Measurement 6.3. The Pulse-Echo
More information(2) 2. (3) 2 Using equation (3), the material time derivative of the Green-Lagrange strain tensor can be obtained as: 1 = + + +
LAGRANGIAN FORMULAION OF CONINUA Review of Continuum Kinematics he reader is referred to Belytscho et al. () for a concise review of the continuum mechanics concepts used here. he notation followed here
More 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 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 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 informationNon-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, Politecnico di Milano, February 17, 2017, Lesson 5
Non-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, attilio.frangi@polimi.it Politecnico di Milano, February 17, 2017, Lesson 5 1 Politecnico di Milano, February 17, 2017, Lesson 5 2 Outline
More information7. Hierarchical modeling examples
7. Hierarchical modeling examples The objective of this chapter is to apply the hierarchical modeling approach discussed in Chapter 1 to three selected problems using the mathematical models studied in
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationFREE VIBRATION OF AXIALLY LOADED FUNCTIONALLY GRADED SANDWICH BEAMS USING REFINED SHEAR DEFORMATION THEORY
FREE VIBRATION OF AXIALLY LOADED FUNCTIONALLY GRADED SANDWICH BEAMS USING REFINED SHEAR DEFORMATION THEORY Thuc P. Vo 1, Adelaja Israel Osofero 1, Marco Corradi 1, Fawad Inam 1 1 Faculty of Engineering
More informationAnalysis of thin plate structures using the absolute nodal coordinate formulation
345 Analysis of thin plate structures using the absolute nodal coordinate formulation K Dufva 1 and A A Shabana 2 1 Department of Mechanical Engineering, Lappeenranta University of echnology, Lappeenranta,
More informationContinuum mechanism: Stress and strain
Continuum mechanics deals with the relation between forces (stress, σ) and deformation (strain, ε), or deformation rate (strain rate, ε). Solid materials, rigid, usually deform elastically, that is the
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 informationDynamics of Ice Sheets and Glaciers
Dynamics of Ice Sheets and Glaciers Ralf Greve Institute of Low Temperature Science Hokkaido University Lecture Notes Sapporo 2004/2005 Literature Ice dynamics Paterson, W. S. B. 1994. The Physics of
More informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems. Prof. Dr. Eleni Chatzi Lecture ST1-19 November, 2015
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Prof. Dr. Eleni Chatzi Lecture ST1-19 November, 2015 Institute of Structural Engineering Method of Finite Elements II 1 Constitutive
More informationSymmetric Bending of Beams
Symmetric Bending of Beams beam is any long structural member on which loads act perpendicular to the longitudinal axis. Learning objectives Understand the theory, its limitations and its applications
More informationReview of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis
uke University epartment of Civil and Environmental Engineering CEE 42L. Matrix Structural Analysis Henri P. Gavin Fall, 22 Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods
More informationChapter 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 informationExternal Work. When a force F undergoes a displacement dx in the same direction i as the force, the work done is
Structure Analysis I Chapter 9 Deflection Energy Method External Work Energy Method When a force F undergoes a displacement dx in the same direction i as the force, the work done is du e = F dx If the
More information