DamageLawNoSoftening.h

Go to the documentation of this file.

#pragma once

#include "DamageLaw.h"
#include "SofteningMaterial.h"

namespace NuTo
{
namespace Constitutive
{


class DamageLawNoSoftening : public DamageLaw
{
public:
    DamageLawNoSoftening(double kappa0)
        : mKappa0(kappa0)
    {
    }
    DamageLawNoSoftening(Material::Softening m)
        : mKappa0(m.ft / m.E)
    {
    }

    double Damage(double kappa) const override
    {
        if (kappa < mKappa0)
            return 0.;
        return 1. - mKappa0 / kappa;
    }

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

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