void
TrabBoneNL3D :: giveRealStressVector_3d(FloatArray &answer, GaussPoint *gp,
const FloatArray &totalStrain, TimeStep *tStep)
{
TrabBoneNL3DStatus *nlStatus = static_cast< TrabBoneNL3DStatus * >( this->giveStatus(gp) );
this->initTempStatus(gp);
double tempDam;
FloatArray effStress, totalStress, densStress;
performPlasticityReturn(gp, totalStrain, tStep);
tempDam = computeDamage(gp, tStep);
effStress = nlStatus->giveTempEffectiveStress();
totalStress = ( 1 - tempDam ) * effStress;
for ( int i = 1; i <= 6; i++ ) {
if ( sqrt( totalStress.at(i) * totalStress.at(i) ) < 1e-8 ) {
totalStress.at(i) = 0.;
}
}
computePlasStrainEnerDensity(gp, totalStrain, totalStress);
if ( densCrit != 0. ) {
computeDensificationStress(densStress, gp, totalStrain, tStep);
answer.add(densStress);
}
answer = totalStress;
nlStatus->setTempDam(tempDam);
nlStatus->letTempStrainVectorBe(totalStrain);
nlStatus->letTempStressVectorBe(answer);
}
void
TrabBoneEmbed :: giveRealStressVector(FloatArray &answer, MatResponseForm form, GaussPoint *gp,
const FloatArray &totalStrain,
TimeStep *atTime)
{
double tempDam, tempTSED;
FloatArray newTotalDef, plasDef;
FloatArray totalStress;
FloatMatrix compliance, elasticity;
this->constructIsoComplTensor(compliance, eps0, nu0);
elasticity.beInverseOf(compliance);
TrabBoneEmbedStatus *status = ( TrabBoneEmbedStatus * ) this->giveStatus(gp);
this->initGpForNewStep(gp);
performPlasticityReturn(gp, totalStrain);
tempDam = computeDamage(gp, atTime);
plasDef.resize(6);
totalStress.beProductOf(elasticity, totalStrain);
tempTSED = 0.5 * totalStrain.dotProduct(totalStress);
answer.resize(6);
answer = totalStress;
status->setTempDam(tempDam);
status->letTempStrainVectorBe(totalStrain);
status->letTempStressVectorBe(answer);
status->setTempTSED(tempTSED);
}
void
MisesMatNl :: giveRealStressVector(FloatArray &answer, GaussPoint *gp,
const FloatArray &totalStrain, TimeStep *tStep)
{
MisesMatNlStatus *nlStatus = static_cast< MisesMatNlStatus * >( this->giveStatus(gp) );
this->initGpForNewStep(gp);
double tempDam;
performPlasticityReturn(gp, totalStrain);
tempDam = this->computeDamage(gp, tStep);
answer.beScaled(1.0 - tempDam, nlStatus->giveTempEffectiveStress());
nlStatus->setTempDamage(tempDam);
nlStatus->letTempStrainVectorBe(totalStrain);
nlStatus->letTempStressVectorBe(answer);
}
void
MisesMatNl :: giveRealStressVector(FloatArray &answer, MatResponseForm form, GaussPoint *gp,
const FloatArray &totalStrain, TimeStep *atTime)
{
MisesMatNlStatus *nlStatus = ( MisesMatNlStatus * ) this->giveStatus(gp);
this->initGpForNewStep(gp);
double tempDam;
FloatArray tempEffStress, totalStress;
MaterialMode mode = gp->giveMaterialMode();
performPlasticityReturn(gp, totalStrain, mode);
tempDam = this->computeDamage(gp, atTime);
nlStatus->giveTempEffectiveStress(tempEffStress);
answer.beScaled( 1.0 - tempDam, tempEffStress);
nlStatus->setTempDamage(tempDam);
nlStatus->letTempStrainVectorBe(totalStrain);
nlStatus->letTempStressVectorBe(answer);
}
