Prescribes $T = g_{i}(x_i-\bar{x}_i)$ where $ T $ are primary unknown. More...

#include <transportgradientdirichlet.h>

Inheritance diagram for oofem::TransportGradientDirichlet:

Collaboration diagram for oofem::TransportGradientDirichlet:

Public Member Functions
	TransportGradientDirichlet (int n, Domain *d)
	Creates boundary condition with given number, belonging to given domain. More...

virtual	~TransportGradientDirichlet ()
	Destructor. More...

virtual IRResultType	initializeFrom (InputRecord *ir)
	Initializes receiver according to object description stored in input record. More...

virtual void	giveInputRecord (DynamicInputRecord &input)
	Setups the input record string of receiver. More...

virtual void	postInitialize ()
	Performs post initialization steps. More...

virtual double	give (Dof *dof, ValueModeType mode, double time)

virtual bcType	giveType () const

double	domainSize ()

void	computeCoefficientMatrix (FloatMatrix &C)
	Constructs a coefficient matrix for all prescribed unknowns. More...

virtual void	computeField (FloatArray &sigma, TimeStep *tStep)
	Computes the homogenized, macroscopic, field (stress). More...

virtual void	computeTangent (FloatMatrix &tangent, TimeStep *tStep)
	Computes the macroscopic tangent for homogenization problems through sensitivity analysis. More...

void	computeXi ()
	Computes the offset values for "improved" Dirichlet. See class description. More...

virtual void	scale (double s)
	Scales the receiver according to given value. More...

virtual const char *	giveClassName () const

virtual const char *	giveInputRecordName () const

Public Member Functions inherited from oofem::BoundaryCondition
	BoundaryCondition (int i, Domain *d)
	Constructor. More...

virtual	~BoundaryCondition ()
	Destructor. More...

virtual double	give (Dof dof, ValueModeType mode, TimeStep tStep)
	Returns the value of a prescribed unknown, respecting requested mode for given time. More...

void	setPrescribedValue (double s)
	Set prescribed value at the input record string of receiver. More...

virtual contextIOResultType	saveContext (DataStream &stream, ContextMode mode, void *obj=NULL)
	Stores receiver state to output stream. More...

virtual contextIOResultType	restoreContext (DataStream &stream, ContextMode mode, void *obj=NULL)
	Restores the receiver state previously written in stream. More...

Public Member Functions inherited from oofem::GeneralBoundaryCondition
	GeneralBoundaryCondition (int n, Domain *d)
	Constructor. More...

virtual	~GeneralBoundaryCondition ()
	Destructor. More...

int	giveSetNumber ()
	Gives the set number which boundary condition is applied to. More...

virtual int	giveNumberOfInternalDofManagers ()
	Gives the number of internal dof managers. More...

virtual DofManager *	giveInternalDofManager (int i)
	Gives an internal dof manager from receiver. More...

Function *	giveTimeFunction ()

virtual bcValType	giveBCValType () const
	Returns receiver load type. More...

virtual bool	isImposed (TimeStep *tStep)
	Returns nonzero if receiver representing BC is imposed at given time, otherwise returns zero. More...

virtual const IntArray &	giveDofIDs () const
	Array with default dofs which b.c. More...

virtual bcGeomType	giveBCGeoType () const
	Returns geometry character of boundary condition. More...

Public Member Functions inherited from oofem::FEMComponent
	FEMComponent (int n, Domain *d)
	Regular constructor, creates component with given number and belonging to given domain. More...

virtual	~FEMComponent ()
	Virtual destructor. More...

Domain *	giveDomain () const

virtual void	setDomain (Domain *d)
	Sets associated Domain. More...

int	giveNumber () const

void	setNumber (int num)
	Sets number of receiver. More...

virtual void	updateLocalNumbering (EntityRenumberingFunctor &f)
	Local renumbering support. More...

virtual int	checkConsistency ()
	Allows programmer to test some internal data, before computation begins. More...

virtual void	printOutputAt (FILE file, TimeStep tStep)
	Prints output of receiver to stream, for given time step. More...

virtual void	printYourself ()
	Prints receiver state on stdout. Useful for debugging. More...

virtual Interface *	giveInterface (InterfaceType t)
	Interface requesting service. More...

std::string	errorInfo (const char *func) const
	Returns string for prepending output (used by error reporting macros). More...

Protected Attributes
FloatArray	mGradient

FloatArray	mCenterCoord

bool	tractionControl

std::map< int, FloatArray >	xis
	Stores one "psi" value for each node. More...

IntArray	surfSets

Protected Attributes inherited from oofem::BoundaryCondition
FloatArray	values
	Prescribed values for each resp. dof. More...

Protected Attributes inherited from oofem::GeneralBoundaryCondition
int	timeFunction
	Associated load time function. More...

bcValType	valType
	Physical meaning of BC value. More...

IntArray	dofs
	Dofs that b.c. is applied to (relevant for Dirichlet type b.c.s). More...

int	isImposedTimeFunction
	Zero by default - the BC is than always imposed. More...

int	set
	Set number for boundary condition to be applied to. More...

Protected Attributes inherited from oofem::FEMComponent
int	number
	Component number. More...

Domain *	domain
	Link to domain object, useful for communicating with other FEM components. More...

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

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

inline

Creates boundary condition with given number, belonging to given domain.

Parameters

n	Boundary condition number.
d	Domain to which new object will belongs.

Definition at line 90 of file transportgradientdirichlet.h.

virtual oofem::TransportGradientDirichlet::~TransportGradientDirichlet ( )

inlinevirtual

Destructor.

Definition at line 93 of file transportgradientdirichlet.h.

Member Function Documentation

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 172 of file transportgradientdirichlet.C.

References oofem::Dof::__givePrescribedEquationNumber(), oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::GeneralBoundaryCondition::dofs, oofem::FEMComponent::domain, oofem::Domain::giveDofManagers(), oofem::FEMComponent::giveDomain(), oofem::Domain::giveEngngModel(), oofem::Domain::giveNumber(), oofem::EngngModel::giveNumberOfDomainEquations(), oofem::Domain::giveNumberOfSpatialDimensions(), mCenterCoord, oofem::FloatMatrix::resize(), tractionControl, xis, and oofem::FloatMatrix::zero().

Referenced by computeField(), and computeTangent().

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 203 of file transportgradientdirichlet.C.

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

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 225 of file transportgradientdirichlet.C.

References oofem::EngngModel::assemble(), oofem::FloatArray::at(), oofem::FloatMatrix::beTProductOf(), oofem::classFactory, computeCoefficientMatrix(), oofem::ClassFactory::createSparseLinSolver(), oofem::ClassFactory::createSparseMtrx(), oofem::GeneralBoundaryCondition::dofs, oofem::FEMComponent::domain, domainSize(), oofem::Element_remote, oofem::Timer::getUtime(), oofem::Domain::giveDofManagers(), oofem::FEMComponent::giveDomain(), oofem::Domain::giveElement(), oofem::Domain::giveEngngModel(), oofem::Domain::giveNumberOfElements(), oofem::Domain::giveNumberOfSpatialDimensions(), oofem::Element::giveParallelMode(), oofem::Element::giveRotationMatrix(), oofem::Element::isActivated(), oofem::FloatMatrix::isNotEmpty(), mCenterCoord, NM_NoSuccess, OOFEM_ERROR, OOFEM_LOG_INFO, oofem::FloatMatrix::rotatedWith(), oofem::ST_Petsc, oofem::Timer::startTimer(), oofem::Timer::stopTimer(), and oofem::FloatMatrix::times().

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 343 of file transportgradientdirichlet.C.

Referenced by postInitialize().

double oofem::TransportGradientDirichlet::domainSize ( )

Definition at line 142 of file transportgradientdirichlet.C.

References oofem::IntArray::at(), oofem::FEMComponent::domain, oofem::FEInterpolation::evalNXIntegral(), oofem::Set::giveBoundaryList(), oofem::FEMComponent::giveDomain(), oofem::Domain::giveElement(), oofem::Element::giveInterpolation(), oofem::Domain::giveNumberOfSpatialDimensions(), oofem::Domain::giveSet(), oofem::IntArray::giveSize(), surfSets, and tractionControl.

Referenced by computeField(), and computeTangent().

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

virtual

Reimplemented from oofem::BoundaryCondition.

Definition at line 111 of file transportgradientdirichlet.C.

References oofem::FloatArray::add(), oofem::FloatArray::beDifferenceOf(), oofem::FloatArray::dotProduct(), oofem::Function::evaluateAccelerationAtTime(), oofem::Function::evaluateAtTime(), oofem::Function::evaluateVelocityAtTime(), oofem::DofManager::giveCoordinates(), oofem::Dof::giveDofManager(), oofem::FEMComponent::giveNumber(), oofem::FloatArray::giveSize(), oofem::GeneralBoundaryCondition::giveTimeFunction(), mCenterCoord, mGradient, OOFEM_ERROR, and xis.

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

inlinevirtual

Returns: Class name of the receiver.

Reimplemented from oofem::BoundaryCondition.

Definition at line 131 of file transportgradientdirichlet.h.

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

virtual

Setups the input record string of receiver.

Parameters

input Dynamic input record to be filled by receiver.

Reimplemented from oofem::BoundaryCondition.

Definition at line 89 of file transportgradientdirichlet.C.

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

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

inlinevirtual

Returns: Input record name of the receiver.

Reimplemented from oofem::BoundaryCondition.

Definition at line 132 of file transportgradientdirichlet.h.

References _IFT_TransportGradientDirichlet_Name.

virtual bcType oofem::TransportGradientDirichlet::giveType ( ) const

inlinevirtual

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 101 of file transportgradientdirichlet.h.

References oofem::DirichletBT.

IRResultType oofem::TransportGradientDirichlet::initializeFrom ( InputRecord * ir )

virtual

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.