Guccione

cardmech/conf/elasticity_benchmarks/ellipsoid.conf

@@ -57,7 +57,7 @@ StretchFactorNormal = 4.0
 ######################
 # Material Parameters#
 ######################
-ActiveMaterial = Linear
+#ActiveMaterial = Linear
 #ActiveMaterial = NeoHooke
 #ActiveMaterial = Bonet
 ActiveMaterial = Holzapfel
@@ -71,7 +71,7 @@ Density = 0.00
 # Volumetric Penalty
 #------------------------
 Incompressible = false
-QuasiCompressiblePenalty = None
+QuasiCompressiblePenalty = Ciarlet
 VolumetricPenalty = 20
 PermeabilityPenalty = 20
 
@@ -132,7 +132,7 @@ PlotCellPotential = 0
 
 # === Mechanics Solver ============================================
 MechSolver = gmres
-MechPreconditioner = GaussSeidel;
+MechPreconditioner = PS;
 MechEpsilon = 1e-8;
 MechReduction = 1e-12;
 MechSteps = 10000;
@@ -164,7 +164,7 @@ postsmoothing = 3;
 
 
 DGSign=-1.0
-DGPenalty=7.0
+DGPenalty=3000000.0
 Overlap = dG1
 #Overlap = NoOverlap
 Overlap_Distribution = 1

cardmech/conf/m++.conf

@@ -16,7 +16,7 @@ GeoPath = ../geo/
 # === Cardiac Mechanics ================
 #loadconf = coupled.conf
 #loadconf = klotz.conf
-loadconf = elasticity_benchmarks/beam.conf
+#loadconf = elasticity_benchmarks/beam.conf
 loadconf = elasticity_benchmarks/ellipsoid.conf
 loadconf = elasticity.conf
 # ======================================

cardmech/src/elasticity/materials/GuccioneMaterial.hpp

@@ -23,7 +23,7 @@ class GuccioneMaterial : public HyperelasticMaterial {
 
   std::pair<size_t, size_t> GetMaxParameterRange() const {
       // start, length
-      return {0, 4};
+      return {0, 1};
   }
 
   Tensor fiberWeights(const Tensor &E) const {

data/benchmarks/CardiacBeamBonetP1DG90

+start program on 1 procs at Thu Mar  3 00:56:05 2022
+Running on: laura-MacBookPro
+   Config Info:
+     PressurePlotVerbose: ..................... -1
+     PrestressPlotVerbose: .................... -1
+     CoupledSolverVTK: ........................ -1
+     NewtonDamping: ........................... 0.75
+     DynamicSolverVTK: ........................ -1
+     PressureSolverVTK: ....................... -1
+     ElphyVerbose: ............................ -1
+     CellModelVerbose: ........................ -10
+     CoupledSolverVerbose: .................... 2
+     DynamicSolverVerbose: .................... -1
+     NonLinearVerbose: ........................ -5
+     MonodomainVerbose: ....................... -1
+     LinearImplicitVerbose: ................... -5
+     NewtonVerbose: ........................... 1
+     LinearVerbose: ........................... -5
+     BaseSolverVerbose: ....................... -1
+     ReferenceVerbose: ........................ 2
+     AssembleVerbose: ......................... 1
+     ElphySolverVerbose: ...................... -1
+     ConfigVerbose: ........................... 1
+     TimeLevel: ............................... 0
+     DebugLevel: .............................. 0
+     Overlap_Distribution: .................... 1
+     DGPenalty: ............................... 90
+     DynamicPlotVerbose: ...................... -1
+     postsmoothing: ........................... 3
+     Density: ................................. 0.00
+     PressureSolverVerbose: ................... 1
+     ActiveMaterial: .......................... Bonet
+     PrestressSteps: .......................... 1
+     PressureDepth: ........................... 0
+     InterpolateStartVector: .................. false
+     PressureIterations: ...................... 10
+     ClearData: ............................... true
+     Overlap: ................................. dG1
+     StretchFactorSheet: ...................... 1.0
+     MechEpsilon: ............................. 1e-8
+     ElphyLinearSolverVerbose: ................ -1
+     SetDisplacement: ......................... 0
+     NewmarkBeta: ............................. 0.25
+     MechModel: ............................... Static
+     PermeabilityPenalty: ..................... 2
+     h_min: ................................... 1e-6
+     MechPolynomialDegree: .................... 1
+     MechVerbose: ............................. -1
+     Model: ................................... ActiveStrainElasticity
+     MechDiscretization: ...................... DG
+     ProblemDimension: ........................ 3
+     MainMechVerbose: ......................... 1
+     CalculateReference: ...................... 1
+     MechPlot: ................................ -1
+     DGSign: .................................. -1
+     PlotCellPotential: ....................... 0
+     MultigridVerbose: ........................ 0
+     MechLevel: ............................... 0
+     GeoPath: ................................. ../geo/
+     ProblemGeometry: ......................... Tet
+     NewtonEpsilon: ........................... 1e-6
+     ElphySolverVTK: .......................... -1
+     NewmarkGamma: ............................ 0.5
+     DebugPressure: ........................... 1
+     MechProblemVerbose: ...................... 1
+     h0: ...................................... 1e-4
+     NewtonSteps: ............................. 200
+     QuasiCompressiblePenalty: ................ Ciarlet
+     MechProblem: ............................. CardiacBeam
+     ReferenceLevel: .......................... 4
+     MechMeshPart: ............................ vessel_base
+     RVPressure: .............................. 0
+     Incompressible: .......................... false
+     precision: ............................... 10
+     MechPreconditioner: ...................... GaussSeidel
+     MeshVerbose: ............................. 2
+     VolumetricPenalty: ....................... 2
+     MechSteps: ............................... 10000
+     GuccioneMat_bs: .......................... 2
+     GuccioneMat_C: ........................... 2
+     PlotVTK: ................................. -1
+     StretchFactorNormal: ..................... 4.0
+     GuccioneMat_bf: .......................... 8
+     WithPrestress: ........................... false
+     GuccioneMat_bfs: ......................... 4
+     RAPressure: .............................. 0
+     LVPressure: .............................. 0
+     NewtonReduction: ......................... 1e-12
+     LAPressure: .............................. 0.004
+     logfile: ................................. ../data/benchmarks/CardiacBeamBonetP1DG90
+     PlotVertexPotential: ..................... 0
+     PlottingSteps: ........................... 1
+     MechSolver: .............................. gmres
+     LinearMat_Mu: ............................ 2
+     PlotStress: .............................. 0
+     MechReduction: ........................... 1e-12
+     NewtonLinearizationReduction: ............ 1e-12
+     EulerVerbose: ............................ 0
+     NewtonLineSearchSteps: ................... 10
+     PressureSteps: ........................... 10
+     LinearMat_Lambda: ........................ 4
+     presmoothing: ............................ 3
+
+   Starting CardiacMechanics Suite of M++.
+   
+   =================== Running level 0================ 
+
+   Evaluation Info:
+     Time: .................................... 0
+     LV-Volume: ............................... 0
+     RV-Volume: ............................... 0
+     LA-Volume: ............................... 0
+     RA-Volume: ............................... 0
+     FrobeniusNorm: ........................... 0
+
+   Solving Pressure Step Info:
+     Step: .................................... 1 / 10
+     Pressure Scale: .......................... 0.1
+
+   Newton iterations Info:
+     Steps: ................................... 0
+     line search steps: ....................... 0
+
+   Update Pressure
+   Newton: d(0)= 0.0005163977795
+   Newton: line search 1: d(0)= 0.001574759735
+   Newton: line search 2: d(0)= 0.000598251999
+   Newton: d(1)= 0.0004701203621