Linear , Cubic and Quintic Coordinate-Dependent Forces and Kinematic Characteristics of a Spring-Mass System

By combining a pair of linear springs we devise a nonlinear vibrator. For a one dimensional scenario the nonlinear force is composed of a polynomial of odd powers of position-dependent variable greater than or equal three. For a chosen initial condition without compromising the generality of the problem we analyze the problem considering only the leading cubic term. We solve the equation of motion analytically leading to The Jacobi Elliptic Function. To avoid the complexity of the latter, we propose a practical, intuitive-based and easy to use alternative semi-analytic method producing the same result. We demonstrate that our method is intuitive and practical vs. the plug-in Jacobi function. According to the proposed procedure, higher order terms such as quintic and beyond easily may be included in the analysis. We also extend the application of our method considering a system of a three-linear spring. Mathematica [1] is being used throughout the investigation and proven to be an indispensable computational tool.


Motivation and Objectives
Literature articles and text books are flooded with sections describing the characteristics of linear oscillators [2].One quick review of these resources reveals these are confined limited to mechanical systems.In the area of electro and magneto-dynamics, the author thoroughly has investigated scenarios of nonlinear oscillators [3,4].Laboratory setups conducive to these scenarios have been proposed, and for the magneto-dynamic case the validity of the theoretical model is qualified with actual data [5].For the sake of completeness, therefore, it is essential to fill in the gap considering a practical scenario conducive to a nonlinear mechanical oscillator.Although the latter is the main motivation of tackling the issue, however, in the course of analyzing the problem we stumbled on a mathematical observation constituting our secondary objective.
To address the first issue we consider combining two linear springs; the device is shown in Figure 1.The springs are identical in shape and have the same stiffness.They are fastened to two separate supports at one end and jointly hold a marble at the other end.In the absence of gravity, the springs exert a nonlinear force on the marble along the horizontal direction when it is pulled away symmetrically from equilibrium.The non-linear force could be a combination of cubic and quintic coordinate-dependent terms; details are discussed in the next section.
Releasing the marble somewhere off from equilibrium results in oscillations along the horizontal direction.The corresponding equation of motion is a nonlinear differential equation.The issues concerning the solution of the equation will be discussed later.Nonetheless, it is worthwhile mentioning not only we solve the equation utiliz- ing traditional methods, we propose a fresh semi-analytic method resulting in the same output.With these objectives we craft this four-section article.In addition to Motivation and Objectives, in Section 2, along with the physics of the problem we formulate the problem.In Section 3 we present the Analysis along with the associated descriptive graphs.We close the article with conclusions and closing remarks.

Physics of the Problem and Its Formulation
As shown in Figure 1, the marble is pulled from its initial equilibrium position to a point such as G, stretching the linear springs with stiffness k beyond their relaxed length L. The differential length elongation of each spring is For displacement x < L the quantity in the parentheses in Equation ( 1) may be replaced with One realizes this simple mechanical device is capable of exerting nonlinear forces; the forces are a combination of cubic and quintic coordinate-dependent terms.For the sake of "completeness", in order to include a linear coordinate-dependent force we modify the previous design by adding a third identical spring, shown in Figure 2.
For this three-spring system the net force is, The criteria x < L suggests that the impact of the third term and terms with powers higher than the fifth are to be insignificant; this is qualified later on.Nonetheless, the equation of motion of the marble under the influence of the force given by Equation (3) along the x-axis yields, with this equation in hand we analyze the impact of the force terms on the kinematic characteristics of the mobile marble.First, by turning off the linear term, a scenario subject to Figure 1, we focus on the impact of the cubic and quintic terms.Then, we consider a case including the combination of the linear and the cubic terms.

Case 1
According to the scenario shown in Figure 1 and its ac- companied force, Equation ( 2), the equation of motion is, First we consider a case assuming x < L, i.e. the initial displacement of the marble from equilibrium is less than the length of the spring; this drops the last term.Sustaining the cubic term only, the equation of motion is an analytically solvable DE with closed form solutions given by Jacobi Elliptic Functions [6].Alternatively, the same equation can be solved numerically.Assigning parameters to the spring length and stiffness and the marble we solve the equation numerically; output graphs of these two approaches are indistinguishable.To establish the basis for the semi-analytic method we focus on the numeric schematic.The values lists a typical set of parameters in MKS units,   2 3 3.0, 6.0 10 , 10.0 10 k m The 4 × 2 graphic-matrix shown in Figure 3 displays the position of the marble vs. time.The differences between the graphs are due only to various initial displacements of the marble.From the top left to the bottom right these correspond to respectively.The common global feature of these plots indicate irrespective of the initial condition the marble under the influence of the cubic force does oscillate.The period of oscillation is not constant; it depends on the initial value of the displacement.This is one the distinct characteristics of a nonlinear force.An inspection of these plots reveals the absence of retarding forces such as friction enforces the constant amplitude; noting the shorter the amplitude the longer the period.Figure 4 puts the observation in perspective.This graph vividly shows a reciprocal relationship between the amplitude and the corresponding period.
It is instructive to compare the oscillations emanating by a cubic force vs. the linear one.This can be done ei- ther by a direct comparison of the amplitudes or more elegantly by comparing their respective phase diagrams shown in Figure 5.
For the sake of clarity the horizontal scale of Figure 5 is magnified by a factor of 10. Figure 5 shows the distinct differences between the oscillation characteristics of a linear vs. a cubic spring.As one may easily prove the phase diagram of a linear spring is a perfect ellipse.The phase diagram of a cubic spring is a depressed "ellipse" along the velocity axis.Depression severity of the latter is a function of the initial condition.A family of such curves is shown in Figure 6.
Figure 6 displays a series of phase diagrams each of which is a result of applying various initial conditions.For instance the inner depressed ellipse comes about from setting the initial displacement to whereas the outer one is associated with the initial dis- The curves between the inner and the outer ellipses correspond to initial displacements   0 0. x  1 α for α = 2, 3, 4, ..., 9, respectively.

Semi-Analytic Method of Solving DE Motion
As we discussed in Section 3.1 the solution of equation of motion of a cubic force is an oscillating function.Accordingly, the value of the initial condition impacts the period.With these two observations we propose a compatible alternative solution.This solution is neither purely numeric nor purely symbolic; it is somewhere in between, it is semi-analytic.From an analytic point of view, first we consider a solution such as,   amp , from a numeric point of view utilizing either Figure 4 or 5 we hunt for T(amp).For instance, for the latter utilizing These two solutions are indistinguishable.Therefore,  as an option, rather than utilizing the numeric solution that in general lacks the physics insight, alternatively, one may objectively apply the semi-analytic method to obtaining the same result.

Case 2
In this section we present the results pertaining scenario shown in Figure 2. The equation of motion is, Equation ( 6) has an analytic closed solution; these are Jacobi Elliptic Functions [6,7].These functions compar-   posed method.
As in Case 1, solutions are indistinguishable.

Conclusions and Remarks
The objective of this study is to design a mechanical system capable of oscillating under the influence of nonlinear forces.Utilizing a set of two and three linear springs we design two such devices; many more may be devised.These are complementary to our previous studies completing the scope of nonlinear oscillations in three different areas of physics: mechanics, electrodynamics and

Figure 1 .
Figure 1.Two-spring arrangement leading to cubic and quintic oscillations.

Figure 2 .
Figure 2. Three-spring arrangement leading to a combination of linear, cubic and quintic oscillations.

Figure 3 .
Figure 3. Output of numeric solution of equation of motion.From the top left to the bottom right the graphs correspond to the initial values x[0] = 0.1 α l for α = 1, 2, 3, ..., 8, respectively.

Figure 4 .
Figure 4. Collective display of the oscillations shown in Figure 3.

Figure 5 .
Figure 5. Phase diagrams: (a) the ellipse corresponds to a linear spring, and (b) the depressed ellipse corresponds to a cubic spring.

Figure 5 7 .
for the chosen initial amplitude we evaluate the maximum ordinate of the corresponding oscillation.A set of one such coordinates are shown with dots in Figure Then we fit the data with an appropriate continuous function.The fitted curve is shown on the same graph as well.Utilizing this function, e.g.T(amp) we compare the semi-analytic solution vs. the numeric solution.A 4 × 2 graphic-matrix shown in Figure 8, qualifies the accuracy of our proposed method.

Figure 6 .
Figure 6.A family of phase diagrams of the cubic spring.The inner curve corresponds to x[0] = 0.1 while the outer one is due x[0] = 1.

Figure 7 .
Figure 7.The dots are the periods and the corresponding amplitudes are deduced from Figure 5.The solid line is the fitted curve T(amp) = 0.0256907 amp −1 .

Figure 8 .
Figure 8.Comparison of the numeric solutions (black curves) vs. semi-analytic solution (gray curves) shown in Figure 2.

Figure 9 .
Figure 9.The dots are the periods and their corresponding amplitudes are deduced from oscillations subject to Equation (6).The solid line is the fitted curve T(amp) = 0.3649 − 6.8834 amp 1.5327 .