gmpMechanicalPhysics.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_H_
#define _GEMA_PLUGIN_MECHANICAL_PHYSICS_H_

#include "gmpMechanicalConfig.h"
#include <gmpFemPhysicsCommon.h>
#include <gmSimulationData.h>
#include <gmStateDump.h>

#include "gmpMechanicPoint.h"
#include "gmpMechanicalContact.h"

class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpMechanicalPhysics : public GmpFemPhysicsCommon
{
public:

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

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

  // Comments on the base class
  virtual const char* pluginType() const { return _pluginType; }
  // Comments on the base class
  virtual const QVariantMap* physicsMetaDataMap();

  // \brief Tell the world that the mechanical physics supports state dumping.
  virtual bool supportsStateDumping() { return true; }

  // supported bc
  //virtual bool supportsBcType(GmBcType type) const { assert(_attrMap);  return _supportedTypes[type]; }

  virtual unsigned stress_state(void);
  
  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 FemResultType fillElementDataForBc(const GmElement* e, const GmBoundaryCondition* bc, int bcIndex, int bcListIndex,
                                             int border, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);  
  virtual FemResultType fillElementDataForLoads(const GmElement* e, int loadId, GmpFemVectorSet& elemVectors);

  virtual FemResultType fillElementDataPressBc(const GmElement* e, const GmBoundaryCondition* bc, int bcIndex, int bcListIndex,
                                          int border, GmMatrix& elemMat, GmVector& elemVec);
  virtual FemResultType fillElementDataHydPressBc(const GmElement* e, const GmBoundaryCondition* bc, int bcIndex, int bcListIndex,
          int border, GmMatrix& elemMat, GmVector& elemVec);
  virtual FemResultType fillElementDataWaterSheepPressBc(const GmElement* e, const GmBoundaryCondition* bc, int bcIndex, int bcListIndex,
          int border, GmMatrix& elemMat, GmVector& elemVec);
  virtual FemResultType fillContactData(const GmContactBoundaryCondition* cbc, int bcIndex1, int bcIndex2,
    GmMatrixDof& dofMap, GmpFemMatrixSet& matrixSet, GmpFemVectorSet& vecSet);
  // Calculates derived results 
  virtual bool calcDerivedResults(bool nonLinearSolver);
  // Calculate material derived results
  virtual bool calcMaterialDerivedResults() const;
  enum PhysicsAttributeIds
  {
    ISOPARAMETRIC_ID,        
    STRESS_MODE_ID,          
    STRAIN_MODE_ID,          
    P_STRAIN_MODE_ID,        
    BODY_FORCE_ID,           
    CGRAV_ID,                
    GEOMETRIC_ID,            
    RF_MODE_ID,              
    RMAPPINGINFO_ID,         
    K0MODE_ID,               
    STRAIN_RESET_ID,         
    CONTACT_ID,         

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

  enum StateVarIds
  {
    U_ID,         

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

  enum ExternalLoadsIds
  {
    PRESSURE_ID,         
    THERMAL_ID,          

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

  enum NodeAttributeIds
  {
    S_NA_ID,           
    E_NA_ID,           
    RF_NA_ID,          
    CF_NA_ID,          
    CD_NA_ID,          
    CIV_NA_ID,         
    CIVOLD_NA_ID,      

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

  enum GaussAttributeIds
  {
    S_GA_ID,            
    SOLD_GA_ID,         
    RS_GA_ID,           
    E_GA_ID,            
    EOLD_GA_ID,         
    YSR_GA_ID,          
    YSROLD_GA_ID,       
    ST_GA_ID,           
    STOLD_GA_ID,        
    ITER_GA_ID,         
    ITEROLD_GA_ID,      
    SUBSTEP_GA_ID,      
    SUBSTEPOLD_GA_ID,   
    LINESEARCH_GA_ID,   
    LINESEARCHOLD_GA_ID,

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

  enum ElementPropertyIds
  {
    THICKNESS_ID,  
    BF_ID,         
    GRAV_ID,       
    DENS_ID,       
    DENSB_ID,      

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

  enum BoundaryConditionIds
  {
    FIXED_U_BC_ID,         
    FIXED_CF_BC_ID,        
    PL_BC_ID,              
    HL_BC_ID,              
    WP_BC_ID,              

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

  enum BoundaryConditionValueIds
  {     
    BC_UX_ID,     
    BC_UY_ID,     
    BC_UZ_ID,     
    BC_CF_ID,     
    BC_PL_ID,     
    BC_HLR_ID,    
    BC_HL_ID,     
    BC_NRL_ID,    
    BC_SWF_ID,    

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

  enum ContactBCIds
  {
    LOCAL_CBC_ID,          
    GLOBAL_CBC_ID,         

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

  enum ContactBCValueIds
  {
    CBC_KCN_ID,     
    CBC_KSN_ID,     
    CBC_MU_ID,      

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

  virtual double fillBuMatrix(const GmElement* e, const GmShape* shape, const GmVector& ncoord, 
                              const GmMatrix& X, const GmVector& N, const GmMatrix& J, GmMatrix& Bu) = 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 axisymmetricFactor(const GmElement* e, const GmMatrix& X, const GmVector& N) { Q_UNUSED(e); Q_UNUSED(X); Q_UNUSED(N); return 1.0; }
  
  virtual bool isPlaneStress() { return false; }

  virtual bool isAxisymmetric() { return false; }

  GmpMechanicPoint* mechanicPoint() { return _mechanicPoint; }

  GmMatrix& dep() { return _Dep; }

  virtual GmVector timeSim() const;

protected:
  GmpMechanicalContact _mecContact;     
  GmpMechanicPoint*    _mechanicPoint;  
  GmMatrix             _Dep;            
  bool                 _setInitialCond; 

private:
  const char* _pluginType; 
};

#endif