XC Open source finite element analysis program
Steel03.h
1 //----------------------------------------------------------------------------
2 // XC program; finite element analysis code
3 // for structural analysis and design.
4 //
5 // Copyright (C) Luis Claudio Pérez Tato
6 //
7 // This program derives from OpenSees <http://opensees.berkeley.edu>
8 // developed by the «Pacific earthquake engineering research center».
9 //
10 // Except for the restrictions that may arise from the copyright
11 // of the original program (see copyright_opensees.txt)
12 // XC is free software: you can redistribute it and/or modify
13 // it under the terms of the GNU General Public License as published by
14 // the Free Software Foundation, either version 3 of the License, or
15 // (at your option) any later version.
16 //
17 // This software is distributed in the hope that it will be useful, but
18 // WITHOUT ANY WARRANTY; without even the implied warranty of
19 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
20 // GNU General Public License for more details.
21 //
22 //
23 // You should have received a copy of the GNU General Public License
24 // along with this program.
25 // If not, see <http://www.gnu.org/licenses/>.
26 //----------------------------------------------------------------------------
27 /* ****************************************************************** **
28 ** OpenSees - Open System for Earthquake Engineering Simulation **
29 ** Pacific Earthquake Engineering Research Center **
30 ** **
31 ** **
32 ** (C) Copyright 1999, The Regents of the University of California **
33 ** All Rights Reserved. **
34 ** **
35 ** Commercial use of this program without express permission of the **
36 ** University of California, Berkeley, is strictly prohibited. See **
37 ** file 'COPYRIGHT' in main directory for information on usage and **
38 ** redistribution, and for a DISCLAIMER OF ALL WARRANTIES. **
39 ** **
40 ** Developed by: **
41 ** Frank McKenna (fmckenna@ce.berkeley.edu) **
42 ** Gregory L. Fenves (fenves@ce.berkeley.edu) **
43 ** Filip C. Filippou (filippou@ce.berkeley.edu) **
44 ** **
45 ** ****************************************************************** */
46 
47 // $Revision: 1.1 $
48 // $Date: 2005/08/08 21:54:39 $
49 // $Source: /usr/local/cvs/OpenSees/SRC/material/uniaxial/Steel03.h,v $
50 
51 
52 #ifndef Steel03_h
53 #define Steel03_h
54 
55 // File: Steel03.h
56 //
57 // Written: mackie
58 // Created: 06/2005
59 // Revision: A
60 //
61 // Description: This file contains the class definition for
62 // Steel03.h. Steel03 is Steel01 verbatim but with added Giuffre-Menegotto-Pinto
63 // transitions on the loading and unloading loops.
64 // references:
65 // 1.) Menegotto, M., and Pinto, P.E. (1973). Method of analysis of cyclically loaded
66 // RC plane frames including changes in geometry and non-elastic behavior of
67 // elements under normal force and bending. Preliminary Report IABSE, vol 13.
68 // 2.) Dhakal, R.J., and Maekawa, K. (2002). Path-dependent cyclic stress-strain relationship
69 // of reinforcing bar including buckling. Engineering Structures, 24(11): 1383-96.
70 // 3.) Gomes, A., and Appleton, J. (1997). Nonlinear cyclic stress-strain relationship of
71 // reinforcing bars including buckling. Engineering Structures, 19(10): 822-6.
72 
73 //
74 //
75 //
76 // What: "@(#) Steel03.h, revA"
77 
78 
79 #include "material/uniaxial/steel/SteelBase0103.h"
80 
81 namespace XC {
83 //
85 class Steel03 : public SteelBase0103
86  {
87  private:
88  /*** Material Properties ***/
89  double r; // radius of rounded corners
90  double cR1;
91  double cR2;
92 
93  double CbStrain;
94  double CbStress;
95  double CrStrain;
96  double CrStress;
97  double Cplastic;
98 
99  /*** CONVERGED State Variables ***/
100  double CcurR;
101 
102  /*** TRIAL State Variables ***/
103  double TcurR;
104 
105  double TbStrain;
106  double TbStress;
107  double TrStrain;
108  double TrStress;
109  double Tplastic;
110 
111 
112  void determineTrialState(double dStrain);
113  double getR(double x_in);
114  protected:
115  int sendData(CommParameters &);
116  int recvData(const CommParameters &);
117  int setup_parameters(void);
118  public:
119  Steel03(int tag, double fy, double E0, double b, double r, double cR1, double cR2,
120  double a1 = STEEL_0103_DEFAULT_A1, double a2 = STEEL_0103_DEFAULT_A2,
121  double a3 = STEEL_0103_DEFAULT_A3, double a4 = STEEL_0103_DEFAULT_A4);
122  Steel03(int tag);
123  Steel03(void);
124 
125  int setTrialStrain(double strain, double strainRate = 0.0);
126  int setTrial(double strain, double &stress, double &tangent, double strainRate = 0.0);
127 
128  int commitState(void);
129  int revertToLastCommit(void);
130 
131  UniaxialMaterial *getCopy(void) const;
132 
133  int sendSelf(CommParameters &);
134  int recvSelf(const CommParameters &);
135 
136  void Print(std::ostream &s, int flag =0);
137  };
138 } // end of XC namespace
139 
140 #endif
int setup_parameters(void)
Sets all history and state variables to initial values.
Definition: Steel03.cpp:79
Steel 03 uniaxial material.
Definition: Steel03.h:85
double E0
Initial stiffness.
Definition: SteelBase.h:42
double a3
coefficient for isotropic hardening in tension
Definition: SteelBase.h:46
int sendSelf(CommParameters &)
Sends object through the channel being passed as parameter.
Definition: Steel03.cpp:298
int sendData(CommParameters &)
Send object members through the channel being passed as parameter.
Definition: Steel03.cpp:278
Base class for Steel01 and Steel03.
Definition: SteelBase0103.h:45
int revertToLastCommit(void)
Reset material to last committed state.
Definition: Steel03.cpp:260
double a4
coefficient for isotropic hardening in tension
Definition: SteelBase.h:47
int recvSelf(const CommParameters &)
Receives object through the channel being passed as parameter.
Definition: Steel03.cpp:312
double b
Hardening ratio (b = Esh/E0)
Definition: SteelBase.h:43
UniaxialMaterial * getCopy(void) const
Virtual constructor.
Definition: Steel03.cpp:274
double a1
coefficient for isotropic hardening in compression
Definition: SteelBase.h:44
double a2
coefficient for isotropic hardening in compression
Definition: SteelBase.h:45
Communication parameters between processes.
Definition: CommParameters.h:65
================================================================================
Definition: ContinuaReprComponent.h:34
void Print(std::ostream &s, int flag=0)
Imprime el objeto.
Definition: Steel03.cpp:330
double fy
Yield stress.
Definition: SteelBase.h:41
Base class for uniaxial materials.
Definition: UniaxialMaterial.h:88
int recvData(const CommParameters &)
Receives object members through the channel being passed as parameter.
Definition: Steel03.cpp:288