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 DP_YF_H
00034 #define DP_YF_H
00035
00036 #include "YieldFunction.h"
00037 #include <math.h>
00038
00039 class DP_YF : public YieldFunction
00040 {
00041 public:
00042 DP_YF(int a_which_in = -1, int index_a_in = 0,
00043 int k_which_in = -1, int index_k_in = 0,
00044 int alpha_which_in = -1, int index_alpha_in = 0);
00045 ~DP_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 const stresstensor& InTensorDerivative(const stresstensor &Stre,
00060 const MaterialParameter &MaterialParameter_in,
00061 int which) const;
00062
00063 int getNumInternalScalar() const;
00064 int getNumInternalTensor() const;
00065 int getYieldFunctionRank() const;
00066
00067 private:
00068 double geta(const MaterialParameter &MaterialParameter_in) const;
00069 double getk(const MaterialParameter &MaterialParameter_in) const;
00070 const stresstensor& getalpha(const MaterialParameter &MaterialParameter_in) const;
00071
00072 private:
00073 int a_which;
00074 int index_a;
00075 int k_which;
00076 int index_k;
00077 int alpha_which;
00078 int index_alpha;
00079
00080 static stresstensor DPst;
00081 };
00082
00083
00084 #endif
00085