#include <transportgradientdirichlet.h>

Inheritance diagram for oofem::TransportGradientDirichlet:

Collaboration diagram for oofem::TransportGradientDirichlet:

Public Member Functions
	TransportGradientDirichlet (int n, Domain *d)
void	initializeFrom (InputRecord &ir) override
void	giveInputRecord (DynamicInputRecord &input) override
void	postInitialize () override
	Performs post initialization steps. Called after all components are created and initialized.
double	give (Dof *dof, ValueModeType mode, double time) override
bcType	giveType () const override
double	domainSize ()
void	computeCoefficientMatrix (FloatMatrix &C)
virtual void	computeField (FloatArray &sigma, TimeStep *tStep)
virtual void	computeTangent (FloatMatrix &tangent, TimeStep *tStep)
void	computeXi ()
	Computes the offset values for "improved" Dirichlet. See class description.
void	scale (double s) override
const char *	giveClassName () const override
const char *	giveInputRecordName () const override
Public Member Functions inherited from oofem::BoundaryCondition
	BoundaryCondition (int i, Domain *d)
virtual	~BoundaryCondition ()
	Destructor.
virtual double	give (Dof dof, ValueModeType mode, TimeStep tStep)
void	setPrescribedValue (double s)
void	saveContext (DataStream &stream, ContextMode mode) override
void	restoreContext (DataStream &stream, ContextMode mode) override
Public Member Functions inherited from oofem::GeneralBoundaryCondition
	GeneralBoundaryCondition (int n, Domain *d)
virtual	~GeneralBoundaryCondition ()
	Destructor.
int	giveSetNumber () const
virtual int	giveNumberOfInternalDofManagers ()
	Gives the number of internal dof managers.
virtual DofManager *	giveInternalDofManager (int i)
	Gives an internal dof manager from receiver.
Function *	giveTimeFunction ()
int	getIsImposedTimeFunctionNumber () const
void	setIsImposedTimeFunctionNumber (int funcIndx)
virtual bcValType	giveBCValType () const
virtual bool	isImposed (TimeStep *tStep)
virtual const IntArray &	giveDofIDs () const
virtual bcGeomType	giveBCGeoType () const
virtual double	giveProperty (int aProperty, TimeStep *tStep) const
void	saveContext (DataStream &stream, ContextMode mode) override
void	restoreContext (DataStream &stream, ContextMode mode) override
virtual void	updateYourself (TimeStep *tStep)
Public Member Functions inherited from oofem::FEMComponent
	FEMComponent (int n, Domain *d)
virtual	~FEMComponent ()=default
	Virtual destructor.
Domain *	giveDomain () const
virtual void	setDomain (Domain *d)
int	giveNumber () const
void	setNumber (int num)
virtual void	updateLocalNumbering (EntityRenumberingFunctor &f)
virtual void	initializeFrom (InputRecord &ir, int priority)
virtual void	initializeFinish ()
virtual int	checkConsistency ()
virtual void	printOutputAt (FILE file, TimeStep tStep)
virtual void	printYourself ()
	Prints receiver state on stdout. Useful for debugging.
virtual Interface *	giveInterface (InterfaceType t)
std::string	errorInfo (const char *func) const
	Returns string for prepending output (used by error reporting macros).

Protected Attributes
FloatArray	mGradient
FloatArray	mCenterCoord
bool	tractionControl = false
std ::map< int, FloatArray >	xis
	Stores one "psi" value for each node.
IntArray	surfSets
Protected Attributes inherited from oofem::BoundaryCondition
FloatArray	values
	Prescribed values for each resp. dof.
Protected Attributes inherited from oofem::GeneralBoundaryCondition
int	timeFunction
	Associated load time function.
bcValType	valType
	Physical meaning of BC value.
IntArray	dofs
	Dofs that b.c. is applied to (relevant for Dirichlet type b.c.s).
int	isImposedTimeFunction
int	set
	Set number for boundary condition to be applied to.
Protected Attributes inherited from oofem::FEMComponent
int	number
	Component number.
Domain *	domain
	Link to domain object, useful for communicating with other FEM components.

Detailed Description

Prescribes \( T = g_{i}(x_i-\bar{x}_i) \) where \( T \) are primary unknown. This is typical boundary condition in multiscale analysis where \( g = \partial_x T \) would be a macroscopic gradient at the integration point, i.e. this is a boundary condition for prolongation. It is also convenient to use when one wants to test a arbitrary specimen for a given average gradient.

Note: This class has optional support for the "psi" parameter, which is computed to improve the bounds for high permeable inclusions in a matrix. Using this parameter requires additional inputs, edgeSets and surfSets, which should be defined in the particular order. The numbering of the surface should be the same as for the hex cube interpolator as defined in the element reference manual.; Experimental code

Author: Mikael Öhman

Definition at line 73 of file transportgradientdirichlet.h.

Constructor & Destructor Documentation

◆ TransportGradientDirichlet()

oofem::TransportGradientDirichlet::TransportGradientDirichlet	(	int	n,
		Domain *	d )

inline

Definition at line 85 of file transportgradientdirichlet.h.

References oofem::BoundaryCondition::BoundaryCondition().

Member Function Documentation

◆ computeCoefficientMatrix()

void oofem::TransportGradientDirichlet::computeCoefficientMatrix ( FloatMatrix & C )

Constructs a coefficient matrix for all prescribed unknowns. Helper routine for computational homogenization.

Todo: Perhaps this routine should only give results for the dof it prescribes?

Parameters

C	Coefficient matrix to fill.

Definition at line 170 of file transportgradientdirichlet.C.

References oofem::Dof::__givePrescribedEquationNumber(), oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::GeneralBoundaryCondition::dofs, oofem::FEMComponent::domain, oofem::FEMComponent::giveDomain(), mCenterCoord, oofem::FloatMatrix::resize(), tractionControl, xis, and oofem::FloatMatrix::zero().

Referenced by computeField(), and computeTangent().

◆ computeField()

void oofem::TransportGradientDirichlet::computeField	(	FloatArray &	sigma,
		TimeStep *	tStep )

virtual

Computes the homogenized, macroscopic, field (stress).

Parameters

sigma	Output quantity (typically stress).
tStep	Active time step.

Definition at line 201 of file transportgradientdirichlet.C.

References oofem::EngngModel::assembleVector(), oofem::FloatArray::beTProductOf(), computeCoefficientMatrix(), oofem::FEMComponent::domain, domainSize(), oofem::FEMComponent::giveDomain(), oofem::FEMComponent::giveNumber(), oofem::EngngModel::giveNumberOfDomainEquations(), oofem::FloatArray::times(), and oofem::FloatArray::zero().

◆ computeTangent()

void oofem::TransportGradientDirichlet::computeTangent	(	FloatMatrix &	tangent,
		TimeStep *	tStep )

virtual

Computes the macroscopic tangent for homogenization problems through sensitivity analysis.

Parameters

tangent	Output tangent.
tStep	Active time step.

Todo: Get this from engineering model

Todo: This rotation matrix is not flexible enough.. it can only work with full size matrices and doesn't allow for flexibility in the matrixassembler.

Definition at line 223 of file transportgradientdirichlet.C.

References oofem::EngngModel::assemble(), oofem::FloatMatrix::at(), oofem::FloatMatrix::beTProductOf(), oofem::classFactory, computeCoefficientMatrix(), oofem::CR_CONVERGED, oofem::GeneralBoundaryCondition::dofs, oofem::FEMComponent::domain, domainSize(), oofem::Element_remote, oofem::Timer::getUtime(), oofem::FEMComponent::giveDomain(), oofem::FloatMatrix::giveNumberOfColumns(), oofem::FloatMatrix::giveNumberOfRows(), oofem::Element::giveParallelMode(), oofem::Element::giveRotationMatrix(), oofem::Element::isActivated(), oofem::FloatMatrix::isNotEmpty(), mCenterCoord, OOFEM_ERROR, OOFEM_LOG_INFO, oofem::FloatMatrix::resize(), oofem::FloatMatrix::rotatedWith(), oofem::ST_Petsc, oofem::Timer::startTimer(), oofem::Timer::stopTimer(), oofem::FloatMatrix::subtract(), and oofem::FloatMatrix::times().

◆ computeXi()

void oofem::TransportGradientDirichlet::computeXi ( )

Computes the offset values for "improved" Dirichlet. See class description.

Todo: Preallocation(?)

Todo: Can't do this yet, problem with material model interface (doesn't know the GP material mode). This should be changed.

Note: We can use the single constraint, but this assumes flat surfaces

Todo: Preallocation

Todo: Can't do this yet, problem with material model interface (doesn't know the GP material mode). This should be changed.

Definition at line 341 of file transportgradientdirichlet.C.

Referenced by postInitialize().

◆ domainSize()

double oofem::TransportGradientDirichlet::domainSize ( )

Definition at line 140 of file transportgradientdirichlet.C.

References oofem::IntArray::at(), oofem::FEMComponent::domain, oofem::FEInterpolation::evalNXIntegral(), oofem::FEMComponent::giveDomain(), oofem::Element::giveInterpolation(), oofem::IntArray::giveSize(), oofem::GeneralBoundaryCondition::set, surfSets, and tractionControl.

Referenced by computeField(), and computeTangent().

◆ give()

double oofem::TransportGradientDirichlet::give	(	Dof *	dof,
		ValueModeType	mode,
		double	time )

overridevirtual

Reimplemented from oofem::BoundaryCondition.

Definition at line 110 of file transportgradientdirichlet.C.

References oofem::FloatArray::add(), oofem::DofManager::giveCoordinates(), oofem::Dof::giveDofManager(), oofem::FEMComponent::giveNumber(), oofem::GeneralBoundaryCondition::giveTimeFunction(), mCenterCoord, mGradient, OOFEM_ERROR, and xis.

◆ giveClassName()

const char * oofem::TransportGradientDirichlet::giveClassName ( ) const

inlineoverridevirtual

Returns: Class name of the receiver.

Reimplemented from oofem::BoundaryCondition.

Definition at line 123 of file transportgradientdirichlet.h.

◆ giveInputRecord()

void oofem::TransportGradientDirichlet::giveInputRecord ( DynamicInputRecord & input )

overridevirtual

Setups the input record string of receiver.

Parameters

input Dynamic input record to be filled by receiver.

Reimplemented from oofem::BoundaryCondition.

Definition at line 88 of file transportgradientdirichlet.C.

References _IFT_TransportGradientDirichlet_centerCoords, _IFT_TransportGradientDirichlet_gradient, _IFT_TransportGradientDirichlet_surfSets, _IFT_TransportGradientDirichlet_tractionControl, mCenterCoord, mGradient, oofem::DynamicInputRecord::setField(), surfSets, and tractionControl.

◆ giveInputRecordName()

const char * oofem::TransportGradientDirichlet::giveInputRecordName ( ) const

inlineoverridevirtual

Returns: Input record name of the receiver.

Reimplemented from oofem::BoundaryCondition.

Definition at line 124 of file transportgradientdirichlet.h.

References _IFT_TransportGradientDirichlet_Name.

◆ giveType()

bcType oofem::TransportGradientDirichlet::giveType ( ) const

inlineoverridevirtual

Returns: Type of boundary condition. It allows to distinguish BC according its mathematical meaning, ie. like Dirichlet, Neumann, or Newton type.

Reimplemented from oofem::BoundaryCondition.

Definition at line 93 of file transportgradientdirichlet.h.

References oofem::DirichletBT.

◆ initializeFrom()

void oofem::TransportGradientDirichlet::initializeFrom ( InputRecord & ir )

overridevirtual

Initializes receiver according to object description stored in input record. This function is called immediately after creating object using constructor. Input record can be imagined as data record in component database belonging to receiver. Receiver may use value-name extracting functions to extract particular field from record. Note that initializeFrom may be called mutiple times.