OOFEM  2.4
OOFEM.org - Object Oriented Finite Element Solver
rcsdnl.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 // *********************************************************************************************
36 // *** CLASS NONLOCAL ROTATING SMEARED CRACK MODEL WITH TRANSITION TO SCALAR DAMAGE ************
37 // *********************************************************************************************
38 
39 #ifndef rcsdnl_h
40 #define rcsdnl_h
41 
42 #include "rcsde.h"
44 #include "../sm/Materials/structuralmaterial.h"
46 
47 #define _IFT_RCSDNLMaterial_Name "rcsdnl"
48 #define _IFT_RCSDNLMaterial_ft "ft"
49 #define _IFT_RCSDNLMaterial_sdtransitioncoeff "sdtransitioncoeff"
50 #define _IFT_RCSDNLMaterial_sdtransitioncoeff2 "sdtransitioncoeff2"
51 #define _IFT_RCSDNLMaterial_r "r"
52 #define _IFT_RCSDNLMaterial_ef "ef"
53 #define _IFT_RCSDNLMaterial_gf "gf"
54 
55 
56 namespace oofem {
61 {
62 protected:
64 
65 public:
66  RCSDNLMaterialStatus(int n, Domain * d, GaussPoint * g);
67  virtual ~RCSDNLMaterialStatus();
68 
69  virtual void printOutputAt(FILE *file, TimeStep *tStep);
70 
73  void setTempNonlocalStrainVector(FloatArray ls) { tempNonlocalStrainVector = std :: move(ls); }
74 
76  void setLocalStrainVectorForAverage(FloatArray ls) { localStrainVectorForAverage = std :: move(ls); }
77 
78  // definition
79  virtual const char *giveInputRecordName() const { return _IFT_RCSDNLMaterial_Name; }
80  virtual const char *giveClassName() const { return "RCSDNLMaterialStatus"; }
81 
82  virtual void initTempStatus();
83  virtual void updateYourself(TimeStep *);
84 
85  // saves current context(state) into stream
86  virtual contextIOResultType saveContext(DataStream &stream, ContextMode mode, void *obj = NULL);
87  virtual contextIOResultType restoreContext(DataStream &stream, ContextMode mode, void *obj = NULL);
88 
90 };
91 
92 
100 {
101 protected:
110  double R;
115  double ef;
116 
117 public:
118  RCSDNLMaterial(int n, Domain * d);
119  virtual ~RCSDNLMaterial();
120 
121  // identification and auxiliary functions
122  virtual const char *giveClassName() const { return "RCSDNLMaterial"; }
123 
125 
127 
128  virtual void giveRealStressVector(FloatArray &answer, GaussPoint *,
129  const FloatArray &, TimeStep *);
130 
141  virtual void updateBeforeNonlocAverage(const FloatArray &strainVector, GaussPoint *gp, TimeStep *tStep);
142 
143  virtual double computeWeightFunction(const FloatArray &src, const FloatArray &coord);
144  virtual int hasBoundedSupport() { return 1; }
148  virtual void giveSupportRadius(double &radius) { radius = this->R; }
149 
150  virtual int packUnknowns(DataStream &buff, TimeStep *tStep, GaussPoint *ip);
151  virtual int unpackAndUpdateUnknowns(DataStream &buff, TimeStep *tStep, GaussPoint *ip);
152  virtual int estimatePackSize(DataStream &buff, GaussPoint *ip);
153 
155 
156 protected:
157  virtual double giveCharacteristicElementLength(GaussPoint *gp, const FloatArray &) { return 1.0; }
158  virtual double giveMinCrackStrainsForFullyOpenCrack(GaussPoint *gp, int i);
159  virtual double computeStrength(GaussPoint *, double) { return this->Ft; }
160 };
161 } // end namespace oofem
162 #endif // rcsdnl_h
Abstract base class for all nonlocal structural materials.
virtual contextIOResultType saveContext(DataStream &stream, ContextMode mode, void *obj=NULL)
Stores receiver state to output stream.
Definition: rcsdnl.C:506
virtual void updateYourself(TimeStep *)
Update equilibrium history variables according to temp-variables.
Definition: rcsdnl.C:493
FloatArray localStrainVectorForAverage
Definition: rcsdnl.h:63
virtual void giveSupportRadius(double &radius)
Determines the width (radius) of limited support of weighting function.
Definition: rcsdnl.h:148
GaussPoint * gp
Associated integration point.
Class and object Domain.
Definition: domain.h:115
const FloatArray & giveNonlocalStrainVector()
Definition: rcsdnl.h:71
FloatArray nonlocalStrainVector
Definition: rcsdnl.h:63
Domain * domain
Link to domain object, useful for communicating with other FEM components.
Definition: femcmpnn.h:82
The purpose of DataStream abstract class is to allow to store/restore context to different streams...
Definition: datastream.h:54
virtual int hasBoundedSupport()
Determines, whether receiver has bounded weighting function (limited support).
Definition: rcsdnl.h:144
virtual double giveCharacteristicElementLength(GaussPoint *gp, const FloatArray &)
Definition: rcsdnl.h:157
const FloatArray & giveLocalStrainVectorForAverage()
Definition: rcsdnl.h:75
This class implements associated Material Status to RCSDNLMaterial.
Definition: rcsdnl.h:60
virtual ~RCSDNLMaterialStatus()
Definition: rcsdnl.C:449
virtual void initTempStatus()
Initializes the temporary internal variables, describing the current state according to previously re...
Definition: rcsdnl.C:472
virtual Interface * giveInterface(InterfaceType)
Interface requesting service.
Definition: rcsdnl.C:552
This class implements associated Material Status to RCSDEMaterial.
Definition: rcsde.h:58
virtual void printOutputAt(FILE *file, TimeStep *tStep)
Print receiver's output to given stream.
Definition: rcsdnl.C:454
This class implements a Rotating Crack Model with transition to scalar damage for fracture in smeared...
Definition: rcsde.h:118
double ef
Strain at complete failure.
Definition: rcsdnl.h:115
virtual IRResultType initializeFrom(InputRecord *ir)
Initializes receiver according to object description stored in input record.
Definition: matstatus.h:140
void setLocalStrainVectorForAverage(FloatArray ls)
Definition: rcsdnl.h:76
RCSDNLMaterialStatus(int n, Domain *d, GaussPoint *g)
Definition: rcsdnl.C:439
FloatArray tempNonlocalStrainVector
Definition: rcsdnl.h:63
virtual const char * giveClassName() const
Definition: rcsdnl.h:122
Abstract base class representing a material status information.
Definition: matstatus.h:84
Class representing vector of real numbers.
Definition: floatarray.h:82
FloatArray strainVector
Equilibrated strain vector in reduced form.
Definition: structuralms.h:69
IRResultType
Type defining the return values of InputRecord reading operations.
Definition: irresulttype.h:47
This class implements a Nonlocal Rotating Crack Model with transition to scalar damage for fracture i...
Definition: rcsdnl.h:99
virtual MaterialStatus * CreateStatus(GaussPoint *gp) const
Creates new copy of associated status and inserts it into given integration point.
Definition: rcsdnl.h:154
Class representing the general Input Record.
Definition: inputrecord.h:101
Base class for all nonlocal structural material statuses.
Class Interface.
Definition: interface.h:82
double SDTransitionCoeff2
Nondimensional parameter controlling the transition between rc and sd model, with respect to shear st...
Definition: rcsdnl.h:106
long ContextMode
Context mode (mask), defining the type of information written/read to/from context.
Definition: contextmode.h:43
virtual const char * giveInputRecordName() const
Definition: rcsdnl.h:79
InterfaceType
Enumerative type, used to identify interface type.
Definition: interfacetype.h:43
the oofem namespace is to define a context or scope in which all oofem names are defined.
void setTempNonlocalStrainVector(FloatArray ls)
Definition: rcsdnl.h:73
virtual const char * giveClassName() const
Definition: rcsdnl.h:80
Class representing integration point in finite element program.
Definition: gausspoint.h:93
virtual contextIOResultType restoreContext(DataStream &stream, ContextMode mode, void *obj=NULL)
Restores the receiver state previously written in stream.
Definition: rcsdnl.C:529
const FloatArray & giveTempNonlocalStrainVector()
Definition: rcsdnl.h:72
Class representing solution step.
Definition: timestep.h:80
virtual double computeStrength(GaussPoint *, double)
Definition: rcsdnl.h:159
#define _IFT_RCSDNLMaterial_Name
Definition: rcsdnl.h:47
double R
Interaction radius, related to the nonlocal characteristic length of material.
Definition: rcsdnl.h:110

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