The influence of architecture on the deformation localisation of cellular materials
|
|
- Stephany O’Connor’
- 6 years ago
- Views:
Transcription
1 The influence of architecture on the deformation localisation of cellular materials Dunlop John 1, Weinkamer Richard 1, Fratzl Peter 1, and Bréchet Yves 2 (1) Max Planck Institute of Colloids and Interfaces, Department of Biomaterials Potsdam, Germany (2) GPM2, Domaine Universitaire de Grenoble, Saint Martin d Heres Cedex, France john.dunlop@mpikg.mpg.de Abtract : In general cellular materials fail by localisation of deformation. For random cellular materials, the localisation plane due to compressive deformation is found to be perpendicular to the applied stress. For periodic cellular materials, however the localisation occurs in particular directions that are related to the underlying crystallography of the cellular architecture. In this contribution we treat the deformation localisation like a martensitic instability, that is a deformation corresponding to a soft elastic deformation mode. The model of Khachaturyan is applied to different cellular architectures both in two and three dimensions. Résumé : En générale les materiaux cellulaires rompent par localisation de la déformation. Pour les matériaux poreux non-ordonnés, le plan de localisation se trouve perpendiculaire a la contrainte, par contre pour les matériaux périodiques la localisation suit des directions bien particulières liées à la «cristallographie» de l architecture. Dans cette contribution nous analysons la relation entre l architecture et l orientation spatiale des bandes de localisation/instabilité dans un materiau cellulaire périodique. Nous appliquons les méthodes développées par Khachaturyan et nous traitons la localisation de la déformation comme un instabilité de type martensitique, c'est-à-dire correspondant à un mode élastique «mou». Le modèle de Khachaturyan est appliqué sur les structures cellulaires pour les architectures différentes en deux et trois dimensions. Key-words : Cellular Materials; Deformation Localisation; Lattice Dynamics 1 Introduction The interest in cellular materials has developed strongly in the last decades principally due to the possibility of optimising material properties through control of the cellular architecture, that is the size, shape and spatial organisation of the pores. Many different optimisation problems can be posed, including multi-criterion optimisation, some examples being; the optimisation of mechanical strength for minimal mass, maximising vibration damping for minimal mass, maximising fluid flow for catalytic or heat transfer applications or in optimising impact energy absorption (Bendsoe and Triantafyllidis 1990; Evans et al. 2001; Neves et al. 2002; Torquato et al In order to optimise cellular materials design it is necessary to have a good model to describe materials properties as a function of the underlying architecture. These models can also be used to predict the effects of architectural modification due to disease (e.g. osteoporosis) of trabecular bone. The modelling of mechanical properties has focused mainly on the elastic properties and failure stress or strength of cellular materials through the analysis of the behaviour of an elementary cell with periodically boundary conditions (Gibson and Ashby 1997 These approaches use the hypothesis of a macroscopically homogeneous deformation. Other approaches have used finite element methods in order to analyse the post-yield behaviour and the effect of imperfections (e.g. Chen et al In general, cellular materials fail through deformation localisation. For random cellular materials such as metal or polymeric foams, the deformation localisation occurs in a direction that is perpendicular to the applied compressive load (Tu et al. 2001; Dillard et al
2 Natural cellular materials, such as trabecular bone, also localise under compressive deformation however the direct link to architecture or fabric is still not clear (Nazarian et al For periodic cellular solids however the deformation localisation occurs along well defined crystallographic directions that depend on the microstructure (e.g. Côté et al This is a very similar to the deformation of single crystals, where deformation is localised along particular slip-planes. In metal physics the lattice dynamics approach developed by Born and Huang has been used to investigate the orientation dependence of the critical ideal shear stress of single crystals (Born and Huang 1954; Cook and De Fontaine 1969; Khachaturyan 1983; Wang et al. 1993; Fratzl and Penrose 1995 Born developed a set of three criteria for the elastic stability of a crystal (Born and Huang 1954) which for a cubic lattice are : C + C 0, C 0, C C 0, (1) > 44 > > where C ij are the elastic constants. These three criterion correspond respectively to; the existence of a Bulk modulus K = ( C C 12 )/ 3, the existence of a shear modulus C 44 and the existence of a shear modulus against tetragonal shear G' = ( C 11 C12 )/ 2 (Wang et al The third criterion is generally reached first, suggesting that as a first approximation we can just use just the difference between the Young s and shear modulus as a means of determining the onset of instability or shear localisation. In crystalline or polycrystalline materials such an instability is normally not reached under mechanical loading, as typically dislocation slip is activated, but these criteria have been applied successfully to the onset of melting in metals (Wang et al In some metals however such as in zinc or cadmium, with a large anisotropy in elastic constants a shear instability can be initiated. Porous materials with an underlying architecture can have elastic moduli which depend strongly on direction, meaning that shear instabilities can be initiated for particular architectures and cell geometries. What is required is a knowledge of the elastic stiffness matrix and how this varies for different deformation modes, macroscopic or global modes or localised perturbations. In this contribution we apply the mathematical tool developed for crystal lattices to the macroscopic crystal of a periodic cellular solid. We treat the deformation localisation much like a martensitic instability, that is a deformation corresponding to a soft elastic deformation mode (Khachaturyan 1983 A similar approach has been used by Triantafyllidis and coworkers, where they use a bifurcation analysis to determine the critical loading conditions for which a periodic cellular lattice loses stability (Triantafyllidis and Bardenhagen 1996; Triantafyllidis et al The model of Khachaturyan is applied to different cellular architectures both in two and three dimensions. The next section outlines the model development and is followed by initial results on simple structures in two and three dimensions. 2 Model x 2 u 3 (r) x 1 α r u 1 (r) u 2 (r ) u 1 (r ) r u 3 (r ) r u 2 (r) r Un-deformed Local Displacements Local Rotations Figure 1: Notation used after (Kumar and McDowell
3 The general approach of Khachaturyan involves calculating the energy of the lattice as a function of the pair-potential between atoms, and the underlying crystallography (Khachaturyan 1983; Fratzl and Penrose 1996 The pair potential gives the energy of a given pair of atoms as a function of the different degrees of freedom (translation and rotation) of each atom type. The calculation of the lattice energy is greatly simplified by the crystallography, as the underlying periodicity means it is only necessary to calculate the energy of each atom pair within the unit cell. The lattice energy is the sum of all the energy of all atom pairs in the unit cell as a function of deformation and rotation of each lattice point. This energy can then be minimised as a function of the lattice deformation modes in order to determine in which direction the lattice will deform or to look at which direction phase transformations will be preferred. In the following we describe our beam lattice as consisting of joints (atoms) connected by beams which have an associated pair potential related to the deformation energy of the beam. The energy of a single beam W as a function of displacements and rotations is given by the following equation: ( r, ) = u( r, K( r, u( r, r ) W, (2) r are the coordinates of the two ends of the beam, and ( ) where r and ' K r, is the stiffness matrix. u ( r, ) is the deformation/rotation vector of both beam ends, and in two dimensions is given by: u ( r ) ( u ( r) u ( r) u ( r) u ( ) u ( ) u ( ')), r = (3) u i for i = 1, 2 corresponds to the displacements of the beam ends and for i = 3 corresponds to the in-plane rotations of the beam ends (Figure 1 In three dimensions the deformation vector will have twelve components including the 6 degrees of freedom at each beam end. The stiffness matrices K ( r, ) of a beam in 2 and 3 dimensions can be found in standard mechanics textbooks (e.g. Bazant and Cedolin 2003 In general the total lattice energy can be written as: W Lattice u can be written as the sum of deformationwaves: u (, ) u( r, K ( r, u( r, ) = W r = (4) r, r, The heterogeneous deformation field ( r) () r = u( k) k 0 ~ T ikr e, (5) where the sum is over all k except k = 0 (which corresponds to the homogeneous displacement The vector u ~ ( k) is simply the Fourier transform of the deformation vector (including rotations) in real space. Equation (5) can then be substituted into equation (4) and then rearranged: Lattice T W = u ~ Lattice ~ (6) K Lattice k ( k) ( k). K ( k). u( k) is the reciprocal stiffness matrix corresponding to the entire lattice, and will be a 3 n by 3 n matrix in 2 dimensions (or 6 n by 6 n in 3 dimensions), where n is the number of joints (atoms) in the unit cell. Triantafyllidis and Schnaidt (1993) used a stiffness matrix in which the beam energy is written as a function of the applied compressive load in the beam. 3
4 They then calculated the critical load at which the minimum eigenvalue of the stiffness matrix K Lattice ( k) tends to zero as a function of k. All the information concerning the energy of an arbitrary deformation mode is contained within the reciprocal lattice stiffness matrix. The k - dependence of the minimum eigenvalue gives information about what is the minimum energy deformation mode, and thus the potential direction dependence of localisation for arbitrary loading. Assuming shear localisation occurs in a plane, then this corresponds to looking at the eigenvalues along lines from the origin in reciprocal or k -space. The eigenvalue surface close to k = 0 gives information about long-wave deformation modes, and in the limit as k 0, the long wave limit, the slope of the eigenvalue surface gives the bulk elastic response of the lattice. A similar approach was used by Kumar and McDowell (2004) to calculate the elastic matrix for two dimensional periodic lattices. In the following we analyse minimum eigenvalue surfaces for different lattices, square, triangular, a bcc-square lattice in 2 dimensions and a cubic lattice in 3 dimensions. 3 Results for simple lattices in two and three dimensions Figure 2a illustrates the eigenvalue surfaces for a square lattice in two dimensions. In two dimensions three degrees of freedom are allowed, rotation of a beam end within the plane, and displacement in the x and y directions, meaning that there are three eigenvalues of the reciprocal stiffness matrix. We are interested in the minimum surface which will be associated with minimal energy deformation modes. Figures 2b and c shows the minimal eigenvalue surface, which appears to show a peak around the origin, however a zoom around the origin shows a finer structure consisting of valleys oriented along the 10 directions. This can also be more clearly seen in Figure 2d which shows the minimal eigenvalue as a function of the angle θ from 01 direction for two different radii. The eigenvalues of the reciprocal stiffness matrix not only give the minimum energies but also give information about the associated deformation mode; whether the mode is rotational, displacive in x or y or some combination of the three deformation modes. For example in the square lattice the lowest energy eigenmodes close to the origin correspond to pure shear modes along the 01 directions, modes further away from the origin are mixed displacement-rotational. Figure 2a: Surface of the three eigenvalues of the reciprocal stiffness matrix of a square lattice as a function of wave vector k = ( k 1, k 2 ) Figure 2b: Surface of the smallest eigenvalue of a square lattice as a function of wave vector k = ( k 1, k 2 ) For triangular lattices the minimal surface can also be calculated, as illustrated in Figure 3 The surface is peaked with a localised minimum around the origin with a circular symmetry which corresponds to the well known isotropy of the triangular lattice. Lattices containing more than one joint in the unit cell, (or a crystal consisting of interpenetrating Bravais lattices) can also be treated, however care must be taken to assemble the reciprocal stiffness matrix. As an example a 2D bcc lattice has been modelled, that is a 2D square lattice with an extra joint 4
5 located in the centre of the unit cell giving diagonal bracing. The minimal eigenvalue surface of this system is illustrated in Figures 4a and 4b. What is interesting is that although there is local minimum at the origin (Figure 4b) there are global minima in the 01 directions (Figure 4a), suggesting that localised or short-wave deformation modes may be activated in the system. Figure 2c: Zoom of the fine structure of the minimal eigenvalue surface around the origin, showing the low energy valleys in the 10 directions Figure 2d: Radial distribution of minimum eigenvalues, normalised by the Young S modulus and length around k = 0 from Figure 2a. Dotted line is for a radius of 10 6 wave vectors from the origin, and filled line is _ wave vectors to the origin Figure 3: The fine structure of the minimal eigenvalue surface of a triangular lattice, showing the circularly symmetric low eigenvalue valley around the origin Figure 4a: Surface of the minimal eigenvalue of a bcc-lattice (two joints per unit cell) as a function of wave vector k Figure 4b: Zoom of the centre of the eigenvalue surface, showing the low energy minimum valley close to the origin Figure 5: Iso-contour surface of the minimum eigenvalue of a cubic lattice as a function of the wave vector k = ( k1, k2, k3 ). Note the direction dependence of the surface in the 100 directions. In the last example the procedure developed above is applied to a three dimensional cubic beam lattice. Figure 5 illustrates an iso-contour surface of the minimum eigenvalues of the reciprocal stiffness matrix for a cubic lattice. The surface is highly oriented having wings 5
6 extending in the 100 directions, corresponding to the directions of localisation in a cubic lattice. 4 Conclusions This paper has briefly outlined the development of a model to investigate the localisation behaviour of periodic cellular materials. The underlying periodicity allows the direct transfer of techniques developed in solid state physics to describe macroscopic lattice properties. The lattice dynamics approach allows the calculation of eigenvalue surfaces or iso-surfaces which is in turn linked to the minimal energy of the system. This has been applied to several different two and three dimensional lattices and can in principal be applied to arbitrary periodic cellular solids, although the calculations will be limited by the number of degrees of freedom. The full analysis of such surfaces will need to be compared with experiments on real porous solids and future work will be done to do mechanical testing of porous solids with different architectures produced using rapid prototyping, as well as other computational experiments such as finite element analysis. References Bazant, Z. P. and L. Cedolin (2003 Stability of Structures. Mineola, Dover. Bendsoe, M. P. and N. Triantafyllidis (1990 International Journal of Solids and Structures 26(7): Born, M. and K. Huang (1954 Dynamical Theory of Crystal Lattices. Oxford, Oxford University Press. Chen, C., T. J. Lu, et al. (2001 International Journal of Mechanical Sciences 43: Cook, H. E. and D. De Fontaine (1969 Acta Metallurgica 17: Côté, F., V. Deshpande, et al. (2006 International Journal of Solids and Structures 43: Dillard, T., F. N'Guyen, et al. (2005 Philosophical Magazine 85(19): Evans, A. G., J. W. Hutchinson, et al. (2001 Progress in Materials Science 46(3-4): Fratzl, P. and O. Penrose (1995 Acta Metallurgica et Materialia 43(8): Fratzl, P. and O. Penrose (1996 Acta Materialia 44: Gibson, L. J. and M. Ashby (1997 Cellular Solids: Structure and Properties. Cambridge, Cambridge University Press. Khachaturyan, A. G. (1983 Theory of Structural Transformations in Solids. New York, John Wiley and Sons. Kumar, R. S. and D. L. McDowell (2004 International Journal of Solids and Structures 41: Nazarian, A., M. Stauber, et al. (2006 Bone 39: Neves, M. M., O. Sigmund, et al. (2002 International Journal for Numerical Methods in Engineering 54: Torquato, S., S. Hyun, et al. (2003 Journal of Applied Physics 94(9): Triantafyllidis, N. and S. Bardenhagen (1996 Journal Of The Mechanics And Physics Of Solids 44(11): Triantafyllidis, N., M. D. Nestorovic, et al. (2006 Journal of Applied Mechanics 73: Triantafyllidis, N. and W. C. Schnaidt (1993 Journal of the Mechanics and Physics of Solids 41(9): Tu, Z. H., V. P. W. Shim, et al. (2001 International Journal of Solids and Structures 38: Wang, J., S. Yip, et al. (1993 Physical Review Letters 71(25):
Regular, low density cellular structures - rapid prototyping, numerical simulation, mechanical testing
Mat. Res. Soc. Symp. Proc. Vol. 823 2004 Materials Research Society W8.8.1 Regular, low density cellular structures - rapid prototyping, numerical simulation, mechanical testing J. Stampfl 1, M.M. Seyr
More informationFOAM TOPOLOGY BENDING VERSUS STRETCHING DOMINATED ARCHITECTURES
Acta mater. 49 (001) 1035 1040 www.elsevier.com/locate/actamat FOAM TOPOLOGY BENDING VERSUS STRETCHING DOMINATED ARCHITECTURES V. S. DESHPANDE, M. F. ASHBY and N. A. FLECK Cambridge University, Department
More informationNumerical analysis for an interpretation of the pressuremeter test in cohesive soil
Numerical analysis for an interpretation of the pressuremeter test in cohesive soil J. Monnet Joseph Fourier University L3S-R Domaine Universitaire, BP n 53, 3841, Grenoble, Cedex 9, France jmonnet@ujf-grenoble.fr
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 informationMechanics of Irregular Honeycomb Structures
Mechanics of Irregular Honeycomb Structures S. Adhikari 1, T. Mukhopadhyay 1 Chair of Aerospace Engineering, College of Engineering, Swansea University, Singleton Park, Swansea SA2 8PP, UK Sixth International
More informationComputational approach to the mechanics of entangled materials
Computational approach to the mechanics of entangled materials Carine Barbier, Rémy Dendievel & David Rodney SIMAP-GPM2 101, rue de la Physique, Domaine Universitaire, BP 46, 38 402 St-Martin d Hères Cedex
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 informationCellular solid structures with unbounded thermal expansion. Roderic Lakes. Journal of Materials Science Letters, 15, (1996).
1 Cellular solid structures with unbounded thermal expansion Roderic Lakes Journal of Materials Science Letters, 15, 475-477 (1996). Abstract Material microstructures are presented which can exhibit coefficients
More informationLecture contents. Stress and strain Deformation potential. NNSE 618 Lecture #23
1 Lecture contents Stress and strain Deformation potential Few concepts from linear elasticity theory : Stress and Strain 6 independent components 2 Stress = force/area ( 3x3 symmetric tensor! ) ij ji
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 informationTable of Contents. Preface...xvii. Part 1. Level
Preface...xvii Part 1. Level 1... 1 Chapter 1. The Basics of Linear Elastic Behavior... 3 1.1. Cohesion forces... 4 1.2. The notion of stress... 6 1.2.1. Definition... 6 1.2.2. Graphical representation...
More informationModule-4. Mechanical Properties of Metals
Module-4 Mechanical Properties of Metals Contents ) Elastic deformation and Plastic deformation ) Interpretation of tensile stress-strain curves 3) Yielding under multi-axial stress, Yield criteria, Macroscopic
More informationICM11. Homogenization of periodic auxetic materials
Available online at www.sciencedirect.com Procedia Engineering 10 (2011) 1847 1852 ICM11 Homogenization of periodic auxetic materials J. Dirrenberger a,, S. Forest a, D. Jeulin a,b, C. Colin a a Centre
More informationStructural Analysis of Truss Structures using Stiffness Matrix. Dr. Nasrellah Hassan Ahmed
Structural Analysis of Truss Structures using Stiffness Matrix Dr. Nasrellah Hassan Ahmed FUNDAMENTAL RELATIONSHIPS FOR STRUCTURAL ANALYSIS In general, there are three types of relationships: Equilibrium
More informationSupplementary figures
Supplementary figures Supplementary figure 1: The fang model (longitudinal section) used for the stiffness analysis. The base of the fang is considered fixed (i.e. zero displacement constraint) and the
More informationExercise: concepts from chapter 8
Reading: Fundamentals of Structural Geology, Ch 8 1) The following exercises explore elementary concepts associated with a linear elastic material that is isotropic and homogeneous with respect to elastic
More informationELASTOPLASTICITY THEORY by V. A. Lubarda
ELASTOPLASTICITY THEORY by V. A. Lubarda Contents Preface xiii Part 1. ELEMENTS OF CONTINUUM MECHANICS 1 Chapter 1. TENSOR PRELIMINARIES 3 1.1. Vectors 3 1.2. Second-Order Tensors 4 1.3. Eigenvalues and
More informationThe fracture toughness of planar lattices: imperfection sensitivity
1 2 The fracture toughness of planar lattices: imperfection sensitivity 3 Naomi E.R. Romijn, Norman A. Fleck 4 5 Cambridge University Engineering Department, Trumpington Street, Cambridge CB2 1PZ, UK 6
More informationVORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS
The 4 th World Conference on Earthquake Engineering October -7, 008, Beijing, China VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS K.
More informationCharacterisation and Modelling of a Melt-Extruded LDPE Closed Cell Foam
Characterisation and Modelling of a Melt-Extruded LDPE Closed Cell Foam Qusai Hatem Jebur 1,a, Philip Harrrison 1,b, Zaoyang Guo,c, Gerlind Schubert 1,d & Vincent Navez 3,e 1 School of Engineering, University
More informationModule 7: Micromechanics Lecture 29: Background of Concentric Cylinder Assemblage Model. Introduction. The Lecture Contains
Introduction In this lecture we are going to introduce a new micromechanics model to determine the fibrous composite effective properties in terms of properties of its individual phases. In this model
More informationStress-strain response and fracture behaviour of plain weave ceramic matrix composites under uni-axial tension, compression or shear
Xi an 2-25 th August 217 Stress-strain response and fracture behaviour of plain weave ceramic matrix composites under uni-axial tension compression or shear Heyin Qi 1 Mingming Chen 2 Yonghong Duan 3 Daxu
More informationOn the determinacy of repetitive structures
Journal of the Mechanics and Physics of Solids 51 (2003) 383 391 www.elsevier.com/locate/jmps On the determinacy of repetitive structures S.D. Guest a;, J.W. Hutchinson b a Department of Engineering, University
More informationMultiscale analyses of the behaviour and damage of composite materials
Multiscale analyses of the behaviour and damage of composite materials Presented by Didier BAPTISTE ENSAM, LIM, UMR CNRS 8006 151 boulevard de l hôpital l 75013 PARIS, France Research works from: K.Derrien,
More informationWe need to be able to describe planes and directions.
We need to be able to describe planes and directions. Miller Indices & XRD 1 2 Determining crystal structure and identifying materials (B) Plastic deformation Plastic deformation and mechanical properties
More informationComputational Analysis for Composites
Computational Analysis for Composites Professor Johann Sienz and Dr. Tony Murmu Swansea University July, 011 The topics covered include: OUTLINE Overview of composites and their applications Micromechanics
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 informationDYNAMIC RESPONSE OF SYNTACTIC FOAM CORE SANDWICH USING A MULTIPLE SCALES BASED ASYMPTOTIC METHOD
ECCM6-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Seville, Spain, -6 June 4 DYNAMIC RESPONSE OF SYNTACTIC FOAM CORE SANDWICH USING A MULTIPLE SCALES BASED ASYMPTOTIC METHOD K. V. Nagendra Gopal a*,
More informationTENSILE CRACKING VIEWED AS BIFURCATION AND INSTABILITY IN A DISCRETE INTERFACE MODEL
Fracture Mechanics of Concrete Structures Proceeding FRAMCOS-3 AEDIFICATIO Publishers, D-79104 Frei burg, Germany TENSILE CRACKING VIEWED AS BIFURCATION AND INSTABILITY IN A DISCRETE INTERFACE MODEL A.
More information1. Background. is usually significantly lower than it is in uniaxial tension
NOTES ON QUANTIFYING MODES OF A SECOND- ORDER TENSOR. The mechanical behavior of rocks and rock-like materials (concrete, ceramics, etc.) strongly depends on the loading mode, defined by the values and
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 informationStrain localization and cyclic damage of polyurethane foam cylinders: experimental tests and theoretical model
Strain localization and cyclic damage of polyurethane foam cylinders: experimental tests and theoretical model Giampiero Pampolini, Gianpietro Del Piero To cite this version: Giampiero Pampolini, Gianpietro
More informationStructure and Dynamics : An Atomic View of Materials
Structure and Dynamics : An Atomic View of Materials MARTIN T. DOVE Department ofearth Sciences University of Cambridge OXFORD UNIVERSITY PRESS Contents 1 Introduction 1 1.1 Observations 1 1.1.1 Microscopic
More informationMaterials: engineering, science, processing and design, 2nd edition Copyright (c)2010 Michael Ashby, Hugh Shercliff, David Cebon.
Modes of Loading (1) tension (a) (2) compression (b) (3) bending (c) (4) torsion (d) and combinations of them (e) Figure 4.2 1 Standard Solution to Elastic Problems Three common modes of loading: (a) tie
More informationVIBRATION PROBLEMS IN ENGINEERING
VIBRATION PROBLEMS IN ENGINEERING FIFTH EDITION W. WEAVER, JR. Professor Emeritus of Structural Engineering The Late S. P. TIMOSHENKO Professor Emeritus of Engineering Mechanics The Late D. H. YOUNG Professor
More informationThermal buckling and post-buckling of laminated composite plates with. temperature dependent properties by an asymptotic numerical method
hermal buckling and post-buckling of laminated composite plates with temperature dependent properties by an asymptotic numerical method F. Abdoun a,*, L. Azrar a,b, E.M. Daya c a LAMA, Higher School of
More informationA Quasicontinuum for Complex Crystals
A Quasicontinuum for Complex Crystals Ellad B. Tadmor Aerospace Engineering and Mechanics University of Minnesota Collaborators: Previous Work: Current Work: U. V. Waghmare, G. S. Smith, N. Bernstein,
More informationAnalysis of high loss viscoelastic composites
Analysis of high loss viscoelastic composites by C. P. Chen, Ph.D. and R. S. Lakes, Ph.D. Department of Engineering Physics Engineering Mechanics Program; Biomedical Engineering Department Materials Science
More informationMechanical Engineering Ph.D. Preliminary Qualifying Examination Solid Mechanics February 25, 2002
student personal identification (ID) number on each sheet. Do not write your name on any sheet. #1. A homogeneous, isotropic, linear elastic bar has rectangular cross sectional area A, modulus of elasticity
More informationThis document is a preview generated by EVS
INTERNATIONAL STANDARD ISO 18437-1 First edition 2012-08-15 Mechanical vibration and shock Characterization of the dynamic mechanical properties of visco-elastic materials Part 1: Principles and guidelines
More informationNumerical modelling of elastic-plastic deformation at crack tips in composite material under stress wave loading
JOURNAL DE PHYSIQUE IV Colloque C8, supplément au Journal de Physique III, Volume 4, septembre 1994 C8-53 Numerical modelling of elastic-plastic deformation at crack tips in composite material under stress
More informationUsing the Timoshenko Beam Bond Model: Example Problem
Using the Timoshenko Beam Bond Model: Example Problem Authors: Nick J. BROWN John P. MORRISSEY Jin Y. OOI School of Engineering, University of Edinburgh Jian-Fei CHEN School of Planning, Architecture and
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting Lectures & 3, 9/31 Aug 017 www.geosc.psu.edu/courses/geosc508 Discussion of Handin, JGR, 1969 and Chapter 1 Scholz, 00. Stress analysis and Mohr Circles Coulomb Failure
More informationIntroduction to Solid State Physics or the study of physical properties of matter in a solid phase
Introduction to Solid State Physics or the study of physical properties of matter in a solid phase Prof. Germar Hoffmann 1. Crystal Structures 2. Reciprocal Lattice 3. Crystal Binding and Elastic Constants
More informationStructural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.
Outline of Multi-Degree-of-Freedom Systems Formulation of Equations of Motions. Newton s 2 nd Law Applied to Free Masses. D Alembert s Principle. Basic Equations of Motion for Forced Vibrations of Linear
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 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 information! EN! EU! NE! EE.! ij! NN! NU! UE! UN! UU
A-1 Appendix A. Equations for Translating Between Stress Matrices, Fault Parameters, and P-T Axes Coordinate Systems and Rotations We use the same right-handed coordinate system as Andy Michael s program,
More information3D simulations of an injection test done into an unsaturated porous and fractured limestone
3D simulations of an injection test done into an unsaturated porous and fractured limestone A. Thoraval *, Y. Guglielmi, F. Cappa INERIS, Ecole des Mines de Nancy, FRANCE *Corresponding author: Ecole des
More informationPost Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method
9210-220 Post Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method You should have the following for this examination one answer book scientific calculator No
More informationGeometry of Crystal Lattice
0 Geometry of Crystal Lattice 0.1 Translational Symmetry The crystalline state of substances is different from other states (gaseous, liquid, amorphous) in that the atoms are in an ordered and symmetrical
More information2.2 - Screening and ranking for optimal selection. Outline
2 - Ashby Method 2.2 - Screening and ranking for optimal selection Outline Basic steps of selection 1. Translation of design requirements into a material specification 2. Screening out of materials that
More informationA Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials
Dublin, October 2010 A Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials FracMan Technology Group Dr Mark Cottrell Presentation Outline Some Physical
More informationAEROELASTIC ANALYSIS OF SPHERICAL SHELLS
11th World Congress on Computational Mechanics (WCCM XI) 5th European Conference on Computational Mechanics (ECCM V) 6th European Conference on Computational Fluid Dynamics (ECFD VI) E. Oñate, J. Oliver
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 informationCheckerboard instabilities in topological shape optimization algorithms
Checkerboard instabilities in topological shape optimization algorithms Eric Bonnetier, François Jouve Abstract Checkerboards instabilities for an algorithm of topological design are studied on a simple
More informationMechanics of Inflatable Fabric Beams
Copyright c 2008 ICCES ICCES, vol.5, no.2, pp.93-98 Mechanics of Inflatable Fabric Beams C. Wielgosz 1,J.C.Thomas 1,A.LeVan 1 Summary In this paper we present a summary of the behaviour of inflatable fabric
More informationFailure and Failure Theories for Anisotropic Rocks
17th international Mining Congress and Exhibition of Turkey- IMCET 2001, 2001, ISBN 975-395-417-4 Failure and Failure Theories for Anisotropic Rocks E. Yaşar Department of Mining Engineering, Çukurova
More informationEnabling Technologies
Enabling Technologies Mechanical Modelling 1 Key Parameter Mineral System Exploration is reflected in scale-dependent translation A. Gradient in hydraulic potential B. Permeability C. Solubility sensitivity
More informationOn optimisation of structures under stability constraints - a simple example
Rakenteiden Mekaniikka Journal of Structural Mechanics Vol. 49, No 2, 26, pp. 69 77 rmseura.tkk.fi/rmlehti/ c The authors 26. Open access under CC BY-SA 4. license. On optimisation of structures under
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 informationCritical Borehole Orientations Rock Mechanics Aspects
Critical Borehole Orientations Rock Mechanics Aspects By R. BRAUN* Abstract This article discusses rock mechanics aspects of the relationship between borehole stability and borehole orientation. Two kinds
More informationStructural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian
Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian ahmadian@iust.ac.ir Dynamic Response of MDOF Systems: Mode-Superposition Method Mode-Superposition Method:
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 informationLecture 2: Stresses in Pavements
Lecture 2: Stresses in Pavements Stresses in Layered Systems At any point, 9 stresses exist. They are 3 normal stresses (s z, s r, s t ) and 6 shearing stresses ( t rz = t zr, t rt = t tr, and t tz = t
More informationThe acoustic characterization of porous media and its standards
The acoustic characterization of porous media and its standards Luc JAOUEN 1, François-Xavier BECOT, Fabien CHEVILLOTTE Matelys, France ABSTRACT While there is a growing number of methods for the acoustic
More informationA CLOSER LOOK AT THE BRAZILIAN TEST AND ITS MODE OF FAILURE Arvid Landva, GEMTEC Limited, Fredericton, New Brunswick, Canada
A CLOSER LOOK AT THE BRAZILIAN TEST AN ITS MOE OF FAILURE Arvid Landva, GEMTEC Limited, Fredericton, New Brunswick, Canada ABSTRACT Timoshenko (1934) showed that a compressive line load applied perpendicularly
More informationEnhancing Prediction Accuracy In Sift Theory
18 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS Enhancing Prediction Accuracy In Sift Theory J. Wang 1 *, W. K. Chiu 1 Defence Science and Technology Organisation, Fishermans Bend, Australia, Department
More informationPHASE-FIELD SIMULATION OF DOMAIN STRUCTURE EVOLUTION IN FERROELECTRIC THIN FILMS
Mat. Res. Soc. Symp. Proc. Vol. 652 2001 Materials Research Society PHASE-FIELD SIMULATION OF DOMAIN STRUCTURE EVOLUTION IN FERROELECTRIC THIN FILMS Y. L. Li, S. Y. Hu, Z. K. Liu, and L. Q. Chen Department
More informationSea to Sky Geotechnique 2006
UNDERSTANDING SLOPE DEFORMATIONS TRIGGERED BY EARTHQUAKES Muhsin Élie RAHHAL, CLERC, École Supérieure d Ingénieurs de Beyrouth, Université Saint Joseph, PO BOX 54, BEIRUT, LEBANON, emails: muhsin.rahal@fi.usj.edu.lb
More informationInternational Journal of Plasticity
International Journal of lasticity 26 (2010) 141 148 Contents lists available at ScienceDirect International Journal of lasticity journal homepage: www.elsevier.com/locate/ijplas The equivalent axisymmetric
More informationMonte Carlo Simulation of Ferroelectric Domain Structure: Electrostatic and Elastic Strain Energy Contributions
Monte Carlo Simulation of Ferroelectric Domain Structure: Electrostatic and Elastic Strain Energy Contributions B.G. Potter, Jr., B.A. Tuttle, and V. Tikare Sandia National Laboratories Albuquerque, NM
More informationAdvanced Structural Analysis EGF Section Properties and Bending
Advanced Structural Analysis EGF316 3. Section Properties and Bending 3.1 Loads in beams When we analyse beams, we need to consider various types of loads acting on them, for example, axial forces, shear
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 informationAnalytical Methods for Materials
Analytical Methods for Materials Lesson 15 Reciprocal Lattices and Their Roles in Diffraction Studies Suggested Reading Chs. 2 and 6 in Tilley, Crystals and Crystal Structures, Wiley (2006) Ch. 6 M. DeGraef
More information548 Advances of Computational Mechanics in Australia
Applied Mechanics and Materials Online: 2016-07-25 ISSN: 1662-7482, Vol. 846, pp 547-552 doi:10.4028/www.scientific.net/amm.846.547 2016 Trans Tech Publications, Switzerland Geometric bounds for buckling-induced
More informationAbstract. 1 Introduction
Contact analysis for the modelling of anchors in concrete structures H. Walter*, L. Baillet** & M. Brunet* *Laboratoire de Mecanique des Solides **Laboratoire de Mecanique des Contacts-CNRS UMR 5514 Institut
More information3.012 Structure An Introduction to X-ray Diffraction
3.012 Structure An Introduction to X-ray Diffraction This handout summarizes some topics that are important for understanding x-ray diffraction. The following references provide a thorough explanation
More informationTHE TOPOLOGICAL DESIGN OF MATERIALS WITH SPECIFIED THERMAL EXPANSION USING A LEVEL SET-BASED PARAMETERIZATION METHOD
11th. World Congress on Computational Mechanics (WCCM XI) 5th. European Conference on Computational Mechanics (ECCM V) 6th. European Conference on Computational Fluid Dynamics (ECFD VI) July 20-25, 2014,
More informationVibration damping in polygonal plates using the acoustic black hole effect: model based on the image source method
Vibration damping in polygonal plates using the acoustic black hole effect: model based on the image source method Jacques Cuenca a,b, Adrien Pelat a, François Gautier a a. Laboratoire d Acoustique de
More informationChapter 6: Plastic Theory
OHP Mechanical Properties of Materials Chapter 6: Plastic Theory Prof. Wenjea J. Tseng 曾文甲 Department of Materials Engineering National Chung Hsing University wenjea@dragon.nchu.edu.tw Reference: W. F.
More informationComputational Modelling of Vibrations Transmission Loss of Auxetic Lattice Structure
Vibrations in Physical Systems Vol. 27 (2016) Abstract Computational Modelling of Vibrations Transmission Loss of Auxetic Lattice Structure Eligiusz IDCZAK Institute of Applied Mechanics, Poznan University
More informationActuation of kagome lattice structures
Actuation of kagome lattice structures A.C.H. Leung D.D. Symons and S.D. Guest Department of Engineering, University of Cambridge, Trumpington Street, Cambridge, CB2 1PZ, UK The kagome lattice has been
More informationBending Load & Calibration Module
Bending Load & Calibration Module Objectives After completing this module, students shall be able to: 1) Conduct laboratory work to validate beam bending stress equations. 2) Develop an understanding of
More informationProf. B V S Viswanadham, Department of Civil Engineering, IIT Bombay
51 Module 4: Lecture 2 on Stress-strain relationship and Shear strength of soils Contents Stress state, Mohr s circle analysis and Pole, Principal stressspace, Stress pathsin p-q space; Mohr-coulomb failure
More informationENGINEERING MECHANICS
ENGINEERING MECHANICS Engineering Mechanics Volume 2: Stresses, Strains, Displacements by C. HARTSUIJKER Delft University of Technology, Delft, The Netherlands and J.W. WELLEMAN Delft University of Technology,
More informationDr.Vinod Hosur, Professor, Civil Engg.Dept., Gogte Institute of Technology, Belgaum
STRUCTURAL DYNAMICS Dr.Vinod Hosur, Professor, Civil Engg.Dept., Gogte Institute of Technology, Belgaum Overview of Structural Dynamics Structure Members, joints, strength, stiffness, ductility Structure
More informationFinite Element Analysis Lecture 1. Dr./ Ahmed Nagib
Finite Element Analysis Lecture 1 Dr./ Ahmed Nagib April 30, 2016 Research and Development Mathematical Model Mathematical Model Mathematical Model Finite Element Analysis The linear equation of motion
More informationViscoelastic Damping Characteristics of Indium-Tin/SiC Particulate Composites
Viscoelastic Damping Characteristics of Indium-Tin/SiC Particulate Composites HyungJoo Kim, Colby C. Swan Ctr. for Computer-Aided Design, Univ. of Iowa Roderic Lakes Engineering Physics, Univ. of Wisconsin
More informationPlane Strain Test for Metal Sheet Characterization
Plane Strain Test for Metal Sheet Characterization Paulo Flores 1, Felix Bonnet 2 and Anne-Marie Habraken 3 1 DIM, University of Concepción, Edmundo Larenas 270, Concepción, Chile 2 ENS - Cachan, Avenue
More informationCavitation instability in rubber with consideration of failure
JOURNAL OF MATERIALS SCIENCE 36 (2001)1901 1909 Cavitation instability in rubber with consideration of failure W. J. CHANG, J. PAN Mechanical Engineering and Applied Mechanics, The University of Michigan,
More informationGLOBAL AND LOCAL LINEAR BUCKLING BEHAVIOR OF A CHIRAL CELLULAR STRUCTURE
GLOBAL AND LOCAL LINEAR BUCKLING BEHAVIOR OF A CHIRAL CELLULAR STRUCTURE Alessandro Spadoni, Massimo Ruzzene School of Aerospace Engineering Georgia Institute of Technology Atlanta, GA 30332 Fabrizio Scarpa
More informationMicromechanical analysis of FRP hybrid composite lamina for in-plane transverse loading
Indian Journal of Engineering & Materials Sciences Vol. 15, October 2008, pp. 382-390 Micromechanical analysis of FRP hybrid composite lamina for in-plane transverse loading K Sivaji Babu a *, K Mohana
More informationChapter 6: Mechanical Properties of Metals. Dr. Feras Fraige
Chapter 6: Mechanical Properties of Metals Dr. Feras Fraige Stress and Strain Tension Compression Shear Torsion Elastic deformation Plastic Deformation Yield Strength Tensile Strength Ductility Toughness
More informationME 176 Final Exam, Fall 1997
Tuesday, December 16, 5:00 8:00 PM, 1997. Answer all questions for a maximum of 100 points. Please write all answers in the space provided. If you need additional space, write on the back sides. Indicate
More informationModule 5: Failure Criteria of Rock and Rock masses. Contents Hydrostatic compression Deviatoric compression
FAILURE CRITERIA OF ROCK AND ROCK MASSES Contents 5.1 Failure in rocks 5.1.1 Hydrostatic compression 5.1.2 Deviatoric compression 5.1.3 Effect of confining pressure 5.2 Failure modes in rocks 5.3 Complete
More informationFLEXIBILITY METHOD FOR INDETERMINATE FRAMES
UNIT - I FLEXIBILITY METHOD FOR INDETERMINATE FRAMES 1. What is meant by indeterminate structures? Structures that do not satisfy the conditions of equilibrium are called indeterminate structure. These
More informationVIBRATION ANALYSIS OF TIE-ROD/TIE-BOLT ROTORS USING FEM
VIBRATION ANALYSIS OF TIE-ROD/TIE-BOLT ROTORS USING FEM J. E. Jam, F. Meisami Composite Materials and Technology Center Tehran, IRAN jejaam@gmail.com N. G. Nia Iran Polymer & Petrochemical Institute, Tehran,
More informationPoisson s ratio effect of slope stability calculations
Poisson s ratio effect of slope stability calculations Murray Fredlund, & Robert Thode SoilVision Systems Ltd., Saskatoon, SK, Canada ABSTRACT This paper presents the results of a study on the effect of
More informationThe Frictional Regime
The Frictional Regime Processes in Structural Geology & Tectonics Ben van der Pluijm WW Norton+Authors, unless noted otherwise 1/25/2016 10:08 AM We Discuss The Frictional Regime Processes of Brittle Deformation
More information