Reliability

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

Moderators: silvia, selimgunay, mhscott, Moderators

Post Reply
shavilasht
Posts: 1
Joined: Sun Oct 13, 2013 8:48 pm
Location: mazandaran-iran

Reliability

Post by shavilasht »

hi
I need some reliability example to help me how can do correct FORM Analyze in opensees.
please let me some complete example.
thanks.
mhscott
Posts: 880
Joined: Tue Jul 06, 2004 3:38 pm
Location: Corvallis, Oregon USA
Contact:

Re: Reliability

Post by mhscott »

Here is a FORM pushover analysis of a steel frame. It's a little heavy on Tcl lists and arrays.

http://opensees.berkeley.edu/WebSVN/fil ... lFrame.tcl
http://opensees.berkeley.edu/WebSVN/fil ... Member.tcl
luiscelorrio
Posts: 8
Joined: Wed Aug 07, 2013 6:37 am
Location: Logroño Spain
Contact:

Re: Reliability

Post by luiscelorrio »

Hello, I have read the paper: "Shape sensitivities in the reliability analysis of nonlinear frame structures". In this paper two examples are explained. I thanks that Michael had uploaded the first example involving a steel structure. I have run it and it works very well. Can you please upload the second example of a reinforced concrete structure?

I have some questions about section command.
Is possible to use the command section WFSection2d to compute DDM based sensitivity?

I have found that the information about the command Section Aggregator in the user manual and in the Sensitivity Command Manual is different. What is the correct command?

Thanks in advance.
Luis
mhscott
Posts: 880
Joined: Tue Jul 06, 2004 3:38 pm
Location: Corvallis, Oregon USA
Contact:

Re: Reliability

Post by mhscott »

luiscelorrio wrote:
> Hello, I have read the paper: "Shape sensitivities in the reliability analysis of
> nonlinear frame structures". In this paper two examples are explained. I thanks
> that Michael had uploaded the first example involving a steel structure. I have
> run it and it works very well. Can you please upload the second example of a reinforced
> concrete structure?

I'm working on cleaning up the RC frame example from that paper.

> I have some questions about section command.
> Is possible to use the command section WFSection2d to compute DDM based sensitivity?

Yes

> I have found that the information about the command Section Aggregator in the user
> manual and in the Sensitivity Command Manual is different. What is the correct
> command?

I'm not sure what the differences are. I would go with what's in the user manual. Third parties wrote the sensitivity command manual, so it could be out of date.


> Thanks in advance.
> Luis
mgrubisic
Posts: 9
Joined: Tue Jan 31, 2012 6:08 am
Location: University of Osijek, Croatia
Contact:

Re: Reliability

Post by mgrubisic »

Dear Michael Scott,

can you share with us the example of FORM analysis of 2-storey 2-bay RC frame from your paper "Shape sensitivities in the reliability analysis of nonlinear frame structures"? It would be a great help!

Thanks in advance!
Marin
Marin Grubišić
Assist. Professor, Structural & Earthquake Eng.
University of Osijek, Faculty of Civil Engineering and Architecture Osijek
3 Vladimir Prelog Street, HR-31000 Osijek, Croatia
marin.grubisic@gfos.hr | www.github.com/mgrubisic
AnggaSetiawan
Posts: 15
Joined: Tue Apr 07, 2015 10:17 pm
Location: Kyoto University

Re: Reliability

Post by AnggaSetiawan »

Dear Michael Scott

Would you like to give me example tcl file of the structure optimization with SNOPT analysis?

Thank you
Technip
Posts: 1
Joined: Tue Mar 14, 2017 7:50 am

Re: Reliability - GFunEvaluator Error

Post by Technip »

Hi Guys,
I am looking for help on the code I have which include the Reliability Assessment.
the problem I have is invalid command GFunEvaluator when run the tcl code using source command. Sorry I could not find the way to attach the tcl file so I copy paste below. Any help on this appreciated.
Regards
Amer Ali
____________________________________________________________

# CREATE THE RELIABILITY MODEL BUILDER. NO FINITE ELEMENT MODEL BUILDER IS NECESSARY IN THIS EXAMPLE
reliability

#
# --- flowline
# (mm)
set fl_od_nom 311.0

#
set fl_oval_max 0.0075
set fl_oval_mean 0.00375
set fl_oval_sd 0.00125
#
set fl_oval_A_mean [expr $fl_od_nom * $fl_oval_mean]
set fl_oval_A_sd [expr $fl_od_nom * $fl_oval_sd]
#
# --- carrier
# (mm)
set cp_od_nom 434.5
set cp_wt_mean 21.75
set cp_wt_sd 0.725
#
set cp_oval_max 0.01
set cp_oval_mean 0.005
set cp_oval_sd 0.001667
#
set cp_oval_A_mean [expr $cp_od_nom * $cp_oval_mean]
set cp_oval_A_sd [expr $cp_od_nom * $cp_oval_sd]
#
# --- centraliser
#
set ct_wt_mean 34.25
set ct_wt_sd 0.0835
#
# --- weld bead
#
set bd_height_max 2.0
set bd_height_mean 0.667
set bd_height_sd 0.333
#
# --- misalignment
#
set ms_max 1.5
set ms_mean 0.5
set ms_sd 0.25
#
# CREATE RANDOM VARIABLES
# -tag- -type- -mean- -sd-
randomVariable 1 lognormal $cp_wt_mean $cp_wt_sd
randomVariable 2 lognormal $cp_wt_mean $cp_wt_sd
randomVariable 3 normal $fl_oval_A_mean $fl_oval_A_sd
randomVariable 4 normal $cp_oval_A_mean $cp_oval_A_sd
randomVariable 5 lognormal $ct_wt_mean $ct_wt_sd
randomVariable 6 lognormal $ct_wt_mean $ct_wt_sd
randomVariable 7 lognormal $bd_height_mean $bd_height_sd
randomVariable 8 lognormal $bd_height_mean $bd_height_sd
randomVariable 9 lognormal $ms_mean $ms_sd


## ............................ SPECIFY CORRELATION
correlate 5 6 0.5

set gap "$cp_od_nom - $fl_od_nom -{x_1} - {x_2} - {x_3} - {x_4} - {x_5} - {x_6} - {x_7} -{x_8} - {x_9}"
#
# -"expression"-
#
performanceFunction 1 "$gap"
#
# CREATE NECESSARY RELIABILITY ANALYSIS TOOLS
#
probabilityTransformation Nataf -print 1
reliabilityConvergenceCheck Standard -e1 1.0e-3 -e2 1.1e-2 -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 -initial 0.01 2 -print 1
startPoint Mean
findDesignPoint StepSearch -maxNumIter 100

# =============================
# === RUN THE FORM ANALYSIS ===
# =============================

runFORMAnalysis FORM34t.out
Last edited by Technip on Tue Mar 14, 2017 9:58 am, edited 1 time in total.
RDB
Posts: 3
Joined: Wed Mar 22, 2017 3:33 am
Location: TechnipFMC

Re: Reliability

Post by RDB »

Hi,

I'm hoping someone can help me. I am working with some code which was written 10years or so ago which uses the commands

gFunEvaluator Basic
gradGEvaluator FiniteDifference -pert 1000

Having searched the internet it would appear that both of these commands have been updated to distinguish that they are tcl commands. From what I have found, similar commands are now written as:

functionEvaluator Tcl -file "Analyze -1"
gradientEvaluator Implicit (although i've replaced this with "FiniteDifference -pert 1000")

Could someone confirm if this is correct or if not, what are the updated versions of these commands.

Thanks,

Ross
Rachellara
Posts: 1
Joined: Tue Jun 11, 2019 10:29 pm
Location: United States

Re: Reliability

Post by Rachellara »

mhscott wrote:
> Here is a FORM pushover analysis of a steel frame. It's a little heavy on
> Tcl lists and arrays.
>
>
> http://opensees.berkeley.edu/WebSVN/fil ... lFrame.tcl
>
> http://opensees.berkeley.edu/WebSVN/fil ... Member.tcl


Hello,
Thanks for this links i have searching for the steel frame tcl list from the last 2 hours. I hope now i can resolve my problem.
leeee
Posts: 2
Joined: Sun May 31, 2020 6:27 pm

Re: Reliability

Post by leeee »

Dear professor scoot. Can you post the link again, the previous link is no longer valid.I needed this case so badly.Thanks
mhscott wrote: Mon Jan 13, 2014 9:14 am Here is a FORM pushover analysis of a steel frame. It's a little heavy on Tcl lists and arrays.

http://opensees.berkeley.edu/WebSVN/fil ... lFrame.tcl
http://opensees.berkeley.edu/WebSVN/fil ... Member.tcl
mhscott
Posts: 880
Joined: Tue Jul 06, 2004 3:38 pm
Location: Corvallis, Oregon USA
Contact:

Re: Reliability

Post by mhscott »

Sorry, just saw your message. I hope this is what you're looking for:

https://github.com/OpenSees/OpenSees/tr ... eliability
imarketingmy
Posts: 1
Joined: Fri May 06, 2022 7:18 am
Contact:

Re: Reliability

Post by imarketingmy »

good, very helpfull
Post Reply