Tcl Procedure for Reliability

A place to ask questions on how to use OpenSees to perform Finite Element Reliability Analysis

Moderators: silvia, selimgunay, mhscott, Moderators

Post Reply
giorgioirotina
Posts: 8
Joined: Tue Jan 21, 2014 11:39 am

Tcl Procedure for Reliability

Post by giorgioirotina » Fri Mar 20, 2015 4:45 am

I have to perform reliability analysis of several different configurations of the same structure.
Since the problem is always the same I defined a procedure where the only inputs are R and S.
When I call the procedure I get the results once and the second time I get an error.
Here is the procedure

proc betaFORM {R S } {


pause
reliability

randomVariable 1 lognormal -mean $R -stdv [expr $R*.18]
randomVariable 2 gumbel -mean $S -stdv [expr $S*.1]
randomVariable 3 normal -mean 1.1 -stdv [expr 1.1*.06]
randomVariable 4 normal -mean 1. -stdv .05
randomVariable 5 normal -mean 1. -stdv .02
randomVariable 6 normal -mean 1. -stdv .1

parameter 1 randomVariable 1
parameter 2 randomVariable 2
parameter 3 randomVariable 3
parameter 4 randomVariable 4
parameter 5 randomVariable 5
parameter 6 randomVariable 6





performanceFunction 1 "\$par(1)-(\$par(2)*\$par(3)*\$par(4)*\$par(5)*\$par(6))"

functionEvaluator Tcl
gradientEvaluator FiniteDifference -pert 1000
startPoint Mean
runFOSMAnalysis barFOSM.out
probabilityTransformation Nataf
searchDirection iHLRF
reliabilityConvergenceCheck Standard -e1 1.0e-4 -e2 1.0e-4
stepSizeRule Fixed -stepSize 1.0
findDesignPoint StepSearch -maxNumIter 10

runFORMAnalysis barFORM.out

set beta $betaFORM(1)
puts "beta $beta"


return $beta

}

mhscott
Posts: 868
Joined: Tue Jul 06, 2004 3:38 pm
Location: Corvallis, Oregon USA
Contact:

Re: Tcl Procedure for Reliability

Post by mhscott » Tue Jul 21, 2015 12:06 pm

Add the 'wipeReliability' command ... your error is probably due to multiply defined rv tags.

Post Reply