An Object Finite Element Library

Equa_Solid.h

Go to the documentation of this file.

/*==============================================================================

                                    O  F  E  L  I

                           Object  Finite  Element  Library

  ==============================================================================

   Copyright (C) 1998 - 2008 Rachid Touzani

   This program is free software; you can redistribute it and/or modify it under
   the terms of the GNU General Public License as published by the Free
   Software Foundation; Version 2 of the License.

   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 General Public License for more
   details.

   You should have received a copy of the GNU General Public License
   along with this program; if not, write to the :

   Free Software Foundation
   Inc., 59 Temple Place - Suite 330
   Boston, MA  02111-1307, USA

  ==============================================================================

                Definition of abstract class 'Equa_Solid' to
                to solve Solid and Structural Mechanicsr problems

  ==============================================================================*/


#ifndef __EQUA_SOLID_H
#define __EQUA_SOLID_H

#include "Equation.h"

namespace OFELI {

template<class T_, size_t NEN_, size_t NEE_, size_t NSN_, size_t NSE_> class Equa_Solid;

template<class T_, size_t NEN_, size_t NEE_, size_t NSN_, size_t NSE_>
class Equa_Solid : virtual public Equation<T_,NEN_,NEE_,NSN_,NSE_>
{

 public:

   using Equation<T_,NEN_,NEE_,NSN_,NSE_>::setMaterialProperty;
   using Equation<T_,NEN_,NEE_,NSN_,NSE_>::_theElement;
   using Equation<T_,NEN_,NEE_,NSN_,NSE_>::_theSide;

#ifndef DOXYGEN_SHOULD_SKIP_THIS
enum PDE_TERMS {
  MASS                =   1,     
  LUMPED_MASS         =   2,     
  STIFFNESS           =   4,     
  DEVIATORIC          =   8,     
  DILATATION          =  12,     
  LOAD                =  20,     
  TRACTION            =  32      
};
#endif /* DOXYGEN_SHOULD_SKIP_THIS */

    Equa_Solid() { }

    virtual ~Equa_Solid() { }

    virtual void LMassToLHS(double coef=1) { _coef = coef; }

    virtual void LMassToRHS(double coef=1) { _coef = coef; }

    virtual void MassToLHS(double coef=1) { _coef = coef; }

    virtual void MassToRHS(double coef=1) { _coef = coef; }

    void setLumpedMass() 
    {
       LMassToLHS(1./AbsEqua<T_>::_time_step);
       LMassToRHS(1./AbsEqua<T_>::_time_step);
    }

    void setMass() 
    {
       MassToLHS(1./AbsEqua<T_>::_time_step);
       MassToRHS(1./AbsEqua<T_>::_time_step);
    }

    virtual void Deviator(double coef=1) { coef=1; }

    virtual void DeviatorToRHS(double coef=1) { coef=1; }

    virtual void Stiffness(double coef=1) { coef=1; }

    virtual void StiffnessToRHS(double coef=1) { coef=1; }

    void setDeviator()
    {
       if (AbsEqua<T_>::time_scheme_==AbsEqua<T_>::STEADY_STATE || 
           AbsEqua<T_>::time_scheme_==AbsEqua<T_>::FORWARD_EULER || 
           AbsEqua<T_>::time_scheme_==AbsEqua<T_>::BACKWARD_EULER)
          Deviator(1.);
       else if (AbsEqua<T_>::time_scheme_==AbsEqua<T_>::CRANK_NICOLSON)
          Deviator(0.5);
       if (AbsEqua<T_>::time_scheme_==AbsEqua<T_>::CRANK_NICOLSON)
          DeviatorToRHS(0.5);
    }

    void setStiffness()
    {
       if (AbsEqua<T_>::_time_scheme==STEADY_STATE || AbsEqua<T_>::_time_scheme==BACKWARD_EULER)
          Stiffness(1.);
       else if (AbsEqua<T_>::_time_scheme==FORWARD_EULER)
          StiffnessToRHS();
       else;
    }


 protected:

    void Young(const double &E) { _E = E; }

    void Poisson(const double &nu) { _nu = nu; }

    void Density(const double &rho) { _rho = rho; }

    void Young(const char *exp) { _E = setMaterialProperty(exp,"Young Modulus"); }

    void Poisson(const char *exp) { _nu = setMaterialProperty(exp,"Poisson Coefficient"); }

    void Density(const char *exp) { _rho = setMaterialProperty(exp,"Density"); }

    void setMaterial()
    {
       _rho = _theElement->getMaterial()->Density();
       _nu  = _theElement->getMaterial()->PoissonRatio();
       _E   = _theElement->getMaterial()->YoungModulus();
    }
#ifndef DOXYGEN_SHOULD_SKIP_THIS
   double  _E, _nu, _lambda, _G, _rho, _coef;
#endif /* DOXYGEN_SHOULD_SKIP_THIS */
};


} /* namespace OFELI */

#endif

Copyright © 1998-2007 Rachid Touzani