From 3da9557a1c4a75cd8ccf1bc5a887bfe61cc49f64 Mon Sep 17 00:00:00 2001
From: jleclerc <julien.leclerc@ulg.ac.be>
Date: Wed, 13 Mar 2019 10:31:39 +0100
Subject: [PATCH] clen impl of transition
---
.../materialLaw/mlawNonLocalPorousCoupled.cpp | 20 +++++++++++++++++++
1 file changed, 20 insertions(+)
diff --git a/NonLinearSolver/materialLaw/mlawNonLocalPorousCoupled.cpp b/NonLinearSolver/materialLaw/mlawNonLocalPorousCoupled.cpp
index de2feaa64..ef3dfaa92 100644
--- a/NonLinearSolver/materialLaw/mlawNonLocalPorousCoupled.cpp
+++ b/NonLinearSolver/materialLaw/mlawNonLocalPorousCoupled.cpp
@@ -577,6 +577,21 @@ void mlawNonLocalPorousCoupledLaw::checkCoalescence(IPNonLocalPorosity* q1, cons
else
{
// In the bulk:
+/*
+ // Only growth model is allowed
+ q1Thom->getRefToCoalescenceOnsetFlag() = false;
+ // Reset onset variables to their default value
+ q1Thom->getRefToPorosityAtCoalescenceOnset() = 0.;
+ q1Thom->getRefToLigamentRatioAtCoalescenceOnset() = 0.;
+ q1Thom->getRefToAspectRatioAtCoalescenceOnset() = 0.;
+ q1Thom->getRefToShapeFactorAtCoalescenceOnset() = 0.;
+
+ q1Thom->getRefToCoalescenceActiveFlag() = false;
+ q1Thom->getRefToCrackOffsetOnCft() = 1.;
+ q1Thom->getRefToAccelerateRate() = 1.;
+ q1Thom->getRefToYieldOffset() = 0.;
+
+*/
if (!q0Thom->getCoalescenceOnsetFlag()){
double yieldThomason = _mlawCoales->yieldFunction(kcorEq,pcor,R,yieldfV,q0,q1,T);
@@ -630,6 +645,11 @@ void mlawNonLocalPorousCoupledLaw::checkCoalescence(IPNonLocalPorosity* q1, cons
q1Thom->getRefToAccelerateRate() = q0Thom->getAccelerateRate();
}
}
+
+
+
+ // In the bulk: end
+
}
else{
// If crack transition is not used: no distinction between bulk and interface ipvs
--
GitLab