RELATED APPLICATIONS

[0001] This application claims the benefit of the following four U.S. Provisional Patent Applications, all of which are incorporated by reference in their entirety:

[0002] 1) Serial No. 62/372,570, filed on August 9, 2016;

[0003] 2) Serial No. 62/438,637, filed on December 23, 2016;

[0004] 3) Serial No. 62/517,406, filed on June 9, 2017; and

[0005] 4) Serial No. 62/538,064, filed on July 28, 2017.

FIELD OF THE DISCLOSURE

[0006] The present disclosure relates generally to improved techniques for the use of metamaterials and acoustic lenses along with acoustic transducers in haptic- based systems.

BACKGROUND

[0007] A continuous distribution of sound energy, referred to as an "acoustic field", may be used for a range of applications including haptic feedback in midair, parametric audio and the levitation of objects.

[0008] By defining one or more control points in space, the acoustic field can be controlled. Each point can be assigned a value equating to a desired amplitude at the control point. A physical set of transducers can then be controlled to create an acoustic field exhibiting the desired amplitude at the control points.

[0009] For this control to be effective there must be as many transducers as possible in order to create enough degrees of freedom for the acoustic field to asymptotically approach the desired configuration. However, from a commercial perspective, using individually actuated transducers is an expensive solution, prone to transducer element failure and variability. From this perspective, reducing the number of elements is preferred, leading to reductions in acoustic capabilities. Being able to reduce the number and complexity of the electrically actuated components, while being able to retain acoustic versatility, is therefore advantageous.

[0010] One set of solutions to achieve these goals is the use of metamaterial, which is a material that has artificially designed structural properties that change its physical behavior in ways that may be generally unavailable in existing materials. Metamaterials are used to control and manipulate light, sound, and many other physical phenomena. The properties of metamaterials are derived both from the inherent properties of their constituent materials and from the geometrical arrangement of those materials.

[0011] Furthermore, to create haptic feedback, an array of transducers is actuated to produce an ultrasonic acoustic field that then induces a haptic effect. It has been discovered through simulations and testing that for best effect the transducers must be positioned in non-uniform or irregular patterns. While this can be to some extent achieved for positioning during printed circuit board (PCB) assembly, arbitrary positioning raises the cost of the assembly process and reduces the transducer density possible. If there is a method to achieve this positioning of acoustic sources without resorting to complex and costly assembly processes, then there is a commercial advantage in such a method. By producing the metamaterial such that an irregular pattern is expressed, this can be achieved without requiring this extra assembly cost.

[0012] Furthermore, by focusing at a point, acoustic waves with increased power are created and can be used to generate a constant force. The force generated is small but naturally non-contact. Further applications of the non-contact forces may also lead to commercially viable devices.

[0013] Current designed metamaterials generally incorporate only static structure. By including structures that can be dynamically switched between states through time, the properties of the metamaterial can be modified on the fly These changes can be affected such that the metamaterial and transducer system can behave as though the input signal has changed, without actually changing the input signal. This is equivalent to effectively modifying the acoustic wave after, rather than before, leaving the transducer element.

[0014] This allows for a further division of labor and the resultant efficiency and potential cost savings. This is achieved by using a single or few highly efficient and powerful sources to generate a (generally monochromatic) acoustic wave that is then modified by switching metamaterial elements which are integrated into the material. This behaves similarly to the a transducer array. In this way, together the acoustic source and each element of the metamaterial may work together as a transducing element, only now the waveform modifications are expressed by the actable metamaterial layer.

[0015] Accordingly, it is desirable to implement one or more of the above techniques to achieve the goal of more effective and efficient use of transducers in haptic systems.

BRIEF DESCRIPTION OF THE FIGURES

[0016] The accompanying figures, where like reference numerals refer to identical or functionally similar elements throughout the separate views, together with the detailed description below, are incorporated in and form part of the specification, and serve to further illustrate embodiments of concepts that include the claimed invention, and explain various principles and advantages of those embodiments.

[0017] Figure 1 is a metamaterial apparatus consisting of a composition of multiple metamaterials, in accordance with some embodiments.

[0018] Figure 2 is a two-dimensional simulation of an acoustic plane wave interacting with a metamaterial structure, in accordance with some embodiments.

[0019] Figure 3 is a diagram of a stacked meta-material configuration, in accordance with some embodiments. [0020] Figure 4 is a graph of 4 layers showing a differing non-linear discretization for the amplitude of a signal, in accordance with some embodiments.

[0021] Figure 5 is a graph, showing the results of sampling the complex plane of possible wave states with 8 layers of metamaterial, in accordance with some embodiments.

[0022] Figure 6 is a cut-away view of a completed model of a structure that addresses the issues of field of view and the removal or reduction of extra unwanted maxima or spurious images, in accordance with some embodiments.

[0023] Figure 7 is a further cut-away view of a completed model of a structure that addresses the issues of field of view and the removal of spurious images, in accordance with some embodiments.

[0024] Figure 8 is a view of a completed model of a structure that addresses the issues of field of view and the removal of spurious images, in accordance with some embodiments.

[0025] Figures 9A and 9B are two-dimensional mock-ups of a section of an acousto-optical system, in accordance with some embodiments.

[0026] Figures 10A, 10B and IOC are series of interlocking geared elements designed to tessellate space, in accordance with some embodiments.

[0027] Figures 11 A and 11B are cross-sections of a geared element showing an acoustic meta-material element, in accordance with some embodiments.

[0028] Figure 12 is another series of interlocking geared elements designed to tessellate space, in accordance with some embodiments

[0029] Figure 13 is yet another series of interlocking geared elements designed to tessellate space, in accordance with some embodiments.

[0030] Figures 14A, 14B and 14C are yet another series of interlocking geared elements designed to tessellate space, in accordance with some embodiments

[0031] Figure 15 A, 15B and 15C are yet another series of interlocking geared elements designed to tessellate space, in accordance with some embodiments. [0032] Skilled artisans will appreciate that elements in the figures are illustrated for simplicity and clarity and have not necessarily been drawn to scale. For example, the dimensions of some of the elements in the figures may be exaggerated relative to other elements to help to improve understanding of embodiments of the present invention.

[0033] The apparatus and method components have been represented where appropriate by conventional symbols in the drawings, showing only those specific details that are pertinent to understanding the embodiments of the present invention so as not to obscure the disclosure with details that will be readily apparent to those of ordinary skill in the art having the benefit of the description herein.

DETAILED DESCRIPTION

[0034] Described herein are certain techniques for improving the performance of haptic-based systems. These techniques may be used either alone or in various combinations.

[0035] I. DYNAMIC CONTROL OF METAM ATERI AL S

[0036] A. Dynamic Metamaterials

[0037] Though there are many structures that qualify as metamaterials, the most common is that of an arrangement of elements whose size and spacing is much smaller relative to the scale of spatial variation of the exciting field. With this limit, the responses of the individual elements, as well as their interactions, can often be incorporated (or homogenized) into continuous, effective material parameters; the collection of discrete elements is thus replaced conceptually by a hypothetical continuous material The advantage of this homogenization procedure is that sophisticated and complex materials can be engineered with properties beyond what is naturally available. For these purposes (creating acoustic fields) the propagation of the monochromatic wave through the material is the crucial property that requires manipulation. By changing the

monochromatic wave, the phase offset of the acoustic wave is changed.

[0038] One such example of an acoustic metamaterial is a tortuous path material shown in Figure 1. Shown in Figure 1 is a metamaterial apparatus 100 consisting of a composition of multiple metamaterials 130. Figure 1 further shows an inset 170 of an exemplary metamaterial cell 160 showing a tortuous or elongated path that an acoustic signal may take. By forcing the wave to either follow a shorter or longer path through the metamaterial cell 160, the propagation of the wave through the metamaterial cell 160 may be allowed straight through or impeded through the introduction of the longer path. As an example, the overall width of the metamaterial cell 160 may be 13 mm wide 110, with each of three sub-cells being 4.4 mm wide 150, 4.3 mm wide 120, and 4.4 mm wide 140.

[0039] Generally, acoustic metamaterial cells may be built to introduce different amplitude attenuations and phase offsets in the waves travelling through them. This can be interpreted as the multiplication of the wave phasor by a complex- valued transmission coefficient.

[0040] Figure 2 shows a two-dimensional simulation of an acoustic plane wave interacting with a metamaterial structure where the overall effect is to generate a negative refractive index. A plane wave 215, interacts with the acoustic material 240, to produce an exiting acoustic wave 220 in the opposite direction to that which would be expected of a material whose refractive index is positive. This is at odds with any natural bulk materials, and is due to the tortuous path

metamaterial shown in which the acoustic waves follow different length paths to produce a subdivided plane wave which exits the material with different phase offsets. These interfere to produce the unnaturally refracted wave.

[0041] In a further step the acoustic metamaterial may then be made

reconfigurable, by mechanically, electrically or magnetically actuating structural elements inside the cells to dynamically change their behavior. This results in different complex transmission coefficients being dynamically expressed in one or more alternate switchable states. In the example of the tortuous path metamaterial, the "teeth" in the cells may be actuated; moved to make the path shorter or longer. This results in a dynamically reconfigurable phase change while adding little in the way of amplitude change. An element produced in this way would change phase without changing amplitude in a significant way. A cell containing two layers of mesh that cross to provide a larger or smaller proportion of open space may describe a simple acoustic metamaterial that implements a different element that results in a controlled amplitude drop. Notably due to physical limitations, a manufactured implementation of such an element will likely result in change in both phase and amplitude for any given input frequency, regardless of the initial configuration of the material.

[0042] In discrete components, a Helmholtz resonator is a container of gas (usually air) with an open hole (or neck or port). A volume of air in and near the open hole vibrates because the elasticity of the air inside behaves in a way analogous to a spring. A common example is an empty bottle: the air inside vibrates when blown across the top, Similarly, the metamaterial may also be constructed with a Helmholtz resonator, which could be actuated by opening a hole in a chamber producing effectively an amplitude increase.

[0043] These actuated metamaterials enable the one or few basic transducer element(s) to be simple, with no amplitude or phase change required. The dynamic metamaterial achieves the modification of the incident acoustic wave into the phase offset and amplitude required to drive control points, points in the acoustic field that provide controlled pressure changes as the result of interfering waves. This frees the transducer element or elements to be electrically driven in a simple way, potentially even being very highly resonant. This is important because resonant devices are highly efficient, but resist changes in output vibration and would be slow to respond and difficult to control. In this case however, control is achieved externally to the transducer, rendering this difficult control unnecessary. As a result, the transducer may be built to be as resonant as possible, as the effect of increasingly long response times becomes less important.

[0044] B. Designing a Metamaterial [0045] In order to create these effects, a monochromatic acoustic source is used, with the phase of the wavefront of the acoustic field manipulated. In existing devices this is achieved by manipulating the signals to individual transducer elements, which involves complicated processing and involved electronic design and fabrication. A method to apply the necessary phase manipulation without involved transducer signaling is commercially desirable.

[0046] Acoustic metamaterials may be designed using structure that is of the same scale as the features of the waves it is interacting with to take advantage of diffraction and interference behavior. These are then designed in order to create novel controlled effects. Ultrasonic waves in air are about the millimeter scale, which puts acoustic metamaterials in this regime easily in the range of simple manufacturing processes through plastics or other materials. These may then be used to influence the acoustic field around the object. In these situations, it is not necessary to place the acoustic source next to the metamaterial. The metamaterial may thus collect and focus incoming acoustic radiation passively into a region of space in order to do useful work such as a power source.

[0047] To design an effective metamaterial, the impedance of the material structure and its efficiency (the proportion of energy that is transmitted versus energy that is reflected) must be considered. One example of a metamaterial that has been designed and successfully fabricated is the tortuous path or labyrinthine structure. This provides close impedance matching to bulk air, which allows for high efficiency in the acoustic output.

[0048] C. Creating Foci

[0049] Acoustic metamaterials may be designed using these tortuous paths or labyrinthine features to manipulate the phase of incoming monochromatic sound. If the monochromatic sound is ultrasound and is amplitude-modulated at a frequency or group of frequencies that may be detected haptically and if it is focused to one or more points in space at once, these can be perceived by the sense of touch. If it is focused and amplitude-modulated at a frequency or group of frequencies that can be detected audibly, it can also be perceived by the sense of hearing as parametric audio. If it is focused and the time-averaged forces involve control, it can be used to levitate objects. A control point may allow for any combination of the foregoing. The acoustic metamaterials may change the phase of the wavefront and the frequency of the ultrasound is greater than any that might be contained in the modulated signal (so that the time lag involved in the metamaterial leaves them mostly unaffected). Then, the acoustic metamaterial may be designed in such a way as to create a focusing of the sound field, which may be off-axis and not symmetrical. This then simplifies the design of the emission of the sound field to a single element that can produce any wave using a single simple signal, which can then be corrected through the application of the metamaterial to create such a focus.

[0050] Alternatively, the array that controls individual transducer elements may be kept and the metamaterials then designed to accomplish the same function remotely. This modifies the acoustic field in the volume efficiently and separately from the individual signal control of the array. Further, if the array control is aware of the behavior of the metamaterial, the array may be configured to actuate the metamaterial remotely. This is much like redirecting a shining torch to point its beam at a lens. This may allow a single transducer array to create a series of foci in different positions or enhance the efficiency of an array in a given region of space.

[0051] D. Actuated Metamaterials

[0052] Acoustic metamaterials can be actuated through mechanistically or electrically modifying their structure in order to dynamically change the properties of the acoustic field. By changing the configuration of 'teeth' in a tortuous path metamaterial, the phase change incurred on a transported acoustic wave can be modified. In this way, multiple states can be encoded into the metamaterial in order to produce, for instance, a focus that moves from one place to another, or foci that change position such as moving forward and backward.

[0053] E. Stacking Metamaterials [0054] These acoustic metamaterial s can be stacked in order to modify the wavefront by passing the sound through two or more metamaterials. For example, the first metamaterial may provide focusing, while the second provides what form of phase singularity to use to levitate and impart angular momentum onto an object, or levitate and trap an object such that it is stationary. A further example might be a metamaterial that generates a focus, and then a further metamaterial that then 'cloaks' the focus by returning the wavefront to a planar or other configuration, which may also involve focusing or levitation.

[0055] F. Stacked Dynamic Metamaterials and Digital Phase Control

[0056] Metamaterial cells may be formed into a plate or "meta-surface" that may then be individually designed to generate a particular complex transmission coefficient when a cell is actuated. Through the stacking of these plates into layers, channels are created through the three-dimensionally actuated system that dynamically controls phases and amplitudes of the waves moving through them. This may be modelled by essentially taking the combined product of the complex- valued transmission coefficients on each layer.

[0057] These layers may be arranged in such a way as to build a digital system with each potentially being a binary "bit", although more states are also possible. On each layer there is an "unconfigured" or "reference" state which is described with unit amplitude and the other state or states are then measured in reference to this "reference" state. For the sake of generality, this is the state with the least transmission loss or most transmission gain. It should be noted that this by no means implies that the transmission through this layer is lossless or that this is a state that does not apply phase change, it simply means that this is the state which is the closest to lossless transmission in this particular design. Thus, the capabilities of the device in terms of dynamic range make the most sense to be described with respect to this, so this state is 'normalized' and assigned a coefficient of unity. By expressing each layer in an 'on' or 'off manner, bitwise digital control of phase becomes plausible through the stacking of plates that embody phase differences corresponding to complex values that are the different principal roots of unity for the powers of two. This is an exponential of a purely imaginary value, e ^l . A potential desired configuration is shown in Figure 2. Arbitrarily, many of the more precise layers may be omitted to provide lower complexity at the expense of reduced control.

[0058] Figure 3 shows a diagram 300 of a potential stacked meta-material configuration. Along the top 301, 303, 305, 307, 309, 311, 313, 315, 317, 319, 321, 323, 325, 327, 329, 331 are the possible complex valued transmission coefficients expressed in polar (r, Θ) addressable from a single switchable layer of amplitude attenuation, when followed by several layers of phase modification underneath. For each, the r corresponds to the change in amplitude of the exiting wave with respect to the incoming wave, while the Θ describes the phase angle offset of the exiting wave with respect to the incoming wave. The layers of phase modification can be expressed as added phase delay or added angle, or alternatively here, values that interact multiplicatively in the domain of complex numbers to give a variety of different resultant signals 302, 304, 306, 308, 310, 312, 314, 316, 318, 320, 322, 324, 326, 328, 330, 332. The numbers at the side 341, 343, 345, 347 give the complex-valued attenuation coefficients of each layer in turn, in both the Off and On' configurations. In the Off configuration, notionally the wave does not experience a change. Nonetheless, although this may be relative as the wave travels through one of these layers, the complex value representing the wave is multiplied by the left value, which is set to unity on all of the layers illustrated. In the On' configuration, the complex number representing the wave state is multiplied by the complex number labelling each layer on the right. These modify the relative phase angle or amplitude, modifying the state of the wave as it exits the layer. In this configuration, embodying binary phase control, the phase modification layers may each be expressed as the principal roots of unity of some power of two (although if the metamaterial plate may be switched among more than two states this may extend to the non-principal roots of unity of values beyond two, the roots of those roots of unity and so on). In an actualized design, however, an objectively defined modulus of the complex transmission coefficient may be less than one. This is because in a physical implementation, the absorption and reflection coefficients are likely to be nonzero and so the transmitted wave will have reduced energy and thus root-mean- square pressure with respect to the incoming wave.

[0059] These stacks may be specifically tailored to an application such as a focusing system that might only have metamaterial plates that implement phase modification. Or a dynamically focusing parametric speaker may have metamaterial plates that implement both phase and amplitude modification. If this form of acoustic metamaterial actuation is not suitable for changing the amplitude responsively, a hybrid system may be designed wherein elements may be changed in amplitude quickly (such as with the sliding mesh implementation of an acoustic metamaterial plate) and the phase component controlled at a slower rate by the mechanically-actuated tortuous path acoustic metamaterial. This would be advantageous in the case of haptic effects, in which moving the focus (which requires phase change and does not require high speed may be implemented) in a device which switches or updates slowly, while amplitude change (which requires high speed oscillation to generate low frequencies in the vibrotactile range) must be actuated at high speed.

[0060] G. Digital Amplitude Control

[0061] Amplitudes may also be dynamically controlled through the stacking of layers with discrete states, but the approach to find the optimal stacking arrangement has added complexity. This is because with each new amplitude attenuation layer, the transmitted amplitude is multiplied. (This differs from the phase angle where due to the complex multiplication the angles can be summed, and is thus simple to compute.)

[0062] Nonetheless, this multiplication can be viewed as a linear combination in a logarithmic space. A summation in logarithmic space for a unit complex value or exponential of a purely imaginary is the same as a linear summation of phase angles, as can be seen from Euler's formula, e' ^e = cos θ+i sin Θ. Considering the case of a ^-layered real- valued only metamaterial, each layer of which expresses one binary digit of amplitude control, there is an optimal set of transmission values that will produce the largest range of amplitudes. This optimal solution for layer m = 0, ... , n— 1 of n layers with binary states can be described by a real- valued exponential e-^ ² ™ ^-1) * = 1 — e ^~x where x can be found by solving _e -(2 ⁿ -i)x _ i _ _e -x^ f _{or rea} j _so i _u ti _on Thi _s determines the x value for which the bit string that describes the layer configuration that represents only the least significant reduction in amplitude with respect to the reference state; (00..01) is the same distance from unity as the layer configuration that represents the most significant reduction in amplitude (11..11) is from zero. These bit strings then correspond to the Off or On' configuration state of the metamaterial layers or plates described previously. By turning each to either the Off or On' state in the prescribed manner, the amplitude level can be tuned to be as expressive as possible given the number of layers available. This expressiveness results from the set of x values, found from the above equation, which are designed into the metamaterial plates for an optimum distribution of metamaterial controllable amplitude. An equivalent equation can be created for any combination of layers with any number of states. This gives the set of layer coefficients that for a given n best cover the unit interval. For example, when n=4 each filter would express a real- valued transmission coefficient corresponding to e ^~x , e ^~2x , e ^~4x and e ^~8x where the value of x can be found by solving e ^~ ^ ² * ^~1>x = 1 — e ^~x . This gives the optimal four acoustic metamaterials with real transmission coefficients values of 0.8720, 0.7603, 0.5781 and 0.3341 respectively. The solution to this equation for x ensures that the range of amplitudes producible by the system is configured to be the optimal possible spread between 0 and 1, in this case with a maximum gap of 0.1280 between 1 and the largest value and between 0 and the smallest value. This yields the most expressive possible stacked metamaterial system for a given design complexity.

[0063] Turning to Figure 4, shown is a graph 400 of how different numbers of layers («=1,2,3 and 4) each give a differing non-linear discretization for the amplitude. The x-axis 404 is the normalized bit pattern value, whose individual bits represent the On' or Off switching state of each of the metamaterial layers. It is also a non-linear distribution of amplitude values; "non-linear" due to the multiplication of the attenuation coefficients on each layer. The y-axis 402 is the output amplitude. This is correlated to the input bit pattern to the layer state, which may also be interpreted as the inversion of an input intensity value that would then map onto output amplitudes in a non-linear way. The coefficient of these layers of amplitude attenuation have in each case been chosen by evaluating the above equation for 1, 2, 3 and 4 layer(s). The curves in the figure show the inversely correlated relationship between the input bit pattern interpreted as a number, and the resulting possible levels of output amplitude, for one layer 406, two layers 408, three layers 410 and finally four layers 412 of actuated amplitude attenuation metamaterials.

[0064] H. Digital Complex Control

[0065] The complex transmission coefficients generated in the actuated and un- actuated positions of the acoustic metamaterial do not have to be purely real or purely imaginary. They may very likely be a combination of real and imaginary parts due to material limitations, acoustic impedance and physical effects that cannot be completely eliminated. In this case, an acoustic metamaterial can be designed to have a fully complex transmission coefficient with respect to the reference state. To achieve this, the approach to generate the real exponential can be merged with the approach to generate the imaginary exponential. For any real k, after having solved for the real part as before, for each layer m a complex

₂ m( _x ! ^2nki )

coefficient of e V 2™ ^y ean be found. These fully complex exponential solutions correspond to values that can be substituted in exactly the same way as before to create a «-ary hierarchy of complex exponentials which can be used to produce a spiral of sample points across the complex plane, sampling uniformly both amplitude and phase in combination. Alternatively, the phase can be sampled using an irrational step by using a coefficient like e ^-2 ™ ( ⁺ ( ^3- ^ ^π _{? in} which each step is separated by the golden angle. [0066] Turning to Figure 5, this less structured approach to sampling are compared with the uniform sampling for n=8. Shown is a graph 500 with the x- axis 504 is the real part of the complex transmission coefficient and the y-axis 502 is the imaginary part of the complex transmission coefficient. Shown on this graph are the results of sampling the complex plane of possible wave states with 8 layers of metamaterial. The input bit interpretation is not shown, but is implied by the positions of the points, as the reference point where all are zero or in the Off state always has a product which produces a complex coefficient of 1. Two uniform samplings of the phase and amplitude space are shown, one in which the system traverses the polar 2π radians in angle once 508 and the other involves traversing the polar 2π radians eight times (16π) 510. The last traversal involves jumping in phase by an irrational angular value 506 in order to better sample the area of the complex plane. This results in increased sampling efficiency of the space of possible exit waves at the expense of a complicated lookup table or function. Such a result occurs when implementing the complex value to bit string conversion to actuate the device into the state required to produce a wave with a given complex coefficient.

[0067] By creating a set of layers using one of these complex solutions for x, non- unity transmission and other losses in efficiency can be factored into the acoustic metamaterials. This occurs while continuing to result in an efficient sampling of the complex wavefunction space by the digitally actuated materials.

[0068] I. Specific applications for metamaterials in haptic systems

[0069] Thus, actuated stacked metamaterial layers may accomplish one or more of the following:

[0070] 1) Mechanically actuated teeth in tortuous path designs to turn on/off phase shifts;

[0071] 2) Mechanically actuated meshes which enable the control of effect open area leading to amplitude control; and [0072] 3) Mechanically opening of holes to control a Helmholtz resonator in a layer.

[0073] These layers do not have to encode a phase shift or an amplitude change - they may encode both in a complex phasor.

[0074] Alternatively the layers may be arranged such that they sample a particular line or curve in the complex plane - this can include both purely phase and purely amplitude curves.

[0075] Alternatively the layers may be arranged such that that they sample the complex plane of potential phasors efficiently using a logarithmic spiral of coefficients reachable through digital control.

[0076] Alternatively a system may conduct a sampling of the complex plane through an irrational step such as the golden angle in order to find the most efficient packing of coefficients in the complex plane.

[0077] Alternatively the system disclosed above may be used to create a haptic effect and/or parametric audio.

[0078] Alternatively the system disclosed above may be used for single or multiple control points using an acoustic lens for creating haptic effects and/or for creating parametric audio and/or for causing objects to levitate.

[0079] Alternatively the system disclosed above may be used for non axis- symmetric focusing or manipulations with acoustic lenses.

[0080] Alternatively the system disclosed above may be used for the design of impedance-matched metamaterial to maximize transmitted energy using tortuous paths and labyrinthine structures

[0081] Alternatively the system disclosed above may use an array underneath to provide input to metamaterial structure. The metamaterial structure provides acoustic "lenses" and "mirrors" in space that can be structurally separate to the reconfigurable acoustic source of a transducer array. Both "lenses" and "mirrors" can be reconfigurable through phase and amplitude changes. [0082] Alternatively, the system disclosed above may use metamaterials that are stacked. These metamaterials do not have to be physically touching for their properties to be concatenated. The layers in such stackings may have or be configured to have, distinct properties, such as: imparting angular momentum through phase manipulation; focusing; and gripping through phase singularity.

[0083] II. ACOUSTIC WIDE-ANGLE LENSES AND TUNABLE

RESONATORS IN HAPTIC SYSTEMS

[0084] To create haptic feedback, an array of transducers is actuated to produce an ultrasonic acoustic field that then induces a haptic effect. It has been discovered through theory, simulations and testing that the transducers must be positioned in non-uniform or irregular patterns for best effect. While this can be to some extent achieved for positioning during PCB assembly, arbitrary positioning raises the cost of the assembly process and reduces the transducer density possible. If there is a method to achieve this positioning of acoustic sources without resorting to complex and costly assembly processes, then there is a commercial advantage in such a method.

[0085] The angles and directions of the acoustic transducers also play a large part in determining the haptic interaction space, which is a function of the carrier frequency and the field of view of the transducer elements. During PCB assembly, precision angling is unrealistic, constraining a straightforward design to have only a limited interaction space that will become smaller as the ultrasonic carrier frequency increases. An approach to enable this to happen at low cost would be advantageous.

[0086] Further simulations have shown that Helmholtz resonators placed over the transducers can increase the output amplitude or efficiency, but only for a fixed frequency dependent on the size and shape of the resonating chamber. Since during the normal creation of haptic effects this frequency can vary in some band around the ultrasonic carrier, it could be advantageous to be able to modify dynamically the resonant chamber to effect substantial power savings. [0087] A. Array Coverings

[0088] To modify the configuration of a transducer array which has been cost- effectively fabricated in a uniform, potentially rectilinear grid or other repeating pattern on a PCB, a structure may be overlaid that has been built that specifically reconfigures the transducer array into a new configuration that is not constrained by the PCB layout. This structure consists of a series of pipes which can redirect the acoustic waves to be emitted from a new given position and in a new given direction. Each pipe must be a close to straight as possible and, in the cases where it is not straight, must be smoothly curved. It must also have prescribed entry and exit points, as well as entry and exit tangents appropriately normal to the wavefront normal direction of the entering and exiting sound waves.

[0089] In order to accomplish this, the footprint of the rectilinear grid of transducers on the PCB to modify must first be decided or found. This necessarily informs the number and distribution of desired acoustic output points and directions that it is possible to create. This is then synthesized with the intention of reconfiguring the array from a rectilinear grid on the PCB (for example) into a different configuration with more desirable properties, such as increased field of view, reduced extra maxima or both. Also, it may be costly, difficult or impossible to produce this in a PCB design. Next, the center point of each physical transducer element with a emission direction is determined, which is then paired with a corresponding output element location in the structure. Next, a squared minimum distance criterion is used to evaluate the fitness of transducer- output pairs, swapping the outputs if this is found to increase fitness (decrease the summed squared minimum distances between the physical transducer input and output). Then a curvature-minimizing spline or other curvature-minimizing curve is solved for and is drawn between the transducer center point and the output point, wherein the endpoints are constrained by the tangent vectors of the transducer emission direction and acoustic output direction respectively. As the pipes are required to be a certain diameter, this is expressed as a minimum distance constraint among the set of created curves and a further step may try to enforce a minimum distance if possible. Finally, the curves and their pipe diameters are carved geometrically from a solid shape appropriate to the exit position and directions given to produce a design for fabrication.

[0090] Turning to Figure 6, shown is a cut-away view of a completed model of a structure 600 that addresses the issues of field of view and the removal or reduction of extra unwanted maxima or spurious images. Inside the 3D model of the structure the pipes 610 that serve to direct the acoustic output are visible. In the lower part of the image the entrance holes to the pipes are visible, taking output from a grid of transducers which have been arranged or fabricated rectilinearly 620.

[0091] Turning to Figure 7, shown is a further cut-away view of a completed model of a structure 700 that addresses the issues of field of view and the removal of spurious images. Now visible are the exit holes from the pipes 710. These are arranged such that their non-regular arrangement generated by irrational arithmetic serves to scatter energy away from any particular location by replacing co-linear arrangements of acoustic output locations. The pipes in the cut-away remain visible 720 and the line of entry points into the structure are also visible in the lower part of the image 730.

[0092] Thus, Figure 7 shows the use of non-uniform output location

arrangements; figures 6 and 7 show the tangents on the curvature minimizing pipes; and figure 6 shows the rectilinear uniform arrangement of the pipe inputs.

[0093] Turning to Figure 8, shown is shown is a view of a completed model of a structure 800 that addresses the issues of field of view and the removal of spurious images. The completed irregular structure of outlet holes is shown 810 following the surface of a hemisphere, although this should not preclude a less extreme example of a smaller adjustment to the addressable field of view (in the form of an array covering of this kind with reduced pipe length, thickness and curvature) from the transducers in which control point can be generated.

[0094] B. Tunable Helmholtz resonators [0095] As can be seen from Figure 6, there is a portion of the tube 610 that moves immediately away from the PCB. This can be readily converted into a Helmholtz resonator with a little extra material. The frequency that a Helmholtz resonator amplifies is tied closely to its volume and the size of the open holes in its surface. By increasing or decreasing the area of the open holes, the frequency that the resonance amplifies is changed to be higher or lower respectively. The same effect can also be applied to the volume contained in the Helmholtz resonator, by changing the volume the resonance can be changed. One way to achieve this is to add an extra piece of moving structure to mechanically open and close the holes. To keep the system cost-effective however, the volume of the chamber can instead be modified by actuating further membranes bounding the resonator separately to the ultrasonic transducer. As the frequency shift required for considerable efficiency savings when creating a haptic effect is of the order of a few percent in either higher or lower frequency, this generates a much lower power design since the amount of energy required by the extra transducer is much smaller.

[0096] III. ACOUSTIC RADIATION FORCE FOR OPTICAL MODULATION THROUGH PHASED ARRAYS

[0097] In an acoustic field, control points may be defined in the field that lead to the creation of tightly controlled pressure at these locations.

[0098] These effects are mediated through the creation of small acoustic radiation forces, which is a non-linear effect generated by acoustic waves. For example, lenses and mirrors may be made of a variety of different materials. Each of these has mechanical properties that governs its response to incident forces. Due to ease of shaping, lenses and mirrors may be made of liquids and gels on top of more traditional materials. Creating a non-contact acoustic radiation force provides such a small incident force that it can be used for the purposes of adaptive optics.

Modifying the force over time may allow lenses or mirrors to change over time, allowing optics to be synchronized to a clock or in a way that captures more data about the scene of interest. [0099] For example, light field cameras can use an array of small lenses configured differently to obtain information about the directions of incoming light. Using an acoustic radiation force to quickly modify the lens in a camera in an undulating pattern would enable light directions to be extracted as the sum of component light pixels recorded from predictable positions and directions.

Computational photography techniques may then be used to extract the resulting data.

[00100] Using acoustic radiation force in this way allows gel and liquid lenses to be reshaped in situ into a variety of different forms suitable for a range of different tasks. They can be made convex or concave depending on the distribution of control points and places where the pressure (and thus the acoustic radiation force) is controlled in the field. By placing these into a distribution, the lens surface may be totally reshaped in a way that enables complete control over the surface geometry. In gels or liquid lenses dominated by surface tension, this would enable further freedom over existing techniques for adjusting lens geometries. Further, given that the acoustic phased array is a separate and completely contactless device, it simplifies the design and manufacture of the lenses or mirrors in question. This enables the phased array to be serviced or replaced separately to the optical elements, especially as technological progress may render obsolete either component over time.

[00101] One such design of a system may place the lens in the center of an annular array or, alternatively, inside a tube of transducers. Shown in Figure 9A is a two-dimensional mock-up of a section of the acousto-optical system 900.

Shown is acoustic field 910, actuated transducers 930 a lens 940 and a single control point 920. The single control point 920 does not interact with the lens 940. Shown in Figure 9B is a two-dimensional mock-up of a section of the acousto-optical system 950. Shown is acoustic field 960, actuated transducers 980 a lens 990 and a single control point 970. Here the single control point 920 containing high pressure is lowered towards the lens 990. As a result the lens 990 deforms slightly (effect greatly exaggerated). This can be used to modify the curvature of a lens or mirror surface without touching them allowing for precise non-contact geometry control in complex systems in which more traditional approaches would involve localized electronics, adding unwanted complexity to existing assemblies of optical elements.

[00102] Further, existing designs of lens systems could recess acoustic elements that can work together in a similar configuration in any system of lenses, such as telescopes, microscopes and other optical systems to provide further non-contact calibration.

[00103] IV TIME-HARMONICS MET AM ATERI AL S

[00104] A. Introduction

[00105] The optimal conditions for producing an acoustic field of a single frequency can be realized by assigning activation coefficients to represent the initial state of each transducer. But to create haptic feedback, the field can be modulated with a signal of a potentially lower frequency. For example, an acoustic field of 40 kHz may be modulated with a 200 Hz frequency to achieve a 200 Hz vibrotactile effect.

[00106] Meta-materials are materials that have artificially designed structural properties that change their physical behavior in ways that may be generally unavailable in existing materials. They can be artificially structured to control and manipulate light, sound and many other physical phenomena. Specifically, meta- materials of interest here are those pertaining to the manipulation of acoustic waves. They may be engineered to take a monochromatic acoustic wave and apply phase shifts or delays through modification to the internal structure, enabling such a material to generate a pre-defined focus or other modification to the resulting acoustic field.

[00107] Modifications to meta-materials may be actuated electronically to provide full re-configurability. This enables the device to be further organized by splitting the production of the field into separated generation and control components, potentially yielding benefits in efficiency. But even this requires a per-channel actuation or more, necessitating an extensive and thus expensive array of electronics and components.

[00108] As the haptic effects that are to be reproduced in mid-air consist primarily of synchronized repeating structures of phase and amplitude through time, this suggests that a meta-material may be designed to fall in between a static meta- material and a fully actuated and re-configurable meta-material. In other words, this is a material that implements phase and amplitude change through purely mechanical means. This would be beneficial because it could embody cyclical haptic effects without requiring complex and potentially expensive electronics.

[00109] B. Design

[00110] A time harmonic meta-material is repeated cyclically in the same manner while changing the acoustic field state continuously. Since time harmonic systems can be parameterized in angle, it is therefore reasonable that this achieved by rotation. Since this must also be synchronized, the angular position of the elements with respect to each other must be actuated only in lockstep, potentially with a fiducial to show the zero angle location in the meta-material to aid assembly. This leads to the consideration of a design involving a plethora of interlocking and rotating structures.

[00111] In some instances, such as in intentionally non-uniform configurations of acoustic elements designed for the reduction of unwanted geometrical artifacts, it is preferable not to pack the structures too tightly. But since the elements must touch and thus cannot be too sparse, this results in a roughly dendritic structure in order to propagate the mechanical motion. To build an efficient metamaterial, similarly to the electronically driven transducers needed to build a traditional array, these rotating structures must ideally be packed as closely together as possible. When close packing of rotating mechanical elements is required, it is reasonable to consider designing them as interlocking geared structures. As these geared structures must interlock, it is prudent to consider alternating clockwise and anti-clockwise rotating structures necessary to tessellate a two-dimensional space without generating a jammed structure that cannot rotate freely. This structure may also instead tessellate space while being non-uniform, through the utilization of a quasi-crystal structure. Each of these structures must have the phases and amplitudes that are to be produced encoded into the structure such that it is expressed as the structure rotates. This suggests that multiple copies of the phase and amplitude behavior may be encoded in the 2π angle rotation around the individual meta-material element.

[00112] A further restriction on the interlocking elements is therefore that the rotation of these elements must parameterize time using an angle. This way, the speed of angular rotations must— in addition to the requirement of unjammed interlocking gear elements— also have an integer ratio to the periodic acoustic pattern. This implies that in the case of only one copy of the periodic field behavior being encoded into the structure, only one angular speed and thus circumference is permissible in order that the structures stay synchronized in their time harmonic behavior.

[00113] But if the behavior of the acoustic field is encoded multiple times around the edge of such an element, this unlocks the possibility of having different circumferences. This works if the ratios of the circumference is matched by the number of times that the periodic behavior is encoded in the 2π angle rotation of any individual meta-material element. This implies that the diameter of the interlocking circular or in three dimensions cylindrical elements must be either the same circumference or have circumferences that are integer ratios.

[00114] Shown in Figure lOA is a series 1000 of interlocking geared elements 1010, 1020, 1030, 1040 so designed as to tessellate space. These gears must travel in opposite directions (compare 1010 and 1020 with 1030 and 1040). This is the simplest design but it causes a large percentage of open area.

[00115] Figure 10B shows a further configurations 2010 that tessellate two- dimensional space using larger interlocking gears 2012, 2016, 2030, 2036 and smaller interlocking gears 2014, 2018, 2020, 2022, 2024, 2026, 2028, 2032, 2034. This arrangement will operate properly if the circumferences of the larger interlocking gears 2012, 2016, 2030, 2036 and smaller interlocking gears 2014, 2018, 2020, 2022, 2024, 2026, 2028, 2032, 2034 have integer ratios. Further, the larger interlocking gears 2012, 2016, 2030, 2036 travel in opposite directions from smaller interlocking gears 2014, 2018, 2020, 2022, 2024, 2026, 2028, 2032, 2034.

[00116] Figure IOC shows a further configurations 2050 that tessellate two- dimensional space using larger interlocking gears 2052, 2058, 2068, 2074, 2080 and smaller interlocking gears 2054, 2056, 2060, 2062, 2064, 2066, 2070, 2072. This arrangement will operate properly if the circumferences of the larger interlocking gears 2052, 2058, 2068, 2074, 2080 and smaller interlocking gears 2054, 2056, 2060, 2062, 2064, 2066, 2070, 2072 have integer ratios. Further, the larger interlocking gears 2052, 2058, 2068, 2074, 2080 travel in opposite directions from smaller interlocking gears 2054, 2056, 2060, 2062, 2064, 2066, 2070, 2072.

[00117] The acoustic elements themselves are composed of a series of stationary teeth on the inner edge of the central spindle, and the rotating teeth structure moving with the mechanism around the outside. As the system rotates, the channel in which the acoustic waves can travel increase and decreases in length as well as changing in amplitude according to the meshing of the internal geometry of the meta-material with the static structure. This gives rise to changing phase and amplitude yielding the desirable properties of the time harmonic meta-material. This may also be achieved through a cam so designed as to have an effect similar to that of the teeth. The order of motion may also change in an implementation of this system; the inner spindle may be rotating while the outer part remains stationary.

[00118] Shown in Figure 11 A is a diagram 1150 of a cross-section 1153 of geared element 1158 highlighting a stationary central spindle 1152 mounted on the common plate 1154 across the array. The acoustic waves travel 1192, 1194 adjacent to the highlighted standardized surface 1152 that the rotational element moves past via rotational movement 1159, 1190. This is designed to allow the rotating spindle options to deflect the acoustic waves differently using different patterns to produce a variety of phase and amplitude changes. The mesh 1165, 1156 is an optional acoustically transparent material.

[00119] Shown in Figure 1 IB is a diagram 1200 of a cross-section 1153 of a geared element 1158 highlighting outer meta-material structure 1160a, 1160b. The acoustic waves travel 1192, 1194 adjacent to the highlighted meta-material structures 1160a, 1160b that the rotational element moves past via rotational movement 1159, 1190. Each gear may have differently shaped meta-material structures 1160a, 1160b that serves to encode a meta-material surface In concert with the stationary inner edge, this provides a path for the acoustic waves that may be longer or shorter. Alternatively, this may attenuate the waves to produce known and controlled changes in phase and amplitude throughout the rotational period of this element. The mesh 1170a, 1170b, 1170c, 1170d shown is an optional acoustically transparent material that also functions as structural support in some embodiments.

[00120] Shown in Figure 12 is the design 1100 of Figure 10A, while also showing the open portions 1110, 1120, 1130, 1140 at the base of the structures 1010, 1020, 1030, 1040. As shown in Figures 11 A and 1 IB, a monochromatic acoustic wave is intended to enter the meta-material, pass through the metamaterial structure on the interior of the geared elements, and leaving the structure with a modified phase. The motion of the structure 1015, 1025, 1035, 1045 generates a small frequency shift in the acoustic wave which translates into the desired motion in the air. The system may use a larger or smaller opening depending on the permissible smoothing of the function of amplitude and phase. This is due to parts of the wave moving through meta-material structure intended for adjacent points of the function in time. While this appears to be a small opening, the meta-material may (due to the only requirement being a monochromatic input wave) be sitting on a reverberant cavity that does not permit the energy to leave except through the openings indicated.

[00121] In Figure 13, shown is a series 3000 of interlocking geared elements 3010, 3020, 3030, 3040 so designed as to tessellate space. Each of the geared elements 3010, 3020, 3030, 3040 have an opening 3080, 3082, 3084, 3086 for the monochromatic acoustic wave. The first geared element 3010 has four radial hatchings 3012a, 3012b, 3012c, 3012d, which corresponds to four different kinds of metamaterials to produce different parts of the harmonic period. The second geared element 3020 has four radial hatchings 3022a, 3022b, 3022c, 3022d, which corresponds to four different kinds of metamaterials to produce different parts of the harmonic period. To produce a different part of the acoustic field in space, these may be different from the first geared element 3010. The third geared element 3030 has four radial hatchings 3032a, 3032b, 3032c, 3032d, which corresponds to four different kinds of metamaterials to produce different parts of the harmonic period. Again, to produce a different part of the acoustic field, these may be different from the first 3010 and second 3020 geared elements. The fourth geared element 3040 has four radial hatchings 3034a, 3034b, 3034c, 3034d, which corresponds to four different kinds of metamaterials to produce different parts of the harmonic period. As this is to produce a fourth part of the acoustic field differing spatially from any of the aforementioned, this may again result in a differing internal meta-material structure. Each of the first radial hatchings 3012a, 3022a, 3032a, 3042a of the four geared elements 3010, 3020, 3030, 3040 represent metamaterial structure whose structurally encoded phase and amplitude are substantially expressed as the resulting acoustic field at time to + nT. Each of the second radial hatchings 3012b, 3022b, 3032b, 3042b of the four geared elements 3010, 3020, 3030, 3040 represent metamaterial structure which is expressed at time to + (n + ¼)T. Each of the third radial hatchings 3012c, 3022c, 3032c, 3042c of the four geared elements 3010, 3020, 3030, 3040 represent metamaterial structure expressed at time to + (n + ½)T. Each of the fourth radial hatchings 3012d, 3022d, 3032d, 3042d of the four geared elements 3010, 3020, 3030, 3040 represent metamaterial structure expressed at time to + (n + 3/4)T. In each case to represents a start time, T the period of each revolution of the gear structures and n any positive integer. [00122] As the structure rotates, different portions of the prebuilt meta-material are exposed to the monochromatic acoustic source and therefore produce different effects that are periodic in time. In this structure, the first and second geared elements 3010, 3020 rotate counterclockwise 3090, 3092 and the third and fourth geared elements 3030, 3040 rotate clockwise 3094, 3096. Thus, the respective radial hatchings (and corresponding metamaterials) 3012a, 3012b, 3012c, 3012d of the first geared element 3010 and the radial hatchings (and corresponding metamaterials) 3022a, 3022b, 3022c, 3022d of the second geared element 3020 are placed in opposite order than the respective radial hatchings (and

corresponding metamaterials) 3032a, 3032b, 3032c, 3032d of the third geared element 3030 and the radial hatchings (and corresponding metamaterials) 3042a, 3042b, 3042c, 3042d of the fourth geared element 3040.

[00123] Although four geared elements with four hatchings each are shown, there may be more or less geared elements and hatchings (and corresponding metamaterials) to produce a similar effect.

[00124] Figure 14A shows a further configuration 3100 that tessellate two- dimensional space using larger interlocking gears 3152, 3154, 3156, 3158, 3160 and smaller interlocking gears 3102, 3104, 3106, 3108, 3110, 3112, 3114, 3116. In this example, the larger interlocking gears 3152, 3154, 3156, 3158, 3160 have twice the circumference of the smaller interlocking gears 3102, 3104, 3106, 3108, 3110, 3112, 3114, 3116.

[00125] Figure 14B shows a detail view 3200 of a typical smaller interlocking gear 3208. The smaller interlocking gear 3208 rotates counterclockwise 3204 about spindle or hub 3206 and has one acoustic opening 3202. In contrast, figure 14C shows a detail view 3300 of a typical large interlocking gear 3310. The larger interlocking gear 3110 rotates clockwise 3304 about spindle or hub 3306 and has two acoustic openings 3302, 3308. The goal of this arrangement is to have the same acoustic field to be generated by both gear types. Here, the angular velocity of the larger interlocking gear 3310 is halved (since it has twice the circumference of the smaller interlocking gear 3208). Accordingly, the larger interlocking gear 3310 has two acoustic openings 3302, 3306 while the smaller interlocking gear 3208 has one acoustic opening 3302.

[00126] Figure 15A shows a further configuration 4000 that tessellate two- dimensional space using larger interlocking gears 4002, 4008, 4012, 4016, 4020 and smaller interlocking gears 4004, 4006, 4010, 4014, 4018, 4022, 4024, 4026. In this example, the larger interlocking gears 4002, 4008, 4012, 4016, 4020 have twice the circumference of the smaller interlocking gears 4004, 4006, 4010, 4014, 4018, 4022, 4024, 4026. Similar to figure 13, each of these interlocking gears incorporates radial hatchings (and corresponding metamaterials).

[00127] Figure 15B shows a detail view 4100 of a typical smaller interlocking gear 4106. The smaller interlocking gear 4106 rotates counterclockwise 4104a about spindle or hub 4102 and has one acoustic opening 4108a. The small interlocking gear 4106 incorporates four radial hatchings (and corresponding metamaterials) 4110a, 4120a, 4130a, 4140a. In contrast, figure 14C shows a detail view 4200 of a typical large interlocking gear 4202. The larger interlocking gear 4202 rotates clockwise 4104b about spindle or hub 4204 and has two acoustic openings 4108b, 4108c. The larger interlocking gear 4202 incorporates four radial hatchings (and corresponding metamaterials), each of which is split into two parts 4110b, 4110c, 4120b, 4120c, 4130b, 4130c, 4140b, 4140c. This is because-unlike between individual gears— within this larger interlocking gear like radial hatching must substantially create a like acoustic field and therefore a like meta-material configuration.

[00128] The goal of this arrangement is to have synchronized contributions to the same acoustic field to be generated repeatedly by both gear types. Here, the angular velocity of the larger interlocking gear 4202 is halved (since it has twice the circumference of the smaller interlocking gear 4106). Accordingly, the larger interlocking gear 4202 has two acoustic openings 4108b, 4108c while the smaller interlocking gear 4106 has one acoustic opening 4108a.

[00129] Further, each of the radial hatching (and corresponding metamaterial) in the smaller interlocking gear 4106 has two counterparts on the larger interlocking gear 4202, each on opposite sides of the other. The first radial hatching (and corresponding metamaterial) 4110a in the smaller interlocking gear 4106 is different to, yet expressed at the same time as the two first radial hatchings 4110b, 4110c (and corresponding metamaterials) in the larger interlocking gear 4202, which belong to the same gear element and thus must be equivalent. The second radial hatching (and corresponding metamaterial) 4120a in the smaller

interlocking gear 4106 is different to, yet expressed at the same time as the two second radial hatchings 4120b, 4120c (and corresponding metamaterials) in the larger interlocking gear 4202, which belong to the same gear element and thus must be equivalent. The third radial hatching (and corresponding metamaterial) 4130a in the smaller interlocking gear 4106 is different to, yet expressed at the same time as the two third radial hatchings 4130b, 4130c (and corresponding metamaterials) in the larger interlocking gear 4202, which belong to the same gear element and thus must be equivalent. The fourth radial hatching (and

corresponding metamaterial) 4140a in the smaller interlocking gear 4106 is different to, yet expressed at the same time as the two fourth radial hatchings 4140b, 4140c (and corresponding metamaterials) in the larger interlocking gear 4202, which belong to the same gear element and thus must be equivalent.

[00130] Further, each of the radial hatching 4110a, 4120a, 4130a, 4140a (and corresponding metamaterial) in the smaller interlocking gear 4106 is laid out in a set order. The counterpart radial hatchings 4110b, 4110c, 4120b, 4120c, 4130b, 4130c, 4140b, 4140c on the larger interlocking gear 4202 (and corresponding metamaterial) are laid out in the opposite order. This is because the smaller interlocking gear 4106 rotates 4104a in one direction and the larger interlocking gear 4204 rotates in the opposite direction 4104b.

[00131] By setting up the larger interlocking gear 4202 and the smaller interlocking gear 4204 in this manner (and then generalizing this to all the gears in Figure 15 A), the same acoustic field may be contributed to by both gear types. [00132] The foregoing embodiments may be generalized to any size gears with integer circumference ratios, any arrangement of such gears to tessellate space, and any amount of differing metamaterials.

[00133] V. CONCLUSION

[00134] The various features of the foregoing embodiments may be selected and combined to produce numerous variations of improved haptic systems.

[00135] In the foregoing specification, specific embodiments have been described. However, one of ordinary skill in the art appreciates that various modifications and changes can be made without departing from the scope of the invention as set forth in the claims below. Accordingly, the specification and figures are to be regarded in an illustrative rather than a restrictive sense, and all such modifications are intended to be included within the scope of present teachings.

[00136] The benefits, advantages, solutions to problems, and any element(s) that may cause any benefit, advantage, or solution to occur or become more pronounced are not to be construed as a critical, required, or essential features or elements of any or all the claims. The invention is defined solely by the appended claims including any amendments made during the pendency of this application and all equivalents of those claims as issued.

[00137] Moreover, in this document, relational terms such as first and second, top and bottom, and the like may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. The terms

"comprises," "comprising," "has", "having," "includes", "including," "contains", "containing" or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises, has, includes, contains a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. An element proceeded by "comprises ... a", "has ... a", "includes ... a", "contains ... a" does not, without more constraints, preclude the existence of additional identical elements in the process, method, article, or apparatus that comprises, has, includes, contains the element. The terms "a" and "an" are defined as one or more unless explicitly stated otherwise herein. The terms "substantially", "essentially", "approximately", "about" or any other version thereof, are defined as being close to as understood by one of ordinary skill in the art. The term "coupled" as used herein is defined as connected, although not necessarily directly and not necessarily mechanically. A device or structure that is "configured" in a certain way is configured in at least that way, but may also be configured in ways that are not listed.

[00138] The Abstract of the Disclosure is provided to allow the reader to quickly ascertain the nature of the technical disclosure. It is submitted with the understanding that it will not be used to interpret or limit the scope or meaning of the claims. In addition, in the foregoing Detailed Description, it can be seen that various features are grouped together in various embodiments for the purpose of streamlining the disclosure. This method of disclosure is not to be interpreted as reflecting an intention that the claimed embodiments require more features than are expressly recited in each claim. Rather, as the following claims reflect, inventive subject matter lies in less than all features of a single disclosed embodiment. Thus the following claims are hereby incorporated into the Detailed Description, with each claim standing on its own as a separately claimed subject matter.