OOFEM  2.4
OOFEM.org - Object Oriented Finite Element Solver
b3solidmat.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 b3solidmat_h
36 #define b3solidmat_h
37 
38 #include "kelvinChM.h"
39 
41 
42 #define _IFT_B3SolidMaterial_Name "b3solidmat"
43 #define _IFT_B3SolidMaterial_emodulimode "emodulimode"
44 #define _IFT_B3SolidMaterial_microprestress "microprestress"
45 #define _IFT_B3SolidMaterial_c0 "c0"
46 #define _IFT_B3SolidMaterial_c1 "c1"
47 #define _IFT_B3SolidMaterial_ksh "ksh"
48 #define _IFT_B3SolidMaterial_finalhumidity "finalhumidity"
49 #define _IFT_B3SolidMaterial_initialhumidity "initialhumidity"
50 #define _IFT_B3SolidMaterial_lambda0 "lambda0"
51 
52 
53 namespace oofem {
58 {
59 protected:
63 
64 public:
65  B3SolidMaterialStatus(int n, Domain *d, GaussPoint *g, int nunits);
66  virtual ~B3SolidMaterialStatus() { }
67 
68  virtual void updateYourself(TimeStep *tStep);
69 
70  virtual contextIOResultType saveContext(DataStream &stream, ContextMode mode, void *obj = NULL);
71  virtual contextIOResultType restoreContext(DataStream &stream, ContextMode mode, void *obj = NULL);
72 
73  double giveMPS() const { return microprestress_old; }
74  void setMPS(double src) { microprestress_new = src; }
75 
76  // definition
77  virtual const char *giveClassName() const { return "B3SolidMaterialStatus"; }
78 };
79 
80 
87 {
88 protected:
89  double t0;
90  double w, E28, q1, q2, q3, q4, q5; // predicted data
91 
93  double lambda0;
94 
95  enum b3ShModeType { B3_NoShrinkage, B3_AverageShrinkage, B3_PointShrinkage } shMode;
97  double EpsSinf, kt, ks, vs, hum;
99  double es0, r, rprime, at;
100  // Additional parameters for sorption isotherm (used to compute relative humidity from water content)
101  double w_h;
102  double n;
103  double a;
104  double EspringVal;
105 
115  double c0;
116  double c1;
117  double tS0;
118  double kSh;
119 
120 
121 public:
123  shMode = B3_NoShrinkage;
124  }
125  virtual ~B3SolidMaterial() { }
126 
127  virtual void giveRealStressVector(FloatArray &answer, GaussPoint *gp,
128  const FloatArray &reducedStrain, TimeStep *tStep);
129  //virtual void updateYourself(GaussPoint *gp, TimeStep *tStep);
130 
131  virtual void giveShrinkageStrainVector(FloatArray &answer, GaussPoint *gp, TimeStep *tStep, ValueModeType mode);
132 
133  virtual const char *giveClassName() const { return "B3SolidMaterial"; }
134  virtual const char *giveInputRecordName() const { return _IFT_B3SolidMaterial_Name; }
136 
137  virtual MaterialStatus *CreateStatus(GaussPoint *gp) const;
138 
140  virtual double computeCreepFunction(double t, double t_prime, GaussPoint *gp, TimeStep *tStep);
141 
142 protected:
143  virtual int hasIncrementalShrinkageFormulation() { return 1; }
144 
145  void computeTotalAverageShrinkageStrainVector(FloatArray &answer, GaussPoint *gp, TimeStep *tStep);
146 
148  void computePointShrinkageStrainVector(FloatArray &answer, GaussPoint *gp, TimeStep *tStep);
149 
150  void predictParametersFrom(double, double, double, double, double, double, double);
151 
153  double computeSolidifiedVolume(TimeStep *tStep);
154 
156  double computeFlowTermViscosity(GaussPoint *gp, TimeStep *tStep);
157 
158  double inverse_sorption_isotherm(double w);
159 
161  virtual void computeCharCoefficients(FloatArray &answer, double tPrime, GaussPoint *gp, TimeStep *tStep);
162 
164  virtual void updateEparModuli(double tPrime, GaussPoint *gp, TimeStep *tStep);
165 
166  virtual void computeCharTimes();
167 
168  virtual double giveEModulus(GaussPoint *gp, TimeStep *tStep);
169 
170  virtual void giveEigenStrainVector(FloatArray &answer, GaussPoint *gp, TimeStep *tStep, ValueModeType mode);
171 
179  double computeMicroPrestress(GaussPoint *gp, TimeStep *tStep, int option);
180 
182  double giveInitMicroPrestress(void);
183 
185  double giveHumidity(GaussPoint *gp, TimeStep *tStep);
186 
188  double giveHumidityIncrement(GaussPoint *gp, TimeStep *tStep);
189 };
190 } // end namespace oofem
191 #endif // b3solidmat_h
virtual const char * giveInputRecordName() const
Definition: b3solidmat.h:134
GaussPoint * gp
Associated integration point.
Class and object Domain.
Definition: domain.h:115
This class implements a solidifying Kelvin chain model describing a viscoelastic material.
Definition: kelvinChM.h:67
virtual const char * giveClassName() const
Definition: b3solidmat.h:133
The purpose of DataStream abstract class is to allow to store/restore context to different streams...
Definition: datastream.h:54
void setMPS(double src)
Definition: b3solidmat.h:74
virtual ~B3SolidMaterial()
Definition: b3solidmat.h:125
ValueModeType
Type representing the mode of UnknownType or CharType, or similar types.
Definition: valuemodetype.h:78
double microprestress_old
Microprestresses.
Definition: b3solidmat.h:61
This class implements the B3 model for concrete creep and shrinkage based on the solidification theor...
Definition: b3solidmat.h:86
double giveMPS() const
Definition: b3solidmat.h:73
int MicroPrestress
If 1, computation exploiting Microprestress solidification theory is done.
Definition: b3solidmat.h:114
virtual contextIOResultType restoreContext(DataStream &stream, ContextMode mode, void *obj=NULL)
Restores the receiver state previously written in stream.
Definition: b3solidmat.C:846
B3SolidMaterial(int n, Domain *d)
Definition: b3solidmat.h:122
double c1
MPS constant c1 (=C1*R*T/M)
Definition: b3solidmat.h:116
virtual void updateYourself(TimeStep *tStep)
Update equilibrium history variables according to temp-variables.
Definition: b3solidmat.C:816
int EmoduliMode
If 0, analysis of retardation spectrum is used for evaluation of Kelvin units moduli (default)...
Definition: b3solidmat.h:109
virtual int hasIncrementalShrinkageFormulation()
If only incremental shrinkage strain formulation is provided, then total shrinkage strain must be tra...
Definition: b3solidmat.h:143
double kSh
MPS shrinkage parameter. Either this or inithum and finalhum must be given in input record...
Definition: b3solidmat.h:118
virtual IRResultType initializeFrom(InputRecord *ir)
Initializes receiver according to object description stored in input record.
Definition: matstatus.h:140
double a
Constant (obtained from experiments) A [Pedersen, 1990].
Definition: b3solidmat.h:103
This class implements associated Material Status to KelvinChainMaterial.
Definition: kelvinChM.h:44
This class implements associated Material Status to B3SolidMaterial.
Definition: b3solidmat.h:57
double n
Constant-exponent (obtained from experiments) n [Pedersen, 1990].
Definition: b3solidmat.h:102
Abstract base class representing a material status information.
Definition: matstatus.h:84
Class representing vector of real numbers.
Definition: floatarray.h:82
double EspringVal
elastic modulus of the aging spring (first member of Kelvin chain if retardation spectrum is used) ...
Definition: b3solidmat.h:104
double lambda0
constant equal to one day in time units of analysis (eg. 86400 if the analysis runs in seconds) ...
Definition: b3solidmat.h:93
IRResultType
Type defining the return values of InputRecord reading operations.
Definition: irresulttype.h:47
double c0
MPS constant c0 [MPa^-1 * day^-1].
Definition: b3solidmat.h:115
Class representing the general Input Record.
Definition: inputrecord.h:101
FloatArray * giveShrinkageStrainVector()
Definition: rheoChM.h:104
long ContextMode
Context mode (mask), defining the type of information written/read to/from context.
Definition: contextmode.h:43
virtual contextIOResultType saveContext(DataStream &stream, ContextMode mode, void *obj=NULL)
Stores receiver state to output stream.
Definition: b3solidmat.C:826
the oofem namespace is to define a context or scope in which all oofem names are defined.
virtual const char * giveClassName() const
Definition: b3solidmat.h:77
double w_h
Constant water content (obtained from experiments) w_h [Pedersen, 1990].
Definition: b3solidmat.h:101
Class representing integration point in finite element program.
Definition: gausspoint.h:93
Class representing solution step.
Definition: timestep.h:80
#define _IFT_B3SolidMaterial_Name
Definition: b3solidmat.h:42
B3SolidMaterialStatus(int n, Domain *d, GaussPoint *g, int nunits)
Definition: b3solidmat.C:812
double tS0
MPS tS0 - necessary for the initial value of microprestress (age when the load is applied) ...
Definition: b3solidmat.h:117

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