## Abstract

Soft phononic crystals have the advantages over their stiff counterparts of being flexible and reconfigurable. Normally, the band gaps of soft phononic crystals will be modified after deformation due to both geometric and constitutive nonlinearity. Indeed these are important properties that can be exploited to tune the dynamic properties of the material. However, in some instances, it may be that one wishes to deform the medium while retaining the band gap structure. A special class of soft phononic crystals is described here with band gaps that are independent or almost-independent of the imposed mechanical deformation, which enables the design of phononic crystals with robust performance. This remarkable behaviour originates from transformation elasticity theory, which leaves the wave equation and the eigenfrequencies invariant after deformation. The necessary condition to achieve such a property is that the Lagrangian elasticity tensor of the hyperelastic material should be constant, i.e. independent of deformation. It is demonstrated that incompressible neo-Hookean materials exhibit such a unique property. Semilinear materials also possess this property under special loading conditions. Phononic crystals composed of these two materials are studied theoretically and the predictions of invariance, or the manner in which the response deviates from invariance, are confirmed via numerical simulation.

## 1. Introduction

Phononic crystals (PCs) are periodic structures that can control the propagation of acoustic or elastic waves via wave filtering in specific frequency ranges [1–10]. Potential applications of PCs include waveguides and filters, sensors, redirectivity devices and many more. Early studies focused predominantly on *stiff* PCs, which have negligible deformation under general external loads. Recently, however, attention has switched to the study of *compliant* or *soft* phononic crystals (SPCs) due to their potential for flexibility, tunability and multifunctionality. To date, almost all work on SPCs is related to their tunability, i.e. band gaps are tuned via mechanical deformation or some other mechanism such as an imposed electrical or magnetic field.

A variety of tuning methods using mechanical deformation have been proposed. Simple one-dimensional models were studied initially [11,12] and since then Bertoldi and co-workers have conducted extensive studies on tuning the band gaps of hyperelastic PCs with the aid of structural instability [13–15]. Rudykh & Boyce [16] discovered that the wrinkling of thin stiff layers embedded in a compliant matrix could tune band gaps and more general wave propagation behaviour. Mousanezhad *et al.* [17] designed flexible honeycomb structures with tunable band gaps under compression and buckling. Most of these aforementioned works are related to body waves. Recently, however, Li *et al.* [18] demonstrated that surface wave band gaps can also be tuned by controlling the surface wrinkling patterns of soft materials. In addition to exploiting instability, a number of other methods have been proposed to tune the dynamic behaviour of PCs. Tang *et al.* [19] designed super-stretchable structures with cut hinges to achieve both tunable band gaps and enhanced strength simultaneously. Barnwell *et al.* [20,21] proposed the idea of tuning band gaps of soft crystals by imposing local pre-deformation without inducing global deformation. Galich *et al.* studied band gap tunability in layered SPCs [22]. Magnetic and electrical effects to tune band gaps have also been described [23–27] having the advantage that contact with the structure is not required, although the constituent materials themselves are generally more complex. Significant progress has therefore been made in the design of tunable PCs by employing mechanical deformation and related techniques.

Given that the majority of the work thus far published has centred on the notion that SPCs are employed to *tune* the band gaps of the medium, an interesting question arises as to *whether it is possible to design SPCs with invariant band gaps even after large deformation?* If this is answered in the affirmative, one would be able to design PCs with robust frequency response while maintaining flexibility and multifunctionality. The primary aim of this work therefore is to explore the conditions required to achieve this apparently abnormal behaviour. It will be shown that this remarkable phenomenon is closely related to the transformation elasticity theory used to design elastodynamic cloaks [28–30].

In §2, the notion of band gap invariance in one-dimensional structures is first discussed before proceeding to the general theory for two-dimensional structures in §3. The properties of the required materials are discussed in §4 and simulations illustrating the influence of deformation on wave propagation for such special materials are then described in §§5 and 6. Conclusion and a general summary are given in §7.

## 2. One-dimensional phononic crystals under pre-stretch

In order to illustrate the effect of pre-deformation on the band structure and wave propagation characteristics of an SPC, consider a simple one-dimensional medium as depicted in figure 1*a*.

The unit cell of this PC is composed of a compliant layer and a stiff layer, having initial lengths *a*_{1} and *a*_{2}, where the subscript indicates the respective phase of the medium. The initial lattice parameter of the PC is therefore *a*=*a*_{1}+*a*_{2}. Under uniaxial stretch perpendicular to the faces of the layers, if the elongation of the stiff layer is negligible, the unit cell of the PC deforms to a length *Λa*, where *Λ* is the overall stretch ratio. Either a Lagrangian or Eulerian scheme can be used to describe the deformation and equilibrium of the PC. Assume that the per-unit-length stiffness and mass density are denoted by *C*_{E} and *M*_{E}, respectively, where the subscript ‘E’ or ‘L’ indicates the Eulerian or Lagrangian scheme, respectively. Therefore, whenever the Lagrangian scheme is used, the corresponding stiffness and mass density are obtained as *C*_{L}=*C*_{E}/λ and *M*_{L}=*λM*_{E}, respectively, where λ is the stretch ratio of any material point. In this case, the wave equation for a linear wave superimposed on the pre-deformation is
*u* is the displacement, *X* and *x* indicate the reference and spatial coordinates (figure 1*a*), respectively, *t* is time and *f*_{,x} indicates the derivative of the function *f* with respect to *x*.

The dispersion relation, i.e. the eigenfrequency *ω*(*K*) as a function of the Lagrangian wavevector *K*, can be readily obtained by assuming a Bloch wave solution
*Λ* of the structure. Hence, the band structure of an SPC will usually be tuned under pre-stretch. An unusual but intriguing phenomenon would arise when
*ω*(*K*) *is independent of deformation*. Therefore, the band structure of the PC will be independent of the uniaxial stretch if equation (2.4) is satisfied.

The longitudinal wave speed is *v*_{L} is independent of the stretch if *C*_{L} is a constant. In this case, the Eulerian wave speed *v*_{E}=*λv*_{L} is proportional to the stretch ratio λ, i.e.
*c*, the wavefront will arrive at the same material point during a given time interval, no matter how large the pre-stretch ratio is. By contrast, the spatial wave speed is proportional to the stretch ratio according to equation (2.5).

An example is provided in order to illustrate the band structure and eigenmodes of an SPC with deformation-independent band gaps. The geometric and material constants are taken as *a*_{1}=2*a*/3, *a*_{2}=*a*/3, *v*_{L2}=10*v*_{L1}, respectively. As illustrated in figure 1*b*, the band structure, if plotted in the Lagrangian scheme as *ω*-*K*, is independent of the overall stretch ratio *Λ*. However, the wave mode *u*(*x*) will be stretched in the spatial domain, as shown in figure 1*c*, which depicts the wave mode for one unit cell (*K*=0.315/*a*, *ω*=9.425*v*_{L1}/*a*). Therefore, the amplitude, phase and shape of the waves are actually retained during the deformation due to the constant Lagrangian wave speed *v*_{L} in the constituent materials. In fact, the material points do not sense any changes of their own motion or neighbourhood in this case. This is the essential reason causing such an intriguing phenomenon with deformation-independent band gaps in one-dimensional SPCs. One could go on to consider the type of material required to produce this type of effect in this one-dimensional scenario. However, given that the purpose of this section was primarily to illustrate the notion of invariance this is bypassed and more complex, realistic configurations are now considered.

## 3. Soft phononic crystals in higher dimensions

In §2, it has been shown that the band gap structure of an SPC in one dimension is independent of stretch and deformation if the Lagrangian wave speed is constant. Does this phenomenon exist in higher dimensions? In order to address this question, the fundamental theory for SPCs will be introduced in this section and a necessary condition for deformation-independent band structures will be discussed.

### (a) Soft phononic crystals with pre-deformation

As shown in figure 2, a two-dimensional PC is considered in this work, but the proposed theory and mechanism also applies for SPCs in general. Although more general designs are possible as will be described below, consider for now the case when the SPC comprises stiff cylinders embedded on a square lattice in a compliant matrix. Therefore, it is possible to deform the PC in figure 2*a* in a finite manner, into the current configuration as depicted in figure 2*b* via the global affine lattice deformation. Correspondingly, the lattice vectors **A**_{i} (*i*=1,2,3) in the reference material space will transform to lattice vectors **a**_{i} in the current material space via
**F**_{L} designates the affine deformation gradient [31–33] of the lattice points. Note that the reference position **X** and current position **x** of an arbitrary material point will not in general obey the affine transformation, i.e. **x**≠**F**_{L}**X** as illustrated in figure 2*a*,*b*, due to the periodic non-affine deformation **F**_{p}(**X**) within a unit cell and hence the total deformation gradient [33] is **F**(**X**)=**F**_{L}**F**_{p}(**X**).

Corresponding to the material spaces, the reference reciprocal space (figure 2*c*) and current reciprocal space (figure 2*d*) can also be related by a mapping. More specifically, the mapping between these two spaces is affine since the reciprocal space is derived from the lattice points in the reference material space, regardless of the constituent materials and interior geometry of the material unit cell. The lattice vectors in the reference and current reciprocal spaces are denoted as (**B**_{1},**B**_{2},**B**_{3}) and (**b**_{1},**b**_{2},**b**_{3}), respectively, which are defined by **A**_{i}⋅**B**_{j}=**a**_{i}⋅**b**_{j}≡2*πδ*_{ij}, where *δ*_{ij} is the Kronecker delta tensor. Therefore, the affine mapping between the two reciprocal spaces is derived as
**K** and **k**, respectively, they can be related via the simple relationship

### (b) Incremental dynamic equation

Linear elastic wave propagation in SPCs with large pre-deformation is described by incremental elastodynamic equations [34], derived via what is often called the theory of *small-on-large*. This theory has been widely used to study elastic wave propagation and scattering in soft materials with pre-deformation [13,14,16,18,20,22,28,29,35,36]. There are two equivalent sets of incremental equations that can be employed to study wave propagation in the pre-stressed state: the *Lagrangian* and *Eulerian* forms [34]. The former specifies equations in the reference (undeformed) material space while the latter specifies the current (deformed) material space. The Lagrangian form is often convenient since boundary conditions can be imposed in a more straightforward manner. The Lagrangian form is adopted in this paper; this gives
*f*_{,i}=∂*f*/∂*X*_{i}, *u*_{i} is the *i*th component of displacement **u**, resolved along Cartesian coordinates, *ρ*_{0} is the material density in its undeformed state and **F** and its components are defined by
*W*(**F**) is the strain energy function of the hyperelastic material under deformation [34,37].

The Bloch wave solution is taken as *ω*(**K**) of the PCs are actually obtained by solving the eigenvalue problem for equation (3.5) by prescribing the Bloch condition in equation (3.7) for a given wavevector **K**.

### (c) Transformation elasticity and eigenfrequency invariance

Transformation elasticity theory [29] provides a scenario where the incremental wave equation (3.5) is invariant to pre-deformation. This is achieved when the ET is independent of the deformation, i.e.
**F**. Similar to the one-dimensional example, an intriguing phenomenon arises when equation (3.8) is satisfied. In this case, the incremental equations (3.5) further simplify to

## 4. Required hyperelastic strain energy functions

### (a) Materials with constant elasticity tensor

What kind of hyperelastic material will exhibit an ET **S**=∂*W*/∂**E** is the second Piola–Kirchhoff stress, **E**=(**F**^{T}**F**−**I**)/2 is the Green–Lagrange strain tensor. The first term on the right-hand side of (4.2) depends on the material stretch only, while the second term depends on both material stretch *and* rotation. Therefore, in order to obtain a constant ET **S** should be constant and the instantaneous modulus *W* must be a linear function of the first invariant of **E**, i.e.
*μ*_{0} is the shear modulus and **E**=*E*_{ii}. The strain energy function in equation (4.2) will inevitably result in a hydrostatic stress term in the undeformed state and it does not appear possible to modify this form for compressible materials in a manner consistent with the imposition that *incompressible* materials however a Lagrange multiplier term *p* can be added, i.e. the strain energy function becomes

### (b) Materials with almost-constant elasticity tensors

As has been illustrated above, it is very challenging to find materials with deformation-independent elastic tensors when subjected to general deformations. Therefore, attention is now restricted to materials with almost-constant ETs or where specific modes of deformation permit invariance of specific wave types. The explicit form of **U** or left stretch tensor **V**. Denoting such a representation as *W*_{i}=∂*W*/∂λ_{i},*W*_{ij}=∂^{2}*W*/∂λ_{i}∂λ_{j} and, importantly, the summation rule does not apply to repeated indices. Note that the ET *is* a possibility.

The formulae in equations (4.5)–(4.8) imply that for invariance, the strain energy function *W* should be a quadratic polynomial function of the principal stretches λ_{1},λ_{2},λ_{3} at least. Any higher powers will in general lead to stretch-dependence. In addition, the constitutive theory requires that *W* be expressed in terms of the invariants of the right stretch tensor **U** for isotropic hyperelastic materials. Hence, the following polynomial form of strain energy function is chosen
*c*_{i} (*i*=0,1,2,3) are unknown coefficients to be determined, and *i*_{1}=*tr* **U**, *i*_{2}=(*tr*^{2}**U**−*tr* **U**^{2})/2 are the first two invariants of **U**. Two of the coefficients *c*_{i} can be eliminated by using the energy-free (*W*=0) and stress-free (∂*W*/∂λ_{i}=0) conditions in the undeformed state (λ_{i}=1). It follows that the strain energy function required is the semilinear form, i.e. [29]
_{0} and *μ*_{0} are the Lamé constant and shear modulus, respectively. The first Piola–Kirchhoff stress **P** of the semilinear model is derived from the relation **R**^{T}**P**=∂*W*/∂**U** [38], as
**R** is the rotation tensor of **F**.

The ET of the semilinear model is derived by substituting equation (4.10) in equations (4.5)–(4.8), as

For special cases, the components can be further simplified. A special case is the situation when *W*_{i}+*W*_{j}=0 which renders the right-hand side of (4.13) zero. As was determined in equation (4.14) of [29], for in-plane incremental displacements superposed on a finite homogeneous deformation, the following condition achieves this:
*ν*_{0} is Poisson's ratio. In this case then, the ET related to the *in-plane* deformation degenerates to a constant tensor, i.e.
_{1},λ_{2}) nor the in-plane rotation. This is therefore a special case when the semilinear material exhibits an invariant ET and this will be exploited to design an invariant SPC of type (ii) as defined in §3c.

A straightforward case to consider, but which will inevitably lead to stretch-dependence of _{3}=1). In this situation, equation (4.13) reduces to the following form:
_{0}+*μ*_{0} in equation (4.17) should be minimized to reduce the effect of stretch on *ν*_{0} because λ_{0}+*μ*_{0}=*μ*_{0}/(1−2*ν*_{0}). If Poisson's ratio

In addition to the effect of stretch, as described above, the transformed ET

## 5. Wave speed invariance

In §2, it was shown that band gap invariance in one-dimensional SPCs is closely related to Lagrangian wave speed invariance. Therefore, the effect of pre-deformation on the wave speed of plane waves in *homogeneous* neo-Hookean and semilinear materials is studied in this section. A homogeneous sample is stretched uniformly along the *X*_{1} direction with stretch ratio λ_{1}=*Λ*. The medium is constrained along the *X*_{3} direction, i.e. λ_{3}=1 is imposed. The condition associated with stretch in the *X*_{2} direction will be chosen for each example to be considered. For neo-Hookean materials, the focus shall be on the horizontally polarized (SH mode, polarized in the *X*_{3} direction) and vertically polarized (SV mode, polarized in the *X*_{2} direction) shear waves, while for the semilinear material the additional longitudinal (P mode) wave is considered.

First, consider shear waves (SV and SH modes) propagating in the *X*_{1} direction in an incompressible neo-Hookean material under uniform pre-stretch. The stretch ratios are taken to be λ_{1}=*Λ*,λ_{2}=*Λ*^{−1},λ_{3}=1 so the deformation is isochoric. The Lagrangian wave speed for plane waves can be derived from equation (3.5) as
*i*=1,2,3 for the P, SV and SH wave mode, respectively. The ET components can be derived from equation (4.4), as

The analysis of the semilinear material is more complicated than the neo-Hookean material. Suppose that the medium is stress-free in the *X*_{2} direction, i.e. imposing *P*_{22}=0, so that the pre-stretch ratios can be shown to be (using (4.11))
*v*_{L1} of the P wave is therefore independent of the stretch ratio *Λ*, while the shear wave speeds *v*_{L2} or *v*_{L3} of the shear wave modes depend on the stretch ratio *Λ*. As shown in figure 3, the Lagrangian wave speeds of the SV and SH waves in the semilinear material are slightly affected by the stretching ratio *Λ*, and the stretch-dependent effect is reduced when Poisson's ratio *ν*_{0} decreases. In addition, the SH wave speed is found to be independent of stretch when Poisson's ratio *ν*_{0}=−1, which can be observed in figure 3 and verified from equation (5.5).

Now consider a third case; the special scenario associated with the imposition of equation (4.15), which gives the invariant in-plane ET as in equation (4.16). The stretch ratios are taken as λ_{1}=*Λ*, λ_{2}=2−*Λ* and λ_{3}=1. The Lagrangian wave speeds of the in-plane P-SV waves are not affected by the pre-deformation. However, the components of the ET corresponding to SH wave propagation will be affected by the pre-stress. Suppose that an SH wave propagates in a direction defined by an angle *θ* subtended from the *X*_{1} axis. The corresponding Lagrangian wave speed is then
*Λ*. This effect of induced anisotropy for the SH wave is illustrated in figure 4.

To summarize, the analysis above indicates that the Lagrangian wave speed of shear waves in an incompressible neo-Hookean material is independent of the deformation. Hence, the neo-Hookean material has potential in designing PCs with deformation-independent band gaps for shear wave modes. On the other hand, the semilinear material cannot exhibit such total invariant behaviour. Having said this, the Lagrangian wave speed of plane waves in a deformed semilinear material are unaffected (for P waves) or slightly affected (for SV and SH waves) by pre-deformation. Hence, one should expect that band gaps in semilinear PCs can be almost-independent of the pre-deformation in general. For a specific case, the semilinear material may exhibit a constant in-plane ET for the special pre-deformation mode as in equation (4.15). This gives rise to the possibility of designing SPCs with deformation-independent in-plane (P/SV) band gaps. These hypotheses will be verified in the next section where wave propagation in various configurations of SPCs will be simulated numerically.

## 6. Simulation and results

### (a) Simulation model

Based on the above analysis, a PC comprising an incompressible neo-Hookean material will have deformation-independent band gaps for shear waves, while one composed of a semilinear material will have band gaps that are almost-invariant for all wave modes. Full in-plane invariance of the semilinear material is achievable for the special deformation case given in equation (4.15), which leads to invariant P-SV band gaps. The effect is illustrated below by considering two examples of SPCs of type (iii) and one example of type (ii) as classified at the end of §3c. An infinitely thick PC (figure 2*a*) is considered with cylinders organized on a square lattice. For the medium of type (iii), the inclusion phase is stiff enough so that its ET is independent of deformation since the induced strain is infinitesimal. For the type (ii) material the inclusion is a void. The matrix phase is neo-Hookean in the first example, and semilinear in subsequent examples. The phononic band structures are calculated by using the finite-element package ABAQUS (v. 6.13-1) and adopting the method proposed by Åberg & Gudmundson [39,13]. One layer of linear brick elements C3D8R are built for the simulation model and the Bloch condition is imposed on the boundaries of the unit cell. A unique feature of this model is that the in-plane mode (P-SV waves) and antiplane mode (SH waves) can be decoupled. The band structures are therefore plotted separately for these two scenarios.

### (b) Phononic crystals of neo-Hookean type

For the first example, a compressible neo-Hookean material is employed so that both the P-SV and SH waves can propagate, although the latter is of specific interest. The corresponding strain energy function is taken as
*I*_{1}= **U**^{2}, *J*= det **F**, the shear modulus *μ*_{m}=25.9 MPa, bulk modulus *κ*_{m}=10*μ*_{m} and density *ρ*_{m}=1000 kg m^{−3}. The cylindrical inclusions are taken to be aluminium with shear modulus *μ*_{c}=25.9 GPa, Poisson's ratio *ν*_{c}=0.33, and density *ρ*_{c}=2700 kg m^{−3}. Additionally in this example, uniaxial stretching (*ϵ*=15% elongation with the lateral sides left free, where *ϵ* is the engineering strain calculated over the unit cell) is applied to the PC, as illustrated in figure 5*b*. For the stretch deformation, the lattice deformation gradient components are *F*_{L11}=1+*ϵ*, *F*_{L12}=*F*_{L21}=0 and *F*_{L22} is determined by the stress-free condition on the lateral sides. The band structures for the undeformed and stretched PCs are shown in figure 5*a*,*b*, respectively, with dispersion curves *ω*_{n}(**K**) calculated along the path *Γ*-*G*_{1}-M-*Γ*-*G*_{2}-M in the Brillouin zone (reference reciprocal space). The dimensionless frequency is defined as
*ω* is the angular frequency and *a*=0.1 *m* is the characteristic length of the unit cell. Additionally, the radius of the cylindrical inclusions is 0.03 m. By comparing the band gaps in figure 5, it is observed that the single-band gap (located in the vicinity of *x*_{1}−*x*_{2} mode) disappears after stretching. This is because in this case the ET components associated with in-plane waves are deformation-dependent. The stretch deformation induces symmetry breaking of the Brillouin zone so that the *G*_{1} and *G*_{2} points are no longer equivalent. This is seen by considering the third band *ω*_{3} associated with in-plane wave modes in figure 5*b*. By contrast, the position of the two band-gaps associated with antiplane waves (the *x*_{3} mode) remain unchanged after stretching since the neo-Hookean material with strain energy function defined in (6.1) behaves very similar to the invariant incompressible medium in response to antiplane shear waves. This confirms the result associated with neo-Hookean media derived theoretically above. Interestingly, invariance of antiplane band gaps has previously been observed in a neo-Hookean SPC [13], but the origin of this phenomenon has not been explained until now. The neo-Hookean material retains the total invariance of transformation elasticity [29]. Given that the antiplane band gaps of the neo-Hookean-type PCs are completely invariant to arbitrary pre-deformation, no further deformation examples for materials of this type will be discussed. Instead attention will focus on the possibility of invariance for in-plane waves.

### (c) Phononic crystals of semilinear type

Consider an SPC with semilinear hyperelastic matrix medium with Lamé constant λ_{m}=50.4 MPa, shear modulus *μ*_{m}=25.9 MPa and density *ρ*_{m}=1000 kg m^{−3}. The material configuration is the same as for the neo-Hookean material example considered in the previous section (aluminium cylinders), except that now the compliant matrix is semilinear. Because the semilinear medium is not completely invariant to deformation, two different deformation modes, stretch and shear, are considered for this example. The resulting band structures of the semilinear SPC are shown in figure 6 in the undeformed and deformed cases. For the shear deformation, the lattice deformation gradient components are *F*_{L12}=*γ*, *F*_{L21}=0, and *F*_{L11} and *F*_{L22} are determined from the free normal stress conditions on the four sides.

The ET *a*, only one band gap is observed for in-plane waves while two band gaps are determined for antiplane waves in the undeformed state. In the deformed state, the band gaps are slightly tuned, as shown in figure 6*b* for the stretching case (*ϵ*=0.15) and figure 6*c* for the shearing case (*γ*=0.3), respectively. The dependence of the pre-deformation on the band gaps is more clearly shown in figure 6*d*,*e* for the stretch and shear deformation cases, respectively. It is observed that the stretch deformation will affect the band gaps slightly, while the influence of the shear deformation on the band gaps is negligible. One explanation for this difference is that the maximum material stretch in the case of shearing is less than that in the case of uniaxial stretching. It is thus concluded that the band gaps in semilinear SPCs are slightly dependent on the pre-deformation and the influence also depends on the deformation mode considered. The advantage of using materials of semilinear type in SPCs then is that they exhibit relatively robust band gaps for both P-SV and SH waves. This more general case has broader significance and potential applications than using a neo-Hookean material, for which SPCs exhibit invariance only for the antiplane wave mode.

Another interesting phenomenon observed in figure 6 is the wave mode degeneracy at the high symmetry point M right above the second band gap of the *x*_{3} modes. It is found that the two dispersion curve branches, which originally coincide at point M in figure 6*a*, become separated in the corresponding band structures in figure 6*b*,*c*. This is actually induced by the symmetry breaking of the PC once the square lattice is distorted due to the deformation [33,40], e.g. the fourfold rotation and reflection symmetries may not be preserved in the deformed configuration. In theory, this kind of degeneracy change should not occur for an SPC comprising only neo-Hookean phases due to the invariance of the wave equations for the *x*_{3} modes. However, these two branches are slightly separated at point M in figure 5*b* due to the existence of the stiff elastic cylinders. Detailed analysis of the symmetry breaking is quite complicated and will be discussed in future work.

The final SPC to be considered is a material of type (ii), i.e. solid–void, with semilinear matrix. Fully invariant P-SV band gaps are possible when the medium is stretched in the *X*_{3} direction (the direction of the axes of the cylindrical voids) while the other two directions are left to be load-free. In this case, the in-plane deformation is homogeneous everywhere (λ_{1}=λ_{2}=1+*ν*_{m}−*ν*_{m}λ_{3}) and it can be verified that equation (4.15) is satisfied automatically. Therefore, the in-plane ET components are always constant (invariant) as shown in equation (4.16), while the ET components related to SH waves will depend on the axial stretch λ_{3}. Hence, while one should expect that the P-SV band gaps are invariant, the SH band gaps will be tuned under pre-deformation. The simulation results are shown in figure 7, which includes the band structures of the undeformed and pre-stretched SPCs. It is observed that the band structures of the P-SV wave modes are indeed independent of the axial stretch, despite the geometrical deformation, while the band gaps of antiplane waves are affected by the pre-deformation.

In summary, a material of type (iii) (neo-Hookean matrix and aluminium cylindrical inclusions) has been shown to exhibit invariant band gaps for SH waves (for arbitrary deformations), whereas a medium of type (ii) (semilinear matrix with cylindrical voids) exhibits invariant band gaps for P-SV waves (under axial stretch only).

## 7. Conclusion

It has been demonstrated here that it is theoretically possible to design SPCs with deformation-independent band gaps or band gaps that are close to being deformation-independent. This was achieved by appealing to hyperelastic transformation theory. This apparently unusual phenomenon has particular utility in engineering applications requiring both robust frequency response and structural flexibility. In order to achieve this non-trivial property, the key issue is to find materials with deformation-independent ETs so that the Lagrangian-formulated incremental equations possess the property of being invariant to deformation, which is the key notion of transformation elasticity theory. It has been proved that only incompressible neo-Hookean materials satisfy such a constraint exactly (i.e. under all deformations). Semilinear hyperelastic materials exhibit almost-constant ETs for the deformation modes considered here and in-plane invariance for specific materials under axial stretch. It was indicated that the invariance of the band gaps is associated with the invariance of the Lagrangian wave speeds in the medium, which are derived from transformation elasticity. Numerical simulations illustrate that a specified SPC with semilinear hyperelastic matrix and aluminium stiff cylindrical inclusions does indeed have band gaps that are almost-deformation-independent for wave modes in the frequency ranges considered. An example that considered the same cylindrical inclusions in an almost-incompressible neo-Hookean medium on the other hand illustrated that only SH wave band gaps are deformation independent, in accordance with the theoretical predictions. An example involving a solid–void PC of the semilinear type illustrated that deformation-independent P-SV band gaps are possible, but only when the structure is stretched axially. This can be thought of as the counterpart to the invariant SH band gap case for neo-Hookean materials, although while the neo-Hookean case holds for arbitrary deformation the semilinear SPC is only fully invariant under axial stretch. Although the deformation-independence is verified for the Bragg-type band gaps in this work, other types of band gaps, e.g. the one caused by local resonance, are expected to be deformation-independent as long as the Lagrangian-formulated wave equations are invariant.

Future work could be directed towards exploring SPCs of other types, as well as the more practical issue of the design and fabrication of polymeric materials that behave as per the suggested strain energy functions. For example, many soft rubbers such as silicone and neoprene [41] behave according to the neo-Hookean model when the stretch deformation is intermediate. A semilinear material could be envisaged as being realized in one dimension via, e.g. a linear spring-mass chain, but would certainly require significant efforts to synthesize in three dimensions. We anticipate that progress in polymer engineering is required to solve this challenge by designing polymers with the required microstructure. Additionally, it would be interesting to design SPCs with band gaps that are partially deformation-independent, namely, the band gaps are invariant for certain deformation modes or wave modes. Finally, it would seem appropriate for the sake of potential applications, to study other aspects of invariance of acoustic and elastic waves to deformation, whether induced by mechanical means or otherwise, e.g. electrical or magnetic.

## Data accessibility

There are no experimental data associated with this work. The simulation input files can be accessed through https://sites.google.com/site/puzhangsite/code.

## Authors' contributions

P.Z. and W.J.P. proposed the idea behind the work and designed the study; P.Z. carried out the research and drafted the manuscript; W.J.P. helped draft and edit the manuscript. Both authors gave final approval for publication.

## Competing interests

We declare we have no competing interests.

## Funding

The authors are grateful to the Engineering and Physical Sciences Research Council (EPSRC) for financial support via grant no. EP/L018039/1.

- Received November 23, 2016.
- Accepted February 24, 2017.

- © 2017 The Authors.

Published by the Royal Society under the terms of the Creative Commons Attribution License http://creativecommons.org/licenses/by/4.0/, which permits unrestricted use, provided the original author and source are credited.