types.h

Go to the documentation of this file.

//********************************************************************************************************
// code:                          ###  ###  #####   ####  ##  ##
//                       ##  ##   ########  ##     ##  ## ##  ##
//                       ##  ##   ## ## ##  ###    ##     ######
//                       ##  ##   ##    ##  ##     ##  ## ##  ##
//                       #####    ##    ##  #####   ####  ##  ##
//                       ##
//
// name:           types.h
// author(s):      Jan Novak
// 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_TYPES_H
#define MUMECH_ESHELBY_SOLU_TYPES_H

// There is conflict name XMLDocument in msxml and tinyxml.
// Solution 1:
// This define have to be before first win32 include.
#define WIN32_LEAN_AND_MEAN

//included libraries
#include <complex>

#include "gelib.h"


#define VERSION 1.34

//#define TOM_DEBUG

#define _4PI_  12.56637061435917295384


//local constant and macros definitions
#define _SELF_BALANCE_NORM_LIMIT_       1e-6

using namespace gelibspace;

namespace mumech {


enum classID {
  classVoid,
  
  classProblem,
  classProblemFEM
};

enum HomogenizationType { HT_Void, HT_RegGrid, HT_Dilute, HT_MoriTanaka, HT_SelfConsistent, HT_ReussBound, HT_VoightBound, HT_DifferentialScheme };


enum rotation { CLOCKWISE = 1, COUNTER_CLOCKWISE = 2 };

//********************************************************************************************************
enum SBM { SBM_VOID, SBM_SB, SBM_SB1, SBM_SB0 };

//********************************************************************************************************
//enum AngleType { AT_Void, AT_DEG, AT_RAD };

//********************************************************************************************************
enum LCmode { LCM_VOID, LCM_SINGLE, LCM_MULTIPLE };

//********************************************************************************************************
enum SBAtype {
  SBAT_VOID = 0,      
  SBAT_ORIGINAL,      
  SBAT_STANDA,        
  SBAT_MESHLESS       
};

//********************************************************************************************************
enum PFCmode { PFCM_VOID = 0, PFCM_FULL = 1, PFCM_OPTIMIZED = 2 };

//********************************************************************************************************
enum PointPositionFlag { PPF_VOID = 0, PPF_INT_POINT = 1, PPF_EXT_POINT = 2 };

//********************************************************************************************************
// termitovo: 2d oproti 3d ma jine napeti ve smyku, udelej priklad z zkontroluj to, jestli tam nefaliruje dvojka
// termitovo: strain a stress a tim padem matice tuhosti se ukladaji v STRN_THEORETICAL_FEEP notaci
//            matice tuhosti se uklada v 36to12 redukovane forme, ale to je mozne jen za predpokladu, ze matice asi isotropni, tj. cleny C_14 atd. jsou nulove
//            ovsem je otazka, jak je to napr. s matici D, nevim v jake je notaci a jestli se da taky takto redukovat, ale snad ano, kdyz to chodi
//
enum STRNotation { // Symetric Tensor Reduction Notation
  STRN_VOID = 0,

  STRN_THEORETICAL_ROW,

  STRN_THEORETICAL_FEEP,

  STRN_VOIGT_FEEP,

  STRN_MANDEL_FEEP

};


//********************************************************************************************************
enum InclusionGeometry {
  IS_VOID = 0,
  IS_ELLIPSOID = 1,        IS_SPHERE = 2,
  IS_OBLATE_SPHEROID = 8,  IS_PROLATE_SPHEROID = 9,
  
  //IS_FLAT_ELLIPSOID = 7,    IS_PENNY = 5,
  //IS_ELLIPTIC_CYLINDER = 3, IS_CYLINDER = 4,
  
  IS_ELLIPSE = 31, IS_CIRCLE = 32
  //IS_STRIP = 33, IS_CRACK = 34
};

inline const char* IST_e2s (InclusionGeometry ig)
{
  switch (ig) {
  case IS_ELLIPSOID:         return "ellipsoid";     break;
  case IS_SPHERE:        return "sphere";            break;
  case IS_OBLATE_SPHEROID:   return "oblate spheroid";   break;
  case IS_PROLATE_SPHEROID:  return "prolate spheroid";  break;
  //case IS_ELLIPTIC_CYLINDER: return "elliptic cylinder"; break;
  //case IS_CYLINDER:        return "cylinder";          break;
  //case IS_PENNY:       return "penny";             break;
  //case IS_FLAT_ELLIPSOID:    return "flat ellipsoid";    break;
  case IS_ELLIPSE:       return "ellipse";           break;
  case IS_CIRCLE:        return "circle";            break;
  //case IS_STRIP:       return "strip";             break;
  //case IS_CRACK:             return "crack";             break;
  default: _errorr ("Unknown inclusion shape");
  }
  return NULL;
}



//Elliptic integrals (potentials) components
//********************************************************************************************************
typedef enum{ //do not change, the complex derivatives depend directly on individual integer values
  _I1_  = 1,  _I2_  = 2,  _I3_  = 3,
  _I11_ = 11, _I12_ = 12, _I13_ = 13,
  _I21_ = 21, _I22_ = 22, _I23_ = 23,
  _I31_ = 31, _I32_ = 32, _I33_ = 33,
  _I_,
}ellipticIntegralComponent;

//definition of derivative directions
//********************************************************************************************************
typedef enum{
  //first derivatives
  _x1_ = 1, _x2_ = 2, _x3_ = 3, _empty_,
  //seccond derivatives
  _x1x1_, _x1x2_, _x1x3_,
  _x2x1_, _x2x2_, _x2x3_,
  _x3x1_, _x3x2_, _x3x3_,
  //mixed or homogeneous directions of 2nd derivative
  _MIXED_, _HOMOGENEOUS_,
}derivativeDirection;

//Components of the Eshelby tensor
//********************************************************************************************************
typedef enum{
  _S1111_, _S1122_, _S1133_, _S1123_, _S1113_, _S1112_,
  _S2211_, _S2222_, _S2233_, _S2223_, _S2213_, _S2212_,
  _S3311_, _S3322_, _S3333_, _S3323_, _S3313_, _S3312_,
  _S2311_, _S2322_, _S2333_, _S2323_, _S2313_, _S2312_,
  _S1311_, _S1322_, _S1333_, _S1323_, _S1313_, _S1312_,
  _S1211_, _S1222_, _S1233_, _S1223_, _S1213_, _S1212_,
}EshelbyTensComponent;

//Components of the perturbation tensor
//********************************************************************************************************
typedef enum{
  //strain and stress tensor perturbation component s
  _D1111_ = 1111, _D1122_ = 1122, _D1133_ = 1133, _D1123_ = 1123, _D1113_ = 1113, _D1112_ = 1112,
  _D2211_ = 2211, _D2222_ = 2222, _D2233_ = 2233, _D2223_ = 2223, _D2213_ = 2213, _D2212_ = 2212,
  _D3311_ = 3311, _D3322_ = 3322, _D3333_ = 3333, _D3323_ = 3323, _D3313_ = 3313, _D3312_ = 3312,
  _D2311_ = 2311, _D2322_ = 2322, _D2333_ = 2333, _D2323_ = 2323, _D2313_ = 2313, _D2312_ = 2312,
  _D1311_ = 1311, _D1322_ = 1322, _D1333_ = 1333, _D1323_ = 1323, _D1313_ = 1313, _D1312_ = 1312,
  _D1211_ = 1211, _D1222_ = 1222, _D1233_ = 1233, _D1223_ = 1223, _D1213_ = 1213, _D1212_ = 1212,
  //displacement field perturbation components
  _L111_ = 111, _L122_ = 122, _L133_ = 133,
  _L112_ = 112, _L121_ = 121, _L123_ = 123, _L132_ = 132, _L113_ = 113, _L131_ = 131,
  _L211_ = 211, _L222_ = 222, _L233_ = 233,
  _L212_ = 212, _L221_ = 221, _L223_ = 223, _L232_ = 232, _L213_ = 213, _L231_ = 231,
  _L311_ = 311, _L322_ = 322, _L333_ = 333,
  _L312_ = 312, _L321_ = 321, _L323_ = 323, _L332_ = 332, _L313_ = 313, _L331_ = 331,
}perturbTensComponent;

//macros for rotating the Eshelby tensor S_ijkl
//according the ellipsoidal semiaxes' dimensions (default:a1>a2>a3)
//********************************************************************************************************
/* 1st row */
#define _11_         0
#define _12_         1
#define _13_         2
/* 2nd row */
#define _21_         3
#define _22_         4
#define _23_         5
/* 3rd row */
#define _31_         6
#define _32_         7
#define _33_         8
/* 4th row */
#define _44_         9
/* 5th row */
#define _55_         10
/* 6th row */
#define _66_         11

//the same principle as for the perturbation strain
/* 1st row */
#define _S1111swap_       0
#define _S1122swap_       1
#define _S1133swap_       2
/* 2nd row */
#define _S2211swap_       3
#define _S2222swap_       4
#define _S2233swap_       5
/* 3rd row */
#define _S3311swap_       6
#define _S3322swap_       7
#define _S3333swap_       8
/* 4th row */
#define _S1212swap_       9
/* 5th row */
#define _S2323swap_      10
/* 6th row */
#define _S1313swap_      11

//the same as _S__swap_, i.e. components of isotropic tensors having 12 non-zero components
/* 1st row */
#define __1111_           0
#define __1122_           1
#define __1133_           2
/* 2nd row */
#define __2211_           3
#define __2222_           4
#define __2233_           5
/* 3rd row */
#define __3311_           6
#define __3322_           7
#define __3333_           8
/* 4th row */
#define __1212_           9
/* 5th row */
#define __2323_          10
/* 6th row */
#define __1313_          11

//macros for rotating the strain perturbation tensor D_ijkl
//according the ellipsoidal semiaxes' dimensions (default:a1>a2>a3)
//*************************************************************************************************************
/* 1st row */
#define _1111_       0
#define _1122_       1
#define _1133_       2
#define _1112_       3
#define _1123_       4
#define _1113_       5
/* 2nd row */
#define _2211_       6
#define _2222_       7
#define _2233_       8
#define _2212_       9
#define _2223_       10
#define _2213_       11
/* 3rd row */
#define _3311_       12
#define _3322_       13
#define _3333_       14
#define _3312_       15
#define _3323_       16
#define _3313_       17
/* 4th row */
#define _1211_       18
#define _1222_       19
#define _1233_       20
#define _1212_       21
#define _1223_       22
#define _1213_       23
/* 5th row */
#define _2311_       24
#define _2322_       25
#define _2333_       26
#define _2312_       27
#define _2323_       28
#define _2313_       29
/* 6th row */
#define _1311_       30
#define _1322_       31
#define _1333_       32
#define _1312_       33
#define _1323_       34
#define _1313_       35

//macros for rotating the displacement perturbation tensor L_ijk
//according the ellipsoidal semiaxes' dimensions (default:a1>a2>a3)
//*************************************************************************************************************
/* 1st row */
#define _111_       0
#define _122_       1
#define _133_       2
#define _112_       3
#define _123_       4
#define _113_       5
/* 2nd row */
#define _211_       6
#define _222_       7
#define _233_       8
#define _212_       9
#define _223_       10
#define _213_       11
/* 3rd row */
#define _311_       12
#define _322_       13
#define _333_       14
#define _312_       15
#define _323_       16
#define _313_       17




// toto se pouzivalo pro closed penny, ted uz je to asi nepotrebne
// //small numbers and multipliers
// #ifndef _MUCH_BIGGER_SCALE_
// #define _MUCH_BIGGER_SCALE_   0.001
// #endif

#define _MIN_AXES_DIFF 1e-7
#define _MAX_AXES_DIFF 1e-6





//components' positions in Point.eInt[13] array
//********************************************************************************************************
#define _Int_        0
//--------------------
#define _Int1_       1
#define _Int2_       2
#define _Int3_       3
//--------------------
#define _Int11_      4
#define _Int12_      5
#define _Int13_      6
//--------------------
#define _Int21_      7
#define _Int22_      8
#define _Int23_      9
//--------------------
#define _Int31_      10
#define _Int32_      11
#define _Int33_      12
//--------------------



//various useful constants for manipulating InclusionRecord objects
//********************************************************************************************************
//macros vrom readVTK2.h


//Ferers-Dysons integrals' macros
//*******************************
#define _J               J[0]

#define J1               J[1]
#define J2               J[2]
#define J3               J[3]

#define J11              J[4]
#define J12              J[5]
#define J13              J[6]

#define J21              J[7]
#define J22              J[8]
#define J23              J[9]

#define J31              J[10]
#define J32              J[11]
#define J33              J[12]

//Derivatives of Ferers-Dysons' integrals
//***************************************
//first derivatives:
//------------------
//Ii_i:
#define J1_1             dJi[0]
#define J1_2             dJi[1]
#define J1_3             dJi[2]

#define J2_1             dJi[3]
#define J2_2             dJi[4]
#define J2_3             dJi[5]

#define J3_1             dJi[6]
#define J3_2             dJi[7]
#define J3_3             dJi[8]

//Iij_1:
#define J11_1            dJij[0]
#define J12_1            dJij[1]
#define J13_1            dJij[2]

#define J21_1            dJij[3]
#define J22_1            dJij[4]
#define J23_1            dJij[5]

#define J31_1            dJij[6]
#define J32_1            dJij[7]
#define J33_1            dJij[8]

//Iij_2:
#define J11_2            dJij[9]
#define J12_2            dJij[10]
#define J13_2            dJij[11]

#define J21_2            dJij[12]
#define J22_2            dJij[13]
#define J23_2            dJij[14]

#define J31_2            dJij[15]
#define J32_2            dJij[16]
#define J33_2            dJij[17]

//Iij_3:
#define J11_3            dJij[18]
#define J12_3            dJij[19]
#define J13_3            dJij[20]

#define J21_3            dJij[21]
#define J22_3            dJij[22]
#define J23_3            dJij[23]

#define J31_3            dJij[24]
#define J32_3            dJij[25]
#define J33_3            dJij[26]

//second derivatives:
//-------------------
//I1_ij:
#define J1_11            ddJi[0]
#define J1_12            ddJi[1]
#define J1_13            ddJi[2]

#define J1_21            ddJi[3]
#define J1_22            ddJi[4]
#define J1_23            ddJi[5]

#define J1_31            ddJi[6]
#define J1_32            ddJi[7]
#define J1_33            ddJi[8]

//I2_ij:
#define J2_11            ddJi[9]
#define J2_12            ddJi[10]
#define J2_13            ddJi[11]

#define J2_21            ddJi[12]
#define J2_22            ddJi[13]
#define J2_23            ddJi[14]

#define J2_31            ddJi[15]
#define J2_32            ddJi[16]
#define J2_33            ddJi[17]

//I3_ij:
#define J3_11            ddJi[18]
#define J3_12            ddJi[19]
#define J3_13            ddJi[20]

#define J3_21            ddJi[21]
#define J3_22            ddJi[22]
#define J3_23            ddJi[23]

#define J3_31            ddJi[24]
#define J3_32            ddJi[25]
#define J3_33            ddJi[26]

//I11_ij:
#define J11_11           ddJij[0]
#define J11_12           ddJij[1]
#define J11_13           ddJij[2]

#define J11_21           ddJij[3]
#define J11_22           ddJij[4]
#define J11_23           ddJij[5]

#define J11_31           ddJij[6]
#define J11_32           ddJij[7]
#define J11_33           ddJij[8]

//I12_ij:
#define J12_11           ddJij[9]
#define J12_12           ddJij[10]
#define J12_13           ddJij[11]

#define J12_21           ddJij[12]
#define J12_22           ddJij[13]
#define J12_23           ddJij[14]

#define J12_31           ddJij[15]
#define J12_32           ddJij[16]
#define J12_33           ddJij[17]

//I13_ij:
#define J13_11           ddJij[18]
#define J13_12           ddJij[19]
#define J13_13           ddJij[20]

#define J13_21           ddJij[21]
#define J13_22           ddJij[22]
#define J13_23           ddJij[23]

#define J13_31           ddJij[24]
#define J13_32           ddJij[25]
#define J13_33           ddJij[26]

//I21_ij:
#define J21_11           ddJij[27]
#define J21_12           ddJij[28]
#define J21_13           ddJij[29]

#define J21_21           ddJij[30]
#define J21_22           ddJij[31]
#define J21_23           ddJij[32]

#define J21_31           ddJij[33]
#define J21_32           ddJij[34]
#define J21_33           ddJij[35]

//I22_ij:
#define J22_11           ddJij[36]
#define J22_12           ddJij[37]
#define J22_13           ddJij[38]

#define J22_21           ddJij[39]
#define J22_22           ddJij[40]
#define J22_23           ddJij[41]

#define J22_31           ddJij[42]
#define J22_32           ddJij[43]
#define J22_33           ddJij[44]

//I23_ij:
#define J23_11           ddJij[45]
#define J23_12           ddJij[46]
#define J23_13           ddJij[47]

#define J23_21           ddJij[48]
#define J23_22           ddJij[49]
#define J23_23           ddJij[50]

#define J23_31           ddJij[51]
#define J23_32           ddJij[52]
#define J23_33           ddJij[53]

//I31_ij:
#define J31_11           ddJij[54]
#define J31_12           ddJij[55]
#define J31_13           ddJij[56]

#define J31_21           ddJij[57]
#define J31_22           ddJij[58]
#define J31_23           ddJij[59]

#define J31_31           ddJij[60]
#define J31_32           ddJij[61]
#define J31_33           ddJij[62]

//I32_ij:
#define J32_11           ddJij[63]
#define J32_12           ddJij[64]
#define J32_13           ddJij[65]

#define J32_21           ddJij[66]
#define J32_22           ddJij[67]
#define J32_23           ddJij[68]

#define J32_31           ddJij[69]
#define J32_32           ddJij[70]
#define J32_33           ddJij[71]

//I31_ij:
#define J33_11           ddJij[72]
#define J33_12           ddJij[73]
#define J33_13           ddJij[74]

#define J33_21           ddJij[75]
#define J33_22           ddJij[76]
#define J33_23           ddJij[77]

#define J33_31           ddJij[78]
#define J33_32           ddJij[79]
#define J33_33           ddJij[80]

//Transformation matrix flags according the order of Euller rotations
typedef enum{
  _NONE_, _XZX_, _XYX_, _YXY_, _YZY_, _ZYZ_, _ZXZ_, _XZY_, _XYZ_, _YXZ_, _YZX_, _ZYX_, _ZXY_,
}EullerRotations;

//h - equidistant for complex diferentiation of Ferers-Dyson's integrals
#define _DIFF_H_              1.0e-7 // Used with complex differentiation
//
//#define _NDIFF_H _             1.0e-4 // Used with real differentioation  - toto je puvodni
//#define _NDIFF_H _MULTIPLIER_  0.008  // Used with real differentioation  - toto pridal standa

// Types of derivatives of elliptical integrals
enum DiffTypes { DT_Void, DT_ANALITICAL, DT_NUMERICAL, DT_COMPLEX };


//structure of various lambda-related parameters for complex differentiation
typedef struct{
  std::complex<double> cx1[3];      //Coordinates of a point in complex plane
  std::complex<double> cx2[3];      //Coordinates of a point in complex plane
  std::complex<double> cx3[3];      //Coordinates of a point in complex plane
  std::complex<double> cla1;        //lambda parameter of a point in complex plane (this is not the Lame lambda !!!)
  std::complex<double> cla2;        //lambda parameter of a point in complex plane (this is not the Lame lambda !!!)
  std::complex<double> cla3;        //lambda parameter of a point in complex plane (this is not the Lame lambda !!!)
  std::complex<double> cDla1;       //Delta lambda parameter in complex plane
  std::complex<double> cDla2;       //Delta lambda parameter in complex plane
  std::complex<double> cDla3;       //Delta lambda parameter in complex plane
  std::complex<double> cdla_11;     //complex derivatives of Lambda according _x1_ in projection to _x1_
  std::complex<double> cdla_12;     //complex derivatives of Lambda according _x1_ in projection to _x2_
  std::complex<double> cdla_13;     //complex derivatives of Lambda according _x1_ in projection to _x3_
  std::complex<double> cdla_22;     //complex derivatives of Lambda according _x2_ in projection to _x2_
  std::complex<double> cdla_23;     //complex derivatives of Lambda according _x2_ in projection to _x3_
  std::complex<double> cdla_33;     //complex derivatives of Lambda according _x3_ in projection to _x3_
  double mult;                //mult = -2. * PI * a[0] * a[1] * a[2]
}cLambda;

//structure of various lambda-related parameters for analytic differentiation
typedef struct{
  double la;        //lambda parameter of a point (this is not the Lame lambda !!!)
  double Dla;       //Delta lambda parameter
  double dla[4];    //first derivatives of Lambda according _x1_, _x2_, _x3_
  double ddla[10];  //second derivatives of Lambda
  double mult;      //mult = -2. * PI * a[0] * a[1] * a[2]
}aLambda;

//structure of various lambda-related parameters for numerical differentiation
typedef struct{
  double la;          //lambda parameter of the point [x1,x2,x3] (this is not the Lame lambda !!!)
  double Dla;         //Delta lambda parameter of the point [x1,x2,x3]
  double dla;
  double x1h[3];      //Coordinates of the point [x1+h,x2,x3]
  double x2h[3];      //Coordinates of the point [x1,x2+h,x3]
  double x3h[3];      //Coordinates of the point [x1,x2,x3+h]
  double la1;         //lambda parameter of the point [x1+h,x2,x3] (this is not the Lame lambda !!!)
  double la2;         //lambda parameter of the point [x1,x2+h,x3] (this is not the Lame lambda !!!)
  double la3;         //lambda parameter of the point [x1,x2,x3+h] (this is not the Lame lambda !!!)
  double Dla1;        //Delta lambda parameter of the point [x1+h,x2,x3]
  double Dla2;        //Delta lambda parameter of the point [x1,x2+h,x3]
  double Dla3;        //Delta lambda parameter of the point [x1,x2,x3+h]
  double dla_11;
  double dla_12;
  double dla_13;
  double dla_22;
  double dla_23;
  double dla_33;
/*   double dla[4];      //first derivatives of Lambda according _x1_, _x2_, _x3_ */
/*   double dla_1[4];    //first derivatives of Lambda in projection to _x1_ */
/*   double dla_2[4];    //first derivatives of Lambda in projection to _x2_ */
/*   double dla_3[4];    //first derivatives of Lambda in projection to _x3_ */
  double mult;        //mult = -2. * PI * a[0] * a[1] * a[2]
}nLambda;

// Number of perturbated points used to differentiate the Eliptic integrals.
#define NUM_PERTURB 19 // Number of all points in 3x3x3 3D grid without the points in corners.
// Define all perturbations needed to compute numerical derivatives of 1st and 2nd order in xyz space.
typedef enum{
  _CENTER_ = 0,
  _Xp_ = 1,
  _Xm_ = 2,
  _Yp_ = 3,
  _Ym_ = 4,
  _Zp_ = 5,
  _Zm_ = 6,
  _XpYp_ = 7,
  _XpYm_ = 8,
  _XmYp_ = 9,
  _XmYm_ = 10,
  _XpZp_ = 11,
  _XpZm_ = 12,
  _XmZp_ = 13,
  _XmZm_ = 14,
  _YpZp_ = 15,
  _YpZm_ = 16,
  _YmZp_ = 17,
  _YmZm_ = 18
}Perturb;



} // end of namespace mumech

#endif

/*end of file*/