ViscousShip.h

#ifndef CViscousShip_h__
#define CViscousShip_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;}
        // amplitude horizon as a fraction of natural period ~B
        void SetPhiAmplitudeHorizon(double horizon){phiAmplitudeHorizon = horizon;}
    
    
        void GetViscousAddedMass(double u, double v, double r, Eigen::Matrix3d &Mat);
    
        double GetViscousSurge(double u, double v, double r);
        double GetViscousSway(double u, double v, double r);
        double GetViscousYaw(double u, double v, double r);
    
        double GetViscousRoll(double phiDot, double u = 0);
        double EddieDamping( double phiDot, double u = 0);
        double BilgeKeelDamping( double phiDot, double u = 0);
        double FricitonDamping( double phiDot, double u = 0);
        double LiftDamping( double phiDot, double u = 0);
        double SkegDamping( double phiDot, double u = 0);
 
        double SkegLiftDamping( double phiDot, double u );
 
        void AddToRollTimeHistory( const double T, double RollAngle, const bool bIsMajorTimeStep );
        void MakeBildgeKeel( double m_BilgeKeelStartPosition, double m_BilgeKeelEndPosition, double m_BilgeKeelHeight );
        void MakeSkeg();
    
        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);
    
        double ManeuverReferenceSpeed() const { return m_maneuverReferenceSpeed; }
        void ManeuverReferenceSpeed(double val) { m_maneuverReferenceSpeed = val; }
    protected:
 
        double deltaWingKP(double A_ratio);
        double deltaWingKV(double A_ratio);
        double deltaWingCL(double alfa, double Kp, double Kv);
 
        double m_maneuverReferenceSpeed;
        double m_viscousDampingBlender;
    
        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_bilgeKeelStartSection;
        int m_bilgeKeelEndSection;
        double m_startSectionFractionBK;
        double m_endSectionFractionBK;  
 
        int rollDampingTheory;
 
        int m_withSkeg;
        int SkegStartSection;
        int SkegEndSection;
        int *SkegStartPoints;
        int *SkegEndPoints;
    
    protected:
        void ReadMainDimensionsVeres( std::string fname );
        void ReadMainDimensionsWamit( std::string fname );
        HullData *myHullData;
    
        double surfaceOfHull;
        double kinematicViscosity;
    
        double *SectionalNormalBilgeKeelForce;
        double *SectionalPressureOnHullFromBilgeKeels;
        double *SectionalEddieForce;
        double *bilgeKeelHeight;
        
        double* Ds;
        double* Bs;
        double* H0;
        double* SigmaS;
        double* A1;
        double* A3;
        double* Ms;
    
        double Lpp,B,T;
        double Cb,VCG;
    
        double Lo,Lr,Kn;
    
        double* SectionLength;
        std::deque<double> phi_history;
        std::deque<double> time_history;
        double phiAmplitude;
        double phiAmplitudeTime;
        double m_rf;
        double OG;
        double phiAmplitudeHorizon;
        int numberOfStations;
        bool m_skegLift;
        bool m_hullLift;
    };
}
#endif // CViscousShip_h__