anisodam.h

Finite element program for mechanical problem. It can solve various problem in solid problem

#ifndef ANISODAM_H
#define ANISODAM_H

#include "iotools.h"
#include "alias.h"
#include "strretalg.h"
struct matrix;
struct vector;
struct atsel;

/**
  This class defines scalar isotropic damage material model.
  The different type of norms for the computing parameters of
  the damage function can be used.
  The order of internal variables is following :
  0 - actual value of param a
  1 - actual value of param At
  2 - actual value of param Ac
  3 - previous value of the damage driving force for volumetric damage
  4..6 - previous value of the damage driving forces in principal directions for deviatoric damage for tension
  7..9 - previous value of the damage driving forces in principal directions for deviatoric damage for compression
  10 - reached value of volumetric damage parameter
  11..13 - reached values of deviatoric damage parameters for tension
  14..16 - reached values of deviatoric damage parameters for compression
*/
class anisodam
{
 public:
  anisodam (void);
  ~anisodam (void);
  void   read (XFILE *in);
  double damdrvforce_vol(long ipp, vector &peps);
  void   damdrvforce_dev(long ipp, vector &peps, matrix &t, vector &pyt, vector &pyc);
  double loadfuncvol(long ipp, vector &peps, double d, double aa);
  void   loadfuncdev(long ipp, vector &peps, matrix &t, vector &damt, vector &damc, double aat, double aac, vector &lft, vector &lfc);
  double daminc_vol(long ipp, double y, double dy, double aa, double lf);
  void   pdaminc_dev(long ipp, vector &pyc, vector &pyt, vector &dyt, vector &dyc, double aat, double aac, vector &lft, vector &lfc, vector &dpdamt, vector &dpdamc);
  double dam_vol(long ipp, double y, double dy, double aa, double lf);
  void   pdam_dev(long ipp, vector &pyc, vector &pyt, vector &dyt, vector &dyc, double aat, double aac, vector &lft, vector &lfc, vector &pdamt, vector &pdamc);
  void   give_actual_param_a(long ipp, long ido, double &aa, double &aac, double &aat);
  void   initvalues(long ipp, long ido);
  void   matstiff (matrix &d,long ipp,long ido);
  void   elmatstiff (matrix &d,long ipp);
  void   nlstresses (long ipp, long im, long ido);
  void   updateval (long ipp,long im,long ido);


  /// correction of disipated energy switch
  corr_disip_en cde;
  /// material parameter A for volumetric damage
  double a;  
  /// material parameter B for volumetric damage
  double b;  
  /// initial treshold for damage driving force Y0 for volumetric damage
  double y0;  
  /// material parameter At for tension
  double at;  
  /// material parameter Bt for tension
  double bt;  
  /// initial treshold for damage driving force Y0 for tension
  double y0t;  
  /// material parameter Ac for compression
  double ac;  
  /// material parameter Bc for compression
  double bc;  
  /// initial treshold for damage driving force Y0 for compression
  double y0c;  
  /// fracture energy of volumetric damage 
  double gf;
  /// fracture energy of deviatoric damage for tension
  double gft;
  /// fracture energy of deviatoric damage for compression
  double gfc;
  
};

#endif