lemMechanicalPhysics.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_LEM_MECHANICAL_PHYSICS_H_
#define _GEMA_LEM_MECHANICAL_PHYSICS_H_

#include <gmpFemPhysicsCommon.h>
#include <gmSimulationData.h>
#include "lemMechanicPoint.h"


class  LemMechanicalPhysics : public GmpFemPhysicsCommon
{
public:

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

  // Comments on the base class
  virtual const char* pluginName() const { return "Lem"; }

  // Comments on the base class
  virtual const char* pluginType() const { return _pluginType; }
  
  virtual const QVariantMap* physicsMetaDataMap();
  
  virtual FemResultType fillElementData(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors) = 0;
  virtual bool fixedNodalForcesBc(QList<int>& nodes, QList<int>&dof, QList<double>& values) const;
  virtual bool fixedNodalDofsBc  (QList<int>& nodes, QList<int>&dof, QList<double>& values, bool* constantValues) const;

  virtual bool calcDerivedResults(bool nonLinearSolver);

  enum PhysicsAttributeIds
  {
    ISOPARAMETRIC_ID,        

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

  enum StateVarIds
  {
    U_ID,         
    W_ID,         

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

  enum NodeAttributeIds
  {
    S_NA_ID,           
    E_NA_ID,           

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

  enum GaussAttributeIds
  {
    S_GA_ID,           
    SOLD_GA_ID,        
    E_GA_ID,           
    EOLD_GA_ID,        
    ST_GA_ID,          
    STOLD_GA_ID,       

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

  enum ElementPropertyIds
  {
    AREA_ID,  

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

  enum BoundaryConditionIds
  {
    FIXED_U_BC_ID,         
    FIXED_CF_BC_ID,        

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

  enum BoundaryConditionValueIds
  {     
    BC_UX_ID,     
    BC_UY_ID,     
    BC_UZ_ID,     
    BC_CF_ID,     

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

  virtual double fillBuMatrix(const GmElement* e, const GmShape* shape, const GmVector& ncoord, const GmMatrix& X, const GmVector& N, const GmMatrix& J, GmMatrix& Bp) = 0;

  void fillElementDisplacements(const GmElement* e, GmVector& ue);
  
  // fills middle coordinates
  virtual void fillMiddleCoordinates(const GmElement* e, GmMatrix Xnode, GmMatrix& MX) const { Q_UNUSED(e); Q_UNUSED(MX); }
  
  // computes external load contribution at specific integration point
  //virtual bool externalLoadContribution(const GmElement* e, GmVector& load, QString sVar, const GmVector& ncoord);

  virtual double axisSymetricFactor(const GmElement* e, const GmMatrix& X, const GmVector& N) { Q_UNUSED(e); Q_UNUSED(X); Q_UNUSED(N); return 1.0; }
  
  virtual bool isPlaneStress() const;

  virtual unsigned stress_state(void);

  GmMatrix& dep() { return _Dep; }

protected:

  lemMechanicPoint* _mechanicPoint; 
  GmMatrix          _Dep;           

private:
  const char*   _pluginType;        
};

#endif