DamageLawExponential.h

Go to the documentation of this file.

#pragma once

#include <cmath>
#include "DamageLaw.h"
#include "SofteningMaterial.h"

namespace NuTo
{
namespace Constitutive
{

class DamageLawExponential : public DamageLaw
{
public:
    DamageLawExponential(double kappa0, double beta, double alpha)
        : mKappa0(kappa0)
        , mBeta(beta)
        , mAlpha(alpha)
    {
    }
    DamageLawExponential(Material::Softening m)
        : mKappa0(m.ft / m.E)
        , mBeta(m.ft / m.gf)
        , mAlpha(1 - m.fMin / m.ft)
    {
    }

    double Damage(double kappa) const override
    {
        if (kappa < mKappa0)
            return 0.;
        return 1 - mKappa0 / kappa * (1 - mAlpha + mAlpha * std::exp(mBeta * (mKappa0 - kappa)));
    }

    double Derivative(double kappa) const override
    {
        if (kappa < mKappa0)
            return 0.;
        return mKappa0 / kappa *
               ((1 / kappa + mBeta) * mAlpha * std::exp(mBeta * (mKappa0 - kappa)) + (1 - mAlpha) / kappa);
    }

private:
    double mKappa0;
    double mBeta;
    double mAlpha;
};
} /* Constitutive */
} /* NuTo */