A Simple Frequency Formulation for the Tangent Oscillator

: The frequency of a nonlinear vibration system is nonlinearly related to its amplitude, and this relationship is critical in the design of a packaging system and a microelectromechanical system (MEMS). This paper proposes a straightforward frequency prediction method for nonlinear oscillators with arbitrary initial conditions. The tangent oscillator, the hyperbolic tangent oscillator, a singular oscillator, and a MEMS oscillator are chosen to elucidate the simple solving process. The results, when compared with those obtained by the homotopy perturbation method, exhibit a good agreement. This paper introduces a very convenient procedure for attaining quick and accurate insight into the vibration property of a nonlinear vibration system. Contributions: Conceptualization, J.-H.H. and C.-H.H.; methodology, Q.Y. and C.-H.H.; software, Q.Y. and Y.K.; validation, and analysis, J.-H.H.; investigation,


Introduction
Vibration absorption and vibration attenuation are two critical factors in designing a nonlinear vibration system; for example, a low amplitude is always considered in the design of the packaging system [1][2][3][4] and the seismic design of architecture [5]. Active vibration control [6] has received a lot of attention in the industrial and academic communities, and many mathematicians are working to predict the periodic property of a practical vibration system. Generally, a nonlinear vibration equation is written as where w is the displacement, m is the mass, h is the nonlinear restoring force, and a and b are constants. For a linear vibration system, one can choose h(w) = kw, k is the spring coefficient, this is the well-known harmonic oscillation. When h(w) = k tan w, we have the well-known tangent oscillator arising in packaging systems [1][2][3][4]. The amplitude is determined by the initial conditions and the frequency of the system. There are numerous analytical methods for solving Equation (1), such as the Li-He method or its modifications [4,7,8], He's variational approach [9], He's Energy Balance Method [10], the iteration perturbation method [11], the exact solution method [12], the homotopy perturbation method (HMP) [13,14], the Gamma function method [15], and the variational method [16][17][18]. The goal of this paper is to suggest a simple method [19,20] for gaining a timely and efficient glimpse into the frequency-amplitude relationship of Equation (1) with arbitrary initial conditions. The comparison with the aforementioned existing methods shows that this work would be greatly challenging for nonlinear vibration theory.

The Frequency Formulation
In Equation (1), h(w)/w > 0 is required for a periodic solution. Equation (1) can be rewritten as follows: In terms of frequency formulation [19,20] is where A is the amplitude, it can be approximated calculated as In Equation (3) N was recommended as √ 3/2 in [19] for non-singular oscillator and 0.8 in [20] for singular oscillators. The frequency can be calculated using Equations (3) and (4). In the literature [21][22][23], Equation (3) is referred to as He's frequency formulation, and there have been numerous modifications [24][25][26][27]; the original formulation began with an ancient Chinese algorithm [28], which led to He's frequency formulation [19,20] for nonlinear vibration systems and Chun-Hui He's algorithm [28,29] in numerical methods.

Tangent Oscillator
The tangent oscillator [3] is One can choose N = √ 3/2 in Equation (3). We have where A is given in Equation (4). For given parameters of ω 0 , a, and b, the frequency can be calculated easily from Equation (6). To obtain an explicit formulation, one can consider the case of A << 1. In that case, we have Equation (6) becomes In view of Equation (4), we can convert Equation (8) to the form Solving ω from Equation (10) and ignoring the meaningless root gives The approximate solution reads where ϕ is determined by the initial conditions In order to reveal the accuracy of the approximate solution, we resolve the problem by the homotopy perturbation method [19]. For small amplitude, Equation (5) is approximated as The homotopy equation is where h is the homotopy parameter, 0 ≤ h ≤ 1. When h = 1, Equation (15) becomes Equation (14). We assume that the solution is written as Following the method's standard steps, we have Equation (17) is a linear differential equation; its solution is where A and ϕ are given, respectively, in Equations (4) and (13). Now Equation (18) is updated as or The solution of Equation (21) is where C is an integral constant. In Equation (22), the term of t sin(ωt + ϕ) is not periodic when t tends to infinity, so its coefficient has to be zero to have a periodic solution, that is or This result can also be obtained by the multiple scales method [24], and it is exactly the same as that given in Equation (9). Considering the initial conditions, w 1 (0) = 0, w 0 (0) = 0, Equation (22) becomes The first-order approximate solution is obtained by setting h = 1 in Equation (16), that is Comparison of the approximate solution, w = A cos(ωt + ϕ), with ω given in Equation (11), with the exact solution for various cases is illustrated in Figure 1, where ω 0 = 1.
For small amplitude, tanh w can be approximated as Equation (28) becomes: In view of Equation (4), Equation (30) turns out to be Solving ω from Equation (32) and ignoring the meaningless root gives This is the same as that obtained by the homotopy perturbation method [1]. Comparison of the approximate solution with the exact solution for various cases is illustrated in Figure 2, where ω 0 = 1.

Singular Oscillator
Now we consider a singular oscillator [9] w + 1 For the singular oscillator, we choose N = 0.8 [20], that is Simplifying Equation (35) gives Solving ω from Equation (36) leads to the result When b = 0 and N = 0.8, we have The exact frequency for b = 0 is [9] ω exact = 1.2533k −1/2 a −1 The relative error is 0.26%. Figure 3 illustrates the accuracy of the approximate solution.

MEMS Oscillator
The fast development of nanotechnology and material science have led to skyrocketing interest in MEMS systems for the last decade [30][31][32]. We consider the following MEMS oscillator [33] where y is the dimensionless displacement, y < 1, and b is constant. In order to use the above frequency formulation, we introduce a transformation: where A is the amplitude, Equation (40) becomes where f (x) is the restoring force, In view of Equation (44), we can rewrite Equation (43) in the form This leads to the following formulation Solving ω from Equation (48), we obtain Finally, we obtain the following approximate periodic solution: By the inverse transformation of Equation (41), we obtain where ω is given in Equation (49). Figure 4 shows a high accuracy of Equation (51)

Conclusions
He's frequency formulation could be extended to fractal oscillators [34] as well as nonconservative oscillators [35][36][37]. In engineering, a simple calculation is always appreciated; the simpler the calculation, the better. This paper proposes possibly the simplest method for quickly inspecting the frequency property of a nonlinear oscillator; the one-step solution yields a highly accurate result, which is quite palatable.