BCs_RVELagrange_Disp Struct Reference

#include <users_modules/homogenisation/src/BCs_RVELagrange_Disp.hpp>

Inheritance diagram for BCs_RVELagrange_Disp:

Collaboration diagram for BCs_RVELagrange_Disp:

Classes
struct	CommonFunctions
struct	MyTriFE
struct	OpDmatRhs
struct	OpRVEBCsLhs
	\biref operator to calculate the LHS for the RVE bounary conditions More...
struct	OpRVEBCsRhs
	\biref operator to calculate the RHS of the constrain for the RVE boundary conditions More...
struct	OpRVEBCsRhsForceExternal_CU
	\biref operator to calculate and assemble CU More...
struct	OpRVEBCsRhsHomoC
	\biref operator to calculate the RHS for the calculation of the homgoensied stiffness matrix More...
struct	OpRVEHomoStress
	\biref operator to calculate the RVE homogenised stress More...
struct	RVEBC_Data

Public Member Functions
MyTriFE &	getLoopFeRVEBCLhs ()
MyTriFE &	getLoopFeRVEBCRhs ()
MyTriFE &	getLoopFeRVEBCRhsResidual ()
MyTriFE &	getLoopFeRVEBCStress ()
MyTriFE &	getLoopFeRVEBCRhsHomoC ()
	BCs_RVELagrange_Disp (MoFEM::Interface &m_field)
PetscErrorCode	addLagrangiangElement (const string element_name, const string field_name, const string lagrang_field_name, const string mesh_nodals_positions)
PetscErrorCode	setRVEBCsOperatorsNonlinear (string field_name, string lagrang_field_name, string mesh_nodals_positions, Mat aij, vector< Vec > &fvec, Vec f, VectorDouble given_strain)
PetscErrorCode	setRVEBCsOperators (string field_name, string lagrang_field_name, string mesh_nodals_positions, Mat aij, vector< Vec > &fvec)
PetscErrorCode	setRVEBCsHomoStressOperators (string field_name, string lagrang_field_name, string mesh_nodals_positions, Vec Stress_Homo)

Public Attributes
boost::ptr_vector< MethodForForceScaling >	methodsOp
MoFEM::Interface &	mField
MyTriFE	feRVEBCRhs
MyTriFE	feRVEBCLhs
MyTriFE	feRVEBCStress
MyTriFE	feRVEBCRhsResidual
MyTriFE	feRVEBCRhsHomoC
map< int, RVEBC_Data >	setOfRVEBC
	maps side set id with appropriate FluxData More...
CommonFunctions	common_functions

Detailed Description

Definition at line 18 of file BCs_RVELagrange_Disp.hpp.

Constructor & Destructor Documentation

BCs_RVELagrange_Disp()

BCs_RVELagrange_Disp::BCs_RVELagrange_Disp ( MoFEM::Interface &  m_field )

inline

Definition at line 44 of file BCs_RVELagrange_Disp.hpp.

                                               :
  mField(m_field),
  feRVEBCRhs(m_field),
  feRVEBCLhs(m_field),
  feRVEBCRhsResidual(m_field),
  feRVEBCStress(m_field),
  feRVEBCRhsHomoC(m_field){
  }

Member Function Documentation

addLagrangiangElement()

PetscErrorCode BCs_RVELagrange_Disp::addLagrangiangElement ( const string  element_name, const string  field_name, const string  lagrang_field_name, const string  mesh_nodals_positions  )

inline

Definition at line 58 of file BCs_RVELagrange_Disp.hpp.

    {
    PetscFunctionBegin;
 
    
    
 
    ierr = mField.add_finite_element(element_name,MF_ZERO); CHKERRQ(ierr);
 
    //============================================================================================================
    //C row as Lagrange_mul_disp and col as DISPLACEMENT
    ierr = mField.modify_finite_element_add_field_row(element_name,lagrang_field_name); CHKERRQ(ierr);
    ierr = mField.modify_finite_element_add_field_col(element_name,field_name); CHKERRQ(ierr);
    //CT col as Lagrange_mul_disp and row as DISPLACEMENT
    ierr = mField.modify_finite_element_add_field_col(element_name,lagrang_field_name); CHKERRQ(ierr);
    ierr = mField.modify_finite_element_add_field_row(element_name,field_name); CHKERRQ(ierr);
    //data
    ierr = mField.modify_finite_element_add_field_data(element_name,lagrang_field_name); CHKERRQ(ierr);
    ierr = mField.modify_finite_element_add_field_data(element_name,field_name); CHKERRQ(ierr);
    //============================================================================================================
 
    if(mField.check_field(mesh_nodals_positions)) { //for high order geometry
      ierr = mField.modify_finite_element_add_field_data(element_name,mesh_nodals_positions); CHKERRQ(ierr);
    }
 
    //this is alternative method for setting boundary conditions, to bypass bu in cubit file reader.
    //not elegant, but good enough
    for(_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(mField,SIDESET,it)) {
      if(it->getName().compare(0,12,"AllBoundSurf") == 0 || it->getMeshsetId() == 103) {
        //          cout<<"Hi from addLagrangiangElement "<<endl;
        rval = mField.get_moab().get_entities_by_type(it->meshset,MBTRI,setOfRVEBC[it->getMeshsetId()].tRis,true); CHKERRQ_MOAB(rval);
        ierr = mField.add_ents_to_finite_element_by_type(setOfRVEBC[it->getMeshsetId()].tRis,MBTRI,element_name); CHKERRQ(ierr);
      }
    }
 
    PetscFunctionReturn(0);
  }

getLoopFeRVEBCLhs()

MyTriFE& BCs_RVELagrange_Disp::getLoopFeRVEBCLhs ( )

inline

Definition at line 35 of file BCs_RVELagrange_Disp.hpp.

{ return feRVEBCLhs; }

getLoopFeRVEBCRhs()

MyTriFE& BCs_RVELagrange_Disp::getLoopFeRVEBCRhs ( )

inline

Definition at line 36 of file BCs_RVELagrange_Disp.hpp.

{ return feRVEBCRhs; }

getLoopFeRVEBCRhsHomoC()

MyTriFE& BCs_RVELagrange_Disp::getLoopFeRVEBCRhsHomoC ( )

inline

Definition at line 41 of file BCs_RVELagrange_Disp.hpp.

{ return feRVEBCRhsHomoC; }

getLoopFeRVEBCRhsResidual()

MyTriFE& BCs_RVELagrange_Disp::getLoopFeRVEBCRhsResidual ( )

inline

Definition at line 38 of file BCs_RVELagrange_Disp.hpp.

{ return feRVEBCRhsResidual; }

getLoopFeRVEBCStress()

MyTriFE& BCs_RVELagrange_Disp::getLoopFeRVEBCStress ( )

inline

Definition at line 39 of file BCs_RVELagrange_Disp.hpp.

{ return feRVEBCStress; }

setRVEBCsHomoStressOperators()

PetscErrorCode BCs_RVELagrange_Disp::setRVEBCsHomoStressOperators ( string  field_name, string  lagrang_field_name, string  mesh_nodals_positions, Vec  Stress_Homo  )

inline

Definition at line 671 of file BCs_RVELagrange_Disp.hpp.

                                                {
    PetscFunctionBegin;
    
    bool ho_geometry = false;
    if(mField.check_field(mesh_nodals_positions)) {
      ho_geometry = true;
    }
    
    map<int,RVEBC_Data>::iterator sit = setOfRVEBC.begin();
    for(;sit!=setOfRVEBC.end();sit++) {
      feRVEBCStress.getOpPtrVector().push_back(new OpRVEHomoStress(field_name, lagrang_field_name, Stress_Homo, sit->second, ho_geometry));
    }
    
    PetscFunctionReturn(0);
  }

setRVEBCsOperators()

PetscErrorCode BCs_RVELagrange_Disp::setRVEBCsOperators ( string  field_name, string  lagrang_field_name, string  mesh_nodals_positions, Mat  aij, vector< Vec > &  fvec  )

inline

Definition at line 580 of file BCs_RVELagrange_Disp.hpp.

        {
    PetscFunctionBegin;
    
    bool ho_geometry = false;
    if(mField.check_field(mesh_nodals_positions)) {
      ho_geometry = true;
    }
    map<int,RVEBC_Data>::iterator sit = setOfRVEBC.begin();
    for(;sit!=setOfRVEBC.end();sit++) {
      
      //Lhs (calculate LHS)
      feRVEBCLhs.getOpPtrVector().push_back(
      new OpRVEBCsLhs(field_name,lagrang_field_name, aij, sit->second,ho_geometry)
      );
 
      //Rhs (F[6] used to calculate Homgenised stiffness matrix)
      feRVEBCRhs.getOpPtrVector().push_back(
       new OpRVEBCsRhsHomoC(field_name,lagrang_field_name,fvec,sit->second,ho_geometry)
       );
     }
    PetscFunctionReturn(0);
  }

setRVEBCsOperatorsNonlinear()

PetscErrorCode BCs_RVELagrange_Disp::setRVEBCsOperatorsNonlinear ( string  field_name, string  lagrang_field_name, string  mesh_nodals_positions, Mat  aij, vector< Vec > &  fvec, Vec  f, VectorDouble  given_strain  )

inline

Definition at line 535 of file BCs_RVELagrange_Disp.hpp.

    {
    PetscFunctionBegin;
 
    bool ho_geometry = false;
    if(mField.check_field(mesh_nodals_positions)) {
      ho_geometry = true;
    }
    map<int,RVEBC_Data>::iterator sit = setOfRVEBC.begin();
    for(;sit!=setOfRVEBC.end();sit++) {
      
      //Lhs (calculate LHS)
      feRVEBCLhs.getOpPtrVector().push_back(
        new OpRVEBCsLhs(field_name,lagrang_field_name, aij, sit->second,ho_geometry)
      );
      
      
      //Rhs (calculate RHS)
      feRVEBCRhs.getOpPtrVector().push_back(
        new OpRVEBCsRhs(field_name,lagrang_field_name,f,given_strain, methodsOp, sit->second,ho_geometry)
      );
      
      
      //Rhs (F[6] used to calculate Homgenised stiffness matrix)
      feRVEBCRhsHomoC.getOpPtrVector().push_back(
         new OpRVEBCsRhsHomoC(field_name,lagrang_field_name,fvec,sit->second,ho_geometry)
       );
 
      //for internal forces (i.e. RHS)
      //calculate and assemble CU and CT Lam
      feRVEBCRhsResidual.getOpPtrVector().push_back(
        new OpRVEBCsRhsForceExternal_CU(field_name,lagrang_field_name,f,common_functions,sit->second,ho_geometry)
      );
      
      
     }
 
    PetscFunctionReturn(0);
  }

Member Data Documentation

common_functions

CommonFunctions BCs_RVELagrange_Disp::common_functions

Definition at line 434 of file BCs_RVELagrange_Disp.hpp.

feRVEBCLhs

MyTriFE BCs_RVELagrange_Disp::feRVEBCLhs

Definition at line 30 of file BCs_RVELagrange_Disp.hpp.

feRVEBCRhs

MyTriFE BCs_RVELagrange_Disp::feRVEBCRhs

Definition at line 29 of file BCs_RVELagrange_Disp.hpp.

feRVEBCRhsHomoC

MyTriFE BCs_RVELagrange_Disp::feRVEBCRhsHomoC

Definition at line 33 of file BCs_RVELagrange_Disp.hpp.

feRVEBCRhsResidual

MyTriFE BCs_RVELagrange_Disp::feRVEBCRhsResidual

Definition at line 32 of file BCs_RVELagrange_Disp.hpp.

feRVEBCStress

MyTriFE BCs_RVELagrange_Disp::feRVEBCStress

Definition at line 31 of file BCs_RVELagrange_Disp.hpp.

methodsOp

boost::ptr_vector<MethodForForceScaling> BCs_RVELagrange_Disp::methodsOp

Definition at line 25 of file BCs_RVELagrange_Disp.hpp.

mField

MoFEM::Interface& BCs_RVELagrange_Disp::mField

Definition at line 28 of file BCs_RVELagrange_Disp.hpp.

setOfRVEBC

map<int,RVEBC_Data> BCs_RVELagrange_Disp::setOfRVEBC

maps side set id with appropriate FluxData

Definition at line 56 of file BCs_RVELagrange_Disp.hpp.

---

The documentation for this struct was generated from the following file:

• users_modules/homogenisation/src/BCs_RVELagrange_Disp.hpp