diff --git a/Plot-Scripts/WoodBilayer_expComparison_r049.py b/Plot-Scripts/WoodBilayer_expComparison_r049.py index f09281bfe4018346f1f5f7cf94466297dfa5bf87..818b705a2b1f52a97cdeb9eaa7da946a7ba3d946 100644 --- a/Plot-Scripts/WoodBilayer_expComparison_r049.py +++ b/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 diff --git a/experiment/wood-bilayer/wood_test.py b/experiment/wood-bilayer/wood_test.py index 5844cf4cf1d9941a0e05a2fc980bd13eb654bbd2..a82df9a1e7be8527e6f3c4a8109ce8e7a2a2004d 100644 --- a/experiment/wood-bilayer/wood_test.py +++ b/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