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