.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "examples/ex17_torsion-gif.py"
.. LINE NUMBERS ARE GIVEN BELOW.

.. only:: html

    .. note::
        :class: sphx-glr-download-link-note

        :ref:`Go to the end <sphx_glr_download_examples_ex17_torsion-gif.py>`
        to download the full example code.

.. rst-class:: sphx-glr-example-title

.. _sphx_glr_examples_ex17_torsion-gif.py:


Animated torsional-shear loading
--------------------------------

.. topic:: Save an animation of a custom boundary condition.

   * assign non-homogeneous displacements to a boundary condition

   * create an animation

   * export a GIF-movie

   * evaluate the reaction moment on a boundary condition

This example demonstrates how to create a :class:`~felupe.Boundary` for torsional
loading. This is somewhat more complicated than a simple boundary condition with
identical mesh-point values because the values of the mesh-points are different.
However, this is not a problem. FElupe supports arrays to be passed as the ``value``-
argument of a :class:`~felupe.Boundary`. This is even possible in a
:class:`~felupe.Step` with the ``ramp``-argument.

.. GENERATED FROM PYTHON SOURCE LINES 22-49

.. code-block:: Python


    import matplotlib.pyplot as plt
    import numpy as np

    import felupe as fem

    mesh = fem.mesh.Triangle(n=3).expand(n=6)
    field = fem.FieldContainer([fem.Field(fem.RegionHexahedron(mesh), dim=3)])

    boundaries = {
        "bottom": fem.Boundary(field[0], fz=0),
        "top": fem.Boundary(field[0], fz=1),
    }
    solid = fem.SolidBody(umat=fem.NeoHookeCompressible(mu=1, lmbda=2), field=field)

    angles_deg = fem.math.linsteps([0, -30, 0], num=10)
    move = []
    for phi in angles_deg:
        top = mesh.points[boundaries["top"].points]
        top_rotated = fem.math.rotate_points(
            points=top,
            angle_deg=phi,
            axis=2,
            center=[0, 0, 1],
        )
        move.append(top_rotated - top)








.. GENERATED FROM PYTHON SOURCE LINES 51-55

The reaction moment on the centerpoint of the right end face is tracked by a
``callback()`` function when we :meth:`~felupe.Job.evaluate` the :class:`~felupe.Job`.
During the callback, the animated deformed body is recorded in a GIF-file. After all
frames are recorded, it is important to ``close()`` the plotter.

.. GENERATED FROM PYTHON SOURCE LINES 55-85

.. code-block:: Python

    moment = []

    plotter = field.plot(
        "Principal Values of Logarithmic Strain", clim=[0, 0.2], off_screen=True
    )
    plotter.open_gif("result.gif", fps=5)


    def record(stepnumber, substepnumber, substep, plotter):
        "Update the mesh-points and the scalars of the plotter."
        if substepnumber in np.arange(len(move), step=2):
            name = plotter.mesh.active_scalars_info.name
            data = substep.x.evaluate.log_strain(tensor=False).mean(-2)[-1]

            plotter.mesh.points[:] = mesh.points + field[0].values
            plotter.mesh[name] = data

            plotter.write_frame()

        # evaluate the reaction moment at the centerpoint of the right end face
        forces = substep.fun
        M = fem.tools.moment(field, forces, boundaries["top"], centerpoint=[0, 0, 1])
        moment.append(M)


    step = fem.Step(items=[solid], ramp={boundaries["top"]: move}, boundaries=boundaries)
    job = fem.Job(steps=[step], callback=record, plotter=plotter).evaluate()

    plotter.close()





.. image-sg:: /examples/images/sphx_glr_ex17_torsion-gif_001.gif
   :alt: ex17 torsion gif
   :srcset: /examples/images/sphx_glr_ex17_torsion-gif_001.gif
   :class: sphx-glr-single-img







.. GENERATED FROM PYTHON SOURCE LINES 86-87

Finally, let's plot the reaction moment vs. torsion angle curve.

.. GENERATED FROM PYTHON SOURCE LINES 87-91

.. code-block:: Python

    fig, ax = plt.subplots()
    ax.plot(abs(angles_deg), abs(np.array(moment)[:, 2]), "o-")
    ax.set_xlabel(r"Torsion Angle $|\phi_3|$ in deg $\rightarrow$")
    ax.set_ylabel(r"Torsion Moment $|M_3|$ in Nmm $\rightarrow$")



.. image-sg:: /examples/images/sphx_glr_ex17_torsion-gif_002.png
   :alt: ex17 torsion gif
   :srcset: /examples/images/sphx_glr_ex17_torsion-gif_002.png
   :class: sphx-glr-single-img






.. rst-class:: sphx-glr-timing

   **Total running time of the script:** (0 minutes 2.284 seconds)


.. _sphx_glr_download_examples_ex17_torsion-gif.py:

.. only:: html

  .. container:: sphx-glr-footer sphx-glr-footer-example

    .. container:: sphx-glr-download sphx-glr-download-jupyter

      :download:`Download Jupyter notebook: ex17_torsion-gif.ipynb <ex17_torsion-gif.ipynb>`

    .. container:: sphx-glr-download sphx-glr-download-python

      :download:`Download Python source code: ex17_torsion-gif.py <ex17_torsion-gif.py>`

    .. container:: sphx-glr-download sphx-glr-download-zip

      :download:`Download zipped: ex17_torsion-gif.zip <ex17_torsion-gif.zip>`


.. only:: html

 .. rst-class:: sphx-glr-signature

    `Gallery generated by Sphinx-Gallery <https://sphinx-gallery.github.io>`_