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