b2element_klt_shell.H

//------------------------------------------------------------------------
// b2element_klt_shell.H --
//
//     Shell element implementing Kirchhoff-Love theory, continuum
//     element approach.
//
// written by Thomas Ludwig
//            Thomas Blome <thomas.blome@dlr.de>
//
// Copyright (c) 2016 SMR Engineering & Development SA
//                    2502 Bienne, Switzerland
//
// Copyright (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_KLT_SHELL_H_
#define _B2ELEMENT_KLT_SHELL_H_
 
#include "model/b2element.H"
#include "utils/b2linear_algebra.H"
 
namespace b2000 {
 
using namespace klt;
using namespace b2linalg;
using namespace node;
using namespace coordof;
 
class ElementShellKLTBase : public TypedElement<double> {
public:
    virtual ElementProperty* get_property() = 0;
 
    virtual int get_num_edges() = 0;
 
    virtual void compute_nodes_interpolation_edge_g1(
          const int edge, const double xi, Vector<double, Vdense>& N) const = 0;
 
    virtual int get_edge_g1_order(const int edge) const = 0;
 
    virtual void get_xi_all_from_xi_edge(const int edge, const double xi, double xi_all[2]) = 0;
 
    virtual void get_dof_numbering_edge_g1(const int edge, Index& dof_numbering) const = 0;
 
    virtual void compute_dof_shape_functions(
          const double xi[2], Matrix<double, Mrectangle>& N) const = 0;
 
    virtual void compute_dof_shape_functions_edge_g1(
          const int edge, const double xi, Matrix<double, Mrectangle>& N) const = 0;
 
    virtual void compute_geom(
          const double thickness, const double offset, const double xi[2], double pos[3],
          double Ar[3], double As[3]) const = 0;
 
    virtual void compute_geom(
          const double thickness, const double offset, const double xi[2], double pos[3],
          double Ar[3], double As[3], double At[3], double Arr[3], double Ars[3], double Ass[3],
          double Atr[3], double Ats[3]) const = 0;
};
 
}  // namespace b2000
 
#endif /* _B2ELEMENT_KLT_SHELL_H_ */