.. _newton_cradle:

*************
Newton Cradle
*************

Product: Sierra/SolidMechanics - Explicit Analysis

The files required for this example can be downloaded here, or found at ``$TESTS/adagio_rtest/ExampleProblemsManual/ContactExampleProblems/NewtonCradle`` (where ``TESTS=/projects/sierra/tests`` in a CEE environment, or your Sierra test distribution otherwise).



.. list-table:: 

 * - .. figure:: newton_cradle/figs/newtoncradle_final_1balldrop.png
		   
           1 ball

   - .. figure:: newton_cradle/figs/newtoncradle_final_2balldrop.png

           2 ball

   - .. figure:: newton_cradle/figs/newtoncradle_final_3balldrop.png

           3 ball

* :download:`1 Ball Drop Mesh <newton_cradle/newtoncradle_final_1balldrop.g>`  
* :download:`2 Ball Drop Mesh <newton_cradle/newtoncradle_final_2balldrop.g>`  
* :download:`3 Ball Drop Mesh <newton_cradle/newtoncradle_final_3balldrop.g>`  
* :download:`Input Deck <newton_cradle/newton_cradle.i>`

.. Problem Description
.. ===================

.. This example demonstrates the conservation of momentum and kinetic energy (through basic Newtonian mechanics) within an explicit dynamic analysis with contact. There are currently three geometric cases available for running: the dropping of one, two, and three balls. These cases can be tested through commenting/uncommenting the Genesis files of interest referenced inside the input deck. These configurations can be seen in Figure  :numref:`fig-threeballs`. The 5 ball-chain system contains a default initial configuration of one raised ball.

.. The balls are given an initial rotational displacement through their geometrical location specified in Cubit. The at rest balls touch at an initial position as fine as the mesh generated. 

.. Loading and Boundary Conditions
.. ===============================

.. The pendulum wires are defined as a truss section with each of the five balls containing an inner rigid body core. Connected in a v-shaped manner, the uppermost 'ceiling' node of the wires define the reference location to each of the five inner rigid body volumes. These rigid bodies are only allowed to rotate along the z-axis and translate along the x-axis and y-axis. Controlling the rigid body displacements in this manner allows one to bypass making the spherical nodeset rigid, hence restricting movement.

.. A constant and uniform gravitational load is applied to the system. 

.. Material Model
.. ==============

.. Assuming small deformations and simple linear elastic behavior, the outermost material for each respective sphere is defined in an elastic model. Modulus of elasticity values were picked out of convenience, i.e., stiff enough for minimal elastic deformation and compliant enough to minimize computational expense for the explicit analysis. Iterations performed to optimize these parameters conceded the outer sphere's Young's modulus at approximately :math:`10^4` times less than that of steel. The inner spheres of the model were defined as rigid bodies. The material properties can be found in :numref:`tab-newton_mats` below. Metric units are used: displacement is in meters, mass is in kilograms, time is in seconds, force is in Newtons, and temperature is in Kelvin.

.. .. %Each block uses an elastic material model where the values were picked for convenience, i.e.,
.. .. %stiff enough to make the moving block act as a rigid body and compliant enough to not drive the
.. .. %time step too small for the explicit analysis.

.. .. _tab-newton_mats:

.. .. csv-table:: Newton Cradle Materials
..    :align: center
..    :delim: &

..    Young's Modulus: Outer Sphere & E & :math:`200\times10^5`
..    Young's Modulus: Core & E & Rigid
..    Young's Modulus: Wire/String & E & :math:`100`
..    Poisson's Ratio & :math:`\nu` & 0.3
..    Density & :math:`\rho` & :math:`7.48\times10^3`

.. Finite Element Model
.. ====================

.. To eliminate rigid body contact of the outer spheres, mesh creation of this model required independent nodeset specification for inner and outer spherical volumes (see figure  :numref:`fig-RB`). The overlap of volumes was eliminated by first creating inner and outer spheres of radii 1 and 0.25 units respectively, followed by elimination/recreation of the inner sphere volumes. Both inner and outer spheres were then meshed using a standard four noded tetrahedral mesh. 

.. The wires were given a truss element type. 

.. Feature Tested
.. ==============

.. Explicit contact and rigid bodies.

.. Results and Discussion
.. ======================

.. In an idealized Newton Cradle exhibiting perfectly elastic collisions, the Total Kinetic Energy before and after a collisions would be equal. In reality, the balls are elastic, so some energy is stored in them as elastic Strain Energy of the balls themselves. Moreover, artificial bulk viscosity is used by default to damp out the high frequency responses from the system. Artificial viscosity is dissipative and so removes Total Energy from the system, hence Kinetic Energy is not conserved.

.. The Internal Energy is calculated as the sum of the nodal internal forces times the nodal velocities, integrated over time.  In contrast, the Strain Energy is calculated as the stress times the strain rate integrated over time.  The Internal Energy includes dissipated energy due to artificial viscosity and stored energy associated with hourglass modes which are not accounted for in the Strain Energy.  The difference in Internal Energy and Strain Energy normalized by the Maximum Initial Potential Energy can be seen in :numref:`fig-normdis-energyhist`\ (a). At sphere-to-sphere contact initiation, steep rates of change in Internal Energy and Strain Energy are present, followed by an approximately constant period when the opposing ball rebounds. This stair-step cycle continues with subsequent oscillations of the Newton Cradle. System energy conservation, including Kinetic, Internal, and Potential Energy, is presented in :numref:`fig-normdis-energyhist`\ (b). The increasing Internal Energy over the simulation can be attributed to bulk viscosity.

.. At 20 seconds of test run time under the configured geometries and loading, the Newton Cradle will swing for approximately 2 periods. The run time for energy dissipation of the system can be carried out as desired through input deck specification. See :ref:`Appendix A <input-newton_cradle>` for the input deck.

.. .. _tab-newton_abbrv:

.. .. csv-table:: Abbreviations Used in Results
..    :align: center
..    :delim: &

..    Kinetic Energy & KE
..    Internal Energy & IE
..    Potential Energy & PE
..    Strain Energy & SE
..    Total Initial Energy & TIE

.. .. _fig-threeballs:

.. .. subfigure:: ABC
..    :subcaptions: below
..    :align: center

..    .. image:: ../../../../tests/adagio_rtest/ExampleProblemsManual/ContactExampleProblems/NewtonCradle/figs/newtoncradle_final_1balldrop.png
..       :scale: 10 %

..    .. image:: ../../../../tests/adagio_rtest/ExampleProblemsManual/ContactExampleProblems/NewtonCradle/figs/newtoncradle_final_2balldrop.png
..       :scale: 10 %

..    .. image:: ../../../../tests/adagio_rtest/ExampleProblemsManual/ContactExampleProblems/NewtonCradle/figs/newtoncradle_final_2balldrop.png
..       :scale: 10 %

..    Initial Configurations

.. .. _fig-RB:

.. .. figure:: ../../../../tests/adagio_rtest/ExampleProblemsManual/ContactExampleProblems/NewtonCradle/figs/newtoncradle_final_rigid.png
..    :align: center
..    :scale: 20 %

..    Rigid Body Highlight

.. .. _fig-normdis-energyhist:

.. .. subfigure:: AB
..    :subcaptions: below
..    :align: center

..    .. image:: ../../../../tests/adagio_rtest/ExampleProblemsManual/ContactExampleProblems/NewtonCradle/figs/Normalized_Diss.png
..       :alt: Normalized Dissipation Energy
..       :scale: 60 %

..    .. image:: ../../../../tests/adagio_rtest/ExampleProblemsManual/ContactExampleProblems/NewtonCradle/figs/Energy_HistoryKJ.png
..       :alt: Energy History
..       :scale: 70 %


.. %\begin{figure}[h]
.. %\centering{\includegraphics[width=.95\linewidth]{figs/newtoncradle_final_1balldrop_KEVsTime}}
.. %\caption{Newton Cradle: Kinetic Energy}
.. %\label{fig-NCKE}
.. %\end{figure}
.. %
.. %
.. %\begin{figure}[h]
.. %\centering{\includegraphics[width=.95\linewidth]{figs/newtoncradle_final_1balldrop_MomVsTime}}
.. %\caption{Newton Cradle: Momentum}
.. %\label{fig-NCM}
.. %\end{figure}

.. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

.. %\begin{figure}
.. %\centering{\includegraphics[width=0.5\linewidth]{figs/NewtonCradle1}}
.. %\caption{Newton Cradle: Rigid Inner Sphere}
.. %\end{figure}


.. %\begin{figure}[h]
.. %\centering
.. %\mbox{\subfloat[KE]{\epsfig{figure=figs/newtoncradle_final_1balldrop_MomVsTime,width=.35\textwidth}}\quad
.. %      \subfloat[Momentum]{\epsfig{figure=figs/newtoncradle_final_1balldrop_KEVsTime,width=.35\textwidth}}}\quad
.. %\caption{Results}
.. %\label{fig-subs}
.. %\end{figure}


.. %\begin{minipage}[c]{\linewidth}
.. %    \centering
.. %    \begin{minipage}{0.32\linewidth}
.. %        \begin{figure}[H]
.. %            \includegraphics[width=\linewidth]{figs/newtoncradle_final_1balldrop.png}
.. %            \caption{Clamped Edges: n=3, m=1}
.. %        \end{figure}
.. %    \end{minipage}
.. %    \hspace{0.1\linewidth}
.. %    \begin{minipage}{0.32\linewidth}
.. %        \begin{figure}[H]
.. %            \includegraphics[width=\linewidth]{figs/newtoncradle_final_2balldrop.png}
.. %            \caption{Clamped Edges: n=3, m=1}
.. %        \end{figure}
.. %    \end{minipage}
.. %\end{minipage}
.. %\begin{minipage}[c]{\linewidth}
.. %    \centering
.. %    \begin{minipage}{0.32\linewidth}
.. %        \begin{figure}[H]
.. %            \includegraphics[width=\linewidth]{figs/newtoncradle_final_3balldrop.png}
.. %            \caption{Unclamped Edges: n=3, m=1}
.. %        \end{figure}
.. %    \end{minipage}
.. %    \hspace{0.1\linewidth}


.. %\begin{minipage}{1.0in}\begin{center}\includegraphics[width=0.5in]{figures/Achtung.eps}\end{center}\end{minipage}%

.. %\begin{figure}
.. %        \centering
.. %        \begin{subfloat}[b]{0.3\textwidth}
.. %                \centering
.. %                \includegraphics[width=\textwidth]{figs/newtoncradle_final_1balldrop}
.. %                \caption{Newton Cradle: 1-Ball Swing}
.. %                \label{fig-NewtonCradle1}
.. %        \end{subfloat}%
.. %       ~ %add desired spacing between images, e. g. ~, \quad, \qquad etc.
..           %(or a blank line to force the subfloat onto a new line)
.. %        \begin{subfloat}[b]{0.3\textwidth}
.. %                \centering
.. %                \includegraphics[width=\textwidth]{figs/newtoncradle_final_2balldrop}
.. %                \caption{Newton Cradle: 2-Ball Swing}
.. %                \label{fig-NewtonCradle2}
.. %        \end{subfloat}
.. %       ~ %add desired spacing between images, e. g. ~, \quad, \qquad etc.
..           %(or a blank line to force the subfloat onto a new line)
.. %        \begin{subfloat}[b]{0.3\textwidth}
.. %                \centering
.. %                \includegraphics[width=\textwidth]{figs/newtoncradle_final_2balldrop}
.. %                \caption{Newton Cradle: 3-Ball Swing}
.. %                \label{fig-NewtonCradle3}
.. %        \end{subfloat}
.. %        \caption{Initial Configurations}\label{fig-NewtonCradlePositions}
.. %\end{figure}



.. %\subsection{Reference}
