Georg_Brantegger/Python-DT_Slot_3
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

added second (easier) test volume

Browse Source
This commit is contained in:
Brantegger Georg
2022-05-17 10:42:17 +02:00
parent 2227599fba
commit 82780978d8
2 changed files with 260073 additions and 260053 deletions
Download Patch File Download Diff File Expand all files Collapse all files

103
functions/volume_change.py
View File

@@ -1,5 +1,4 @@
import numpy as np
# Testvolume # Testvolume
# Depth of the whole structure is constant and given by the variable d # Depth of the whole structure is constant and given by the variable d
# #
@@ -57,61 +56,85 @@ import numpy as np
def test_1_parameters():
h_1 = 10
h_2 = 5 + h_1
h_3 = 5 + h_2
h_1 = 10 x_1 = 100
h_2 = 5 + h_1 x_2 = 101
h_3 = 5 + h_2 x_3 = 30
x_1 = 100 d = 5
x_2 = 101
x_3 = 30
d = 5 vol_1 = x_1*d*h_1
vol_2 = x_1*d*(h_2-h_1)+(x_2-x_1)*d*(h_2-h_1)**2/(2*(h_2-h_1)) + vol_1
vol_3 = x_2*d*(h_3-h_2)+(x_3-x_2)*d*(h_3-h_2)**2/(2*(h_3-h_2)) + vol_2
V_1 = x_1*d*h_1 a_2 = 0.5*((x_2-x_1)*d)/(h_2-h_1)
V_2 = x_1*d*(h_2-h_1)+(x_2-x_1)*d*(h_2-h_1)**2/(2*(h_2-h_1)) + V_1 a_3 = 0.5*((x_3-x_2)*d)/(h_3-h_2)
V_3 = x_2*d*(h_3-h_2)+(x_3-x_2)*d*(h_3-h_2)**2/(2*(h_3-h_2)) + V_2
def V_h(h): b_2 = x_1*d-((x_2-x_1)*d*h_1)/(h_2-h_1)
b_3 = x_2*d-((x_3-x_2)*d*h_2)/(h_3-h_2)
c_2 = ((x_2-x_1)*d*h_1**2)/(2*(h_2-h_1))-h_1*x_1*d
c_3 = ((x_3-x_2)*d*h_2**2)/(2*(h_3-h_2))-h_2*x_2*d
return h_1,h_2,h_3,x_1,x_2,x_3,d,vol_1,vol_2,vol_3,a_2,a_3,b_2,b_3,c_2,c_3
def V_h_test_1(h):
h_1,h_2,h_3,x_1,x_2,x_3,d,vol_1,vol_2,vol_3,a_2,a_3,b_2,b_3,c_2,c_3 = test_1_parameters()
if h <= h_1: if h <= h_1:
V = x_1*d*h V = x_1*d*h
elif (h_1 < h) and (h <= h_2): elif (h_1 < h) and (h <= h_2):
V = x_1*d*(h-h_1)+(x_2-x_1)*d*(h-h_1)**2/(2*(h_2-h_1)) + V_1 V = x_1*d*(h-h_1)+(x_2-x_1)*d*(h-h_1)**2/(2*(h_2-h_1)) + vol_1
elif (h_2 < h) and (h <= h_3): elif (h_2 < h) and (h <= h_3):
V = x_2*d*(h-h_2)+(x_3-x_2)*d*(h-h_2)**2/(2*(h_3-h_2)) + V_2 V = x_2*d*(h-h_2)+(x_3-x_2)*d*(h-h_2)**2/(2*(h_3-h_2)) + vol_2
elif (h_3 < h): elif (h_3 < h):
V = x_3*d*(h-h_3) + V_3 V = x_3*d*(h-h_3) + vol_3
return V return V
a_2 = 0.5*((x_2-x_1)*d)/(h_2-h_1) def h_V_test_1(V):
a_3 = 0.5*((x_3-x_2)*d)/(h_3-h_2) h_1,h_2,h_3,x_1,x_2,x_3,d,vol_1,vol_2,vol_3,a_2,a_3,b_2,b_3,c_2,c_3 =test_1_parameters()
if V <= vol_1:
b_2 = x_1*d-((x_2-x_1)*d*h_1)/(h_2-h_1)
b_3 = x_2*d-((x_3-x_2)*d*h_2)/(h_3-h_2)
c_2 = ((x_2-x_1)*d*h_1**2)/(2*(h_2-h_1))-h_1*x_1*d
c_3 = ((x_3-x_2)*d*h_2**2)/(2*(h_3-h_2))-h_2*x_2*d
def h_V(V):
if V <= V_1:
h = V/(x_1*d) h = V/(x_1*d)
elif (V_1 < V) and (V <= V_2): elif (vol_1 < V) and (V <= vol_2):
h = (-b_2+(b_2**2-4*a_2*(c_2-(V-V_1)))**0.5)/(2*a_2) h = (-b_2+(b_2**2-4*a_2*(c_2-(V-vol_1)))**0.5)/(2*a_2)
elif (V_2 < V) and (V <= V_3): elif (vol_2 < V) and (V <= vol_3):
h = (-b_3+(b_3**2-4*a_3*(c_3-(V-V_2)))**0.5)/(2*a_3) h = (-b_3+(b_3**2-4*a_3*(c_3-(V-vol_2)))**0.5)/(2*a_3)
elif (V_3 < V): elif (vol_3 < V):
h = (V-V_3)/(x_3*d)+h_3 h = (V-vol_3)/(x_3*d)+h_3
return h return h
def show_parameters(): def test_2_parameters():
print('h_1: ', h_1) x = 10
print('h_2: ', h_2) d = 10
print('h_3: ', h_3) return x,d
print('x_1: ', x_1)
print('x_2: ', x_2) def V_h_test_2(h):
print('x_3: ', x_3) x,d = test_2_parameters()
return x*d*h
def h_V_test_2(V):
x,d = test_2_parameters()
return V/(x*d)
def show_parameters(test_version):
h_1,h_2,h_3,x_1,x_2,x_3,d,vol_1,vol_2,vol_3,a_2,a_3,b_2,b_3,c_2,c_3 = test_1_parameters()
x,d = test_2_parameters()
if test_version == 1:
print('h_1: ', h_1)
print('h_2: ', h_2)
print('h_3: ', h_3)
print('x_1: ', x_1)
print('x_2: ', x_2)
print('x_3: ', x_3)
elif test_version == 2:
print('x: ', x)
print('d: ', d)

520023
functions/volume_change_nb.ipynb
View File

File diff suppressed because it is too large Load Diff