OOFEM  2.4
OOFEM.org - Object Oriented Finite Element Solver
intelline1PF.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 intelline1pf_h
36 #define intelline1pf_h
37 
40 #define _IFT_IntElLine1PF_Name "IntElLine1PF"
41 #define _IFT_IntElLine1PF_prescribedDamage "prescribeddamage"
42 
43 namespace oofem {
44 class FEI2dLineLin;
45 
53 {
54 protected:
56 
57 public:
58  IntElLine1PF(int n, Domain *d);
59  virtual ~IntElLine1PF() { }
60 
61  virtual FEInterpolation *giveInterpolation() const;
62 
63  virtual int computeNumberOfDofs() { return 12; }
64  virtual void giveDofManDofIDMask(int inode, IntArray &answer) const;
65 
66  virtual double computeAreaAround(GaussPoint *gp);
67  virtual void computeTransformationMatrixAt(GaussPoint *gp, FloatMatrix &answer);
68  virtual void computeCovarBaseVectorAt(GaussPoint *gp, FloatArray &G);
69 
70  virtual int testElementExtension(ElementExtension ext) { return 0; }
71 
72  //virtual Interface *giveInterface(InterfaceType) { return NULL; }
73 
74  // definition & identification
75  virtual const char *giveInputRecordName() const { return _IFT_IntElLine1PF_Name; }
76  virtual const char *giveClassName() const { return "IntElLine1PF"; }
78 
79  virtual void giveEngTraction(FloatArray &answer, GaussPoint *gp, const FloatArray &jump, TimeStep *tStep)
80  {
81  this->giveInterfaceCrossSection()->giveEngTraction_2d(answer, gp, jump, tStep);
82  }
83 
85  {
86  this->giveInterfaceCrossSection()->give2dStiffnessMatrix_Eng(answer, rMode, ip, tStep);
87  }
88 
89  virtual void computeStiffnessMatrix(FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep);
90 
91 
92 
93  // PF
94  int giveNumberOfDofManagers() const { return 4; }
96  virtual void giveDofManDofIDMask_u(IntArray &answer);
97  virtual void giveDofManDofIDMask_d(IntArray &answer);
98 
102  //virtual void computeStiffnessMatrix_du(FloatMatrix &, MatResponseMode, TimeStep *);
103 
104  virtual void giveInternalForcesVector(FloatArray &answer, TimeStep *tStep, int useUpdatedGpRecord = 0);
105  virtual void giveInternalForcesVectorUD(FloatArray &fu, FloatArray &fd, TimeStep *tStep, int useUpdatedGpRecord = 0);
106  virtual double computeDamageAt(GaussPoint *gp, ValueModeType valueMode, TimeStep *stepN);
107  virtual void computeLocationArrayOfDofIDs( const IntArray &dofIdArray, IntArray &answer );
108  virtual void computeBd_vectorAt(GaussPoint *gp, FloatArray &B);
109  virtual void computeNd_vectorAt(const FloatArray &lCoords, FloatArray &N);
110  virtual double computeFreeEnergy(GaussPoint *gp, TimeStep *tStep);
111 
112  double computeOldG(GaussPoint *gp, ValueModeType valueMode, TimeStep *stepN);
113  double neg_MaCauley(double par)
114  {
115  return 0.5 * ( abs(par) - par );
116  }
117 
118  double neg_MaCauleyPrime(double par)
119  {
120  return 0.5 * ( abs(par)/(par + 1.0e-12) - 1.0 ); // 0.5*(sign - 1) taken from Ragnars code
121  }
122 
123  double MaCauley(double par)
124  {
125  return 0.5 * ( abs(par) + par );
126  }
127 
128  double MaCauleyPrime(double par)
129  {
130  return 0.5 * ( abs(par)/(par + 1.0e-12) + 1.0 ); // 0.5*(sign + 1) taken from Ragnars code
131  }
132 
133 
134  void solveForLocalDamage(FloatMatrix &answer, TimeStep *tStep);
135  void computeGMatrix(FloatMatrix &answer, const double damage, GaussPoint *gp, ValueModeType valueMode, TimeStep *stepN);
136 
137 protected:
138 
142 
147 
148  virtual void computeNmatrixAt(GaussPoint *gp, FloatMatrix &answer);
149  virtual void computeGaussPoints();
150 
151  virtual int giveApproxOrder() { return 1; }
152  Element_Geometry_Type giveGeometryType() const { return EGT_quad_1_interface; };
153 };
154 } // end namespace oofem
155 #endif
virtual void computeCovarBaseVectorAt(GaussPoint *gp, FloatArray &G)
Definition: intelline1PF.C:102
virtual double computeDamageAt(GaussPoint *gp, ValueModeType valueMode, TimeStep *stepN)
Definition: intelline1PF.C:574
double neg_MaCauley(double par)
Definition: intelline1PF.h:113
Class and object Domain.
Definition: domain.h:115
Element_Geometry_Type
Enumerative type used to classify element geometry Possible values are: EGT_point - point in space EG...
virtual IRResultType initializeFrom(InputRecord *ir)
Initializes receiver according to object description stored in input record.
Definition: intelline1PF.C:133
void solveForLocalDamage(FloatMatrix &answer, TimeStep *tStep)
Definition: intelline1PF.C:448
virtual void computeStiffnessMatrix_ud(FloatMatrix &, MatResponseMode, TimeStep *)
Definition: intelline1PF.C:271
Element_Geometry_Type giveGeometryType() const
Returns the element geometry type.
Definition: intelline1PF.h:152
void computeGMatrix(FloatMatrix &answer, const double damage, GaussPoint *gp, ValueModeType valueMode, TimeStep *stepN)
Definition: intelline1PF.C:539
ValueModeType
Type representing the mode of UnknownType or CharType, or similar types.
Definition: valuemodetype.h:78
Abstract class for phase field formulation.
FloatArray deltaAlpha
Definition: intelline1PF.h:143
#define _IFT_IntElLine1PF_Name
Definition: intelline1PF.h:40
virtual void giveDofManDofIDMask_u(IntArray &answer)
Definition: intelline1PF.C:743
virtual ~IntElLine1PF()
Definition: intelline1PF.h:59
FloatArray unknownVectorD
Definition: intelline1PF.h:140
StructuralInterfaceElement * giveElement()
Definition: intelline1PF.h:95
double computeOldG(GaussPoint *gp, ValueModeType valueMode, TimeStep *stepN)
Definition: intelline1PF.C:732
virtual void computeNd_vectorAt(const FloatArray &lCoords, FloatArray &N)
Definition: intelline1PF.C:779
Class implementing an array of integers.
Definition: intarray.h:61
MatResponseMode
Describes the character of characteristic material matrix.
virtual FEInterpolation * giveInterpolation() const
Definition: intelline1PF.C:171
This class implements a two dimensional interface element.
Definition: intelline1PF.h:52
virtual void computeNmatrixAt(GaussPoint *gp, FloatMatrix &answer)
Computes modified interpolation matrix (N) for the element which multiplied with the unknowns vector ...
Definition: intelline1PF.C:72
int giveNumberOfDofManagers() const
Definition: intelline1PF.h:94
double neg_MaCauleyPrime(double par)
Definition: intelline1PF.h:118
virtual double computeAreaAround(GaussPoint *gp)
Definition: intelline1PF.C:119
Class representing a general abstraction for finite element interpolation class.
Definition: feinterpol.h:132
virtual void giveInternalForcesVectorUD(FloatArray &fu, FloatArray &fd, TimeStep *tStep, int useUpdatedGpRecord=0)
Definition: intelline1PF.C:626
FloatArray deltaUnknownVectorD
Definition: intelline1PF.h:141
virtual void computeBd_vectorAt(GaussPoint *gp, FloatArray &B)
Definition: intelline1PF.C:788
virtual int testElementExtension(ElementExtension ext)
Tests if the element implements required extension.
Definition: intelline1PF.h:70
ElementExtension
Type representing element extension.
virtual void giveStiffnessMatrix_Eng(FloatMatrix &answer, MatResponseMode rMode, IntegrationPoint *ip, TimeStep *tStep)
Definition: intelline1PF.h:84
virtual void giveEngTraction(FloatArray &answer, GaussPoint *gp, const FloatArray &jump, TimeStep *tStep)
Definition: intelline1PF.h:79
#define N(p, q)
Definition: mdm.C:367
Class representing a 2d line with linear interpolation.
Definition: fei2dlinelin.h:46
virtual void computeStiffnessMatrix_uu(FloatMatrix &, MatResponseMode, TimeStep *)
Definition: intelline1PF.C:218
double MaCauley(double par)
Definition: intelline1PF.h:123
Class representing vector of real numbers.
Definition: floatarray.h:82
Implementation of matrix containing floating point numbers.
Definition: floatmatrix.h:94
IntElLine1PF(int n, Domain *d)
Definition: intelline1PF.C:56
IRResultType
Type defining the return values of InputRecord reading operations.
Definition: irresulttype.h:47
virtual int computeNumberOfDofs()
Computes or simply returns total number of element's local DOFs.
Definition: intelline1PF.h:63
virtual void computeGaussPoints()
Initializes the array of integration rules member variable.
Definition: intelline1PF.C:91
void give2dStiffnessMatrix_Eng(FloatMatrix &answer, MatResponseMode rMode, GaussPoint *gp, TimeStep *tStep)
Class representing the general Input Record.
Definition: inputrecord.h:101
double MaCauleyPrime(double par)
Definition: intelline1PF.h:128
virtual void computeTransformationMatrixAt(GaussPoint *gp, FloatMatrix &answer)
Definition: intelline1PF.C:155
virtual const char * giveClassName() const
Definition: intelline1PF.h:76
virtual const char * giveInputRecordName() const
Definition: intelline1PF.h:75
virtual double computeFreeEnergy(GaussPoint *gp, TimeStep *tStep)
Definition: intelline1PF.C:711
virtual void giveInternalForcesVector(FloatArray &answer, TimeStep *tStep, int useUpdatedGpRecord=0)
Returns equivalent nodal forces vectors.
Definition: intelline1PF.C:599
virtual int giveApproxOrder()
Definition: intelline1PF.h:151
FloatArray unknownVectorU
Definition: intelline1PF.h:139
Abstract base class for all structural interface elements.
void giveEngTraction_2d(FloatArray &answer, GaussPoint *gp, const FloatArray &jump, TimeStep *tStep)
the oofem namespace is to define a context or scope in which all oofem names are defined.
static FEI2dLineLin interp
Definition: intelline1PF.h:55
virtual void computeLocationArrayOfDofIDs(const IntArray &dofIdArray, IntArray &answer)
Definition: intelline1PF.C:758
virtual void computeStiffnessMatrix_dd(FloatMatrix &, MatResponseMode, TimeStep *)
Definition: intelline1PF.C:352
Class representing integration point in finite element program.
Definition: gausspoint.h:93
virtual void giveDofManDofIDMask_d(IntArray &answer)
Definition: intelline1PF.C:751
Class representing solution step.
Definition: timestep.h:80
virtual void computeStiffnessMatrix(FloatMatrix &answer, MatResponseMode rMode, TimeStep *tStep)
Computes the stiffness/tangent matrix of receiver.
Definition: intelline1PF.C:180
virtual void giveDofManDofIDMask(int inode, IntArray &answer) const
Returns dofmanager dof mask for node.
Definition: intelline1PF.C:149
StructuralInterfaceCrossSection * giveInterfaceCrossSection()

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