Imposes weak periodicity on the doftype of choice. More...

#include <weakperiodicbc.h>

Inheritance diagram for oofem::WeakPeriodicBoundaryCondition:

Collaboration diagram for oofem::WeakPeriodicBoundaryCondition:

Public Member Functions
	WeakPeriodicBoundaryCondition (int n, Domain *d)

virtual	~WeakPeriodicBoundaryCondition ()

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

basisType	giveBasisType ()

virtual void	assemble (SparseMtrx &answer, TimeStep *tStep, CharType type, const UnknownNumberingScheme &r_s, const UnknownNumberingScheme &c_s, double scale=1.0)
	Assembles B.C. More...

virtual void	assembleVector (FloatArray &answer, TimeStep tStep, CharType type, ValueModeType mode, const UnknownNumberingScheme &s, FloatArray eNorm=NULL)
	Assembles B.C. More...

void	giveInternalForcesVector (FloatArray &answer, TimeStep tStep, CharType type, ValueModeType mode, const UnknownNumberingScheme &s, FloatArray eNorm=NULL)

void	giveExternalForcesVector (FloatArray &answer, TimeStep *tStep, CharType type, ValueModeType mode, const UnknownNumberingScheme &s)

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

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

virtual void	addElementSide (int elem, int side)
	Adds element for active boundary condition. More...

virtual const char *	giveClassName () const

virtual const char *	giveInputRecordName () const

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

virtual	~ActiveBoundaryCondition ()
	Destructor. More...

virtual void	giveLocationArrays (std::vector< IntArray > &rows, std::vector< IntArray > &cols, CharType type, const UnknownNumberingScheme &r_s, const UnknownNumberingScheme &c_s)
	Gives a list of location arrays that will be assembled. More...

virtual bool	requiresActiveDofs ()
	Checks to see if active boundary condition requires special DOFs. More...

virtual bool	isPrimaryDof (ActiveDof *dof)
	Checks to see if the dof is a primary DOF. More...

virtual double	giveBcValue (Dof dof, ValueModeType mode, TimeStep tStep)
	Returns the prescribed value of a dof (if any). More...

virtual bool	hasBc (Dof dof, TimeStep tStep)
	Returns the prescribed value of a dof (if any). More...

virtual int	giveNumberOfMasterDofs (ActiveDof *dof)
	Allows for active boundary conditions to handle their own special DOF. More...

virtual Dof *	giveMasterDof (ActiveDof *dof, int mdof)
	Give the pointer to master dof belonging to active DOF. More...

virtual void	computeDofTransformation (ActiveDof *dof, FloatArray &masterContribs)

virtual double	giveUnknown (PrimaryField &field, ValueModeType mode, TimeStep tStep, ActiveDof dof)
	Computes the value of the dof. More...

virtual double	giveUnknown (ValueModeType mode, TimeStep tStep, ActiveDof dof)
	Computes the value of the dof. 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...

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 bcType	giveType () const

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

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

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

virtual void	giveInputRecord (DynamicInputRecord &input)
	Setups 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::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 Member Functions
void	computeElementTangent (FloatMatrix &answer, Element e, int boundary, TimeStep tStep)

Private Member Functions
void	giveEdgeNormal (FloatArray &answer, int element, int side)

void	updateSminmax ()

void	updateDirection ()

double	computeBaseFunctionValue (int baseID, FloatArray coordinate)

double	computeBaseFunctionValue1D (int baseID, double coordinate)

double	computeBaseFunctionValue2D (int baseID, FloatArray coordinate)

double	factorial (int n)

double	binomial (double n, int k)

void	getExponents (int n, int &i, int &j)
	Compute exponent for term n. More...

void	computeOrthogonalBasis ()
	Compute orthogonal polynomial basis using Gram-Smidth process. More...

double	computeProjectionCoefficient (int vIndex, int uIndex)

void	computeDeformationGradient (FloatMatrix &answer, Element e, FloatArray lcoord, TimeStep *tStep)

Private Attributes
basisType	useBasisType

int	bcID

int	orderOfPolygon

FloatArray	g
	Contains prescribed gradient. More...

int	direction
	Direction of normal. More...

IntArray	surfaceIndexes
	Keeps info on which coordinates varies over the surface. More...

FloatArray	smax

FloatArray	smin

bool	doUpdateSminmax

int	ngp
	Number of Gausspoints used when integrating along the element edges. More...

int	ndof
	Number of degrees of freedom (number of terms) More...

int	posSet
	Set containing positive side. More...

int	negSet
	Set containing negative side. More...

IntArray	dofids
	ID of dofs on which weak periodicity is imposed. More...

signed int	sideSign [2]
	sideSign is the sign of the normal for each side More...

std::vector< int >	side [2]
	side[] keeps track of which side of the triangle is located along the boundary. More...

std::vector< int >	element [2]

std::vector< double >	directions
	Keeps track of which coordinate(s) are changing on the surface/edge. More...

std::unique_ptr< Node >	gammaDman

IntArray	gamma_ids

FloatMatrix	gsMatrix
	gsMatrix contains coefficients for the Gram-Schmidt polynomials More...

int	tcount
	Number of terms in polynomial. More...

int	ndofids
	Number of dofIDs. More...

bool	nlgeo
	Use finite strains? More...

Additional Inherited Members
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

Imposes weak periodicity on the doftype of choice.

2D. It is required that the two edges are parallel and either horizontal or vertical.

Author: Carl Sandström

Definition at line 75 of file weakperiodicbc.h.

Constructor & Destructor Documentation

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

Definition at line 63 of file weakperiodicbc.C.

References doUpdateSminmax, oofem::monomial, and useBasisType.

oofem::WeakPeriodicBoundaryCondition::~WeakPeriodicBoundaryCondition ( )

virtual

Definition at line 69 of file weakperiodicbc.C.

Member Function Documentation

void oofem::WeakPeriodicBoundaryCondition::addElementSide	(	int	elem,
		int	side
	)

virtual

Adds element for active boundary condition.

Parameters

elem	Element number.
side	Side number.

Reimplemented from oofem::ActiveBoundaryCondition.

Definition at line 361 of file weakperiodicbc.C.

References oofem::FloatArray::at(), direction, element, giveEdgeNormal(), and side.

void oofem::WeakPeriodicBoundaryCondition::assemble	(	SparseMtrx &	answer,
		TimeStep *	tStep,
		CharType	type,
		const UnknownNumberingScheme &	r_s,
		const UnknownNumberingScheme &	c_s,
		double	scale = `1.0`
	)

virtual

Assembles B.C.

contributions to specified matrix.

Parameters

[in,out]	answer	Matrix to assemble to.
	tStep	Active time step.
	tStep	Active time step.
	type	Type of matrix to assemble.
	r_s	Row numbering scheme.
	c_s	Column numbering scheme.
	scale	Scaling factor.

Reimplemented from oofem::ActiveBoundaryCondition.

Definition at line 491 of file weakperiodicbc.C.

void oofem::WeakPeriodicBoundaryCondition::assembleVector	(	FloatArray &	answer,
		TimeStep *	tStep,
		CharType	type,
		ValueModeType	mode,
		const UnknownNumberingScheme &	s,
		FloatArray *	eNorms = `NULL`
	)

virtual

Assembles B.C.

contributions to specified vector.

Parameters

[in,out]	answer	Vector to assemble to.
	tStep	Active time step.
	type	Type of matrix to assemble.
	mode	Mode of value.
	s	Numbering scheme.
	eNorms	Norms for each dofid.

Returns: Equivalent of the sum of the element norm (squared) of assembled vector.

Reimplemented from oofem::ActiveBoundaryCondition.

Definition at line 656 of file weakperiodicbc.C.

References giveExternalForcesVector(), and giveInternalForcesVector().

double oofem::WeakPeriodicBoundaryCondition::binomial	(	double	n,
		int	k
	)

private

Definition at line 881 of file weakperiodicbc.C.

Referenced by computeBaseFunctionValue1D().

double oofem::WeakPeriodicBoundaryCondition::computeBaseFunctionValue	(	int	baseID,
		FloatArray	coordinate
	)

private

Definition at line 598 of file weakperiodicbc.C.

References oofem::IntArray::at(), oofem::FloatArray::at(), computeBaseFunctionValue1D(), computeBaseFunctionValue2D(), oofem::FEMComponent::domain, oofem::Domain::giveNumberOfSpatialDimensions(), and surfaceIndexes.

Referenced by assemble(), computeElementTangent(), giveExternalForcesVector(), and giveInternalForcesVector().

double oofem::WeakPeriodicBoundaryCondition::computeBaseFunctionValue1D	(	int	baseID,
		double	coordinate
	)

private

Definition at line 626 of file weakperiodicbc.C.

References oofem::FloatArray::at(), binomial(), oofem::FloatArray::giveSize(), oofem::legendre, M_PI, oofem::monomial, oofem::FloatArray::resize(), smax, smin, oofem::trigonometric, and useBasisType.

Referenced by computeBaseFunctionValue().

double oofem::WeakPeriodicBoundaryCondition::computeBaseFunctionValue2D	(	int	baseID,
		FloatArray	coordinate
	)

private

Definition at line 607 of file weakperiodicbc.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), getExponents(), gsMatrix, oofem::legendre, oofem::monomial, ndof, and useBasisType.

Referenced by computeBaseFunctionValue().

void oofem::WeakPeriodicBoundaryCondition::computeDeformationGradient	(	FloatMatrix &	answer,
		Element *	e,
		FloatArray *	lcoord,
		TimeStep *	tStep
	)

private

Definition at line 392 of file weakperiodicbc.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::FloatMatrix::beInverseOf(), oofem::FloatMatrix::beMatrixForm(), oofem::FloatArray::beProductOf(), oofem::FloatMatrix::beTranspositionOf(), oofem::Element::computeVectorOf(), dofids, oofem::FEInterpolation::evaldNdx(), oofem::Element::giveInterpolation(), oofem::FloatMatrix::giveNumberOfRows(), oofem::FloatMatrix::resize(), and oofem::FloatMatrix::zero().

Referenced by assemble(), computeElementTangent(), and giveInternalForcesVector().

void oofem::WeakPeriodicBoundaryCondition::computeElementTangent	(	FloatMatrix &	answer,
		Element *	e,
		int	boundary,
		TimeStep *	tStep
	)

protected

void oofem::WeakPeriodicBoundaryCondition::computeOrthogonalBasis ( )

private

Compute orthogonal polynomial basis using Gram-Smidth process.

Definition at line 158 of file weakperiodicbc.C.

References oofem::FloatArray::at(), oofem::FloatMatrix::at(), computeProjectionCoefficient(), gsMatrix, ndof, oofem::FloatArray::resize(), oofem::FloatMatrix::resize(), oofem::FloatArray::times(), and oofem::FloatMatrix::zero().

Referenced by updateSminmax().

double oofem::WeakPeriodicBoundaryCondition::computeProjectionCoefficient	(	int	vIndex,
		int	uIndex
	)

private

Definition at line 188 of file weakperiodicbc.C.

References oofem::IntArray::at(), oofem::FloatArray::at(), oofem::FloatMatrix::at(), oofem::FEInterpolation::boundaryGiveTransformationJacobian(), oofem::FEInterpolation::boundaryLocal2Global(), oofem::FEMComponent::domain, element, getExponents(), oofem::FEInterpolation::giveBoundaryIntegrationRule(), oofem::Domain::giveElement(), oofem::Element::giveInterpolation(), gsMatrix, ndof, side, and surfaceIndexes.

Referenced by computeOrthogonalBasis().

double oofem::WeakPeriodicBoundaryCondition::factorial ( int n )

private

Definition at line 872 of file weakperiodicbc.C.

void oofem::WeakPeriodicBoundaryCondition::getExponents	(	int	n,
		int &	i,
		int &	j
	)

private

Compute exponent for term n.

Exponents i and j are x^i*y^j

Definition at line 890 of file weakperiodicbc.C.

Referenced by computeBaseFunctionValue2D(), and computeProjectionCoefficient().

basisType oofem::WeakPeriodicBoundaryCondition::giveBasisType ( )

inline

Definition at line 166 of file weakperiodicbc.h.

Referenced by oofem::MatlabExportModule::doOutputSpecials().

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

inlinevirtual

Returns: Class name of the receiver.

Implements oofem::FEMComponent.

Definition at line 189 of file weakperiodicbc.h.

void oofem::WeakPeriodicBoundaryCondition::giveEdgeNormal	(	FloatArray &	answer,
		int	element,
		int	side
	)

private

Definition at line 239 of file weakperiodicbc.C.

References oofem::FEInterpolation::boundaryEvalNormal(), dofids, oofem::FEMComponent::domain, oofem::Domain::giveElement(), oofem::Element::giveInterpolation(), oofem::Domain::giveNumberOfSpatialDimensions(), and oofem::FloatArray::resize().

Referenced by addElementSide(), updateDirection(), and updateSminmax().

void oofem::WeakPeriodicBoundaryCondition::giveExternalForcesVector	(	FloatArray &	answer,
		TimeStep *	tStep,
		CharType	type,
		ValueModeType	mode,
		const UnknownNumberingScheme &	s
	)

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

inlinevirtual

Returns: Input record name of the receiver.

Implements oofem::FEMComponent.

Definition at line 190 of file weakperiodicbc.h.

References _IFT_WeakPeriodicBoundaryCondition_Name.

DofManager * oofem::WeakPeriodicBoundaryCondition::giveInternalDofManager ( int i )

virtual

Gives an internal dof manager from receiver.

Reimplemented from oofem::GeneralBoundaryCondition.

Definition at line 863 of file weakperiodicbc.C.

References gammaDman.

Referenced by oofem::MatlabExportModule::doOutputSpecials().

void oofem::WeakPeriodicBoundaryCondition::giveInternalForcesVector	(	FloatArray &	answer,
		TimeStep *	tStep,
		CharType	type,
		ValueModeType	mode,
		const UnknownNumberingScheme &	s,
		FloatArray *	eNorm = `NULL`
	)

Definition at line 669 of file weakperiodicbc.C.

References oofem::FloatArray::assemble(), oofem::FloatArray::assembleSquared(), oofem::FloatMatrix::at(), oofem::FloatMatrix::beInverseOf(), oofem::FloatMatrix::beProductOf(), oofem::FloatMatrix::beTProductOf(), oofem::FloatMatrix::beTranspositionOf(), oofem::FEInterpolation::boundaryEvalN(), oofem::FEInterpolation::boundaryGiveNodes(), oofem::FEInterpolation::boundaryGiveTransformationJacobian(), oofem::FEInterpolation::boundaryLocal2Global(), computeBaseFunctionValue(), oofem::Element::computeBoundaryVectorOf(), computeDeformationGradient(), dofids, oofem::FEMComponent::domain, element, gamma_ids, gammaDman, oofem::FEInterpolation::giveBoundaryIntegrationRule(), oofem::Element::giveBoundaryLocationArray(), oofem::Domain::giveElement(), oofem::Element::giveInterpolation(), oofem::IntArray::giveSize(), oofem::FloatArray::giveSize(), oofem::FEInterpolation::global2local(), N, ndof, ndofids, nlgeo, orderOfPolygon, oofem::FloatArray::plusProduct(), oofem::FloatArray::resize(), oofem::FloatMatrix::resize(), side, sideSign, tcount, updateSminmax(), oofem::FloatArray::zero(), and oofem::FloatMatrix::zero().

Referenced by assembleVector().

int oofem::WeakPeriodicBoundaryCondition::giveNumberOfInternalDofManagers ( )

virtual

Gives the number of internal dof managers.

Reimplemented from oofem::GeneralBoundaryCondition.

Definition at line 858 of file weakperiodicbc.C.

Referenced by oofem::MatlabExportModule::doOutputSpecials().

IRResultType oofem::WeakPeriodicBoundaryCondition::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.