beam_element.hpp

Go to the documentation of this file.

#pragma once
 
#include <array>
#include <span>
#include <vector>
 
#include "beam_section.hpp"
 
namespace kynema {
 
struct BeamElement {
    size_t ID;                                      // Element identifier
    std::vector<size_t> node_ids;                   // Element node identifiers
    std::vector<BeamSection> sections;              // Element mass/stiffness in material frame
    std::vector<std::array<double, 2>> quadrature;  // Element quadrature points and weights
    std::array<double, 6> mu;  // Element stiffness-proportional damping coefficients
 
    BeamElement(
        size_t id, std::span<const size_t> n, std::span<const BeamSection> s,
        std::span<const std::array<double, 2>> q,
        const std::array<double, 6>& mu_ = {0., 0., 0., 0., 0., 0.}
    )
        : ID(id), mu(mu_) {
        node_ids.assign(std::begin(n), std::end(n));
        sections.assign(std::begin(s), std::end(s));
        quadrature.assign(std::begin(q), std::end(q));
    }
};
 
}  // namespace kynema