OOFEM  2.4
OOFEM.org - Object Oriented Finite Element Solver
concretefcm.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 - 2016 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 concretefcm_h
36 #define concretefcm_h
37 
38 #include "../fcm.h"
39 #include "randommaterialext.h"
40 
42 
43 #define _IFT_ConcreteFCM_Name "concretefcm"
44 #define _IFT_ConcreteFCM_softType "softtype"
45 #define _IFT_ConcreteFCM_shearType "sheartype"
46 #define _IFT_ConcreteFCM_shearStrengthType "shearstrengthtype"
47 #define _IFT_ConcreteFCM_gf "gf"
48 #define _IFT_ConcreteFCM_ft "ft"
49 #define _IFT_ConcreteFCM_beta "beta"
50 #define _IFT_ConcreteFCM_sf "sf"
51 #define _IFT_ConcreteFCM_fc "fc"
52 #define _IFT_ConcreteFCM_ag "ag"
53 #define _IFT_ConcreteFCM_lengthScale "lengthscale"
54 #define _IFT_ConcreteFCM_soft_w "soft_w"
55 #define _IFT_ConcreteFCM_soft_function_w "soft(w)"
56 #define _IFT_ConcreteFCM_soft_eps "soft_eps"
57 #define _IFT_ConcreteFCM_soft_function_eps "soft(eps)"
58 #define _IFT_ConcreteFCM_beta_w "beta_w"
59 #define _IFT_ConcreteFCM_beta_function "beta(w)"
60 #define _IFT_ConcreteFCM_H "h"
61 #define _IFT_ConcreteFCM_eps_f "eps_f"
62 
63 
64 namespace oofem {
69 {
70 public:
71  ConcreteFCMStatus(int n, Domain *d, GaussPoint *g);
72  virtual ~ConcreteFCMStatus();
74  virtual void printOutputAt(FILE *file, TimeStep *tStep);
75 
76  // definition
77  virtual const char *giveClassName() const { return "ConcreteFCMStatus"; }
78 
79  virtual void initTempStatus();
80  virtual void updateYourself(TimeStep *tStep);
81 
83 
84  // saves current context(state) into stream
85  virtual contextIOResultType saveContext(DataStream &stream, ContextMode mode, void *obj = NULL);
86  virtual contextIOResultType restoreContext(DataStream &stream, ContextMode mode, void *obj = NULL);
87 };
88 
89 
90 
99 {
100 public:
101  ConcreteFCM(int n, Domain *d);
102  virtual ~ConcreteFCM() {
103  delete linearElasticMaterial;
104  }
105 
106  // identification and auxiliary functions
108  virtual int hasNonLinearBehaviour() { return 1; }
109  virtual const char *giveClassName() const { return "ConcreteFCM"; }
110  virtual const char *giveInputRecordName() const { return _IFT_ConcreteFCM_Name; }
111 
112  virtual MaterialStatus *CreateStatus(GaussPoint *gp) const { return new ConcreteFCMStatus(1, domain, gp); }
113 
114  virtual double give(int aProperty, GaussPoint *gp);
115 
116  virtual int giveIPValue(FloatArray &answer, GaussPoint *gp, InternalStateType type, TimeStep *tStep);
117 
118  virtual MaterialStatus *giveStatus(GaussPoint *gp) const;
119 
120 
121 protected:
122 
124  double Gf;
126  double Ft;
128  double beta;
130  double sf;
131 
132  // 3 parameters for collins aggregate interlock:
134  double fc;
136  double ag;
138  double lengthScale;
139 
141  FloatArray soft_w, soft_function_w;
145  FloatArray beta_w, beta_function;
146 
148  double H;
150  double eps_f;
151 
153  virtual double giveTensileStrength(GaussPoint *gp) { return this->give(ft_strength, gp); }
154 
156  virtual double giveFractureEnergy(GaussPoint *gp) { return this->give(gf_ID, gp); }
157 
159  virtual double giveCrackingModulus(MatResponseMode rMode, GaussPoint *gp, int i);
160 
162  virtual double computeEffectiveShearModulus(GaussPoint *gp, int i);
163 
165  virtual double computeD2ModulusForCrack(GaussPoint *gp, int icrack);
166 
168  virtual double giveNormalCrackingStress(GaussPoint *gp, double eps_cr, int i);
169 
171  virtual double maxShearStress(GaussPoint *gp, int i);
172 
174  virtual void checkSnapBack(GaussPoint *gp, int crack);
175 
177  enum SofteningType { ST_NONE, ST_Exponential, ST_Linear, ST_Hordijk, ST_UserDefinedCrack, ST_LinearHardeningStrain, ST_UserDefinedStrain, ST_Unknown };
179 
181  enum ShearRetentionType { SHR_NONE, SHR_Const_ShearRetFactor, SHR_Const_ShearFactorCoeff, SHR_UserDefined_ShearRetFactor, SHR_Unknown };
183 
185  enum ShearStrengthType { SHS_NONE, SHS_Const_Ft, SHS_Collins_Interlock, SHS_Unknown };
187 };
188 } // end namespace oofem
189 #endif // concretefcm_h
InternalStateType
Type representing the physical meaning of element or constitutive model internal variable.
double H
hardening modulus
Definition: concretefcm.h:148
GaussPoint * gp
Associated integration point.
Class and object Domain.
Definition: domain.h:115
Domain * domain
Link to domain object, useful for communicating with other FEM components.
Definition: femcmpnn.h:82
This class implements associated Material Status to FCMMaterial (fixed crack material).
Definition: fcm.h:68
#define ft_strength
Definition: matconst.h:91
The purpose of DataStream abstract class is to allow to store/restore context to different streams...
Definition: datastream.h:54
This class implements a Fixed Crack Model for fracture (after initiation the crack directions cannot ...
Definition: fcm.h:197
This class manages the status of ConcreteFCM.
Definition: concretefcm.h:68
virtual Interface * giveInterface(InterfaceType it)
Interface requesting service.
Definition: concretefcm.C:857
MatResponseMode
Describes the character of characteristic material matrix.
virtual int hasNonLinearBehaviour()
Returns nonzero if receiver is non linear.
Definition: concretefcm.h:108
ShearStrengthType shearStrengthType
Definition: concretefcm.h:186
SofteningType softType
Definition: concretefcm.h:178
virtual contextIOResultType saveContext(DataStream &stream, ContextMode mode, void *obj=NULL)
saves current context(state) into stream
Definition: concretefcm.C:869
virtual ~ConcreteFCM()
Definition: concretefcm.h:102
ConcreteFCMStatus(int n, Domain *d, GaussPoint *g)
Definition: concretefcm.C:814
double beta
shear retention factor
Definition: concretefcm.h:128
virtual double giveFractureEnergy(GaussPoint *gp)
returns fracture energy (can be random)
Definition: concretefcm.h:156
virtual IRResultType initializeFrom(InputRecord *ir)
Initializes receiver according to object description stored in input record.
Definition: matstatus.h:140
virtual const char * giveClassName() const
Definition: concretefcm.h:109
This class implements a ConcreteFCM material in a finite element problem.
Definition: concretefcm.h:98
#define gf_ID
Definition: matconst.h:85
SofteningType
type of post-peak behavior in the normal direction to the crack plane
Definition: concretefcm.h:177
Abstract base class for all random materials.
virtual MaterialStatus * CreateStatus(GaussPoint *gp) const
Creates new copy of associated status and inserts it into given integration point.
Definition: concretefcm.h:112
FloatArray beta_w
user-defined shear retention factor (with respect to crack opening)
Definition: concretefcm.h:145
ShearStrengthType
defines the maximum value of shear stress
Definition: concretefcm.h:185
double ag
Collins' aggregate interlock: aggregate diameter in appropriate units (same as FE mesh) ...
Definition: concretefcm.h:136
ShearRetentionType shearType
Definition: concretefcm.h:182
Abstract base class for all random constitutive model statuses.
double eps_f
strain at failure
Definition: concretefcm.h:150
double Ft
Tensile strenght.
Definition: concretefcm.h:126
Abstract base class representing a material status information.
Definition: matstatus.h:84
Class representing vector of real numbers.
Definition: floatarray.h:82
IRResultType
Type defining the return values of InputRecord reading operations.
Definition: irresulttype.h:47
virtual void printOutputAt(FILE *file, TimeStep *tStep)
Writes information into the output file.
Definition: concretefcm.C:825
virtual void initTempStatus()
initializes temporary status
Definition: concretefcm.C:835
double lengthScale
Collins' aggregate interlock: 1 for meter, 1000 for analysis in mm.
Definition: concretefcm.h:138
#define _IFT_ConcreteFCM_Name
Definition: concretefcm.h:43
Class representing the general Input Record.
Definition: inputrecord.h:101
Class Interface.
Definition: interface.h:82
ShearRetentionType
type of reduction of the shear stiffness caused by cracking
Definition: concretefcm.h:181
virtual contextIOResultType restoreContext(DataStream &stream, ContextMode mode, void *obj=NULL)
restores context(state) from stream
Definition: concretefcm.C:886
long ContextMode
Context mode (mask), defining the type of information written/read to/from context.
Definition: contextmode.h:43
InterfaceType
Enumerative type, used to identify interface type.
Definition: interfacetype.h:43
virtual void updateYourself(TimeStep *tStep)
replaces equilibrated values with temporary values
Definition: concretefcm.C:846
virtual const char * giveClassName() const
Definition: concretefcm.h:77
the oofem namespace is to define a context or scope in which all oofem names are defined.
Class representing integration point in finite element program.
Definition: gausspoint.h:93
virtual const char * giveInputRecordName() const
Definition: concretefcm.h:110
Class representing solution step.
Definition: timestep.h:80
double sf
shear factor
Definition: concretefcm.h:130
FloatArray soft_w
user-defined softening (traction-COD)
Definition: concretefcm.h:141
double fc
Collins' aggregate interlock: compressive strength in MPa.
Definition: concretefcm.h:134
FloatArray soft_function_eps
Definition: concretefcm.h:143
double Gf
Fracture energy.
Definition: concretefcm.h:124
virtual double giveTensileStrength(GaussPoint *gp)
returns tensile strength (can be random)
Definition: concretefcm.h:153

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