GlobalFractureEnergyIntegrator.h

Go to the documentation of this file.

//
// Created by Thomas Titscher on 10/28/16.
//
#pragma once
#include <string>
#include <Eigen/Core>

namespace NuTo
{
namespace Tools
{

class GlobalFractureEnergyIntegrator
{
public:
    GlobalFractureEnergyIntegrator(const Eigen::VectorXd& rForces, const Eigen::VectorXd& rDisplacements);

    GlobalFractureEnergyIntegrator(std::string rFileForces, std::string rFileDisplacements, int rColumn);

    double IntegrateSofteningCurve(double rArea) const;

    double IntegrateSofteningCurve(double rArea, double rForceThreshold) const;

private:
    double IntegrateIndexToForceThreshold(int rIndexEnd, double rForceThreshold) const;

    int FindIndexJustAboveForceThreshold(double rForceThreshold) const;

    double IntegrateSofteningCurveInternal(int rEnd) const;

    void CheckData() const;

    Eigen::VectorXd mForce;
    Eigen::VectorXd mDispl;
};

} // namespace Tools
} // namespace NuTo