Cable.h

 
#ifndef CABLE_H
#define CABLE_H
#include "SimObject.h"
class CEnvironment;
class C3DLine;
 
class CCable : public SimObject
{
   public:
   CCable(const string& simObjectName, ISimObjectCreator* const creator);
 
   ~CCable(); 
 
   void OdeFcn(const double T, const double* const X, double* const XDot, const bool bIsMajorTimeStep);
 
 
   void InitialConditionSetup(const double T, const double* const currentIC, double* const updatedIC, ISimObjectCreator* const creator);
 
   void FinalSetup(const double T, const double* const X, ISimObjectCreator* const creator); 
 
   const double* ForceA(const double T, const double* const X); 
   const double* ForceB(const double T, const double* const X); 
 
   void computeEndForces(const double T, const double* const X); 
 
#ifdef FH_VISUALIZATION
   void RenderInit(Ogre::Root* const ogreRoot, ISimObjectCreator* const creator); 
   void RenderUpdate(const double T, const double* const X); 
#endif
 
   protected:
   void computeSegmentForces(const double T, const double* X, const double* pA, const double* pB, const double* vA, const double* vB, double* FA, double* FB);
 
 
   void computeVirtualPoint(const double* const target, const double* const point, const double* const targetVel, const double* const pointVel, double* virtualPoint, double* virtualPointVel, double residual, double dResidual);
 
   void testPositiveParam(double value, ISimObjectCreator* creator, string valueName); 
   void computeElementTension(const double* prev, const double* curr, double* force); 
 
   double m_cableLength; 
   double m_cableWeight; 
   double m_Emodulus; 
   double m_tensileCrossSectionArea; 
   double m_displacementCrossSectionArea; 
   double m_diameter; 
   double m_hydrodynamicDiameter; 
   double m_submergedWeight; 
   double m_internalDamping; 
   int m_numElements; 
   double m_elementMass; 
 
   double m_elementMassInverse; 
   double m_segmentMass; 
   double m_segmentSubmergedMass; 
   double m_segmentLength; 
   double m_segmentStiffness; 
   int* m_elementPositionIndex; 
   int* m_elementVelocityIndex; 
 
   ISignalPort* m_PositionA; 
   ISignalPort* m_PositionB; 
   ISignalPort* m_VelocityA; 
   ISignalPort* m_VelocityB; 
   ISignalPort* m_retractedLengthPortA; 
   ISignalPort* m_retractedLengthPortB; 
   ISignalPort* m_retractedSpeedPortA; 
   ISignalPort* m_retractedSpeedPortB; 
 
   double m_endForceA[3]; 
   double m_endForceB[3]; 
   double m_nodeForceA[3]; 
   double m_nodeForceB[3]; 
   double m_retractedLengthA; 
   double m_retractedLengthB; 
   int m_retractedNodesA; 
   int m_retractedNodesB; 
   double m_virtualPointA[3]; 
   double m_virtualPointB[3]; 
   double m_virtualPointVelocityA[3]; 
   double m_virtualPointVelocityB[3]; 
   double m_residualA; 
   double m_residualB; 
   ICommonComputation* m_computeEndForces; 
 
   const double* m_posA; 
   const double* m_posB; 
   const double* m_velA; 
   const double* m_velB; 
 
   CEnvironment* m_environment; 
 
#ifdef FH_VISUALIZATION
   C3DLine* m_lines1; 
   C3DLine* m_lines2; 
#endif
};
#endif