geom_elem.h

Go to the documentation of this file.

/*
 * MAST: Multidisciplinary-design Adaptation and Sensitivity Toolkit
 * Copyright (C) 2013-2020  Manav Bhatia and MAST authors
 *
 * This library is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public
 * License as published by the Free Software Foundation; either
 * version 2.1 of the License, or (at your option) any later version.
 *
 * This library is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public
 * License along with this library; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
 */


#ifndef __mast_geom_elem_h__
#define __mast_geom_elem_h__

// MAST includes
#include "base/mast_data_types.h"

// libMesh includes
#include "libmesh/elem.h"
#include "libmesh/fe_type.h"


namespace MAST {
    
    // Forward declerations
    class FEBase;
    class BoundaryConditionBase;
    class SystemInitialization;
    
    class GeomElem {
        
    public:
        GeomElem();
        
        virtual ~GeomElem();
        
        virtual const libMesh::Elem&
        get_reference_elem() const;

        virtual const libMesh::Elem& get_reference_local_elem() const;

        virtual const libMesh::Elem& get_quadrature_elem() const;

        virtual const libMesh::Elem& get_quadrature_local_elem() const;

        unsigned int dim() const;

        unsigned int n_sides_quadrature_elem() const;
        
        libMesh::FEType get_fe_type(unsigned int i) const;

        void set_local_y_vector(const RealVectorX& y_vec);
        
        void set_bending(bool onoff);
        
        virtual void init(const libMesh::Elem& elem,
                          const MAST::SystemInitialization& sys_init);

        virtual std::unique_ptr<MAST::FEBase>
        init_fe(bool init_grads,
                bool init_second_order_derivative,
                int extra_quadrature_order = 0) const;


        virtual std::unique_ptr<MAST::FEBase>
        init_side_fe(unsigned int s,
                     bool init_grads,
                     bool init_second_order_derivative,
                     int extra_quadrature_order = 0) const;

        
        virtual void
        external_side_loads_for_quadrature_elem
        (std::multimap<libMesh::boundary_id_type, MAST::BoundaryConditionBase*>& bc,
         std::map<unsigned int, std::vector<MAST::BoundaryConditionBase*>>& loads) const;
        
        virtual void
        get_boundary_ids_on_quadrature_elem_side
        (unsigned int s, std::vector<libMesh::boundary_id_type>& bc_ids) const;

        bool use_local_elem() const;


        const RealVectorX&
        domain_surface_normal() const;
        
        void
        transform_point_to_global_coordinate(const libMesh::Point& local_pt,
                                             libMesh::Point& global_pt) const;

        void
        transform_vector_to_global_coordinate(const libMesh::Point& local_vec,
                                              libMesh::Point& global_vec) const;

        void
        transform_vector_to_local_coordinate(const libMesh::Point& global_vec,
                                             libMesh::Point& local_vec) const;

        void
        transform_vector_to_global_coordinate(const RealVectorX& local_vec,
                                              RealVectorX& global_vec) const;
        
        void
        transform_vector_to_local_coordinate(const RealVectorX& global_vec,
                                             RealVectorX& local_vec) const;

        const RealMatrixX& T_matrix() const;
        
    protected:
        
        
        void _init_local_elem();

        void _init_local_elem_1d();

        void _init_local_elem_2d();

        const MAST::SystemInitialization  *_sys_init;

        bool                               _use_local_elem;
        
        const libMesh::Elem*               _ref_elem;
        
        libMesh::Elem*                     _local_elem;

        RealVectorX                        _local_y;

        RealVectorX                        _domain_surface_normal;

        std::vector<libMesh::Node*>        _local_nodes;

        RealMatrixX                        _T_mat;
        
        bool                                _bending = true;
    };
}

#endif  // __mast_geom_elem_h__