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
00032
00033 #ifndef RMC_PF_H
00034 #define RMC_PF_H
00035
00036 #include "PlasticFlow.h"
00037 #include <math.h>
00038
00039 class RMC_PF : public PlasticFlow
00040 {
00041 public:
00042 RMC_PF(int dilatant_which_in = -1, int index_dilatant_in = 0,
00043 int r_which_in = -1, int index_r_in = 0);
00044
00045 ~RMC_PF();
00046
00047 PlasticFlow* newObj();
00048
00049 const straintensor& PlasticFlowTensor(const stresstensor &Stre,
00050 const straintensor &Stra,
00051 const MaterialParameter &MaterialParameter_in) const;
00052
00053 private:
00054
00055 double getdilatant(const MaterialParameter &MaterialParameter_in) const;
00056 double getr(const MaterialParameter &MaterialParameter_in) const;
00057
00058 double RoundedFunctionf1(double s, double r) const;
00059 double RoundedFunctionf2(double s, double r) const;
00060 double RoundedFunctiondf1(double s, double r) const;
00061 double RoundedFunctiondf2(double s, double r) const;
00062
00063 private:
00064
00065 int dilatant_which;
00066 int index_dilatant;
00067 int r_which;
00068 int index_r;
00069
00070 static straintensor RMCm;
00071 };
00072
00073
00074 #endif
00075