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
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
|
#!/usr/bin/env python
import openturns as ot
import persalys
Study_0 = persalys.Study("Study_0")
persalys.Study.Add(Study_0)
# variables
z0 = persalys.Input("z0", 100, "")
v0 = persalys.Input("v0", 55, "")
m = persalys.Input("m", 80, "")
c = persalys.Input("c", 15, "")
z = persalys.Output("z", "")
z2 = persalys.Output("z2", "")
inputs = [z0, v0, m, c]
outputs = [z]
# mesh model
meshModel = persalys.GridMeshModel(ot.Interval(0.0, 12.0), [20])
# Symbolic model
formula = [
"max(0, z0 - (m * 9.81 / c) * t + (m / c) * (v0 + (m * 9.81 / c)) * (1 - exp(-t * c / m)))"
]
model1 = persalys.SymbolicFieldModel("model1", meshModel, inputs, outputs, formula)
Study_0.add(model1)
f = model1.getPointToFieldFunction()
print(f([100, 55, 81, 17]))
print(f([[100, 55, 81, 17], [100, 55, 81, 18]]))
# Symbolic model : two outputs
outputs = [z, z2]
formula = [
"max(0, z0 - (m * 9.81 / c) * t + (m / c) * (v0 + (m * 9.81 / c)) * (1 - exp(-t * c / m)))"
] * 2
model2 = persalys.SymbolicFieldModel("model2", meshModel, inputs, outputs, formula)
Study_0.add(model2)
f = model2.getPointToFieldFunction()
print(f([100, 55, 81, 17]))
print(f([[100, 55, 81, 17], [100, 55, 81, 18]]))
# script
script = Study_0.getPythonScript()
print(script)
# exec(script)
|
