.. _linear_elasticity-linear_elastic_probes:

linear_elasticity/linear_elastic_probes.py
==========================================

**Description**


This example shows how to use the post_process_hook and probe_hook options.

Use it as follows (assumes running from the sfepy directory; on Windows, you
may need to prefix all the commands with "python " and remove "./"):

1. solve the problem:

   ./simple.py examples/linear_elasticity/linear_elastic_probes.py

2. optionally, view the results:

   ./postproc.py cylinder.h5 -b

3. optionally, convert results to VTK, and view again:

   ./extractor.py -d cylinder.h5
   ./postproc.py cylinder.vtk -b

4. probe the data:

   ./probe.py examples/linear_elasticity/linear_elastic_probes.py cylinder.h5

Find :math:`\ul{u}` such that:

.. math::
    \int_{\Omega} D_{ijkl}\ e_{ij}(\ul{v}) e_{kl}(\ul{u})
    = 0
    \;, \quad orall \ul{v} \;,

where

.. math::
    D_{ijkl} = \mu (\delta_{ik} \delta_{jl}+\delta_{il} \delta_{jk}) +
    \lambda \ \delta_{ij} \delta_{kl}
    \;.


.. image:: /../gallery/images/linear_elasticity-linear_elastic_probes.png


:download:`source code </../examples/linear_elasticity/linear_elastic_probes.py>`

.. literalinclude:: /../examples/linear_elasticity/linear_elastic_probes.py

