beam_input.hpp

Go to the documentation of this file.

#pragma once
 
#include <array>
#include <cstdint>
#include <stdexcept>
#include <vector>
 
namespace kynema::interfaces::components {
 
struct ReferenceAxis {
    std::vector<double> coordinate_grid;
 
    std::vector<std::array<double, 3>> coordinates;
 
    std::vector<double> twist_grid;
 
    std::vector<double> twist;
};
 
struct Section {
    double location;
 
    std::array<std::array<double, 6>, 6> mass_matrix;
 
    std::array<std::array<double, 6>, 6> stiffness_matrix;
 
    Section(
        double loc, std::array<std::array<double, 6>, 6> m, std::array<std::array<double, 6>, 6> k
    )
        : location(loc), mass_matrix(m), stiffness_matrix(k) {
        // Check that the section location is in range [0,1]
        if (loc < 0. || loc > 1.) {
            throw std::invalid_argument("Section location must be in range [0, 1]");
        }
    }
};
 
struct Root {
    bool prescribe_root_motion{false};
 
    std::array<double, 7> position{0., 0., 0., 1., 0., 0., 0.};
 
    std::array<double, 6> velocity{0., 0., 0., 0., 0., 0.};
 
    std::array<double, 6> acceleration{0., 0., 0., 0., 0., 0.};
};
 
struct BeamInput {
    enum class QuadratureRule : std::uint8_t {
        GaussLegendre,
        GaussLobatto
    };
 
    enum class QuadratureStyle : std::uint8_t {
        Segmented,
        WholeBeam
    };
 
    size_t element_order{10};
 
    QuadratureRule quadrature_rule{QuadratureRule::GaussLobatto};
 
    QuadratureStyle quadrature_style{QuadratureStyle::Segmented};
 
    size_t section_refinement{0};
 
    ReferenceAxis ref_axis;
 
    Root root;
 
    std::vector<Section> sections;
};
 
}  // namespace kynema::interfaces::components