#include "users_modules/mofem_um_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

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.

35{ return feRVEBCLhs; }

◆ getLoopFeRVEBCRhs()

MyTriFE & BCs_RVELagrange_Disp::getLoopFeRVEBCRhs ( )

inline

Definition at line 36 of file BCs_RVELagrange_Disp.hpp.

36{ return feRVEBCRhs; }

◆ getLoopFeRVEBCRhsHomoC()

MyTriFE & BCs_RVELagrange_Disp::getLoopFeRVEBCRhsHomoC ( )

inline

Definition at line 41 of file BCs_RVELagrange_Disp.hpp.

41{ return feRVEBCRhsHomoC; }

◆ getLoopFeRVEBCRhsResidual()

MyTriFE & BCs_RVELagrange_Disp::getLoopFeRVEBCRhsResidual ( )

inline

Definition at line 38 of file BCs_RVELagrange_Disp.hpp.

38{ return feRVEBCRhsResidual; }

◆ getLoopFeRVEBCStress()

MyTriFE & BCs_RVELagrange_Disp::getLoopFeRVEBCStress ( )

inline

Definition at line 39 of file BCs_RVELagrange_Disp.hpp.

39{ 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 640 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 556 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 511 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 416 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/mofem_um_homogenisation/src/BCs_RVELagrange_Disp.hpp

Classes

Public Member Functions

Public Attributes

Detailed Description

Constructor & Destructor Documentation

◆ BCs_RVELagrange_Disp()

Member Function Documentation

◆ addLagrangiangElement()

◆ getLoopFeRVEBCLhs()

◆ getLoopFeRVEBCRhs()

◆ getLoopFeRVEBCRhsHomoC()

◆ getLoopFeRVEBCRhsResidual()

◆ getLoopFeRVEBCStress()

◆ setRVEBCsHomoStressOperators()

◆ setRVEBCsOperators()

◆ setRVEBCsOperatorsNonlinear()

Member Data Documentation

◆ common_functions

◆ feRVEBCLhs

◆ feRVEBCRhs

◆ feRVEBCRhsHomoC

◆ feRVEBCRhsResidual

◆ feRVEBCStress

◆ methodsOp

◆ mField

◆ setOfRVEBC