Dear Silvia,
I have one question about truss element axialForce reliability analysis.
Here is the file:
model basic -ndm 2 -ndf 2
reliability
node 1 0.0 0.0
node 2 10.0 1.0
node 3 20.0 0.0
set E 200.0
set P 9.0
set A0 250.0
set fy 400.0
#uniaxialMaterial Elastic 1 $E
uniaxialMaterial Steel01 1 $fy $E 0.001
element truss 1 1 2 $A0 1
element truss 2 2 3 $A0 1
fix 1 1 1
fix 2 0 0
fix 3 1 1
pattern Plain 1 Linear {
load 2 0.0 $P
}
recorder Element -file ele1local.out -time -ele 1 axialForce
constraints Plain
numberer RCM
test NormUnbalance 1.0e-6 25 0
integrator LoadControl 1 1 1 1
algorithm Newton
system ProfileSPD
sensitivityIntegrator -static
sensitivityAlgorithm -computeAtEachStep
analysis Static
randomVariablePositioner 1 -createRV3 normal $E [expr 0.05*$E] -element 1 -material E
randomVariablePositioner 2 -createRV3 normal $E [expr 0.05*$E] -element 2 -material E
randomVariablePositioner 3 -createRV3 normal $fy [expr 0.05*$fy] -element 1 -material sigmaY
randomVariablePositioner 4 -createRV3 normal $fy [expr 0.05*$fy] -element 2 -material sigmaY
randomVariablePositioner 5 -createRV3 normal $P [expr 0.5*$P] -loadPattern 1 -loadAtNode 2 2
correlateGroup 1 2 0.9999
correlateGroup 3 4 0.9999
performanceFunction 1 "50.0-{rec_element_1_axialForce}"
randomNumberGenerator CStdLib
probabilityTransformation Nataf -print 0
reliabilityConvergenceCheck Standard -e1 1.0e-3 -e2 1.0e-3 -print 1
#gFunEvaluator OpenSees -file try_2_aux.tcl
gFunEvaluator OpenSees -analyze 1
#gradGEvaluator FiniteDifference
gradGEvaluator OpenSees
searchDirection iHLRF
meritFunctionCheck AdkZhang -multi 2.0 -add 10.0 -factor 0.5
stepSizeRule Armijo -maxNum 50 -base 0.5 -initial 1.0 2 -print 1 -sphere 1000.0 1.0 1.0
startPoint Mean
findDesignPoint StepSearch -maxNumIter 100 -printDesignPointX designPointX.out
runFORMAnalysis FORMoutput.out
Truss element axialForce reliability question
Moderators: silvia, selimgunay, Moderators