OOFEM  2.4
OOFEM.org - Object Oriented Finite Element Solver
qplanestressgrad.C
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 #include "../sm/Elements/PlaneStress/qplanestressgrad.h"
36 #include "fei2dquadlin.h"
37 #include "gausspoint.h"
38 #include "gaussintegrationrule.h"
39 #include "floatmatrix.h"
40 #include "floatarray.h"
41 #include "intarray.h"
42 #include "crosssection.h"
43 #include "classfactory.h"
44 
45 namespace oofem {
46 REGISTER_Element(QPlaneStressGrad);
47 
48 FEI2dQuadLin QPlaneStressGrad :: interpolation_lin(1, 2);
49 
51  // Constructor.
52 {
53  nPrimNodes = 8;
54  nPrimVars = 2;
55  nSecNodes = 4;
56  nSecVars = 1;
61 }
62 
63 
64 void
66 {
67  if ( inode <= 4 ) {
68  answer = {D_u, D_v, G_0};
69  } else {
70  answer = {D_u, D_v};
71  }
72 }
73 
74 
77 {
79 }
80 
81 
82 void
84 {
85  if ( integrationRulesArray.size() == 0 ) {
86  integrationRulesArray.resize( 1 );
87  integrationRulesArray [ 0 ].reset( new GaussIntegrationRule(1, this, 1, 3) );
89  }
90 }
91 
92 void
94 {
96 }
97 
98 void
100 {
101  FloatMatrix dnx;
103  answer.beTranspositionOf(dnx);
104 }
105 }
CrossSection * giveCrossSection()
Definition: element.C:495
Class and object Domain.
Definition: domain.h:115
virtual void computeBkappaMatrixAt(GaussPoint *gp, FloatMatrix &answer)
QPlaneStressGrad(int n, Domain *d)
Class implementing an array of integers.
Definition: intarray.h:61
Abstract class for gradient formulation of coupled damage-plasticity model(GradDp).
Definition: graddpelement.h:48
virtual void computeGaussPoints()
Initializes the array of integration rules member variable.
This class implements an Quadratic isoparametric eight-node quadrilateral plane- stress elasticity fi...
Definition: qplanstrss.h:52
REGISTER_Element(LSpace)
virtual IRResultType initializeFrom(InputRecord *ir)
Wrapper around element definition to provide FEICellGeometry interface.
Definition: feinterpol.h:95
virtual int setupIntegrationPoints(IntegrationRule &irule, int npoints, Element *element)
Sets up integration rule for the given element.
Definition: crosssection.C:54
int numberOfGaussPoints
Number of integration points as specified by nip.
Definition: element.h:188
Class representing vector of real numbers.
Definition: floatarray.h:82
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
Class representing the general Input Record.
Definition: inputrecord.h:101
void beTranspositionOf(const FloatMatrix &src)
Assigns to the receiver the transposition of parameter.
Definition: floatmatrix.C:323
std::vector< std::unique_ptr< IntegrationRule > > integrationRulesArray
List of integration rules of receiver (each integration rule contains associated integration points a...
Definition: element.h:170
the oofem namespace is to define a context or scope in which all oofem names are defined.
virtual IRResultType initializeFrom(InputRecord *ir)
Initializes receiver according to object description stored in input record.
virtual double evaldNdx(FloatMatrix &answer, const FloatArray &lcoords, const FEICellGeometry &cellgeo)
Evaluates the matrix of derivatives of interpolation functions (shape functions) at given point...
Definition: fei2dquadlin.C:75
Class representing integration point in finite element program.
Definition: gausspoint.h:93
static FEI2dQuadLin interpolation_lin
virtual void evalN(FloatArray &answer, const FloatArray &lcoords, const FEICellGeometry &cellgeo)
Evaluates the array of interpolation functions (shape functions) at given point.
Definition: fei2dquadlin.C:59
const FloatArray & giveNaturalCoordinates()
Returns coordinate array of receiver.
Definition: gausspoint.h:138
virtual void giveDofManDofIDMask(int inode, IntArray &answer) const
Returns dofmanager dof mask for node.
virtual void computeNkappaMatrixAt(GaussPoint *gp, FloatArray &answer)
Class representing Gaussian-quadrature integration rule.

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