MomentumBalance.h

Go to the documentation of this file.

#pragma once

#include "nuto/mechanics/constitutive/MechanicsInterface.h"
#include "nuto/mechanics/dofs/DofType.h"
#include "nuto/mechanics/dofs/DofVector.h"
#include "nuto/mechanics/dofs/DofMatrix.h"
#include "nuto/mechanics/interpolation/TypeDefs.h"
#include "nuto/mechanics/cell/CellIpData.h"
#include "nuto/mechanics/cell/CellData.h"

namespace NuTo
{
namespace Integrands
{

template <int TDim>
class MomentumBalance
{
public:
    MomentumBalance(DofType dofType, const Laws::MechanicsInterface<TDim>& law)
        : mDofType(dofType)
        , mLaw(law)
    {
    }

    DofVector<double> Gradient(const CellIpData& cellIpData, double deltaT)
    {
        DofVector<double> gradient;

        Eigen::MatrixXd B = cellIpData.B(mDofType, Nabla::Strain());
        gradient[mDofType] =
                B.transpose() * mLaw.Stress(cellIpData.Apply(mDofType, Nabla::Strain()), deltaT, cellIpData.Ids());

        return gradient;
    }

    DofMatrix<double> Hessian0(const CellIpData& cellIpData, double deltaT)
    {
        DofMatrix<double> hessian0;

        Eigen::MatrixXd B = cellIpData.B(mDofType, Nabla::Strain());
        hessian0(mDofType, mDofType) =
                B.transpose() * mLaw.Tangent(cellIpData.Apply(mDofType, Nabla::Strain()), deltaT, cellIpData.Ids()) * B;

        return hessian0;
    }

protected:
    DofType mDofType;

private:
    const Laws::MechanicsInterface<TDim>& mLaw;
};
} /* Integrand */
} /* NuTo */