element_property_card_base.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__element_property_card_base__
#define __mast__element_property_card_base__


// MAST includes
#include "base/function_set_base.h"
#include "elasticity/bending_operator.h"


namespace MAST
{
    // Forward decleration
    class MaterialPropertyCardBase;
    class ElementBase;
    class GeomElem;
    template <typename ValType> class FieldFunction;
    
    
    enum StrainType {
        LINEAR_STRAIN,
        NONLINEAR_STRAIN
    };
    
    
    
    class ElementPropertyCardBase:
    public MAST::FunctionSetBase {
        
    public:
        ElementPropertyCardBase():
        MAST::FunctionSetBase(),
        _strain_type(MAST::LINEAR_STRAIN),
        _diagonal_mass(false)
        { }
        
        virtual ~ElementPropertyCardBase() { }
        
        virtual MAST::BendingOperatorType
        bending_model(const MAST::GeomElem& elem) const = 0;
        
        virtual int extra_quadrature_order(const MAST::GeomElem& elem) const = 0;
        

        virtual std::unique_ptr<MAST::FieldFunction<RealMatrixX> >
        stiffness_A_matrix(const MAST::ElementBase& e) const = 0;

        virtual std::unique_ptr<MAST::FieldFunction<RealMatrixX> >
        stiffness_B_matrix(const MAST::ElementBase& e) const = 0;

        virtual std::unique_ptr<MAST::FieldFunction<RealMatrixX> >
        stiffness_D_matrix(const MAST::ElementBase& e) const = 0;

        virtual std::unique_ptr<MAST::FieldFunction<RealMatrixX> >
        damping_matrix(const MAST::ElementBase& e) const = 0;

        virtual std::unique_ptr<MAST::FieldFunction<RealMatrixX> >
        inertia_matrix(const MAST::ElementBase& e) const = 0;

        virtual std::unique_ptr<MAST::FieldFunction<RealMatrixX> >
        thermal_expansion_A_matrix(const MAST::ElementBase& e) const = 0;

        virtual std::unique_ptr<MAST::FieldFunction<RealMatrixX> >
        thermal_expansion_B_matrix(const MAST::ElementBase& e) const = 0;

        virtual std::unique_ptr<MAST::FieldFunction<RealMatrixX> >
        transverse_shear_stiffness_matrix(const MAST::ElementBase& e) const = 0;

        virtual std::unique_ptr<MAST::FieldFunction<RealMatrixX> >
        prestress_A_matrix( MAST::ElementBase& e) const = 0;

        virtual std::unique_ptr<MAST::FieldFunction<RealMatrixX> >
        prestress_B_matrix( MAST::ElementBase& e) const = 0;
        
        virtual std::unique_ptr<MAST::FieldFunction<RealMatrixX> >
        thermal_conductance_matrix(const MAST::ElementBase& e) const = 0;

        virtual std::unique_ptr<MAST::FieldFunction<RealMatrixX> >
        thermal_capacitance_matrix(const MAST::ElementBase& e) const= 0;

        virtual const MAST::FieldFunction<Real>*
        section(const MAST::ElementBase& e) const = 0;

        
        virtual bool if_isotropic() const = 0;
        
        
        virtual const MAST::MaterialPropertyCardBase& get_material() const {
            libmesh_error();
        }
        
        
        virtual const MAST::MaterialPropertyCardBase& set_material(MAST::MaterialPropertyCardBase& mat) const {
            libmesh_error_msg("Not Implemented, this needs to be reimplemented for individual card type; In " << __PRETTY_FUNCTION__ << " in " << __FILE__ << " at line " << __LINE__);
        }
        
        
        virtual RealVectorX& y_vector() {
            libmesh_error_msg("Not Implemented, this needs to be reimplemented for individual card type; In " << __PRETTY_FUNCTION__ << " in " << __FILE__ << " at line " << __LINE__);
        }
        

        virtual const RealVectorX& y_vector() const {
            libmesh_error_msg("Not Implemented, this needs to be reimplemented for individual card type; In " << __PRETTY_FUNCTION__ << " in " << __FILE__ << " at line " << __LINE__);
        }
        
        virtual void init() {
            libmesh_error_msg("Not Implemented, this needs to be reimplemented for individual card type; In " << __PRETTY_FUNCTION__ << " in " << __FILE__ << " at line " << __LINE__);
        }
        
        
        virtual unsigned int dim() const = 0;
        
        
        void set_strain(MAST::StrainType strain) {
            _strain_type = strain;
        }
        
        
        const MAST::StrainType strain_type() const {
            return _strain_type;
        }
        
        
        virtual void set_bending_model(MAST::BendingOperatorType b)  {
            libmesh_error_msg("Not implemented, this needs to be reimplemented for individual card type; In " << __PRETTY_FUNCTION__ << " in " << __FILE__ << " at line " << __LINE__);
        }
        
        
        void set_diagonal_mass_matrix(bool m) {
            _diagonal_mass = m;
        }
        
        
        bool if_diagonal_mass_matrix() const {
            return _diagonal_mass;
        }
        
        
        virtual bool if_prestressed() const {
            return this->contains("prestress");
        }
        
        
        virtual bool get_warping_only() const
        {
            return _warping_only;
        }
        
    protected:
        
        MAST::StrainType _strain_type;
        
        bool _diagonal_mass;
        
        bool _warping_only = false;
    };
    
}



#endif // __mast__element_property_card_base__