melement.h

Go to the documentation of this file.

//********************************************************************************************************
// code:                          ###  ###  #####   ####  ##  ##
//                       ##  ##   ########  ##     ##  ## ##  ##
//                       ##  ##   ## ## ##  ###    ##     ######
//                       ##  ##   ##    ##  ##     ##  ## ##  ##
//                       #####    ##    ##  #####   ####  ##  ##
//                       ##
//
// name:           melement.h
// author(s):      Ladislav Svoboda
// language:       C, C++
// license:        This program is free software; you can redistribute it and/or modify
//                 it under the terms of the GNU Lesser General Public License as published by
//                 the Free Software Foundation; either version 2 of the License, or
//                 (at your option) any later version.
//
//                 This program is distributed in the hope that it will be useful,
//                 but WITHOUT ANY WARRANTY; without even the implied warranty of
//                 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//                 GNU Lesser General Public License for more details.
//
//                 You should have received a copy of the GNU Lesser General Public License
//                 along with this program; if not, write to the Free Software
//                 Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
//********************************************************************************************************
#ifndef MUMECH_MELEMENT_H
#define MUMECH_MELEMENT_H

#include "types.h"


namespace mumech {

class Mesh;

class mElement
{
 public:
  long id;                
  const Mesh *M;          
  
  long nnodes;            
  long *nodes;            
 private:
  long type;              
  double volume;          
  
 public:
  long region;            
  
  double centroids[3];    
  
  double ***displc;        
  double ***strain;        
  double ***stress;        
  
  double **energy;         
  double **eqstress;       

  // termitovo
  // mesh predelat, aby nebyl pod problemem, protoze kazdy problem muze byt jine reseni, a ja chci na meshich zkoumat vice reseni
  // pak udelat mesh, udelat fci ktera do nej nacte mumech reseni, udelat fci kteeera nacte ansys reseni, spocitat chybu ...

  double **error   ;       




 public:
  mElement (long i, const Mesh *m);
  mElement (long i, const Mesh *m, const mElement *src);
  virtual ~mElement();
  
  long give_type (void) const { return type; }
  void set_type (long t);
  
  void compute_centroids(void);
  
  void allocate_fields (int pid);
  
  double give_volume(void);
  double give_energy (int pid, int lc);
  
  void save_error (double val, int pid, int lc);
  
  void giveReducedStiffMatrix (double *C, int pid) const;
  
  void give_strain (double *strain, int pid, int lc) const;
  void give_stress (double *stress, int pid, int lc) const;
  
  //double give_equiv_stress (int lc);
  
  
};//end of class declaration

} // end of namespace mumech

#endif // end of MUMECH_MELEMENT_H

/*end of file*/