I require example for Finite Element Reliability Analysis

[Shahrakih]

Hi friends
I require an input example.tcl for Finite Element Reliability Analysis for RC moment frames .

[Shahrakih]

Hi Freinds
I've executed, for test, the first example proposed by Terje Haukas in his doctor thesis:


reliability
randomVariable 1 lognormal 500.0 100.0 500.0
randomVariable 2 lognormal 2000.0 400.0 2000.0
randomVariable 3 uniform 5.0 0.5 5.0
correlate 1 2 0.3
correlate 1 3 0.2
correlate 2 3 0.2
set a "{x_2}/(1000.0*{x_3})"
set b "{x_1}/(200.0*{x_3})"
performanceFunction 1 "1.0-\$a-\$b*\$b"
probabilityTransformation Nataf -print 0
reliabilityConvergenceCheck Standard -e1 1.0e-3 -e2 1.0e-3 -print 1
gFunEvaluator Basic
gradGEvaluator FiniteDifference -pert 1000
searchDirection iHLRF
meritFunctionCheck AdkZhang -multi 2.0 -add 10.0 -factor 0.5
stepSizeRule Armijo -maxNum 50 -base 0.5 -print 0
startPoint Mean
findDesignPoint StepSearch -maxNumIter 100
runFORMAnalysis FORMoutput.out


but the reliability index = 0
why ???????

[asalnaz]

by asalnaz » Fri Nov 04, 2011 9:47 pm

hi
i do this work but it is not true in results yet!i have index reliability but i dont have eta or delta!
i change gFunEvaluator and the results was beter.
can anyone help me?thank u.


# CREATE THE RELIABILITY MODEL BUILDER. NO FINITE ELEMENT MODEL BUILDER IS NECESSARY IN THIS EXAMPLE
reliability
randomVariable 1 lognormal 500.0 100.0 500.0
randomVariable 2 lognormal 2000.0 400.0 2000.0
randomVariable 3 uniform 5.0 0.5 5.0
correlate 1 2 0.3
correlate 1 3 0.2
correlate 2 3 0.2
#set a "\$xrv(2) / (1000.0*\$xrv(3))"
#set b "\$xrv(1) / (200.0*\$xrv(3))"
performanceFunction 1 "1.0 - (\$xrv(2) / (1000.0*\$xrv(3))) - ((\$xrv(1) / (200.0*\$xrv(3)))*(\$xrv(1) / (200.0*\$xrv(3))))"
probabilityTransformation Nataf -print 0
reliabilityConvergenceCheck Standard -e1 1.0e-3 -e2 1.0e-3 -print 1
gFunEvaluator Tcl -file "set a 1"
#gFunEvaluator Basic
gradGEvaluator FiniteDifference -pert 1000
searchDirection iHLRF
meritFunctionCheck AdkZhang -multi 2.0 -add 10.0 -factor 0.5
stepSizeRule Armijo -maxNum 50 -base 0.5 -print 0
startPoint Mean
findDesignPoint StepSearch -maxNumIter 100
runFORMAnalysis FORMoutput.out

[Shahrakih]

Hi Dear Asalnaz
you must substitute the last line with :
runFORMAnalysis FORMoutput.out -relSens 1
coodluck

[asalnaz]

Hi Dear Shahrakin
i try your sugestion, but it is not similar to answer of Terje Haukas example!!!!!

[asalnaz]

i want example about fosm and mcs analysis. is there any example about fosm (first order second moment) and mcs(mont carlo simulation) analysis?
thank u.

[linguan118]

I think the wiki manu has many examples contributed by Gu Quan
http://opensees.berkeley.edu/wiki/index ... l_Commands:

[Shahrakih]

Dear linguan
I checked this link :
http://opensees.berkeley.edu/wiki/index ... l_Commands
but I didn't find any example , do you run Finite Element Reliability Analysis
( FERA ) in OPENSEES ?