double TrabBoneEmbed :: computeDamage(GaussPoint *gp, TimeStep *atTime) { double tempAlpha; computeCumPlastStrain(tempAlpha, gp, atTime); double tempDam = computeDamageParam(tempAlpha, gp); // double dam=0.0; return tempDam; }
double MisesMat :: computeDamage(GaussPoint *gp, TimeStep *tStep) { double tempKappa, dam; MisesMatStatus *status = static_cast< MisesMatStatus * >( this->giveStatus(gp) ); dam = status->giveDamage(); computeCumPlastStrain(tempKappa, gp, tStep); double tempDam = computeDamageParam(tempKappa); if ( dam > tempDam ) { tempDam = dam; } return tempDam; }