inline double
M4Material :: FN(double en, double sv, double k1, double c1, double c2, double c3, double c4,
double E, double Ev)
/*normal boundary */
{
return ( E * k1 * c1 * exp( -macbra(en - c1 * c2 * k1) / ( k1 * c3 + macbra( -c4 * ( sv / Ev ) ) ) ) );
}
inline double
M4Material :: FT(double sn, double ev, double k1, double k2, double c10,
double c11, double c12, double Et)
/*shear boundary */
{
double a, sn0;
//oprava sn0 fabs<=>macbra
sn0 = Et * k1 * c11 / ( 1 + c12 * macbra(ev) );
a = macbra(-sn + sn0);
return ( Et * k1 * k2 * c10 * a / ( Et * k1 * k2 + c10 * a ) );
}
inline double
M4Material :: FDplus(double ed, double k1, double c5, double c6, double c7, double c20,
double E)
/*positive deviatoric bondary */
{
double a;
a = ( macbra(ed - c5 * c6 * k1) / ( k1 * c20 * c7 ) );
return ( E * k1 * c5 / ( 1 + a * a ) );
}
inline double
M4Material :: FDminus(double ed, double k1, double c7, double c8, double c9,
double E)
/*negative deviatoric boundary */
{
double a;
a = macbra(-ed - c8 * c9 * k1) / ( k1 * c7 );
return ( -E * k1 * c8 / ( 1 + a * a ) );
}
double
RCSDMaterial :: computeCurrEquivStrain(GaussPoint *gp, const FloatArray &reducedTotalStrainVector, double e, TimeStep *tStep)
{
FloatArray effStress, princEffStress, fullEffStress;
FloatMatrix De;
double answer = 0.0;
linearElasticMaterial->giveStiffnessMatrix(De, TangentStiffness, gp, tStep);
effStress.beProductOf(De, reducedTotalStrainVector);
StructuralMaterial :: giveFullSymVectorForm( fullEffStress, effStress, gp->giveMaterialMode() );
this->computePrincipalValues(princEffStress, fullEffStress, principal_stress);
for ( int i = 1; i <= 3; i++ ) {
answer = max( answer, macbra( princEffStress.at(i) ) );
}
return answer / e;
}
void
IsoInterfaceDamageMaterial :: computeEquivalentStrain(double &kappa, const FloatArray &strain, GaussPoint *gp, TimeStep *atTime)
{
kappa = macbra(strain.at(1));
}
inline double
M4Material :: FVplus(double ev, double k1, double c13, double c14, double c15, double Ev)
/*positive volumetric boundary */
{
return ( Ev * k1 * c13 / ( 1 + ( c14 / k1 ) * macbra(ev - c13 * c15 * k1) ) );
}
