Analytical Mechanics: Elastic Deformation
|
|
- Jason Eaton
- 5 years ago
- Views:
Transcription
1 Analytical Mechanics: Elastic Deformation Shinichi Hirai Dept. Robotics, Ritsumeikan Univ. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation / 59
2 Agenda Agenda Two-dimensional Elastic Deformation Stress-strain Relationship Piecewise Linear Approximation Approximating Potential Energy Dynamic Deformation 2 Three-dimensional Elastic Deformation Stress-strain Relationship Piecewise Linear Approximation Approximating potential energy 3 Summary Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 2 / 59
3 Stress-strain Relationship Displacement vector P(ξ, η) u(ξ, η) v(ξ, η) arbitrary point inside object displacement of point P(ξ, η) along ξ-axis displacement of point P(ξ, η) along η-axis displacement vector u = [ u v Note: vector u depends on ξ and η. ] Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 3 / 59
4 Stress-strain Relationship Cauchy strain B P(ξ,η) A + v η P u η + B u ξ A v ξ Deformation of small square region u v = extension along ξ-axis, = extension along η-axis ξ η u v = shear rotation, = shear + rotation η ξ Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 4 / 59
5 Stress-strain Relationship Cauchy strain normal strain component along ξ-axis at point P: ε ξξ = u ξ normal strain component along η-axis at point P: shear strain at point P: ε ηη = v η 2ε ξη = u η + v ξ Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 5 / 59
6 Stress-strain Relationship Stress-strain relationship pseudo strain vector (strain vector): ε = ε ξξ ε ηη 2ε ξη σ ξξ σ ηη σ ξη normal stress component along ξ-axis at point P normal stress component along η-axis at point P a shear stress component at point P σ = σ ξξ σ ηη σ ξη Note: σ T ε represents energy density (energy per unit volume) Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 6 / 59
7 Stress-strain Relationship Stress-strain relationship in linear elastic material Stress-strain relationship: σ = σ(ε) Stress-strain relationship in linear elastic material: σ = Dε where 3 3 matrix D is referred to as elasticity matrix Elasticity matrix of an linear isotropic material λ + 2µ λ 0 D = λ λ + 2µ µ where λ and µ denote Lamé s constants Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 7 / 59
8 Stress-strain Relationship Stress-strain relationship in linear elastic material Elasticity matrix of an linear isotropic material: D = λi λ + µi µ where I λ = , I µ = Lamé s constants characterize linear isotropic elasticity: λ = νe ( + ν)( 2ν), µ = E 2( + ν) where Young s modulus E and Poisson s ratio ν Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 8 / 59
9 Stress-strain Relationship Potential energy potential energy density of linear elatic material: S 2 σt ε = 2 εt Dε potential energy stored in an elastic object: U = 2 σt ε h ds = 2 εt Dε h ds where h denotes the constant thickness of the object Note: h ds denotes volume of small region S Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 9 / 59
10 Piecewise Linear Approximation Cover of region by triangles region S cover by triangles P i nodal point of a triangle [ ξ i, η i ] T coordinates of point P i P i P j P k a triangle consisting of nodal points P i, P j, P k Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 0 / 59
11 Piecewise Linear Approximation Potential energy U i,j,k : Potential energy stored in P i P j P k U = all triangles U i,j,k How to compute U i,j,k? Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation / 59
12 Piecewise Linear Approximation Signed area Signed area positive if the triangular loop is counter clockwise negative if the loop is clockwise Signed area of triangle OP i P j : OP i P j = 2 ξ i η i ξ j η j = 2 (ξ iη j η i ξ j ) Signed area of triangle P i P j P k : P i P j P k = 2 ξ j ξ i ξ k ξ i η j η i η k η i = 2 {(ξ iη j ξ j η i ) + (ξ j η k ξ k η j ) + (ξ k η i ξ i η k )} Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 2 / 59
13 Piecewise Linear Approximation Shape functions two-dimensional shape function N i,j,k (ξ, η) on triangle P i P j P k : N i,j,k (ξ, η) = (η j η k )ξ (ξ j ξ k )η + (ξ j η k ξ k η j ) 2 P i P j P k Note that { at point Pi N i,j,k (ξ, η) = 0 at point P j and P k Note: P(ξ, η) be an arbitrary point within the triangle: N i,j,k (ξ, η) = PP jp k P i P j P k assume that function N i,j,k (ξ, η) vanishes outside P i P j P k Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 3 / 59
14 Piecewise Linear Approximation Shape functions Any function u(ξ, η) can be linearly approximated inside P i P j P k : u(ξ, η) = u i N i,j,k (ξ, η) + u j N j,k,i (ξ, η) + u k N k,i,j (ξ, η) where u i = u(ξ i, η i ), u j = u(ξ j, η j ), u k = u(ξ k, η k ) Partial derivatives of N i,j,k (ξ, η): N i,j,k ξ = η j η k 2 P i P j P k, N i,j,k η = (ξ j ξ k ) 2 P i P j P k. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 4 / 59
15 Piecewise Linear Approximation Approximating displacement vector displacement vector u in region P i P j P k : u = N i,j,k u i + N j,k,i u j + N k,i,j u k collective vector: u i,j,k = u i u j u k collective vectors: γ u = u i u j u k, γ v = v i v j v k, γ = [ γu γ v ] Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 5 / 59
16 Piecewise Linear Approximation Strain components partial derivatives inside = P i P j P k : u ξ = a T γ u, u η = b T γ u, v ξ = a T γ v, v η = b T γ v, where a = 2 η j η k η k η i η i η j, b = 2 normal and shear strain components inside : ξ j ξ k ξ k ξ i ξ i ξ j ε ξξ = a T γ u, ε ηη = b T γ v 2ε ξη = b T γ u + a T γ v Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 6 / 59
17 Approximating Potential Energy Potential energy density of potential energy of an isotropic linear elastic material: 2 εt Dε = 2 εt (λi λ + µi µ )ε = 2 λ(ε ξξ + ε ηη ) µ { 2ε 2 ξξ + 2ε 2 ηη + (2ε ξη ) 2} potential energy stored in element P i P j P k : where G λ = 2 G µ = 2 U i,j,k = λg λ + µg µ (ε ξξ + ε ηη ) 2 h ds P i P j P k P i P j P k 2 ( ε 2 ξξ + ε 2 ηη ) h ds + 2 P i P j P k (2ε ξη ) 2 h ds Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 7 / 59
18 Approximating Potential Energy Potential energy computing G λ : where L aa = L ab = G λ = 2 γt Lγ = 2 are constant matrices. [ γ T u γ T v ] [ L aa L ab L ba L bb ] [ γu aa T h ds = aa T h L bb = bb T h ds = bb T h ab T h ds = ab T h, L ba = ba T h ds = ba T h γ v ] Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 8 / 59
19 Approximating Potential Energy Potential energy computing G µ : where G µ = 2 γt Mγ = 2 [ γ T u γ T v ] [ M aa M ab M ba M bb ] [ γu γ v ] M aa = 2L aa + L bb, M bb = 2L bb + L aa M ab = L ba, M ba = L ab are constant matrices. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 9 / 59
20 Approximating Potential Energy Permutation matrix permutation converting u i,j,k into γ: u i u j u k v i v j v k = γ = Pu i,j,k u i v i u j v j u k v k P : permutation matrix Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 20 / 59
21 Approximating Potential Energy Potential energy G λ = 2 γt Lγ = 2 (Pu i,j,k) T L(Pu i,j,k ) = 2 ut i,j,kj λ i,j,ku i,j,k G µ = 2 γt Mγ = 2 (Pu i,j,k) T M(Pu i,j,k ) = 2 ut i,j,kj µ i,j,k u i,j,k where Ji,j,k λ = P T LP, J µ i,j,k = P T MP partial connection matrices Potential energy U i,j,k = λg λ + µg µ = 2 ut i,j,kk i,j,k u i,j,k where K i,j,k = λj λ i,j,k + µj µ i,j,k stiffness matrix Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 2 / 59
22 Approximating Potential Energy Nodal elastic forces a set of nodal forces applied to P i, P j, P k : where f i,j,k = K i,j,k u i,j,k K i,j,k = λj λ i,j,k + µj µ i,j,k λ and µ are physical parameters. Ji,j,k λ and J µ i,j,k are geometric; they include no physical parameters. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 22 / 59
23 Approximating Potential Energy Example η P3 P4 P5 T T3 T0 T2 P0 P P2 ξ square divided into 4 triangles: T 0 = P 0 P P 3, T = P P 4 P 3, T 2 = P P 2 P 4, T 3 = P 2 P 5 P 4 thickness h is constantly equal to 2 Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 23 / 59
24 Approximating Potential Energy Example T 0 = P 0 P P 3 : a = [,, 0 ] T, b = [, 0, ] T L = J 0,,3 λ = Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 24 / 59
25 Approximating Potential Energy Example T 0 = P 0 P P 3 : a = [,, 0 ] T, b = [, 0, ] T M = Jµ 0,,3 = Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 25 / 59
26 Approximating Potential Energy Example T = P P 4 P 3 : a = [ 0,, ] T, b = [,, 0 ] T J,4,3 λ = J,4,3 µ = Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 26 / 59
27 Approximating Potential Energy Example T 2 = P P 2 P 4 : a = [,, 0 ] T, b = [, 0, ] T J,2,4 λ = J 0,,3 λ, Jµ,2,4 = Jµ 0,,3 T 3 = P 2 P 5 P 4 : a = [ 0,, ] T, b = [,, 0 ] T J 2,5,4 λ = J,4,3 λ, Jµ 2,5,4 = Jµ,4,3 Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 27 / 59
28 Approximating Potential Energy Example adding the contribution of matrix J 0,,3 λ to connection matrix J λ : J λ = Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 28 / 59
29 Approximating Potential Energy Example adding the contribution of matrix J,4,3 λ to connection matrix J λ : J λ = Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 29 / 59
30 Approximating Potential Energy Example adding the contribution of matrix J,2,4 λ to connection matrix J λ : J λ = Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 30 / 59
31 Approximating Potential Energy Example adding the contribution of matrix J 2,5,4 λ to connection matrix J λ : J λ = Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 3 / 59
32 Approximating Potential Energy Example computing J µ yields: J µ = Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 32 / 59
33 Approximating Potential Energy Example potential energy: u N = u 0 u u 2 u 3 u 4 u 5 U = 2 ut NKu N where K = λj λ + µj µ Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 33 / 59
34 Approximating Potential Energy Example P3 P0 T0 T P4 P T2 P5 T3 p P2 edge P 0 P 3 is fixed on a rigid wall uniform pressure p = [ p ξ, p η ] T is exerted over edge P 2 P 5 Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 34 / 59
35 Approximating Potential Energy Example work done by the pressure W = [ (P2 P 5 h/2) p (P 2 P 5 h/2) p ] T [ u2 u 5 ] = f T extu N where f ext = 0 0 (P 2 P 5 h/2) p 0 0 (P 2 P 5 h/2) p Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 35 / 59
36 Approximating Potential Energy Example u 0 = 0 and u 3 = 0 are integrated into A T u N = 0, where I 2 2 O O O A = O O O I 2 2 = O O O O Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 36 / 59
37 Approximating Potential Energy Example Variational principle in statics under geometric constraints minimize J(u N, λ A ) = 2 u N T Ku N f T extu N λ T A A T u N where λ A is a collective vector consisting of four Lagrange multipliers [ K J u N = Ku N f ext Aλ A = 0, J λ A = A T u N = 0 A A T ] [ un λ A ] = [ fext 0 ] Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 37 / 59
38 Dynamic Deformation Kinetic energy ρ: density of an object at point P(ξ, η) kinetic energy inside = P i P j P k : T i,j,k = P i P j P k 2 ρ ut u h ds total kinetic energy of the object: T = P i P j P k T i,j,k velocity of any point within the triangle: u = u i N i,j,k + u j N j,k,i + u k N k,i,j. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 38 / 59
39 Dynamic Deformation Kinetic energy assume that density ρ is constant: T i,j,k = [ ] u T 2 i u j T u k T Mi,j,k u i u j u k where total kinetic energy: M i,j,k = ρh P ip j P k 2 2I 2 2 I 2 2 I 2 2 I 2 2 2I 2 2 I 2 2 I 2 2 I 2 2 2I 2 2 T = 2 u N T M u N matrix M is referred to as an inertia matrix Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 39 / 59
40 Dynamic Deformation Dynamic equation Lagrangian under geometric constraints L = T U + W + λ T A A T u N = 2 u N T M u N 2 u N T Ku N + f T extu N + λ T A A T u N, where λ A is a set of Lagrange multipliers a set of Lagrange equations of motion: L d L = 0 u N dt u N a set of motion equations of nodal points: Ku N + f ext + Aλ A Mü N = 0 Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 40 / 59
41 Dynamic Deformation Dynamic equation Equation for stabilizing constraints: (A T ü N ) + 2α(A T u N ) + α 2 (A T u N ) = 0 introducing velocity vector v N = u N : u N = v N, M v N Aλ A = Ku N + f ext, A T v N = A T (2αv N + α 2 u N ) u N = v N, [ ] [ ] [ ] M A vn Ku A T = N + f ext A T (2αv N + α 2 u N ) λ A Note: the second linear equation is solvable, implying that we can compute v N. Thus, we can sketch u N and v N using any ODE solver. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 4 / 59
42 Dynamic Deformation Example (E = 30, ν = 0.35 and c = 20, ν vis = 0.35) s 0 s s 30 s Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 42 / 59
43 Dynamic Deformation Question Report #5 due date : Jan. (Wed.) Simulate the deformation of a rectangular viscoelastic object shown in the figure. The bottom surface is fixed to the ground. Uniform pressure is applied to the middle of the top surface downward for a while, then the pressure is released. Use appropriate values of geometrical and physical parameters of the object. p P2 P3 P4 P8 P9 P0 P4 P5 P6 P5 P P7 P0 P P2 P3 Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 43 / 59
44 Three-dimensional Elastic Deformation Stress-strain Relationship 3D Cauchy strain C C B B A P A P Deformation of small cubic region Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 44 / 59
45 Three-dimensional Elastic Deformation Stress-strain Relationship 3D Cauchy strain u v w = ext. along ξ, = ext. along η, = ext. along ζ, ξ η ζ u v = shear in ξη rot ard ζ, = shear in ξη + rot ard ζ, η ξ v w = shear in ηζ rot ard ξ, = shear in ηζ + rot ard ξ, ζ η w u = shear in ζξ rot ard η, = shear in ζξ + rot ard η ξ ζ Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 45 / 59
46 Three-dimensional Elastic Deformation Stress-strain Relationship 3D Cauchy strain normal strain components: ε ξξ = u ξ, shear strain components: 2ε ηζ = v ζ + w η, strain vector: ε ηη = v η, 2ε ζξ = w ξ + u ζ, ε = ε ξξ ε ηη ε ζζ 2ε ηζ 2ε ζξ 2ε ξη ε ζζ = w ζ 2ε ξη = u η + v ξ Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 46 / 59
47 Three-dimensional Elastic Deformation Stress-strain Relationship Stress-strain relationship stress vector: σ = σ ξξ σ ηη σ ζζ σ ηζ σ ζξ σ ξη stress-strain relationship in linear elastic material: σ = Dε Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 47 / 59
48 Three-dimensional Elastic Deformation Stress-strain Relationship Stress-strain relationship elasticity matrix of an isotropic material: D = λi λ + µi µ where [ O I λ = O O ] [ 2I O, I µ = O I ] and =, I = , O = Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 48 / 59
49 Three-dimensional Elastic Deformation Piecewise Linear Approximation Signed volume signed volume of tetrahedron OP i P j P k : OP i P j P k = ξ i ξ j ξ k 6 η i η j η k ζ i ζ j ζ k signed volume of tetrahedron P i P j P k P l : P i P j P k P l = ξ j ξ i ξ k ξ i ξ l ξ i 6 η j η i η k η i η l η i ζ j ζ i ζ k ζ i ζ l ζ i Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 49 / 59
50 Three-dimensional Elastic Deformation Piecewise Linear Approximation Shape functions 3D shape function N i,j,k,l (ξ, η, ζ) on tetrahedron P i P j P k P l : N i,j,k,l (ξ, η, ζ) = PP { jp k P l at point Pi = P i P j P k P l 0 within P j P k P l linear approximation of function u(ξ, η, ζ) inside P i P j P k P l : u(ξ, η, ζ) = u i N i,j,k,l (ξ, η, ζ) + u j N j,k,l,i (ξ, η, ζ) + u k N k,l,i,j (ξ, η, ζ) + u l N l,i,j,k (ξ, η, ζ) Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 50 / 59
51 Three-dimensional Elastic Deformation Piecewise Linear Approximation Shape functions N i,j,k,l ξ N i,j,k,l η where PP j P k P l = = a j,k,l, P i P j P k P l ξ 6 P i P j P k P l = N i,j,k,l b j,k,l, = c j,k,l 6 P i P j P k P l ζ 6 P i P j P k P l a j,k,l = η j ζ j b j,k,l = ζ j ξ j c j,k,l = ξ j η j η k ζ k ζ k ξ k ξ k η k + η k ζ k + ζ k ξ k + ξ k η k η l ζ l ζ l ξ l ξ l η l + η l ζ l + ζ l ξ l + ξ l η l η j ζ j ζ j ξ j ξ j η j Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 5 / 59
52 Three-dimensional Elastic Deformation Approximating potential energy Potential energy potential energy in = P i P j P k P l : U i,j,k,l = P i P j P k P l 2 εt Dε dv collective vectors: and γ u = u i u j u k u l, γ v = u i,j,k,l = v i v j v k v l, γ w u i u j u k u l = w i w j w k w l, γ = γ u γ v γ w Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 52 / 59
53 Three-dimensional Elastic Deformation Approximating potential energy Potential energy partial derivatives: where a = 6 u ξ = a T γ u, u η = b T γ u, u ζ = c T γ u v ξ = a T γ v, v η = b T γ v, v ζ = c T γ v w ξ = a T γ w, w η = b T γ w, w ζ = c T γ w a j,k,l a k,l,i a l,i,j a i,j,k, b = 6 b j,k,l b k,l,i b l,i,j b i,j,k, c = 6 c j,k,l c k,l,i c l,i,j c i,j,k Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 53 / 59
54 Three-dimensional Elastic Deformation Approximating potential energy Potential energy normal strain components: shear strain components: ε ξξ = a T γ u, ε ηη = b T γ v, ε ζζ = c T γ w 2ε ηζ = c T γ v + b T γ w, 2ε ζξ = a T γ w + c T γ u, 2ε ξη = b T γ u + a T γ v density of potential energy of an isotropic linear elastic material: 2 εt Dε = 2 λ(ε ξξ + ε ηη + ε ζζ ) µ { 2ε 2 ξξ + 2ε 2 ηη + 2ε 2 ζζ + (2ε ηζ ) 2 + (2ε ζξ ) 2 + (2ε ξη ) 2} Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 54 / 59
55 Three-dimensional Elastic Deformation Approximating potential energy Potential energy potential energy in = P i P j P k P l : where G λ = 2 G µ = 2 U i,j,k,l = λg λ + µg µ. {ε ξξ + ε ηη + ε ζζ } 2 dv P i P j P k P l { 2ε 2 ξξ + 2ε 2 ηη + 2ε 2 ζζ + (2ε ηζ ) 2 + (2ε ζξ ) 2 + (2ε ξη ) 2} dv Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 55 / 59
56 Three-dimensional Elastic Deformation Approximating potential energy Potential energy G λ = 2 γt Lγ = 2 [ γ T u γ T v γ T w ] L aa L ab L ac L ba L bb L bc L ca L cb L cc γ u γ v γ w L aa = aa T, L bb = bb T, L cc = cc T L ab = L T ba = ab T, L bc = L T cb = bc T, L ca = L T ac = ca T G µ = 2 γt Mγ = [ ] M aa M ab M ac γ u γ T 2 u γv T γw T M ba M bb M bc γ v M ca M cb M cc γ w M aa = 2L aa + L bb + L cc, M bb = 2L bb + L cc + L aa, M cc = 2L cc + L aa + L bb, M ab = L ba, M ba = L ab, M bc = L cb, M cb = L bc, M ca = L ac, M ac = L ca Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 56 / 59
57 Three-dimensional Elastic Deformation Approximating potential energy Permutation matrix permutation converting u i,j,k,l into γ: γ = Pu i,j,k,l u i u j u k u l v i v j v k v l w i w j w k w l = u i v i w i u j v j w j u k v k w k u l v l w l Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 57 / 59
58 Three-dimensional Elastic Deformation Approximating potential energy Potential energy potential energy G λ = 2 ut i,j,kj λ i,j,k u i,j,k, J λ i,j,k = P T LP G µ = 2 ut i,j,kj µ i,j,k u i,j,k, J µ i,j,k = P T MP where U i,j,k,l = 2 ut i,j,k,l K i,j,k,l u i,j,k,l K i,j,k,l = λj λ i,j,k,l + µj µ i,j,k,l Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 58 / 59
59 Summary Summary 2D elastic deformation stress-strain relationship and potential energy piecewise linear approximation of potential energy formulating static deformation based on variational principle formulating dynamic deformation based on variational principle 3D elastic deformation piecewise linear approximation in 3D space Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 59 / 59
Analytical 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 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 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 information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationCIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen
CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall
More informationLecture notes Models of Mechanics
Lecture notes Models of Mechanics Anders Klarbring Division of Mechanics, Linköping University, Sweden Lecture 7: Small deformation theories Klarbring (Mechanics, LiU) Lecture notes Linköping 2012 1 /
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 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 informationUNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich
UNIVERSITY OF SASKATCHEWAN ME 313.3 MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich A CLOSED BOOK EXAMINATION TIME: 3 HOURS For Marker s Use Only LAST NAME (printed): FIRST
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 informationStress analysis of a stepped bar
Stress analysis of a stepped bar Problem Find the stresses induced in the axially loaded stepped bar shown in Figure. The bar has cross-sectional areas of A ) and A ) over the lengths l ) and l ), respectively.
More informationSimple derivation of the parameter formulas for motions and transfers
Einfache Ableitung der Parameterformeln für Bewegungen and Umlegungen Rend. del Circ. mat. di Palermo (1) 9 (1910) 39-339. Simple derivation of the parameter formulas for motions and transfers By A. Schoenflies
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 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 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 informationMidterm Examination. Please initial the statement below to show that you have read it
EN75: Advanced Mechanics of Solids Midterm Examination School of Engineering Brown University NAME: General Instructions No collaboration of any kind is permitted on this examination. You may use two pages
More information16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations
6.2 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive quations Constitutive quations For elastic materials: If the relation is linear: Û σ ij = σ ij (ɛ) = ρ () ɛ ij σ ij =
More 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 informationNONLINEAR STRUCTURAL DYNAMICS USING FE METHODS
NONLINEAR STRUCTURAL DYNAMICS USING FE METHODS Nonlinear Structural Dynamics Using FE Methods emphasizes fundamental mechanics principles and outlines a modern approach to understanding structural dynamics.
More informationProblem " Â F y = 0. ) R A + 2R B + R C = 200 kn ) 2R A + 2R B = 200 kn [using symmetry R A = R C ] ) R A + R B = 100 kn
Problem 0. Three cables are attached as shown. Determine the reactions in the supports. Assume R B as redundant. Also, L AD L CD cos 60 m m. uation of uilibrium: + " Â F y 0 ) R A cos 60 + R B + R C cos
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 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 informationModeling of a soft link continuum formulation
Modeling of a soft link continuum formulation Stanislao Grazioso Thursday 12 th April, 2018 Stanislao Grazioso Geometric Theory of Soft Robots Thursday 12 th April, 2018 1 / 38 Introduction Soft arm: 2D
More informationBasic concepts to start Mechanics of Materials
Basic concepts to start Mechanics of Materials Georges Cailletaud Centre des Matériaux Ecole des Mines de Paris/CNRS Notations Notations (maths) (1/2) A vector v (element of a vectorial space) can be seen
More informationThe 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 informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
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 informationMHA042 - Material mechanics: Duggafrågor
MHA042 - Material mechanics: Duggafrågor 1) For a static uniaxial bar problem at isothermal (Θ const.) conditions, state principle of energy conservation (first law of thermodynamics). On the basis of
More 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 informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More 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 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 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 informationCHAPTER 5. Beam Theory
CHPTER 5. Beam Theory SangJoon Shin School of Mechanical and erospace Engineering Seoul National University ctive eroelasticity and Rotorcraft Lab. 5. The Euler-Bernoulli assumptions One of its dimensions
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 informationElmer :Heat transfert with phase change solid-solid in transient problem Application to silicon properties. SIF file : phasechange solid-solid
Elmer :Heat transfert with phase change solid-solid in transient problem Application to silicon properties 3 6 1. Tb=1750 [K] 2 & 5. q=-10000 [W/m²] 0,1 1 Ω1 4 Ω2 7 3 & 6. α=15 [W/(m²K)] Text=300 [K] 4.
More informationIntroduction to structural dynamics
Introduction to structural dynamics p n m n u n p n-1 p 3... m n-1 m 3... u n-1 u 3 k 1 c 1 u 1 u 2 k 2 m p 1 1 c 2 m2 p 2 k n c n m n u n p n m 2 p 2 u 2 m 1 p 1 u 1 Static vs dynamic analysis Static
More informationNon-linear and time-dependent material models in Mentat & MARC. Tutorial with Background and Exercises
Non-linear and time-dependent material models in Mentat & MARC Tutorial with Background and Exercises Eindhoven University of Technology Department of Mechanical Engineering Piet Schreurs July 7, 2009
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 informationSpectral Element Methods (SEM) May 12, D Wave Equations (Strong and Weak form)
Spectral Element Methods (SEM) May 12, 2009 1 3-D Wave Equations (Strong and Weak form) Strong form (PDEs): ρ 2 t s = T + F in (1) Multiply a global test function W(x) to both sides and integrate over
More informationUniversity of Groningen
University of Groningen Nature-inspired microfluidic propulsion using magnetic actuation Khaderi, S. N.; Baltussen, M. G. H. M.; Anderson, P. D.; Ioan, D.; den Toonder, J.M.J.; Onck, Patrick Published
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 information3.2 Hooke s law anisotropic elasticity Robert Hooke ( ) Most general relationship
3.2 Hooke s law anisotropic elasticity Robert Hooke (1635-1703) Most general relationship σ = C ε + C ε + C ε + C γ + C γ + C γ 11 12 yy 13 zz 14 xy 15 xz 16 yz σ = C ε + C ε + C ε + C γ + C γ + C γ yy
More informationKinetic energy. Analytical Mechanics: Link Mechanisms. Kinetic energy. Agenda. Kinematics of two link open mechanism.
Analtical Mechanics: Link Mechanisms Shinichi Hirai Dept. Robotics, Ritsumeikan Univ. Kinetic energ velocit of the center of mass of link : angular velocit of link : kinetic energ of link : ẋ c l c θ S
More informationChapter 4-b Axially Loaded Members
CIVL 222 STRENGTH OF MATERIALS Chapter 4-b Axially Loaded Members AXIAL LOADED MEMBERS Today s Objectives: Students will be able to: a) Determine the elastic deformation of axially loaded member b) Apply
More informationFigure 1: General Plane Motion (Translation and Rotation)
STRIN ND TH TRNSFORMTION OF STRIN INTRODUCTION - DFORMBL BODY MOTION ) Rigid Bod Motion T T Translation Rotation Figure : General Plane Motion (Translation and Rotation) Figure shows the general plane
More informationUsing MATLAB and. Abaqus. Finite Element Analysis. Introduction to. Amar Khennane. Taylor & Francis Croup. Taylor & Francis Croup,
Introduction to Finite Element Analysis Using MATLAB and Abaqus Amar Khennane Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Croup, an informa business
More information16.20 HANDOUT #2 Fall, 2002 Review of General Elasticity
6.20 HANDOUT #2 Fall, 2002 Review of General Elasticity NOTATION REVIEW (e.g., for strain) Engineering Contracted Engineering Tensor Tensor ε x = ε = ε xx = ε ε y = ε 2 = ε yy = ε 22 ε z = ε 3 = ε zz =
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 informationOptimal thickness of a cylindrical shell under dynamical loading
Optimal thickness of a cylindrical shell under dynamical loading Paul Ziemann Institute of Mathematics and Computer Science, E.-M.-A. University Greifswald, Germany e-mail paul.ziemann@uni-greifswald.de
More informationSymmetry and Properties of Crystals (MSE638) Stress and Strain Tensor
Symmetry and Properties of Crystals (MSE638) Stress and Strain Tensor Somnath Bhowmick Materials Science and Engineering, IIT Kanpur April 6, 2018 Tensile test and Hooke s Law Upto certain strain (0.75),
More informationHomework 7-8 Solutions. Problems
Homework 7-8 Solutions Problems 26 A rhombus is a parallelogram with opposite sides of equal length Let us form a rhombus using vectors v 1 and v 2 as two adjacent sides, with v 1 = v 2 The diagonals of
More informationCH.1. THERMODYNAMIC FOUNDATIONS OF CONSTITUTIVE MODELLING. Computational Solid Mechanics- Xavier Oliver-UPC
CH.1. THERMODYNAMIC FOUNDATIONS OF CONSTITUTIVE MODELLING Computational Solid Mechanics- Xavier Oliver-UPC 1.1 Dissipation approach for constitutive modelling Ch.1. Thermodynamical foundations of constitutive
More informationEnergy Considerations
Physics 42200 Waves & Oscillations Lecture 4 French, Chapter 3 Spring 2016 Semester Matthew Jones Energy Considerations The force in Hooke s law is = Potential energy can be used to describe conservative
More informationME FINITE ELEMENT ANALYSIS FORMULAS
ME 2353 - FINITE ELEMENT ANALYSIS FORMULAS UNIT I FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PROBLEMS 01. Global Equation for Force Vector, {F} = [K] {u} {F} = Global Force Vector [K] = Global Stiffness
More informationNon-linear dynamic modeling using component mode synthesis
University of New Mexico UNM Digital Repository Mechanical Engineering ETDs Engineering ETDs 2-9-2010 Non-linear dynamic modeling using component mode synthesis Jamey Bond Follow this and additional works
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 informationRegent College. Maths Department. Core Mathematics 4. Vectors
Regent College Maths Department Core Mathematics 4 Vectors Page 1 Vectors By the end of this unit you should be able to find: a unit vector in the direction of a. the distance between two points (x 1,
More informationLecture 15: Revisiting bars and beams.
3.10 Potential Energy Approach to Derive Bar Element Equations. We place assumptions on the stresses developed inside the bar. The spatial stress and strain matrices become very sparse. We add (ad-hoc)
More information3. Numerical integration
3. Numerical integration... 3. One-dimensional quadratures... 3. Two- and three-dimensional quadratures... 3.3 Exact Integrals for Straight Sided Triangles... 5 3.4 Reduced and Selected Integration...
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 informationNONLINEAR WAVE EQUATIONS ARISING IN MODELING OF SOME STRAIN-HARDENING STRUCTURES
NONLINEAR WAE EQUATIONS ARISING IN MODELING OF SOME STRAIN-HARDENING STRUCTURES DONGMING WEI Department of Mathematics, University of New Orleans, 2 Lakeshore Dr., New Orleans, LA 7148,USA E-mail: dwei@uno.edu
More informationNONLINEAR CONTINUUM FORMULATIONS CONTENTS
NONLINEAR CONTINUUM FORMULATIONS CONTENTS Introduction to nonlinear continuum mechanics Descriptions of motion Measures of stresses and strains Updated and Total Lagrangian formulations Continuum shell
More informationName (Print) ME Mechanics of Materials Exam # 1 Date: October 5, 2016 Time: 8:00 10:00 PM
Name (Print) (Last) (First) Instructions: ME 323 - Mechanics of Materials Exam # 1 Date: October 5, 2016 Time: 8:00 10:00 PM Circle your lecturer s name and your class meeting time. Gonzalez Krousgrill
More informationOBJECTIVES & ASSUMPTIONS
OBJECTIVES & ASSUMPTIONS Z Project Objectives: Construct a template for the 5-node pyramid element Conduct higher-order patch tests Verify element formulation Determine the optimal element in bending X
More informationChapter 4 Analysis of a cantilever
Chapter 4 Analysis of a cantilever Before a complex structure is studied performing a seismic analysis, the behaviour of simpler ones should be fully understood. To achieve this knowledge we will start
More informationGeneral elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
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 informationFurther Linear Elasticity
Torsion of cylindrical bodies Further Linear Elasticity Problem Sheet # 1. Consider a cylindrical body of length L, the ends of which are subjected to distributions of tractions that are statically equivalent
More informationMMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM
B Course Content: A INTRODUCTION AND OVERVIEW Numerical method and Computer-Aided Engineering; Phsical problems; Mathematical models; Finite element method;. B Elements and nodes, natural coordinates,
More informationTHE USE OF DYNAMIC RELAXATION TO SOLVE THE DIFFERENTIAL EQUATION DESCRIBING THE SHAPE OF THE TALLEST POSSIBLE BUILDING
VII International Conference on Textile Composites and Inflatable Structures STRUCTURAL MEMBRANES 2015 E. Oñate, K.-U.Bletzinger and B. Kröplin (Eds) THE USE OF DYNAMIC RELAXATION TO SOLVE THE DIFFERENTIAL
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 informationINTRODUCTION TO STRAIN
SIMPLE STRAIN INTRODUCTION TO STRAIN In general terms, Strain is a geometric quantity that measures the deformation of a body. There are two types of strain: normal strain: characterizes dimensional changes,
More informationSupplementary Information. for. Origami based Mechanical Metamaterials
Supplementary Information for Origami based Mechanical Metamaterials By Cheng Lv, Deepakshyam Krishnaraju, Goran Konjevod, Hongyu Yu, and Hanqing Jiang* [*] Prof. H. Jiang, C. Lv, D. Krishnaraju, Dr. G.
More informationLevel 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method
9210-203 Level 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method You should have the following for this examination one answer book No additional data is attached
More informationBHAR AT HID AS AN ENGIN E ERI N G C O L L E G E NATTR A MPA LL I
BHAR AT HID AS AN ENGIN E ERI N G C O L L E G E NATTR A MPA LL I 635 8 54. Third Year M E C H A NICAL VI S E M ES TER QUE S T I ON B ANK Subject: ME 6 603 FIN I T E E LE ME N T A N A L YSIS UNI T - I INTRODUCTION
More informationElements of Rock Mechanics
Elements of Rock Mechanics Stress and strain Creep Constitutive equation Hooke's law Empirical relations Effects of porosity and fluids Anelasticity and viscoelasticity Reading: Shearer, 3 Stress Consider
More informationParametric Study Of The Material On Mechanical Behavior Of Pressure Vessel
Parametric Study Of The Material On Mechanical Behavior Of Pressure Vessel Dr. Mohammad Tariq Assistant Professor Mechanical Engineering Department SSET, SHIATS-Deemed University, Naini, Allahabad, U.P.,
More informationNumerical analyses of cement-based piezoelectric smart composites
Numerical analyses of cement-based piezoelectric smart composites *Jan Sladek 1, Pavol Novak 2, Peter L. Bishay 3, and Vladimir Sladek 1 1 Institute of Construction and Architecture, Slovak Academy of
More informationElasticity in two dimensions 1
Elasticity in two dimensions 1 Elasticity in two dimensions Chapters 3 and 4 of Mechanics of the Cell, as well as its Appendix D, contain selected results for the elastic behavior of materials in two and
More informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Thermomechanics
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Thermomechanics Prof. Dr. Eleni Chatzi Dr. Giuseppe Abbiati, Dr. Konstantinos Agathos Lecture 13-14 December, 2017 1 / 30 Forewords
More informationChapter 3 Variational Formulation & the Galerkin Method
Institute of Structural Engineering Page 1 Chapter 3 Variational Formulation & the Galerkin Method Institute of Structural Engineering Page 2 Today s Lecture Contents: Introduction Differential formulation
More informationInstabilities and Dynamic Rupture in a Frictional Interface
Instabilities and Dynamic Rupture in a Frictional Interface Laurent BAILLET LGIT (Laboratoire de Géophysique Interne et Tectonophysique) Grenoble France laurent.baillet@ujf-grenoble.fr http://www-lgit.obs.ujf-grenoble.fr/users/lbaillet/
More informationFig. 1. Circular fiber and interphase between the fiber and the matrix.
Finite element unit cell model based on ABAQUS for fiber reinforced composites Tian Tang Composites Manufacturing & Simulation Center, Purdue University West Lafayette, IN 47906 1. Problem Statement In
More informationEarthquake response analysis of rock-fall models by discontinuous deformation analysis
c Earthquake response analysis of rock-fall models by discontinuous deformation analysis T. Sasaki, I. Hagiwara & K. Sasaki Rock Engineering Lab., Suncoh Consultants Co. Ltd., Tokyo, Japan R. Yoshinaka
More informationTopic 5: Finite Element Method
Topic 5: Finite Element Method 1 Finite Element Method (1) Main problem of classical variational methods (Ritz method etc.) difficult (op impossible) definition of approximation function ϕ for non-trivial
More informationENGN 2340 Final Project Report. Optimization of Mechanical Isotropy of Soft Network Material
ENGN 2340 Final Project Report Optimization of Mechanical Isotropy of Soft Network Material Enrui Zhang 12/15/2017 1. Introduction of the Problem This project deals with the stress-strain response of a
More informationDiscrete Analysis for Plate Bending Problems by Using Hybrid-type Penalty Method
131 Bulletin of Research Center for Computing and Multimedia Studies, Hosei University, 21 (2008) Published online (http://hdl.handle.net/10114/1532) Discrete Analysis for Plate Bending Problems by Using
More informationBasic Energy Principles in Stiffness Analysis
Basic Energy Principles in Stiffness Analysis Stress-Strain Relations The application of any theory requires knowledge of the physical properties of the material(s) comprising the structure. We are limiting
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 informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 13
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras (Refer Slide Time: 00:25) Module - 01 Lecture - 13 In the last class, we have seen how
More informationConstitutive models. Constitutive model: determines P in terms of deformation
Constitutive models Constitutive model: determines P in terms of deformation Elastic material: P depends only on current F Hyperelastic material: work is independent of path strain energy density function
More informationSEISMOLOGY I. Laurea Magistralis in Physics of the Earth and of the Environment. Elasticity. Fabio ROMANELLI
SEISMOLOGY I Laurea Magistralis in Physics of the Earth and of the Environment Elasticity Fabio ROMANELLI Dept. Earth Sciences Università degli studi di Trieste romanel@dst.units.it 1 Elasticity and Seismic
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 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 informationA unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation
A unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation S. Bordère a and J.-P. Caltagirone b a. CNRS, Univ. Bordeaux, ICMCB,
More informationPART I. Basic Concepts
PART I. Basic Concepts. Introduction. Basic Terminology of Structural Vibration.. Common Vibration Sources.. Forms of Vibration.3 Structural Vibration . Basic Terminology of Structural Vibration The term
More informationAnalytical Mechanics: Variational Principles
Analytical Mechanics: Variational Principles Shinichi Hirai Dept. Robotics, Ritsumeikan Univ. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Variational Principles 1 / 71 Agenda
More information202 Index. failure, 26 field equation, 122 force, 1
Index acceleration, 12, 161 admissible function, 155 admissible stress, 32 Airy's stress function, 122, 124 d'alembert's principle, 165, 167, 177 amplitude, 171 analogy, 76 anisotropic material, 20 aperiodic
More informationPURE BENDING. If a simply supported beam carries two point loads of 10 kn as shown in the following figure, pure bending occurs at segment BC.
BENDING STRESS The effect of a bending moment applied to a cross-section of a beam is to induce a state of stress across that section. These stresses are known as bending stresses and they act normally
More information