shefplast.h

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

#ifndef SHEFPLAST_H
#define SHEFPLAST_H

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

/**
   This class defines Sheffield plastic material model

   14.1.2004
   Frederic Dufour, Tomas Koudelka
*/

class shefplast
{
 public:
  shefplast (void);
  ~shefplast (void);
  void read (XFILE *in);

  void compute_khat(matrix &sig, vector &q);
  void compute_rhoc(matrix &sig);
  void compute_r(matrix &sig);
  void compute_classvar (matrix &sig, vector &q);
  double yieldfunction (matrix &sig, vector &q);
  void deryieldfsigma (long ipp,matrix &sig, vector &q, matrix &dfds,long ido);
  void dderyieldfsigma (matrix &ddfds, strastrestate ssst);
  void derqdsigma (matrix &sig, matrix &drds);
  void derpotsigma (long ipp,matrix &sig, vector &q, matrix &dgds,long ido);
  void deryieldfq(matrix &sig, vector &q, vector &dq);
  void derkhatds (long ipp,matrix &sig, vector &q, matrix &dkhatds,long ido);
  //  void der_q_gamma(long ipp, matrix &sig, vector &eps, vector &qtr, vector &dqdg);
  //  double plasmodscalar(long ipp, matrix &sig, vector &eps, vector &qtr);
  void updateq(long ipp, vector &eps, matrix &sig, vector &q,long ido);
  //  void updateq(vector &epsp, strastrestate ssst, vector &q);
  void matstiff (matrix &d, long ipp, long ido);
  void tangentstiff (matrix &d, matrix &td, long ipp,long ido);
  void nlstresses (long ipp,long ido);
  void updateval (long ipp, long ido);
  double give_consparam (long ipp, long ido);

  void stress_return (long ipp,double &lambda,vector &k,vector &eps,vector &epsp, long ido);
  void compzeta();
  double hardening(long ipp,matrix &sigtens,vector &q,long ido);
  void numdiff_dfdsds(long ipp,matrix &sigtens,vector &q,matrix &dfdsds,long ido);
  void numdiff_dfdsdsc(long ipp,matrix &sigtens,vector &q,matrix &dfdsds,long ido);
  void numdiff_dfdsdk(long ipp,matrix &sigtens,vector &q,vector &dfdsdk,long ido);
  void numdiff_dfdsdkc(long ipp,matrix &sigtens,vector &q,vector &dfdsdk,long ido);
  void numdiff_dhds(long ipp,matrix &sigtens,vector &q,vector &dhds,long ido);
  void numdiff_dhdsc(long ipp,matrix &sigtens,vector &q,vector &dhds,long ido);
  void numdiff_dhdk(long ipp,matrix &sigtens,vector &q,double &dhdk,long ido);
  void numdiff_dhdkc(long ipp,matrix &sigtens,vector &q,double &dhdk,long ido);

  double maxim (double a,double b);

  /// 
  double rc;
  ///
  double rt;
  ///
  double gamma;
  ///
  double p;
  ///
  double a;
  ///
  double k0;
  /// 
  double alpha;
  ///
  double ah, bh, ch;
  /// initial value of hardening parameter
  double kh0;

  // Internal temporary variables

  /// stress deviator
  matrix dev;
  /// the first invariant of the stress tensor
  double i1s;
  /// the second invariant of the stress deviator
  double j2s;
  /// the third invariant of the stress deviator
  double j3s;
  double khat, rhoc, r, m, xi, ft;
  double b0, b1, c0, c1, c2, c3, c4, c5, theta, d0, d1, d2;
  
  
  double zeta;

  ///  stress return algorithm
  strretalg sra;
};

#endif