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& cellTable, 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) const;
 
  // Fills the Logarithm strain vector form
  void LogarithmStrain(const GmMatrix& F, GmVector& elog, unsigned sc) const;
 
  void NLStressStateMatrix(const GmVector& S, GmMatrix& SvM, const GmMatrix& F, int d) const;
 
  // Fills the Cauchy stress vector form
  void CauchyStress(const GmVector& S, GmVector& sv, const GmMatrix& F, int d) const;
 
  //void FillDofIndexPosition(const GmElement * e, QList<int>& iUe);
  
  GmTLS<GmMatrix, true>     _F;            
};
 
#endif