Supplementary Information. for. Origami based Mechanical Metamaterials
|
|
- Melanie Sparks
- 5 years ago
- Views:
Transcription
1 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. Konjevod School for ngineering of Matter, Transport and nergy Arizona State University, Tempe, AZ 8587 (USA) -mail: [*] Prof. H. Yu School of arth and Space xploration School of lectrical, Computer and nergy ngineering Arizona State University, Tempe, AZ 8587 (USA)
2 Change of Length W The length W is given by q. () of the main text. Figure Sa shows the derivatives of W s two terms, i.e., ( n ) cos β η φ cos ( φ / ) and cos / φ, along with W, as a function of φ, for n = =, η = / and β = 78.5, the same parameters used in Fig. b. It can be seem that these n 5 two derivatives work against each other with ( n ) cos β η > 0 φ cos ( φ / ) to decrease W, while with cos / < 0 φ to increase W, from a planar state to a collapsed state. Therefore, the one among these two derivatives with larger absolute value dominates the change of W. It is apparent that when one folds a Miura-ori from its planar state to a collapsed state, ( n ) cos β η > 0 φ cos ( φ / ) dominates firstly to decrease W. Once the stationary point is reached, increase W. cos / < 0 φ starts to dominate and Out-of-plane Poisson s Ratio Using q. () in the main text, the out-of-plane Poisson s ratios, related to in-plane strains upon height change, are given by β ε cot sin sin 3 = = ε cos β ε 0 ν ν ( ) ( n ) ( ) ( n ) + ( ) β φ η β φ ε sin sin cos cos 3 = = ε cos β ηcos β cos φ ε 0 (S) dh where ε = is the strain in x 3 -direction. Figure Sb shows that ν 3 is positive across the entire H range of φ, due to φ 0,β, and monotonically decreases from to 0 as the Miura-ori varies from
3 its completely collapsed state ( φ = 0 ) to the planar state ( φ β = ), i.e., [ ] ν. For an extreme 3 0, case when β 90 or φ 0, ν3. Figure Sc shows ν 3 for n = 3 and η = 0.5. It is observed that ν 3 can be both negative and positive, separated by a boundary defined by ( n ) η β cos = cos which is shown as the white dashed line. Figure Sc also shows that ν 3 [, ]. Nonlinearity of the lastic nergy with respect to the Coordinates of Vertices Depending on the type of rigid origami and the number of dihedral angles per unit cell, the elastic energy can be always expressed by U T ( α α ), n = k, (S) total i i i eq i= α i where k i are the stiffness constants of the dihedral angles α i with i, eq α as the equilibrium angle, and T n is the number of types of dihedral angles. To obtain the stiffness matrix K and non-equilibrium force P, we need to express U total in terms of the coordinates of vertices. Considering the case of Miura-ori where T α n =, the two dihedral angles are given as cos sin β sin = α = β φ ( ) cos cot tan, (S3) where φ is the projection angle between two ridges. β 0,90 and φ 0,β. Using the distances between the three vertices,, 3 (Fig. c) that form this angle, φ can be determined by using the cosine rule, R R3 R3 φ = cos +, (S4) RR3 3
4 where R ij is the distance between vertices i and j. By combining qs. (S) to (S4), a relationship between the elastic energy and the coordinates of vertices can be obtained. Clearly, this relationship is nonlinear, i.e., U total is nonlinear with respect to the coordinates of vertices. Bulk dulus of Miura-ori The bulk modulus K of Miura-ori can be defined by K θ =, (S5) p p= 0 where p is the hydrostatic pressure. Using the principle of superposition, when only p( = σ ) is applied and σ σ 0 = =, we have ( ) θ = ν ν ε, where σ, σ σ are normal stresses in 3 x, x, and x 3 directions, respectively. Using a similar approach for σ and σ, the bulk modulus K K ε ε ε = ν ν + ν ν + ν ν, where σ = σ = σ = p. is given by ( ) ( ) ( ) σ σ σ For vanishing nominal stress σ = σ = σ 0, the tensile moduli are given by ε ε ε σ σ σ =, =, =. (S6) σ = 0 σ = 0 σ = 0 Thus the bulk modulus K is given by ν ν ν ν ν ν K = + +. (S7) Here the tensile moduli are the tangential moduli of the stress-strain curve, i.e., dσ dσ dσ =, =, =. dε dε d 0 dε ε = dε = 0 dε = 0 4
5 Work Conjugate Relation Stress and duli for Miura-ori The elastic energy density W tot for a (n, n ) Miura-ori is given by W tot = V ( n )( n ) k ( α α, eq ) + ( n )( n ) k ( α α, eq ), (S8) where ( ) ( ) ( ) β φ β φ sin V = ( n ) ab n acos + bcos sin sin cos (S9) is the volume of this (n, n ) Miura-ori. The work conjugate relation provides stress by taking derivatives of W tot with respect to strains, i.e., W W W σ =, σ =, σ =. tot tot tot ε ε ε (S0) Since the Miura-ori is a periodic structure and for a given Miura-ori (i.e., fixed (n, n ), a, b, and β), the deformation can be solely determined by a single parameter φ, these derivatives can be implemented by taking derivatives with respect to φ, i.e., W / φ W / φ W / φ σ =, σ =, σ =. ε / φ ε / φ ε / φ tot tot tot (S) The strains are explicitly given by dl ε = = cot dφ L dw tan cos cos ε = = dφ. (S) W n a b ( n ) a β b ( ) cos β + cos tan cos ( β) dφ ( β) dh ε = = H sin sin Thus the stresses are obtained as 5
6 tan / σ = κ V σ σ ( n ) a β + b V ( n ) acos β bcos ( φ / ) β cot / cos cos / = κ, (S3) cot / sin sin / = cos β κ where cos / κ = ( n )( n ) k ( α α, eq ) sin β sin /. (S4) cos β dv + ( n )( n ) k ( α α, eq ) Wtot cos / sin sin / d φ β The moduli are given by dσ dσ / dφ = = dε dε / dφ ε = 0 dφ= 0 dσ dσ / dφ = = dε dε / dφ ε = 0 dφ= 0 dσ dσ / dφ = = dε dε / dφ ε = 0 dφ= 0. (S5) Implementation of q. (S5) leads to k = ab ζ ξ tan k ζ = ξ η β φ ( ) ( ) n ( ) ηcos β + cos φ φ ( n ) cos cos ( ) ab cot k ζ = ab ξ β cot sin sin 4 cos β (S6) where = ( n )( n ) ( ) + ( n )( n ) ζ 4 cos φ cos β 4 and ( n ) sin ( ) ( n ) cos cos ( ) sin sin ( ) 3 ξ = φ η β + φ β φ. Here k is the spring constant of the hinges for dihedral angles for Miura-ori. 6
7 Range of Tensile and Bulk dulus Given that n [ 3, ], [ 3, ] n, and φ 0,β collapsed state ( φ = 0 ) and for the planar state ( φ β, [ ] 0, = ), [ ] 0, with 0 for the completely with 0 for the planar state ( φ = β ) and for the completely collapsed state ( φ = 0 ). varies from a finite positive value (depending on n, n, and φ) to at both the planar and completely collapsed states. Figures Sa-c show the tensile moduli,, and normalized by k / ab as a function of φ for a few representative n and n, and η = /, β = 45. Above discussed trends are observed. Now we study the bulk modulus using q. (S7). Since some extreme values (e.g., 0,, and ) present in either tensile moduli or Poisson's ratios, it is interesting to study the extreme values of K. At ν ν3 φ 0, ν ν, 3 ν3 ν3 0, and 0, thus the bulk modulus K 0. ν ν3 ν ν3 At φ β, 0, ν3 ν3 0, and, thus the bulk modulus K 0. Another interesting point is at a particular state (φ) for the prescribed n, n and β, the right hand of q. (S7) vanishes, which provides an infinity bulk modulus. The vanishing of the right hand of q. (S7) is determined by the numerators, specifically ν ν3, ν ν3, and ν3 ν3, and again the condition of vanishing these three terms is given by their numerators. Based on q. (4) in the main text and q. (S), the numerators of ν ν3, ν ν3, and ν3 ν3 happen to be identical, 6 4 ( + β) ( n ) ( )( ) ( n ) cos cos cos. (S7) + ηcos φ cos β + cos β ηcos β 3 3 The vanishing of the expression (S7) provides a particular state of folding characterized by the angle φ and dependent on n to reach an infinity bulk modulus. Figure Sd shows the bulk modulus K 7
8 normalized by k / ab as a function of φ for a few representative n and n, and / η =, β = 45, where the signature of changing from 0 to and then 0 is represented. Non-local Interactions in the Ron Resch Pattern Figure S3a shows the planar state of a Ron Resch pattern, which features some equilateral triangles connected by some right triangles. The insets of Fig. S3a show two different folded states of Ron Resch patterns with the upper left one for a dome shape and the upper right one for a completely collapsed state or namely, a Ron Resch plate. Three dihedral angles, β, β, and β 3 are required to describe this rigid origami folding (Fig. S3b). When β = β = β3 = 80, i.e., all triangles are in the same plane, it represents a planar state (e.g., Fig. S3a). When β 0,80, β 0,80, and β 3 0,80, it corresponds to a curved state, illustrated by the upper left inset of Fig. S3a as an example. When β =, β = 0, β 3 = 90, it describes another planar but more compact state (illustrated by the upper 0 right inset of Fig. S3a), by the name of a Ron Resch plate. It is noticed that there are two types of vertices in a Ron Resch pattern, specifically, the centroids of the equilateral triangles (e.g., the vertex marked by a solid blue dot in Fig. S3a) and vertices between the right and equilateral triangles (e.g., the vertex marked by an open blue dot in Fig. S3a). The non-local feature can be similarly observed from these two vertices. For example for the solid blue vertex, it is seen that its motion influences its nearest-neighbor vertices (i.e., the ones marked by solid red dots) through dihedral angles β and β, and its second-neighbor vertices (i.e., the ones marked by solid green dots) through dihedral angles β 3. Buckling Analysis of a Ron Resch Plate and a Six-Fold Supporting Structure These two structures are all periodic so that only the unit cells are utilized to conduct the buckling analysis. Figure S6 shows the unit cell of these two structures. Same thickness, height of supporting, elastic modulus and Poisson s ratio are assigned to the two models. The finite element package 8
9 ABAQUS is used, where the eigenvalues of the different modes can be calculated by using the built-in buckling module. Critical load then can be obtained by using th i cr P th i = Pλ (S8) where th i cr P is the critical load for the i th mode, P is the infinitesimal load applied in the simulation, th i λ is the i th eigenvalue. For the Ron Resch plate, 84,58 R3 (3-node triangular shell) elements are used, with the fixed displacement boundary conditions along the in-plane directions of the plate and the spike. Contact at the spike between the ground plane and the Ron Resch plate is considered. A very small concentrated load is applied at the centroid of the plate. For the six-fold supporting structure, 49,68 S4 (4-node doubly curved shell) elements are used. Same boundary conditions and loads are applied. The cross-sectional properties and material properties of these two structures are the same. Here only the first buckling mode is concerned. The buckling analysis shows that the critical load for the Ron Resch plate is 57% higher than that for the six-fold supporting structure. 9
10 Figures (a) (b) (c) Supplementary Figure S. Geometry characteristics of Miura-ori. (a) Change of size W. Here the change rate of W s two terms with respective to φ is also shown. (b) Contour plot of out-of-plane Poisson's ratio ν 3 as a function of φ and β. (c) Contour plot of out-of-plane Poisson's ratio ν 3 as a function of φ and β. 0
11 (a) (b) (c) (d) Supplementary Figure S. Tensile and bulk moduli for the Miura-ori. (a) Tensile modulus as a function of φ for (3,3) and (3,3) Miura-ori. (b) Tensile modulus as a function of φ for (3,3) and (3,3) Miura-ori. (c) Tensile modulus as a function of φ for (3,3) and (3,3) Miura-ori. (d) Bulk modulus K as a function of φ for (3,3) and (3,3) Miura-oris. Here η = /, β = 45 and all moduli are normalized by k / ab.
12 (a) (b) (c) (d) Supplementary Figure S3. Ron Resch pattern. (a) The planar state of a Ron Resch pattern, where the solid lines are for mountain creases and the dashed lines are for valley creases. Insets are two different folded states. On the upper left is a dome shape and the upper right is a completely collapsed state. (b) Three dihedral angles β, β, and β 3 are needed to describe a Ron Resch pattern. (c) One type of non-local element for the Ron Resch pattern with the centroid of the equilateral triangle as the central vertex. (d) Another type of non-local element for the Ron Resch pattern with the intersections between pleated triangles as the central vertex.
13 Supplementary Figure S4. Shear deformation of a Miura-ori. Deformation of a (3,3) Miura-ori under shear loading along the negative x direction. It is observed that the opposite relationship between shear loading and shear deformation. 3
14 (a) (b) (c) Supplementary Figure S5. Histograms of the three dihedral angles (a) β, (b) β, and (c) β 3 for three Ron Resch patterns, namely a Ron Resch dome, a tube and a stingray. 4
15 (a) (b) Supplementary Figure S6. Unit cell for buckling analysis. (a) Completed collapsed Ron Resch plate. (b) Six-fold supporting structure. 5
16 References. B. Liu, Y. Huang, H. Jiang, S. Qu, K. C. Hwang, Computer Methods in Applied Mechanics and ngineering 93, 849 (004).. B. Liu et al., Physical Review B 7, 8 (Jul, 005). 6
A 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 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 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 informationTheoretical and Finite Element Analysis of. Origami and Kirigami Based Structures. Cheng Lv
Theoretical and Finite Element Analysis of Origami and Kirigami Based Structures by Cheng Lv A Dissertation Presented in Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy Approved
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 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 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 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 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 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 informationNonlinear Mechanics of Monolayer Graphene Rui Huang
Nonlinear Mechanics of Monolayer Graphene Rui Huang Center for Mechanics of Solids, Structures and Materials Department of Aerospace Engineering and Engineering Mechanics The University of Texas at Austin
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 informationBending of Simply Supported Isotropic and Composite Laminate Plates
Bending of Simply Supported Isotropic and Composite Laminate Plates Ernesto Gutierrez-Miravete 1 Isotropic Plates Consider simply a supported rectangular plate of isotropic material (length a, width b,
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 information46th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference April 2005 Austin, Texas
th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics & Materials Conference - April, Austin, Texas AIAA - AIAA - Bi-stable Cylindrical Space Frames H Ye and S Pellegrino University of Cambridge, Cambridge,
More informationEE C247B ME C218 Introduction to MEMS Design Spring 2017
247B/M 28: Introduction to MMS Design Lecture 0m2: Mechanics of Materials CTN 2/6/7 Outline C247B M C28 Introduction to MMS Design Spring 207 Prof. Clark T.- Reading: Senturia, Chpt. 8 Lecture Topics:
More informationSupplementary Material
Mangili et al. Supplementary Material 2 A. Evaluation of substrate Young modulus from AFM measurements 3 4 5 6 7 8 Using the experimental correlations between force and deformation from AFM measurements,
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 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 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 informationARTICLE A-8000 STRESSES IN PERFORATED FLAT PLATES
ARTICLE A-8000 STRESSES IN PERFORATED FLAT PLATES Delete endnote 18, which says "Express metric values in exponential form" A-8100 INTRODUCTION A-8110 SCOPE (a) This Article contains a method of analysis
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 informationEMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading
MA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading MA 3702 Mechanics & Materials Science Zhe Cheng (2018) 2 Stress & Strain - Axial Loading Statics
More informationLocalization in Undrained Deformation
Localization in Undrained Deformation J. W. Rudnicki Dept. of Civil and Env. Engn. and Dept. of Mech. Engn. Northwestern University Evanston, IL 621-319 John.Rudnicki@gmail.com Fourth Biot Conference on
More informationSupplementary Information to the article entitled: Internally Architectured Materials with Directionally Asymmetric Friction
Supplementary Information to the article entitled: Internally Architectured Materials with Directionally Asymmetric Friction Ehsan Bafekrpour 1,2, Arcady Dyskin 3, Elena Pasternak 4, Andrey Molotnikov
More informationCollapse behaviour and simplified modeling of triangular cross-section columns
Indian Journal of ngineering & Materials Sciences Vol. 16, April 2009, pp. 71-78 Collapse behaviour and simplified ing of triangular cross-section columns Yucheng Liu* Department of Mechanical ngineering,
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 informationIntroduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams.
Outline of Continuous Systems. Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams. Vibrations of Flexible Strings. Torsional Vibration of Rods. Bernoulli-Euler Beams.
More informationLecture 4 Honeycombs Notes, 3.054
Honeycombs-In-plane behavior Lecture 4 Honeycombs Notes, 3.054 Prismatic cells Polymer, metal, ceramic honeycombs widely available Used for sandwich structure cores, energy absorption, carriers for catalysts
More informationUnderstand basic stress-strain response of engineering materials.
Module 3 Constitutive quations Learning Objectives Understand basic stress-strain response of engineering materials. Quantify the linear elastic stress-strain response in terms of tensorial quantities
More informationDetermining the Elastic Modulus and Hardness of an Ultrathin Film on a Substrate Using Nanoindentation
Determining the Elastic Modulus and Hardness of an Ultrathin Film on a Substrate Using Nanoindentation The Harvard community has made this article openly available. Please share how this access benefits
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 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 informationCAAM 335 Matrix Analysis Planar Trusses
CAAM 5 Matrix Analysis Planar Trusses September 1, 010 1 The Equations for the Truss We consider trusses with m bars and n nodes. Each node can be displaced in horizontal and vertical direction. If the
More informationChapter 2: Elasticity
OHP 1 Mechanical Properties of Materials Chapter 2: lasticity Prof. Wenjea J. Tseng ( 曾文甲 ) Department of Materials ngineering National Chung Hsing University wenjea@dragon.nchu.edu.tw Reference: W.F.
More information20. Rheology & Linear Elasticity
I Main Topics A Rheology: Macroscopic deformation behavior B Linear elasticity for homogeneous isotropic materials 10/29/18 GG303 1 Viscous (fluid) Behavior http://manoa.hawaii.edu/graduate/content/slide-lava
More informationFEA A Guide to Good Practice. What to expect when you re expecting FEA A guide to good practice
FEA A Guide to Good Practice What to expect when you re expecting FEA A guide to good practice 1. Background Finite Element Analysis (FEA) has transformed design procedures for engineers. Allowing more
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 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 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 informationPrediction of Elastic Constants on 3D Four-directional Braided
Prediction of Elastic Constants on 3D Four-directional Braided Composites Prediction of Elastic Constants on 3D Four-directional Braided Composites Liang Dao Zhou 1,2,* and Zhuo Zhuang 1 1 School of Aerospace,
More informationCHAPTER 6 MECHANICAL PROPERTIES OF METALS PROBLEM SOLUTIONS
CHAPTER 6 MECHANICAL PROPERTIES OF METALS PROBLEM SOLUTIONS Concepts of Stress and Strain 6.1 Using mechanics of materials principles (i.e., equations of mechanical equilibrium applied to a free-body diagram),
More informationChapter 2 Examples of Optimization of Discrete Parameter Systems
Chapter Examples of Optimization of Discrete Parameter Systems The following chapter gives some examples of the general optimization problem (SO) introduced in the previous chapter. They all concern the
More informationSupplementary Information for Saddles, Twists, and Curls: Shape Transitions in Freestanding Nanoribbons
Supplementary Information for Saddles, Twists, and Curls: Shape Transitions in Freestanding Nanoribbons Hailong Wang and Moneesh Upmanyu Group for Simulation and Theory of Atomic-Scale Material Phenomena
More informationStress-Strain Behavior
Stress-Strain Behavior 6.3 A specimen of aluminum having a rectangular cross section 10 mm 1.7 mm (0.4 in. 0.5 in.) is pulled in tension with 35,500 N (8000 lb f ) force, producing only elastic deformation.
More informationMECHANICS OF 2D MATERIALS
MECHANICS OF 2D MATERIALS Nicola Pugno Cambridge February 23 rd, 2015 2 Outline Stretching Stress Strain Stress-Strain curve Mechanical Properties Young s modulus Strength Ultimate strain Toughness modulus
More informationStrength of Material. Shear Strain. Dr. Attaullah Shah
Strength of Material Shear Strain Dr. Attaullah Shah Shear Strain TRIAXIAL DEFORMATION Poisson's Ratio Relationship Between E, G, and ν BIAXIAL DEFORMATION Bulk Modulus of Elasticity or Modulus of Volume
More informationPractice Final Examination. Please initial the statement below to show that you have read it
EN175: Advanced Mechanics of Solids Practice Final Examination School of Engineering Brown University NAME: General Instructions No collaboration of any kind is permitted on this examination. You may use
More informationComb resonator design (2)
Lecture 6: Comb resonator design () -Intro Intro. to Mechanics of Materials School of Electrical l Engineering i and Computer Science, Seoul National University Nano/Micro Systems & Controls Laboratory
More information2/28/2006 Statics ( F.Robilliard) 1
2/28/2006 Statics (.Robilliard) 1 Extended Bodies: In our discussion so far, we have considered essentially only point masses, under the action of forces. We now broaden our considerations to extended
More informationMechanics of Materials II. Chapter III. A review of the fundamental formulation of stress, strain, and deflection
Mechanics of Materials II Chapter III A review of the fundamental formulation of stress, strain, and deflection Outline Introduction Assumtions and limitations Axial loading Torsion of circular shafts
More informationCritical Strain of Carbon Nanotubes: An Atomic-Scale Finite Element Study
X. Guo A. Y. T. Leung 1 e-mail: bcaleung@cityu.edu.hk Department of Building and Construction, City University of Hong Kong, Hong Kong, China H. Jiang Department of Mechanical and Aerospace Engineering,
More informationAvailable online at ScienceDirect. Procedia IUTAM 13 (2015 ) 82 89
Available online at www.sciencedirect.com ScienceDirect Procedia IUTAM 13 (215 ) 82 89 IUTAM Symposium on Dynamical Analysis of Multibody Systems with Design Uncertainties The importance of imperfections
More informationUnit 18 Other Issues In Buckling/Structural Instability
Unit 18 Other Issues In Buckling/Structural Instability Readings: Rivello Timoshenko Jones 14.3, 14.5, 14.6, 14.7 (read these at least, others at your leisure ) Ch. 15, Ch. 16 Theory of Elastic Stability
More informationInfluence of residual stresses in the structural behavior of. tubular columns and arches. Nuno Rocha Cima Gomes
October 2014 Influence of residual stresses in the structural behavior of Abstract tubular columns and arches Nuno Rocha Cima Gomes Instituto Superior Técnico, Universidade de Lisboa, Portugal Contact:
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 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 information: APPLIED MECHANICS & STRENGTH OF MATERIALS COURSE CODE : 4021 COURSE CATEGORY : A PERIODS/ WEEK : 5 PERIODS/ SEMESTER : 75 CREDIT : 5 TIME SCHEDULE
COURSE TITLE : APPLIED MECHANICS & STRENGTH OF MATERIALS COURSE CODE : 4021 COURSE CATEGORY : A PERIODS/ WEEK : 5 PERIODS/ SEMESTER : 75 CREDIT : 5 TIME SCHEDULE MODULE TOPIC PERIODS 1 Simple stresses
More informationCRITERIA FOR SELECTION OF FEM MODELS.
CRITERIA FOR SELECTION OF FEM MODELS. Prof. P. C.Vasani,Applied Mechanics Department, L. D. College of Engineering,Ahmedabad- 380015 Ph.(079) 7486320 [R] E-mail:pcv-im@eth.net 1. Criteria for Convergence.
More informationXI. NANOMECHANICS OF GRAPHENE
XI. NANOMECHANICS OF GRAPHENE Carbon is an element of extraordinary properties. The carbon-carbon bond possesses large magnitude cohesive strength through its covalent bonds. Elemental carbon appears in
More informationAircraft Stress Analysis and Structural Design Summary
Aircraft Stress Analysis and Structural Design Summary 1. Trusses 1.1 Determinacy in Truss Structures 1.1.1 Introduction to determinacy A truss structure is a structure consisting of members, connected
More informationChapter 12. Static Equilibrium and Elasticity
Chapter 12 Static Equilibrium and Elasticity Static Equilibrium Equilibrium implies that the object moves with both constant velocity and constant angular velocity relative to an observer in an inertial
More informationMechanical Properties
Mechanical Properties Elastic deformation Plastic deformation Fracture I. Elastic Deformation S s u s y e u e T I II III e For a typical ductile metal: I. Elastic deformation II. Stable plastic deformation
More informationLearning Units of Module 3
Module Learning Units of Module M. tress-train oncepts in Three- Dimension M. Introduction to Anisotropic lasticity M. Tensorial oncept and Indicial Notations M.4 Plane tress oncept Action of force (F)
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 informationSupplementary Figure 1 Extracting process of wetting ridge profiles. a1-4, An extraction example of a ridge profile for E 16 kpa.
Supplementary Figure 1 Extracting process of wetting ridge profiles. a1-4, An extraction example of a ridge profile for E 16 kpa. An original image (a1) was binarized, as shown in a2, by Canny edge detector
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 informationLecture 15 Strain and stress in beams
Spring, 2019 ME 323 Mechanics of Materials Lecture 15 Strain and stress in beams Reading assignment: 6.1 6.2 News: Instructor: Prof. Marcial Gonzalez Last modified: 1/6/19 9:42:38 PM Beam theory (@ ME
More informationELASTICITY (MDM 10203)
ELASTICITY () Lecture Module 3: Fundamental Stress and Strain University Tun Hussein Onn Malaysia Normal Stress inconstant stress distribution σ= dp da P = da A dimensional Area of σ and A σ A 3 dimensional
More informationSamantha Ramirez, MSE. Stress. The intensity of the internal force acting on a specific plane (area) passing through a point. F 2
Samantha Ramirez, MSE Stress The intensity of the internal force acting on a specific plane (area) passing through a point. Δ ΔA Δ z Δ 1 2 ΔA Δ x Δ y ΔA is an infinitesimal size area with a uniform force
More informationRHEOLOGY & LINEAR ELASTICITY
GG303 Lecture 20 10/25/09 1 RHEOLOGY & LINEAR ELASTICITY I Main Topics A Rheology: Macroscopic deformation behavior B Importance of fluids and fractures in deformation C Linear elasticity for homogeneous
More informationN = Shear stress / Shear strain
UNIT - I 1. What is meant by factor of safety? [A/M-15] It is the ratio between ultimate stress to the working stress. Factor of safety = Ultimate stress Permissible stress 2. Define Resilience. [A/M-15]
More informationStrain Transformation equations
Strain Transformation equations R. Chandramouli Associate Dean-Research SASTRA University, Thanjavur-613 401 Joint Initiative of IITs and IISc Funded by MHRD Page 1 of 8 Table of Contents 1. Stress transformation
More informationThe effect of plasticity in crumpling of thin sheets: Supplementary Information
The effect of plasticity in crumpling of thin sheets: Supplementary Information T. Tallinen, J. A. Åström and J. Timonen Video S1. The video shows crumpling of an elastic sheet with a width to thickness
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 informationDynamic Analysis of Laminated Composite Plate Structure with Square Cut-Out under Hygrothermal Load
Dynamic Analysis of Laminated Composite Plate Structure with Square Cut-Out under Hygrothermal Load Arun Mukherjee 1, Dr. Sreyashi Das (nee Pal) 2 and Dr. A. Guha Niyogi 3 1 PG student, 2 Asst. Professor,
More informationContent. Department of Mathematics University of Oslo
Chapter: 1 MEK4560 The Finite Element Method in Solid Mechanics II (January 25, 2008) (E-post:torgeiru@math.uio.no) Page 1 of 14 Content 1 Introduction to MEK4560 3 1.1 Minimum Potential energy..............................
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 informationCOURSE TITLE : APPLIED MECHANICS & STRENGTH OF MATERIALS COURSE CODE : 4017 COURSE CATEGORY : A PERIODS/WEEK : 6 PERIODS/ SEMESTER : 108 CREDITS : 5
COURSE TITLE : APPLIED MECHANICS & STRENGTH OF MATERIALS COURSE CODE : 4017 COURSE CATEGORY : A PERIODS/WEEK : 6 PERIODS/ SEMESTER : 108 CREDITS : 5 TIME SCHEDULE MODULE TOPICS PERIODS 1 Simple stresses
More informationProperties of a chiral honeycomb with a Poisson's ratio -1
Properties of a chiral honeycomb with a Poisson's ratio -1 D. Prall, R. S. Lakes Int. J. of Mechanical Sciences, 39, 305-314, (1996). Abstract A theoretical and experimental investigation is conducted
More informationAPPLICATION OF ISOTENSOID-BASED CROSS SECTIONS TO FILAMENT-WOUND TOROIDAL PRESSURE VESSELS
APPLICATION OF ISOTENSOID-BASED CROSS SECTIONS TO FILAMENT-WOUND TOROIDAL PRESSURE VESSELS L. Zu, S. Koussios and A. Beukers Design and Production of Composite Structures, Faculty of Aerospace Engineering
More informationAuthor(s) Senda, Kei; Petrovic, Mario; Nakani. Citation Acta Astronautica (2015), 111:
TitleWrinkle generation in shear-enforce Author(s) Senda, Kei; Petrovic, Mario; Nakani Citation Acta Astronautica (2015), 111: 110- Issue Date 2015-06 URL http://hdl.handle.net/2433/196845 2015 IAA. Published
More informationNumerical Modelling of Dynamic Earth Force Transmission to Underground Structures
Numerical Modelling of Dynamic Earth Force Transmission to Underground Structures N. Kodama Waseda Institute for Advanced Study, Waseda University, Japan K. Komiya Chiba Institute of Technology, Japan
More informationUniversity of Sheffield The development of finite elements for 3D structural analysis in fire
The development of finite elements for 3D structural analysis in fire Chaoming Yu, I. W. Burgess, Z. Huang, R. J. Plank Department of Civil and Structural Engineering StiFF 05/09/2006 3D composite structures
More informationAnalysis of a portal steel frame subject to fire by use of a truss model
Analysis of a portal steel frame subject to fire by use of a truss model P. G. Papadopoulos & A. Mathiopoulou Department of Civil Engineering, Aristotle University of Thessaloniki, Greece Abstract A plane
More informationComb Resonator Design (2)
Lecture 6: Comb Resonator Design () -Intro. to Mechanics of Materials Sh School of felectrical ti lengineering i and dcomputer Science, Si Seoul National University Nano/Micro Systems & Controls Laboratory
More informationCourse: US01CPHY01 UNIT 1 ELASTICITY I Introduction:
Course: US0CPHY0 UNIT ELASTICITY I Introduction: If the distance between any two points in a body remains invariable, the body is said to be a rigid body. In practice it is not possible to have a perfectly
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 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 informationHomework Problems. ( σ 11 + σ 22 ) 2. cos (θ /2), ( σ θθ σ rr ) 2. ( σ 22 σ 11 ) 2
Engineering Sciences 47: Fracture Mechanics J. R. Rice, 1991 Homework Problems 1) Assuming that the stress field near a crack tip in a linear elastic solid is singular in the form σ ij = rλ Σ ij (θ), it
More informationRe-entrant origami-based metamaterials with negative Poisson s ratio and bistability
UPPLEENTAL ATERIAL Re-entrant origami-based metamaterials with negative Poisson s ratio and bistability April 3 rd, 5 H. Yasuda and J. Yang * Department of echanical Engineering, University of Washington,
More informationREPORT TO D-FLEX LTD APRIL 2009 REPORT. (to D-Flex Ltd, NWDA Innovation Vouchers Award, 28 April 2009)
REPORT (to D-Flex Ltd, NWDA Innovation Vouchers Award, 28 April 2009) Finite Element Modelling of Rigidity of Fishing Rods with Various Cross-sections Report to D-Flex Ltd by Z Guan Department of Engineering,
More informationA circular tunnel in a Mohr-Coulomb medium with an overlying fault
MAP3D VERIFICATION EXAMPLE 9 A circular tunnel in a Mohr-Coulomb medium with an overlying fault 1 Description This example involves calculating the stresses and displacements on a fault overlying a 5 m
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 informationMAAE 2202 A. Come to the PASS workshop with your mock exam complete. During the workshop you can work with other students to review your work.
It is most beneficial to you to write this mock final exam UNDER EXAM CONDITIONS. This means: Complete the exam in 3 hours. Work on your own. Keep your textbook closed. Attempt every question. After the
More informationME 383S Bryant February 17, 2006 CONTACT. Mechanical interaction of bodies via surfaces
ME 383S Bryant February 17, 2006 CONTACT 1 Mechanical interaction of bodies via surfaces Surfaces must touch Forces press bodies together Size (area) of contact dependent on forces, materials, geometry,
More informationLINEAR AND NONLINEAR BUCKLING ANALYSIS OF STIFFENED CYLINDRICAL SUBMARINE HULL
LINEAR AND NONLINEAR BUCKLING ANALYSIS OF STIFFENED CYLINDRICAL SUBMARINE HULL SREELATHA P.R * M.Tech. Student, Computer Aided Structural Engineering, M A College of Engineering, Kothamangalam 686 666,
More informationBUCKLING OF CYLINDRICAL COMPOSITE SHELLS UNDER DEFORMATION-RELATED LOADS BY USING FINITE STRIP METHOD
BUCKLING OF CYLINDRICAL COMPOSITE SHELLS UNDER DEFORMATION-RELATED LOADS BY USING FINITE STRIP METHOD Arman Amini Safir-E-Andisheh Language Academy, Iran Abstract Current advances in today s industries
More informationTopics in Ship Structures
Topics in Ship Structures 8 Elastic-lastic Fracture Mechanics Reference : Fracture Mechanics by T.L. Anderson Lecture Note of Eindhoven University of Technology 17. 1 by Jang, Beom Seon Oen INteractive
More informationEsben Byskov. Elementary Continuum. Mechanics for Everyone. With Applications to Structural Mechanics. Springer
Esben Byskov Elementary Continuum Mechanics for Everyone With Applications to Structural Mechanics Springer Contents Preface v Contents ix Introduction What Is Continuum Mechanics? "I Need Continuum Mechanics
More information