dkt.h

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

#ifndef DKT_H
#define DKT_H

#include "alias.h"
struct matrix;
struct vector;
struct ivector;

class dktelem
{
 public:
  dktelem (void);
  ~dktelem (void);
  void eleminit (long eid);
  void geom_matrix (matrix &gm,vector &x,vector &y,vector &l);
  void transf_matrix (ivector &nodes,matrix &tmat);
  void dbmat (matrix &d,matrix &db,double t);
  void res_stiffness_matrix (long eid,matrix &sm);
  void stiffness_matrix (long eid,long ri,long ci,matrix &sm,vector &x,vector &y);
  void nodecoord (vector &xi,vector &eta);
  void appval (vector &l, long fi,long nc,vector &eps,double **val);

  void res_mainip_strains (long lcid,long eid);
  void mainip_strains (long lcid,long eid,long ri,long ci,vector &x,vector &y,vector &r);
  void nod_strains (long lcid,long eid,long ri,long ci);
  void elem_strains (double **stra,long lcid,long eid,long ri,long ci);
  void appstrain (long lcid,long eid,vector &l,long fi,long ncomp,vector &eps);
  void allip_strains (double **stra,long lcid,long eid,long ri,long ci);
  void strains (long lcid,long eid,long ri,long ci);

  void res_allip_stresses (long lcid,long eid);
  void mainip_stresses (long lcid,long eid,long ri,long ci);
  void nod_stresses (long lcid,long eid,long ri,long ci);
  void elem_stresses (double **stra,double **stre,long lcid,long eid,long ri,long ci);
  void appstress (long lcid,long eid,vector &l,long fi,long ncomp,vector &sig);
  void allip_stresses (double **stre,long lcid,long eid,long ri,long ci);
  void stresses (long lcid,long eid,long ri, long ci);
  void nodeforces (long eid,long *le,double *nv,vector &nf);
  void areaforces (long eid,double *nv,vector &nf);
  void inicipval(long eid, long ri, long ci, matrix &nodval, inictype *ictn);

  //  number of DOFs on the element
  long ndofe;
  //  number of nodes on one element
  long nne;
  ///  total number of integration points on element
  long tnip;
  //  array containing numbers of components of stress and strain tensors
  long *ncomp;
  long *cncomp;
  //  total number of components of stress and strain tensors
  long tncomp;
  //  number of approximated functions on the element
  long napfun;
  //  number of edges on one element
  long ned;
  //  number of nodes on one edge
  long nned;
  //  order of integration of stiffness matrix
  long **intordsm;
  //  array of numbers of integration points in sets
  long **nip;
  //  number of blocks
  long nb;
  //  order of integration for mass matrix
  long intordmm;
  //  stress/strain state
  strastrestate ssst;
};

#endif