esuf_EllipticCylinder.h

Go to the documentation of this file.

//********************************************************************************************************
// code:                          ###  ###  #####   ####  ##  ##
//                       ##  ##   ########  ##     ##  ## ##  ##
//                       ##  ##   ## ## ##  ###    ##     ######
//                       ##  ##   ##    ##  ##     ##  ## ##  ##
//                       #####    ##    ##  #####   ####  ##  ##
//                       ##
//
// name:           eshelbySoluUniformEllipticCylinder.h
// description:    header file of the function returning the Eshelby solution of an elliptic-cylindrical inclusion
//                 loaded by the uniform strain/stress field
// author(s):      Jan Novak, Tomas Janda

// last edit:      14-1-2014
// 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_ESHELBY_SOLU_UNIFORM_FIELD_ELLIPTIC_CYLINDER_H
#define MUMECH_ESHELBY_SOLU_UNIFORM_FIELD_ELLIPTIC_CYLINDER_H

#include "esuf.h"
#include "types.h"


namespace mumech {

class eshelbySoluUniformFieldEllipticCylinder : public eshelbySoluUniformField
{
 public:
  eshelbySoluUniformFieldEllipticCylinder (const InclusionRecord3D *i) : eshelbySoluUniformField (i) { }
  virtual ~eshelbySoluUniformFieldEllipticCylinder () {}
  
  //******************************************************************************************************
  //Function gives a component of the Eshelby tensor
  double eshelbyTensCompUniformField( const double sort_a[3], const double eInt[13], double nu,
                              EshelbyTensComponent flag );
  //******************************************************************************************************
  //Function initialize an 1D array by Eshelby tensor components
  void eshelbyTensUniformField( double eshTens[12], const double sort_a[3], const double eInt[13]);
  //******************************************************************************************************
  //Function gives D_ijkl(x) generalized Eshelby's tensor. Function requiers "constructor2"
  void giveDijkl( double D[36], const double S[12], const double J[13], const double dJi[9], const double dJij[27],
                  const double ddJi[27], const double ddJij[81], const double sort_a[3],  const double x[3]);

 protected:
  
  //******************************************************************************************************
  //Function gives the displacement perturbation tensor of internal fields.
  //Function requiers "constructor2"
  void giveLijkINT( double Lint[18], const double J[13], const double sort_a[3], const double x[3] );
  //******************************************************************************************************
  //Function gives the displacement perturbation tensor of external fields.
  //Function requiers "constructor2"
  void giveLijkEXT(double Lext[18], const double Lint[18], const double dJi[9], const double dJij[27], const double sort_a[3],
                    const double x[3]);
  //******************************************************************************************************
  //Function gives S_ijkl(x) position dependent Eshelby's tensor. (This is not
  //the strain perturbation tensor! Just its part ). The returned members of the
  //S_ijkl matrix are in tensorial (theoretical) notation!
  //Function requiers "constructor2"
  void giveSijkl( double S[36], const double J[13], const double sort_a[3], double nu, bool newFormulation );
  //******************************************************************************************************

};//end of class declaration

} // end of namespace mumech

#endif

/*end of file*/