arxiv: v2 [cond-mat.mtrl-sci] 28 Sep 2014

Transcription

1 Mechanics of elastic networks Andrew N. Norris Mechanical and Aerospace Engineering, Rutgers University, Piscataway, NJ , USA arxiv:49.848v [cond-mat.mtrl-sci] 8 Sep 4 Abstract We consider a periodic lattice structure in d = or 3 dimensions with unit cell comprising Z thin elastic members emanating from a similarly situated central node. A general theoretical approach provides an algebraic formula for the effective elasticity of such frameworks. The method yields the effective cubic elastic constants for 3D space-filling lattices with Z = 4, 6, 8, and 4, the latter being the stiffest lattice proposed by []. The analytical expressions provide explicit formulas for the effective properties of pentamode materials, both isotropic and anisotropic, obtained from the general formulation in the stretch dominated limit for Z = d+.. Introduction Space frames, or periodic lattice structures of rods and joints, have long been of interest to engineers, architects, materials scientists and others. The octet truss, for instance, which is common in modern large scale structures because of its load bearing capacity may be attributed to Alexander Graham Bell s interest in tetrahedral cells for building man-carrying kites []. Recent fabrication of micro-architectured materials have used the octet truss tetrahedral cell design to achieve ultralight and ultrastiff structures [3]. An even stiffer structure comprising tetrakaidecahedral unit cells was proposed by [4, ], see Fig.. Unlike the octet truss which has cubic elastic symmetry, the tetrakaidecahedral structure can display isotropic effective elastic properties. At the other end of the stiffness spectrum for elastic lattice structures are pentamode materials PMs with five easy modes of deformation [5] see also page 666 of [6]. The range of such material properties, including high stiffness, strength and fracture toughness, exhibited by low density micro-architectured materials is reviewed in [7]. The response of low density lattice structures depends on whether the deformation under load is dominated by stretching versus bending. This in turn depends upon the coordination number, Z, the number of nearest neighboring joints in the unit cell, see Fig. for several examples ranging from pentamodal Z = 4 to stiffest Z = 4. Maxwell [8] described the necessary although not sufficient condition for a d-dimensional d =,3 space frame of b struts and j pin joints to be just rigid: b 3 = j 3d. For an infinite periodic structure, b jz/, Maxwell s condition becomes Z = d. Structures with Z = d, known as isostatic lattices, are at the threshold of mechanical stability [9]. A closer examination of the issue taking into account the degrees of freedom in the applied strain field, dd+/, leads to the conclusion that the necessary and sufficient condition for rigidity of frameworks with coordination number Z is Z dd + []. The octet truss lattice Z = is an example of a 3D lattice which satisfies the rigidity condition []. Three dimensional frameworks with Z < admit soft modes; thus, as we will see in 44.6, a cubic framework with Z = 6 has 3 soft modes. Zero frequency modes, floppy modes, that occur for Z < d correspond to collapse mechanisms, a topic also examined by [] for truss-like D lattices. Three dimensional elasticity is characterized by 6 positive eigenvalues [3]. A pentamode material PM in 3D is the special case of elasticity with five zero eigenvalues, hence penta. An inviscid compressible fluid like water serves as a useful reference material for PMs since it has a single bulk modulus but zero shear rigidity, the elastic stiffness tensor is C = K I I C ijkl = K δ ij δ kl, where K is the bulk modulus and C ijkl are the components of C. This form of the elastic moduli corresponds to a rank one 6 6 matrix Voigt matrix [C IJ ] with single non-zero eigenvalue 3K. PMs can therefore be thought of as

2 elastic generalizations of water but without the ability to flow; however, unlike water, for which the stress is isotropic, PMs can display anisotropy. Recent interest in PMs has increased after the observation that they provide the potential for realizing transformation acoustics [4]. Pentamode materials can be realized from specific microstructures with tetrahedral-like unit cells [5, 5]. These types of PM lattice structures are related to low density materials such as foams in which the low density is a consequence of the low filling fraction of the solid phase, see [6] for a review of mechanical properties of low density materials. Here we consider specific microstructures and find explicit values for the elastic moduli for isotropic and anisotropic PMs. Thepurposeofthe presentpaperistwo-fold. First, wefillthe need forageneraltheoreticalapproachthat provides a simple means to estimate the effective elasticity of frameworks with nodes which are all similarly situated. Nodes are similarly situated if the framework appears the same when viewed from any one of the nodes [7]; the unit cell must therefore be space-filling, as are the cases in Fig.. Specific homogenization methods have been proposed for lattice structures; e.g. [8] use a mix of analytical and finite element methods, while [9, ] provide a general mathematical scheme that is not easy to implement in practice. More general micropolar elasticity theories have also been considered for two-dimensional frameworks, e.g. by applying force and moment balances on the unit cell [], or alternatively, using energy based methods []. The method proposed here derives the elastic tensor relating the symmetric stress to the strain. It does not assume micropolar theories, although the solution involves a local rotation within the unit cell required for balancing the moments, see 3. In contrast to prior works, the present method is explicit and practical; it provides for instance, the effective cubic elastic constants for all the examples in Fig., see The second objective is to provide analytical expressions for the effective properties of pentamode materials, both isotropic and anisotropic. The general theory derived here is perfectly suited to this goal. We show in 44. that the minimal coordination number necessary for a fully positive definite elasticity tensor is Z = d + d = or 3, the pentamode limit therefore follows by taking the stretch dominated limit for Z = d+. The paper proceeds as follows: The lattice model is introduced and the main results for the effective properties are summarized in. The detailed derivation is presented in 3. In 4 some properties of the effective moduli are described, including the stretch dominated limit, and examples of 5 different lattice structures are given. Pentamode materials, which arise as a special case of the stretch dominated limit when the coordination number is d +, are discussed in 5. The two dimensional case is presented in 6 and conclusions are given in 7.. Lattice model The structural unit cell in d-dimensions d =,3 comprises Z d + rods and has volume V. Let denote the position of the single junction in the unit cell with the cell edges at the midpoints of the rods, located at R i for i =,...,Z. Under the action of a static loading the relative position of the vertex initially located at R i moves to r i. The angle between members i and j before and after deformation is Ψ ij = cos R i R j /R i R j and ψ ij = cos r i r j /r i r j, respectively, where R i = R i, r i = r i. The end displacement r i = r i R i is decomposed as r i = r i + r i. In the linear approximation assumed here r i r ie i where r i = r i R i and the unit axial vector is e i = R i /R i e i =. The change in angle between members i and j is ψ ij ψ ij Ψ ij, j i. The transverse displacement r i can include a contribution r rot i e i r rot i = caused by rigid body rotation of the unit cell. We therefore define r b i = r i r rot i, the transverse displacement associated with flexural bending. Vectors perpendicular to e i are used to define transverse bending forces: The unit vector e ij is perpendicular to e i and lies in the plane spanned by e i and e j with e ij e j <, that is e ij = e i e i e j ei e j = cosψ ij e i e j /sinψ ij, i j Z. The unit vectors e α i, α = : d, are such that {e i, e α i } form an orthonormalset of d-vectors. Summation on repeated lower case Greek superscripts is understood and only relevant for 3D. Define P i = e i e i, P i = e α i e α i, If e i = e j we consider a slight perturbation so that ψ ij π.

3 so that P i + P i = I, the unit matrix in d-dimensions. The axial tensor of a vector v is defined by its action on a vector w as axvw = v w. Finally, although the derivation will be mostly coordinate free, for the purpose of defining examples and the components of the effective stiffness tensor, we will use the orthonormal basis {a q } q = : d. a Z = 4 b Z = 6 c Z = 8 d Z = e Z = 4 Figure : Unit cell for some lattices considered, see Table. The stretch dominated Z = 4 lattice with node at the center of the tetrakaidecahedral unit cell has maximal stiffness []. Colours are used to illustrate the structure for Z = and to differentiate the shorter blue members from the longer onesgreen for Z = 4... Forces on individual members The members interact in the static limit via combined axial forces directed along the members, and bending moments, associated with axial deformation and transverse flexure, respectively. We also include the possibility of nodal bending stiffness, associated with torsional spring effects at the junction The strain energy can then be represented [] H = H s +H b +H n for stretch, bending and nodal deformation, respectively. Later, we consider the limit in which the contributions from bending, H b and H n, are small, and the deformation may be approximated by axial forces only, 3

4 which is the stretch-dominated limit. Physically, this corresponds to slender members with small thickness to length ratio. We assume strain energy of the form H s = r i M i, H b = r b i N i, H n = R i R j N ij ψ ij 3 where M i are the axial compliances, N i the bending compliances and N ij arethe nodal bending compliances. The force acting at the end of member i i =,...,Z is f i = f s i +fb i +fn i 4 where f s i = r i /M i, acting parallel to the member, is associated with stretching. The perpendicular component of the force acting on the member s end is comprised of a shear force f b i = r b i /N i caused by the bending of the member, plus a shear force f n i = R j ψ ij e ij /N ij associated with the node compliance. The axial and bending compliances can be related to the member properties via M i = Ri dx Ri x dx, N i =, i Z 5 E i A i E i I i where E i x, A i x, I i x are the Young s modulus, cross-sectional area and moment of inertia, with x = at the nodal junction. We assume circular or square cross-section in 3D so that only a single bending compliance is required for each member, otherwise the results below involving N i are not generally valid although they could be amended with necessary analytical complication. The nodal bending compliances N ij are arbitrary and satisfy the symmetry N ij = N ji which ensures that the sum of the moments of the node bending forces are zero... Effective stress and moduli We consider the forces on the members of the unit cell responding to an applied macroscopic loading. The forces acting at the node of the unit cell are equilibrated, as are the moments, f i =, R i f i =. 6 Treating the volume of the cell as a continuum with equilibrated stress σ, integrating divx σ = σ over V and identifying the tractions as the point forces f i acting on the cell boundary, implies the well-known connection Z σ = V R i f i. 7 The symmetry of the stress, σ = σ T, is guaranteed by the moment balance 6. Our aim is to derive the effective elastic moduli defined by the fourth order tensor C which relates the stress to the macroscopic strain ǫ according to σ = Cǫ. 8 The elements of the elastic stiffness C when expressed in an orthonormal basis possess the symmetries C ijkl = C jikl and C ijkl = C klij, and the elements can also be represented in terms of the Voigt notation via C ijkl C IJ = C JI. 4

5 .3. Summary of the main result for the effective elastic stiffness We first introduce the vectors d i, d α i, d ij = d ji, the second order symmetric tensors D i, D α i, D ij = D ji and the L L matrix with elements P ij, where L = Zd+ZZ /: d i = e i Mi, d α i = eα i, α = : d d ij = Ni D i = R ip i, D α i = R i VMi VNi e i e α i +eα i e i, D ij = R i R j eij + e ji, 9a N ij R i R j R i R j ei e ij +e j e ji 9b VN ij {u k } L k= = {d i, d α i, d ij }, {U k } L k= = {D i, D α i, D ij }, α = : d 9c L P ij = δ ij u i u k u k uj, i, j = : L, k= then, under some general assumptions applicable to the 3D structures in Fig., eq., the effective moduli can be be written L C = P ij U i U j. i,j= These results are derived in the next section and implications are discussed in 4, including a simple expression 8 for the elastic moduli represented in 6 6 Voigt notation. The general structure of eqs. 9 holds for d = without requiring the zero rotation conditions of eq., as discussed in 6. 9d 3. Derivation of the effective elasticity tensor 3.. Affine deformation Strain is introduced through the affine kinematic assumption that the effect of deformation is to cause the cell edges to displace in a linear manner proportional to the local deformation gradient F. Edge points originally located at R i are translated to FR i. In addition to the affine motion, we include two d-vectors, introduced to satisfy the equilibrium conditions 6. Following [3] we assume that the junction moves from the origin to χ. An additional rotation Q SOd is introduced, so that the vector defining the edge relative to the vertex is r i = QFR i χ. The linear approximation for the deformation is F = I + ǫ + ω with ǫ = ǫ T and ω = ω T. We take Q = e Γ = I + Γ+OΓ where the skew symmetric matrix Γ is defined by the d-vector γ as Γ = axγ. Hence, r i r i R i = ǫ+ω +ΓR i χ. In the linear approximationr i can equally well be taken along R i as far as second order terms are concerned. Thus, r i = R i P i : ǫ e i χ e i, r i = P i R i ǫe i χ+r i ω +Γe i, ψ ij = e i ǫe ij +e j ǫe ji 3 R i e ij +R j e ji χ. The tangential displacement governing the shear bending force is, after removing the affine rigid body rotation, r b i = r i R i ωe i. 4 5

6 Note that we retain the unknown rotation Γ in order to satisfy the moment equilibrium condition 6. Hence, in the linear approximation 4 becomes f i = M i R i P i : ǫe i +N i R i P i ǫe i +γ e i + N ij R j ei ǫe ij +e j ǫe ji eij M i P i +N i P i + N ij R je ij R i e ij +R j e ji χ. 5 This explicit expression for the forces allows us to determine the vectors χ and γ, next. 3.. Solution of the equilibrium equations Consider first the moment balance condition 6. Of the three terms comprising the force in eq. 4 only the bending shear forces fi b does not automatically yield zero moment. Equilibrium of the moments therefore reduces to R i fi b =. 6 Substituting f b i = rb i /N i and using eqs. 4 and 6 allows us to find γ in the form γ = B g χ j= The force on member i becomes, using eq. 5, where f i = R i M i P i : ǫe i + R i N i P i ǫe i +axe i B Rj Z R i, Z e j ǫe j where B = P R i i g = e i. 7 N j N i N i j= Rj R j axe j ǫe j + ei ǫe ij +e j ǫe ji eij N j N ij P i + P i + R i axe i Baxg+ R j e ij e ij + e ji χ. 8 M i N i N i N ij R i R j The equilibrium condition 6 can then be solved for χ as A = χ = A Z Ri M i P i : ǫe i + R i N i P i +axgbr i axe i ǫe i + R j 9 ei ǫe ij +e j ǫe ji eij N ij P i + P i + R i R j e ij e ij + e ji +axgbaxg. M i N i N ij R i R i R j Equations 7, 8 and 9 provide the desired linear relation between the strain and the stress from which one can derive the effective elastic moduli A simplification While eqs. 7, 8, 8- provide all of the necessary ingredients for the most general situation we assume for the remainder of the paper that the unit cell rotation vanishes, implying γ =. Hence, the vector g and the second term in the expression for γ in eq. 7 vanish. The latter identity is equivalent to 6

7 Dv v = v where D = Z R i N i e i e i. This implies that D must be proportional to the identity, hence the zero rotation condition may be written R i N i e i = and R i N i ei e i d I = zero cell rotation. The identities hold for the examples considered later. Note that the assumption of zero rotation is not necessary for stretch dominated lattices in which bending effects are negligible Effective stiffness In order to arrive at an explicit expression for the elastic stiffness tensor we first write the stress in terms of strain, using eqs. 7, 5-9, σ = V V Z k= R i P i M P i : ǫ+ R i e i e α i i N eα i ǫe i + i Z R i Mi P i d i + R i Ni e i e α i dα i dk R k Mk P i : ǫ+dα k + Z R k e α k ǫe k + Nk k= l k R i R j N ij e i e ij ei ǫe ij +e j ǫe ji R i R j N ij e i e ij d ij A d kl Rk R l N kl ek ǫe kl +e l ǫe lk. It follows from, the symmetry of the stress and strain and from the definition of the second order symmetric tensors D i, D ij in 9b that the elastic moduli can be expressed C = Di D i +D α i D α i Z Di d i +D α i d α i Z k= + + dk D k +d α k D α k + D ij D ij D ij d ij A d kl D kl. 3 k= l k Finally, we note, based on the definitions of the vectors in 9a, that A = di d i +d α i dα i + d ij d ij = L u i u i. 4 The sets {u k } and {U k } defined in 9c combine the Z vectors/tensors associated with stretch, the d Z vectors/tensors associated with shear, and the ZZ / vectors/tensors associated with nodal bending into sets of L = dz +ZZ / elements in terms of which 3 becomes C = L U i U i L L U i u i u j u j L u k U k. 5 It then follows from the definition of P in 9d that C can be expressed in the form. 7 j= k=

8 4. Properties of the effective moduli 4.. Generalized Kelvin form The L L symmetric matrix P with elements P ij defined in eq. 9d has the crucial properties P = P, rankp = L d, 6 i.e. P is a projector, and the dimension of its projection space is trp = L d. Hence, the summation in is essentially the sum of L d tensor products of second order tensors. This is to be compared with the Kelvin form for the elasticity tensor [3] C = 3d 3 λ i S i S i where λ i >, trs i S j = δ ij. 7 The second order symmetric tensors are eigenvectors {S i } that diagonalize the elasticity tensor, with eigenvectors λ i known as the Kelvin stiffnesses. Equation provides a non-diagonal representation for C. Note that L L s + L b + L n where L s = Z is associated with stretch, L b = d Z with bending shear and L n = ZZ / with nodal bending. A necessary although not sufficient condition for positive definiteness of C is that the rank of P, which is L d, exceed 3d 3. Ignoring nodal bending L = L s +L b this is satisfied if Z d + for d = and 3. The requirement is stricter in the stretch dominated limit L = L s : Z 6 in D and Z in 3D matrix in 3 dimensions The main result of eq. implies a simple representation for the 6 6 matrix of elastic moduli [C IJ ] based on the compact Voigt notation C ijkl C IJ in the orthonormalbasis {a, a, a 3 }. Let [u] 3 L denote the L vectors {u k } and let [U] 6 L denote the L second order tensors {U k } according to U Ik = a i U k a j with the standard correspondence I {,,3,4,5,6} ij {,,33,3,3,}. Then eq. becomes C C C 3 C 4 C 5 C 6 C C C 3 C 4 C 5 C 6 C 3 C 3 C 33 C 34 C 35 C 36 C 4 C 4 C 34 C 44 C 45 C 46 = [U][U] T [U][u] T [u][u] T [u][u] T. 8 C 5 C 5 C 35 C 45 C 55 C 56 C 6 C 6 C 36 C 46 C 56 C Bulk modulus If the effective medium has isotropic or cubic symmetry then a strain ǫ = εi produces strain σ = dkεi where K is the d dimensional bulk modulus. More generally, whether or not the symmetry is cubic or isotropic, we can define K = d C iijj. The bulk modulus follows from eqs. 9 and as K = d V i,j= This simplifies further under the broad assumption that certainly true of all the examples of Fig. considered in 44.6, so that K = d V P ij R i R j Mi M j. 9 e i M =, 3 8 R i M i. 3

9 Note that the bulk modulus depends only on the axial stiffness of the members. Assume the members are the same material E i = E, and each has constant cross-section area or width A i, then according to eqs. 5 and 3, K = φ d E where φ = V A i R i 3 is the volume fraction of solid material in the lattice. The scaling of bulk modulus with volume fraction, K φe, is well known, e.g. [4, eq..] for d =, [5], [6] for tetrakaidecahedral unit cells see below, and [] Model simplification While the model considered is quite general, in practice there is little information on the form of the nodal compliances for practical situations. For the remainder of the paper we concentrate on just the stretch and shear bending effects, so that L = L s +L b +L n L s +L b = dz. The stress-strain relation is then σ = [ R i X i I Ri Z Z ] X k X j R j : ǫ where X i = P i + P i. 33 M i N i k= j= A further simplifications is obtained by ignoring shear bending effects, i.e. L L s = Z, the stretch dominated limit, considered next Stretch dominated limit In this limit the forces f i have no transverse components. Physically, this corresponds to infinite bending compliances, /N i =, /N ij =, and may be achieved approximately by long slender members. By ignoring shear and nodal bending the expression for C reduces to C = V i,j= R i R j Mi M j P ij P i P j, P ij = δ ij e i Mi Z k= P k M k e j Mj. 34 It follows from eq. 6 that the Z Z projection matrix P with elements P ij has rank Z d Examples in 3D: Z = 4,6,8,,4 Allexamplesdisplaycubicsymmetry,withthreeindependentelasticmoduli: C, C andc 44. Introduce the fourth order tensors I, J and D with components I ijkl = δ ikδ jl +δ il δ jk, J ijkl = 3 δ ijδ kl, and D ijkl = δ i δ j δ k δ l +δ i δ j δ k δ l +δ i3 δ j3 δ k3 δ l3. A solid of cubic symmetry has elasticity of the form C = 3KJ+µ I D +µ D J. 35 The isotropic tensor J and the tensors of cubic symmetry I D and D J are positive definite [7], so the requirement of positive strain energy is that K, µ and µ are positive. These three parameters, called the principal elasticities by Kelvin [3], can be related to the standard Voigt stiffness notation: K = C +C /3, µ = C 44 and µ = C C /. The bulk modulus follows from eq. 3 as K = φ ZR E Z; K = 9 9VM 4a, Z 4; K = +, Z = V M M where for Z = 4 M, M arethe axial compliances of the two different types of members. It may be checked that K 4 = K 6 + K 8 where K Z denotes the bulk modulus for coordination number Z. The shear moduli are given in Table. Note that the effective compliance, relating strain to stress by ǫ = C σ is simply 9

10 Table : 3D lattice structures considered. They display cubic elastic symmetry with C = K µ, C = K 3 µ and C 44 = µ where K is given by eq. 36. All cases except Z = 4 have uniform rod length R and compliances M, N. The boundary of the tetrakaidecahedral Z = 4 unit cell has 36 edges each of length a and the cell comprises members of two types: 6 of length R = a and 8 of length R = 3 a, the average length of the members is R =.36a. The associated compliances are M, N and M, N. The volume fraction φ in all cases is based on cylindrical rods of uniform radius b. Note that the volume fraction increases with coordination number Z. Z cell V φ {e i } not normalized 64 4 diamond 3 3 R3. b R ± 6 simple 8R 3.36 b ± R cubic ± 8 BCC R3 4.8 b R FCC octet truss 4 tetrakaidecahedral 4 R b R 8 a b R ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± µ K 9M 4M+N 3M N + M N M 4N 3 M + N M + M µ K 3M N 3 3M N M 8N N + M + 3 N M + M C = 3K J+µ I D +µ D J. The ratio M/N may also be expressed in terms of the volume fraction φ since the rods are assumed to be solid circular so that M N = 3 4R. 37 Hence, Table indicates that µ =Oφ for Z = 4, 6 and µ =Oφ for Z = 4, 8; otherwise µ, µ =Oφ. Z = 4: the tetrakaidecahedral unit cell The tetrakaidecahedron is a truncated octahedron with all edges of the same length a V = 8 a 3. Rods extend from the centerto allfaces ofthe Kelvin cell as shownin Fig.. Note the functional dependence µ = µ M,N, µ = µ M,N,N. Isotropy µ = µ is achieved if in which case the effective Poisson s ratio is ν = b 3 M 3 N = M N 38 M N 4M N +N. 39 In the stretch dominated limit /N, /N the 6 6 Voigt matrix of effective elastic moduli is C 4 = C 6 +C 8 where C 6 = φ 6 I 3 E, C 8 = φ 8 J 9 E, 4 I all elements of the 3 3 matrix J are unity, and the volume fractions φ 6 = 6R A /V, φ 8 = 8R A /V satisfy φ = φ 6 +φ 8. The three moduli follow from 4 as K = φ 9 E, µ = φ 8 9 E, µ = φ 6 E. 4 6

11 Isotropy, µ = µ µ, is achieved if 3φ 6 = φ 8, i.e. A = A µ = φ 5 E 4 in which case the effective Poisson s ratio is ν = 4, in agreement with 39. This effective solid is the 3D isotropic optimal material introduced by []. 5. Pentamode lattices 5.. Z = d+ and the pentamode limit As discussed in 44., Z = d+ is the minimal coordination number necessary for a fully positive definite elasticity tensor. We now examine this case in particular in the limit of stretch dominant deformation. Given that a PM is an elastic solid with a single Kelvin modulus the elastic stiffness must be of the form C = λs S, λ >, S Sym. 43 Note that the parameter λ is somewhat arbitrary since it can be replaced by unity by subsuming it into the definition of S. Since rankp = Z d it follows that the single non-zero eigenvalue of P of 34 is unity, i.e. there exists a d+-vector b such that Hence, 34 yields the moduli explicitly in the form 43 with P = bb T where b T b =. 44 d+ λ =, S = V / R i M / i b i P i. 45 The eigenvalue property Pb = b implies that b satisfies d+ b iu i =, i.e. it is closely related with the fact that the d+ vectors u i are necessarily linearly dependent. Alternatively, b follows by assuming that C of eq. 34 has PM form C = S S, then use CI = StrS and I : CI = trs, from which we deduce that the moduli have the form 43 with λ = d+, d+ V γ k S = k= γ i P i where γ i = R i R d+ i M i M i k= P d+ k R j. 46 M k M j= j Equations 43, 45 and 46 provide two alternative and explicit formulas for the PM moduli. It is interesting to note that either of the above formulas for C leads to an expression for the axial force in member i based on equations 7 and 8. Thus, using eq. 46 gives f i = VλS : ǫri γ i e i. It may be checked from the definition of γ i that the forces are equilibrated, since d+ R i γ i e i =. 47 This identity implies that γ i = for some member i only if but not iff the remaining d members are linearly dependent. When this unusual circumstance occurs the member i bears no load since f i = for any applied strain. For instance, if two members are collinear in D, say members and, then the third member is not load bearing only if R γ = R γ. When d of the members span a d -plane the remaining member is non-load bearing if it is orthogonal to the plane. Writing S in terms of its principal directions and eigenvalues, S = s q q + s q q + s 3 q 3 q 3, where {q,q,q 3 } is an orthonormal triad, it follows that the elastic moduli in this basis are C IJ = λs I s J if I,J {,,3}, otherwise. 48

12 .5 C D.5 C 3D.5 C 3D C D θ deg Figure : The elastic moduli for D and 3D PM lattices with rods of equal length R = R and stiffness M = M as a function of the junction angle θ. Note that the D 3D moduli are identical at the isotropy angle The axial stiffness C vanishes at θ = π. Since C = C C it follows that C also vanishes at θ = π. The material symmetry displayed by PMs is therefore isotropic, transversely isotropic or orthotropic, the lowest symmetry, depending as the triplet of eigenvalues {s,s,s 3 } has one, two or three distinct members. The five easy pentamode strains correspond to the 5-dimensional space S : ǫ =. Three of the easy strains are pure shear: q i q j +q j q i, i j and the other two are s q q s q q and s q 3 q 3 s 3 q q. Any other zero-energy strain is a linear combination of these. 5.. Poisson s ratio of a PM In practice there must be some small but finite rigidity that makes C full rank, the material is unstable otherwise. The five soft modes of the PM are represented by < {µ i, i =,...,5} K where the set of generalized shear moduli must be determined as part of the full elasticity tensor. A measurable quantity that depends upon the soft moduli is the Poisson s ratio: for a given pair of directions defined by the orthonormal vectors n and m the Poisson s ratio ν nm is the ratio of the contraction in the m-direction to the extension in the n-direction for a uniaxial applied stress along n, i.e. ν nm = mm : Mnn / nn : Mnn where M = C is the fourth order tensor of elastic compliance. As an example, consider the diamond-like structure of Fig. a with shear moduli given by Table, Z = 4. In the pentamode limit K µ = 3µ, with n i, m i as the components in the principal axes, we obtain see e.g. [8] ν nm = n m n m n 3m 3 n 4 +n4 +n4 3 [, ]. 49 The actual values of the soft moduli {µ i, i =,...,5} are sensitive to features such as junction strength and might not be easily calculated in comparison with the pentamode stiffness. An estimate of the Poisson effect can be obtained by assuming the five soft moduli equal, in which case C C of eq. 43 is modified to Cµ C+µ I λtrs C, 5 which is invertible and positive definite for µ >. Using M = C µ define ν nm µ, then the limit exists as the shear modulus is reduced to zero: ν nm ν nm where ν nm = m Smn Sn S : S n Sn. 5 For the example of Figure a S = I and eq. 5 gives ν nm = /. Generally, the values of ν nm from eq. 5 associated with the principal axes of S see 48 are ν ij = s i s j /s j + s k, i j k i. If s > s > s 3 > then the largest and smallest values are ν > and ν 3 <, respectively. Compare this

13 with the Poisson s ratio of an incompressible isotropic elastic material: ν =. Negative values of Poisson s ratio occur if the principal values of S are simultaneously positive and negative. a θ = 5 b θ = 6 c θ = 7 Figure 3: Each of these two-dimensional PM lattices have isotropic quasi-static properties. The ratio of the R red to R blue is determined by 55. The pure honeycomb structure is θ = 6..5 ν ij.5 3D D θ deg a ν, ν b D: θ = 5, Figure 4: a The solid curves show Poisson s ratio ν for the same configuration as Figure R = R, M = M. ν describes the lateral contraction for loading along the axial e-direction. The related Poisson s ratio ν = ν / + ν is shown by the dashed curves. b The D lattice for θ = 5 top and θ = Transversely isotropic PM lattice Assume the unit cell has symmetry consistent with macroscopic transverse isotropy. It comprises two types of rods: i = with R, M in direction e = e, and i =,...,d + with R, M in directions e i symmetrically situated about e with e e i = cosθ. Let c = cosθ, s = sinθ. We find after some simplification, that 43 and 46 give the PM elastic stiffness as C = ds 4 R I+β e e I+β e e Vd dc M +M 5 3

14 C ii C C 3 ν.5 ν ij.5 ν p p Figure 5: The principal stiffnesses left and Poisson s ratios right for a diamond lattice with the center atom shifted along the cube diagonal. The four vertices of the unit cell at,,,, and the center junction atom lies at ppp. Isotropy is p =. where the non-dimensional parameter β and the unit cell volume V are β = d c { R +cr s, V = sr d 4, d =, R +cr R 6 3, d = Note that the elasticity of the rods enters only through the combination dc M +M. The nondimensional geometrical parameter β defines the anisotropy of the pentamode material, with isotropy iff β =. If β > the PM is stiffer along the axial or preferred direction e than in the orthogonal plane, and conversely it is stiffer in the plane if < β <. The axial stiffness vanishes if β = which is possible if θ = π. The unit cell becomes re-entrant if θ > π c <. If c < then β < and the principal values of S are simultaneously positive and negative with the negative value associated with the axialdirection. Note thatr +cr must bepositivesince theunit cellvolumev is positive. AsR +cr the members criss-cross and the infinite lattice becomes stacked in a slab of unit thickness, hence the volume per cell tends to zero V. Let e, the axis of transverse isotropy, be in the -direction. A transversely isotropic elastic solid d = 3 has 5 independent moduli: C, C = C 33, C = C 3, C 3 and C 66 = C 55 with C 44 = C C 3. The PM has C 66 = and C 3 = C C 44 = and C C = C, which are consistent with rankc =. The D version, technically of orthotropic symmetry, is defined by 4 independent moduli C, C, C and C 66, which in the PM limit satisfy C 66 = and C C = C. In either case the non-zero moduli are C C β d cs R R +cr = K C C β where K = d VM +dc M. 54 The PM is isotropic for β =, i.e. when the angle θ and R /R are related by R R = dcosθ d cosθ isotropy β =. 55 Hence, isotropy can be obtained if θ [cos d, π ] with the proper ratio of lengths, see Fig. 3. At the limiting angles R R as θ cos d θ π. If the lengths are equal R = R isotropy is obtained for cosθ = d, i.e. θ =6, 7.53, for d =, 3, corresponding to hexagonal and tetrahedral unit cells, respectively. Some examples of isotropic PMs and their properties are illustrated in Fig. 3. Transversely isotropic PMs are considered in Figs ThestiffnessparameterK of 54isthebulkmodulusoftheisotropicPM.NotethatK isnotequivalent to K of 3 since the latter is consequent upon the condition 3 which is not assumed here. Instead, eqs. 4

15 5, 54 imply that the isotropic PM bulk modulus for uniform members is K = Kf, f = d s 4[ d + A A dc dc ] [ d + A A dc dc ] 56 where A, A are the cross-sectional areas strut thicknesses for d =. For a given θ and d, f with equality iff A A = dc. Hence the maximum possible isotropic effective bulk modulus for a given volume fraction φ is precisely K of 3. This result agrees with [4, eq..] for d =, and with the bulk modulus for a regular lattice with tetrakaidecahedral unit cells [5], [6], i.e. an open Kelvin foam, see Fig. 6. The latter structure, comprising joints with 4 struts and a unit cell of 4 faces 6 squares and 8 hexagons, has cubic symmetry; however the two shear moduli are almost equal so that the structure is almost isotropic. In fact, if the struts are circular and have Poisson s ratio equal to zero then the effective material is precisely isotropic with shear modulus µ = 4 9π φ E [5]. Figure 6: The tetrakaidecahedral open foam unit cell [5] has low density PM behavior similar to the diamond lattice. Note that [9] considered a tetrahedral unit cell of four identical half-struts that join at equal angles and found K = φe 8 not φe 9 ; the difference arises from taking the cell volume for the tetrahedron, but since the tetrahedron is not a space filling polyhedron, this is not the correct unit volume to use. 6. Two dimensions: a special case 6.. Shear force as a nodal bending force For d = the total force 4 on member i can be simplified as f i = M i r i e i + N ij Rj ψ ij e ij 57 with generalized nodal compliance N ij given by N ij = N ij + N b ij where N b ij N in j Rk. 58 R i R j N k k Hence, the shear force can be considered as an equivalent nodal bending force. Significantly, the moments of the shear forces are now automatically equilibrated due to the symmetry N ij = N ji. Equation 57 follows by first noting that the vector moment of the shear force is in the direction perpendicular to the plane of the lattice, say a 3. Define the angle of deflection associated with flexural bending: θ i a 3 e i r b i /R i. The moment of the shear force is R i f b i = R i /N iθ i a 3, and the moment equilibrium condition 7 becomes i Ri θ i = θ i = N i k 5 Rk N k R j N j θ i θ j. 59

16 However, θ i θ j = ± ψ ij more precisely θ i θ j = ψ ij a 3 e j e i / e j e i, therefore eq. 59 allows us to express the single shear force acting on member i as the sum of nodal bending forces with compliances N b ij, from which eq. 57 follows. The significance of eq. 57 is that it allows us to express the effective moduli for d = as follows: Define d i = e i Mi, D i = R i e i e i VMi, d ij = D ij = R i R j VN ij ei e ij +e j e ji where N ij R i R j eij N ij + e ji, 6a R i R j = + R ir j Z R k, N ij N i N j N k k= 6b {u k } L k= = {d i, d ij }, {U k } L k= = {D i, D ij }, L = ZZ +/, 6c N L P ij = δ ij u i u k u k uj C = P ij U i U j. 6d k= Note that this result is valid for any similarly situated D lattice structure; in particular it does not require the zero rotation assumption. 6.. Example: honeycomb lattice As an application of eq. 6 we consider the transversely isotropic lattice of 55.3 in -dimensions Z = 3, now including the effects of the bending compliances of the individual members, N and N. Using the same notation as in 55.3 we find C C = cs β N +s c M, c C M +M N +s M M β N +s M +c M, N M, 6 C 66 = sr R +cr s R N +R N M. + cr +R These are in agreement with the in-plane moduli found by [3]. Note that the moduli of eq. 6 reduce to the PM moduli 54 as the bending compliance N, independent of the bending compliance N. 7. Conclusions Our main result, eq. 9, is that the effective moduli of the lattice structure can be expressed C = L i,j= P iju i U j where L = Zd + ZZ /, U i are second order tensors, and P ij are elements of a L L projection matrix P of rank L d. Explicit forms for the parameters {U i, P ij } have been derived in terms of the cell volume, and the length, orientation, axial and bending stiffness of each of the Z rods. This Kelvin-like representation for the elasticity tensor implies as a necessary although not sufficient condition for positive definiteness of C that the rank of P exceed 3d 3, which is satisfied if the coordination number satisfies Z d +. The L second order tensors {U i } are split into Z stretch dominated and ZZ / bending dominated elements. The latter contribute little to the stiffness in the limit of very thin members, in which case the elastic stiffness is stretch dominated and, at most, of rank Z d. The formulation developed here is applicable to the entire range of stiffness possible in similarly situated lattice frameworks, from the Z = 4 structure proposed by [] with full rank C to pentamode materials corresponding to coordination number Z = d+, with C of rank one. i,j= Acknowledgment Thanks to Adam Nagy and Xiaoshi Su for discussions and graphical assistance. This work was supported under ONR MURI Grant No. N4363 6

17 References [] Gurtner, G & Durand, M, 4 Stiffest elastic networks. Proc. R. Soc. A 47, 3,6+. doi:.98/rspa.3.6. [] Bell, AG, 93 Tetrahedral principle in kite structures. Nat. Geogr. Mag. 6, 9 5. [3] Zheng, X, Lee, H, Weisgraber, TH, Shusteff, M, Deotte, JR, Duoss, E, Kuntz, JD, Biener, MM, Kucheyev, SO, Ge, Q, Jackson, J, Fang, NX, & Spadaccini, CM, 4 Ultra-light, ultra-stiff mechanical metamaterials. Science 344, doi:.6/science.59. [4] Gurtner, G & Durand, M, 9 Structural properties of stiff elastic networks. Europhys. Lett., 4,+doi:.9/ /87/4. [5] Milton, GW & Cherkaev, AV, 995 Which elasticity tensors are realizable? J. Eng. Mater. Tech. 7, doi:.5/ [6] Milton, GW, The Theory of Composites. Cambridge University Press, st edition. [7] Fleck, NA, Deshpande, VS, & Ashby, MF, Micro-architectured materials: past, present and future. Proc. R. Soc. A 466, doi:.98/rspa..5. [8] Maxwell, JC, 864 On the calculation of the equilibrium and stiffness of frames. Phil. Mag. 7, [9] Sun, K, Souslov, A, Mao, X, & Lubensky, TC, Surface phonons, elastic response, and conformal invariance in twisted kagome lattices. Proc. Nat. Acad. Sci. 9,,369,374. doi:.73/pnas [] Deshpande, VS, Fleck, NA, & Ashby, MF, Effective properties of the octet-truss lattice material. J. Mech. Phys. Solids 49, doi:.6/s [] Deshpande, VS, Ashby, MF, & Fleck, NA, Foam topology: bending versus stretching dominated architectures. Acta Materialia 49, doi:.6/s [] Hutchinson, RG & Fleck, NA, 6 The structural performance of the periodic truss. J. Mech. Phys. Solids 54, doi:.6/j.jmps [3] Thomson, W, 856 Elements of a mathematical theory of elasticity. Phil. Trans. R. Soc. Lond. 46, [4] Norris, AN, 8 Acoustic cloaking theory. Proc. R. Soc. A 464, doi:.98/rspa [5] Kadic, M, Bückmann, T, Schittny, R, & Wegener, M, 3 On anisotropic versions of three-dimensional pentamode metamaterials. New J. Phys. 5, 3,9+. doi:.88/367-63/5//39. [6] Christensen, RM, Mechanics of cellular and other low-density materials. Int. J. Solids Struct. 37, doi:.6/s [7] Nye, JF, 985 Physical Properties of Crystals. Clarendon Press, Oxford. [8] Tollenaere, H & Caillerie, D, 998 Continuous modeling of lattice structures by homogenization. Adv. Engng. Software 9, doi:.6/s [9] Martinsson, PG & Babuška, I, 7 Homogenization of materials with periodic truss or frame micro-structures. Math. Models Methods Appl. Sci. 7, doi:.4/s857x. [] Martinsson, PG & Babuška, I, 7 Mechanics of materials with periodic truss or frame micro-structures. Arch. Rat. Mech. Anal. 85, 34. doi:.7/s [] Wang, XL & Stronge, WJ, 999 Micropolar theory for two-dimensional stresses in elastic honeycomb. Proc. R. Soc. A 455, 9 6. doi:.98/rspa [] Kumar, RS & McDowell, DL, 4 Generalized continuum modeling of -D periodic cellular solids. Int. J. Solids Struct. 4, doi:.6/j.ijsolstr [3] Wang, Y & Cuitino, A, Three-dimensional nonlinear open-cell foams with large deformations. J. Mech. Phys. Solids 48, doi:.6/s [4] Christensen, RM, 995 The hierarchy of microstructures for low density materials. In J Casey & M Crochet, editors, Theoretical, Experimental, and Numerical Contributions to the Mechanics of Fluids and Solids. Birkhäuser Basel, pages doi:.7/ \ 7. [5] Warren, WE & Kraynik, AM, 997 Linear elastic behavior of a low-density Kelvin foam with open cells. J. Appl. Mech. ASME 64, doi:.5/ [6] Zhu, HX, Knott, JF, & Mills, NJ, 997 Analysis of the elastic properties of open-cell foams with tetrakaidecahedral cells. J. Mech. Phys. Solids 45, doi:.6/s [7] Walpole, LJ, 984 Fourth rank tensors of the thirty-two crystal classes: multiplication tables. Proc. R. Soc. A 39, doi:.98/rspa [8] Norris, AN, 6 Poisson s ratio in cubic materials. Proc. R. Soc. A 46, doi:.98/rspa [9] Warren, WE & Kraynik, AM, 988 The linear elastic properties of open-cell foams. J. Appl. Mech. ASME 55, doi:.5/ [3] Kim, HS & Al-Hassani, S, 3 Effective elastic constants of two-dimensional cellular materials with deep and thick cell walls. Int. J. Mech. Sc. 45, doi:.6/j.ijmecsci