OOFEM  2.4
OOFEM.org - Object Oriented Finite Element Solver
tr_shell02.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 tr_shell02_h
36 #define tr_shell02_h
37 
38 #include "../sm/Elements/structuralelement.h"
39 #include "zznodalrecoverymodel.h"
40 #include "../sm/ErrorEstimators/zzerrorestimator.h"
41 #include "../sm/Elements/Plates/dkt3d.h"
42 #include "../sm/Elements/PlaneStress/trplanestressrotallman3d.h"
43 #include "spatiallocalizer.h"
44 
45 #include <memory>
46 
47 #define _IFT_TR_SHELL02_Name "tr_shell02"
48 
49 namespace oofem {
57 {
58 protected:
60  std :: unique_ptr< DKTPlate3d > plate;
62  std :: unique_ptr< TrPlanestressRotAllman3d > membrane;
68  std :: unique_ptr< IntegrationRule > compositeIR;
69 
72 
73 public:
75  TR_SHELL02(int n, Domain * d);
77  virtual ~TR_SHELL02() {}
78 
79  virtual FEInterpolation *giveInterpolation() const { return plate->giveInterpolation(); }
80 
81  virtual int computeNumberOfDofs() { return 18; }
82  virtual void giveDofManDofIDMask(int inode, IntArray &answer) const
83  { plate->giveDofManDofIDMask(inode, answer); }
84  // definition & identification
85  virtual const char *giveInputRecordName() const { return _IFT_TR_SHELL02_Name; }
86  virtual const char *giveClassName() const { return "TR_SHELL02"; }
88 
89  virtual void giveCharacteristicVector(FloatArray &answer, CharType mtrx, ValueModeType mode, TimeStep *tStep);
90  virtual void giveCharacteristicMatrix(FloatMatrix &answer, CharType mtrx, TimeStep *tStep);
91  virtual double computeVolumeAround(GaussPoint *gp);
92  virtual bool giveRotationMatrix(FloatMatrix &answer);
93 
94  virtual void updateYourself(TimeStep *tStep);
95  virtual void updateInternalState(TimeStep *tStep);
96  virtual void printOutputAt(FILE *file, TimeStep *tStep);
97  virtual contextIOResultType saveContext(DataStream &stream, ContextMode mode, void *obj = NULL);
98  virtual contextIOResultType restoreContext(DataStream &stream, ContextMode mode, void *obj = NULL);
99  virtual void postInitialize();
101  void setCrossSection(int csIndx);
102 #ifdef __OOFEG
103  virtual void drawRawGeometry(oofegGraphicContext &gc, TimeStep *tStep);
104  virtual void drawDeformedGeometry(oofegGraphicContext &gc, TimeStep *tStep, UnknownType type);
105  virtual void drawScalar(oofegGraphicContext &gc, TimeStep *tStep);
106 #endif
107  // the membrane and plate irules are same (chacked in initializeFrom)
108  virtual int giveDefaultIntegrationRule() const { return plate->giveDefaultIntegrationRule(); }
109  virtual IntegrationRule *giveDefaultIntegrationRulePtr() { return plate->giveDefaultIntegrationRulePtr(); }
110  virtual IntegrationRule *giveIntegrationRule(int i) { return plate->giveIntegrationRule(i); }
111  virtual Element_Geometry_Type giveGeometryType() const { return EGT_triangle_1; }
113  virtual MaterialMode giveMaterialMode() { return _Unknown; }
114 
116  virtual int giveIPValue(FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep);
117 
119  InternalStateType type, TimeStep *tStep);
120 
122  virtual void ZZErrorEstimatorI_computeLocalStress(FloatArray &answer, FloatArray &sig);
123 
124  // SpatialLocalizerI
125  virtual void SpatialLocalizerI_giveBBox(FloatArray &bb0, FloatArray &bb1);
126 
127 
128  virtual int computeGlobalCoordinates(FloatArray &answer, const FloatArray &lcoords) {
129  return this->plate->computeGlobalCoordinates(answer, lcoords);
130  }
131 
132  virtual bool computeLocalCoordinates(FloatArray &answer, const FloatArray &gcoords) {
133  return this->plate->computeLocalCoordinates(answer, gcoords);
134  }
135 
136 protected:
137  virtual void computeBmatrixAt(GaussPoint *, FloatMatrix &, int = 1, int = ALL_STRAINS)
138  { OOFEM_ERROR("calling of this function is not allowed"); }
139  virtual void computeNmatrixAt(const FloatArray &iLocCoord, FloatMatrix &)
140  { OOFEM_ERROR("calling of this function is not allowed"); }
141 
143 protected:
144  virtual void computeGaussPoints()
145  {
146  this->membrane->computeGaussPoints();
147  this->plate->computeGaussPoints();
148  }
149  virtual void computeStressVector(FloatArray &answer, const FloatArray &strain, GaussPoint *gp, TimeStep *tStep)
150  { OOFEM_ERROR("calling of this function is not allowed"); }
152  { OOFEM_ERROR("calling of this funciton is not allowed"); }
153  virtual void computeBodyLoadVectorAt(FloatArray &answer, Load *forLoad, TimeStep *tStep, ValueModeType mode);
154 
155 public:
156  virtual void computeStiffnessMatrix(FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep)
157  { OOFEM_ERROR("calling of this function is not allowed"); }
158  virtual void computeMassMatrix(FloatMatrix &answer, TimeStep *tStep)
159  { OOFEM_ERROR("calling of this function is not allowed"); }
160  virtual void giveInternalForcesVector(FloatArray &answer, TimeStep *tStep, int useUpdatedGpRecord)
161  { OOFEM_ERROR("calling of this function is not allowed"); }
162 };
163 } // end namespace oofem
164 #endif
InternalStateType
Type representing the physical meaning of element or constitutive model internal variable.
The element interface required by NodalAvergagingRecoveryModel.
std::unique_ptr< IntegrationRule > compositeIR
Element integraton rule (plate and membrane parts have their own integration rules) this one used to ...
Definition: tr_shell02.h:68
virtual void ZZErrorEstimatorI_computeLocalStress(FloatArray &answer, FloatArray &sig)
Returns stress vector in global c.s.
Definition: tr_shell02.C:351
integrationDomain
Used by integrator class to supply integration points for proper domain to be integrated (Area...
std::unique_ptr< TrPlanestressRotAllman3d > membrane
Pointer to membrane (plane stress) element.
Definition: tr_shell02.h:62
virtual contextIOResultType saveContext(DataStream &stream, ContextMode mode, void *obj=NULL)
Stores receiver state to output stream.
Definition: tr_shell02.C:309
Class and object Domain.
Definition: domain.h:115
#define _IFT_TR_SHELL02_Name
Definition: tr_shell02.h:47
virtual IntegrationRule * giveDefaultIntegrationRulePtr()
Access method for default integration rule.
Definition: tr_shell02.h:109
Element_Geometry_Type
Enumerative type used to classify element geometry Possible values are: EGT_point - point in space EG...
virtual void computeStiffnessMatrix(FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep)
Computes numerically stiffness matrix of receiver.
Definition: tr_shell02.h:156
The element interface required by ZZNodalRecoveryModel.
virtual double computeVolumeAround(GaussPoint *gp)
Returns volume related to given integration point.
Definition: tr_shell02.C:218
The purpose of DataStream abstract class is to allow to store/restore context to different streams...
Definition: datastream.h:54
virtual void giveDofManDofIDMask(int inode, IntArray &answer) const
Returns dofmanager dof mask for node.
Definition: tr_shell02.h:82
virtual void giveCharacteristicVector(FloatArray &answer, CharType mtrx, ValueModeType mode, TimeStep *tStep)
Computes characteristic vector of receiver of requested type in given time step.
Definition: tr_shell02.C:114
ValueModeType
Type representing the mode of UnknownType or CharType, or similar types.
Definition: valuemodetype.h:78
virtual void computeBodyLoadVectorAt(FloatArray &answer, Load *forLoad, TimeStep *tStep, ValueModeType mode)
Computes the load vector due to body load acting on receiver, at given time step. ...
Definition: tr_shell02.C:224
oofem::oofegGraphicContext gc[OOFEG_LAST_LAYER]
virtual integrationDomain giveIntegrationDomain() const
Returns integration domain for receiver, used to initialize integration point over receiver volume...
Definition: tr_shell02.h:112
virtual void computeNmatrixAt(const FloatArray &iLocCoord, FloatMatrix &)
Computes interpolation matrix for element unknowns.
Definition: tr_shell02.h:139
virtual contextIOResultType restoreContext(DataStream &stream, ContextMode mode, void *obj=NULL)
Restores the receiver state previously written in stream.
Definition: tr_shell02.C:325
virtual const char * giveClassName() const
Definition: tr_shell02.h:86
MaterialMode
Type representing material mode of integration point.
Definition: materialmode.h:89
virtual int computeGlobalCoordinates(FloatArray &answer, const FloatArray &lcoords)
Computes the global coordinates from given element&#39;s local coordinates.
Definition: tr_shell02.h:128
void setCrossSection(int csIndx)
Sets the cross section model of receiver.
Definition: tr_shell02.C:104
Class implementing an array of integers.
Definition: intarray.h:61
MatResponseMode
Describes the character of characteristic material matrix.
virtual void NodalAveragingRecoveryMI_computeNodalValue(FloatArray &answer, int node, InternalStateType type, TimeStep *tStep)
Computes the element value in given node.
Definition: tr_shell02.C:252
virtual int giveDefaultIntegrationRule() const
Returns id of default integration rule.
Definition: tr_shell02.h:108
virtual void computeStressVector(FloatArray &answer, const FloatArray &strain, GaussPoint *gp, TimeStep *tStep)
Computes the stress vector of receiver at given integration point, at time step tStep.
Definition: tr_shell02.h:149
Abstract base class representing integration rule.
virtual void computeGaussPoints()
Definition: tr_shell02.h:144
virtual void drawRawGeometry(oofegGraphicContext &gc, TimeStep *tStep)
Definition: tr_shell02.C:447
virtual void drawDeformedGeometry(oofegGraphicContext &gc, TimeStep *tStep, UnknownType type)
Definition: tr_shell02.C:481
virtual void SpatialLocalizerI_giveBBox(FloatArray &bb0, FloatArray &bb1)
Creates a bounding box of the nodes and checks if it includes the given coordinate.
Definition: tr_shell02.C:408
Class representing a general abstraction for finite element interpolation class.
Definition: feinterpol.h:132
virtual IRResultType initializeFrom(InputRecord *ir)
Initializes receiver according to object description stored in input record.
Definition: tr_shell02.C:65
Abstract base class for all "structural" finite elements.
virtual void computeBmatrixAt(GaussPoint *, FloatMatrix &, int=1, int=ALL_STRAINS)
Computes the geometrical matrix of receiver in given integration point.
Definition: tr_shell02.h:137
virtual IntegrationRule * giveIntegrationRule(int i)
Definition: tr_shell02.h:110
The element interface corresponding to ZZErrorEstimator.
static IntArray loc_membrane
Definition: tr_shell02.h:71
#define OOFEM_ERROR(...)
Definition: error.h:61
static IntArray loc_plate
Definition: tr_shell02.h:70
virtual int computeNumberOfDofs()
Computes or simply returns total number of element&#39;s local DOFs.
Definition: tr_shell02.h:81
UnknownType
Type representing particular unknown (its physical meaning).
Definition: unknowntype.h:55
virtual void updateYourself(TimeStep *tStep)
Updates element state after equilibrium in time step has been reached.
Definition: tr_shell02.C:191
virtual FEInterpolation * giveInterpolation() const
Definition: tr_shell02.h:79
virtual void printOutputAt(FILE *file, TimeStep *tStep)
Prints output of receiver to stream, for given time step.
Definition: tr_shell02.C:262
virtual Element_Geometry_Type giveGeometryType() const
Returns the element geometry type.
Definition: tr_shell02.h:111
virtual bool giveRotationMatrix(FloatMatrix &answer)
Transformation matrices updates rotation matrix between element-local and primary DOFs...
Definition: tr_shell02.C:150
virtual void giveCharacteristicMatrix(FloatMatrix &answer, CharType mtrx, TimeStep *tStep)
Computes characteristic matrix of receiver of requested type in given time step.
Definition: tr_shell02.C:132
#define ALL_STRAINS
virtual const char * giveInputRecordName() const
Definition: tr_shell02.h:85
virtual void postInitialize()
Performs post initialization steps.
Definition: tr_shell02.C:87
virtual MaterialMode giveMaterialMode()
Returns material mode for receiver integration points.
Definition: tr_shell02.h:113
virtual void updateInternalState(TimeStep *tStep)
Updates element state after equilibrium in time step has been reached.
Definition: tr_shell02.C:183
TR_SHELL02(int n, Domain *d)
Constructor.
Definition: tr_shell02.C:57
Class representing vector of real numbers.
Definition: floatarray.h:82
virtual Interface * giveInterface(InterfaceType it)
Interface requesting service.
Definition: tr_shell02.C:201
virtual void drawScalar(oofegGraphicContext &gc, TimeStep *tStep)
Definition: tr_shell02.C:515
Implementation of matrix containing floating point numbers.
Definition: floatmatrix.h:94
IRResultType
Type defining the return values of InputRecord reading operations.
Definition: irresulttype.h:47
virtual ~TR_SHELL02()
Destructor.
Definition: tr_shell02.h:77
CharType
Definition: chartype.h:87
Class representing the general Input Record.
Definition: inputrecord.h:101
virtual IntegrationRule * ZZErrorEstimatorI_giveIntegrationRule()
Returns element integration rule used to evaluate error.
Definition: tr_shell02.C:341
Class Interface.
Definition: interface.h:82
The spatial localizer element interface associated to spatial localizer.
long ContextMode
Context mode (mask), defining the type of information written/read to/from context.
Definition: contextmode.h:43
virtual void computeMassMatrix(FloatMatrix &answer, TimeStep *tStep)
Computes mass matrix of receiver.
Definition: tr_shell02.h:158
InterfaceType
Enumerative type, used to identify interface type.
Definition: interfacetype.h:43
Load is base abstract class for all loads.
Definition: load.h:61
the oofem namespace is to define a context or scope in which all oofem names are defined.
virtual void giveInternalForcesVector(FloatArray &answer, TimeStep *tStep, int useUpdatedGpRecord)
Returns equivalent nodal forces vectors.
Definition: tr_shell02.h:160
virtual int giveIPValue(FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep)
Returns the integration point corresponding value in full form.
Definition: tr_shell02.C:230
virtual bool computeLocalCoordinates(FloatArray &answer, const FloatArray &gcoords)
Computes the element local coordinates from given global coordinates.
Definition: tr_shell02.h:132
Class representing integration point in finite element program.
Definition: gausspoint.h:93
Class representing solution step.
Definition: timestep.h:80
virtual void computeConstitutiveMatrixAt(FloatMatrix &answer, MatResponseMode rMode, GaussPoint *gp, TimeStep *tStep)
Computes constitutive matrix of receiver.
Definition: tr_shell02.h:151
std::unique_ptr< DKTPlate3d > plate
Pointer to plate element.
Definition: tr_shell02.h:60
void updateLocalNumbering(EntityRenumberingFunctor &f)
Local renumbering support.
Definition: tr_shell02.C:97
This class implements an triangular three-node shell finite element, composed of dkt3d and trplanestr...
Definition: tr_shell02.h:56

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:32 for OOFEM by doxygen 1.8.11 written by Dimitri van Heesch, © 1997-2011