DamageLawLinear.h

Go to the documentation of this file.

#pragma once

#include <algorithm> // std::min
#include "DamageLaw.h"
#include "SofteningMaterial.h"

namespace NuTo
{
namespace Constitutive
{

class DamageLawLinear : public DamageLaw
{
public:
    DamageLawLinear(double kappa0, double kappaC, double omegaMax)
        : mKappa0(kappa0)
        , mKappaC(kappaC)
        , mOmegaMax(omegaMax)
    {
    }

    DamageLawLinear(Material::Softening m)
        : mKappa0(m.ft / m.E)
        , mKappaC(mKappa0 + 2. * m.gf / m.ft)
        , mOmegaMax(1. - m.fMin / m.ft) // this is not totally right...
    {
    }

    double Damage(double kappa) const override
    {
        if (kappa < mKappa0)
            return 0.;
        double damage = mKappaC / kappa * (kappa - mKappa0) / (mKappaC - mKappa0);
        return std::min(damage, mOmegaMax);
    }

    double Derivative(double kappa) const override
    {
        if (kappa < mKappa0)
            return 0.;
        if (Damage(kappa) < mOmegaMax)
            return mKappaC * mKappa0 / (kappa * kappa * (mKappaC - mKappa0));
        return 0;
    }

private:
    double mKappa0;
    double mKappaC;
    double mOmegaMax;
};
} /* Constitutive */
} /* NuTo */