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 #ifndef FeapMaterial_h
00032 #define FeapMaterial_h
00033
00034 #include <NDMaterial.h>
00035
00036 class FeapMaterial : public NDMaterial
00037 {
00038 public:
00039 FeapMaterial(int tag, int classTag, int numHV, int numData,
00040 double rho = 0.0);
00041 FeapMaterial(int classTag);
00042 virtual ~FeapMaterial();
00043
00044 virtual const char *getClassType(void) const {return "FeapMaterial";};
00045
00046 virtual int setTrialStrain(const Vector &strain);
00047 virtual const Vector &getStrain(void);
00048 virtual const Vector &getStress(void);
00049 virtual const Matrix &getTangent(void);
00050 virtual double getRho(void);
00051
00052 virtual int commitState(void);
00053 virtual int revertToLastCommit(void);
00054 virtual int revertToStart(void);
00055
00056 virtual NDMaterial *getCopy(void);
00057 virtual NDMaterial *getCopy(const char *type);
00058 virtual const char *getType(void) const;
00059 virtual int getOrder(void) const;
00060
00061 virtual int sendSelf(int commitTag, Channel &theChannel);
00062 virtual int recvSelf(int commitTag, Channel &theChannel,
00063 FEM_ObjectBroker &theBroker);
00064
00065 virtual void Print(OPS_Stream &s, int flag = 0);
00066
00067 protected:
00068
00069 virtual int invokeSubroutine(int isw);
00070 virtual int fillDArray(void);
00071
00072 double *ud;
00073 double *hstv;
00074
00075 static double d[];
00076
00077 double rho;
00078
00079 private:
00080 int numHV;
00081 int numData;
00082
00083 double eps[6];
00084 static double sig[6];
00085 static double dd[36];
00086
00087 static Vector strain3;
00088 static Vector strain4;
00089 static Vector strain6;
00090
00091 static Vector sigma3;
00092 static Vector sigma4;
00093 static Vector sigma6;
00094
00095 static Matrix tangent3;
00096 static Matrix tangent4;
00097 static Matrix tangent6;
00098
00099 enum Formulation{Unknown, ThreeDimensional, PlaneStrain, AxiSymmetric};
00100 int myFormulation;
00101 void setType(Formulation form);
00102 };
00103
00104 #endif