#include "users_modules/mofem_um_homogenisation/src/BCs_RVELagrange_Periodic.hpp"

Inheritance diagram for BCs_RVELagrange_Periodic:

Collaboration diagram for BCs_RVELagrange_Periodic:

Classes
struct	CommonDataPeriodic

struct	CommonFunctionsPeriodic

struct	MyPrismFE

struct	OpDmatRhs
	\biref operator to calculate the RHS of the constrain for the RVE boundary conditions More...

struct	OpRVEBCsPeriodicCalAssemCmat
	\biref operator to calculate and assemble Cmat More...

struct	OpRVEBCsPeriodicCalCTLam

struct	OpRVEBCsPeriodicCalCU

struct	OpRVEBCsPeriodicCalDispAtGaussPts

struct	OpRVEBCsPeriodicCalLagMulAtGaussPts

struct	OpRVEBCsPeriodicColInd

struct	OpRVEBCsPeriodicRhs

struct	OpRVEBCsPeriodicRhs_givenStrain

struct	OpRVEBCsPeriodicRowInd

struct	OpRVEHomoStress
	\biref operator to calculate the RVE homogenised stress More...

struct	RVEBC_Data_Periodic

Public Member Functions
MyPrismFE &	getLoopFeRVEBCLhs ()

MyPrismFE &	getLoopFeRVEBCRhs ()

MyPrismFE &	getLoopFeRVEBCRhsResidual ()

MyPrismFE &	getLoopFeRVEBCStress ()

	BCs_RVELagrange_Periodic (MoFEM::Interface &m_field)

PetscErrorCode	addLagrangiangElement (const string element_name, const string field_name, const string lagrang_field_name, const string mesh_nodals_positions, Range &periodic_prisms)

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 > &f)

PetscErrorCode	setRVEBCsHomoStressOperators (string field_name, string lagrang_field_name, string mesh_nodals_positions, Vec stress_homo)

Public Member Functions inherited from BCs_RVELagrange_Disp
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

map< int, RVEBC_Data_Periodic >	setOfRVEBCPrisms
	maps side set id with appropriate FluxData

MyPrismFE	feRVEBCRhs

MyPrismFE	feRVEBCLhs

MyPrismFE	feRVEBCStress

MyPrismFE	feRVEBCRhsResidual

CommonFunctionsPeriodic	commonFunctionsPeriodic

CommonDataPeriodic	commonDataPeriodic

Public Attributes inherited from BCs_RVELagrange_Disp
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 19 of file BCs_RVELagrange_Periodic.hpp.

Constructor & Destructor Documentation

◆ BCs_RVELagrange_Periodic()

BCs_RVELagrange_Periodic::BCs_RVELagrange_Periodic ( MoFEM::Interface & m_field )

inline

Definition at line 44 of file BCs_RVELagrange_Periodic.hpp.

                                                   :
  BCs_RVELagrange_Disp(m_field),
  feRVEBCRhs(m_field),
  feRVEBCLhs(m_field),
  feRVEBCRhsResidual(m_field),
  feRVEBCStress(m_field) {
  }

Member Function Documentation

◆ addLagrangiangElement()

PetscErrorCode BCs_RVELagrange_Periodic::addLagrangiangElement	(	const string	element_name,
		const string	field_name,
		const string	lagrang_field_name,
		const string	mesh_nodals_positions,
		Range &	periodic_prisms
	)

inline

Definition at line 52 of file BCs_RVELagrange_Periodic.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);
    }
 
    //Adding Prisims to Element Lagrange_elem (to loop over these prisims)
    setOfRVEBCPrisms[1].pRisms = periodic_prisms;
    ierr = mField.add_ents_to_finite_element_by_type(periodic_prisms,MBPRISM,element_name); CHKERRQ(ierr);
 
    
    Range tris;
    rval = mField.get_moab().get_adjacencies(periodic_prisms,2,false,tris,moab::Interface::UNION); CHKERRQ_MOAB(rval);
    setOfRVEBC[1].tRis = tris.subset_by_type(MBTRI);
 
    PetscFunctionReturn(0);
  }

◆ getLoopFeRVEBCLhs()

MyPrismFE & BCs_RVELagrange_Periodic::getLoopFeRVEBCLhs ( )

inline

Definition at line 38 of file BCs_RVELagrange_Periodic.hpp.

38{ return feRVEBCLhs; }

◆ getLoopFeRVEBCRhs()

MyPrismFE & BCs_RVELagrange_Periodic::getLoopFeRVEBCRhs ( )

inline

Definition at line 39 of file BCs_RVELagrange_Periodic.hpp.

39{ return feRVEBCRhs; }

◆ getLoopFeRVEBCRhsResidual()

MyPrismFE & BCs_RVELagrange_Periodic::getLoopFeRVEBCRhsResidual ( )

inline

Definition at line 41 of file BCs_RVELagrange_Periodic.hpp.

41{ return feRVEBCRhsResidual; }

◆ getLoopFeRVEBCStress()

MyPrismFE & BCs_RVELagrange_Periodic::getLoopFeRVEBCStress ( )

inline

Definition at line 42 of file BCs_RVELagrange_Periodic.hpp.

42{ return feRVEBCStress; }

◆ setRVEBCsHomoStressOperators()

PetscErrorCode BCs_RVELagrange_Periodic::setRVEBCsHomoStressOperators	(	string	field_name,
		string	lagrang_field_name,
		string	mesh_nodals_positions,
		Vec	stress_homo
	)

inline

Definition at line 1312 of file BCs_RVELagrange_Periodic.hpp.

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

◆ setRVEBCsOperators()

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

inline

Definition at line 1176 of file BCs_RVELagrange_Periodic.hpp.

    {
    PetscFunctionBegin;
 
    bool ho_geometry = false;
    if(mField.check_field(mesh_nodals_positions)) {
      ho_geometry = true;
    }
    //      cout<<"Hi 1 from setRVEBCsOperators periodic "<<endl;
    map<int,RVEBC_Data_Periodic>::iterator sit = setOfRVEBCPrisms.begin();
    for(;sit!=setOfRVEBCPrisms.end();sit++) {
 
      //        LHS
      //Col indices and Col N
      feRVEBCLhs.getOpPtrVector().push_back(
        new OpRVEBCsPeriodicColInd(
          field_name,sit->second,commonDataPeriodic,commonFunctionsPeriodic,ho_geometry
        )
      );
 
 
      //Row indices and Row N
      feRVEBCLhs.getOpPtrVector().push_back(
        new OpRVEBCsPeriodicRowInd(
          lagrang_field_name,sit->second,commonDataPeriodic,commonFunctionsPeriodic,ho_geometry
        )
      );
 
      //Calculate and assemble Cmat
      feRVEBCLhs.getOpPtrVector().push_back(
        new OpRVEBCsPeriodicCalAssemCmat(
          field_name,lagrang_field_name,aij,sit->second,commonFunctionsPeriodic, commonDataPeriodic,ho_geometry
        )
      );
 
 
 
      //RHS
      feRVEBCRhs.getOpPtrVector().push_back(
        new OpRVEBCsPeriodicRhs(
          field_name,lagrang_field_name,f,sit->second,commonDataPeriodic,commonFunctionsPeriodic,ho_geometry
        )
      );
 
    }
    PetscFunctionReturn(0);
  }

◆ setRVEBCsOperatorsNonlinear()

PetscErrorCode BCs_RVELagrange_Periodic::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 1091 of file BCs_RVELagrange_Periodic.hpp.

    {
    PetscFunctionBegin;
 
    bool ho_geometry = false;
    if(mField.check_field(mesh_nodals_positions)) {
      ho_geometry = true;
    }
 
    map<int,RVEBC_Data_Periodic>::iterator sit = setOfRVEBCPrisms.begin();
    for(;sit!=setOfRVEBCPrisms.end();sit++) {
 
      //Col indices and Col N
      feRVEBCLhs.getOpPtrVector().push_back(
        new OpRVEBCsPeriodicColInd(
          field_name,sit->second,commonDataPeriodic,commonFunctionsPeriodic,ho_geometry
        )
      );
 
 
      //Row indices and Row N
      feRVEBCLhs.getOpPtrVector().push_back(
        new OpRVEBCsPeriodicRowInd(
          lagrang_field_name,sit->second,commonDataPeriodic,commonFunctionsPeriodic,ho_geometry
        )
      );
 
      //Calculate and assemble Cmat
      feRVEBCLhs.getOpPtrVector().push_back(
        new OpRVEBCsPeriodicCalAssemCmat(
          field_name,lagrang_field_name,aij,sit->second,commonFunctionsPeriodic, commonDataPeriodic,ho_geometry
        )
      );
 
      //RHS
      feRVEBCRhs.getOpPtrVector().push_back(
        new OpRVEBCsPeriodicRhs_givenStrain(
          field_name,lagrang_field_name,f,given_strain, methodsOp, sit->second,commonDataPeriodic,commonFunctionsPeriodic,ho_geometry
        )
      );
 
 
      //Calculate displacement at gauss point
      feRVEBCRhsResidual.getOpPtrVector().push_back(
        new OpRVEBCsPeriodicCalDispAtGaussPts(
          field_name,sit->second,commonDataPeriodic,commonFunctionsPeriodic,ho_geometry
        )
      );
 
      //Calculate and assemble CU = NT * u
      feRVEBCRhsResidual.getOpPtrVector().push_back(
        new OpRVEBCsPeriodicCalCU(
          lagrang_field_name,sit->second,f,commonDataPeriodic,commonFunctionsPeriodic,ho_geometry
        )
      );
 
 
      //Calculate lagrange multipliers at gauss point
      feRVEBCRhsResidual.getOpPtrVector().push_back(
        new OpRVEBCsPeriodicCalLagMulAtGaussPts(
          lagrang_field_name,sit->second,commonDataPeriodic,commonFunctionsPeriodic,ho_geometry
        )
      );
 
      //Calculate and assemble CT Lam = NT * Lam
      feRVEBCRhsResidual.getOpPtrVector().push_back(
        new OpRVEBCsPeriodicCalCTLam(
          field_name,sit->second,f,commonDataPeriodic,commonFunctionsPeriodic,ho_geometry
        )
      );
 
//      cout<<"Hi 1 from setRVEBCsOperators periodic "<<endl;
//      string aaa;
//      cin>>aaa;
    }
    PetscFunctionReturn(0);
  }