prandtl_reuss.h

Go to the documentation of this file.

/*
 * Copyright (c) 2014-2015:  G-CSC, Goethe University Frankfurt
 * Author: Raphael Prohl
 * 
 * This file is part of UG4.
 * 
 * UG4 is free software: you can redistribute it and/or modify it under the
 * terms of the GNU Lesser General Public License version 3 (as published by the
 * Free Software Foundation) with the following additional attribution
 * requirements (according to LGPL/GPL v3 §7):
 * 
 * (1) The following notice must be displayed in the Appropriate Legal Notices
 * of covered and combined works: "Based on UG4 (www.ug4.org/license)".
 * 
 * (2) The following notice must be displayed at a prominent place in the
 * terminal output of covered works: "Based on UG4 (www.ug4.org/license)".
 * 
 * (3) The following bibliography is recommended for citation and must be
 * preserved in all covered files:
 * "Reiter, S., Vogel, A., Heppner, I., Rupp, M., and Wittum, G. A massively
 *   parallel geometric multigrid solver on hierarchically distributed grids.
 *   Computing and visualization in science 16, 4 (2013), 151-164"
 * "Vogel, A., Reiter, S., Rupp, M., Nägel, A., and Wittum, G. UG4 -- a novel
 *   flexible software system for simulating pde based models on high performance
 *   computers. Computing and visualization in science 16, 4 (2013), 165-179"
 * 
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
 * GNU Lesser General Public License for more details.
 */
 
#ifndef PRANDTL_REUSS_H_
#define PRANDTL_REUSS_H_
 
namespace ug{
namespace SmallStrainMechanics{
 
 
template <typename TDomain>
class PrandtlReuss
  : public IMaterialLaw<TDomain>
{
  private:
    typedef IMaterialLaw<TDomain> base_type;
 
    typedef PrandtlReuss<TDomain> this_type;
 
  public:
    static const int dim = base_type::dim;
 
    typedef typename base_type::TBaseElem TBaseElem;
 
  public:
    PrandtlReuss();
 
    ~PrandtlReuss(){};
 
    void set_bulk_modulus(const number bulkModulus)
    {matConsts.kappa = bulkModulus;
    std::stringstream ss; ss << m_materialConfiguration << " bulk modulus kappa = " << bulkModulus << "\n";
    m_materialConfiguration = ss.str();}
 
    void set_shear_modulus(const number shearModulus)
    {matConsts.mu = shearModulus;
    std::stringstream ss; ss << m_materialConfiguration << " shear modulus mu = " << shearModulus << "\n";
    m_materialConfiguration = ss.str();}
 
    void set_initial_flow_stress(const number initialFlowStress)
    {matConsts.K_0 = initialFlowStress;
    std::stringstream ss; ss << m_materialConfiguration << " initial flow-stress K_0 = " << initialFlowStress << "\n";
    m_materialConfiguration = ss.str();}
 
    void set_residual_flow_stress(const number resFlowStress)
    {matConsts.K_inf = resFlowStress;
    std::stringstream ss; ss << m_materialConfiguration << " residual flow-stress (saturation stress) K_inf = " << resFlowStress << "\n";
    m_materialConfiguration = ss.str();}
 
    void set_hardening_modulus(const number hardModulus)
    {matConsts.Hard = hardModulus; m_bHardModulus = true;
    std::stringstream ss; ss << m_materialConfiguration << " hardening modulus = " << hardModulus << "\n";
    m_materialConfiguration = ss.str();}
 
    void set_hardening_exponent(const number hardExponent)
    {matConsts.omega = hardExponent; m_bHardExp = true;
    std::stringstream ss; ss << m_materialConfiguration << " hardening exponent = " << hardExponent << "\n";
    m_materialConfiguration = ss.str();}
 
    void set_hardening_behavior(int hard);
 
    void set_tangent_precision(const number tanAccur)
    {m_tangentAccur = tanAccur;
    std::stringstream ss; ss << m_materialConfiguration << " accuracy of the tangent approximation = " << tanAccur << "\n";
    m_materialConfiguration = ss.str();}
 
 
  public:
  // INTERFACE-METHODS
    void init();
 
    void stressTensor(MathMatrix<dim,dim>& stressTens, const size_t ip,
        const MathMatrix<dim, dim>& GradU);
 
    SmartPtr<MathTensor4<TDomain::dim,TDomain::dim,TDomain::dim,TDomain::dim> >
      elasticityTensor(const size_t ip, const MathMatrix<dim, dim>& GradU);
 
    virtual bool needs_to_add_jac_m(){return false;}
 
 
    virtual void init_internal_vars(TBaseElem* elem, const size_t numIP);
    virtual void internal_vars(TBaseElem* elem);
    virtual void update_internal_vars(const size_t ip, const MathMatrix<dim, dim>& GradU);
 
    virtual void write_data_to_console(const number t);
 
 
    virtual SmartPtr<MathMatrix<dim, dim> > inelastic_strain_tensor(const size_t ip){
      MathMatrix<dim, dim>& inelastStrain = m_pElemData->internalVars[ip].strain_p_old_t;
      SmartPtr<MathMatrix<dim, dim> > spInelasticStrain (new MathMatrix<dim, dim>(inelastStrain));
      return spInelasticStrain;
    }
    virtual number hardening_parameter(const size_t ip){
      return m_pElemData->internalVars[ip].alpha;
    }
    virtual number plastic_multiplier(const size_t ip, const MathMatrix<dim, dim>& GradU);
 
 
    virtual void attach_internal_vars(typename TDomain::grid_type& grid)
    {
  //todo: Move this to an initialization routine. Best would probably to make
  //    set_domain(...) of the base class virtual and to attach m_aElemData there.
 
      grid.template attach_to<TBaseElem>(m_aElemData);
      m_aaElemData.access(grid, m_aElemData);
    }
 
    virtual void clear_attachments(typename TDomain::grid_type& grid)
    {
      if(grid.template has_attachment<TBaseElem>(m_aElemData)){
        grid.template detach_from<TBaseElem>(m_aElemData);
        m_aaElemData.invalidate();
      }
    }
 
  private:
    void strainTensor(MathMatrix<dim,dim>& strainTens, const MathMatrix<dim, dim>& GradU);
 
    void Flowrule(MathMatrix<dim, dim>& strain_p_new,
        MathMatrix<dim, dim>& strain,
        number& gamma,
        MathMatrix<dim, dim>& strial,
        MathMatrix<dim, dim>& normal,
        const MathMatrix<dim, dim>& GradU,
        const MathMatrix<dim, dim>& strain_p_old_t,
        const number alpha);
 
    void ConstLaw(MathMatrix<dim, dim>& stressTens, const MathMatrix<dim, dim>& strainTens,
        const MathMatrix<dim, dim>& strial, const number& gamma,
        const MathMatrix<dim, dim>& normal);
 
    void StressTensor(MathMatrix<dim,dim>& stressTens, const MathMatrix<dim, dim>& GradU,
            const MathMatrix<dim, dim>& strain_p_old_t, const number alpha);
 
    number Hardening(const number alpha);
 
    number Hardening_d(const number alpha);
 
    number PerfectPlasticity(const number flowcondtrial);
 
    number LinearHardening(const number flowcondtrial);
 
    number ExponentialHardening(const number strialnorm, const number alpha);
 
    void Update_internal_vars(MathMatrix<dim, dim>& strain_p_new,
        number& alpha,
        const MathMatrix<dim, dim>& GradU,
        const MathMatrix<dim, dim>& strain_p_old_t);
 
  private:
    struct InternalVars {
      //  plastic part of strain tensor wrt to last time step
      MathMatrix<dim, dim> strain_p_old_t;
 
      //  hardening variable
      number alpha;
    };
 
  //  std-vector of InternalVars
    struct ElemData{
      std::vector<InternalVars> internalVars;
    };
 
    ElemData* m_pElemData;
 
  //  attachment type: attachment of ElemData
    typedef Attachment<ElemData> AElemData;
    AElemData m_aElemData;
  //  the attachment accessor
    typedef Grid::AttachmentAccessor<TBaseElem, AElemData>  ElemDataAccessor;
    ElemDataAccessor m_aaElemData;
 
    int m_hardening;
    size_t m_MaxHardIter;
    number m_HardAccuracy;
 
    struct MaterialConstants{
      number mu;    //  shear modulus
      number kappa; //  bulk modulus
      number K_0;   //  initial flow-stress:
      number K_inf;   //  residual flow-stress (saturation stress)
      number Hard;  //  linear hardening modulus
      number omega; //  hardening exponent
    }matConsts;
 
    number m_tangentAccur;
 
    bool m_bHardModulus;
    bool m_bHardExp;
 
    number m_max_k_tan;
    number m_min_k_tan;
    size_t m_plasticIPs;
 
  public:
    using base_type::m_materialConfiguration;
};
 
}// end of namespace SmallStrainMechanics
}// end of namespace ug
 
#include "prandtl_reuss_impl.h"
 
#endif /* PRANDTL_REUSS_H_ */