[ot-users] Open Turns FORM

Douglas Long douglaslong11 at gmail.com
Fri Aug 4 19:38:43 CEST 2017


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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://openturns.org/pipermail/users/attachments/20170804/a3405179/attachment.html>


More information about the users mailing list