void
RCM2Material :: updateCrackStatus(GaussPoint *gp, const FloatArray &crackStrain)
//
// updates gp records and MatStatus due to cracking.
//
// Updates MatStatus (respective it's temporary variables) to current
// reached status during integrating incremental constitutive relations
// Temporary variables are used, because we may integrate constitutive
// realtions many times for different strainIncrement in order to
// reach equilibrium state, so we don't want to change other variables
// which describing previously reached equlibrium. After a new equilibrium
// is reached, this->updateYourself is called which invokes matStatus->
// updateYourself(), which copies temporary variables to variables
// describing equilibrium.
//
//
{
int i;
double minCrackStrainsForFullyOpenCrack;
RCM2MaterialStatus *status = ( RCM2MaterialStatus * ) this->giveStatus(gp);
IntArray crackMap;
status->giveCrackMap(crackMap);
// check if material previously cracked
// and compute possible crack planes
// or if newer cracked, so we compute principal stresses
// (for iso mat. coincide with princ strains)
// and compare them with reduced tension strength
// as a criterion for crack initiation
// principalStrain = status->givePrincipalStrainVector();
// transform stresses to principal directions of strains
//
// local stress is updated according to reached local crack strain
//
for ( i = 1; i <= 3; i++ ) { // loop over each possible crack plane
if ( ( crackMap.at(i) != 0 ) &&
( this->giveStressStrainComponentIndOf(FullForm, gp->giveMaterialMode(), i) ) ) {
if ( status->giveTempMaxCrackStrain(i) < crackStrain.at(i) ) {
status->setTempMaxCrackStrain( i, crackStrain.at(i) );
}
minCrackStrainsForFullyOpenCrack = this->giveMinCrackStrainsForFullyOpenCrack(gp, i);
//if ((crackStrain.at(i) >= status->giveMinCrackStrainsForFullyOpenCrack(i)) &&
if ( ( crackStrain.at(i) >= minCrackStrainsForFullyOpenCrack ) &&
( crackStrain.at(i) >= status->giveTempMaxCrackStrain(i) ) ) {
//
// fully open crack
//
status->setTempCrackStatus(i, pscm_OPEN);
} else if ( crackStrain.at(i) >= status->giveTempMaxCrackStrain(i) ) {
//
// further softening of crack
//
status->setTempCrackStatus(i, pscm_SOFTENING);
// status->giveTempReachedSofteningStress()->at(i) = localStress->at(i);
} else if ( crackStrain.at(i) <= 0. ) {
if ( status->giveTempCrackStatus(i) != pscm_NONE ) {
// previously active crack becomes closed
status->setTempCrackStatus(i, pscm_CLOSED);
}
} else {
//
// crack unloading or reloading
//
status->setTempCrackStatus(i, pscm_UNLOADING);
}
} // end for possible direction
} // end loop over prin directions
}
