QuasistaticSolver.h

Go to the documentation of this file.

#pragma once

#include "nuto/mechanics/solver/Solve.h"
#include "nuto/mechanics/tools/TimeDependentProblem.h"
#include "nuto/mechanics/constraints/Constraints.h"
#include <iosfwd>

namespace NuTo
{
class QuasistaticSolver
{
public:
    QuasistaticSolver(TimeDependentProblem& equations, DofType dof);

    QuasistaticSolver(TimeDependentProblem& equations, std::vector<DofType> dofs);

    void SetConstraints(Constraint::Constraints constraints);

    void SetGlobalTime(double globalTime);

    DofVector<double> TrialState(double newGlobalTime, const NuTo::ConstrainedSystemSolver& solver);

    void UpdateHistory(const DofVector<double>& x);

    DofVector<double> Residual(const DofVector<double>& u);

    DofMatrixSparse<double> Derivative(const DofVector<double>& u);


    double Norm(const DofVector<double>& residual) const;


    void Info(int i, const DofVector<double>& x, const DofVector<double>& r) const;


    double mTolerance = 1.e-10;


    int DoStep(double newGlobalTime, std::string solverType = "EigenSparseLU");

    void WriteTimeDofResidual(std::ostream& out, DofType dofType, std::vector<int> dofNumbers);

private:
    DofVector<double> mX;

    TimeDependentProblem& mProblem;
    Constraint::Constraints mConstraints;

    std::vector<DofType> mDofs;
    DofMatrixSparse<double> mCmatUnit;

    double mGlobalTime = 0;
    double mTimeStep = 0;
};

} /* NuTo */