Skip to content
Snippets Groups Projects
Commit 7e8a2100 authored by Ludovic Noels's avatar Ludovic Noels
Browse files

modicy local variabel locatin

parent 0cac8ac5
No related branches found
No related tags found
No related merge requests found
......@@ -229,8 +229,6 @@ void mlawNonLocalDamageIsotropicElasticity::constitutive(
}
// modify ipv value
q1->getRefToLocalEffectiveStrains() = sqrt(equivalentStrainsSquared);
localVariable(0)=q1->getRefToLocalEffectiveStrains();
// - Elastic energy release rate
double effEner = deformationEnergy(epsilon_strains,Peff);
......@@ -308,6 +306,7 @@ void mlawNonLocalDamageIsotropicElasticity::constitutive(
}
}
}
localVariable(0)=q1->getRefToLocalEffectiveStrains();
......
......@@ -321,7 +321,6 @@ void mlawNonLocalDamageJ2FullyCoupledThermoMechanics::predictorCorrector(
const STensor3 &Fp = ipvcur->getConstRefToFp();
STensorOperation::inverseSTensor3(Fp,Fpinv);
STensorOperation::multSTensor3(Fn,Fpinv,Fe);
localVariable(0) = ipvcur->getConstRefToEquivalentPlasticStrain();
double ene = ipvcur->defoEnergy();
......@@ -409,6 +408,7 @@ void mlawNonLocalDamageJ2FullyCoupledThermoMechanics::predictorCorrector(
else{
irrEneg = 0.;
}
localVariable(0) = ipvcur->getConstRefToEquivalentPlasticStrain();
if(stiff)
......
......@@ -164,7 +164,6 @@ void mlawNonLocalDamageJ2Hyper::constitutive(
const STensor3 &Fp = ipvcur->getConstRefToFp();
STensorOperation::inverseSTensor3(Fp,Fpinv);
STensorOperation::multSTensor3(Fn,Fpinv,Fe);
localVariable(0) = ipvcur->getConstRefToEquivalentPlasticStrain();
double ene = ipvcur->defoEnergy();
......@@ -248,6 +247,7 @@ void mlawNonLocalDamageJ2Hyper::constitutive(
else{
irrEneg = 0.;
}
localVariable(0) = ipvcur->getConstRefToEquivalentPlasticStrain();
if(stiff)
......
......@@ -137,7 +137,6 @@ void mlawNonLocalDamageJ2SmallStrain::constitutive(
const STensor3& epsPl = ipvcur->_j2lepsp;
const STensor3& epsEl = ipvcur->_Ee;
localVariable(0) = ipvcur->getConstRefToEquivalentPlasticStrain();
double ene = ipvcur->defoEnergy();
if (ipvcur->dissipationIsBlocked()){
......@@ -226,6 +225,7 @@ void mlawNonLocalDamageJ2SmallStrain::constitutive(
{
irrEneg = 0.;
}
localVariable(0) = ipvcur->getConstRefToEquivalentPlasticStrain();
if(stiff)
{
......
......@@ -306,7 +306,7 @@ check.equal(0.000000e+00,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSm
#check.equal(3.399153e-04,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 4),1.e-5)
check.equal(-3.484000e-03,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfInd, 4),1.e-2)
check.equal(3.144224e-03,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfVacOut, 4),1.e-5)
check.equal(-6.847045e-04,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 5),1.e-5)
check.equal(-6.847148e-04,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfSmp, 5),1.e-5)
check.equal(3.720132e-01,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfInd, 5),1.e-5)
check.equal(5.949466e+03,mysolver.getArchivedForceOnPhysicalGroup("Face", SurfVacOut, 5),1.e-5)
check.equal(0.0,mysolver.getArchivedNodalValue(SmpRefPoint,0,mysolver.displacement),1.e-5)
......
SetFactory("OpenCASCADE");
Mesh.Optimize = 1;
cm = 1e-2; // Unit
pp = "Input/10Geometric dimensions/0";
close_menu = 1;
colorpp = "Ivory";
// Smp
RintSMP = 0.95*cm;
RextSMP = 1*cm;
H = 2*cm;
DefineConstant[
RintSMP = {0.95*cm, Name StrCat[pp, "4Sample inner radius [m]"], Highlight Str[colorpp]}
RextSMP = {1*cm, Name StrCat[pp, "5Sample outer radius [m]"], Highlight Str[colorpp]}
H = {2*cm, Name StrCat[pp, "6Sample height [m]"], Highlight Str[colorpp]}
];
ECORE = 1000;
SKINECORE = 11110;
REFCORE = 11111;
REFCOREBOTTOM = 11115;
REFCOREPOINT = 11112;
REFCOREOUT = 11113;
REFCOREIN = 11114;
// Smp
vE()+=newv; Cylinder(newv) = {0,-H/2,0, 0,H,0, RintSMP};
vE()+=newv; Cylinder(newv) = {0,-H/2,0, 0,H,0, RextSMP};
vCoreE()+=newv; BooleanDifference(newv) = { Volume{vE(1)}; Delete; }{ Volume{vE(0)}; Delete; };
lc2 = RextSMP/4; // smp
Characteristic Length { PointsOf{ Volume{vCoreE()}; } } = lc2;
Recursive Color SteelBlue {Volume{vCoreE()};}
Physical Volume(ECORE) = vCoreE();
bnd_vCoreE() = CombinedBoundary{Volume{vCoreE()};};
Physical Surface(SKINECORE) = bnd_vCoreE();
Physical Surface(REFCORE)= bnd_vCoreE(1);
Physical Surface(REFCOREBOTTOM)= bnd_vCoreE(2);
Physical Surface(REFCOREOUT)= bnd_vCoreE(0);
Physical Surface(REFCOREIN)= bnd_vCoreE(3);
pts() = PointsOf {Surface { bnd_vCoreE(1) }; };
Physical Point(REFCOREPOINT)= pts(0);
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment