Nonlinear resonators can have advantages over linear designs including increased sensitivity towards changes in their physical properties and environment, and high quality factors which make them attractive in applications such as mass/chemical sensors or signal filters. Designing nonlinear structures, however, requires much understanding of nonlinear behavior characteristics of structures. Similarly, the proliferation of 3D or additive manufacturing/printing capabilities has opened the doors to deploying nonlinear resonators on scales not possible earlier. However, to obtain consistent nonlinear dynamic performance the designer must perform a careful analysis to explore the existence and repeatability of desired nonlinear behavior. Also, the use of 3D printing with the associated substrate material properties poses its own challenges in regards to device simulation in view of the fact that most of the traditional literature on nonlinear resonators assumes linear material stiffness. In this chapter, the authors discuss computational design methods for structural design, and specifically study the case of 1:2 internal resonances in resonators made of nonlinear (hyperelastic) materials. The design methods allow for development of large number of candidate resonator designs without a required significant nonlinear structural design experience, and the study of the dynamic response of the resonators provides a glimpse in to the 1:2 nonlinear internal resonance exhibited by the candidate resonators.
- nonlinear dynamics
- internal resonances
- hyperelastic materials
- 3D printing
- topology optimization
The development of the micro- and nano-electronics industry coupled with the advances in semiconductor manufacturing techniques led to an interest in developing and applying micro- and nano-electromechanical systems (MEMS and NEMS) [1, 2]. Due to the small length scales of such devices and the popular modes of actuation employed by designers for MEMS or NEMS, such as electrostatic actuation, it led to the observations that nonlinear effects in many of these devices were the norm rather than the exception. This realization led to the study of effects of nonlinearities on the dynamic response of MEMS and NEMS in their various modes of operation, the effects in some cases being detrimental to linear design performance, and in some cases being beneficial to performance . The study of nonlinear dynamics is an area of research with long history in structural and mechanical systems . Several attempts have been made to incorporate the nonlinear dynamic effects and use the plethora of associated phenomena into operating mechanisms of MEMS devices particularly as mass and/or chemical sensors and filters . While such nonlinear devices have several advantages over linear ones with the same functionality in terms of measurement resolution, there have been challenges in making sure that the designed topologies are “in-tune” with the semiconductor manufacturing processes. At a conceptual level, researchers in nonlinear dynamics often work with lumped-parameter models  and many interesting applications have been considered [5, 6] for systems characterized by a single degree of freedom. For systems with more than one degree of freedom, one of the most compact representations is a system exhibiting a nonlinear 1:2 internal resonance, the spring-pendulum system [4, 7]. While this system lends itself very conveniently to a systematic analytical study , it is a relatively hard task to replicate it physically at micro- or nano-scales. This is even more so when the structural components fabricated involve two- or three-dimensional elastic structures.
3D printing or additive manufacturing offers several appealing advantages in terms of building devices based on nonlinear dynamic principles. The fabrication processes and dimensions are such that there is better repeatability with complex mechanical designs, as well as initial prototyping and low volume production costs come without the need for significant capital expenditure [9, 10]. Recently studies have been reported for prototyping nonlinear vibratory components made with 3D printers having a feature size of less than 1 mm . While the dimensions of polymeric resonators created by 3D printers may limit their measurement resolution and frequency range of operations, the manufacturing process itself is more repeatable for certain kinds of resonators. However, using 3D printing for producing nonlinear resonators also comes with its own challenges. While micro- and nano-resonators operate in an environment with many kinds of nonlinearities, 3D printed structures have to largely rely on only two sources of nonlinearities, namely, geometric and material nonlinearities unless controlled material variations or composite structures are explicitly introduced in fabrication. Fortunately, as was demonstrated by Tripathi and Bajaj [12, 13], both geometric nonlinearities due to finite deformations and material nonlinearities due to nonlinear hyperelastic properties of the 3D printed material can produce nonlinear dynamic effects such as 1:2 internal resonance.
A 1:2 internal resonance is a popular mechanism exhibited and employed by many nonlinear dynamics based resonators . Internal resonance in a structure refers to the energy transfer that occurs between two modes of the structure when their natural frequencies are almost commensurable and the structure has some appropriate nonlinearity. For example, for 1:2 internal resonance, if the two modes of a structure have their natural frequencies close to the ratio 1:2 and harmonic excitation of the higher mode is above a certain threshold, energy can be transferred from the resonant response of the higher mode to the lower frequency mode in the presence of quadratic nonlinearities. The mathematical description of a 1:2 internal resonance and the dynamics is well established [4, 7, 8]. For the purposes of evaluating the suitability of using 3D printing to produce resonators exhibiting 1:2 internal resonances, it is important to demonstrate that the dynamic response equations for the resonators exhibit the same mathematical characteristics as the cardinal examples.
In the context of elastic structures exhibiting various internal resonances, the present work focuses on elastic plate-type structures [4, 15]. A few representative works on different aspects of nonlinear vibrations of rectangular plates with internal resonances are [16, 17, 18]. In general, isotropic plates with different simple boundary conditions do not exhibit any commensurate frequencies unless there exists some type of symmetry of the structure. Some works have considered optimization of geometry and material distribution to affect frequency distributions as well as internal resonances [19, 20]. A systematic approach is based on the concepts in “topology optimization” . Applications of topology or shape optimization have now appeared in the literation on nonlinear dynamics as well, with the works in [22, 23, 24] focusing on general one-dimensional elastic systems whereas the works in [12, 13] focusing on plate structures with internal resonances. The overall goal is to tailor the system’s dynamic response to some desired form for appropriate external excitations.
In the present study, particular classes of resonator designs consisting of rectangular plates with cutouts which can be easily fabricated using 3D printing are analyzed for their nonlinear dynamic response. To obtain a suitable resonator design with commensurable (1:2) natural frequencies, a parametric optimization process which varies the sizes of the cutouts is employed. The natural frequencies themselves are computed using linear finite element analysis (FEA). The resonators are assumed to be made of a Mooney-Rivlin hyperelastic material  which is anticipated to provide the material nonlinearity necessary to produce 1:2 internal resonances. Once the optimization process is able to provide a candidate structure, the mode shapes obtained by the finite element analysis are used to build a reduced order model of the resonator displacements. This displacement field can then be used to derive the kinetic and strain energies of the structure which can provide the system Lagrangian. This Lagrangian is then averaged and subjected to the Euler-Lagrange conditions to derive the slow-amplitude equations of motion of the structure that provide the dynamic steady state response.
This work has two following sections: Section 2 describes the design and optimization process which leads to a desired candidate structure. It discusses the aspects of the hyperelastic material model as well as the use of mode shapes to construct the reduced order model. Section 3 describes the development of the structure’s Lagrangian, the extraction of the nonlinear equations of motion for the modal amplitudes, and the steady state dynamic response of the system under harmonic excitation of the higher frequency mode. Section 4 contains some concluding remarks for this work.
2. Candidate structure synthesis
The principal objective of the structural synthesis proves is to obtain a resonator design with commensurable natural frequencies. As it is difficult to come up with such a structure by just relying on the researcher’s experience, a computational optimization method is proposed to design the candidate resonators. For 1:2 internal resonances, the desired frequency relation between the two modes taking part in the energy transfer can be expressed as:
Thus, the optimization process attempts to minimize the deviation of the two natural frequencies from the perfect 1:2 natural frequency ratio. Solving the optimization problem posed by Eq. (2) would lead to a structure with two of its natural frequencies close to the ratio of 1:2 which is a major requirement for resonators exhibiting 1:2 internal resonance. In this study, two methods for solving the optimization problem are discussed. The first method is a topology optimization method based on simple isotropic material with penalization (SIMP) model and the method of moving asymptotes (MMA) . The second method is a parametric optimization method in which a starting parameterized base structure is chosen whose topology is similar to the final desired candidate structure. Then this base structure is optimized by a nonlinear quadratic programming process to produce a viable candidate structure.
2.1 Topology optimization with SIMP
Topology optimization techniques have been widely used to solve a broad range of structural optimization problems. While quite versatile, an occasional drawback against topology optimization generated optimal topologies has been the difficulty of their reproduction using conventional manufacturing processes. In this regard 3D printing is eminently suited to produce topologically optimized design as both techniques are adept at producing extruded structures with complex planar patterns. Topology optimization methods are based on finite element discretization of the design spaces. In the context of designing candidate hyperelastic resonators for 1:2 internal resonance, the design space can be discretized with finite elements and the density and material stiffness of the
where, is the material density, is the material stiffness and is the design variable which varies between 0 and 1. and are infinitesimal constants to prevent numerical singularities in case becomes equal to zero. The exponents,
As an example of the topology optimization based resonator generation process, consider the structure shown in Figure 1. This structure is a rectangular plate which is constrained at its bottom edge. This plate is assigned Mooney-Rivlin material properties and is meshed with four node planar elements as it is assumed that the plate is undergoing vibrations in its plane.
The ratio of the first two planar natural frequencies of the base structure shown in Figure 1 was 3.3. The aim of the topology optimization process is to fill the central cavity of the starting structure so that its first two planar natural frequencies are in the ratio close to 1:2. Thus, the design space is discretized with finite elements and the optimization problem posed by Eq. (2) is solved using the method moving asymptotes (MMA) which yields the optimized structure shown in Figure 2.
The ratio of the first two planar natural frequencies of the optimized structure shown in Figure 2 was 1.99. Thus, the topology optimization process was successful in bringing the natural frequencies of interest close to the ratio of 1:2. As the optimization process uses finite elements to compute the natural frequencies of the structure, the mode shapes of the optimized structure also become available and are shown in Figure 3.
2.2 Parametric optimization
Parametric optimization process is a simple but powerful tool which can also be used to generate various candidate structures for 1:2 internal resonance. As an example to illustrate the essential aspects of this procedure, consider the base structure shown in Figure 4. This base structure consists of a rectangular cantilever plate with two cutouts.
This plate can be assigned Mooney-Rivlin material properties and meshed with four node shell elements. In this study, Abaqus software is used to compute the natural frequencies of the base structure with the frequencies of interest being the second and third natural frequencies respectively. For the base structure shown in Figure 4, the ratio between the natural frequencies of the higher and lower mode of interest (third and second natural frequencies, respectively) was computed as 2.4. This base structure was then subjected to an optimization process with the objective function being described by Eq. (2). The design parameters for this optimization were the cut-out size and positions on the cantilever plate and the optimization was performed by a sequential quadratic programming method. The optimized structure obtained by applying the optimization process on the base structure is shown in Figure 5.
The ratio between the third and second natural frequencies of the optimized structure shown in Figure 5 was 2.0. Thus, the optimization method was able to successfully bring the natural frequencies of the structure close to the desired ratio of 1:2. The mode shapes of the optimized structure also become available from the finite element model and are shown in Figure 6. These mode shapes can then be used to construct a reduced order model for the system which will be used to develop the nonlinear dynamic response for the candidate structure.
The method of parametric optimization allows for development a wide range of topologies which can each potentially exhibit 1:2 internal response. The optimal topology obtained depends on the starting structure, reflecting the local optimal nature of the solution. For example, consider the starting structure shown in Figure 7, the ratio between the natural frequencies of the higher and lower mode of interest (third and second natural frequencies, respectively) was computed as 1.6. Also note that the boundary conditions in this case involve fixing the resonator along both of its vertical sides. In this particular case, the optimization parameters were the size and location of the three circular cutouts. After performing the shape optimization process again using the sequential quadratic programming method, the optimized structure obtained is shown in Figure 8. The ratio between the third and second natural frequencies of the optimized structure shown in Figure 8 was 2.0. Thus, the examples of Figures 5 and 8 illustrate the possibilities of generating a large number of examples with different topologies as candidate resonators for 1:2 internal resonance. The mode shapes of the optimized structure shown in Figure 8 are shown in Figure 9.
3. Nonlinear dynamic response
For the development of the nonlinear dynamic response of a 3D printed structure, consider the structure shown in Figure 10. This structure was designed using the simple iterative optimization procedure detailed in Section 2, and the resulting structure’s modes 2 and 3 are in near internal resonance of 1:2. Thus, the frequency ratio achieved was 2:0005. The resonator was then fabricated using 3D printing machine Stratsys Dimension 1200es. This structure has the ratio of its second and third natural frequencies as ∼2.0. The two transverse modes of interest for this candidate structure are shown in Figure 11. Using the mode shapes shown in Figure 11, assuming that the structure is subjected to a base excitation, and using the Kirchhoff plate theory , the displacement at any point on this rectangular plate can we written as:
where the dot (.) represents the derivative of the displacements with respect to time and
where, and are the first and second deviatoric invariants, respectively, of the Left Cauchy Green deformation tensor
Then the deformation gradient
The left Cauchy Green deformation tensor is given by
The deviatoric strain invariants of the left Cauchy Green deformation tensor
This Lagrangian from Eq. (19) will be a nonlinear function of the modal amplitudes owing to the nonlinear nature of the strain energy potential given in Eq. (9). The base excitation of the structure is now assumed to be of the form
Similarly, another mistuning parameter, the internal mistuning
To further study the nonlinear dynamic response of the structure for small nonlinear motions, and to formulate the application of the method of averaging , the modal amplitudes (for Eqs. (5)–(7)) can be written as
where a prime () denotes a derivative with respect to the slow time
Subjecting the averaged Lagrangian shown in Eq. (27) to the Euler-Lagrange conditions (, ,
where a prime () denotes a derivative with respect to the slow time
The steady-state solutions for the system of equations described by Eqs. (34)–(37) can be obtained by setting and . These equations can be solved for steady-state solutions to give single-mode (only second modal amplitude
As is clear from Eqs. (34)–(37), the modal amplitudes depend upon many parameters. Some of the more interesting of these are the internal mistuning
Figure 14 shows the effect of damping coefficients on the nonlinear response curves obtained using Eqs. (34)–(37). Increasing damping coefficients and has interesting effects on the overall nonlinear response. As can be observed from Figure 14a, increasing damping leads to a reduction in the frequency range in which the nonlinear coupled-mode response is observed.
4. Summary and conclusions
This work explored the possibility of synthesizing 3D printed hyperelastic plate structure exhibiting 1:2 internal resonances. 3D printing occupies a potential sweet spot in terms of dimensional capabilities and repeatability to produce nonlinear resonators which can be used as vibration absorbers, sensors, or for signal processing applications. The synthesis methodology allows for designing a large set of designs meeting the desired internal resonance conditions resulting in complex modal coupling and energy transfer between modes of the structure.
While the nonlinear dynamical response studied here was focused on 3D printed cantilever plates that exhibited 1:2 internal resonances on account of material nonlinearities, the methodology can be easily applied to other boundary conditions and internal resonances, as well as for structures with geometric nonlinearities caused by finite deformations of plates.
Conflict of interest
The authors declare no conflict of interest.