Skip to content
Snippets Groups Projects
Commit 487c54a7 authored by Julien Leclerc's avatar Julien Leclerc
Browse files

working version of cft offset method 2

parent b9186830
Branches
Tags
1 merge request!128Coal jl8
...@@ -601,6 +601,9 @@ void mlawNonLocalPorousCoupledLaw::checkFailed(IPNonLocalPorosity* q1, const IPN ...@@ -601,6 +601,9 @@ void mlawNonLocalPorousCoupledLaw::checkFailed(IPNonLocalPorosity* q1, const IPN
// If coalescence plastic mode is currently used: criterion based on Cft // If coalescence plastic mode is currently used: criterion based on Cft
double Cft = q1Thom->getConcentrationFactor(); double Cft = q1Thom->getConcentrationFactor();
double CftOffest = q1Thom->getCrackOffsetOnCft(); double CftOffest = q1Thom->getCrackOffsetOnCft();
if (_CfTOffsetMethod==2){
if (!q0->getConstRefToIPCoalescence().getCoalescenceOnsetFlag()){CftOffest = 1.;};
};
if (Cft*CftOffest < (1.-_failedTol)*5.){ if (Cft*CftOffest < (1.-_failedTol)*5.){
q1->setFailed(true); q1->setFailed(true);
Msg::Info("mlawNonLocalPorousCoupledLaw::checkFailed: an ip has failed at Cft = %e", Cft); Msg::Info("mlawNonLocalPorousCoupledLaw::checkFailed: an ip has failed at Cft = %e", Cft);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment