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 [email protected]
_______________________________________________ OpenTURNS users mailing list [email protected] http://openturns.org/mailman/listinfo/users
