FlukeAnchor.h

#ifndef CFlukeAnchor_H
#define CFlukeAnchor_H
 
#include <fhsim/simobject/ISimObjectCreator.h>
 
#include <vector>
 
namespace rigidbody
{
class CRigidCompositeBody;
namespace soilforces
{
class SoilElement;
}
}; // namespace rigidbody
namespace environment
{
class EnvironmentProvider;
}
 
class FlukeAnchor : public SimObject
{
   public:
   FlukeAnchor(const std::string& simObjectName, ISimObjectCreator* const creator);
   ~FlukeAnchor();
   void OdeFcn(const double T, const double* const X, double* const XDot) const;
 
   void FinalSetup(const double T, const double* const X, ISimObjectCreator* const creator);
 
#ifdef FH_VISUALIZATION
   void RenderInit(Ogre::Root* const ogreRoot, ISimObjectCreator* const creator);
   void RenderUpdate(const double T, const double* const X);
#endif
 
   const double* Position(const double T, const double* const X);
   const double* Velocity(const double T, const double* const X);
   const double* Quater(const double T, const double* const X);
   const double* Omega(const double T, const double* const X);
 
   // Jacobian support (drag and added-mass derivatives omitted; none while the anchor is in contact with the seabed)
   bool HasJacobians() const override;
   void OdeJacobian(double T, const double* X, double* J, int nStates) override;
   int GetJacobianSparsity(int nStates, int* rowPtr, int* colIdx) override;
 
   protected:
   rigidbody::CRigidCompositeBody* m_rigidBody;
 
 
   double m_MaterialRho; // Total mass of anchor [kg]
   double m_ShankLength; // Total volume displacement of anchor [m^3]
   double m_ShankDiameter; // Relative added mass (to displaced fluid volume) [-]
   double m_ShankThickness; // Total length of anchor [m]
   double m_ShankFlukeAngle; // Length from anchor-tip to fluke [m]
   double m_FlukeWidth; // Fluke width times fluke height [m^2]
   double m_FlukeSideLength; // Main diameter of shank [m]
   double m_FlukePointLength;
   double m_FlukeThickness; // Anchor drag coefficient [-]
 
   ISignalPort* m_Force;
   int m_PositionIndex;
   int m_VelocityIndex;
   int m_QuaterIndex;
   int m_OmegaIndex;
 
   environment::EnvironmentProvider* m_environment;
   std::vector<rigidbody::soilforces::SoilElement*> m_soilElements; 
 
#ifdef FH_VISUALIZATION
   Ogre::SceneNode* m_RenderNode; 
   Ogre::SceneNode* m_seabed;
#endif
};
 
 
#endif