/home/hauberg/Dokumenter/Capture/humim-tracker-0.1/src/OpenTissue/OpenTissue/dynamics/fem/fem_initialize_plastic.h

Go to the documentation of this file.

#ifndef OPENTISSUE_DYNAMICS_FEM_FEM_INITIALIZE_PLASTIC_H
#define OPENTISSUE_DYNAMICS_FEM_FEM_INITIALIZE_PLASTIC_H
//
// OpenTissue Template Library
// - A generic toolbox for physics-based modeling and simulation.
// Copyright (C) 2008 Department of Computer Science, University of Copenhagen.
//
// OTTL is licensed under zlib: http://opensource.org/licenses/zlib-license.php
//
#include <OpenTissue/configuration.h>

#include <OpenTissue/dynamics/fem/fem_compute_b.h>
#include <OpenTissue/dynamics/fem/fem_compute_isotropic_elasticity.h>

namespace OpenTissue
{
  namespace fem
  {
    namespace detail
    {
      template < typename tetrahedron_iterator,typename real_type >
      inline void initialize_plastic(
        tetrahedron_iterator begin
        , tetrahedron_iterator end
        , real_type const & c_yield
        , real_type const & c_creep
        , real_type const & c_max
        )
      {
        assert(c_yield>=0 || !"initialize_plastic(): yield must be non-negative");
        assert(c_creep>=0 || !"initialize_plastic(): creep must be non-negative");
        assert(c_max>=0   || !"initialize_plastic(): max must be non-negative");

        for (tetrahedron_iterator T = begin; T != end; ++T)
        {
          T->m_yield = c_yield;
          T->m_creep = c_creep;
          T->m_max = c_max;
          for(int i=0;i<6;++i)
            T->m_plastic[i] = 0;

          compute_B(T->m_e10, T->m_e20, T->m_e30, T->m_V, T->m_B);
          compute_isotropic_elasticity_vector (T->m_young, T->m_poisson, T->m_D);
        }
      }

    } // namespace detail
  } // namespace fem
} // namespace OpenTissue

//OPENTISSUE_DYNAMICS_FEM_FEM_INITIALIZE_PLASTIC_H
#endif