isotropic_element_property_card_3D.h

Go to the documentation of this file.

/*
 * MAST: Multidisciplinary-design Adaptation and Sensitivity Toolkit
 * Copyright (C) 2013-2019  Manav Bhatia
 *
 * 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__isotropic_element_property_card_3D__
#define __mast__isotropic_element_property_card_3D__

// MAST includes
#include "property_cards/element_property_card_base.h"



namespace MAST
{
    
    class IsotropicElementPropertyCard3D:
    public MAST::ElementPropertyCardBase {
        
    public:
        IsotropicElementPropertyCard3D():
        MAST::ElementPropertyCardBase(),
        _material(nullptr)
        { }
        
        virtual ~IsotropicElementPropertyCard3D() { }
        
        
        virtual unsigned int dim() const {
            return 3;
        }

        virtual bool if_isotropic() const {
            return true;
        }

        
        virtual void set_material(MAST::MaterialPropertyCardBase& mat) {
            _material = &mat;
        }
        
        
        const MAST::MaterialPropertyCardBase& get_material() const {
            libmesh_assert(_material); // make sure it has already been set
            return *_material;
        }
        

        virtual MAST::BendingOperatorType
        bending_model(const MAST::GeomElem& elem) const {
            libmesh_assert(false);
        }
        
        virtual int extra_quadrature_order(const MAST::GeomElem& elem) const {
            libmesh_assert(false);
            return 0;
        }

        virtual bool depends_on(const MAST::FunctionBase& f) const;

        
        virtual std::unique_ptr<MAST::FieldFunction<RealMatrixX> >
        stiffness_A_matrix(const MAST::ElementBase& e) const;
        
        virtual std::unique_ptr<MAST::FieldFunction<RealMatrixX> >
        stiffness_B_matrix(const MAST::ElementBase& e) const;
        
        virtual std::unique_ptr<MAST::FieldFunction<RealMatrixX> >
        stiffness_D_matrix(const MAST::ElementBase& e) const;
        
        virtual std::unique_ptr<MAST::FieldFunction<RealMatrixX> >
        damping_matrix(const MAST::ElementBase& e) const;
        
        virtual std::unique_ptr<MAST::FieldFunction<RealMatrixX> >
        inertia_matrix(const MAST::ElementBase& e) const;
        
        virtual std::unique_ptr<MAST::FieldFunction<RealMatrixX> >
        thermal_expansion_A_matrix(const MAST::ElementBase& e) const;
        
        virtual std::unique_ptr<MAST::FieldFunction<RealMatrixX> >
        thermal_expansion_B_matrix(const MAST::ElementBase& e) const;
        
        virtual std::unique_ptr<MAST::FieldFunction<RealMatrixX> >
        transverse_shear_stiffness_matrix(const MAST::ElementBase& e) const;
        
        virtual std::unique_ptr<MAST::FieldFunction<RealMatrixX> >
        prestress_A_matrix( MAST::ElementBase& e) const;
        
        virtual std::unique_ptr<MAST::FieldFunction<RealMatrixX> >
        prestress_B_matrix( MAST::ElementBase& e) const;
        
        virtual std::unique_ptr<MAST::FieldFunction<RealMatrixX> >
        thermal_conductance_matrix(const MAST::ElementBase& e) const;
        
        virtual std::unique_ptr<MAST::FieldFunction<RealMatrixX> >
        thermal_capacitance_matrix(const MAST::ElementBase& e) const;
        
        virtual const MAST::FieldFunction<Real>&
        section(const MAST::ElementBase& e) const {
            libmesh_error();
            MAST::FieldFunction<Real>* ptr = nullptr;
            return *ptr;
        }

    protected:
        
        MAST::MaterialPropertyCardBase* _material;
    };
}


#endif // __mast__isotropic_element_property_card_3D__