00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00035
00036
00037
00039
00040 #ifndef CloughDamage_H
00041 #define CloughDamage_H
00042
00043 #include <UniaxialMaterial.h>
00044 #include <Vector.h>
00045 #include <stdio.h>
00046
00047 class DamageModel;
00048
00049 class CloughDamage : public UniaxialMaterial
00050 {
00051 public:
00052 CloughDamage();
00053 CloughDamage(int tag, Vector inputParam, DamageModel *strength, DamageModel *stiffness,DamageModel *accelerated, DamageModel *capping);
00054 virtual ~CloughDamage();
00055
00056 const char *getClassType(void) const {return "CloughDamage";};
00057
00058 int setTrialStrain(double d, double strainRate = 0.0);
00059 double getStrain(void);
00060
00061 double getStress(void);
00062 double getTangent(void);
00063 double getInitialTangent(void);
00064 int commitState(void);
00065 int revertToLastCommit(void);
00066 int revertToStart(void);
00067
00068
00069 UniaxialMaterial *getCopy(void);
00070
00071 int sendSelf(int commitTag, Channel &theChannel);
00072 int recvSelf(int commitTag, Channel &theChannel,FEM_ObjectBroker &theBroker);
00073
00074 void Print(OPS_Stream &s, int flag =0);
00075
00076 protected:
00077 void envelPosCap(double fy, double alphaPos, double alphaCap,
00078 double cpDsp, double d, double *f, double *ek );
00079
00080 void envelNegCap(double fy, double alphaNeg, double alphaCap,
00081 double cpDsp, double d, double *f, double *ek);
00082
00083 void recordInfo(int cond =0);
00084
00085
00086 private:
00087
00088 double elstk,fyieldPos,fyieldNeg,alpha,Resfac;
00089 double capSlope,capDispPos,capDispNeg;
00090
00091
00092 double dyieldPos,dyieldNeg;
00093 DamageModel *StrDamage;
00094 DamageModel *StfDamage;
00095 DamageModel *AccDamage;
00096 DamageModel *CapDamage;
00097
00098 double hsTrial[24], hsCommit[24], hsLastCommit[24];
00099
00100 };
00101
00102 #endif