contactmodelmechanical.h

Go to the documentation of this file.

#pragma once
// contactmodelmechanical.h
#include "contactmodel.h"
#include "contactmodel/interface/icontactmodelmechanical.h"

#pragma warning(push)
#pragma warning(disable:4251)

namespace itascaxd
{
    class IContactMechanical;
    class IContactThermal;
}

namespace cmodelsxd
{
    using namespace itasca;
    using namespace itascaxd;

    class ContactModelThermalState;
    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) { }


        virtual ~ContactModelMechanicalState() { };
        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();
        
        virtual ~ContactModelMechanical();
        
        static const char *  getPluginPrefix() { return "contactmodelmechanical"; }
        
        static const char *  getPluginDirectory() { return "contactmodels/mechanical"; }

        //  IContactModel overrides

        virtual bool    isOKToDelete() const { return true; }

        // IContactModelMechanical functions.

        virtual ContactModel * getContactModel() {return this;}
        
        virtual const ContactModel * getContactModel() const {return this;}
        
        virtual bool    checkActivity(const double &) {return false;}
        
        virtual bool    isSliding() const { return false; }
        
        virtual bool    isBonded() const { return false; }

        virtual void    unbond() { }

        virtual bool    hasNormal() const { return false; };
        
        virtual DVect3  getNormal() const { return DVect3(0.0); }
        
        virtual void    propagateStateInformation(IContactModelMechanical *,const CAxes &n=CAxes(),const CAxes &l=CAxes()) {n; l;}
        
        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 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);}

        virtual void    getSphereList(const IContact *con,std::vector<DVect> *pos,std::vector<double> *rad,std::vector<double> *val) { con; pos; rad; val; }
#ifdef THREED

        virtual void    getDiskList(const IContact *con,std::vector<DVect> *pos,std::vector<DVect> *normal,std::vector<double> *radius,std::vector<double> *val) { con; pos; normal; radius; val; }
#endif

        virtual void    getCylinderList(const IContact *con,std::vector<DVect> *bot,std::vector<DVect> *top,std::vector<double> *radlow,std::vector<double> *radhi,std::vector<double> *val) { con; bot; top; radlow; radhi; val; }

        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); }

        virtual ContactModelMechanical *clone() const { throw Exception("Must be reimplemented."); }

    };
} // namespace cmodelxd

#pragma warning(pop)


// EoF