generalized Prestress in MainElasticity

cardmech/src/elasticity/MainElasticity.cpp

@@ -105,8 +105,8 @@ void MainElasticity::initAssemble() {
   config.get("WithPrestress", prestress);
   if (prestress) {
     StaticPressureProblem prestressPressure(pressureProblem->InitialPressure());
-    LagrangeElasticity pEla(*elasticityProblem, prestressPressure, *meshes, degree,isStatic);
-    IterativePressureSolver prestressSolver(pEla, prestressPressure);
+    prestressedEla = CreateFiniteElasticityAssemble(*elasticityProblem, prestressPressure, *meshes, degree, isStatic);
+    IterativePressureSolver prestressSolver(*prestressedEla, prestressPressure);
     Vector p(0.0, *displacement);
     prestressSolver.Solve(p, prestressSteps);
     mechA->InitializePrestress(p);

cardmech/src/elasticity/MainElasticity.hpp

@@ -23,6 +23,7 @@ class MainElasticity : public MainCardMech {
   std::unique_ptr<IDiscretization> scalarDisc = nullptr;
 
   std::unique_ptr<FiniteElasticity> mechA = nullptr;
+  std::unique_ptr<FiniteElasticity> prestressedEla = nullptr;
 
   std::unique_ptr<TimeSeries> tSeries = nullptr;