CNetStructureArray.h

#pragma once
 
// Includes
#include <string>
#include "SimObject.h"
#include "CNetCableElement.h"
#include "CNetElementDef.h"
#include "CNormalForceCalc.h"
#include "IBasicBottom.h"
#include "CNetBottomContactFilter.h"
#include "CInputReader.h"
#include "ICurrentVel.h"
#include "SimObject.h"
 
#ifdef FH_VISUALIZATION
#include <sfh/ogre/NMeshGenerators.h>
#endif
 
//Class definition
class CNetStructureArray : public SimObject
{
public:
    CNetStructureArray(std::string sSimObjectName, ISimObjectCreator *pCreator);
    virtual ~CNetStructureArray();
    virtual void OdeFcn(const double dT, const double *const adX, double *const adXDot, const bool bIsMajorTimeStep);
    virtual void InitialConditionSetup(const double dT, const double *const adCurrentIC, double *const adUpdatedIC, ISimObjectCreator *const pCreator);
    virtual void FinalSetup(const double dT, const double *const adX, ISimObjectCreator *const pCreator);
#ifdef FH_VISUALIZATION
    virtual void RenderInit(Ogre::Root *const pOgreRoot, ISimObjectCreator *const pCreator);
    virtual void RenderUpdate(const double dT, const double *const adX);
#endif
 
    void SetAddedDragPerMPS(double AddedDragPerMPS);
    virtual const double * NumExternalPos(const double dT, const double *const adX);
    virtual const double * NodePos(const double dT, const double *const adX, int iExtNode);
    virtual const double * NodeVel(const double dT, const double *const adX, int iExtNode);
    virtual const double * BottomContact(const double dT, const double *const adX);
protected:
    virtual void CalcStateDerivatives(const double *const adX, double *const adXDot);
    virtual void AddExternalForces(const double dT, const double *const adX);
    virtual void DeleteAllNodeForces();
    virtual void ReadInputFile(ISimObjectCreator *pCreator);
    virtual void AddNetForces(const double dT, const double *const adX);
    virtual void AddCableForces(const double dT, const double *const adX);
    void FindNodesToDivide(int iNodesToDivide[2], double *dLength, double *dLengthFactor);
    void RefineMesh();
    void CheckNodeNumber(int iNodeNumber, ISimObjectCreator *const pCreator, std::string sNodeDescription = "");
    void CheckNodes(ISimObjectCreator *pCreator);
    unsigned short AddCalcNode(int aiNodes[2]);
    virtual double AddNodeBottomForce(const double *const adX, int iNode, double dDiam, double adForceInOut[3]);
    virtual void AddBottomForces(const double *const adX);
    virtual void AddInputPorts(ISimObjectCreator *pCreator);
 
    sfh::timers::StopWatch m_OdeTimer;
 
    // Node indices
    int                     m_iNumNodeInputsUsed;   
    int                     m_iNumNodeOutputsUsed;  
 
    int*                    m_aIExternalConn;           
    long*               m_aiExternalNodes;          
 
    // Input ports
    ISignalPort**           m_apInExternalForces;           
 
    // Output ports
    double                  m_dNumExternalPos;          
 
    // States
    int*                    m_aIStatePos;           
    int*                    m_aIStateVel;           
 
    // Calculation variables
    int                     m_iNumPanelsCalc;       
    int                     m_iNumPanelsGlobal;     
    int                     m_iNumCablesCalc;       
    int                     m_iNumCablesGlobal;     
    int                     m_iNumNodesOriginal;    
    int                     m_iNumNodesCalc;        
    double                  m_dFluidRho;            
    std::vector<CNetElementDef>     m_VSPanelData;          
    std::vector<CNetCableElement*>  m_VCableElements;       
    CNetBottomContactFilter m_BottomContactFilter;
    double                  m_adBottomContact[6];
    // Variables for controlling the calculation mesh refinement.
    double                  m_dMaxSideLengthFactor; 
    double                  m_dMinSideLength;       
    int                     m_iMaxNumNodes;         
    int                     m_iMaxNumPanels;        
    int                     m_iMaxNumCables;        
    double                          m_dScale;               
    CNetCableElement**              m_apCableElements;      
    double*                         m_adInertia;            
    double*                         m_adNodeForce;          
    std::string                     m_sNetName;             
    IBasicBottom **m_ppBottom;
    ICurrentVel**                   m_ppCurrentVel;         
    // Net element properties.
    int *                   m_aiNetConn;
    // Cable properties.
    int **m_aaiCableConn;
    // Variables controlling the simulation performance and adaptation.
    double                          m_dMaxStepLength;       
    double                          m_dAdaptationPeriod;    
    double                          m_dStepSafetyFactor;    
    double                          m_dMinCableLength;      
    double                          m_dMeanCableTension;    
    double                          m_dMaxCableTension;     
    double                          m_dMassFactor;          
    double                          m_dMaxAcceleration;     
    double                          m_dMaxVelocity;         
    double                          m_dAddedDragPerMPS;     
 
    int m_iNumNodeOutputs;
    VSPairStrInt        m_VsNodesOutput;        
    VSPairStrInt        m_VsNodesInput;     
    VSPairStrInt        m_VsNodesAvailable;     
 
    bool m_bUseTagForExternalNodes;
    int CreateExternalNodeMap(VSPairStrInt& VsNodes, std::string sTag, ISimObjectCreator* pCreator);
    void CreateNetPanels( ISimObjectCreator* pCreator);
    void CreateNetCables( ISimObjectCreator* pCreator);
    void CreateDataStructures();
    void CreateOutputPorts(ISimObjectCreator* pCreator);
 
 
    // Net element data structures
   // double*       m_adNetElementNodeForces;
    double*         adNodesDistMesh;
    double*         adMesdDet;
    double*         adEATwine;
    double*         adMeshBarLength_unstretched;
    double*         adKnotDiameter;
    double*         adKnotMomStiff_contact;
    double*         adKKnot;
    double*         adTwineDragNormal;
    double*         adTwineDiameter;
    double*         adElementWaterVel;
    double*         adDampingCoeff;
    double*         adElementNodeWeight;
 
    void PRINTALLNODEFORCES(std::string spec, double dT);
};