gmpThermoFemPhysics.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_THERMO_FEM_PHYSICS_H_
#define _GEMA_PLUGIN_THERMO_FEM_PHYSICS_H_
 
#include <gmpFemPhysicsCommon.h>
 
#include "gmpThermoConfig.h"
 
class GMP_THERMO_PHYSICS_API_EXPORT GmpThermoFemPhysics: public GmpFemPhysicsCommon
{
public:
  GmpThermoFemPhysics(const char* pluginType, GmSimulationData* simulation, QString id, QString description, const GmLogCategory& logger);
  virtual ~GmpThermoFemPhysics();
 
  // Comments on the base class
  virtual const char* pluginName() const { return "ThermoFemPhysics"; }
 
  // Comments on the base class
  //virtual const char* pluginType() const { return "heat"; }
  // Comments on the base class
  virtual const char* pluginType() const { return _pluginType; }
 
  virtual const QVariantMap* physicsMetaDataMap();
 
  virtual bool beforeElementStiffnessLoop(const GmpFemMatrixSet& elemMatrices, const GmpFemVectorSet& elemVectors);
  virtual FemResultType fillElementData(const GmElement* e, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);
 
  virtual FemResultType fillElementDataForBc(const GmElement* e, const GmBoundaryCondition* bc, int bcIndex, int bcListIndex,
                                             int border, GmpFemMatrixSet& elemMatrices, GmpFemVectorSet& elemVectors);
 
  virtual bool fixedNodalDofsBc(QVector<int>& nodes, QVector<int>&dof, QVector<double>& values, bool* constantValues) const;
 
  virtual bool calcDerivedResults(bool nonLinearSolver);
  
  virtual bool supportsStateDumping() { return true; }
 
  void fillConductivityMatrix(GmMatrix& K, const GmElement* e, const GmVector* ip, int ipIndex);
 
  enum PhysicsAttributeIds
  {
    ISOPARAMETRIC_ID,           
    HAS_ELEMENT_GENERATION_ID,  
    HF_MODE_ID,                 
    HF_DIRECT_ID,               
  };
 
  enum StateVarIds
  {
    T_ID,   
  };
 
  enum NodeAttributeIds
  {
    Q_NA_ID,  
  };
 
  enum GaussAttributeIds
  {
    Q_GA_ID,     
  };
 
  enum ElementPropertyIds
  {
    K_ID,          
    THICKNESS_ID,  
    G_ID,          
    RHO_ID,        
    CP_ID,         
 
    // --- NO ADDING BELOW THIS LINE
    NUM_PROPER_IDS,  
  };
 
  enum BoundaryConditionIds
  {
    FIXED_T_BC_ID,          
    FIXED_HF_BC_ID,         
    CONVECTIVE_BC_ID,       
    POINT_GENERATION_BC_ID, 
  };
 
  enum BoundaryConditionValueIds
  {
    BC_T_ID,       
    BC_Q_ID,       
    BC_QDIR_ID,    
    BC_H_ID,       
    BC_TA_ID,      
    BC_G_ID,       
    BC_POS_ID,     
  };
 
  
protected:
  virtual void calcElementDerivedResultAtPoints(const GmElement* e, int resultId, const GmMatrix& evalPoints, bool evalAtIp, GmMatrix& result);
  virtual bool isAxisymmetric() { return false; }
 
  virtual double axisymmetricFactor(const GmElement* e, const GmMatrix& X, const GmVector& N) const = 0;
 
private:
  const char* _pluginType; 
  bool _hasElementGeneration;  
};
 
#endif