Commit 3c125fda authored by thomaskummer's avatar thomaskummer
Browse files

intrepid lifev update

parent 0b54d1f9
......@@ -887,29 +887,30 @@ int main (int argc, char** argv)
//============================================
if ( 0 == comm->MyPID() ) circulationSolver.exportSolution( circulationOutputFile );
//============================================
// Power computations
//============================================
Real leftVentPower = heartSolver.externalPower(disp, dispPre, dETFESpace, p("lv"), dt_mechanics, 454);
Real rightVentPower = heartSolver.externalPower(disp, dispPre, dETFESpace, p("rv"), dt_mechanics, 455);
Real patchPower1 = heartSolver.externalPower(disp, dispPre, dETFESpace, p("lv"), dt_mechanics, 454);
Real patchPower2 = heartSolver.externalPower(disp, dispPre, dETFESpace, p("lv"), dt_mechanics, 454);
Real dispNormInf;
disp.normInf(dispNormInf);
if ( 0 == comm->MyPID() )
{
std::cout << "\n******************************************";
std::cout << "\nLeft ventricular power is " << leftVentPower;// << " " << p("lv");
std::cout << "\nRight ventricular power is " << dispNormInf;
std::cout << "\n******************************************\n\n";
}
dispPre = disp;
}
//============================================
// Power computations
//============================================
Real leftVentPower = heartSolver.externalPower(disp, dispPre, dETFESpace, p("lv"), dt_mechanics, 454);
Real rightVentPower = heartSolver.externalPower(disp, dispPre, dETFESpace, p("rv"), dt_mechanics, 455);
Real patchPower1 = heartSolver.externalPower(disp, dispPre, dETFESpace, p("lv"), dt_mechanics, 454);
Real patchPower2 = heartSolver.externalPower(disp, dispPre, dETFESpace, p("lv"), dt_mechanics, 454);
Real dispNormInf;
disp.normInf(dispNormInf);
if ( 0 == comm->MyPID() )
{
std::cout << "\n******************************************";
std::cout << "\nLeft ventricular power is " << leftVentPower;// << " " << p("lv");
std::cout << "\nRight ventricular power is " << dispNormInf;
std::cout << "\n******************************************\n\n";
}
dispPre = disp;
//============================================
// Export FE-solution
......
......@@ -266,6 +266,16 @@ public:
M_emSolver.tensionEstimator().vonMisesStressPtr(),
UInt (0) );
m_exporter->addVariable ( ExporterData<RegionMesh<LinearTetra> >::ScalarField,
"Principal Stress",
M_emSolver.electroSolverPtr()->feSpacePtr(),
M_emSolver.tensionEstimator().principal stresses(),
UInt (0) );
// export stresses
// export global eigenvalue vector
m_exporter->addVariable ( ExporterData<RegionMesh<LinearTetra> >::VectorField,
"Fibers",
M_emSolver.structuralOperatorPtr()->dispFESpacePtr(),
......
......@@ -1007,7 +1007,6 @@ WallTensionEstimator<Mesh >::analyzeTensionsRecoveryVonMisesStress ()
*M_sigmaVonMises += stressComponent1 * stressComponent1
+ stressComponent2 * stressComponent2
+ stressComponent3 * stressComponent3;
*M_sigmaVonMises *= 6.;
//Diagonal elements
exportStressXX ( stressComponent1 );
......
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