add wood

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()   
 

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

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"

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]):

experiment/wood-bilayer/#readme.txt#

+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

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)

experiment/wood-bilayer/results/0/BMatrix.txt

+1 1 1.22028443964317246
+1 2 -0.2845146939877648
+1 3 -1.41110998730663744e-33

experiment/wood-bilayer/results/0/QMatrix.txt

+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

experiment/wood-bilayer/results/0/output.txt

+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 & 
+---------------------------------------------------------------------------------------------------------------------------------------

experiment/wood-bilayer/results/0/parameter.txt

+r = 0.22
+h = 0.0053
+omega_flat = 17.17547062
+omega_target = 14.72680026

experiment/wood-bilayer/results/1/BMatrix.txt

+1 1 1.7573834621025497
+1 2 -0.42040684090086522
+1 3 -1.93627631123640932e-32

experiment/wood-bilayer/results/1/QMatrix.txt

+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

experiment/wood-bilayer/results/1/output.txt

+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 & 
+---------------------------------------------------------------------------------------------------------------------------------------

experiment/wood-bilayer/results/1/parameter.txt

+r = 0.22
+h = 0.0053
+omega_flat = 17.17547062
+omega_target = 13.64338887

experiment/wood-bilayer/results/2/BMatrix.txt

+1 1 2.36523001516716569
+1 2 -0.581569592422200854
+1 3 2.33065121544418776e-34

experiment/wood-bilayer/results/2/QMatrix.txt

+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