scale_element_hexahedral Module Reference

module FElib / Element / hexahedron More...

Data Types
type	hexahedralelement
	Derived type representing a hexahedral element. More...

Functions/Subroutines
subroutine	hexhedralelement_init (elem, elemorder_h, elemorder_v, lumpedmassmatflag)
	Initialize an object to manage a hexahedral element.

Detailed Description

module FElib / Element / hexahedron

Description

A module for a hexahedral finite element

Author

Yuta Kawai, Team SCALE

Function/Subroutine Documentation

hexhedralelement_init()

subroutine scale_element_hexahedral::hexhedralelement_init	(	class(hexahedralelement), intent(inout)	elem,
		integer, intent(in)	elemorder_h,
		integer, intent(in)	elemorder_v,
		logical, intent(in)	lumpedmassmatflag )

Initialize an object to manage a hexahedral element.

Parameters

elem	Object of finite element
elemOrder_h	Polynomial order with 1D horizontal direction
elemOrder_v	Polynomial order with vertical direction
LumpedMassMatFlag	Flag whether mass lumping is considered

Definition at line 52 of file scale_element_hexahedral.F90.

    implicit none
    
    class(HexahedralElement), intent(inout) :: elem
    integer, intent(in) :: elemOrder_h
    integer, intent(in) :: elemOrder_v
    logical, intent(in) :: LumpedMassMatFlag
 
    !-----------------------------------------------------------------------------
    
    elem%PolyOrder_h = elemorder_h
    elem%PolyOrder_v = elemorder_v
    elem%Nnode_h1D = elemorder_h + 1
    elem%Nnode_v = elemorder_v + 1
 
    elem%Nv = 8
    elem%Nfaces_h = 4
    elem%Nfaces_v = 2  
    elem%Nfaces = elem%Nfaces_h + elem%Nfaces_v
 
    elem%Nfp_h = elem%Nnode_h1D*elem%Nnode_v
    elem%Nfp_v = elem%Nnode_h1D**2
    elem%NfpTot = elem%Nfp_h*elem%Nfaces_h + elem%Nfp_v*elem%Nfaces_v
 
    elem%Np = elem%Nfp_v * elem%Nnode_v
    
    call elementbase3d_init(elem, lumpedmassmatflag)
    call construct_element(elem)
 
    return

References scale_element_base::elementbase3d_final(), scale_element_base::elementbase3d_init(), scale_element_base::elementbase_construct_liftmat(), scale_element_base::elementbase_construct_massmat(), scale_element_base::elementbase_construct_stiffmat(), scale_linalgebra::linalgebra_inv(), scale_polynomial::polynomial_gendlagrangepoly_lglpt(), scale_polynomial::polynomial_gendlegendrepoly(), scale_polynomial::polynomial_gengausslegendrept(), scale_polynomial::polynomial_gengausslegendreptintweight(), scale_polynomial::polynomial_gengausslobattopt(), scale_polynomial::polynomial_gengausslobattoptintweight(), scale_polynomial::polynomial_genlagrangepoly(), and scale_polynomial::polynomial_genlegendrepoly().