Add dataset for ratio r=0.49 and plot script

Plot-Scripts/WoodBilayer_expComparison_r049.py

@@ -50,11 +50,13 @@ fig, ax = plt.subplots(figsize=(width,height))
 
 
 
-input = [14.72680026, 13.64338887, 12.41305478, 11.66482931, 11.09781471, 9.435795985, 8.959564147]  #moisture-content (in %)
+input = [15.30614414,14.49463867,13.46629742,12.78388234,12.23057715,10.21852839,9.341730605] #moisture-content (in %)
+
+
+kappa_sim = [0.59118236, 0.87174349, 1.24248497, 1.49298597, 1.70340681, 2.46492986, 2.79559118]         #kappa from Simulation
+kappa_exp = [0.357615902, 0.376287785, 0.851008627, 0.904475291, 1.039744708, 1.346405241, 1.566568558]   #kappa from Experiment
 
 
-kappa_sim = [1.20401338 , 1.72575251, 2.28762542 , 2.64882943 , 2.92976589 , 3.73244147 , 3.93311037 ]    #kappa from Simulation
-kappa_exp = [1.058078122,1.544624544, 2.317033799, 2.686043143, 2.967694189, 3.913528418, 4.262750825]    #kappa from Experiment
 
 # compute difference:
 relative_error = (np.array(kappa_sim) - np.array(kappa_exp)) / np.array(kappa_exp)
@@ -97,7 +99,7 @@ line_2 = ax.plot(input, kappa_exp,                     # data
 # Fix_value = 7.674124
 # l3 = plt.axhline(y = Fix_value, color = 'black', linewidth=0.75, linestyle = 'dashed')
 # --------------- Set Axes  -----------------------
-ax.set_title(r"ratio $r = 0.22$")   # Plot - Title
+ax.set_title(r"ratio $r = 0.49$")   # Plot - Title
 # plt.xscale('log') # Use Logarithmic-Scale
 # plt.yscale('log')
 ax.set_xlabel(r"Wood moisture content $\omega (\%)$", labelpad=4)
@@ -130,5 +132,5 @@ frame.set_edgecolor('black')
 
 # ---------- Output Figure as pdf:
 fig.set_size_inches(width, height)
-fig.savefig('WoodBilayer_expComparison.pdf')
+fig.savefig('WoodBilayer_expComparison_r049.pdf')
 # plt.show()
\ No newline at end of file

experiment/wood-bilayer/wood_test.py

@@ -124,14 +124,26 @@ gamma = 1.0
 # experimental_kappa = curvature measure in experiment
 
 #First Experiment:
+# Ratio r = 0.22
+# materialFunctionParameter=[
+#    [0.22, 0.0053, 17.17547062, 14.72680026, 0, 1.058078122],
+#    [0.22, 0.0053, 17.17547062, 13.64338887, 0, 1.544624544],
+#    [0.22, 0.0053, 17.17547062, 12.41305478, 0, 2.317033799],
+#    [0.22, 0.0053, 17.17547062, 11.66482931, 0, 2.686043143],
+#    [0.22, 0.0053, 17.17547062, 11.09781471, 0, 2.967694189],
+#    [0.22, 0.0053, 17.17547062, 9.435795985, 0, 3.913528418],
+#    [0.22, 0.0053,  17.17547062, 8.959564147, 0, 4.262750825]
+# ]
+
+# Ratio r = 0.49
 materialFunctionParameter=[
-   [0.22, 0.0053, 17.17547062, 14.72680026, 0, 1.058078122],
-   [0.22, 0.0053, 17.17547062, 13.64338887, 0, 1.544624544],
-   [0.22, 0.0053, 17.17547062, 12.41305478, 0, 2.317033799],
-   [0.22, 0.0053, 17.17547062, 11.66482931, 0, 2.686043143],
-   [0.22, 0.0053, 17.17547062, 11.09781471, 0, 2.967694189],
-   [0.22, 0.0053, 17.17547062, 9.435795985, 0, 3.913528418],
-   [0.22, 0.0053,  17.17547062, 8.959564147, 0, 4.262750825]
+   [0.49, 0.008,  17.01520754, 15.30614414, 0, 0.357615902],
+   [0.49, 0.008,  17.01520754, 14.49463867, 0, 0.376287785],
+   [0.49, 0.008,  17.01520754, 13.46629742, 0, 0.851008627],
+   [0.49, 0.008,  17.01520754, 12.78388234, 0, 0.904475291],
+   [0.49, 0.008,  17.01520754, 12.23057715, 0, 1.039744708],
+   [0.49, 0.008,  17.01520754, 10.21852839, 0, 1.346405241],
+   [0.49, 0.008,  17.01520754, 9.341730605, 0, 1.566568558]
 ]
 
 #Second Experiment: Rotate "active" bilayer phase