Cellular solid structures with unbounded thermal expansion. Roderic Lakes. Journal of Materials Science Letters, 15, (1996).

Similar documents
Piezoelectric composite lattices with high sensitivity

Properties of a chiral honeycomb with a Poisson's ratio -1

Three-dimensional stiff cellular structures with negative Poisson's ratio

Analysis of high loss viscoelastic composites

Lecture 15 Strain and stress in beams

Lecture 4 Honeycombs Notes, 3.054

STRESSES WITHIN CURVED LAMINATED BEAMS OF DOUGLAS-FIR

Bending of Simply Supported Isotropic and Composite Laminate Plates

A short review of continuum mechanics

Curvature of a Cantilever Beam Subjected to an Equi-Biaxial Bending Moment. P. Krulevitch G. C. Johnson

TABLE OF CONTENTS. Mechanics of Composite Materials, Second Edition Autar K Kaw University of South Florida, Tampa, USA

Fig. 1. Circular fiber and interphase between the fiber and the matrix.

Module 7: Micromechanics Lecture 34: Self Consistent, Mori -Tanaka and Halpin -Tsai Models. Introduction. The Lecture Contains. Self Consistent Method

Choi, J. B. and Lakes, R. S., "Nonlinear properties of polymer cellular materials with a negative Poisson's. ratio",

ELECTROMECHANICAL RESPONSE OF PIEZOELECTRIC FOAMS

Laminated Beams of Isotropic or Orthotropic Materials Subjected to Temperature Change

Macroscopic Elastic Constitutive Relationship of Cast-in-Place Hollow-Core Slabs

ADVANCES IN NEGATIVE POISSON'S RATIO MATERIALS. Roderic Lakes adapted from Advanced Materials, 5, , (1993). Diagrams at bottom.

Large Thermal Deflections of a Simple Supported Beam with Temperature-Dependent Physical Properties

ENGN 2340 Final Project Report. Optimization of Mechanical Isotropy of Soft Network Material

Module 7: Micromechanics Lecture 29: Background of Concentric Cylinder Assemblage Model. Introduction. The Lecture Contains

Measurement of deformation. Measurement of elastic force. Constitutive law. Finite element method

Stresses in Curved Beam

Materials: engineering, science, processing and design, 2nd edition Copyright (c)2010 Michael Ashby, Hugh Shercliff, David Cebon.

Composite materials: mechanical properties

3D Elasticity Theory

Stability of Cosserat solids: size effects, ellipticity and waves

Mechanics of Irregular Honeycomb Structures

Design of a fastener based on negative Poisson's ratio foam adapted from

EE C247B / ME C218 INTRODUCTION TO MEMS DESIGN SPRING 2014 C. Nguyen PROBLEM SET #4

Mechanics PhD Preliminary Spring 2017

PURE 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.

The fracture toughness of planar lattices: imperfection sensitivity

Presented By: EAS 6939 Aerospace Structural Composites

Advanced Structural Analysis EGF Section Properties and Bending

16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations

Fundamentals of Linear Elasticity

PROPAGATION OF CURVED CRACKS IN HOMOGENEOUS AND GRADED MATERIALS

Computational Analysis for Composites

Lecture 7, Foams, 3.054

Reduced warp in torsion of reticulated foam due to Cosserat elasticity: experiment

Strain Gages. Approximate Elastic Constants (from University Physics, Sears Zemansky, and Young, Reading, MA, Shear Modulus, (S) N/m 2

The Rotating Inhomogeneous Elastic Cylinders of. Variable-Thickness and Density

CHEM-C2410: Materials Science from Microstructures to Properties Composites: basic principles

7.4 The Elementary Beam Theory

A Note on Suhir s Solution of Thermal Stresses for a Die-Substrate Assembly

LAMINATED COMPOSITE PLATES

Chapter 2: Elasticity

Mechanical Engineering Ph.D. Preliminary Qualifying Examination Solid Mechanics February 25, 2002

Forced Response of Plate with Viscoelastic Auxetic Dampers

Supplementary Figures

Composite Structures. Indian Institute of Technology Kanpur

Prediction of Elastic Constants on 3D Four-directional Braided

THE INFLUENCE OF IN-PLANE DENSITY VARIATION ON ENGINEERING PROPERTIES OF ORIENTED STRANDBOARD: A FINITE ELEMENT SIMULATION

MATERIAL ELASTIC ANISOTROPIC command

6. Bending CHAPTER OBJECTIVES

Buckling Behavior of 3D Randomly Oriented CNT Reinforced Nanocomposite Plate

CHAPTER -6- BENDING Part -1-

BENCHMARK LINEAR FINITE ELEMENT ANALYSIS OF LATERALLY LOADED SINGLE PILE USING OPENSEES & COMPARISON WITH ANALYTICAL SOLUTION

Effect of Growth Direction on Twin Formation in GaAs Crystals Grown by the Vertical Gradient Freeze Method

Characterisation and Modelling of a Melt-Extruded LDPE Closed Cell Foam

The aims of this experiment were to obtain values for Young s modulus and Poisson s ratio for

Example 3.7 Consider the undeformed configuration of a solid as shown in Figure 3.60.

Keywords: Re-entrant; cellular solid; Negative Poisson s ratio; Voronoi; Energy absorption; Compression

USING A HOMOGENIZATION PROCEDURE FOR PREDICTION OF MATERIAL PROPERTIES AND THE IMPACT RESPONSE OF UNIDIRECTIONAL COMPOSITE

Theories of Straight Beams

Lecture 7: The Beam Element Equations.

Supplementary figures

4. BEAMS: CURVED, COMPOSITE, UNSYMMETRICAL

Basic Equations of Elasticity

PEAT SEISMOLOGY Lecture 2: Continuum mechanics

Elasticity in two dimensions 1

Frequently Asked Questions

Regular, low density cellular structures - rapid prototyping, numerical simulation, mechanical testing

DYNAMIC RESPONSE OF SYNTACTIC FOAM CORE SANDWICH USING A MULTIPLE SCALES BASED ASYMPTOTIC METHOD

The New Boundary Condition Effect on The Free Vibration Analysis of Micro-beams Based on The Modified Couple Stress Theory

MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MATERIALS SCIENCE AND ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139

MICROMECHANICAL ANALYSIS OF FRP COMPOSITES SUBJECTED TO LONGITUDINAL LOADING

Consider an elastic spring as shown in the Fig.2.4. When the spring is slowly

UNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich

Strong Cosserat elastic effects in a unidirectional composite

CHAPTER OBJECTIVES Use various methods to determine the deflection and slope at specific pts on beams and shafts: 2. Discontinuity functions

THE TOPOLOGICAL DESIGN OF MATERIALS WITH SPECIFIED THERMAL EXPANSION USING A LEVEL SET-BASED PARAMETERIZATION METHOD

Viscoelastic Properties of Foamed Thermoplastic Vulcanizates and their Dependence on Void Fraction

MODIFIED HYPERBOLIC SHEAR DEFORMATION THEORY FOR STATIC FLEXURE ANALYSIS OF THICK ISOTROPIC BEAM

CHAPTER 5. Beam Theory

CHEM-E2200: Polymer blends and composites Fibre architecture and principles of reinforcement

Elements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004

Lecture 8. Stress Strain in Multi-dimension

On the visco-elastic properties of open-cell polyurethane foams in uniaxial compression

Strain Gages. Approximate Elastic Constants (from University Physics, Sears Zemansky, and Young, Reading, MA, 1979

3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1

FOAM TOPOLOGY BENDING VERSUS STRETCHING DOMINATED ARCHITECTURES

INFLUENCE OF TEMPERATURE ON BEHAVIOR OF THE INTERFACIAL CRACK BETWEEN THE TWO LAYERS

Mechanical and Thermal Properties of Coir Fiber Reinforced Epoxy Composites Using a Micromechanical Approach

International Journal of Scientific & Engineering Research Volume 8, Issue 12, December-2017 ISSN

ASSESSMENT OF MIXED UNIFORM BOUNDARY CONDITIONS FOR PREDICTING THE MACROSCOPIC MECHANICAL BEHAVIOR OF COMPOSITE MATERIALS

Module 7: Micromechanics Lecture 25: Strength of Materials Approach. Introduction. The Lecture Contains. Effective Transverse Modulus

International Journal of Pure and Applied Sciences and Technology

Viscoelastic Damping Characteristics of Indium-Tin/SiC Particulate Composites

Transcription:

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 of thermal expansion larger than that of either constituent. Thermal expansion of materials is of considerable practical interest since materials in service may experience a considerable range of operating temperature. Materials with microstructure, including composites, honeycombs, and foams (which are composites with one void phase) are finding increasing application for many purposes. It is therefore of interest to explore the range of thermal expansion coefficient attainable in materials with microstructure. The constitutive equation for a linear isotropic thermoelastic continuum is [1]: ε ij = 1 + ν E σ ij - [ ν E σ ii - α T]δ ij, (1) in which ε ij is strain, σ ij is stress, ν is Poisson's ratio, E is Young's modulus, α is the thermal expansion coefficient, δ ij is Kronecker's delta, and the usual index notation is used. Thermal expansion of crystalline solids is considered to arise from the anharmonicity of the interatomic potential. Thermal expansion coefficients for solids are on the order of 3 x 10-6 / C for ceramics, 10-5 / C for metals, and 10-4 / C for polymers [2]; many specific values are given in reference [3]. Thermal expansion of foams was discussed briefly by Gibson and Ashby [2] who pointed out that for a foam viewed as a framework (of homogeneous ribs), the thermal expansion coefficient is the same as the solid from which it is made. An exception is a closed cell foam of very small solid volume fraction, for which thermal expansion is enhanced by the expansion of the gas within the cells [2]. An upper bound for this effect can be estimated as follows. If the foam is very compliant, expansion of the gas occurs at approximately constant pressure. Considering the ideal gas law PV = nrt, the expansion of such a compliant foam at room temperature (T = 300 K) can be no more than α = 1.1 x 10-3 / C, greater than that of a solid polymer. For arbitrary two-phase composite materials, bounds have been developed for the thermal expansion coefficient α of the composite in terms of that of each constituent [4,5]; the upper bound is a rule of mixtures. In deriving these bounds it was assumed that there are only two phases, and there is no empty space in the composite. In the present work, it is shown that if one relaxes the

2 assumption of no empty space, it is possible to generate cellular micro-structures which exhibit a thermal expansion coefficient much larger than that of either constituent. Consider a cellular solid consisting of bi-layer ribs. The rib element in the cellular solid consists of a laminate of two layers of differing thermal expansion coefficient α. The layers have Young's moduli E 1 and E 2, Poisson's ratio ν 1 and ν 2 ; thermal expansion coefficient α 1 and α 2 ; and thicknesses h 1 and h 2. The curvature κ (the inverse of the radius of curvature ρ) of a thin bilayer elastic beam is given in terms of the temperature change T by Timoshenko [6], neglecting end effects [7]. The curvature may be written κ = 6(α 2 - α 1 ) T (1 + h 1 h 2 ) 2 (h 1 + h 2 ) 1 3((1 + h 1 h ) 2 ) + (1 + h 1 E 1 2 h 2 E ) (( h 1 2 h ) 2 + h 2E 2 2 h 1 E ) 1. (2) The curvature can be made arbitrarily large by making the strip thickness h small. The ratio of Young's moduli has little effect on the curvature [6]. We determine the thermal expansion coefficient of cellular solids with two-layer rib elements. Consider first the deformations and strains caused by rib curvature due to strain. The transverse displacement δ at the midpoint of an initially straight rib of length l is [6] δ = κl2 8. (3) Cellular solids which exhibit global thermal expansion due to such transverse rib displacements must have an unusual connectivity, such as that shown in Fig. 1. It is easier to make use of longitudinal displacement of rib elements as determined as follows. Referring to Fig. 2, consider now a curved rib element of arc length l arc, radius of curvature r, and included angle θ, in radians. The rib length is given by l = l arc 2 θ sin θ 2. (4) The strain of a unit cell due to end to end length changes in the rib is ε = dl l. (5) The effect of a curvature increment dκ due to bending is ε = [ 1 2 ctn θ 2-1 θ ] l arcdκ. (6)

3 Eq. 2 gives the curvature of a straight rib, or the increment in curvature of an initially curved rib. Eq. 6 shows that there must be some initial curvature or no finite length change occurs. Cellular solids which expand due to longitudinal displacement by bending of ribs with initial curvature are shown in Fig. 3. The thermal expansion coefficient of the cellular solid structure is α = ε/ T (at constant stress, see Eq. 1), so α = l arc (h 1 + h 2 ) 6(α 2 - α 1 )(1 + h 1 h 2 ) 2 3((1 + h 1 h ) 2 ) + (1 + h 1 E 1 2 h 2 E ) (( h 1 2 h ) 2 + h 2E 2 2 h 1 E ) 1 [ 1 2 ctn θ 2-1 θ ] (7) The ratio of rib arc length to rib total thickness can be made arbitrarily large, by making the ribs very slender. The initial curvature of the ribs, as quantified by the included angle θ, should also be made large in order to maximize the magnitude of α given by Eq. 7 (Fig. 4). The thermal expansion coefficient of the cellular solid structure is therefore unbounded. It can be made much larger in magnitude than the thermal expansion coefficient of either solid constituent. The sign of the thermal expansion coefficient of the cellular solid is governed by the placement of the constituents within each rib. Referring to Fig. 3, if the constituent with the higher thermal expansion coefficient is on its concave side, an increase in temperature will cause the rib to straighten, giving rise to a positive thermal expansion coefficient for the cellular solid. Conversely, if the constituent with the higher thermal expansion coefficient is on its convex side, an increase in temperature will cause the rib to curve more and become shorter, giving rise to a negative thermal expansion coefficient for the cellular solid. The structure shown in Fig. 3 has two dimensional cubic symmetry giving rise to anisotropic elasticity but isotropic thermal expansion. Hexagonal structures which are elastically isotropic in plane are also possible. The examples presented are two-dimensional honeycombs. They could be produced by a co-extrusion process, by progressive lamination, or by lithography. One may also envisage threedimensional cellular solids, or foams, containing bi-layer ribs, which may be sintered together. Materials based on polymer bi-layers can deliver higher thermal expansion than those based on metallic bi-layers. These structures provide a further illustration of the importance of void space in determining physical properties. Inclusion of void space permits one to substantially exceed bounds for thermal expansion developed under the assumption of no void space. Control of the shape of void space has led to materials with a negative Poisson's ratio [8] and lightweight hierarchical materials with high ratios of strength to weight [9]. More developments in materials with unusual or extreme properties are anticipated.

4 To conclude, arbitrarily high values of thermal expansion can be attained by preparing cellular solids with ribs consisting of a bi-layer containing materials of two different expansion coefficients. References 1 I. S. Sokolnikoff, Mathematical theory of elasticity, (Krieger, Malabar, FL, 2nd ed., 1983). 2 L. J. Gibson and M.F. Ashby, Cellular Solids (Pergamon press, Oxford, 1987). 3 Staff, "Comparison of materials; coefficient of thermal expansion", Materials Engineering 103 (1986) 25-7 4 J. L. Cribb, "Shrinkage and thermal expansion of a two-phase material", Nature, 220 (1963) 576-577. 5 S. Torquato, "Bounds on the thermoelastic properties of suspensions of spheres", J. Appl. Phys. 67, (1990), 7223-7227. 6 S. P. Timoshenko, "Analysis of bi-metal thermostats", J. Optical Soc. America, 11 (1925) 233-355. 7 E. Suhir, "Stresses in bi-metal thermostats", J. Appl. Mech., 53 (1986) 657-660. 8 R.S. Lakes, "Foam structures with a negative Poisson's ratio", Science, 235 (1987) 1038-1040. 9 R.S. Lakes, "Materials with structural hierarchy", Nature, 361, (1993) 511-515.

5 Figures 1 Cellular solid which undergoes thermal expansion via lateral bending displacement of ribs. 2 Rib element geometry with initial curvature. l arc = θ l l/2 θ r θ /2 r

6 3 Cellular solid with curved ribs. It undergoes thermal expansion via longitudinal bending displacement of ribs.

7 4 Normalized thermal expansion and rib aspect ratio l arc h α α 2 - α 1 of cellular solid vs. rib curve angle θ (in radians) for h 1 = h 2 = h and E 1 = E 2.