b2element_beam_cross_section.H

//------------------------------------------------------------------------
// b2element_beam_cross_section.H --
//
// written by Mathias Doreille
//            Thomas Blome <thomas.blome@dlr.de>
//
// Copyright (c) 2011-2012 SMR Engineering & Development SA
//                         2502 Bienne, Switzerland
//
// (c) 2023 Deutsches Zentrum für Luft- und Raumfahrt (DLR) e.V.
//          Linder Höhe, 51147 Köln
//
// All Rights Reserved.  Proprietary source code.  The contents of
// this file may not be disclosed to third parties, copied or
// duplicated in any form, in whole or in part, without the prior
// written permission of SMR.
//------------------------------------------------------------------------
 
#ifndef B2ELEMENT_BEAM_CROSS_SECTION_H_
#define B2ELEMENT_BEAM_CROSS_SECTION_H_
 
#include "model/b2element.H"
#include "utils/b2linear_algebra.H"
 
namespace b2000 {
 
class BeamCrossSectionBaseElement : public TypedElement<double> {
public:
    virtual void add_value_0(
          Model* model, double& area, double neutral_point[2], double& mass,
          double inertia_point[2]) = 0;
 
    virtual void add_value_01(
          Model* model, const double center[2], double neutral_axis[3],
          double inertia_moment[3]) = 0;
 
    virtual void add_value_1(
          Model* model, const double neutral_point[2], const double inetria_axis[2][2],
          const size_t number_of_disp_dof, b2linalg::Index& dof_numbering,
          b2linalg::Matrix<double, b2linalg::Mpacked>& K) = 0;
 
    virtual void add_value_2(
          Model* model, const double neutral_point[2], const double neutral_axis[2][2],
          const b2linalg::Matrix<double, b2linalg::Mrectangle_constref>& v,
          const b2linalg::Matrix<double, b2linalg::Mrectangle_constref>& lambda,
          const bool beta_cd[6][6], const int beta_cs, const int beta_c[6], const int beta_ds,
          const int beta_d[6], b2float128 beta[6][6]) = 0;
 
    virtual void add_value_3(
          Model* model, const double neutral_point[2], const double neutral_axis[2][2],
          const size_t number_of_disp_dof,
          const b2linalg::Matrix<double, b2linalg::Mrectangle_constref>& v,
          const b2linalg::Matrix<double, b2linalg::Mrectangle_constref>& lambda,
          const b2float128 beta[6][6], b2linalg::Matrix<double>& u) = 0;
 
    // for the free bc stress
    virtual void add_value_4(
          Model* model, const double neutral_point[2], const double neutral_axis[2][2],
          const b2linalg::Matrix<double, b2linalg::Mrectangle_constref>& v,
          const b2linalg::Matrix<double, b2linalg::Mrectangle_constref>& lambda,
          const b2float128 beta[6][6], std::vector<size_t>& node_weight,
          b2linalg::Matrix<double> strain[6], b2linalg::Matrix<double> stress[6]) = 0;
 
    // for the boundary condition
    virtual void add_value_5(
          Model* model, const double neutral_point[2], const double neutral_axis[2][2],
          double f[6]) = 0;
 
    virtual void add_value_6(
          Model* model, const double neutral_point[2], const double neutral_axis[2][2],
          const size_t number_of_disp_dof, const double f[6],
          const b2linalg::Matrix<double, b2linalg::Mrectangle_constref>& u,
          const b2linalg::Matrix<double> stress[6], b2linalg::Vector<double>& value) = 0;
 
    // for the bc stress contribution
    virtual void add_value_7(
          Model* model, const double neutral_point[2], const double neutral_axis[2][2],
          const double f[6], const b2linalg::Matrix<double, b2linalg::Mrectangle_constref>& u,
          const b2linalg::Vector<double, b2linalg::Vdense_constref>& vd,
          const b2linalg::Vector<double, b2linalg::Vdense_constref>& lambdad,
          b2linalg::Matrix<double>& straind, b2linalg::Matrix<double>& stressd) = 0;
 
    virtual void add_value_8(
          Model* model, const double neutral_point[2], const double neutral_axis[2][2],
          const b2linalg::Matrix<double, b2linalg::Mrectangle_constref>& u,
          const b2linalg::Vector<double, b2linalg::Vdense_constref>& vd,
          const b2linalg::Matrix<double> stress[6], const b2linalg::Matrix<double>& stressd,
          double betad[6]) = 0;
};
 
}  // namespace b2000
 
#endif /* B2ELEMENT_BEAM_CROSS_SECTION_H_ */