microsim.h

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

#ifndef MICROSIM_H
#define MICROSIM_H

#include "iotools.h"
#include "matrix.h"
#include "vector.h"
#include "intpoints.h"



class microSIM
{
 public:
  microSIM (void);
  ~microSIM (void);
  void read (XFILE *in);
  void matstiff (matrix &d,long ipp, long ido);
  void nlstresses (long ipp, long ido);
  void updateval(long ipp, long im, long ido);

  //material parameters
  long numberOfMicroplanes;
  double e,nu;

 protected:
  vector  microplaneWeights; //(numberOfMicroplanes);
  matrix projN;
  vector kronecker; //(6);
  double k1,k2,k3,k4,c3,c20,c1,c2,c4,c5,c6,c7,c8,c9,c10,c11,c12,c13,c14,c15,c16,c17,c18,c19,mu,ev,ed,et;
  void initializeData (long numberOfMicroplanes );
  inline double macbra(double x);
  inline double FVplus(double epsV);
  inline double FVminus (double epsV);
  inline double FDminus(double epsD);
  inline double FDplus(double epsD);
  inline double FN(double epsN,double sigmaV);
  inline double maxim (double a,double b);
  inline double minim (double a,double b);
};

#endif