esuf_Sphere.h

Go to the documentation of this file.

//********************************************************************************************************
// code:                          ###  ###  #####   ####  ##  ##
//                       ##  ##   ########  ##     ##  ## ##  ##
//                       ##  ##   ## ## ##  ###    ##     ######
//                       ##  ##   ##    ##  ##     ##  ## ##  ##
//                       #####    ##    ##  #####   ####  ##  ##
//                       ##
//
// name:           eshelbySoluUniformFieldSphere.h
// author(s):      Jan Novak, Lukasz Kaczmarczyk
// 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_SPHERE_H
#define MUMECH_ESHELBY_SOLU_UNIFORM_FIELD_SPHERE_H

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

namespace mumech {

class eshelbySoluUniformFieldSphere : public eshelbySoluUniformField
{
 public:
  eshelbySoluUniformFieldSphere (const InclusionRecord3D *i) : eshelbySoluUniformField (i) { }
  virtual ~eshelbySoluUniformFieldSphere () {}
  
  
  //******************************************************************************************************
  // Function gives a component of the Eshelby tensor
  double giveEshelbyTensCompUniformField( EshelbyTensComponent flag, double nu, double a );
  //******************************************************************************************************
  // Function initialize an 1D array by Eshelby tensor components
  void eshelbyTensUniformField (double eshTens[12], const double sort_a[3], const double eSInt[13]);

  //******************************************************************************************************
  // Function gives the inverse of the Eshelby tensor
  void giveEshelbyTensorInverse(double SInv[12], const double S[] );
  //******************************************************************************************************
  // Function gives D_ijkl(x) generalized Eshelby's tensor. Function requiers "constructor1".
  void giveDijkl(double D[12], const double S[], const double J[], const double dJi[], const double dJij[],
                  const double ddJi[], const double ddJij[], const double sort_a[3],  const double x[]);
  
 protected:
  
  //******************************************************************************************************
  //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 the displacement perturbation tensor of internal fields.
  //Function requiers "constructor2"!
  void giveLijkINT(double Lint[18], const double J[13], const double 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!
  void giveSijkl(double S[36], const double J[], const double sort_a[], double nu , bool newFormulation);
  //******************************************************************************************************

};

} // end of namespace mumech

#endif

/*end of file*/