doc/CNetElement3NWakeEffect_8h_source.html

#ifndef C_NET_ELEMENT_3N_WAKE_EFFECT

#define C_NET_ELEMENT_3N_WAKE_EFFECT


#include <math.h>

#include "CNetElementDef.h"

#include "CPrintDuringExec.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 CNetElement3NWakeEffect

{

public:

    CNetElement3NWakeEffect(

        const std::string       sNetName,

        const unsigned long     iID,

        const CNetElementDef*   pNetPanel,

        const int               NBarsInWakeEffectCalc,

        const double            dNy_water,

        const double            dRho_water,

        const double            dKnotDiameterRatio,

        const double            dKnotMomStiff,

        const double            dKknotMomStiff_contact,

        const double            dCnLinearLimitAngle,

        const double            dCnLinear,

        const double            dCnTurb,

        const double            dCnLam,

        const double            dTurbLimit,

        const double            dLamLimit,

        const double            dCt_nominal,

        const double            dCnKnots_nominal,

        const double            dDampingRatio,

        const bool              internalWakeEffectsActive,

        const bool              useActualSizeInWakeEffect,

        const bool              useReynoldsDependentDrag);


    ~CNetElement3NWakeEffect(void);


    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          dUseScreenModel,

        const bool          dUseDragCoeffMoeFore2021,

        const double        dSnImg,

        const double        dRhoWater,

        double              dHydroForcesRatio = 1.0,

        double              dAddedLinearDrag = 0.0

        );


    double CalcNodeInertia(int iNode);


    //****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* pOgreRoot);


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


    // 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;

    double  *m_ExternalWakeEffectPolynomial;

    int     m_orderOfInternalPolynomial;

    int     m_orderOfExternalPolynomial;

    bool    m_InternalWakeEffectsActive;

    bool    m_isInShadowedRegionOfCage;

    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_iID;

    std::string m_sNetName;

    double m_ny_water;

    double m_adNodeInertia[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_dBarD;

    double m_dKnotD;

    double m_dBarL0;

    double m_dKnotMomStiff;

    double m_dKnotMomStiffContact;

    double m_adNodeAPos_mesh[2];

    double m_adNodeBPos_mesh[2];

    double m_adNodeCPos_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_adNodesDistAB_mesh[2];

    double m_adNodesDistAC_mesh[2];

    double m_adNodesDistBC_mesh[2];

    double m_dMeshDet;

    double m_dNumMeshes;

    double m_dNumUBars;

    double m_dNumVBars;

    double m_dNumKnots;

    double m_dWeight;

    double m_dEATwine;


    double m_dMeshOpeningAngle;

    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_dDampingCoeff;

#ifdef FILTER_FORCES

    double m_nodeAforce_ned_last[3];

    double m_nodeBforce_ned_last[3];

    double m_nodeCforce_ned_last[3];

#endif


    bool                    m_isCoveredByROV;


#ifdef FH_VISUALIZATION

    // Visualization specifics

    Ogre::Root*             m_pOgreRoot;

    Ogre::SceneManager*     m_pSceneMgr;

    Ogre::Entity*           m_apRenderElement[3];

    Ogre::Entity*           m_apRenderSphere[3];

    Ogre::SceneNode*        m_pRenderNodeMesh;

    Ogre::SceneNode*        m_apRenderNodeSphere[3];

    Ogre::SceneNode*        m_pRenderNodeOutline;

    Ogre::SceneNode*        m_pRenderNodeMaterial;

#ifdef USE_DYNAMIC_LINES

    CDynamicLines*          m_pLines;

#else

    C3DLine*                m_pLines;

#endif

    CDynamicLines*          m_pLineOutline;

    Ogre::ManualObject *    m_pNetRenderObj;


    // Visualization variables

    unsigned short*         m_aiSortU;

    unsigned short*         m_aiSortV;

    unsigned long               m_iNumTwines;

    unsigned long               m_iNumUTwines;

    unsigned long               m_iNumVTwines;

    double*                 m_adInterpCoeffAu;

    double*                 m_adInterpCoeffAv;

    double*                 m_adInterpCoeffBu;

    double*                 m_adInterpCoeffBv;

#endif


};


#endif

CNetElement3NWakeEffect
Definition: CNetElement3NWakeEffect.h:42

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

CNetElement3NWakeEffect::CNetElement3NWakeEffect
CNetElement3NWakeEffect(const std::string sNetName, const unsigned long iID, const CNetElementDef *pNetPanel, const int NBarsInWakeEffectCalc, const double dNy_water, const double dRho_water, const double dKnotDiameterRatio, const double dKnotMomStiff, const double dKknotMomStiff_contact, const double dCnLinearLimitAngle, const double dCnLinear, const double dCnTurb, const double dCnLam, const double dTurbLimit, const double dLamLimit, const double dCt_nominal, const double dCnKnots_nominal, const double dDampingRatio, const bool internalWakeEffectsActive, const bool useActualSizeInWakeEffect, const bool useReynoldsDependentDrag)
The constructor.

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

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