[ot-users] Open Turns FORM

Douglas Long douglaslong11 at gmail.com
Fri Aug 4 19:48:23 CEST 2017


my apologies. there was a type in the distributions list.

import openturns as ot


myFunction = ot.NumericalMathFunction(['P', 'L', 'W', 'T'], ['d'],
                                      ['W*T - P*L/4'])

distributions_list = [
ot.Normal(0.0001, 0.00002),
ot.Normal(600000, 100000),
ot.Normal(10, 2),
ot.Normal(8, 0.0001)
]
copula = ot.IndependentCopula(4)
compose_distribution = ot.ComposedDistribution(distributions_list, copula)


vect = ot.RandomVector(compose_distribution)
output = ot.RandomVector(myFunction, vect)
myEvent = ot.Event(output, ot.Less(), 0)

myCobyla = ot.Cobyla()
myAlgo = ot.FORM(myCobyla, myEvent, [0.0001,600000,10,8])
myAlgo.run()
result = myAlgo.getResult()

print(result.getPhysicalSpaceDesignPoint())
print(result.getHasoferReliabilityIndex())
print(result.getEventProbability())
print(result.getLimitStateVariable())
print(result.getStandardSpaceDesignPoint())
print(result.getImportanceFactors())


On Fri, Aug 4, 2017 at 11:38 AM, Douglas Long <douglaslong11 at gmail.com>
wrote:

> Hey Folks,
>
> I am attempting to recreate a FORM example.
>
>    - http://www2.mae.ufl.edu/haftka/stropt/Lectures/FORM.pdf
>    - I am trying to get BETA = 2.48 as in the example linked above.
>
>
>
> here is my code but my results are different. I have tried many different
> solutions.
> any help would be greatly appreciated.
>
> Thanks,
> Doug
>
> import openturns as ot
>
>
> myFunction = ot.NumericalMathFunction(['P', 'L', 'W', 'T'], ['d'],
>                                       ['W*T - P*L/4'])
>
> distributions_list = [
> ot.Normal(0.0001, 2),
> ot.Normal(600000, 0.1),
> ot.Normal(10, 0.00002),
> ot.Normal(8, 100000)
> ]
> copula = ot.IndependentCopula(4)
> compose_distribution = ot.ComposedDistribution(distributions_list, copula)
>
>
> vect = ot.RandomVector(compose_distribution)
> output = ot.RandomVector(myFunction, vect)
> myEvent = ot.Event(output, ot.Less(), 0)
>
> myCobyla = ot.Cobyla()
> myAlgo = ot.FORM(myCobyla, myEvent, [0.0001,600000,10,8])
> myAlgo.run()
> result = myAlgo.getResult()
>
> print(result.getPhysicalSpaceDesignPoint())
> print(result.getHasoferReliabilityIndex())
> print(result.getEventProbability())
> print(result.getLimitStateVariable())
> print(result.getStandardSpaceDesignPoint())
> print(result.getImportanceFactors())
>
>
> --
> Douglas Long
> douglaslong11 at gmail.com
>



-- 
Douglas Long
douglaslong11 at gmail.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://openturns.org/pipermail/users/attachments/20170804/10527fbd/attachment.html>


More information about the users mailing list