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) const;
  
  virtual FemResultType fillElementData(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors) = 0;
  virtual bool fixedNodalForcesBc(QVector<int>& nodes, QVector<int>&dof, QVector<double>& values) const;
  virtual bool fixedNodalDofsBc  (QVector<int>& nodes, QVector<int>&dof, QVector<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,              
    STRAIN_T_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,          
    TREF_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,
    ET_GA_ID,           
 
    // ------    NO ADDING BELOW THIS LINE
    NUM_GA_IDS,  
  };
 
  enum ElementPropertyIds
  {
    THICKNESS_ID,  
    BF_ID,         
    GRAV_ID,       
    DENS_ID,       
    DENSB_ID,      
    ALPHA_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 nodal thermal increment
  void fillElementThermalIncrement(const GmElement* e, GmVector& dTe);
 
  // 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) const { 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; }
 
  virtual bool evaluatesThermalStrain(const GmElement* e, GmVector& strainT, GmVector& NT, GmVector& dT, double alpha) const;
 
  inline GmMatrix& dep() { return _Dep.localData(); }
 
protected:
  GmpMechanicPoint*    _mechanicPoint;  
  GmTLS<GmMatrix, true>      _Dep;            
  bool                 _setInitialCond; 
 
private:
  const char*          _pluginType;     
};
 
#endif