diff --git a/Plot-Scripts/WoodBilayer_ExperimentComparison.py b/Plot-Scripts/WoodBilayer_ExperimentComparison.py
index 2a2d68d0700f6464ca2d763ce82df78b573f3c75..81963fc0ce292266694cb0ccb936e02fd342c4a1 100644
--- a/Plot-Scripts/WoodBilayer_ExperimentComparison.py
+++ b/Plot-Scripts/WoodBilayer_ExperimentComparison.py
@@ -55,42 +55,65 @@ for dataset_number in dataset_numbers:
# fig = plt.figure() #main
fig, ax = plt.subplots(figsize=(width,height))
+
+ # sampled with N=500 and length=5 via 'PolaPlotLocalEnergy.py'
+
dataset = [
[ # Dataset Ratio r = 0.12
[14.70179844, 13.6246 , 12.42994508, 11.69773413,11.14159987, 9.500670278, 9.005046347] , # input (moisture-content in %)
- [1.29323308, 1.83458647 , 2.40601504 , 2.76691729 , 3.03759398, 3.81954887 , 4.03007519], # curvature kappa from Simulation
- [1.140351217, 1.691038688, 2.243918105, 2.595732726, 2.945361006,4.001528043, 4.312080261] # curvature kappa from Experiment
+ # [1.30260521, 1.83366733, 2.40480962, 2.76553106, 3.0260521, 3.80761523, 4.03807615], # curvature kappa from Simulation GridLevel=3
+ # [1.30260521, 1.83366733, 2.41482966, 2.76553106, 3.03607214, 3.81763527,4.04809619], # curvature kappa from Simulation GridLevel=4
+ [1.28256513, 1.79358717, 2.36472946, 2.70541082, 2.9759519, 3.73747495, 3.96793587], # curvature kappa from Simulation GridLevel=5
+ [1.140351217, 1.691038688, 2.243918105, 2.595732726, 2.945361006,4.001528043, 4.312080261] # curvature kappa from Experiment
],
[ # Dataset Ratio r = 0.17
- [14.75453569,13.71227639,12.54975012,11.83455959,11.29089521,9.620608917,9.101671742], # input (moisture-content in %)
- [1.26315789, 1.77443609, 2.34586466, 2.70676692, 2.97744361, 3.78947368,4.03007519], # curvature kappa from Simulation
- [1.02915975,1.573720805,2.407706364,2.790518802,3.173814476,4.187433094,4.511739072] # curvature kappa from Experiment
+ [14.75453569,13.71227639,12.54975012,11.83455959,11.29089521,9.620608917,9.101671742], # input (moisture-content in %)
+ # [1.27254509, 1.78356713, 2.36472946, 2.70541082, 2.9759519, 3.78757515, 4.02805611] # curvature kappa from Simulation GridLevel=3
+ # [1.29258517, 1.81362725, 2.39478958, 2.74549098, 3.01603206, 3.83767535, 4.08817635], # curvature kappa from Simulation GridLevel=4
+ [1.28256513, 1.80360721, 2.38476954, 2.73547094, 3.00601202, 3.81763527,4.06813627], # curvature kappa from Simulation GridLevel=5
+ [1.02915975,1.573720805,2.407706364,2.790518802,3.173814476,4.187433094,4.511739072] # curvature kappa from Experiment
],
+
+
[ # Dataset Ratio r = 0.22
- [14.72680026, 13.64338887, 12.41305478, 11.66482931, 11.09781471, 9.435795985, 8.959564147], # input (moisture-content in %)
- [1.20401338 , 1.72575251 , 2.28762542 , 2.64882943 , 2.92976589 , 3.73244147 , 3.93311037 ], # curvature kappa from Simulation
- [1.058078122, 1.544624544, 2.317033799, 2.686043143, 2.967694189, 3.913528418, 4.262750825] # curvature kappa from Experiment
+ [14.72680026, 13.64338887, 12.41305478, 11.66482931, 11.09781471, 9.435795985, 8.959564147], # input (moisture-content in %)
+ # [1.19238477, 1.71342685, 2.30460922, 2.66533066, 2.93587174, 3.71743487,3.94789579], # curvature kappa from Simulation GridLevel=3
+ # [1.20240481, 1.73346693, 2.3246493, 2.68537074, 2.95591182, 3.74749499,3.97795591], # curvature kappa from Simulation GridLevel=4
+ [1.19238477, 1.72344689, 2.31462926, 2.6753507, 2.94589178, 3.72745491,3.95791583], # curvature kappa from Simulation GridLevel=5
+ [1.058078122, 1.544624544, 2.317033799, 2.686043143, 2.967694189, 3.913528418, 4.262750825] # curvature kappa from Experiment
],
+
[ # Dataset Ratio r = 0.34
- [14.98380876,13.97154915,12.77309253,12.00959929,11.42001731,9.561447179,8.964704969], # input (moisture-content in %)
- [0.87218045, 1.26315789, 1.7443609, 2.07518797, 2.28571429, 3.03759398, 3.27819549], # curvature kappa from Simulation
- [0.789078472,1.1299263,1.738136936,2.159520896,2.370047499,3.088299431,3.18097558] # curvature kappa from Experiment
+ [14.98380876,13.97154915,12.77309253,12.00959929,11.42001731,9.561447179,8.964704969], # input (moisture-content in %)
+ # [0.87174349, 1.27254509, 1.75350701, 2.06412826, 2.29458918, 3.03607214,3.26653307], # curvature kappa from Simulation GridLevel=3
+ # [0.87174349, 1.28256513, 1.76352705, 2.0741483, 2.31462926, 3.05611222,3.28657315], # curvature kappa from Simulation GridLevel=4
+ [0.87174349, 1.27254509, 1.75350701, 2.06412826, 2.29458918, 3.03607214,3.26653307], # curvature kappa from Simulation GridLevel=5
+ [0.789078472,1.1299263,1.738136936,2.159520896,2.370047499,3.088299431,3.18097558] # curvature kappa from Experiment
],
+
+
+
[ # Dataset Ratio r = 0.43
- [15.11316339,14.17997082,13.05739844,12.32309209,11.74608518,9.812372466,9.10519385 ], # input (moisture-content in %)
- [0.63157895, 0.93233083, 1.29323308, 1.53383459, 1.71428571, 2.34586466,2.55639098], # curvature kappa from Simulation
- [0.577989364,0.829007544,1.094211707,1.325332511,1.400455154,1.832325697,2.047483977] # curvature kappa from Experiment
+ [15.11316339,14.17997082,13.05739844,12.32309209,11.74608518,9.812372466,9.10519385 ], # input (moisture-content in %)
+ # [0.63126253, 0.93186373, 1.29258517, 1.53306613, 1.71342685, 2.33466934,2.55511022], # curvature kappa from Simulation GridLevel=3
+ # [0.61122244, 0.90180361, 1.25250501, 1.48296593, 1.66332665, 2.26452906, 2.48496994], # curvature kappa from Simulation GridLevel=4
+ [0.6012024, 0.89178357, 1.23246493, 1.46292585, 1.63326653, 2.2244489,2.44488978], # curvature kappa from Simulation GridLevel=5
+ [0.577989364,0.829007544,1.094211707,1.325332511,1.400455154,1.832325697,2.047483977] # curvature kappa from Experiment
],
+
+
[ # Dataset Ratio r = 0.49
- [15.30614414,14.49463867,13.46629742,12.78388234,12.23057715,10.21852839,9.341730605], # input (moisture-content in %)
- [0.60150376, 0.87218045, 1.23308271, 1.5037594 , 1.71428571, 2.46616541,2.79699248], # curvature kappa from Simulation
- [0.357615902,0.376287785,0.851008627,0.904475291,1.039744708,1.346405241,1.566568558] # curvature kappa from Experiment
+ [15.30614414,14.49463867,13.46629742,12.78388234,12.23057715,10.21852839,9.341730605], # input (moisture-content in %)
+ # [0.59118236, 0.87174349, 1.24248497, 1.49298597, 1.70340681, 2.46492986,2.79559118], # curvature kappa from Simulation GridLevel=3
+ # [0.54108216, 0.80160321, 1.14228457, 1.37274549, 1.56312625, 2.26452906,2.5751503 ], # curvature kappa from Simulation GridLevel=4
+ [0.56112224, 0.84168337, 1.19238477, 1.43286573, 1.63326653, 2.36472946,2.68537074], # curvature kappa from Simulation GridLevel=5
+ [0.357615902,0.376287785,0.851008627,0.904475291,1.039744708,1.346405241,1.566568558] # curvature kappa from Experiment
]
]
diff --git a/experiment/wood-bilayer/PolarPlotLocalEnergy.py b/experiment/wood-bilayer/PolarPlotLocalEnergy.py
index 0abd7e9f0ca9863047aa806b6c3672f3da33b38c..ae5fe253542151f516001ba2078032a9c53ca6d7 100644
--- a/experiment/wood-bilayer/PolarPlotLocalEnergy.py
+++ b/experiment/wood-bilayer/PolarPlotLocalEnergy.py
@@ -55,8 +55,8 @@ def ReadEffectiveQuantities(QFilePath, BFilePath):
number=7
show_plot = False
-# dataset_numbers = [0, 1, 2, 3, 4, 5]
-dataset_numbers = [0]
+dataset_numbers = [0, 1, 2, 3, 4, 5]
+# dataset_numbers = [0]
for dataset_number in dataset_numbers:
@@ -79,8 +79,8 @@ for dataset_number in dataset_numbers:
B=np.transpose([B])
#
- N=400
- length=12
+ N=500
+ length=5
r, theta = np.meshgrid(np.linspace(0,length,N),np.radians(np.linspace(0, 360, N)))
E=np.zeros(np.shape(r))
for i in range(0,N):