From 41316745787e8250fc10509ad918782f4be2143b Mon Sep 17 00:00:00 2001
From: Stefan Neukamm <stefan.neukamm@tu-dresden.de>
Date: Wed, 13 Sep 2023 09:16:09 +0200
Subject: [PATCH] add wood
---
.../rotation-test/PolarPlotLocalEnergy.py | 8 +--
experiment/rotation-test/cellsolver.parset | 4 +-
.../isotrop_orthotrop_rotation.py | 5 +-
experiment/rotation-test/rotation_test.py | 5 +-
experiment/wood-bilayer/#readme.txt# | 12 +++++
experiment/wood-bilayer/.gitignore | 0
.../wood-bilayer/PolarPlotLocalEnergy.py | 2 +
.../elasticity_toolbox.cpython-310.pyc | Bin 0 -> 3564 bytes
.../wood_european_beech.cpython-310.pyc | Bin 0 -> 2977 bytes
experiment/wood-bilayer/readme.txt | 0
experiment/wood-bilayer/results/0/BMatrix.txt | 3 ++
experiment/wood-bilayer/results/0/QMatrix.txt | 9 ++++
experiment/wood-bilayer/results/0/output.txt | 49 ++++++++++++++++++
.../wood-bilayer/results/0/parameter.txt | 4 ++
experiment/wood-bilayer/results/1/BMatrix.txt | 3 ++
experiment/wood-bilayer/results/1/QMatrix.txt | 9 ++++
experiment/wood-bilayer/results/1/output.txt | 49 ++++++++++++++++++
.../wood-bilayer/results/1/parameter.txt | 4 ++
experiment/wood-bilayer/results/2/BMatrix.txt | 3 ++
experiment/wood-bilayer/results/2/QMatrix.txt | 9 ++++
experiment/wood-bilayer/results/2/output.txt | 49 ++++++++++++++++++
.../wood-bilayer/results/2/parameter.txt | 4 ++
experiment/wood-bilayer/results/3/BMatrix.txt | 3 ++
experiment/wood-bilayer/results/3/QMatrix.txt | 9 ++++
experiment/wood-bilayer/results/3/output.txt | 49 ++++++++++++++++++
.../wood-bilayer/results/3/parameter.txt | 4 ++
experiment/wood-bilayer/results/4/BMatrix.txt | 3 ++
experiment/wood-bilayer/results/4/QMatrix.txt | 9 ++++
experiment/wood-bilayer/results/4/output.txt | 49 ++++++++++++++++++
.../wood-bilayer/results/4/parameter.txt | 4 ++
experiment/wood-bilayer/results/5/BMatrix.txt | 3 ++
experiment/wood-bilayer/results/5/QMatrix.txt | 9 ++++
experiment/wood-bilayer/results/5/output.txt | 49 ++++++++++++++++++
.../wood-bilayer/results/5/parameter.txt | 4 ++
experiment/wood-bilayer/results/6/BMatrix.txt | 3 ++
experiment/wood-bilayer/results/6/QMatrix.txt | 9 ++++
experiment/wood-bilayer/results/6/output.txt | 49 ++++++++++++++++++
.../wood-bilayer/results/6/parameter.txt | 4 ++
experiment/wood-bilayer/results/auswertung.py | 19 +++++--
.../wood-bilayer/results/kappa_simulation.txt | 2 +
.../wood-bilayer/wood_european_beech.py | 49 ++++++++----------
experiment/wood-bilayer/wood_test.py | 12 +++--
42 files changed, 527 insertions(+), 46 deletions(-)
create mode 100644 experiment/wood-bilayer/#readme.txt#
create mode 100644 experiment/wood-bilayer/.gitignore
create mode 100644 experiment/wood-bilayer/__pycache__/elasticity_toolbox.cpython-310.pyc
create mode 100644 experiment/wood-bilayer/__pycache__/wood_european_beech.cpython-310.pyc
create mode 100644 experiment/wood-bilayer/readme.txt
create mode 100644 experiment/wood-bilayer/results/0/BMatrix.txt
create mode 100644 experiment/wood-bilayer/results/0/QMatrix.txt
create mode 100644 experiment/wood-bilayer/results/0/output.txt
create mode 100644 experiment/wood-bilayer/results/0/parameter.txt
create mode 100644 experiment/wood-bilayer/results/1/BMatrix.txt
create mode 100644 experiment/wood-bilayer/results/1/QMatrix.txt
create mode 100644 experiment/wood-bilayer/results/1/output.txt
create mode 100644 experiment/wood-bilayer/results/1/parameter.txt
create mode 100644 experiment/wood-bilayer/results/2/BMatrix.txt
create mode 100644 experiment/wood-bilayer/results/2/QMatrix.txt
create mode 100644 experiment/wood-bilayer/results/2/output.txt
create mode 100644 experiment/wood-bilayer/results/2/parameter.txt
create mode 100644 experiment/wood-bilayer/results/3/BMatrix.txt
create mode 100644 experiment/wood-bilayer/results/3/QMatrix.txt
create mode 100644 experiment/wood-bilayer/results/3/output.txt
create mode 100644 experiment/wood-bilayer/results/3/parameter.txt
create mode 100644 experiment/wood-bilayer/results/4/BMatrix.txt
create mode 100644 experiment/wood-bilayer/results/4/QMatrix.txt
create mode 100644 experiment/wood-bilayer/results/4/output.txt
create mode 100644 experiment/wood-bilayer/results/4/parameter.txt
create mode 100644 experiment/wood-bilayer/results/5/BMatrix.txt
create mode 100644 experiment/wood-bilayer/results/5/QMatrix.txt
create mode 100644 experiment/wood-bilayer/results/5/output.txt
create mode 100644 experiment/wood-bilayer/results/5/parameter.txt
create mode 100644 experiment/wood-bilayer/results/6/BMatrix.txt
create mode 100644 experiment/wood-bilayer/results/6/QMatrix.txt
create mode 100644 experiment/wood-bilayer/results/6/output.txt
create mode 100644 experiment/wood-bilayer/results/6/parameter.txt
create mode 100644 experiment/wood-bilayer/results/kappa_simulation.txt
diff --git a/experiment/rotation-test/PolarPlotLocalEnergy.py b/experiment/rotation-test/PolarPlotLocalEnergy.py
index 67429fda..5cb9d8de 100644
--- a/experiment/rotation-test/PolarPlotLocalEnergy.py
+++ b/experiment/rotation-test/PolarPlotLocalEnergy.py
@@ -48,7 +48,7 @@ def ReadEffectiveQuantities(QFilePath, BFilePath):
B = np.array([X[0][2], X[1][2], X[2][2]])
return Q, B
-number=2
+number=1
kappa=np.zeros(number)
for n in range(0,number):
# Read from Date
@@ -62,8 +62,8 @@ for n in range(0,number):
B=np.transpose([B])
#
- N=500
- length=5
+ N=300
+ 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):
@@ -87,7 +87,7 @@ for n in range(0,number):
plt.colorbar(pcm, extend='max', ticks=colorbarticks, pad=0.1)
plt.show()
-f = open("./experiment/wood-bilayer/results/kappa_simulation.txt", "w")
+f = open("./experiment/rotation-test/results/kappa_simulation.txt", "w")
f.write(str(kappa))
f.close()
diff --git a/experiment/rotation-test/cellsolver.parset b/experiment/rotation-test/cellsolver.parset
index e6745496..065f1d15 100644
--- a/experiment/rotation-test/cellsolver.parset
+++ b/experiment/rotation-test/cellsolver.parset
@@ -4,7 +4,7 @@
# --------------------------------------------------------
# Path for results and logfile
-outputPath=./experiment/rotation-test/results/0
+outputPath=./experiment/rotation-test/results/3
# Path for material description
geometryFunctionPath =experiment/rotation-test/
@@ -24,7 +24,7 @@ geometryFunctionPath =experiment/rotation-test/
## {start,finish} computes on all grid from 2^(start) to 2^finish refinement
#----------------------------------------------------
-numLevels=4 4
+numLevels=3 3
#numLevels = 1 1 # computes all levels from first to second entry
#numLevels = 2 2 # computes all levels from first to second entry
#numLevels = 1 3 # computes all levels from first to second entry
diff --git a/experiment/rotation-test/isotrop_orthotrop_rotation.py b/experiment/rotation-test/isotrop_orthotrop_rotation.py
index d4809134..1669e604 100644
--- a/experiment/rotation-test/isotrop_orthotrop_rotation.py
+++ b/experiment/rotation-test/isotrop_orthotrop_rotation.py
@@ -25,7 +25,7 @@ prestrain_wood=np.array([[-0.050821460301886785, 0, 0],
[0, -2.134501332679245, 0],
[0, 0, -0.8824453561509432]])
# rotation angle
-param_theta = 0
+param_theta = 1.5707963267948966
# untere Schicht: isotrop, ohne prestrain
@@ -49,7 +49,8 @@ materialParameters_phase1 = np.dot(np.dot(N,compliance_wood),N.T)
materialParameters_phase1 = (0.5*(materialParameters_phase1.T+materialParameters_phase1)).tolist()
# rotation of strain
def prestrain_phase1(x):
- return elast.voigt_to_strain(np.dot(elast.rotation_matrix_compliance(2,param_theta),np.dot(elast.strain_to_voigt(np.array(prestrain_wood)),N.T))).tolist()
+ R=elast.rotation_matrix(2,param_theta)
+ return np.dot(R,np.dot(np.array(prestrain_wood),R.T)).tolist()
# --- PHASE 2
phase2_type="isotropic"
diff --git a/experiment/rotation-test/rotation_test.py b/experiment/rotation-test/rotation_test.py
index 34061707..07d889e3 100644
--- a/experiment/rotation-test/rotation_test.py
+++ b/experiment/rotation-test/rotation_test.py
@@ -107,13 +107,14 @@ outputPath = Path + '/results/'
class ParameterSet:
pass
+#ParameterSet.materialFunction = "isotrop_orthotrop_rotation"
ParameterSet.materialFunction = "isotrop_orthotrop_rotation"
ParameterSet.gamma=1.0
-ParameterSet.numLevels=4
+ParameterSet.numLevels=3
# ----- Define Parameters for Material Function --------------------
# Liste mit Drehwinkel theta
-materialFunctionParameter=[0,np.pi/12, 2*np.pi/12,3*np.pi/12, 4*np.pi/12, 5*np.pi/12]
+materialFunctionParameter=[0, 2*np.pi/12, 4*np.pi/12, 6*np.pi/12]
# ------ Loops through Parameters for Material Function -----------
for i in range(0,np.shape(materialFunctionParameter)[0]):
diff --git a/experiment/wood-bilayer/#readme.txt# b/experiment/wood-bilayer/#readme.txt#
new file mode 100644
index 00000000..81861ace
--- /dev/null
+++ b/experiment/wood-bilayer/#readme.txt#
@@ -0,0 +1,12 @@
+Zum Starten im Verzeichnis dune-microstructure:
+python ./experimenta/wood_text.py
+(sonst kommt das script mit den Verzeichnissen durcheinander)
+
+Definition des Komposits:
+wood_europea_beech.py
+
+Zum Plotten:
+PolarPlotLocalEnergy.py
+
+Toolbox
+elasticity_toolboy.py
\ No newline at end of file
diff --git a/experiment/wood-bilayer/.gitignore b/experiment/wood-bilayer/.gitignore
new file mode 100644
index 00000000..e69de29b
diff --git a/experiment/wood-bilayer/PolarPlotLocalEnergy.py b/experiment/wood-bilayer/PolarPlotLocalEnergy.py
index fae4dbb6..8a7cc06c 100644
--- a/experiment/wood-bilayer/PolarPlotLocalEnergy.py
+++ b/experiment/wood-bilayer/PolarPlotLocalEnergy.py
@@ -48,12 +48,14 @@ def ReadEffectiveQuantities(QFilePath, BFilePath):
B = np.array([X[0][2], X[1][2], X[2][2]])
return Q, B
+# Number of experiments / folders
number=7
kappa=np.zeros(number)
for n in range(0,number):
# Read from Date
print(str(n))
DataPath = './experiment/wood-bilayer/results/'+str(n)
+ DataPath = './experiment/wood-bilayer/results/'+str(n)
QFilePath = DataPath + '/QMatrix.txt'
BFilePath = DataPath + '/BMatrix.txt'
Q, B = ReadEffectiveQuantities(QFilePath,BFilePath)
diff --git a/experiment/wood-bilayer/__pycache__/elasticity_toolbox.cpython-310.pyc b/experiment/wood-bilayer/__pycache__/elasticity_toolbox.cpython-310.pyc
new file mode 100644
index 0000000000000000000000000000000000000000..23ea9ef696e037927165d5e533a22f55855c9a95
GIT binary patch
literal 3564
zcmd1j<>g{vU|<M5x;AAMKLf*K5C<7EF)%PVFfcF_`!Fytq%fo~<}gG-XvQd}6s8pB
z7KSM16qXd$7KSL66t)!h7KSL+6pj?m7KSLc6s{ER7KSKxus)D_o)q2|hA55{#$X0b
zzLy|-{E`_#6clqXFff3yGswmU1_p*2hAf5@#%$&y&JxBHre?+_#u|ncW=V!x5X~aV
zP|H-qkisg-P|KXcR>P3OF3C{CkisF!AkHAdP|E@q;{=OwLd96YVq9P`E~pq=3DW}R
z8pax?8s-|78rB-Nbhd?zwd^(Q!3>(*ek++?GBYqR6oH5;ZO^=d(h`N-#FC=S3Wef~
z{L-8hg``x4L<QptV};_9qQuNRsI1;frXrBvz{E<XmtZy<0|UcKrk9{lEdqHo9pdo^
z_A8l+KyqNgl}um?EC4eEq!`8qDb!@S#af(^Sde;)D=aamG}X1JD8J|yQ(nO>*2JQs
z#7a$$Tf7kG#6ukWg_VJU;oR2viV7I3vJ%Fss)DdK$5+Bx6$n=4Exxk+%=D7@y!?{H
zlFa<Pl?+9a3=9mvg7q`<b5r$;OH$Jk^Yl|n^HRZtZf<6BQGRhrQE75XX;G?vYDGb6
zQD$yxUWtBretwE>Qf5wKWonUrYEEKtNoI0pNo9OVetu3;euZ8^<t=`wvrF>h!M1|j
zT`a=Dz`(`C!N|eP!N|cV0)}GDEKEhx3=9kyNr;_+fdQO^vKSZ`N*HPwni;`~h>K7n
zs$ocBlVqr6fh46`R!Cy1g(VQS8rF2?8nzmy8kThC8rB+?8s>E7bY^&B;q-e6ikm7k
zNGd5yO)kkVf~As-#IjTcGlkTgR8UA1D`e&=gn>d;0T!xy;FJQ6bx8b!Qrt@x1_p+g
zpfJ#6xy785npeaJN{sv<LV$sRL6f6M0>l;q5yBt>6jVi`AeI=25C;*0AOaL+;CL&N
zVqjpn#SaeM_>%m1NPvSB7K3DX;Nc7kYYs*paG10GV__*$0EIRv@j{X<0|UsrAPf%e
zE#S~zz*xh$kg=Ajgkb>_n8jQJ%EC-VvLy@)n9~_jm=-anGo&zsNtQ*7;P7L~Qv_w?
zY^EXus46zFJVY%!TooUREW`|u8c-atfaE}#n42MmF`c=FxrPZuLu~}9U@c)-08+(L
z!;-~X!@7{EhGijBFoPzCA2{eWnZXo<0Qm|g0AhnhHJM-j|Ns9#m<=I7Qku+=NCW3^
z5CPMt2nu$v09z3##%?hu=NI2%F3!v=0>#}e=9K)BTZ|!^Y`2&aD{iqS=B4MP-r`Ko
zFODzCNG(ac#aW!02W2t_-eL@b#2z>vK@oL}zX+Nkp?O>fRBUR2vI-*?BL^cFBO4<d
zBga1u78r|*6(qv;kB^lDB==8@xd@~u897RM7#J8pjs-{Q8b(M4s$pzqs$~KND5$Kh
zVM<{FXA>rHS<G6)lmaS?*=m@;WifjV6F3`j)G&cc<XX-eCQzAN%T>da!YRp6%U!~_
zfT@N#g-e*Bh6PNrf=QMdRuP68W+=%b!jR5Z!v>aN2a_CNlD&pQgrSBFO0t7hae`&I
zz$7=A<f`EoVW{DRl3XBFDO@#dATbe!6s{U}5r!HrFbiZFNG&^9EeA+c9ArL7jH8BA
zgaORs5@D!ehw#$bQn+fE;U=-bO#-<dq>BY?5-ZdsW+)G05;I7&h6T+eHn2-TCb5HA
zAd^7uVgdV-735lRP{@G9SV6u3^FV%Jf$-A7br!h3Dgu=lMW9$L;s)0?NQo0vjudf1
zL_swQBn^XF07alO6I|?QLec~{!Ge>iCQp$XsK8(aC(0si5St}1zM|q5h^z#YRUonw
zEK&&;sRWBufkmpoB2`uB$yW+3`Nk*b=N9B-Cgvrlx`0x$3OFT03L#{Yg;5O*y%<4|
z?N^Z{s8|A(mat+86a=802~NXu;9?1s=1Q1rm_RAHgauS3EM#hCtYuDN2xib^^wVSo
z7j9r*-(t&6EXgg+xh0P7z|?}`%$)o@aIFI>k5dbZ!C_OR4+<1UCv+e3f(=N8maq{Z
zPl2+U9Al9k$jP9h0p?_oILO&xCpUq!-2%oGP)Q1E>M%(%)G($n+kl!kERqZ?ARem?
zh-3qcu|vf;z+#+WF)pYWH&~1ZEXE5J;{%KFgT(}(VuD~XA+VS**gg?3TNElU29_5G
zi%CGmB*9`*U@>W^m<(7<7Az(Q6_W>xDS*Wk!DcIg*~(CP6|lT2SWFEnrVbX<0E=lt
z#k9a;+F&squsOP5wjNYoA1rSG7Bhs38G*%&!D1$0J*Hr`8C2dJEN=l8vjodqfkc8C
zG_CzKnTljUg*cc1MHeJ@fzlQvi|K-7!C48E*^2Z)Tu@5Xgf{-5Y!C&>Y_gzA8!8E9
zgD8k3D5=79LfIe+A_;Z_C<ng;=T%TPE&^rUmmupkq3uVIMVd_D#ws}9zXU0RNrEiW
zWGVtt5J`}Bm?X#|O{O9c1(5`M2xQeuP+_78ZImKg0J8WcFT_J2HrQmC1t5W!pvDSJ
zC&*-&1z<@~EmH(a_b)-UA<Wev*TKvHS@sgtc!5bGn*o*twR>QyK{{W8+CVS?WYu6v
zP}>OR1gIp)5w}>YQj7A7HJNVlKpI7mCi5+}lGME7{37&Rjc7SR<YPhk8Wi(fj66(`
zwke22uwg6+4I+7%SeS|eKq1F@izPR)B;yuqUTJPYCAj<n7cY?F2vj?Ri#KqA0oH-!
n8bnKg!zMRBr8Fni4wOcVL4m@-0B$V{afotAaENd)^DqGbVy~?~
literal 0
HcmV?d00001
diff --git a/experiment/wood-bilayer/__pycache__/wood_european_beech.cpython-310.pyc b/experiment/wood-bilayer/__pycache__/wood_european_beech.cpython-310.pyc
new file mode 100644
index 0000000000000000000000000000000000000000..12b7e98a66b84afe9b00cb5b987c3f42eb6e8be7
GIT binary patch
literal 2977
zcmd1j<>g{vU|{$dwl1YdoPpsnh=Yvz7#J8F7#J9er!X@xq%fo~<}gG-XvQc;FrO)k
z2~0CbF@tHAC{{4d7R8>z6vdIk9L1T!62+Cm8pWN$5XF<qoysM_5XGCy9mSW*x`020
zErmUuDTO12Go2}gD}_6qDTOD6H=QYkFNHsyDMcVfFr6tyC`CA(DMchjG@U6$EJZw>
zDV=j6W0XLuK#GJUL#jZEq$EQsp9q6ELzG~uK#CMtOd2jGlq!%S0~V8oiwUO+q{xBA
z<l$l>sRAhqU@=9wm}sg%iV|2%87?N4Dv+WA7E^_biKhyrsDZ`Q;bIc00x23`F-^Fb
zWU6GUkO)Jnurx!IRH}HYpa?^%kTgS-bgBfH4HB10(Tb8y(T<W!(TS2z(M{27VTe*l
z(T`G0F^E!16;2gO6-;N2VUALcQi)Pcl}wdL6$kUxqST`_Ql(O5Ql-Is%_yxX?NqT;
z(NqyIUnfd8N)H_7hN;TZ45^yj45=#845?b&45_Np3@Jt+n9i1}#?1hdPca0`YNN@j
zLu6AmAi7fwk@bR1NHIz=Ox1&!m(B(i0g35?#E?WlVmi_c>1^q23mK#IQ}v?^QjAkf
zS{R}XQw*YvQUyVN1-l&LQ<+rhR7o%&;!}xK@l+`=AL3KdRFPCMFn{p-)Hua7#jJ%P
z$|S`g$~2XIfmw=qD$hcuDDza-1r}gzo@%g=Db+m7av@WcRf<KbRhDI{l{kY4LzGpD
zWr|e`LzHz2V=#lJ^-EBG^-I@`=|A*U<BWa!!t~9rino^9rwjPkEqSeD=#c)p@AY5r
z*NP75#vs7G#36mw%$YN1-t;{H5n0{ske&tx<_FTxo;h<yXU_$w$g;u%5UG1D4(Si8
zMAo?{>K({D#mc~t-ugOlhg_@PflMD328Q$<av;9pf%IEze#PV-y=0$$^z`MOi>3tZ
zPY>vwo+s3N!#;h_CT_OtcN+Jn+w2XRrr3AYKK)VA9;S<voAzrm6)`d}Fn|dt`z1&~
zlc|UaNgT`u5t>XdK@?0qGe`|cK$Gbum<`p-f+P-RL&aH<#KCN+I2)2Um<<(YM-m6K
zq2e4!;$SvZoD(LV?!Ccc&d;kd_UY5ZoFB<~$J?iGG;TeY{ieu18Kf2zbs!9KJILYA
zprTfWfq@}~v6iuhA%#hjA%(ewL4qNhwTK54-6>4L44N!{kWc~Xcn;Fe!0^CclkpaN
zL1Iy2ZhVm@(=En|TWo2G$tC$kD;bLT7#J9Sh3IGG=ceixm!zg8=IN)D=B0uO-Q3LN
zqWt2LqSEA&(xO!T)QW=CqRiaXyb}HL{QMN%q|BVe%G4q-CqA{bD8C>zF)uzTH8nXy
zub}dlU}j!QW^!Ulevw;gUUErheqJQV3!uPYVH99wVPs({5&(q;DEMy)r>Ewn7A5Ay
zC+1}q=a&@a7i1=by^;>`1p@<0sDMHU6b+zI`3$OsQW#4Zq!~&WMHot$L>N+-Qkc@2
zQ<#uMnL(mpF=SB|5r%Z;bmm}&mCVqP1DUMJ3Jb7|TihwBIVFkl`MIg-iMQAja|$vN
z<9%R^5Evt9C1Vi}$TOe-1&3-8KLZ29ErEif)Z&t&#LT?-f{euCRKrpR28Pd|ICWqw
zl48K)t0Ezgal#A?44SM(A|SRXh>!pgVjw~sL`Z@N0(Kg8ffPxD2vGVd0;P>2P~DCZ
zqPz?Y44@D#wqRgjC}CK@Si-b`xrAW>OA2EN>p~_dkFA7h0ecBY4MQ{ILdJ!RC7fAY
z3z-%&Eo5B4y^tZ8L6gZ3Y=tJzE!NbW#Nv`${AKx>=_T<c`SFmTxW$xLaEm!5zvPxw
zQGQ7xC~3s!CYBUsR>UXg=N9B-Cgvrl-r|NN*^-RZlEhp55S<`1z=qvoO)M%(tON()
zEk?gvj3Kw!O7e3ui%W`#2@>P!pri<j1`Ea_4Nw@{++xX1EXlaVmRwR<kXn3;DZls@
zb8%(yE!Mo!+=9v?c~EQ$f&G)2oLN#CUy`4nlayakBnt{xZV&;Aq+0?ISA&yFd|FOo
z$t@vdk&?us^wg3fC6IbheBNROt1MCh882P{%0opZnW@F`$@!^iX_?8Hsd**Ex0qey
zgFr+Gi14|^;vOFq0wR5az+?!>k?|oxV9EzdK?Qxlf<X{okWUdf>3|avIQ@V!;4QWQ
zP|_^E#SKZK@t_a@*<=VJj6j4j*c3w}2nFF8K`BG9pfQv(0=dW#Od1=4*ha>;co4A~
zAACzPH?btOC^Io908}KUmZTOHL-PhCHX*Jsg17>wVxuBZT828=7_7KR5frMbAVL8|
zsDZ*1rvb*`c!y+Fv<$>ylbfGXnv-e=DhrCIFfcH1FmNz&FhU>;52)y3VPs+CU}0fm
z;bCK9V?=}BMEFFQP$3T!4<ipF3kz5uI|~O3hX5l4Dly72$}tLoA=_^jHV#G>c8H7s
SBZeF^3p)>!7$YAe4-){{*V{P&
literal 0
HcmV?d00001
diff --git a/experiment/wood-bilayer/readme.txt b/experiment/wood-bilayer/readme.txt
new file mode 100644
index 00000000..e69de29b
diff --git a/experiment/wood-bilayer/results/0/BMatrix.txt b/experiment/wood-bilayer/results/0/BMatrix.txt
new file mode 100644
index 00000000..84844fde
--- /dev/null
+++ b/experiment/wood-bilayer/results/0/BMatrix.txt
@@ -0,0 +1,3 @@
+1 1 1.22028443964317246
+1 2 -0.2845146939877648
+1 3 -1.41110998730663744e-33
diff --git a/experiment/wood-bilayer/results/0/QMatrix.txt b/experiment/wood-bilayer/results/0/QMatrix.txt
new file mode 100644
index 00000000..2969e15f
--- /dev/null
+++ b/experiment/wood-bilayer/results/0/QMatrix.txt
@@ -0,0 +1,9 @@
+1 1 304.52882446568799
+1 2 19.1058810430896244
+1 3 -1.65436122510605535e-24
+2 1 19.1058932055586403
+2 2 524.281062366142351
+2 3 1.18131731230229265e-23
+3 1 7.66396164928418806e-33
+3 2 1.21058723308231638e-32
+3 3 88.8021141416043207
diff --git a/experiment/wood-bilayer/results/0/output.txt b/experiment/wood-bilayer/results/0/output.txt
new file mode 100644
index 00000000..63688560
--- /dev/null
+++ b/experiment/wood-bilayer/results/0/output.txt
@@ -0,0 +1,49 @@
+Number of Grid-Elements in each direction: [16,16,16]
+Solver-type used: GMRES-Solver
+---------- OUTPUT ----------
+ --------------------
+Corrector-Matrix M_1:
+-0.233954 -1.54165e-27 0
+-1.54165e-27 0.0110659 0
+0 0 0
+
+ --------------------
+Corrector-Matrix M_2:
+-0.0025163 1.89724e-27 0
+1.89724e-27 0.11024 0
+0 0 0
+
+ --------------------
+Corrector-Matrix M_3:
+0 0.0227636 0
+0.0227636 0 0
+0 0 0
+
+ --------------------
+--- Effective moduli ---
+Qeff_:
+304.529 19.1059 -1.65436e-24
+19.1059 524.281 1.18132e-23
+7.66396e-33 1.21059e-32 88.8021
+
+------------------------
+--- Prestrain Output ---
+Bhat_: 366.176 -125.851 -1.19402e-31
+Beff_: 1.22028 -0.284515 -1.41111e-33 (Effective Prestrain)
+------------------------
+size of FiniteElementBasis: 13056
+q1=304.529
+q2=524.281
+q3=88.8021
+q12=19.1059
+q23=1.18132e-23
+q_onetwo=19.105881
+b1=1.220284
+b2=-0.284515
+b3=-0.000000
+mu_gamma=88.802114
+---------------------------------------------------------------------------------------------------------------------------------------
+ Levels | q1 | q2 | q3 | q12 | q23 | b1 | b2 | b3 |
+---------------------------------------------------------------------------------------------------------------------------------------
+ 4 & 3.04529e+02 & 5.24281e+02 & 8.88021e+01 & 1.91059e+01 & 1.18132e-23 & 1.22028e+00 & -2.84515e-01 & -1.41111e-33 &
+---------------------------------------------------------------------------------------------------------------------------------------
diff --git a/experiment/wood-bilayer/results/0/parameter.txt b/experiment/wood-bilayer/results/0/parameter.txt
new file mode 100644
index 00000000..ee7ee08f
--- /dev/null
+++ b/experiment/wood-bilayer/results/0/parameter.txt
@@ -0,0 +1,4 @@
+r = 0.22
+h = 0.0053
+omega_flat = 17.17547062
+omega_target = 14.72680026
diff --git a/experiment/wood-bilayer/results/1/BMatrix.txt b/experiment/wood-bilayer/results/1/BMatrix.txt
new file mode 100644
index 00000000..2b051e95
--- /dev/null
+++ b/experiment/wood-bilayer/results/1/BMatrix.txt
@@ -0,0 +1,3 @@
+1 1 1.7573834621025497
+1 2 -0.42040684090086522
+1 3 -1.93627631123640932e-32
diff --git a/experiment/wood-bilayer/results/1/QMatrix.txt b/experiment/wood-bilayer/results/1/QMatrix.txt
new file mode 100644
index 00000000..d6f7407d
--- /dev/null
+++ b/experiment/wood-bilayer/results/1/QMatrix.txt
@@ -0,0 +1,9 @@
+1 1 315.078030302231809
+1 2 20.6661963902054353
+1 3 -5.16987882845642297e-24
+2 1 20.6662089007290781
+2 2 537.77534277449513
+2 3 -1.27954501004296468e-24
+3 1 9.31699229364431294e-32
+3 2 -1.2996446844110351e-32
+3 3 90.3065313535569913
diff --git a/experiment/wood-bilayer/results/1/output.txt b/experiment/wood-bilayer/results/1/output.txt
new file mode 100644
index 00000000..d135fa65
--- /dev/null
+++ b/experiment/wood-bilayer/results/1/output.txt
@@ -0,0 +1,49 @@
+Number of Grid-Elements in each direction: [16,16,16]
+Solver-type used: GMRES-Solver
+---------- OUTPUT ----------
+ --------------------
+Corrector-Matrix M_1:
+-0.232471 2.20455e-26 0
+2.20455e-26 0.0115767 0
+0 0 0
+
+ --------------------
+Corrector-Matrix M_2:
+-0.0026877 1.24089e-27 0
+1.24089e-27 0.110113 0
+0 0 0
+
+ --------------------
+Corrector-Matrix M_3:
+0 0.0222889 0
+0.0222889 0 0
+0 0 0
+
+ --------------------
+--- Effective moduli ---
+Qeff_:
+315.078 20.6662 -5.16988e-24
+20.6662 537.775 -1.27955e-24
+9.31699e-32 -1.29964e-32 90.3065
+
+------------------------
+--- Prestrain Output ---
+Bhat_: 545.025 -189.766 -1.57938e-30
+Beff_: 1.75738 -0.420407 -1.93628e-32 (Effective Prestrain)
+------------------------
+size of FiniteElementBasis: 13056
+q1=315.078
+q2=537.775
+q3=90.3065
+q12=20.6662
+q23=-1.27955e-24
+q_onetwo=20.666196
+b1=1.757383
+b2=-0.420407
+b3=-0.000000
+mu_gamma=90.306531
+---------------------------------------------------------------------------------------------------------------------------------------
+ Levels | q1 | q2 | q3 | q12 | q23 | b1 | b2 | b3 |
+---------------------------------------------------------------------------------------------------------------------------------------
+ 4 & 3.15078e+02 & 5.37775e+02 & 9.03065e+01 & 2.06662e+01 & -1.27955e-24 & 1.75738e+00 & -4.20407e-01 & -1.93628e-32 &
+---------------------------------------------------------------------------------------------------------------------------------------
diff --git a/experiment/wood-bilayer/results/1/parameter.txt b/experiment/wood-bilayer/results/1/parameter.txt
new file mode 100644
index 00000000..23cacfff
--- /dev/null
+++ b/experiment/wood-bilayer/results/1/parameter.txt
@@ -0,0 +1,4 @@
+r = 0.22
+h = 0.0053
+omega_flat = 17.17547062
+omega_target = 13.64338887
diff --git a/experiment/wood-bilayer/results/2/BMatrix.txt b/experiment/wood-bilayer/results/2/BMatrix.txt
new file mode 100644
index 00000000..06fe4580
--- /dev/null
+++ b/experiment/wood-bilayer/results/2/BMatrix.txt
@@ -0,0 +1,3 @@
+1 1 2.36523001516716569
+1 2 -0.581569592422200854
+1 3 2.33065121544418776e-34
diff --git a/experiment/wood-bilayer/results/2/QMatrix.txt b/experiment/wood-bilayer/results/2/QMatrix.txt
new file mode 100644
index 00000000..b04d2a62
--- /dev/null
+++ b/experiment/wood-bilayer/results/2/QMatrix.txt
@@ -0,0 +1,9 @@
+1 1 327.048205759343205
+1 2 22.5165789852819032
+1 3 -3.72877510502419507e-24
+2 1 22.5165911472273521
+2 2 553.147634470043045
+2 3 1.34416849539866997e-24
+3 1 9.61672195943999387e-33
+3 2 6.28741120991523548e-33
+3 3 92.0137173817475684
diff --git a/experiment/wood-bilayer/results/2/output.txt b/experiment/wood-bilayer/results/2/output.txt
new file mode 100644
index 00000000..27357ee2
--- /dev/null
+++ b/experiment/wood-bilayer/results/2/output.txt
@@ -0,0 +1,49 @@
+Number of Grid-Elements in each direction: [16,16,16]
+Solver-type used: GMRES-Solver
+---------- OUTPUT ----------
+ --------------------
+Corrector-Matrix M_1:
+-0.230885 2.19762e-26 0
+2.19762e-26 0.0121584 0
+0 0 0
+
+ --------------------
+Corrector-Matrix M_2:
+-0.00288395 3.79325e-27 0
+3.79325e-27 0.109978 0
+0 0 0
+
+ --------------------
+Corrector-Matrix M_3:
+0 0.0217664 0
+0.0217664 0 0
+0 0 0
+
+ --------------------
+--- Effective moduli ---
+Qeff_:
+327.048 22.5166 -3.72878e-24
+22.5166 553.148 1.34417e-24
+9.61672e-33 6.28741e-33 92.0137
+
+------------------------
+--- Prestrain Output ---
+Bhat_: 760.449 -268.437 4.05344e-32
+Beff_: 2.36523 -0.58157 2.33065e-34 (Effective Prestrain)
+------------------------
+size of FiniteElementBasis: 13056
+q1=327.048
+q2=553.148
+q3=92.0137
+q12=22.5166
+q23=1.34417e-24
+q_onetwo=22.516579
+b1=2.365230
+b2=-0.581570
+b3=0.000000
+mu_gamma=92.013717
+---------------------------------------------------------------------------------------------------------------------------------------
+ Levels | q1 | q2 | q3 | q12 | q23 | b1 | b2 | b3 |
+---------------------------------------------------------------------------------------------------------------------------------------
+ 4 & 3.27048e+02 & 5.53148e+02 & 9.20137e+01 & 2.25166e+01 & 1.34417e-24 & 2.36523e+00 & -5.81570e-01 & 2.33065e-34 &
+---------------------------------------------------------------------------------------------------------------------------------------
diff --git a/experiment/wood-bilayer/results/2/parameter.txt b/experiment/wood-bilayer/results/2/parameter.txt
new file mode 100644
index 00000000..1b0b796f
--- /dev/null
+++ b/experiment/wood-bilayer/results/2/parameter.txt
@@ -0,0 +1,4 @@
+r = 0.22
+h = 0.0053
+omega_flat = 17.17547062
+omega_target = 12.41305478
diff --git a/experiment/wood-bilayer/results/3/BMatrix.txt b/experiment/wood-bilayer/results/3/BMatrix.txt
new file mode 100644
index 00000000..82579fb3
--- /dev/null
+++ b/experiment/wood-bilayer/results/3/BMatrix.txt
@@ -0,0 +1,3 @@
+1 1 2.73380820761312782
+1 2 -0.68295979603524104
+1 3 -1.51829981430077995e-33
diff --git a/experiment/wood-bilayer/results/3/QMatrix.txt b/experiment/wood-bilayer/results/3/QMatrix.txt
new file mode 100644
index 00000000..ba86d737
--- /dev/null
+++ b/experiment/wood-bilayer/results/3/QMatrix.txt
@@ -0,0 +1,9 @@
+1 1 334.324056037577748
+1 2 23.6828912676120602
+1 3 4.21345124519198472e-24
+2 1 23.6828964275446054
+2 2 562.522258871768145
+2 3 4.30392412468997212e-24
+3 1 -1.16480740815231617e-32
+3 2 -8.12511691239257716e-33
+3 3 93.0513170279158714
diff --git a/experiment/wood-bilayer/results/3/output.txt b/experiment/wood-bilayer/results/3/output.txt
new file mode 100644
index 00000000..cdcdccdb
--- /dev/null
+++ b/experiment/wood-bilayer/results/3/output.txt
@@ -0,0 +1,49 @@
+Number of Grid-Elements in each direction: [16,16,16]
+Solver-type used: GMRES-Solver
+---------- OUTPUT ----------
+ --------------------
+Corrector-Matrix M_1:
+-0.229968 2.88958e-27 0
+2.88958e-27 0.012513 0
+0 0 0
+
+ --------------------
+Corrector-Matrix M_2:
+-0.00300407 -2.3938e-28 0
+-2.3938e-28 0.109899 0
+0 0 0
+
+ --------------------
+Corrector-Matrix M_3:
+0 0.0214569 0
+0.0214569 0 0
+0 0 0
+
+ --------------------
+--- Effective moduli ---
+Qeff_:
+334.324 23.6829 4.21345e-24
+23.6829 562.522 4.30392e-24
+-1.16481e-32 -8.12512e-33 93.0513
+
+------------------------
+--- Prestrain Output ---
+Bhat_: 897.803 -319.436 -1.67574e-31
+Beff_: 2.73381 -0.68296 -1.5183e-33 (Effective Prestrain)
+------------------------
+size of FiniteElementBasis: 13056
+q1=334.324
+q2=562.522
+q3=93.0513
+q12=23.6829
+q23=4.30392e-24
+q_onetwo=23.682891
+b1=2.733808
+b2=-0.682960
+b3=-0.000000
+mu_gamma=93.051317
+---------------------------------------------------------------------------------------------------------------------------------------
+ Levels | q1 | q2 | q3 | q12 | q23 | b1 | b2 | b3 |
+---------------------------------------------------------------------------------------------------------------------------------------
+ 4 & 3.34324e+02 & 5.62522e+02 & 9.30513e+01 & 2.36829e+01 & 4.30392e-24 & 2.73381e+00 & -6.82960e-01 & -1.51830e-33 &
+---------------------------------------------------------------------------------------------------------------------------------------
diff --git a/experiment/wood-bilayer/results/3/parameter.txt b/experiment/wood-bilayer/results/3/parameter.txt
new file mode 100644
index 00000000..3ab41d62
--- /dev/null
+++ b/experiment/wood-bilayer/results/3/parameter.txt
@@ -0,0 +1,4 @@
+r = 0.22
+h = 0.0053
+omega_flat = 17.17547062
+omega_target = 11.66482931
diff --git a/experiment/wood-bilayer/results/4/BMatrix.txt b/experiment/wood-bilayer/results/4/BMatrix.txt
new file mode 100644
index 00000000..a0f25ae1
--- /dev/null
+++ b/experiment/wood-bilayer/results/4/BMatrix.txt
@@ -0,0 +1,3 @@
+1 1 3.01257794449986038
+1 2 -0.761426746467108506
+1 3 -2.00301081070281749e-33
diff --git a/experiment/wood-bilayer/results/4/QMatrix.txt b/experiment/wood-bilayer/results/4/QMatrix.txt
new file mode 100644
index 00000000..7dbc4bb9
--- /dev/null
+++ b/experiment/wood-bilayer/results/4/QMatrix.txt
@@ -0,0 +1,9 @@
+1 1 339.836353364801084
+1 2 24.5874784431021602
+1 3 -9.71937219749807518e-24
+2 1 24.5874898806789268
+2 2 569.639985045040589
+2 3 2.29930360895599411e-23
+3 1 -2.0503720577364613e-32
+3 2 -1.25878986003714669e-32
+3 3 93.8373197975450353
diff --git a/experiment/wood-bilayer/results/4/output.txt b/experiment/wood-bilayer/results/4/output.txt
new file mode 100644
index 00000000..29680bb5
--- /dev/null
+++ b/experiment/wood-bilayer/results/4/output.txt
@@ -0,0 +1,49 @@
+Number of Grid-Elements in each direction: [16,16,16]
+Solver-type used: GMRES-Solver
+---------- OUTPUT ----------
+ --------------------
+Corrector-Matrix M_1:
+-0.229296 1.44303e-25 0
+1.44303e-25 0.012782 0
+0 0 0
+
+ --------------------
+Corrector-Matrix M_2:
+-0.00309544 1.50887e-26 0
+1.50887e-26 0.109842 0
+0 0 0
+
+ --------------------
+Corrector-Matrix M_3:
+0 0.0212263 0
+0.0212263 0 0
+0 0 0
+
+ --------------------
+--- Effective moduli ---
+Qeff_:
+339.836 24.5875 -9.71937e-24
+24.5875 569.64 2.2993e-23
+-2.05037e-32 -1.25879e-32 93.8373
+
+------------------------
+--- Prestrain Output ---
+Bhat_: 1005.06 -359.667 -2.40141e-31
+Beff_: 3.01258 -0.761427 -2.00301e-33 (Effective Prestrain)
+------------------------
+size of FiniteElementBasis: 13056
+q1=339.836
+q2=569.64
+q3=93.8373
+q12=24.5875
+q23=2.2993e-23
+q_onetwo=24.587478
+b1=3.012578
+b2=-0.761427
+b3=-0.000000
+mu_gamma=93.837320
+---------------------------------------------------------------------------------------------------------------------------------------
+ Levels | q1 | q2 | q3 | q12 | q23 | b1 | b2 | b3 |
+---------------------------------------------------------------------------------------------------------------------------------------
+ 4 & 3.39836e+02 & 5.69640e+02 & 9.38373e+01 & 2.45875e+01 & 2.29930e-23 & 3.01258e+00 & -7.61427e-01 & -2.00301e-33 &
+---------------------------------------------------------------------------------------------------------------------------------------
diff --git a/experiment/wood-bilayer/results/4/parameter.txt b/experiment/wood-bilayer/results/4/parameter.txt
new file mode 100644
index 00000000..391a2bf4
--- /dev/null
+++ b/experiment/wood-bilayer/results/4/parameter.txt
@@ -0,0 +1,4 @@
+r = 0.22
+h = 0.0053
+omega_flat = 17.17547062
+omega_target = 11.09781471
diff --git a/experiment/wood-bilayer/results/5/BMatrix.txt b/experiment/wood-bilayer/results/5/BMatrix.txt
new file mode 100644
index 00000000..8e6000c5
--- /dev/null
+++ b/experiment/wood-bilayer/results/5/BMatrix.txt
@@ -0,0 +1,3 @@
+1 1 3.82701535427334427
+1 2 -0.999166807715384686
+1 3 1.36223980815737315e-34
diff --git a/experiment/wood-bilayer/results/5/QMatrix.txt b/experiment/wood-bilayer/results/5/QMatrix.txt
new file mode 100644
index 00000000..c14621ce
--- /dev/null
+++ b/experiment/wood-bilayer/results/5/QMatrix.txt
@@ -0,0 +1,9 @@
+1 1 355.98933444216442
+1 2 27.3425221293822105
+1 3 -1.18907213054497728e-24
+2 1 27.3425338419695017
+2 2 590.572959316512765
+2 3 1.55096364853692689e-25
+3 1 -1.32587679531366614e-33
+3 2 -1.39370368771301563e-33
+3 3 96.1397819755019469
diff --git a/experiment/wood-bilayer/results/5/output.txt b/experiment/wood-bilayer/results/5/output.txt
new file mode 100644
index 00000000..f9656ea8
--- /dev/null
+++ b/experiment/wood-bilayer/results/5/output.txt
@@ -0,0 +1,49 @@
+Number of Grid-Elements in each direction: [16,16,16]
+Solver-type used: GMRES-Solver
+---------- OUTPUT ----------
+ --------------------
+Corrector-Matrix M_1:
+-0.227431 4.59477e-27 0
+4.59477e-27 0.0135725 0
+0 0 0
+
+ --------------------
+Corrector-Matrix M_2:
+-0.00336484 3.60248e-27 0
+3.60248e-27 0.109684 0
+0 0 0
+
+ --------------------
+Corrector-Matrix M_3:
+0 0.0205697 0
+0.0205697 0 0
+0 0 0
+
+ --------------------
+--- Effective moduli ---
+Qeff_:
+355.989 27.3425 -1.18907e-24
+27.3425 590.573 1.55096e-25
+-1.32588e-33 -1.3937e-33 96.1398
+
+------------------------
+--- Prestrain Output ---
+Bhat_: 1335.06 -485.441 9.41494e-33
+Beff_: 3.82702 -0.999167 1.36224e-34 (Effective Prestrain)
+------------------------
+size of FiniteElementBasis: 13056
+q1=355.989
+q2=590.573
+q3=96.1398
+q12=27.3425
+q23=1.55096e-25
+q_onetwo=27.342522
+b1=3.827015
+b2=-0.999167
+b3=0.000000
+mu_gamma=96.139782
+---------------------------------------------------------------------------------------------------------------------------------------
+ Levels | q1 | q2 | q3 | q12 | q23 | b1 | b2 | b3 |
+---------------------------------------------------------------------------------------------------------------------------------------
+ 4 & 3.55989e+02 & 5.90573e+02 & 9.61398e+01 & 2.73425e+01 & 1.55096e-25 & 3.82702e+00 & -9.99167e-01 & 1.36224e-34 &
+---------------------------------------------------------------------------------------------------------------------------------------
diff --git a/experiment/wood-bilayer/results/5/parameter.txt b/experiment/wood-bilayer/results/5/parameter.txt
new file mode 100644
index 00000000..861b07bc
--- /dev/null
+++ b/experiment/wood-bilayer/results/5/parameter.txt
@@ -0,0 +1,4 @@
+r = 0.22
+h = 0.0053
+omega_flat = 17.17547062
+omega_target = 9.435795985
diff --git a/experiment/wood-bilayer/results/6/BMatrix.txt b/experiment/wood-bilayer/results/6/BMatrix.txt
new file mode 100644
index 00000000..72e1ddee
--- /dev/null
+++ b/experiment/wood-bilayer/results/6/BMatrix.txt
@@ -0,0 +1,3 @@
+1 1 4.05964843125063535
+1 2 -1.06932937107897685
+1 3 5.63758516062832998e-32
diff --git a/experiment/wood-bilayer/results/6/QMatrix.txt b/experiment/wood-bilayer/results/6/QMatrix.txt
new file mode 100644
index 00000000..b1b97742
--- /dev/null
+++ b/experiment/wood-bilayer/results/6/QMatrix.txt
@@ -0,0 +1,9 @@
+1 1 360.61715085991392
+1 2 28.1605370711061482
+1 3 -3.51551760335036762e-24
+2 1 28.1605407781263679
+2 2 596.590840822605628
+2 3 9.85508151674505628e-24
+3 1 -1.17988935076666896e-31
+3 2 -8.86266063963928227e-32
+3 3 96.7991406490075548
diff --git a/experiment/wood-bilayer/results/6/output.txt b/experiment/wood-bilayer/results/6/output.txt
new file mode 100644
index 00000000..96a2e7ea
--- /dev/null
+++ b/experiment/wood-bilayer/results/6/output.txt
@@ -0,0 +1,49 @@
+Number of Grid-Elements in each direction: [16,16,16]
+Solver-type used: GMRES-Solver
+---------- OUTPUT ----------
+ --------------------
+Corrector-Matrix M_1:
+-0.226924 7.72095e-26 0
+7.72095e-26 0.0137995 0
+0 0 0
+
+ --------------------
+Corrector-Matrix M_2:
+-0.00344244 -2.16662e-27 0
+-2.16662e-27 0.109641 0
+0 0 0
+
+ --------------------
+Corrector-Matrix M_3:
+0 0.0203866 0
+0.0203866 0 0
+0 0 0
+
+ --------------------
+--- Effective moduli ---
+Qeff_:
+360.617 28.1605 -3.51552e-24
+28.1605 596.591 9.85508e-24
+-1.17989e-31 -8.86266e-32 96.7991
+
+------------------------
+--- Prestrain Output ---
+Bhat_: 1433.87 -523.63 5.07291e-30
+Beff_: 4.05965 -1.06933 5.63759e-32 (Effective Prestrain)
+------------------------
+size of FiniteElementBasis: 13056
+q1=360.617
+q2=596.591
+q3=96.7991
+q12=28.1605
+q23=9.85508e-24
+q_onetwo=28.160537
+b1=4.059648
+b2=-1.069329
+b3=0.000000
+mu_gamma=96.799141
+---------------------------------------------------------------------------------------------------------------------------------------
+ Levels | q1 | q2 | q3 | q12 | q23 | b1 | b2 | b3 |
+---------------------------------------------------------------------------------------------------------------------------------------
+ 4 & 3.60617e+02 & 5.96591e+02 & 9.67991e+01 & 2.81605e+01 & 9.85508e-24 & 4.05965e+00 & -1.06933e+00 & 5.63759e-32 &
+---------------------------------------------------------------------------------------------------------------------------------------
diff --git a/experiment/wood-bilayer/results/6/parameter.txt b/experiment/wood-bilayer/results/6/parameter.txt
new file mode 100644
index 00000000..4dc71462
--- /dev/null
+++ b/experiment/wood-bilayer/results/6/parameter.txt
@@ -0,0 +1,4 @@
+r = 0.22
+h = 0.0053
+omega_flat = 17.17547062
+omega_target = 8.959564147
diff --git a/experiment/wood-bilayer/results/auswertung.py b/experiment/wood-bilayer/results/auswertung.py
index 12e9fc7a..204c9a30 100644
--- a/experiment/wood-bilayer/results/auswertung.py
+++ b/experiment/wood-bilayer/results/auswertung.py
@@ -1,12 +1,23 @@
import matplotlib.pyplot as plt
+from matplotlib.ticker import StrMethodFormatter
+
omega=[
14.72680026, 13.64338887, 12.41305478, 11.66482931, 11.09781471, 9.435795985, 8.959564147]
kappa_sim=[1.20240481, 1.73346693, 2.3246493, 2.68537074, 2.95591182, 3.74749499,
3.97795591]
kappa_exp=[
1.058078122, 1.544624544, 2.317033799, 2.686043143, 2.967694189, 3.913528418, 4.262750825]
-plt.scatter(omega,kappa_sim, label="simulation")
-plt.scatter(omega, kappa_exp, marker='x', label="experiment")
-plt.xlabel="target humidity"
-plt.ylabel="curvature"
+
+# plotting a line plot after changing it's width and height
+f = plt.figure()
+f.set_figwidth(6)
+f.set_figheight(3)
+plt.gca().yaxis.set_major_formatter(StrMethodFormatter('{x:,.1f}')) # 2 decimal places
+plt.gca().xaxis.set_major_formatter(StrMethodFormatter('{x:,.1f}')) # 2 decimal places
+plt.scatter(omega, kappa_exp, marker='o', label="Experiment")
+plt.scatter(omega, kappa_sim, marker='x', label="Simulation")
+plt.xticks(omega)
+plt.ylabel(r"curvature $\kappa$ (1/mm)")
+plt.xlabel(r"wood moisture content $\omega$ (%)")
+plt.legend()
plt.show()
diff --git a/experiment/wood-bilayer/results/kappa_simulation.txt b/experiment/wood-bilayer/results/kappa_simulation.txt
new file mode 100644
index 00000000..ff67e05b
--- /dev/null
+++ b/experiment/wood-bilayer/results/kappa_simulation.txt
@@ -0,0 +1,2 @@
+[1.20240481 1.73346693 2.3246493 2.68537074 2.95591182 3.74749499
+ 3.97795591]
\ No newline at end of file
diff --git a/experiment/wood-bilayer/wood_european_beech.py b/experiment/wood-bilayer/wood_european_beech.py
index 456cb0bc..e30fd721 100644
--- a/experiment/wood-bilayer/wood_european_beech.py
+++ b/experiment/wood-bilayer/wood_european_beech.py
@@ -15,16 +15,20 @@ import elasticity_toolbox as elast
#--- define indicator function
# Parameters of the model
-# -- ratio between upper layer / thickness, thickness [m]
+# -- (thickness upper layer) / (thickness)
param_r = 0.22
+# -- thickness [meter]
param_h = 0.0053
-# -- moisture content in the flat state and target state [%]
+# -- moisture content in the flat state [%]
param_omega_flat = 17.17547062
+# -- moisture content in the target state [%]
param_omega_target = 8.959564147
+# -- Drehwinkel (nicht implementiert)
param_theta = 0
#
#
#
+# -- increment of the moisture content
delta_omega=param_omega_target-param_omega_flat
# moisture content for material law
omega=param_omega_target
@@ -41,40 +45,31 @@ omega=param_omega_target
properties_coefficients=np.array([
# [b_0, b_1]
-[2565.6,-59.7], # E_R [MPa]
-[885.4, -23.4], # E_T [MPa]
-[17136.7,-282.4], # E_L [MPa]
-[667.8, -15.19], # G_RT [MPa]
-[1482, -15.26], # G_RL [MPa]
-[1100, -17.72], # G_TL [MPa]
-[0.2933, -0.001012], # nu_TR [1]
-[0.383, -0.008722], # nu_LR [1]
-[0.3368, -0.009071] # nu_LT [1]
-])
-# Compute actual properties
-#E_R=1500
+ [2565.6,-59.7], # E_R [MPa]
+ [885.4, -23.4], # E_T [MPa]
+ [17136.7,-282.4], # E_L [MPa]
+ [667.8, -15.19], # G_RT [MPa]
+ [1482, -15.26], # G_RL [MPa]
+ [1100, -17.72], # G_TL [MPa]
+ [0.2933, -0.001012], # nu_TR [1]
+ [0.383, -0.008722], # nu_LR [1]
+ [0.3368, -0.009071] # nu_LT [1]
+ ])
+# Compute actual material properties
E_R = properties_coefficients[0,0]+properties_coefficients[0,1]*omega
-# E_T=450
E_T = properties_coefficients[1,0]+properties_coefficients[1,1]*omega
-#E_L=12000
E_L = properties_coefficients[2,0]+properties_coefficients[2,1]*omega
-# G_RT=400
G_RT = properties_coefficients[3,0]+properties_coefficients[3,1]*omega
-# G_RL=1200
G_LR = properties_coefficients[4,0]+properties_coefficients[4,1]*omega
-# G_TL=800
G_LT = properties_coefficients[5,0]+properties_coefficients[5,1]*omega
-# nu_TR=0.28
nu_TR = properties_coefficients[6,0]+properties_coefficients[6,1]*omega
-# nu_LR=0.2125
nu_LR = properties_coefficients[7,0]+properties_coefficients[7,1]*omega
-# nu_LT=0.175
nu_LT = properties_coefficients[8,0]+properties_coefficients[8,1]*omega
-#
+# Compute the remaining Poisson ratios
nu_TL=nu_LT*E_T/E_L
nu_RT=nu_TR*E_R/E_T
nu_RL=nu_LR*E_R/E_L
-
+#
# --- differential swelling strain
# --- relation to swelling strain eps: eps=alpha* delta_omega with delta_omega = change of water content in %
alpha_L=0.00011 # PLOS paper
@@ -114,7 +109,9 @@ compliance_S=np.array([[1/E_1, -nu_21/E_2, -nu_31/E_3, 0.0,
[0.0, 0.0, 0.0, 0.0, 1/G_31, 0.0],
[0.0, 0.0, 0.0, 0.0, 0.0, 1/G_12]]);
materialParameters_phase1 = compliance_S
+
def prestrain_phase1(x):
+ # hB=delta_omega * alpha with delta_omega increment of moisture content and alpha swelling factor.
return [[1/param_h*delta_omega*alpha_L, 0, 0], [0,1/param_h*delta_omega*alpha_T,0], [0,0,1/param_h*delta_omega*alpha_R]]
# --- PHASE 2
@@ -133,8 +130,6 @@ compliance_S=np.array([[1/E_1, -nu_21/E_2, -nu_31/E_3, 0.0,
[0.0, 0.0, 0.0, 0.0, 1/G_31, 0.0],
[0.0, 0.0, 0.0, 0.0, 0.0, 1/G_12]]);
materialParameters_phase2 = compliance_S
-#phase2_type="orthotropic"
-#materialParameters_phase2 = [E_R, E_L, E_T,G_RL, G_TL, G_RT,nu_LR, nu_TR,nu_LT]
def prestrain_phase2(x):
return [[1/param_h*delta_omega*alpha_R, 0, 0], [0,1/param_h*delta_omega*alpha_L,0], [0,0,1/param_h*delta_omega*alpha_T]]
@@ -149,4 +144,4 @@ materialParameters_phase3 = np.dot(np.dot(N,materialParameters_phase1),N.T)
materialParameters_phase3 = 0.5*(materialParameters_phase3.T+materialParameters_phase3)
# rotation of strain
def prestrain_phase3(x):
- return elast.voigt_to_strain(np.dot(elast.rotation_matrix_compliance(2,param_theta),np.dot(elast.strain_to_voigt(np.array(prestrain_phase1(x))),N.T))).tolist()
+ return elast.voigt_to_strain(np.dot(elast.rotation_matrix_compliance(2,param_theta),np.dot(elast.strain_to_voigt(np.array(prestrain_phase1(x))),N.T))).tolist()
\ No newline at end of file
diff --git a/experiment/wood-bilayer/wood_test.py b/experiment/wood-bilayer/wood_test.py
index 56a48fd0..b8194e49 100644
--- a/experiment/wood-bilayer/wood_test.py
+++ b/experiment/wood-bilayer/wood_test.py
@@ -113,6 +113,13 @@ ParameterSet.numLevels=4
# ----- Define Parameters for Material Function --------------------
# [r, h, omega_flat, omega_target, theta, experimental_kappa]
+# r = (thickness upper layer)/(thickness)
+# h = thickness [meter]
+# omega_flat = moisture content in the flat state before drying [%]
+# omega_target = moisture content in the target state [%]
+# theta = rotation angle (not implemented and used)
+# experimental_kappa = curvature measure in experiment
+#
materialFunctionParameter=[
[0.22, 0.0053, 17.17547062, 14.72680026, 0, 1.058078122],
[0.22, 0.0053, 17.17547062, 13.64338887, 0, 1.544624544],
@@ -122,11 +129,6 @@ materialFunctionParameter=[
[0.22, 0.0053, 17.17547062, 9.435795985, 0, 3.913528418],
[0.22, 0.0053, 17.17547062, 8.959564147, 0, 4.262750825]
]
-# materialFunctionParameter=[
-# [0.22, 0.0053, 17.17547062, 8.959564147, 0, 4.262750825],
-# [0.22, 0.0053, 17.17547062, 8.959564147, np.pi/2, 4.262750825]
-# ]
-
# ------ Loops through Parameters for Material Function -----------
for i in range(0,np.shape(materialFunctionParameter)[0]):
print("------------------")
--
GitLab