Tube.hpp

#pragma once
 
// #define _USE_MATH_DEFINES
#include <cmath>
#include <vector>
#include <iostream>
#include <blaze/Math.h>
 
class Tube
{
private:
    static constexpr double pi_64 = M_PI / 64.00;
 
    static constexpr double pi_32 = M_PI / 32.00;
 
    double m_OD;
 
    double m_ID;
 
    double m_E;
 
    double m_I;
 
    double m_G;
 
    double m_J;
 
    blaze::DiagonalMatrix<blaze::StaticMatrix<double, 3UL, 3UL, blaze::rowMajor>> m_K;
 
    double m_ls;
 
    double m_lc;
 
    blaze::StaticVector<double, 3UL> m_u_ast;
 
public:
    Tube();
 
    Tube(double OD, double ID, double E, double G, double ls, double lc, const blaze::StaticVector<double, 3UL> &u_ast);
 
    ~Tube() = default;
 
    Tube(const Tube &rhs);
 
    Tube(Tube &&rhs) noexcept;
 
    Tube &operator=(const Tube &rhs);
 
    Tube &operator=(Tube &&rhs) noexcept;
 
    std::tuple<double, double, double, double, double, blaze::StaticVector<double, 3UL>> getTubeParameters();
 
    void setYoungModulus(double E);
 
    void setShearModulus(double G);
 
    double getTubeLength();
 
    blaze::StaticVector<double, 3UL> get_u_ast();
 
    double get_u_ast(const size_t id);
 
    void set_u_ast(const blaze::StaticVector<double, 3UL> &u_ast);
 
    void set_u_ast(const size_t id, const double u);
 
    double getStraightLen();
 
    double getCurvLen();
 
    void setStraightLen(double ls);
 
    void setCurvLen(double lc);
 
    blaze::DiagonalMatrix<blaze::StaticMatrix<double, 3UL, 3UL, blaze::rowMajor>> getK_Matrix();
 
    double getK(int i);
 
    void setK(const double EI, const double GJ);
 
    void setBendingK(const double EI);
};