.. _xfem-angled-crack-cylinder:

*********************
Angled Crack Cylinder
*********************

Product: Sierra/SolidMechanics - Explicit Analysis

.. figure:: angled_crack_cylinder/figs/cylinderProblemSetup.png
   :align: center
   :width: 600px
	   
   Angled crack cylinder problem set-up.

* :download:`Mesh <angled_crack_cylinder/cylinder_shell.g>`   
* :download:`Input Deck <angled_crack_cylinder/cylinder_shell.i>`

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

.. .. _fig-CylShellProblemSetup:

.. .. figure:: ../../../../tests/adagio_rtest/ExampleProblemsManual/XfemExampleProblems/AngledCrackCylinder/figs/cylinderProblemSetup.png
..    :align: center
..    :scale: 75 %

..    Angled crack cylinder problem setup.

.. The purpose of the following numerical example is to display the X-FEM cutting by prescribed object and planar crack growth capabilities for 2-D shell elements. Consider a hollow cylinder under uniform tensile loading at the rim of the cylinder and an angled prescribed crack as shown in :numref:`fig-CylShellProblemSetup`.  As the force increases as the top and bottom of the cylinder are pulled apart, a stress concentration forms at the crack tip. Once the stress concentration reaches the crack growth stress, the crack grows across that particular element. This problem is run using explicit dynamics.  

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

.. A prescribed displacement is applied to the top and bottom of the rim of the cylinder at a linear rate.  The prescribed object used for cutting is a disk whose midpoint is placed at the end of the cylinder, halfway between the center of the cylinder and the bottom rim.  The crack growth parameter that was used was chosen for this problem is  maximum principal stress.  

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

.. An elastic-plastic material model is used, and the material properties can be found in :numref:`tab-materials` 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.

.. .. _tab-materials:

.. .. csv-table:: Cylinder Material
..    :align: center
..    :delim: &

..    Young's Modulus: & :math:`E` & :math:`1.0\times 10^9`
..    Poisson's Ratio & :math:`\nu` & :math:`0.25`
..    Density & :math:`\rho` & :math:`2.61\times 10^-4`
..    Yield Stress & :math:`\sigma_{yield}` & :math:`36,000`
..    Hardening Modulus & :math:`H` & :math:`0.0`
..    Beta & :math:`\beta` & :math:`1.0`

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

.. The elements that were used for this simulation were Belytschko-Tsay shell elements.

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

.. X-FEM cut by prescribed object, and planar crack growth.

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

.. As can be seen in :numref:`fig-PCG`\ (a) below, the prescribed crack grows when the user-specified maximum principal stress is reached. The crack propagates up the geometry in a planar fashion, and a sliver of the cylinder is cut off. As the displacement is applied, the geometry separates into two separate pieces. In :numref:`fig-PCG`\ (b), the triangles over the geometry are not elements but visualization surfaces; the element used is still the Belytschko-Tsay four node shell.See :ref:`Appendix A <input-angled-crack-cylinder>` for the input deck.

.. .. _fig-PCG:

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

..    .. image:: ../../../../tests/adagio_rtest/ExampleProblemsManual/XfemExampleProblems/AngledCrackCylinder/figs/cylinderCrackGrow.png
..       :alt: Planar crack growth.
..       :scale: 40 %

..    .. image:: ../../../../tests/adagio_rtest/ExampleProblemsManual/XfemExampleProblems/AngledCrackCylinder/figs/cylinderSeparation.png
..       :alt: Piece of cylinder is cut off and separates.
..       :scale: 40 %
