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
00032 #include <stdlib.h>
00033 #include <FedeasConcr3Material.h>
00034
00035 FedeasConcr3Material::FedeasConcr3Material(int tag,
00036 double fc, double ec, double fu, double eu,
00037 double rat, double ft, double epst0,
00038 double ft0, double beta, double epstu):
00039
00040 FedeasMaterial(tag, MAT_TAG_FedeasConcrete3, 2, 10)
00041 {
00042 data[0] = fc;
00043 data[1] = ec;
00044 data[2] = fu;
00045 data[3] = eu;
00046 data[4] = rat;
00047 data[5] = ft;
00048 data[6] = epst0;
00049 data[7] = ft0;
00050 data[8] = beta;
00051 data[9] = epstu;
00052
00053 tangent = 2.0*data[0]/data[1];
00054 tangentP = tangent;
00055 }
00056
00057 FedeasConcr3Material::FedeasConcr3Material(int tag, const Vector &d):
00058
00059 FedeasMaterial(tag, MAT_TAG_FedeasConcrete3, 2, 10)
00060 {
00061 if (d.Size() != numData) {
00062 opserr << "FedeasConcr3Material::FedeasConcr3Material -- not enough input arguments\n";
00063 exit(-1);
00064 }
00065
00066 for (int i = 0; i < numData; i++)
00067 data[i] = d(i);
00068
00069 tangent = 2.0*data[0]/data[1];
00070 tangentP = tangent;
00071 }
00072
00073 FedeasConcr3Material::FedeasConcr3Material(void):
00074 FedeasMaterial(0, MAT_TAG_FedeasConcrete3, 2, 10)
00075 {
00076
00077 }
00078
00079 FedeasConcr3Material::~FedeasConcr3Material(void)
00080 {
00081
00082 }
00083
00084 UniaxialMaterial*
00085 FedeasConcr3Material::getCopy(void)
00086 {
00087 Vector d(data, numData);
00088
00089 FedeasConcr3Material *theCopy = new FedeasConcr3Material(this->getTag(), d);
00090
00091
00092 for (int i = 0; i < 2*numHstv; i++)
00093 theCopy->hstv[i] = hstv[i];
00094
00095 theCopy->epsilonP = epsilonP;
00096 theCopy->sigmaP = sigmaP;
00097 theCopy->tangentP = tangentP;
00098
00099 theCopy->epsilon = epsilonP;
00100 theCopy->sigma = sigmaP;
00101 theCopy->tangent = tangentP;
00102
00103 return theCopy;
00104 }
00105
00106 double
00107 FedeasConcr3Material::getInitialTangent(void)
00108 {
00109
00110 return 2.0*data[0]/data[1];
00111 }