contactmodelmechanical.h

Go to the documentation of this file.

#pragma once
// contactmodelmechanical.h
#include "contactmodel.h"
#include "contactmodel/interface/icontactmodelmechanical.h"
 
#ifndef __GNUC__
#pragma warning(push)
#pragma warning(disable : 4251)
#endif
 
namespace itascaxd {
    class IContactMechanical;
    class IContactThermal;
    class IContactFluid;
} // namespace itascaxd
 
namespace cmodelsxd {
    using namespace itasca;
    using namespace itascaxd;
 
    class ContactModelThermalState;
    class ContactModelFluidState;
    class ContactModelMechanicalState : public ContactModelState {
      public:
        ContactModelMechanicalState()
            : // force_(0.0)
              //, momentOn1_(0.0)
              //, momentOn2_(0.0)
              /*,*/ relativeTranslationalIncrement_(0.0), relativeAngularIncrement_(0.0),
              end1Curvature_(0.0), end2Curvature_(0.0), inertialMass_(0.0), gap_(0.0),
              canFail_(true) {}
 
        ~ContactModelMechanicalState() override{};
        virtual const IContactMechanical* getMechanicalContact() const = 0;
        virtual const IContact* getContact() const = 0;
 
        // DVect       force_;                            ///< current total force  at this contact
        // (computed by the contact model) DAVect      momentOn1_;                        ///<
        // moment applied to end1 DAVect      momentOn2_;                        ///< moment applied
        // to end2
        DVect relativeTranslationalIncrement_; 
        DAVect relativeAngularIncrement_;      
        DVect2 end1Curvature_;                 
        DVect2 end2Curvature_;                 
        double inertialMass_; 
        double gap_;          
        bool canFail_;        
    };
 
    class CONTACTMODEL_EXPORT ContactModelMechanical : public ContactModel,
                                                       public IContactModelMechanical {
      public:
        ContactModelMechanical();
 
        ~ContactModelMechanical() override;
 
        static const char* getPluginPrefix() { return "contactmodelmechanical"; }
 
        static const char* getPluginDirectory() { return "contactmodels/mechanical"; }
 
        //  IContactModel overrides
 
        bool isOKToDelete() const override { return true; }
 
        // IContactModelMechanical functions.
 
        ContactModel* getContactModel() override { return this; }
 
        const ContactModel* getContactModel() const override { return this; }
 
        virtual bool checkActivity(const double&) { return false; }
 
        bool isSliding() const override { return false; }
 
        bool isBonded() const override { return false; }
 
        void unbond() override {}
 
        bool hasNormal() const override { return false; };
 
        DVect3 getNormal() const override { return DVect3(0.0); }
 
        void setPorePressure(const double&) override {
            throw Exception("The contact model does not support pore pressures.");
        }
 
        void propagateStateInformation(IContactModelMechanical*,
            [[maybe_unused]] const CAxes& n = CAxes(),
            [[maybe_unused]] const CAxes& l = CAxes()) override {}
 
        virtual bool endPropertyUpdated(const QString& name, const IContactMechanical* c) = 0;
 
        virtual bool validate(ContactModelMechanicalState* state, const double& timestep) = 0;
 
        virtual bool forceDisplacementLaw(
            ContactModelMechanicalState* state, const double& timestep) = 0;
 
        virtual bool thermalCoupling(ContactModelMechanicalState*, ContactModelThermalState*,
            IContactThermal*, const double&) {
            return false;
        }
        virtual bool coupling(ContactModelMechanicalState*, ContactModelThermalState*,
            ContactModelFluidState*, const double&) {
            return false;
        }
 
        virtual DVect2 getTimestepLimits(ContactModelMechanicalState*, const double&) const {
            return DVect2(0.0, limits<double>::max());
        }
 
        virtual DVect2 getEffectiveTranslationalStiffness() const { return DVect2(0.0, 0.0); }
 
        virtual DAVect getEffectiveRotationalStiffness() const { return DAVect(0.0); }
 
        void getSphereList(const IContact*, std::vector<DVect>*, std::vector<double>*,
            std::vector<double>*) override {}
#ifdef THREED
 
        void getDiskList(const IContact*, std::vector<DVect>*, std::vector<DVect>*,
            std::vector<double>*, std::vector<double>*) override {}
#endif
 
        void getCylinderList(const IContact*, std::vector<DVect>*, std::vector<DVect>*,
            std::vector<double>*, std::vector<double>*, std::vector<double>*) override {}
 
        virtual DVect getForce(const IContactMechanical*) const { return DVect(0.0); }
 
        virtual DAVect getMomentOn1(const IContactMechanical*) const { return DAVect(0.0); }
 
        virtual DAVect getMomentOn2(const IContactMechanical* c) const { return getMomentOn1(c); }
 
        ContactModelMechanical* clone() const override {
            throw Exception("Must be reimplemented.");
        };
    };
} // namespace cmodelsxd
 
#ifndef __GNUC__
#pragma warning(pop)
#endif
 
// EoF