Line data    Source code

       1             : #ifndef MODEL_INTERFACE_CRACK_CONSTANT_H
       2             : #define MODEL_INTERFACE_CRACK_CONSTANT_H
       3             : 
       4             : #include <iostream>
       5             : #include <fstream>
       6             : 
       7             : #include "AMReX.H"
       8             : #include "Crack.H"
       9             : #include "Set/Set.H"
      10             : #include "Util/Util.H"
      11             : 
      12             : namespace Model
      13             : {
      14             : namespace Interface
      15             : {
      16             : namespace Crack
      17             : {
      18             : class Constant : public Crack
      19             : {
      20             : public:
      21           0 :     Constant(){};
      22             :     
      23           0 :     Set::Scalar Gc (Set::Scalar /*theta*/) { return _Gc; }
      24           0 :     Set::Scalar DGc (Set::Scalar /*theta*/) { return 0.0; }
      25           0 :     Set::Scalar DDGc (Set::Scalar /*theta*/) { return 0.0; }
      26           0 :     Set::Scalar Zeta(Set::Scalar /*theta*/) { return zeta; }
      27           0 :     Set::Scalar Mobility (Set::Scalar /*theta*/) {return mobility;}
      28           0 :     Set::Scalar DrivingForceThreshold (Set::Scalar /*theta*/) {return threshold;}
      29             : private:
      30             :     Set::Scalar _Gc = 1.0, zeta = 1.e-2, mobility = 1.e-2, threshold = 0.0;
      31             : 
      32             : public:
      33           0 :     static void Parse(Constant & value, IO::ParmParse & pp)
      34             :     {
      35           0 :         pp_query("G_c",value._Gc); // Fracture energy 
      36           0 :         pp_query("zeta",value.zeta); // Lengthscale regularization
      37           0 :         pp_query("mobility",value.mobility); // Mobility (speed)
      38           0 :         pp_query("threshold", value.threshold); // Threshold
      39             :         
      40           0 :         std::string gtype = "";
      41           0 :         std::string wtype = "";
      42           0 :         pp_query("gtype", gtype); // Type of g function to use {square, multiwell, 4c3, squarep, squarepexp, cubicm}
      43           0 :         pp_query("wtype", wtype); // Type o w function to use {square, multiwell, multiwell2, 4c3}
      44             : 
      45           0 :         Set::Scalar ductile_exponent = 1.;
      46           0 :         if(gtype == "squarepexp"|| gtype == "cubicm") pp_query("exponent",ductile_exponent);
      47             : 
      48           0 :         std::map<std::string,Model::Interface::Crack::Crack::GType>  g_map;
      49           0 :         g_map["square"] = Model::Interface::Crack::Crack::GType::GSQUARE;
      50           0 :         g_map["multiwell"] = Model::Interface::Crack::Crack::GType::GMULTIWELL;
      51           0 :         g_map["4c3"] = Model::Interface::Crack::Crack::GType::GPHI4C3;
      52           0 :         g_map["squarep"] = Model::Interface::Crack::Crack::GType::GSQUAREP;
      53           0 :         g_map["squarepexp"] = Model::Interface::Crack::Crack::GType::GSQUAREPM;
      54           0 :         g_map["cubicm"] = Model::Interface::Crack::Crack::GType::GCUBICM;
      55             : 
      56           0 :         std::map<std::string,Model::Interface::Crack::Crack::WType>  w_map;
      57           0 :         w_map["square"] = Model::Interface::Crack::Crack::WType::WSQUARE;
      58           0 :         w_map["multiwell"] = Model::Interface::Crack::Crack::WType::WMULTIWELL;
      59           0 :         w_map["multiwell2"] = Model::Interface::Crack::Crack::WType::WMULTIWELL2;
      60           0 :         w_map["4c3"] = Model::Interface::Crack::Crack::WType::WPHI4C3;
      61             : 
      62           0 :         value.SetGType(g_map[gtype]);
      63           0 :         value.SetWType(w_map[wtype]);
      64             : 
      65           0 :         value.SetDuctileExponent(ductile_exponent);
      66           0 :     }
      67             : };
      68             : }
      69             : }
      70             : }
      71             : #endif