Auswertung SPZ Steifigkeit nach TP 26 hinzugefügt

src/pytestpavement/functions/citt.py

@@ -8,3 +8,9 @@ def stiffness_tp26(T, f, Emax, Emin, phi, z0, z1, T0=20.0):
     E = Emin + (Emax - Emin) / (1 + np.exp(z0 * x + z1))
 
     return E
+
+
+def calc_nu(T):
+    #TODO: Prüfen ob Formel stimmt!
+    nu = 0.15 + (0.35) / (1 + np.exp(3.1849 - 0.04233 * (9 / 5 * T + 32)))
+    return nu

src/pytestpavement/labtests/base.py

@@ -318,14 +318,22 @@ class DataSineLoad():
                     'Number of load cycles smaller than selectect values')
 
         if not isinstance(self.data, list):
-            df_sel = [
+            if self.number_of_load_cycles_for_analysis > 1:
-                sel_df(self.data, num=self.number_of_load_cycles_for_analysis)
+                df_sel = [
-            ]
+                    sel_df(self.data,
+                           num=self.number_of_load_cycles_for_analysis)
+                ]
+            else:
+                df_sel = [self.data]
 
         else:
             df_sel = []
             for d in self.data:
-                d_sel = sel_df(d, num=self.number_of_load_cycles_for_analysis)
+                if self.number_of_load_cycles_for_analysis > 1:
+                    d_sel = sel_df(d,
+                                   num=self.number_of_load_cycles_for_analysis)
+                else:
+                    d_sel = d
 
                 df_sel.append(d_sel)
 

src/pytestpavement/models/citt.py

@@ -2,6 +2,7 @@ import datetime
 
 from mongoengine import *
 
+from .data import CITTSiffness
 from .material import Material
 
 
@@ -41,40 +42,43 @@ class CITTSiffness(CyclicIndirectTensileTest):
     sigma_set = FloatField()
     T_set = FloatField()
 
+    N_from = IntField()
+    N_to = IntField()
+
+    data = LazyReferenceField(CITTSiffness,
+                              required=True,
+                              reverse_delete_rule=CASCADE)
+
     #results
     stiffness = FloatField()
+    nu = FloatField()
+    phase = FloatField()
     #fit parameter
     ## F
-    fit_a_F = FloatField()
+    F_amp = FloatField()
-    fit_b_F = FloatField()
+    F_freq = FloatField()
-    fit_d_F = FloatField()
+    F_phase = FloatField()
-    fit_e_F = FloatField()
+    F_offset = FloatField()
-    fit_f_F = FloatField()
+    F_slope = FloatField()
-    r2_F = FloatField()
+    F_r2 = FloatField()
     ## S1
-    fit_a_s_hor_1 = FloatField()
+    s_hor_1_amp = FloatField()
-    fit_b_s_hor_1 = FloatField()
+    s_hor_1_freq = FloatField()
-    fit_d_s_hor_1 = FloatField()
+    s_hor_1_phase = FloatField()
-    fit_e_s_hor_1 = FloatField()
+    s_hor_1_offset = FloatField()
-    fit_f_s_hor_1 = FloatField()
+    s_hor_1_slope = FloatField()
-    r2_s_hor_1 = FloatField()
+    s_hor_1_r2 = FloatField()
     ## S2
-    fit_a_s_hor_2 = FloatField()
+    s_hor_2_amp = FloatField()
-    fit_b_s_hor_2 = FloatField()
+    s_hor_2_freq = FloatField()
-    fit_d_s_hor_2 = FloatField()
+    s_hor_2_phase = FloatField()
-    fit_e_s_hor_2 = FloatField()
+    s_hor_2_offset = FloatField()
-    fit_f_s_hor_2 = FloatField()
+    s_hor_2_slope = FloatField()
-    r2_s_hor_2 = FloatField()
+    s_hor_2_r2 = FloatField()
     ## S-Sum
-    fit_a_s_hor_sum = FloatField()
+    s_hor_amp = FloatField()
-    fit_b_s_hor_sum = FloatField()
+    s_hor_freq = FloatField()
-    fit_d_s_hor_sum = FloatField()
+    s_hor_phase = FloatField()
-    fit_e_s_hor_sum = FloatField()
+    s_hor_offset = FloatField()
-    fit_f_s_hor_sum = FloatField()
+    s_hor_slope = FloatField()
-    r2_s_hor_sum = FloatField()
+    s_hor_r2 = FloatField()
-    # data
-    time = ListField(FloatField())
-    F = ListField(FloatField())
-    N = ListField(IntField())
-    s_hor_1 = ListField(FloatField())
-    s_hor_2 = ListField(FloatField())

src/pytestpavement/models/data.py

@@ -31,3 +31,15 @@ class DataSheartest(RawData):
     s_vert_2 = ListField(FloatField())
     s_vert_sum = ListField(FloatField())
     s_piston = ListField(FloatField())
+
+
+class CITTSiffness(RawData):
+
+    # data
+    time = ListField(FloatField())
+    F = ListField(FloatField())
+    N = ListField(IntField())
+    s_hor_1 = ListField(FloatField())
+    s_hor_2 = ListField(FloatField())
+    s_hor_sum = ListField(FloatField())
+    s_piston = ListField(FloatField(), required=False)