contact_spr_error_process.h

Go to the documentation of this file.

// KRATOS    ______            __             __  _____ __                  __                   __
//          / ____/___  ____  / /_____ ______/ /_/ ___// /________  _______/ /___  ___________ _/ /
//         / /   / __ \/ __ \/ __/ __ `/ ___/ __/\__ \/ __/ ___/ / / / ___/ __/ / / / ___/ __ `/ / 
//        / /___/ /_/ / / / / /_/ /_/ / /__/ /_ ___/ / /_/ /  / /_/ / /__/ /_/ /_/ / /  / /_/ / /  
//        \____/\____/_/ /_/\__/\__,_/\___/\__//____/\__/_/   \__,_/\___/\__/\__,_/_/   \__,_/_/  MECHANICS
//
//  License:         BSD License
//                   license: ContactStructuralMechanicsApplication/license.txt
//
//  Main authors:    Anna Rehr
//  Co-author   :    Vicente Mataix Ferrandiz
//
 
#pragma once
 
// System includes
 
// External includes
 
// Project includes
#include "custom_processes/spr_error_process.h"
 
namespace Kratos
{
 
 
    using SizeType = std::size_t;
 
 
 
 
template<SizeType TDim>
class KRATOS_API(CONTACT_STRUCTURAL_MECHANICS_APPLICATION) ContactSPRErrorProcess
    : public SPRErrorProcess<TDim>
{
public:
 
 
    using BaseType = SPRErrorProcess<TDim>;
 
    using NodesArrayType = ModelPart::NodesContainerType;
    using ElementsArrayType = ModelPart::ElementsContainerType;
    using ConditionsArrayType = ModelPart::ConditionsContainerType;
 
    using WeakElementItType = GlobalPointersVector<Element>::iterator;
    using NodeItType = NodesArrayType::iterator;
    using ElementItType = ElementsArrayType::iterator;
 
    using IndexType = std::size_t;
 
    KRATOS_CLASS_POINTER_DEFINITION(ContactSPRErrorProcess);
 
    static constexpr SizeType SigmaSize = (TDim == 2) ? 3 : 6;
 
 
    ContactSPRErrorProcess(
        ModelPart& rThisModelPart,
        Parameters ThisParameters = Parameters(R"({})")
        );
 
    ~ContactSPRErrorProcess() override {}
 
 
    void operator()()
    {
        this->Execute();
    }
 
    const Parameters GetDefaultParameters() const override;
 
 
 
 
 
 
 
    std::string Info() const override
    {
        return "ContactSPRErrorProcess";
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "ContactSPRErrorProcess";
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
    }
 
protected:
 
 
 
 
 
 
 
    void CalculatePatch(
        NodeItType itNode,
        NodeItType itPatchNode,
        const SizeType NeighbourSize,
        Vector& rSigmaRecovered
        ) override;
 
 
 
 
 
 
 
 
private:
 
 
    double mPenaltyNormal;   
    double mPenaltyTangent;  
 
 
 
    void ComputeNormalTangentMatrices(
        BoundedMatrix<double, 1, SigmaSize>& rNk,
        BoundedMatrix<double, 1, SigmaSize>& rTk1,
        BoundedMatrix<double, 1, SigmaSize>& rTk2,
        const array_1d<double, 3>& rNormal
        );
 
 
 
 
 
    ContactSPRErrorProcess& operator=(ContactSPRErrorProcess const& rOther)
    {
        return *this;
    };
 
    //ContactSPRErrorProcess(ContactSPRErrorProcess const& rOther);
 
};// class ContactSPRErrorProcess
 
};// namespace Kratos.