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 #include <DrainClough2Material.h>
00032 #include <Vector.h>
00033
00034 DrainClough2Material::DrainClough2Material(int tag,
00035 double E, double fyp, double fyn, double alpha,
00036 double ecaps, double ecapk, double ecapa, double ecapd,
00037 double cs, double ck, double ca, double cd,
00038 double capSlope, double capDispP, double capDispN, double res,
00039 double b):
00040
00041 DrainMaterial(tag, MAT_TAG_DrainClough2, 19, 16, b)
00042 {
00043 data[0] = E;
00044 data[1] = fyp;
00045 data[2] = fyn;
00046 data[3] = alpha;
00047 data[4] = ecaps;
00048 data[5] = ecapk;
00049 data[6] = ecapa;
00050 data[7] = ecapd;
00051 data[8] = cs;
00052 data[9] = ck;
00053 data[10] = ca;
00054 data[11] = cd;
00055 data[12] = capSlope;
00056 data[13] = capDispP;
00057 data[14] = capDispN;
00058 data[15] = res;
00059
00060
00061 this->revertToStart();
00062 }
00063
00064 DrainClough2Material::DrainClough2Material(int tag, const Vector &input, double b):
00065
00066 DrainMaterial(tag, MAT_TAG_DrainClough2, 19, 16, b)
00067 {
00068 for (int i = 0; i < 16; i++)
00069 data[i] = input(i);
00070
00071
00072 this->revertToStart();
00073 }
00074
00075 DrainClough2Material::DrainClough2Material(void):
00076 DrainMaterial(0, MAT_TAG_DrainClough2, 19, 16)
00077 {
00078
00079 }
00080
00081 DrainClough2Material::~DrainClough2Material(void)
00082 {
00083
00084 }
00085
00086 int
00087 DrainClough2Material::revertToStart(void)
00088 {
00089 double dyp = data[1]/data[0];
00090 double dyn = data[2]/data[0];
00091
00092 hstv[0] = data[0];
00093 hstv[1] = data[0];
00094 hstv[2] = dyp;
00095 hstv[3] = dyn;
00096 hstv[4] = 0.0;
00097 hstv[5] = dyp;
00098 hstv[6] = dyn;
00099 hstv[7] = data[1];
00100 hstv[8] = data[2];
00101 hstv[9] = data[13];
00102 hstv[10] = data[14];
00103 hstv[11] = 0.0;
00104 hstv[12] = 0.0;
00105 hstv[13] = 0.0;
00106 hstv[14] = data[0];
00107 hstv[15] = 0.0;
00108 hstv[16] = 0.0;
00109 hstv[17] = 0.0;
00110 hstv[18] = 0.0;
00111
00112
00113 for (int i = 0; i < 19; i++)
00114 hstv[i+19] = hstv[i];
00115
00116 return 0;
00117 }
00118
00119 UniaxialMaterial*
00120 DrainClough2Material::getCopy(void)
00121 {
00122 Vector input(data, 16);
00123
00124 DrainClough2Material *theCopy =
00125 new DrainClough2Material(this->getTag(), input, beto);
00126
00127 return theCopy;
00128 }