lemMaterialElastic.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_MATERIAL_ELASTIC_H_
#define _GEMA_LEM_MATERIAL_ELASTIC_H_

#include <gmCellAccessor.h>
#include <gmpFemPhysics.h>
#include <gmMathUtils.h>

#include "lemMechanicalMaterial.h"
#include "lemMechanicPoint.h"

class  LemMaterialElastic : public LemMechanicalMaterial
{
protected:
  enum ElementPropertyIds
  {
    E_ID,             
    NU_ID,            

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

public:

  LemMaterialElastic(int typeIndex, QString typeName, const GmLogCategory& logger)
    : LemMechanicalMaterial(typeIndex, typeName, logger) {}

  virtual ~LemMaterialElastic() {}

  static GmpFemPhysicsCommonMaterial* instance(GmSimulationData* simulation, int typeIndex, QString typeName, const GmLogCategory& logger)
  {
    Q_UNUSED(simulation);
    return new LemMaterialElastic(typeIndex, typeName, logger);
  }

  // Returns a map with material associated properties. 
  virtual const QVariantMap* materialMetaDataMap();

  // Returns the stresses according to the material behavior adopted
  virtual bool returnMapping(const GmElement* e, GmMatrix& Dep, const lemMechanicPoint* mp, const GmVector* coord, unsigned nc, unsigned ips) const;

  // Returns the constitutive tangent matrix
  virtual void tangentModulus(const GmElement* e, GmMatrix& Dep, const lemMechanicPoint* mp, const GmVector* coord, unsigned nc, unsigned ips) const;

  virtual bool isIsotropic() const { return true; }

  virtual double elasticModulus(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(E_ID)->scalarValueAt(e, coord, ip);
  }
  virtual double poissonRatio(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(NU_ID)->scalarValueAt(e, coord, ip);
  }

  //virtual void stiffness(const GmElement*, unsigned, GmMatrix&, const GmVector*, unsigned, int, const GmVector* = NULL) const;
  //virtual void flexibility(const GmElement*, unsigned, GmMatrix&, const GmVector*, unsigned, int, const GmVector* = NULL) const;
  virtual double bulkModulus(const GmElement*, const GmVector*, int) const;
  virtual double shearModulus(const GmElement*, const GmVector*, int) const;
};

#endif