beam.hpp

Go to the documentation of this file.

#pragma once
 
#include <array>
#include <vector>
 
#include "interfaces/node_data.hpp"
 
namespace kynema {
class Model;
struct BeamSection;
}  // namespace kynema
 
namespace kynema::interfaces::components {
 
struct BeamInput;
 
class Beam {
public:
    using DeviceType =
        Kokkos::Device<Kokkos::DefaultExecutionSpace, Kokkos::DefaultExecutionSpace::memory_space>;
 
    static constexpr size_t invalid_id{9999999};
 
    static constexpr size_t kMaxGeometryPoints{10};
 
    size_t beam_element_id{invalid_id};
 
    std::vector<NodeData> nodes;
 
    size_t prescribed_root_constraint_id{invalid_id};
 
    std::vector<double> node_xi;
 
    Beam() = default;
 
    Beam(const BeamInput& input, Model& model);
 
    [[nodiscard]] std::vector<double> GetNodeWeights(double s) const;
 
    void AddPointLoad(double s, std::array<double, 6> loads);
 
    void ClearLoads();
 
    void GetMotion(const HostState<DeviceType>& host_state);
 
    void SetLoads(HostState<DeviceType>& host_state) const;
 
private:
    std::vector<std::array<double, 3>> node_coordinates;  
    std::vector<std::array<double, 3>> node_tangents;     
 
    static void ValidateInput(const BeamInput& input);
    void SetupNodeLocations(const BeamInput& input);
 
    void CreateNodeGeometry(const BeamInput& input);
    void CreateBeamElement(const BeamInput& input, Model& model);
    void PositionBladeInSpace(const BeamInput& input, Model& model) const;
 
    void SetupRootNode(const BeamInput& input, Model& model);
 
    void CalcNodeTangents();
    static std::vector<BeamSection> BuildBeamSections(const BeamInput& input);
};
 
}  // namespace kynema::interfaces::components