OOFEM  2.4
OOFEM.org - Object Oriented Finite Element Solver
tet1bubblestokes.h
Go to the documentation of this file.
1 /*
2  *
3  * ##### ##### ###### ###### ### ###
4  * ## ## ## ## ## ## ## ### ##
5  * ## ## ## ## #### #### ## # ##
6  * ## ## ## ## ## ## ## ##
7  * ## ## ## ## ## ## ## ##
8  * ##### ##### ## ###### ## ##
9  *
10  *
11  * OOFEM : Object Oriented Finite Element Code
12  *
13  * Copyright (C) 1993 - 2013 Borek Patzak
14  *
15  *
16  *
17  * Czech Technical University, Faculty of Civil Engineering,
18  * Department of Structural Mechanics, 166 29 Prague, Czech Republic
19  *
20  * This library is free software; you can redistribute it and/or
21  * modify it under the terms of the GNU Lesser General Public
22  * License as published by the Free Software Foundation; either
23  * version 2.1 of the License, or (at your option) any later version.
24  *
25  * This program is distributed in the hope that it will be useful,
26  * but WITHOUT ANY WARRANTY; without even the implied warranty of
27  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
28  * Lesser General Public License for more details.
29  *
30  * You should have received a copy of the GNU Lesser General Public
31  * License along with this library; if not, write to the Free Software
32  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
33  */
34 
35 #ifndef tet1bubblestokes_h
36 #define tet1bubblestokes_h
37 
38 #include "fmelement.h"
39 #include "zznodalrecoverymodel.h"
40 #include "spatiallocalizer.h"
41 #include "elementinternaldofman.h"
42 #include "matresponsemode.h"
43 
44 #include <memory>
45 
46 #define _IFT_Tet1BubbleStokes_Name "tet1bubblestokes"
47 
48 namespace oofem {
49 class FEI3dTetLin;
50 
59 class Tet1BubbleStokes : public FMElement,
62 {
63 protected:
69  static IntArray edge_ordering [ 6 ];
71  static IntArray surf_ordering [ 4 ];
72 
74  std :: unique_ptr< ElementDofManager > bubble;
75  // Coordinates associated with the bubble dofs.
76  //FloatArray bubbleCoord; // Assumed fixed at 0 for now (i.e. only linear geometry)
77 
78 public:
79  Tet1BubbleStokes(int n, Domain * d);
80  virtual ~Tet1BubbleStokes();
81 
82  virtual double computeVolumeAround(GaussPoint *gp);
83 
84  virtual void computeGaussPoints();
85 
86  virtual void giveCharacteristicVector(FloatArray &answer, CharType type, ValueModeType mode, TimeStep *tStep);
87  virtual void giveCharacteristicMatrix(FloatMatrix &answer, CharType type, TimeStep *tStep);
88 
89  void computeInternalForcesVector(FloatArray &answer, TimeStep *tStep);
90  void computeStiffnessMatrix(FloatMatrix &answer, MatResponseMode mode, TimeStep *tStep);
91 
92  void computeExternalForcesVector(FloatArray &answer, TimeStep *tStep);
93  virtual void computeLoadVector(FloatArray &answer, BodyLoad *load, CharType type, ValueModeType mode, TimeStep *tStep);
94  virtual void computeBoundarySurfaceLoadVector(FloatArray &answer, BoundaryLoad *load, int boundary, CharType type, ValueModeType mode, TimeStep *tStep, bool global=true);
95 
96  virtual const char *giveClassName() const { return "Tet1BubbleStokes"; }
97  virtual const char *giveInputRecordName() const { return _IFT_Tet1BubbleStokes_Name; }
98  virtual MaterialMode giveMaterialMode() { return _3dFlow; }
99 
100  virtual int computeNumberOfDofs();
101 
102  virtual int giveNumberOfInternalDofManagers() const { return 1; }
103  virtual DofManager *giveInternalDofManager(int i) const { return bubble.get(); }
104  virtual void giveInternalDofManDofIDMask(int i, IntArray &answer) const;
105 
106  virtual FEInterpolation *giveInterpolation() const;
107  virtual FEInterpolation *giveInterpolation(DofIDItem id) const;
108 
109  virtual void giveDofManDofIDMask(int inode, IntArray &answer) const;
110 
111  virtual void updateYourself(TimeStep *tStep);
112 
114 
115  virtual void computeField(ValueModeType u, TimeStep *tStep, const FloatArray &coords, FloatArray &answer);
116 };
117 } // end namespace oofem
118 #endif // tet1bubblestokes_h
virtual Interface * giveInterface(InterfaceType it)
Interface requesting service.
#define _IFT_Tet1BubbleStokes_Name
Class and object Domain.
Definition: domain.h:115
static IntArray momentum_ordering
Ordering of dofs in element. Used to assemble the element stiffness.
void computeInternalForcesVector(FloatArray &answer, TimeStep *tStep)
virtual void computeField(ValueModeType u, TimeStep *tStep, const FloatArray &coords, FloatArray &answer)
Computes the unknown vector interpolated at the specified local coordinates.
The element interface required by ZZNodalRecoveryModel.
void computeExternalForcesVector(FloatArray &answer, TimeStep *tStep)
ValueModeType
Type representing the mode of UnknownType or CharType, or similar types.
Definition: valuemodetype.h:78
Class implementing element body load, acting over whole element volume (e.g., the dead weight)...
Definition: bodyload.h:49
Class representing implementation of linear tetrahedra interpolation class.
Definition: fei3dtetlin.h:44
virtual int giveNumberOfInternalDofManagers() const
Base class for dof managers.
Definition: dofmanager.h:113
virtual void giveDofManDofIDMask(int inode, IntArray &answer) const
Returns dofmanager dof mask for node.
static IntArray surf_ordering[4]
Ordering of dofs on surfaces. Used to assemble surface loads.
MaterialMode
Type representing material mode of integration point.
Definition: materialmode.h:89
Class implementing an array of integers.
Definition: intarray.h:61
MatResponseMode
Describes the character of characteristic material matrix.
virtual void giveInternalDofManDofIDMask(int i, IntArray &answer) const
Returns internal dofmanager dof mask for node.
Class representing a general abstraction for finite element interpolation class.
Definition: feinterpol.h:132
virtual void updateYourself(TimeStep *tStep)
Updates element state after equilibrium in time step has been reached.
Tet1BubbleStokes(int n, Domain *d)
Abstract base class representing a boundary load (force, momentum, ...) that acts directly on a bound...
Definition: boundaryload.h:110
virtual int computeNumberOfDofs()
Computes or simply returns total number of element&#39;s local DOFs.
virtual MaterialMode giveMaterialMode()
Returns material mode for receiver integration points.
static FEI3dTetLin interp
Interpolation for pressure.
DofIDItem
Type representing particular dof type.
Definition: dofiditem.h:86
virtual void giveCharacteristicVector(FloatArray &answer, CharType type, ValueModeType mode, TimeStep *tStep)
Computes characteristic vector of receiver of requested type in given time step.
virtual double computeVolumeAround(GaussPoint *gp)
Returns volume related to given integration point.
virtual void giveCharacteristicMatrix(FloatMatrix &answer, CharType type, TimeStep *tStep)
Computes characteristic matrix of receiver of requested type in given time step.
static IntArray conservation_ordering
virtual FEInterpolation * giveInterpolation() const
Class representing vector of real numbers.
Definition: floatarray.h:82
void computeStiffnessMatrix(FloatMatrix &answer, MatResponseMode mode, TimeStep *tStep)
This abstract class represent a general base element class for fluid dynamic problems.
Definition: fmelement.h:54
virtual void computeLoadVector(FloatArray &answer, BodyLoad *load, CharType type, ValueModeType mode, TimeStep *tStep)
Computes the contribution of the given body load (volumetric).
Implementation of matrix containing floating point numbers.
Definition: floatmatrix.h:94
std::unique_ptr< ElementDofManager > bubble
The extra dofs from the bubble.
virtual void computeBoundarySurfaceLoadVector(FloatArray &answer, BoundaryLoad *load, int boundary, CharType type, ValueModeType mode, TimeStep *tStep, bool global=true)
Computes the contribution of the given load at the given boundary surface in global coordinate system...
CharType
Definition: chartype.h:87
Class Interface.
Definition: interface.h:82
The spatial localizer element interface associated to spatial localizer.
InterfaceType
Enumerative type, used to identify interface type.
Definition: interfacetype.h:43
virtual DofManager * giveInternalDofManager(int i) const
Returns i-th internal element dof manager of the receiver.
the oofem namespace is to define a context or scope in which all oofem names are defined.
virtual void computeGaussPoints()
Initializes the array of integration rules member variable.
Class representing integration point in finite element program.
Definition: gausspoint.h:93
Class representing solution step.
Definition: timestep.h:80
virtual const char * giveClassName() const
virtual const char * giveInputRecordName() const
Tetrahedral element for Stokes flow using Bubble basis function for stabilization.
static IntArray edge_ordering[6]
Ordering of dofs on edges. Used to assemble edge loads.

This page is part of the OOFEM documentation. Copyright (c) 2011 Borek Patzak
Project e-mail: info@oofem.org
Generated at Tue Jan 2 2018 20:07:31 for OOFEM by doxygen 1.8.11 written by Dimitri van Heesch, © 1997-2011