doc/NetElement3NWakeEffect_8h_source.html

#ifndef C_NET_ELEMENT_3N_WAKE_EFFECT

#define C_NET_ELEMENT_3N_WAKE_EFFECT


#include <math.h>

#include "NetElementDef.h"

#include <fhsim/PrintDuringExec.h>

#include <complex>

//#include "eigen_matrix_defs.h"

//#include "Eigen/Eigen"

#include "Eigen/StdVector"

using namespace std;

using std::vector;

typedef Eigen::Matrix<double,3,1> vec3;

//#include <Eigen/Core>


#ifdef FH_VISUALIZATION

//#define USE_DYNAMIC_LINES

#ifdef USE_DYNAMIC_LINES

#include "sfh/ogre/CDynamicLines.h"

#else

#include "sfh/ogre/C3DLine.h"

#endif

#endif


//#define RENDER_OUTLINE

#define RENDER_MESH

//#define ATTACH_SPHERE

//#define RENDER_MATERIAL

//#define FILTER_FORCES

//#define CNetElement3NWakeEffect_TRANSPARENT


class CDynamicLines; // forward declaration


class NetElement3NWakeEffect

{

public:

    NetElement3NWakeEffect(

        const std::string       netName,

        const unsigned long     iID,

        const NetElementDef*    netPanel,

        const int               NBarsInWakeEffectCalc,

        const double            dNy_water,

        const double            dRho_water,

        const double            knotDiameterRatio,

        const double            knotMomStiff,

        const double            dKknotMomStiff_contact,

        const double            cnLinearLimitAngle,

        const double            cnLinear,

        const double            cnTurb,

        const double            cnLam,

        const double            turbLimit,

        const double            lamLimit,

        const double            dCt_nominal,

        const double            dCnKnots_nominal,

        const double            dampingRatio,

        const bool              internalWakeEffectsActive,

        const bool              useActualSizeInWakeEffect,

        const bool              useReynoldsDependentDrag);


    ~NetElement3NWakeEffect(void);


    // Owns raw buffers (freed in the destructor), so it must not be copied (MARE-0318).

    NetElement3NWakeEffect(const NetElement3NWakeEffect&) = delete;

    NetElement3NWakeEffect& operator=(const NetElement3NWakeEffect&) = delete;


    void AddNodeForces(const double* const adNodePosA_ned,

        const double* const adNodePosB_ned,

        const double* const adNodePosC_ned,

        const double* const adNodeVelA_ned,

        const double* const adNodeVelB_ned,

        const double* const adNodeVelC_ned,

        const double* const adAmbientCurrentVel_ned,

        const double* const adWaveParticleVel_ned,

        const double* const adWaterVel_ned,

        double * const      adNodeForceA_ned,

        double * const      adNodeForceB_ned,

        double * const      adNodeForceC_ned,

        const bool          useScreenModel,

        const bool          dUseDragCoeffMoeFore2021,

        const double        snImg,

        const double        rhoWater,

        double              hydroForcesRatio = 1.0,

        double              addedLinearDrag = 0.0

        );


    double CalcNodeInertia(int node);


    //****NEW STUFF FOR WAKE EFFECT INCLUSION****//

    void    GenerateInternalWakeEffectPolynomial(double* currentVel);

    void    GenerateExternalWakeEffectPolynomial(double* currentVel);

    bool    IsInShadowedRegionOfCage();

    void    SetInShadowedRegionOfCage(bool value);

    void    SetDistanceToWallInFront(double distance);

    double  GetExternalWakeEffect(double* currentVel, double* posA, double* posB, double* posC);

    double  GetInternalWakeEffect(double* currentVel, double* posA, double* posB, double* posC);


    //****END NEW STUFF FOR WAKE EFFECT INCLUSION****//


#ifdef FH_VISUALIZATION

    void RenderInit(Ogre::Root* ogreRoot);


    void RenderUpdate( const double* const adPosA, const double* const adPosB,

        const double* const adPosC );

#endif


    void setElementIsCoveredByROV(bool _covered);


private:

    void CalcTensionForces( double meshBarComp_panel[6], double meshBarLength[2],

        double twineTension[2], double nodeATension_panel[3], double nodeBTension_panel[3],

        double nodeCTension_panel[3]);


    void CalcLocalUVComponents( double nodesABdist_panel[3], double nodesACdist_panel[3],

        double nodesBCdist_panel[3], double meshBarComp_panel[6], double meshBarLength[2]);


    void CalcMeshContactForces( double meshBarComp_panel[6],

        double meshBarLength[2],

        double nodeAMeshOpeningRes_panel[3],

        double nodeBMeshOpeningRes_panel[3],

        double nodeCMeshOpeningRes_panel[3],

        double nodeATwineContactRes_panel[3],

        double nodeBTwineContactRes_panel[3],

        double nodeCTwineContactRes_panel[3]);


    double CalcVelocities(double R_ned2panel[3][3],

        const double* const adNodeVelA_ned,

        const double* const adNodeVelB_ned,

        const double* const adNodeVelC_ned,

        const double* const adWaterVel_ned,

        double              elementVel_ned[3],

        double              hydDynVel_ned[3],

        double              hydDynVel_panel[3]);


    void CalcHydroDynamicForcesScreenModel(

        const double* const adNodePosA_ned,

        const double* const adNodePosB_ned,

        const double* const adNodePosC_ned,

        const double* const adNodeVelA_ned,

        const double* const adNodeVelB_ned,

        const double* const adNodeVelC_ned,

        const double* const adAmbientCurrentVel_ned,

        const double* const adWaveParticleVel_ned,

        const double* const adWaterVel_ned,

        double * const      nodeAforce_ned,

        double * const      nodeBforce_ned,

        double * const      nodeCforce_ned,

        const double        rhoWater,

        const bool          useDragCoeffMoeFore2021,

        const double        SnImg);


    void CalcHydroDynamicForces(double meshBarComp_panel[6],

        double meshBarLength[2],

        double R_ned2panel[3][3],

        const double* const adNodeVelA_ned,

        const double* const adNodeVelB_ned,

        const double* const adNodeVelC_ned,

        const double* const adWaterVel_ned,

        double elementVel_ned[3],

        double totalDragPerNode_panel[3]);


    void CalcDampingForces( double elementVel_ned[3],

        const double* const adNodeVelA_ned,

        const double* const adNodeVelB_ned,

        const double* const adNodeVelC_ned,

        const double* const adWaterVel_ned,

        double adDampingForcesNodeA[3],

        double adDampingForcesNodeB[3],

        double adDampingForcesNodeC[3],

        double addedLinearDrag);


    // FUNCTIONS FOR SCREEN MODEL

    double DragCoeffCircularCylinderFromReynoldsNumber(double Sn, double Urel);

    vec3 CalcRelativeVelocity(

        vec3 vec_Ucurr,

        vec3 vec_Uwave,

        const double* const adNodeVelA_ned,

        const double* const adNodeVelB_ned,

        const double* const adNodeVelC_ned);


    //****NEW STUFF FOR WAKE EFFECT INCLUSION****//


    void    SortPanelNodes(int index, const double * aadNodes[3], double adNodePosA_mesh[2], double adNodePosB_mesh[2], double adNodePosC_mesh[2], unsigned short* aiSort);

    void    FindTwineCrossings(unsigned short index, unsigned short kiNumMeshPerLine, double* CoeffA, double* CoeffB, const double* aadNodes[3]);


    double  GetDragCoefficientFromPolynomial(double nominalReynoldsNumber);

    double  Blevins(double DX, double DY, double Cd);

    double  SquaredSumColumnIn2DDoubleArray(Eigen::MatrixXd input,int colIndex,int length);

    Eigen::VectorXd PolynomialFit(Eigen::VectorXd alpha, int nAlpha, Eigen::VectorXd targetArray, int order);

    double  ComputePolynomialValue(std::string polynomialType, double inputValue);

    int     m_maxUBarDimension;

    int     m_maxVBarDimension;

    double  *m_InternalWakeEffectPolynomial = nullptr;

    double  *m_ExternalWakeEffectPolynomial = nullptr;

    int     m_orderOfInternalPolynomial;

    int     m_orderOfExternalPolynomial;

    bool    m_InternalWakeEffectsActive;

    bool    m_inShadowedRegionOfCage;

    double  m_distanceToWallInFront;

    bool    m_useReynoldsDependentDrag;

    bool    m_UseActualSizeInWakeEff;

    int     m_NBarsInWakeEffectCalc;

    //****END NEW STUFF FOR WAKE EFFECT INCLUSION****//


    double  m_Sn;

    double  m_knotFactor;


    // Input

    unsigned long  m_iD;

    std::string m_netName;

    double m_ny_water;

    double m_nodeInertia[3];

    //double m_rho_water;                   ///< Density of water.

    //double m_rho_twines;              ///< Material density of twines.

    //double m_EModulus;                    ///< E-modulus of twine material.

    double m_barD;

    double m_knotD;

    double m_barL0;

    double m_knotMomStiff;

    double m_knotMomStiffContact;

    double m_nodeAPos_mesh[2];

    double m_nodeBPos_mesh[2];

    double m_nodeCPos_mesh[2];

    double m_CnLinearLimitAngle;

    double m_CnLinear;

    double m_CnTurb;

    double m_CnLam;

    double m_TurbLimit;

    double m_LamLimit;

    double m_Ct_nominal;

    double m_CnKnots_nominal;


    // Calculated element properties

    //double m_twineCrossSecArea;           ///< The cross sectional area of each twine.

    double m_nodesDistAB_mesh[2];

    double m_nodesDistAC_mesh[2];

    double m_nodesDistBC_mesh[2];

    double m_meshDet;

    double m_numMeshes;

    double m_numUBars;

    double m_numVBars;

    double m_numKnots;

    double m_weight;

    double m_eATwine;


    double m_meshOpeningAngle;

    double m_ShadowEffect[3];

    double m_ShadowEffectBias;

    // Calculated constants for increased computational performance.

    double m_Ct;

    double m_kKnot;

    double m_NodeWeight;

    double m_kTwineDragNormal[2];

    double m_nodeAforce_ned[3];

    double m_nodeBforce_ned[3];

    double m_nodeCforce_ned[3];

    double m_dampingCoeff;

#ifdef FILTER_FORCES

    double m_nodeAforce_ned_last[3];

    double m_nodeBforce_ned_last[3];

    double m_nodeCforce_ned_last[3];

#endif


    bool                    m_coveredByROV;


#ifdef FH_VISUALIZATION

    // Visualization specifics

    Ogre::Root*             m_ogreRoot;

    Ogre::SceneManager*     m_sceneMgr;

    Ogre::Entity*           m_renderElement[3];

    Ogre::Entity*           m_renderSphere[3];

    Ogre::SceneNode*        m_renderNodeMesh;

    Ogre::SceneNode*        m_renderNodeSphere[3];

    Ogre::SceneNode*        m_renderNodeOutline;

    Ogre::SceneNode*        m_renderNodeMaterial;

#ifdef USE_DYNAMIC_LINES

    CDynamicLines*          m_lines = nullptr;

#else

    C3DLine*                m_lines = nullptr;

#endif

    CDynamicLines*          m_lineOutline;

    Ogre::ManualObject *    m_netRenderObj;


    // Visualization variables

    unsigned short*         m_sortU = nullptr;

    unsigned short*         m_sortV = nullptr;

    unsigned long               m_numTwines;

    unsigned long               m_numUTwines;

    unsigned long               m_numVTwines;

    double*                 m_interpCoeffAu = nullptr;

    double*                 m_interpCoeffAv = nullptr;

    double*                 m_interpCoeffBu = nullptr;

    double*                 m_interpCoeffBv = nullptr;

#endif


};


#endif

NetElement3NWakeEffect
Definition NetElement3NWakeEffect.h:48

NetElement3NWakeEffect::CalcNodeInertia
double CalcNodeInertia(int node)
Calculates the inertia of a node.

NetElement3NWakeEffect::NetElement3NWakeEffect
NetElement3NWakeEffect(const std::string netName, const unsigned long iID, const NetElementDef *netPanel, const int NBarsInWakeEffectCalc, const double dNy_water, const double dRho_water, const double knotDiameterRatio, const double knotMomStiff, const double dKknotMomStiff_contact, const double cnLinearLimitAngle, const double cnLinear, const double cnTurb, const double cnLam, const double turbLimit, const double lamLimit, const double dCt_nominal, const double dCnKnots_nominal, const double dampingRatio, const bool internalWakeEffectsActive, const bool useActualSizeInWakeEffect, const bool useReynoldsDependentDrag)
The constructor.

NetElement3NWakeEffect::~NetElement3NWakeEffect
~NetElement3NWakeEffect(void)
The destructor.

NetElementDef
Class containing the description of a net panel and methods necessary for taking this over to triangu...
Definition NetElementDef.h:22