add new func

NonLinearSolver/materialLaw/mlawNonLocalPorous.h

@@ -110,6 +110,7 @@ class mlawNonLocalPorosity: public materialLaw
     virtual void setScatterredInitialPorosity(double f0min, double f0max);
     virtual void setCoalescenceLaw(const CoalescenceLaw& added_coalsLaw);
     virtual void setShearPorosityGrowthFactor(const double k);
+    virtual void setStressTriaxialityFunction_kw(const scalarFunction& fct);
     virtual void setLocalRegularizedFunction(const scalarFunction& fct);
     virtual void setCorrectedRegularizedFunction(const scalarFunction& fct);
     virtual void setPostBlockedDissipationBehavior(const int method);
@@ -139,7 +140,7 @@ class mlawNonLocalPorosity: public materialLaw
     J2IsotropicHardening* &getRefToJ2IsotropicHardening() {
     return _j2IH;
     } ;
-    virtual void setStressTriaxialityFunction_kw(const scalarFunction& fct);
+
 //---------------------end
   #ifndef SWIG
     virtual mlawNonLocalPorosity::nonLocalMethod getNonLocalMethod() const;

NonLinearSolver/materialLaw/mlawNonLocalPorousCoupled.cpp

@@ -205,6 +205,12 @@ void mlawNonLocalPorousCoupledLaw::setTemperatureFunction_alp(const scalarFuncti
   _mlawCoales->setTemperatureFunction_alp(Tfunc);
 };
 
+void mlawNonLocalPorousCoupledLaw::setStressTriaxialityFunction_kw(const scalarFunction& fct){
+  mlawNonLocalPorosity::setStressTriaxialityFunction_kw(fct);
+  _mlawGrowth->setStressTriaxialityFunction_kw(fct);
+  _mlawCoales->setStressTriaxialityFunction_kw(fct);
+};
+
 void mlawNonLocalPorousCoupledLaw::createIPState(IPStateBase* &ips,const bool* state_,const MElement *ele, const int nbFF_, const IntPt *GP, const int gpt) const
 {
   //bool inter=true;

NonLinearSolver/materialLaw/mlawNonLocalPorousCoupled.h

@@ -69,6 +69,8 @@ class mlawNonLocalPorousCoupledLaw : public mlawNonLocalPorosity
   virtual void setTemperatureFunction_nu(const scalarFunction& Tfunc);
   virtual void setTemperatureFunction_alp(const scalarFunction& Tfunc);
 
+  virtual void setStressTriaxialityFunction_kw(const scalarFunction& fct);
+
   // Constructors & destructor
   mlawNonLocalPorousCoupledLaw(const mlawNonLocalPorousCoupledLaw &source);
   virtual ~mlawNonLocalPorousCoupledLaw()

dG3D/src/nonLocalDamageDG3DMaterialLaw.cpp

@@ -853,6 +853,10 @@ void NonLocalDamageGursonDG3DMaterialLaw::setShearPorosityGrowthFactor(const dou
   _nldGursonlaw->setShearPorosityGrowthFactor(k);
 };
 
+void NonLocalDamageGursonDG3DMaterialLaw::setStressTriaxialityFunction_kw(const scalarFunction& fct){
+  _nldGursonlaw->setStressTriaxialityFunction_kw(fct);
+};
+
 void NonLocalDamageGursonDG3DMaterialLaw::clearCLengthLaw(){
   _nldGursonlaw->clearCLengthLaw();
 };
@@ -1138,6 +1142,9 @@ void NonLocalPorousThomasonDG3DMaterialLaw::setOnsetTriaxialityForRoundedYieldSu
 void NonLocalPorousThomasonDG3DMaterialLaw::setShearPorosityGrowthFactor(const double k){
   _nlpthomlaw-> setShearPorosityGrowthFactor(k);
 };
+void NonLocalPorousThomasonDG3DMaterialLaw::setStressTriaxialityFunction_kw(const scalarFunction& fct){
+  _nlpthomlaw->setStressTriaxialityFunction_kw(fct);
+};
 
 void NonLocalPorousThomasonDG3DMaterialLaw::setSubStepping(const bool fl, const int maxNumStep){
   _nlpthomlaw->setSubStepping(fl,maxNumStep);
@@ -1366,6 +1373,10 @@ void NonLocalPorousCoupledDG3DMaterialLaw::setShearPorosityGrowthFactor(const do
   _nlpcoupledlaw->setShearPorosityGrowthFactor(k);
 };
 
+void NonLocalPorousCoupledDG3DMaterialLaw::setStressTriaxialityFunction_kw(const scalarFunction& fct){
+  _nlpcoupledlaw->setStressTriaxialityFunction_kw(fct);
+};
+
 void NonLocalPorousCoupledDG3DMaterialLaw::setCrackTransition(const bool fl){
   _nlpcoupledlaw->setCrackTransition(fl);
 }

dG3D/src/nonLocalDamageDG3DMaterialLaw.h

@@ -326,6 +326,7 @@ public:
     void setLocalRegularizedFunction(const scalarFunction& fct);
     void setNonLocalMethod(const int i);
     void setShearPorosityGrowthFactor(const double k);
+    void setStressTriaxialityFunction_kw(const scalarFunction& fct);
     void clearCLengthLaw();
     void setCLengthLaw(const CLengthLaw& clength);
 
@@ -395,6 +396,7 @@ class NonLocalPorousThomasonDG3DMaterialLaw : public dG3DMaterialLaw
     void setSubStepping(const bool fl, const int maxNumStep);
     void setPostBlockedDissipationBehavior(const int method);
     void setShearPorosityGrowthFactor(const double k);
+    void setStressTriaxialityFunction_kw(const scalarFunction& fct);
     void setBalanceEquationType(const int type);
     void clearCLengthLaw();
     void setCLengthLaw(const CLengthLaw& clength);
@@ -462,6 +464,7 @@ protected:
     void setYieldSurfaceExponent(const double newN);
     void setSubStepping(const bool fl, const int maxNumStep);
     void setShearPorosityGrowthFactor(const double k);
+    void setStressTriaxialityFunction_kw(const scalarFunction& fct);
     void setCrackTransition(const bool fl);
     void setPostBlockedDissipationBehavior(const int method);
     void setBalanceEquationType(const int type);