1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
|
Y0 > 20
Y0 < 10
Y0 < 0
Y0 < 15
Y1 < 15
#!/usr/bin/env python
import openturns as ot
import persalys
myStudy = persalys.Study('myStudy')
persalys.Study.Add(myStudy)
dist_X0 = ot.Normal(1, 1)
X0 = persalys.Input('X0', 0, dist_X0, '')
dist_X1 = ot.Normal(1, 1)
X1 = persalys.Input('X1', 0, dist_X1, '')
Y0 = persalys.Output('Y0', '')
Y1 = persalys.Output('Y1', '')
inputs = [X0, X1]
outputs = [Y0, Y1]
formulas = ['sin(X0) + 8*X1', '1 + sin(X0) + 8*X1']
aModelPhys = persalys.SymbolicPhysicalModel('aModelPhys', inputs, outputs, formulas)
myStudy.add(aModelPhys)
aLimitState = persalys.LimitState('aLimitState', aModelPhys, 'Y0', ot.Less(), 10)
myStudy.add(aLimitState)
aLimitState2 = persalys.LimitState('aLimitState2', aModelPhys, 'Y1', ot.Less(), 15)
myStudy.add(aLimitState2)
|
