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 PINCHING_H
00041 #define PINCHING_H
00042
00043 #include <UniaxialMaterial.h>
00044 #include <Vector.h>
00045 #include <stdio.h>
00046
00047 class Pinching : public UniaxialMaterial
00048 {
00049 public:
00050 Pinching();
00051 Pinching(int tag, Vector inputParam );
00052 virtual ~Pinching();
00053
00054 const char *getClassType(void) const {return "Pinching";};
00055 int setTrialStrain(double d, double strainRate = 0.0);
00056 double getStrain(void);
00057
00058 double getStress(void);
00059 double getTangent(void);
00060 double getInitialTangent(void);
00061 int commitState(void);
00062 int revertToLastCommit(void);
00063 int revertToStart(void);
00064
00065
00066 UniaxialMaterial *getCopy(void);
00067
00068 int sendSelf(int commitTag, Channel &theChannel);
00069 int recvSelf(int commitTag, Channel &theChannel,
00070 FEM_ObjectBroker &theBroker);
00071
00072 void Print(OPS_Stream &s, int flag =0);
00073
00074 protected:
00075 void envelPosCap(double fy, double alfaPos, double alfaCap,
00076 double cpDsp, double d, double *f, double *ek );
00077
00078 void envelNegCap(double fy, double alfaNeg, double alfaCap,
00079 double cpDsp, double d, double *f, double *ek);
00080
00081 void recordInfo(int cond =0);
00082
00083
00084 private:
00085
00086
00087 double elstk,fyieldPos,fyieldNeg,alpha,Resfac;
00088 double capSlope,capDispPos,capDispNeg;
00089 double ecaps,ecapk,ecapa,ecapd,cs,ck,ca,cd;
00090 double fpPos,fpNeg,a_pinch;
00091
00092
00093 double dyieldPos,dyieldNeg,Enrgts,Enrgta,Enrgtk,Enrgtd;
00094
00095 double hsTrial[22], hsCommit[22], hsLastCommit[22];
00096
00097 };
00098
00099 #endif