diff --git a/inputs/cellsolver.parset b/inputs/cellsolver.parset index f096251875be01605c9ca6628a47baa316c7eef5..71dc1e54751e09a138e9d63845203c23ac82f321 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 53b9688b190cbd8e3da43f8216d33ee7fb04c08b..3813a62be3e2aee8d2740682ed797c8e06f37174 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();