//////////////////////////////////////////////////////////////////////
// the proceure invoked by the interpreter when limitCurve is invoked
//////////////////////////////////////////////////////////////////////
int
Tcl_AddLimitCurveCommand (ClientData clientData, Tcl_Interp *interp, int argc,
TCL_Char **argv, Domain *theDomain)
{
OPS_ResetInputNoBuilder(clientData, interp, 2, argc, argv, theDomain);
// Make sure there is a minimum number of arguments
if (argc < 8)
{
opserr << "WARNING insufficient number of limit curve arguments\n";
opserr << "Want: limitCurve type? tag? <specific curve args>" << endln;
return TCL_ERROR;
}
// Pointer to a limit curve that will be added to the model builder
LimitCurve *theCurve = 0;
///////////////////////
// Axial Limit Curve
///////////////////////
if (strcmp(argv[1],"Axial") == 0)
{
if (argc != 9 && argc != 12 && argc != 14 && argc != 15)
{
opserr << "WARNING invalid number of arguments\n";
opserr << "Want: limitCurve Axial tag? eleTag? Fsw? Kdeg? Fres? defType? forType?" << endln; //SDK
opserr << "<ndI? ndJ? dof? perpDirn? delta? eleRemove?>" << endln;
return TCL_ERROR;
}
int tag;
int eleTag;
double Fsw;//SDK
double Kdeg;
double Fres;
int defType, forType;
int ndI = 0;
int ndJ = 0;
int dof = 0;
int perpDirn = 0;
int eleRemove = 0;
double delta = 0.0;
if (Tcl_GetInt(interp, argv[2], &tag) != TCL_OK)
{
opserr << "WARNING invalid Axial LimitCurve tag" << endln;
return TCL_ERROR;
}
if (Tcl_GetInt(interp, argv[3], &eleTag) != TCL_OK)
{
opserr << "WARNING invalid element tag for associated beam-column element (eleTag)\n";
opserr << "LimitCurve Axial: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetDouble(interp,argv[4], &Fsw) != TCL_OK)
{
opserr << "WARNING invalid Fsw\n";
opserr << "LimitCurve Axial: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetDouble(interp,argv[5], &Kdeg) != TCL_OK)
{
opserr << "WARNING invalid degrading slope Kdeg\n";
opserr << "LimitCurve Axial: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetDouble(interp,argv[6], &Fres) != TCL_OK)
{
opserr << "WARNING invalid residual capacity Fres\n";
opserr << "LimitCurve Axial: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetInt(interp,argv[7], &defType) != TCL_OK)
{
opserr << "WARNING invalid deformation type defType\n";
opserr << "LimitCurve Axial: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetInt(interp,argv[8], &forType) != TCL_OK)
{
opserr << "WARNING invalid force type forType\n";
opserr << "LimitCurve Axial: " << tag << endln;
return TCL_ERROR;
}
if (defType == 2)
{
if (Tcl_GetInt(interp,argv[9], &ndI) != TCL_OK)
{
opserr << "WARNING invalid node I\n";
opserr << "LimitCurve Axial: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetInt(interp,argv[10], &ndJ) != TCL_OK)
{
opserr << "WARNING invalid node J\n";
opserr << "LimitCurve Axial: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetInt(interp,argv[11], &dof) != TCL_OK)
{
opserr << "WARNING invalid degree of freedom for drift\n";
opserr << "LimitCurve Axial: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetInt(interp,argv[12], &perpDirn) != TCL_OK)
{
opserr << "WARNING invalid direction for column length\n";
opserr << "LimitCurve Axial: " << tag << endln;
return TCL_ERROR;
}
}
if (argc >= 14)
{
if (Tcl_GetDouble(interp,argv[13], &delta) != TCL_OK)
{
opserr << "WARNING invalid shift in drift surface (delta)\n";
opserr << "LimitCurve Axial: " << tag << endln;
return TCL_ERROR;
}
}
if (argc >= 15)
{
if (Tcl_GetInt(interp,argv[14], &eleRemove) != TCL_OK)
{
opserr << "WARNING invalid element removal option\n";
opserr << "LimitCurve Axial: " << tag << endln;
return TCL_ERROR;
}
}
// Parsing was successful, allocate the limit curve
// Subtract one from dof and perpDirn for C indexing
theCurve = new AxialCurve(interp, tag, eleTag, theDomain, Fsw, //SDK
Kdeg, Fres, defType, forType, ndI, ndJ, dof-1, perpDirn-1,
delta, eleRemove);
}
//////////////////////////
// Three Point LimitCurve
//////////////////////////
else if (strcmp(argv[1],"RotationShearCurve") == 0) {
void *theRSC = OPS_RotationShearCurve();
if (theRSC != 0) {
theCurve = (LimitCurve *)theRSC;
} else
return TCL_ERROR;
}
else if (strcmp(argv[1],"ThreePoint") == 0)
{
if (argc < 14 || argc > 18)
{
opserr << "WARNING insufficient arguments\n";
opserr << "Want: limitCurve ThreePoint tag? eleTag? x1? y1? x2? y2? x3? y3?";
opserr << "Kdeg? Fres? defType? forType?" << endln;
opserr << "<ndI? ndJ? dof? perpDirn?>" << endln;
return TCL_ERROR;
}
int tag;
int eleTag;
double Kdeg;
double Fres;
int defType, forType;
double x1, y1;
double x2, y2;
double x3, y3;
int ndI = 0;
int ndJ = 0;
int dof = 0;
int perpDirn = 0;
if (Tcl_GetInt(interp, argv[2], &tag) != TCL_OK)
{
opserr << "WARNING invalid limitCurve ThreePoint tag" << endln;
return TCL_ERROR;
}
if (Tcl_GetInt(interp, argv[3], &eleTag) != TCL_OK)
{
opserr << "WARNING invalid element tag for associated beam-column element (eleTag)\n";
opserr << "LimitCurve ThreePoint: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetDouble(interp, argv[4], &x1) != TCL_OK)
{
opserr << "WARNING invalid x1\n";
opserr << "limitCurve ThreePoint: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetDouble(interp, argv[5], &y1) != TCL_OK)
{
opserr << "WARNING invalid y1\n";
opserr << "limitCurve ThreePoint: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetDouble(interp, argv[6], &x2) != TCL_OK)
{
opserr << "WARNING invalid x2\n";
opserr << "limitCurve ThreePoint: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetDouble(interp, argv[7], &y2) != TCL_OK)
{
opserr << "WARNING invalid y2\n";
opserr << "limitCurve ThreePoint: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetDouble(interp, argv[8], &x3) != TCL_OK)
{
opserr << "WARNING invalid x3\n";
opserr << "limitCurve ThreePoint: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetDouble(interp, argv[9], &y3) != TCL_OK)
{
opserr << "WARNING invalid y3\n";
opserr << "limitCurve ThreePoint: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetDouble(interp,argv[10], &Kdeg) != TCL_OK)
{
opserr << "WARNING invalid degrading slope Kdeg\n";
opserr << "LimitCurve ThreePoint: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetDouble(interp,argv[11], &Fres) != TCL_OK)
{
opserr << "WARNING invalid residual capacity Fres\n";
opserr << "LimitCurve ThreePoint: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetInt(interp,argv[12], &defType) != TCL_OK)
{
opserr << "WARNING invalid deformation type defType\n";
opserr << "LimitCurve ThreePoint: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetInt(interp,argv[13], &forType) != TCL_OK)
{
opserr << "WARNING invalid force type forType\n";
opserr << "LimitCurve ThreePoint: " << tag << endln;
return TCL_ERROR;
}
if (defType == 2)
{
if (Tcl_GetInt(interp,argv[14], &ndI) != TCL_OK)
{
opserr << "WARNING invalid node I\n";
opserr << "LimitCurve ThreePoint: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetInt(interp,argv[15], &ndJ) != TCL_OK)
{
opserr << "WARNING invalid node J\n";
opserr << "LimitCurve ThreePoint: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetInt(interp,argv[16], &dof) != TCL_OK)
{
opserr << "WARNING invalid degree of freedom for drift\n";
opserr << "LimitCurve ThreePoint: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetInt(interp,argv[17], &perpDirn) != TCL_OK)
{
opserr << "WARNING invalid direction for column length\n";
opserr << "LimitCurve ThreePoint: " << tag << endln;
return TCL_ERROR;
}
}
// Parsing was successful, allocate the material
// Subtract one from dof and perpDirn for C indexing
theCurve = new ThreePointCurve(tag, eleTag, theDomain,
x1, y1, x2, y2, x3, y3, Kdeg, Fres, defType, forType,
ndI, ndJ, dof-1, perpDirn-1);
}
/////////////////////////////
// Shear Limit Curve
/////////////////////////////
else if (strcmp(argv[1],"Shear") == 0)
{
if (argc < 14 || argc > 19)
{ //SDK
opserr << "WARNING insufficient arguments\n";
// printCommand(argc,argv); // Commented out by Terje
opserr << "Want: limitCurve Shear tag? eleTag? rho? fc? b? h? d? Fsw? "; //SDK
opserr << "Kdeg? Fres? defType? forType?" << endln;
opserr << "<ndI? ndJ? dof? perpDirn? delta?>" << endln;
return TCL_ERROR;
}
int tag;
int eleTag;
double Kdeg;
double Fres;
int defType, forType;
double rho;
double fc;
double b, h, d;
int ndI = 0;
int ndJ = 0;
int dof = 0;
int perpDirn = 0;
double Fsw = 0.0; //SDK
double delta =0.0;
if (Tcl_GetInt(interp, argv[2], &tag) != TCL_OK)
{
opserr << "WARNING invalid limitCurve Shear tag" << endln;
return TCL_ERROR;
}
if (Tcl_GetInt(interp, argv[3], &eleTag) != TCL_OK)
{
opserr << "WARNING invalid element tag for associated beam-column element (eleTag)\n";
opserr << "LimitCurve Shear: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetDouble(interp, argv[4], &rho) != TCL_OK)
{
opserr << "WARNING invalid trans reinf ratio\n";
opserr << "limitCurve Shear: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetDouble(interp, argv[5], &fc) != TCL_OK)
{
opserr << "WARNING invalid concrete strength\n";
opserr << "limitCurve Shear: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetDouble(interp, argv[6], &b) != TCL_OK)
{
opserr << "WARNING invalid b\n";
opserr << "limitCurve Shear: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetDouble(interp, argv[7], &h) != TCL_OK)
{
opserr << "WARNING invalid h\n";
opserr << "limitCurve Shear: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetDouble(interp, argv[8], &d) != TCL_OK)
{
opserr << "WARNING invalid d\n";
opserr << "limitCurve Shear: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetDouble(interp, argv[9], &Fsw) != TCL_OK)
{ //SDK
opserr << "WARNING invalid Fsw\n";
opserr << "limitCurve Shear: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetDouble(interp,argv[10], &Kdeg) != TCL_OK)
{
opserr << "WARNING invalid degrading slope Kdeg\n";
opserr << "LimitCurve Shear: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetDouble(interp,argv[11], &Fres) != TCL_OK)
{
opserr << "WARNING invalid residual capacity Fres\n";
opserr << "LimitCurve Shear: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetInt(interp,argv[12], &defType) != TCL_OK)
{
opserr << "WARNING invalid deformation type defType\n";
opserr << "LimitCurve Shear: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetInt(interp,argv[13], &forType) != TCL_OK)
{
opserr << "WARNING invalid force type forType\n";
opserr << "LimitCurve Shear: " << tag << endln;
return TCL_ERROR;
}
if (defType == 2)
{
if (Tcl_GetInt(interp,argv[14], &ndI) != TCL_OK)
{
opserr << "WARNING invalid node I\n";
opserr << "LimitCurve Shear: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetInt(interp,argv[15], &ndJ) != TCL_OK)
{
opserr << "WARNING invalid node J\n";
opserr << "LimitCurve Shear: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetInt(interp,argv[16], &dof) != TCL_OK)
{
opserr << "WARNING invalid degree of freedom for drift\n";
opserr << "LimitCurve Shear: " << tag << endln;
return TCL_ERROR;
}
if (Tcl_GetInt(interp,argv[17], &perpDirn) != TCL_OK)
{
opserr << "WARNING invalid direction for column length\n";
opserr << "LimitCurve Shear: " << tag << endln;
return TCL_ERROR;
}
}
if (argc == 19)
{
if (Tcl_GetDouble(interp,argv[18], &delta) != TCL_OK) {
opserr << "WARNING invalid shift in drift surface (delta)\n";
opserr << "LimitCurve Shear: " << tag << endln;
return TCL_ERROR;
}
}
// Parsing was successful, allocate the material
// Subtract one from dof and perpDirn for C indexing
theCurve = new ShearCurve(tag, eleTag, theDomain,
rho, fc, b, h, d, Fsw, Kdeg, Fres, defType, forType, //SDK
ndI, ndJ, dof-1, perpDirn-1, delta);
}
//***********************************************************************************************************************start package commands //**MRL
else
{
//MRL Start Comment to allow OpenSees to look in package
/*opserr << "WARNING unknown type of limitCurve: " << argv[1];
opserr << "\nValid types: Axial, ThreePoint, Shear" << endln;
return TCL_ERROR;
}*/ //MRL End Comment
//////////////////////////////////////
// Curve could be defined in a package
//////////////////////////////////////
if (theCurve == 0)
{
// maybe element in a class package already loaded
// loop through linked list of loaded functions comparing names & if find call it
LimitCurvePackageCommand *limCrvCommands = theLimitCurvePackageCommands;
bool found = false;
while (limCrvCommands != NULL && found == false)
{
if (strcmp(argv[1], limCrvCommands->funcName) == 0)
{
OPS_ResetInputNoBuilder(clientData, interp, 2, argc, argv, theDomain);
theCurve = (LimitCurve *)(*(limCrvCommands->funcPtr))(argc, argv);
found = true;;
}
else
limCrvCommands = limCrvCommands->next;
}
}
if (theCurve == 0)
{
// check to see if element is a procedure
// maybe curve in a routine
char *limCrvType = new char[strlen(argv[1])+1];
strcpy(limCrvType, argv[1]);
limCrvObj *limCrvObject = OPS_GetLimitCurveType(limCrvType, strlen(limCrvType));
delete [] limCrvType;
if (limCrvObject != 0)
{
theCurve = Tcl_addWrapperLimitCurve(limCrvObject, clientData, interp, argc, argv);
if (theCurve == 0)
delete limCrvObject;
}
}
// maybe limit curve class exists in a package yet to be loaded
if (theCurve == 0)
{
void *libHandle;
void * (*funcPtr)(int argc, TCL_Char **argv);
int limCrvNameLength = strlen(argv[1]);
char *tclFuncName = new char[limCrvNameLength+12];
strcpy(tclFuncName, "OPS_");
strcpy(&tclFuncName[4], argv[1]);
int res = getLibraryFunction(argv[1], tclFuncName, &libHandle, (void **)&funcPtr);
delete [] tclFuncName;
if (res == 0)
{
// add loaded function to list of functions
char *limCrvName = new char[limCrvNameLength+1];
strcpy(limCrvName, argv[1]);
LimitCurvePackageCommand *theLimCrvCommand = new LimitCurvePackageCommand;
theLimCrvCommand->funcPtr = funcPtr;
theLimCrvCommand->funcName = limCrvName;
theLimCrvCommand->next = theLimitCurvePackageCommands;
theLimitCurvePackageCommands = theLimCrvCommand;
OPS_ResetInputNoBuilder(clientData, interp, 2, argc, argv, theDomain);
theCurve = (LimitCurve *)(*funcPtr)(argc, argv);
}
}
// if still here the element command does not exist
if (theCurve == 0)
{
opserr << "WARNING could not create LimitCurve " << argv[1] << endln;
return TCL_ERROR;
}
// Now add the limit curve to the modelBuilder
if (OPS_addLimitCurve(theCurve) < 0)
{
opserr << "WARNING could not add LimitCurve to the domain\n";
opserr << *theCurve << endln;
delete theCurve; // invoke the limit curve objects destructor, otherwise mem leak
return TCL_ERROR;
}
//MRL end else for package
}
// Ensure we have created the Material, out of memory if got here and no material
if (theCurve == 0) {
opserr << "WARNING ran out of memory creating limitCurve\n";
opserr << argv[1] << endln;
return TCL_ERROR;
}
// Now add the curve to the modelBuilder
if (OPS_addLimitCurve(theCurve) < 0) {
opserr << "WARNING could not add limitCurve to the domain\n";
opserr << *theCurve << endln;
delete theCurve; // invoke the curve objects destructor, otherwise mem leak
return TCL_ERROR;
}
return TCL_OK;
}
int OPS_LimitCurve()
{
// Make sure there is a minimum number of arguments
if (OPS_GetNumRemainingInputArgs() < 6)
{
opserr << "WARNING insufficient number of limit curve arguments\n";
opserr << "Want: limitCurve type? tag? <specific curve args>" << endln;
return -1;
}
const char* type = OPS_GetString();
// Pointer to a limit curve that will be added to the model builder
LimitCurve *theCurve = 0;
if (strcmp(type, "Axial") == 0) {
opserr << "WARNING to be implemented ...\n";
return -1;
} else if (strcmp(type, "RotationShearCurve") == 0) {
void *theRSC = OPS_RotationShearCurve();
if (theRSC != 0) {
theCurve = (LimitCurve *)theRSC;
} else {
return -1;
}
} else if (strcmp(type,"ThreePoint") == 0) {
void* curve = OPS_RotationShearCurve();
if (curve != 0) {
theCurve = (LimitCurve *)curve;
} else {
return -1;
}
} else if (strcmp(type,"Shear") == 0) {
void* curve = OPS_ShearCurve();
if (curve != 0) {
theCurve = (LimitCurve *)curve;
} else {
return -1;
}
} else {
opserr << "WARNING type of limit curve is unknown\n";
return -1;
}
// Ensure we have created the Material, out of memory if got here and no material
if (theCurve == 0) {
opserr << "WARNING ran out of memory creating limitCurve\n";
return -1;
}
// Now add the material to the modelBuilder
if (OPS_addLimitCurve(theCurve) == false) {
opserr << "WARNING could not add limitCurve to the domain\n";
delete theCurve; // invoke the material objects destructor, otherwise mem leak
return -1;
}
return 0;
}