gmpMaterialSolidifyingKelvinChain.h

Go to the documentation of this file.

/************************************************************************
**
** Copyright (C) 2019 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_MECHANICALMATERIAL_SOLIDIFYINGKELVINCHAIN_H_
#define _GEMA_PLUGIN_MECHANICALMATERIAL_SOLIDIFYINGKELVINCHAIN_H_

#include "gmpMechanicalConfig.h"
#include <gmpFemPhysics.h>
#include <gmMathUtils.h>

#include "gmpMechanicalMaterial.h"
#include "gmpMechanicPoint.h"

class GMP_MECHANICAL_PHYSICS_API_EXPORT GmpMaterialSolidifyingKelvinChain : public GmpMechanicalMaterial
{
protected:
  enum ElementPropertyIds
  {
    E0_ID,             
    NU_ID,             
    E_CHAIN_ID,        
    RT_CHAIN_ID,       
    ALPHA_AGING_ID,    

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

  enum GaussAttributeIds
  {
    ECC_GA_ID,         
    ECC_OLD_GA_ID,     

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

public:

  GmpMaterialSolidifyingKelvinChain(int typeIndex, QString typeName, const GmLogCategory& logger)
    : GmpMechanicalMaterial(typeIndex, typeName, logger) {}

  virtual ~GmpMaterialSolidifyingKelvinChain() {}

  static GmpFemPhysicsCommonMaterial* instance(GmSimulationData* simulation, int typeIndex, QString typeName, const GmLogCategory& logger)
  {
    Q_UNUSED(simulation);
    return new GmpMaterialSolidifyingKelvinChain(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 mechanicalConstitutiveModel(const GmElement* e, GmMatrix& Dep, const GmpMechanicPoint* mp, const GmVector* coord, const GmVector& Time, unsigned sc, bool ips) const;

  virtual bool isIsotropic() const { return true; }

  virtual double elasticModulus(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(E0_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 GmVector elasticModulusChain(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE();
    GmVector data;
    GmCellAccessor* Acc = propertyAc(E_CHAIN_ID);
    if (Acc != NULL) {
      vectorPropertyValue(e, Acc, data, coord, ip);
    }
    return data;
  }

  virtual GmVector retardationTimeChain(const GmElement* e, const GmVector* coord, int ip) const
  {
    S_TRACE();
    GmVector data;
    GmCellAccessor* Acc = propertyAc(RT_CHAIN_ID);
    if (Acc != NULL) {
      vectorPropertyValue(e, Acc, data, coord, ip);
    }
    return data;
  }

  virtual double alphaAging(const GmElement* e, const GmVector* coord, int ip) const
  {
    return propertyAc(ALPHA_AGING_ID)->scalarValueAt(e, coord, ip);
  }

  void strainInteg(double E0, double nu, const GmVector& Ei, const GmVector& tau, double dTime, const GmVector& dStrain, const GmVector& stress, bool ips, double& Cve, GmVector& stress_new, GmMatrix& strainCreep) const;

  // Checks material loaded data. 
  virtual bool checkLoadedData(const GmElement* e) const;
};

#endif