From 78d8e89a7854572473c9d961e8b3f47043431125 Mon Sep 17 00:00:00 2001
From: Klaus <klaus.boehnlein@tu-dresden.de>
Date: Fri, 1 Jul 2022 12:27:56 +0200
Subject: [PATCH] Update Cellsolver.parset
---
inputs/cellsolver.parset | 16 +++++++++-------
src/Cell-Problem.cc | 7 +++++--
2 files changed, 14 insertions(+), 9 deletions(-)
diff --git a/inputs/cellsolver.parset b/inputs/cellsolver.parset
index f0962518..71dc1e54 100644
--- a/inputs/cellsolver.parset
+++ b/inputs/cellsolver.parset
@@ -39,15 +39,13 @@ numLevels = 3 3 # computes all levels from first to second entry
########################################################################################
-# --- Choose Gamma value:
-
-gamma=1.0 #(default)
+# --- Choose scale ratio gamma:
+gamma=1.0 #(default = 1.0)
#gamma=50.0
#gamma=0.01
#gamma=3.0
#gamma=0.5
#gamma = 0.05
-#gamma=2.5
#############################################
# Material / Prestrain parameters and ratios
@@ -68,8 +66,8 @@ alpha = 2.0
#--- Lame-Parameters
mu1=1.0
-#lambda1=0.0
-lambda1=1.0
+lambda1=0.0
+#lambda1=5.0
# ---volume fraction (default value = 1.0/4.0)
@@ -86,9 +84,13 @@ material_prestrain_imp= "parametrized_Laminate"
#material_prestrain_imp= "analytical_Example"
#material_prestrain_imp="isotropic_bilayer"
#material_prestrain_imp= "circle_fiber" #TEST
+#material_prestrain_imp= "matrix_material_circles"
+#material_prestrain_imp= "matrix_material_squares"
+#nF = 10 # Number of Fibers (default = 3)
+#rF = 0.1 # Fiber-Radius (default = 0.5*(width/(2.0*nF)) //half of the max-fiber-radius mrF = (width/(2.0*nF)) )
-# --- (Optional output) write Material / prestrain / Corrector functions to .vtk-Files:
+# --- (Optional output) write Material / prestrain / Corrector functions to .vtk-Files (default=false):
#write_materialFunctions = true
#write_prestrainFunctions = true # VTK norm of B ,
#write_VTK = true
diff --git a/src/Cell-Problem.cc b/src/Cell-Problem.cc
index 53b9688b..3813a62b 100644
--- a/src/Cell-Problem.cc
+++ b/src/Cell-Problem.cc
@@ -1697,7 +1697,8 @@ int main(int argc, char *argv[])
{
using VTKGridType = YaspGrid<dim, EquidistantOffsetCoordinates<double, dim> >;
// VTKGridType grid_VTK({-1.0/2.0, -1.0/2.0, -1.0/2.0},{1.0/2.0, 1.0/2.0, 1.0/2.0},{80,80,80});
- VTKGridType grid_VTK({-1.0/2.0, -1.0/2.0, -1.0/2.0},{1.0/2.0, 1.0/2.0, 1.0/2.0},nElements);
+ VTKGridType grid_VTK({-1.0/2.0, -1.0/2.0, -1.0/2.0},{1.0/2.0, 1.0/2.0, 1.0/2.0},{40,40,40});
+// VTKGridType grid_VTK({-1.0/2.0, -1.0/2.0, -1.0/2.0},{1.0/2.0, 1.0/2.0, 1.0/2.0},nElements);
using GridViewVTK = VTKGridType::LeafGridView;
const GridViewVTK gridView_VTK = grid_VTK.leafGridView();
@@ -1745,7 +1746,9 @@ int main(int argc, char *argv[])
if (write_prestrainFunctions)
{
using VTKGridType = YaspGrid<dim, EquidistantOffsetCoordinates<double, dim> >;
- VTKGridType grid_VTK({-1.0/2.0, -1.0/2.0, -1.0/2.0},{1.0/2.0, 1.0/2.0, 1.0/2.0},{80,80,80});
+// VTKGridType grid_VTK({-1.0/2.0, -1.0/2.0, -1.0/2.0},{1.0/2.0, 1.0/2.0, 1.0/2.0},{80,80,80});
+ VTKGridType grid_VTK({-1.0/2.0, -1.0/2.0, -1.0/2.0},{1.0/2.0, 1.0/2.0, 1.0/2.0},{40,40,40});
+// VTKGridType grid_VTK({-1.0/2.0, -1.0/2.0, -1.0/2.0},{1.0/2.0, 1.0/2.0, 1.0/2.0},2*nElements);
using GridViewVTK = VTKGridType::LeafGridView;
const GridViewVTK gridView_VTK = grid_VTK.leafGridView();
--
GitLab