ViscousShip.h

#ifndef CViscousShip_h__
#define CViscousShip_h__
 
#include <vector>
#include "HelperFunctions.h"
#include <deque>
#include <string>
 
#include <Eigen/Eigen>
 
#include "HullData/HullData.h"
namespace Ship{
    class ViscousShip
    {
    public :
        enum EFileType
        {
            FILE_TYPE_VERES,
            FILE_TYPE_WAMIT
        };
    
        ViscousShip();
    
        void SetupFromFile( std::string fname, EFileType type, HullData *H);
    
        void SetHullData( HullData *HD ){ myHullData = HD; }
        // 0 = ITTC, 1= SeaWay
        void SetTheoryFormulation( int theory ){ rollDampingTheory = theory;}
        // The window [s] of the roll amplitude: the amplitude is the largest |roll| of the last @p horizon seconds. It is also
        // the roll period of the ITTC friction coefficient. Must be positive (set it before the first AddToRollTimeHistory)
        void SetManoeuvreMinSpeed(double speed){m_manoeuvreMinSpeed = speed;}
        void SetPhiAmplitudeHorizon(double horizon){phiAmplitudeHorizon = horizon;}
    
    
        void GetViscousAddedMass(double u, double v, double r, Eigen::Matrix3d &Mat) const;
    
        double GetViscousSurge(double u, double v, double r) const;
        double GetViscousSway(double u, double v, double r) const;
        double GetViscousYaw(double u, double v, double r) const;
    
        double GetViscousRoll(double phiDot, double u = 0) const;
        bool IsSkegStation( int station ) const;
        double EddieDamping( double phiDot, double u = 0) const;
        double BilgeKeelDamping( double phiDot, double u = 0) const;
        double FricitonDamping( double phiDot, double u = 0) const;
        double LiftDamping( double phiDot, double u = 0) const;
        double SkegDamping( double phiDot, double u = 0) const;
 
        double SkegLiftDamping( double phiDot, double u ) const;
 
        double RollAmplitude() const { return phiAmplitude; }
        void AddToRollTimeHistory( const double T, double RollAngle, const bool bIsMajorTimeStep );
        bool MakeBildgeKeel( double m_BilgeKeelStartPosition, double m_BilgeKeelEndPosition, double m_BilgeKeelHeight );
        void SetSkeg( int startSection, int endSection, const std::vector<int>& startPoints, const std::vector<int>& endPoints );
        bool HullHasTheSkegStations() const;
    
        void SetHorizontalDampingBlend(double blendSpeed, double blendFraction=0.5);
        void SetDpDampingX(double Xu, double Xv, double Xr);
        void SetDpDampingY(double Yu, double Yv, double Yr);
        void SetDpDampingN(double Nu, double Nv, double Nr);
    
        const std::vector<double>& GetSectionLength() const { return SectionLength; }
 
        double ManeuverReferenceSpeed() const { return m_maneuverReferenceSpeed; }
        void ManeuverReferenceSpeed(double val) { m_maneuverReferenceSpeed = val; }
    protected:
 
        double deltaWingKP(double A_ratio) const;
        double deltaWingKV(double A_ratio) const;
        double deltaWingCL(double alfa, double Kp, double Kv) const;
 
        double ManoeuvreSpeed(double u) const;
        double m_manoeuvreMinSpeed = 0.1;
        double m_maneuverReferenceSpeed;
        double m_viscousDampingBlender;
    
    public:
        // The coefficients the ship sets from its parameters
        double Yv_man, Yr_man, Yrd_man, Yvd_man;
        double Nv_man, Nr_man, Nrd_man, Nvd_man;
        double Xrr_man, Xud_man, Xvr_man, Xvv_man, Xvvvv_man;
    
        double Xu_dp, Yu_dp, Nu_dp;
        double Xv_dp, Yv_dp, Nv_dp;
        double Xr_dp, Yr_dp, Nr_dp;
        int m_withSkeg;
        bool m_skegLift;
        bool m_hullLift;
    protected:
    
        int m_bilgeKeelStartSection;
        int m_bilgeKeelEndSection;
        double m_startSectionFractionBK;
        double m_endSectionFractionBK;  
 
        int rollDampingTheory;
 
        int SkegStartSection;
        int SkegEndSection;
        std::vector<int> SkegStartPoints;
        std::vector<int> SkegEndPoints;
    
    protected:
        void ReadMainDimensionsVeres( std::string fname );
        void ReadMainDimensionsWamit( std::string fname );
        HullData *myHullData;
    
        double surfaceOfHull;
        double kinematicViscosity;
    
        std::vector<double> bilgeKeelHeight;
        
        std::vector<double> Ds;
        std::vector<double> Bs;
        std::vector<double> H0;
        std::vector<double> SigmaS;
        std::vector<double> A1;
        std::vector<double> A3;
        std::vector<double> Ms;
    
        double Lpp = kNotRead, B = kNotRead, T = kNotRead;
        double Cb = kNotRead, VCG = kNotRead;
    
        double Lo,Lr,Kn;
    
        std::vector<double> SectionLength;
        std::deque<double> phi_history;
        std::deque<double> time_history;
        double phiAmplitude;
        double phiAmplitudeTime;
        double m_rf;
        double OG;
        double phiAmplitudeHorizon = 2.75*3.1622776601683795; 
        int numberOfStations;
    };
}
#endif // CViscousShip_h__