Steel01.hGo to the documentation of this file.00001 /* ****************************************************************** ** 00002 ** OpenSees - Open System for Earthquake Engineering Simulation ** 00003 ** Pacific Earthquake Engineering Research Center ** 00004 ** ** 00005 ** ** 00006 ** (C) Copyright 1999, The Regents of the University of California ** 00007 ** All Rights Reserved. ** 00008 ** ** 00009 ** Commercial use of this program without express permission of the ** 00010 ** University of California, Berkeley, is strictly prohibited. See ** 00011 ** file 'COPYRIGHT' in main directory for information on usage and ** 00012 ** redistribution, and for a DISCLAIMER OF ALL WARRANTIES. ** 00013 ** ** 00014 ** Developed by: ** 00015 ** Frank McKenna (fmckenna@ce.berkeley.edu) ** 00016 ** Gregory L. Fenves (fenves@ce.berkeley.edu) ** 00017 ** Filip C. Filippou (filippou@ce.berkeley.edu) ** 00018 ** ** 00019 ** ****************************************************************** */ 00020 00021 // $Revision: 1.13 $ 00022 // $Date: 2006/09/05 22:27:12 $ 00023 // $Source: /usr/local/cvs/OpenSees/SRC/material/uniaxial/Steel01.h,v $ 00024 00025 00026 #ifndef Steel01_h 00027 #define Steel01_h 00028 00029 // Written: MHS 00030 // Created: 06/99 00031 // Revision: A 00032 // 00033 // Description: This file contains the class definition for 00034 // Steel01.h 00035 // 00036 // 00037 // 00038 // What: "@(#) Steel01.h, revA" 00039 00040 00041 #include <UniaxialMaterial.h> 00042 00043 // Default values for isotropic hardening parameters a1, a2, a3, and a4 00044 #define STEEL_01_DEFAULT_A1 0.0 00045 #define STEEL_01_DEFAULT_A2 55.0 00046 #define STEEL_01_DEFAULT_A3 0.0 00047 #define STEEL_01_DEFAULT_A4 55.0 00048 00049 class Steel01 : public UniaxialMaterial 00050 { 00051 public: 00052 Steel01(int tag, double fy, double E0, double b, 00053 double a1 = STEEL_01_DEFAULT_A1, double a2 = STEEL_01_DEFAULT_A2, 00054 double a3 = STEEL_01_DEFAULT_A3, double a4 = STEEL_01_DEFAULT_A4); 00055 Steel01(); 00056 ~Steel01(); 00057 00058 const char *getClassType(void) const {return "Steel01";}; 00059 00060 int setTrialStrain(double strain, double strainRate = 0.0); 00061 int setTrial (double strain, double &stress, double &tangent, double strainRate = 0.0); 00062 double getStrain(void); 00063 double getStress(void); 00064 double getTangent(void); 00065 double getInitialTangent(void) {return E0;}; 00066 00067 int commitState(void); 00068 int revertToLastCommit(void); 00069 int revertToStart(void); 00070 00071 UniaxialMaterial *getCopy(void); 00072 00073 int sendSelf(int commitTag, Channel &theChannel); 00074 int recvSelf(int commitTag, Channel &theChannel, 00075 FEM_ObjectBroker &theBroker); 00076 00077 void Print(OPS_Stream &s, int flag =0); 00078 00079 // AddingSensitivity:BEGIN ////////////////////////////////////////// 00080 int setParameter(const char **argv, int argc, Parameter ¶m); 00081 int updateParameter (int parameterID, Information &info); 00082 int activateParameter (int parameterID); 00083 double getStressSensitivity (int gradNumber, bool conditional); 00084 double getInitialTangentSensitivity(int gradNumber); 00085 int commitSensitivity (double strainGradient, int gradNumber, int numGrads); 00086 // AddingSensitivity:END /////////////////////////////////////////// 00087 00088 protected: 00089 00090 private: 00091 /*** Material Properties ***/ 00092 double fy; // Yield stress 00093 double E0; // Initial stiffness 00094 double b; // Hardening ratio (b = Esh/E0) 00095 double a1; 00096 double a2; 00097 double a3; 00098 double a4; // a1 through a4 are coefficients for isotropic hardening 00099 00100 /*** CONVERGED History Variables ***/ 00101 double CminStrain; // Minimum strain in compression 00102 double CmaxStrain; // Maximum strain in tension 00103 double CshiftP; // Shift in hysteresis loop for positive loading 00104 double CshiftN; // Shift in hysteresis loop for negative loading 00105 int Cloading; // Flag for loading/unloading 00106 // 1 = loading (positive strain increment) 00107 // -1 = unloading (negative strain increment) 00108 // 0 initially 00109 00110 /*** CONVERGED State Variables ***/ 00111 double Cstrain; 00112 double Cstress; 00113 double Ctangent; 00114 00115 /*** TRIAL History Variables ***/ 00116 double TminStrain; 00117 double TmaxStrain; 00118 double TshiftP; 00119 double TshiftN; 00120 int Tloading; 00121 00122 /*** TRIAL State Variables ***/ 00123 double Tstrain; 00124 double Tstress; 00125 double Ttangent; // Not really a state variable, but declared here 00126 // for convenience 00127 00128 // Calculates the trial state variables based on the trial strain 00129 void determineTrialState (double dStrain); 00130 00131 // Determines if a load reversal has occurred based on the trial strain 00132 void detectLoadReversal (double dStrain); 00133 00134 // AddingSensitivity:BEGIN ////////////////////////////////////////// 00135 int parameterID; 00136 Matrix *SHVs; 00137 // AddingSensitivity:END /////////////////////////////////////////// 00138 }; 00139 00140 #endif
Generated on Mon Oct 23 15:05:22 2006 for OpenSees by |
opensees-support @ berkeley.edu
Supported by the National Science Foundation
©2006, UC Regents