Commit ffe6c6b7 authored by thomaskummer's avatar thomaskummer
Browse files

euler lifev update

parent fa10879d
......@@ -625,15 +625,15 @@ int main (int argc, char** argv)
}
// Update pressure b.c.
modifyPressureBC(preloadPressure(bcValues, i, preloadSteps));
//modifyNaturalPatchBC(i);
//modifyPressureBC(preloadPressure(bcValues, i, preloadSteps));
modifyNaturalPatchBC(i);
// Solve mechanics
solver.bcInterfacePtr() -> updatePhysicalSolverVariables();
solver.solveMechanics();
//solver.saveSolution (i-1);
//heartSolver.exporter()->postProcess(i-1);
heartSolver.exporter()->postProcess(i-1);
}
auto maxI4fValue ( solver.activationModelPtr()->I4f().maxValue() );
......
......@@ -249,7 +249,7 @@ public:
}
void setupExporter(std::string problemFolder = "./", std::string outputFileName = "heartSolution")
void setupExporter(std::string problemFolder = "./", std::string outputFileName = "humanHeartSolution")
{
m_exporter.reset (new exporter_Type());
setupExporter<mesh_Type>(*m_exporter, M_emSolver.localMeshPtr(), M_emSolver.comm(), outputFileName, problemFolder);
......@@ -259,6 +259,12 @@ public:
M_emSolver.structuralOperatorPtr()->dispFESpacePtr(),
M_emSolver.structuralOperatorPtr()->displacementPtr(),
UInt (0) );
// m_exporter->addVariable ( ExporterData<RegionMesh<LinearTetra> >::ScalarField,
// "Von Mises Stress",
// M_emSolver.electroSolverPtr()->feSpacePtr(),
// M_wteTotal.vonMisesStressPtr(),
// UInt (0) );
}
......@@ -276,6 +282,15 @@ public:
}
// void postProcess(const Real& time)
// {
// M_wteTotal.setDisplacement ( M_EMStructuralOperatorPtr -> displacement() );
// M_wteTotal.analyzeTensionsRecoveryVonMisesStress();
//
// m_exporter->postProcess(time);
// }
exporterPtr_Type exporter()
{
return m_exporter;
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment