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_YF_H
00034 #define RMC_YF_H
00035
00036 #include "YieldFunction.h"
00037 #include <math.h>
00038
00039 class RMC_YF : public YieldFunction
00040 {
00041 public:
00042 RMC_YF(int a_which_in = -1, int index_a_in = 0,
00043 int k_which_in = -1, int index_k_in = 0,
00044 int r_which_in = -1, int index_r_in = 0);
00045 ~RMC_YF();
00046
00047 YieldFunction *newObj();
00048
00049 double YieldFunctionValue(const stresstensor &Stre,
00050 const MaterialParameter &MaterialParameter_in) const;
00051
00052 const stresstensor& StressDerivative(const stresstensor &Stre,
00053 const MaterialParameter &MaterialParameter_in) const;
00054
00055 double InScalarDerivative(const stresstensor &Stre,
00056 const MaterialParameter &MaterialParameter_in,
00057 int which) const;
00058
00059 int getNumInternalScalar() const;
00060 int getNumInternalTensor() const;
00061 int getYieldFunctionRank() const;
00062
00063 private:
00064 double geta(const MaterialParameter &MaterialParameter_in) const;
00065 double getk(const MaterialParameter &MaterialParameter_in) const;
00066 double getr(const MaterialParameter &MaterialParameter_in) const;
00067
00068 double RoundedFunctionf1(double s, double r) const;
00069 double RoundedFunctionf2(double s, double r) const;
00070 double RoundedFunctiondf1(double s, double r) const;
00071 double RoundedFunctiondf2(double s, double r) const;
00072 double RoundedFunction(double s, double r) const;
00073
00074 private:
00075 int a_which;
00076 int index_a;
00077 int k_which;
00078 int index_k;
00079 int r_which;
00080 int index_r;
00081
00082 static stresstensor RMCst;
00083 };
00084
00085
00086 #endif
00087