gmpMechanicalLargeDisplacementSolid.h

Go to the documentation of this file.

/************************************************************************
**
** Copyright (C) 2014 by Carlos Augusto Teixera Mendes
** All rights reserved.
**
** This file is part of the "GeMA" software. It's use should respect
** the terms in the license agreement that can be found together 
** with this source code.
** It is provided AS IS, with NO WARRANTY OF ANY KIND, 
** INCLUDING THE WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR 
** A PARTICULAR PURPOSE.
**
************************************************************************/

#ifndef _GEMA_PLUGIN_MECHANICAL_PHYSICS_LARGEDISPLACEMENT_SOLID_H_
#define _GEMA_PLUGIN_MECHANICAL_PHYSICS_LARGEDISPLACEMENT_SOLID_H_

#include "gmpMechanicalSolid.h"
#include <gmpFemPhysicsCommon.h>
#include <gmTrace.h>


class GmGaussAccessor;

class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpMechanicalLargeDisplacementSolid: public GmpMechanicalSolid
{
public:

  // constructor
        GmpMechanicalLargeDisplacementSolid(const char* pluginType, GmSimulationData* simulation, QString id, QString description,
                 const GmpFemPhysicsCommonMaterialFactory* matFactory, const GmLogCategory& logger);

  // destructor
  virtual ~GmpMechanicalLargeDisplacementSolid();

  enum PhysicsAttributeIds
  {
          REFERENCE_LD_ID= GmpMechanicalSolid::NUM_ATTRI_IDS,        


          // ------         NO ADDING BELOW THIS LINE
          NUM_ATTRI_IDS,    
  };

  // set physic map
  virtual const QVariantMap* physicsMetaDataMap();
  virtual FemResultType fillElementData(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);
  
  virtual double fillBuMatrix(const GmElement* e, const GmShape* shape, const GmVector& ncoord, const GmMatrix& X, const GmVector& N, const GmMatrix& J, GmMatrix& Bu);
  
  virtual double fillBuMatrix(const GmElement* e, const GmShape* shape, const GmVector& ncoord, const GmMatrix& X, const GmVector& N, const GmMatrix& J, GmMatrix& Bu, const GmMatrix& F) = 0;

  virtual void fillBnlMatrix(const GmElement* e, const GmShape* shape, const GmVector& ncoord, const GmMatrix& X, const GmVector& N, const GmMatrix& J, GmMatrix& Bnl, const GmMatrix& F) = 0;

  virtual bool DeformationGradient(const GmElement* e, const GmMatrix& X0, const GmMatrix& Xt, const GmVector& N, const GmMatrix& J0, const GmMatrix& Jt, GmMatrix& F) = 0;

protected:
  // See comments on the base class
  virtual bool checkAndLoadAttributeAccessors(LuaTable& nodeTable, LuaTable& gaussTable);

  // Fills the Green-Lagrange strain vector form
  void GreenLagrangeStrain(const GmMatrix& F, GmVector& Ev, unsigned sc);

  // Fills the Euler-Almansi strain vector form
  void EulerAlmansiStrain(const GmMatrix& F, GmVector& ev, unsigned sc);

  // Fills the Logarithm strain vector form
  void LogarithmStrain(const GmMatrix& F, GmVector& elog, unsigned sc);

  void NLStressStateMatrix(const GmVector& S, GmMatrix& SvM, const GmMatrix& F, int d);

  // Fills the Cauchy stress vector form
  void CauchyStress(const GmVector& S, GmVector& sv, const GmMatrix& F, int d);

  //void FillDofIndexPosition(const GmElement * e, QList<int>& iUe);
  
  GmMatrix     _F;            
};

#endif