<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">WJM</journal-id><journal-title-group><journal-title>World Journal of Mechanics</journal-title></journal-title-group><issn pub-type="epub">2160-049X</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/wjm.2013.36027</article-id><article-id pub-id-type="publisher-id">WJM-36191</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Engineering</subject><subject> Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Linear, Cubic and Quintic Coordinate-Dependent Forces and Kinematic Characteristics of a Spring-Mass System
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>aiduke</surname><given-names>Sarafian</given-names></name><xref ref-type="aff" rid="aff1"><sub>1</sub></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>University College, Pennsylvania State University, York, USA</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>has2@psu.edu</email></corresp></author-notes><pub-date pub-type="epub"><day>23</day><month>08</month><year>2013</year></pub-date><volume>03</volume><issue>06</issue><fpage>265</fpage><lpage>269</lpage><history><date date-type="received"><day>June</day>	<month>4,</month>	<year>2013</year></date><date date-type="rev-recd"><day>July</day>	<month>6,</month>	<year>2013</year>	</date><date date-type="accepted"><day>August</day>	<month>1,</month>	<year>2013</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
   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. 
 
</p></abstract><kwd-group><kwd>Linear; Cubic and Quintic Nonlinear Oscillations; Semi-Analytic Solution to Equation of Motion; Mathematica</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Motivation and Objectives</title><p>Literature articles and text books are flooded with sections describing the characteristics of linear oscillators [<xref ref-type="bibr" rid="scirp.36191-ref2">2</xref>]. 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 [<xref ref-type="bibr" rid="scirp.36191-ref5">5</xref>]. 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.</p><p>To address the first issue we consider combining two linear springs; the device is shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>. 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.</p><p>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-</p><p>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.</p></sec><sec id="s2"><title>2. Physics of the Problem and Its Formulation</title><p>As shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>, 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<img src="1-4900219\690e5df0-69af-4dfb-95d6-08a63b51a3cc.jpg" />, yielding the net force</p><disp-formula id="scirp.36191-formula7326"><label>(1)</label><graphic position="anchor" xlink:href="1-4900219\b9ac529c-2d13-4e75-83a5-96f34e479b8b.jpg"  xlink:type="simple"/></disp-formula><p>For displacement x &lt; L the quantity in the parentheses in Equation (1) may be replaced with</p><p><img src="1-4900219\ecfcfef8-7a62-4cc8-a86f-5ab9830624e4.jpg" /></p><disp-formula id="scirp.36191-formula7327"><label>(2)</label><graphic position="anchor" xlink:href="1-4900219\5840b294-df29-498c-9db6-76034be972e1.jpg"  xlink:type="simple"/></disp-formula><p>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 <xref ref-type="fig" rid="fig2">Figure 2</xref>.</p><p>For this three-spring system the net force is,</p><disp-formula id="scirp.36191-formula7328"><label>(3)</label><graphic position="anchor" xlink:href="1-4900219\43a01583-acbb-47ed-9d87-c323264003de.jpg"  xlink:type="simple"/></disp-formula><p>The criteria x &lt; 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,</p><disp-formula id="scirp.36191-formula7329"><label>(4)</label><graphic position="anchor" xlink:href="1-4900219\32d5a4d2-f353-4a95-8fe2-12a6d06e4147.jpg"  xlink:type="simple"/></disp-formula><p>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 <xref ref-type="fig" rid="fig1">Figure 1</xref>, 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.</p></sec><sec id="s3"><title>3. Analysis</title><sec id="s3_1"><title>3.1. Case 1</title><p>According to the scenario shown in <xref ref-type="fig" rid="fig1">Figure 1</xref> and its ac-</p><p>companied force, Equation (2), the equation of motion is,</p><disp-formula id="scirp.36191-formula7330"><label>. (5)</label><graphic position="anchor" xlink:href="1-4900219\3add406e-aa69-4543-adca-032fc93b46d3.jpg"  xlink:type="simple"/></disp-formula><p>First we consider a case assuming x &lt; 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 [<xref ref-type="bibr" rid="scirp.36191-ref6">6</xref>]. 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,</p><p><img src="1-4900219\01d8da97-3938-4ae7-96b2-3b337f185ad9.jpg" />.</p><p>The 4 &#215; 2 graphic-matrix shown in <xref ref-type="fig" rid="fig3">Figure 3</xref> 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 <img src="1-4900219\78a8d8d5-a9c5-40af-aed5-99471abcb399.jpg" /> 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. <xref ref-type="fig" rid="fig4">Figure 4</xref> puts the observation in perspective. This graph vividly shows a reciprocal relationship between the amplitude and the corresponding period.</p><p>It is instructive to compare the oscillations emanating by a cubic force vs. the linear one. This can be done ei-</p><p>ther by a direct comparison of the amplitudes or more elegantly by comparing their respective phase diagrams shown in <xref ref-type="fig" rid="fig5">Figure 5</xref>.</p><p>For the sake of clarity the horizontal scale of <xref ref-type="fig" rid="fig5">Figure 5</xref> is magnified by a factor of 10. <xref ref-type="fig" rid="fig5">Figure 5</xref> 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 <xref ref-type="fig" rid="fig6">Figure 6</xref>.</p><p><xref ref-type="fig" rid="fig6">Figure 6</xref> 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 <img src="1-4900219\782fc3b4-b408-4baa-bc22-b15346f6ebf1.jpg" /> whereas the outer one is associated with the initial displacement<img src="1-4900219\31375dbe-b694-4465-86bb-aebf55ce2018.jpg" />. The curves between the inner and the outer ellipses correspond to initial displacements <img src="1-4900219\2a653f9f-936f-470d-b571-79eb72f9d9d2.jpg" /> α for α &#160;= 2, 3, 4, ..., 9, respectively.</p></sec></sec></body><back><ref-list><title>References</title><ref id="scirp.36191-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Mathematica, “A General Computer Software System and Language Intended for Mathematical and Other Applications,” V9.0, Wolfram Research, 2013.</mixed-citation></ref><ref id="scirp.36191-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Tronton and Marion, “Classical Dynamics of Particles and Systems,” 5th Edition, Cengage Learning, 2008.</mixed-citation></ref><ref id="scirp.36191-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">H. Sarafian, “Static Electric-Spring and Nonlinear Oscillations,” Journal of Electromagnetic Analysis &amp; Applications, Vol. 2, No. 2, 2010, pp. 75-81.</mixed-citation></ref><ref id="scirp.36191-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">H. Sarafian, “Nonlinear Oscillations of a Magneto Static Spring-Mass,” Journal of Electromagnetic Analysis and Applications, Vol. 3, No. 5, 2011, pp. 133-139. 
doi:10.4236/jemaa.2011.35022</mixed-citation></ref><ref id="scirp.36191-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">H. Sarafian, “Dynamic Dipole-Dipole Magnetic Interaction and Damped Nonlinear Oscillations,” Journal of Electromagnetic Analysis &amp; Applications, Vol. 1, No. 4, 2009, pp. 195-204. doi:10.4236/jemaa.2009.14030</mixed-citation></ref><ref id="scirp.36191-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Jacobi Elliptic Functions.  
http://mathworld.wolfram.com/JacobiEllipticFunctions.html</mixed-citation></ref><ref id="scirp.36191-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">J. Richter-Gebert and U. H. Kortenkamp, “The Interactive Geometry Software,” Cinderella 2.0, Springer, 2012.</mixed-citation></ref></ref-list></back></article>