PEER NEES NEESit

• HOME
• USER
• DEVELOPER
• PROJECTS
• SUPPORT
• SITE MAP

• Dev Doc
• API
• Source
• Download
• Builds
• Bug Reports
• Message Board

OpenSeesGradGEvaluator.cpp Go to the documentation of this file. /* ****************************************************************** ** ** OpenSees - Open System for Earthquake Engineering Simulation ** ** Pacific Earthquake Engineering Research Center ** ** ** ** ** ** (C) Copyright 2001, The Regents of the University of California ** ** All Rights Reserved. ** ** ** ** Commercial use of this program without express permission of the ** ** University of California, Berkeley, is strictly prohibited. See ** ** file 'COPYRIGHT' in main directory for information on usage and ** ** redistribution, and for a DISCLAIMER OF ALL WARRANTIES. ** ** ** ** Developed by: ** ** Frank McKenna (fmckenna@ce.berkeley.edu) ** ** Gregory L. Fenves (fenves@ce.berkeley.edu) ** ** Filip C. Filippou (filippou@ce.berkeley.edu) ** ** ** ** Reliability module developed by: ** ** Terje Haukaas (haukaas@ce.berkeley.edu) ** ** Armen Der Kiureghian (adk@ce.berkeley.edu) ** ** ** ** ****************************************************************** */ // $Revision: 1.2 $ // $Date: 2003/10/27 23:45:44 $ // $Source: /usr/local/cvs/OpenSees/SRC/reliability/analysis/sensitivity/OpenSeesGradGEvaluator.cpp,v $ // // Written by Terje Haukaas (haukaas@ce.berkeley.edu) // #include <OpenSeesGradGEvaluator.h> #include <Vector.h> #include <Matrix.h> #include <GradGEvaluator.h> #include <ReliabilityDomain.h> #include <LimitStateFunction.h> #include <RandomVariable.h> #include <tcl.h> #include <string.h> #include <fstream> #include <iomanip> #include <iostream> using std::ifstream; using std::ios; using std::setw; using std::setprecision; using std::setiosflags; OpenSeesGradGEvaluator::OpenSeesGradGEvaluator( Tcl_Interp *passedTclInterp, ReliabilityDomain *passedReliabilityDomain, bool PdoGradientCheck) :GradGEvaluator(passedReliabilityDomain, passedTclInterp) { theReliabilityDomain = passedReliabilityDomain; doGradientCheck = PdoGradientCheck; int nrv = passedReliabilityDomain->getNumberOfRandomVariables(); grad_g = new Vector(nrv); grad_g_matrix = 0; DgDdispl = 0; } OpenSeesGradGEvaluator::~OpenSeesGradGEvaluator() { if (grad_g != 0) delete grad_g; if (DgDdispl != 0) delete DgDdispl; if (grad_g_matrix != 0) delete grad_g_matrix; } Vector OpenSeesGradGEvaluator::getGradG() { return (*grad_g); } Matrix OpenSeesGradGEvaluator::getAllGradG() { if (grad_g_matrix==0) { Matrix dummy(1,1); return dummy; } else { return (*grad_g_matrix); } } int OpenSeesGradGEvaluator::computeGradG(double g, Vector passed_x) { // Zero out the previous result matrix if (DgDdispl != 0) { delete DgDdispl; DgDdispl = 0; } // Call base class method computeParameterDerivatives(g); // Initial declaractions double perturbationFactor = 0.001; // (is multiplied by stdv and added to others...) char tclAssignment[500]; char *dollarSign = "$"; char *underscore = "_"; char lsf_expression[500] = ""; char separators[5] = "}{"; int nrv = theReliabilityDomain->getNumberOfRandomVariables(); RandomVariable *theRandomVariable; char tempchar[100]=""; char newSeparators[5] = "_"; double g_perturbed; int i; double onedudx; Vector dudx(nrv); // Compute gradients if this is a path-INdependent analysis // (This command only has effect if it IS path-independent.) sprintf(tclAssignment,"computeGradients"); Tcl_Eval( theTclInterp, tclAssignment ); // Initialize gradient vector grad_g->Zero(); // "Download" limit-state function from reliability domain int lsf = theReliabilityDomain->getTagOfActiveLimitStateFunction(); LimitStateFunction *theLimitStateFunction = theReliabilityDomain->getLimitStateFunctionPtr(lsf); char *theExpression = theLimitStateFunction->getExpression(); char *lsf_copy = new char[500]; strcpy(lsf_copy,theExpression); // Tokenize the limit-state function and COMPUTE GRADIENTS char *tokenPtr = strtok( lsf_copy, separators); while ( tokenPtr != NULL ) { strcpy(tempchar,tokenPtr); if ( strncmp(tokenPtr, "x",1) == 0) { // Get random variable number int rvNum; sscanf(tempchar,"x_%i",&rvNum); // Perturb its value according to its standard deviation theRandomVariable = theReliabilityDomain->getRandomVariablePtr(rvNum); double stdv = theRandomVariable->getStdv(); sprintf(tclAssignment , "set x_%d %35.20f", rvNum, (passed_x(rvNum-1)+perturbationFactor*stdv) ); Tcl_Eval( theTclInterp, tclAssignment); // Evaluate limit-state function again char *theTokenizedExpression = theLimitStateFunction->getTokenizedExpression(); Tcl_ExprDouble( theTclInterp, theTokenizedExpression, &g_perturbed ); // Make assignment back to its original value sprintf(tclAssignment , "set x_%d %35.20f", rvNum, passed_x(rvNum-1) ); Tcl_Eval( theTclInterp, tclAssignment); // Add gradient contribution (*grad_g)(rvNum-1) += (g_perturbed-g)/(perturbationFactor*stdv); } // If a nodal velocity is detected else if ( strncmp(tokenPtr, "ud", 2) == 0) { // Get node number and dof number int nodeNumber, direction; sscanf(tempchar,"ud_%i_%i", &nodeNumber, &direction); // Keep the original value double originalValue; sprintf(tclAssignment,"$ud_%d_%d", nodeNumber, direction); Tcl_ExprDouble( theTclInterp, tclAssignment, &originalValue); // Set perturbed value in the Tcl workspace double newValue = originalValue*(1.0+perturbationFactor); sprintf(tclAssignment,"set ud_%d_%d %35.20f", nodeNumber, direction, newValue); Tcl_Eval( theTclInterp, tclAssignment); // Evaluate the limit-state function again char *theTokenizedExpression = theLimitStateFunction->getTokenizedExpression(); Tcl_ExprDouble( theTclInterp, theTokenizedExpression, &g_perturbed ); // Compute gradient double onedgdu = (g_perturbed-g)/(originalValue*perturbationFactor); // Make assignment back to its original value sprintf(tclAssignment,"set ud_%d_%d %35.20f", nodeNumber, direction, originalValue); Tcl_Eval( theTclInterp, tclAssignment); // Obtain DDM gradient vector for (int i=1; i<=nrv; i++) { sprintf(tclAssignment , "set sens [sensNodeVel %d %d %d ]",nodeNumber,direction,i); Tcl_Eval( theTclInterp, tclAssignment); sprintf(tclAssignment , "$sens "); Tcl_ExprDouble( theTclInterp, tclAssignment, &onedudx ); dudx( (i-1) ) = onedudx; } // Add gradient contribution (*grad_g) += onedgdu*dudx; } // If a nodal displacement is detected else if ( strncmp(tokenPtr, "u", 1) == 0) { // Get node number and dof number int nodeNumber, direction; sscanf(tempchar,"u_%i_%i", &nodeNumber, &direction); // Keep the original value double originalValue; sprintf(tclAssignment,"$u_%d_%d", nodeNumber, direction); Tcl_ExprDouble( theTclInterp, tclAssignment, &originalValue); // Set perturbed value in the Tcl workspace double newValue = originalValue*(1.0+perturbationFactor); sprintf(tclAssignment,"set u_%d_%d %35.20f", nodeNumber, direction, newValue); Tcl_Eval( theTclInterp, tclAssignment); // Evaluate the limit-state function again char *theTokenizedExpression = theLimitStateFunction->getTokenizedExpression(); Tcl_ExprDouble( theTclInterp, theTokenizedExpression, &g_perturbed ); // Compute gradient double onedgdu = (g_perturbed-g)/(originalValue*perturbationFactor); // Store the DgDdispl in a matrix if (DgDdispl == 0) { DgDdispl = new Matrix(1, 3); (*DgDdispl)(0,0) = (double)nodeNumber; (*DgDdispl)(0,1) = (double)direction; (*DgDdispl)(0,2) = onedgdu; } else { int oldSize = DgDdispl->noRows(); Matrix tempMatrix = *DgDdispl; delete DgDdispl; DgDdispl = new Matrix(oldSize+1, 3); for (i=0; i<oldSize; i++) { (*DgDdispl)(i,0) = tempMatrix(i,0); (*DgDdispl)(i,1) = tempMatrix(i,1); (*DgDdispl)(i,2) = tempMatrix(i,2); } (*DgDdispl)(oldSize,0) = (double)nodeNumber; (*DgDdispl)(oldSize,1) = (double)direction; (*DgDdispl)(oldSize,2) = onedgdu; } // Make assignment back to its original value sprintf(tclAssignment,"set u_%d_%d %35.20f", nodeNumber, direction, originalValue); Tcl_Eval( theTclInterp, tclAssignment); // Obtain DDM gradient vector for (int i=1; i<=nrv; i++) { sprintf(tclAssignment , "set sens [sensNodeDisp %d %d %d ]",nodeNumber,direction,i); Tcl_Eval( theTclInterp, tclAssignment); sprintf(tclAssignment , "$sens "); Tcl_ExprDouble( theTclInterp, tclAssignment, &onedudx ); dudx( (i-1) ) = onedudx; } // Add gradient contribution (*grad_g) += onedgdu*dudx; } else if ( strncmp(tokenPtr, "rec",3) == 0) { } tokenPtr = strtok( NULL, separators); // read next token and go up and check the while condition again } delete [] lsf_copy; if (doGradientCheck) { char myString[100]; ofstream outputFile( "DDMgradients.out", ios::out ); opserr << endln; for (int ddm=0; ddm<grad_g->Size(); ddm++) { opserr << "DDM("<< (ddm+1) << ") = " << (*grad_g)(ddm) << endln; sprintf(myString,"%20.16e ",(*grad_g)(ddm)); outputFile << myString << endln; } outputFile.close(); opserr << "PRESS Ctrl+C TO TERMINATE APPLICATION!" << endln; while(true) { } } return 0; } int OpenSeesGradGEvaluator::computeAllGradG(Vector gFunValues, Vector passed_x) { // Allocate result matrix Vector gradG(passed_x.Size()); if (grad_g_matrix == 0) { grad_g_matrix = new Matrix(passed_x.Size(), gFunValues.Size()); } else { grad_g_matrix->Zero(); } // Loop over performance functions for (int j=1; j<=gFunValues.Size(); j++) { // Set tag of active limit-state function theReliabilityDomain->setTagOfActiveLimitStateFunction(j); this->computeGradG(gFunValues(j-1),passed_x); gradG = this->getGradG(); for (int i=1; i<=passed_x.Size(); i++) { (*grad_g_matrix)(i-1,j-1) = gradG(i-1); } } return 0; } Matrix OpenSeesGradGEvaluator::getDgDdispl() { return (*DgDdispl); } Generated on Mon Oct 23 15:05:25 2006 for OpenSees by  1.5.0

---

opensees-support @ berkeley.edu

Supported by the National Science Foundation

©2006, UC Regents