Added debug log output and a plot variable for outputting the current…

… pressure of an ideal-gas-pressure load.

FEBioMech/FEBioMechModule.cpp

@@ -865,6 +865,7 @@ void FEBioMech::InitModule()
     REGISTER_FECORE_CLASS(FEPlotGrowthRightHencky, "growth right Hencky");
     REGISTER_FECORE_CLASS(FEPlotGrowthLeftHencky, "growth left Hencky");
     REGISTER_FECORE_CLASS(FEPlotGrowthRelativeVolume, "growth relative volume");
+    REGISTER_FECORE_CLASS(FEPlotIdealGasPressure, "ideal gas pressure");
 
 	// beam variables
 	REGISTER_FECORE_CLASS(FEPlotBeamStress      , "beam stress");

FEBioMech/FEBioMechPlot.cpp

@@ -72,6 +72,7 @@ SOFTWARE.*/
 #include <FECore/FEElement.h>
 #include <FEBioMech/FEElasticBeamDomain.h>
 #include <FEBioMech/FEElasticBeamMaterial.h>
+#include "FEIdealGasPressure.h"
 
 //=============================================================================
 //                            N O D E   D A T A
@@ -4660,3 +4661,33 @@ bool FEPlotGrowthRelativeVolume::Save(FEDomain &dom, FEDataStream& a)
     return true;
 }
 
+
+bool FEPlotIdealGasPressure::Init()
+{
+	FEModel* fem = GetFEModel();
+	if (fem == nullptr) return false;
+
+	for (int i = 0; i < fem->ModelLoads(); ++i)
+	{
+		m_load = dynamic_cast<FEIdealGasPressure*>(fem->ModelLoad(i));
+		if (m_load) return true;
+	}
+	return (m_load != nullptr);
+}
+
+bool FEPlotIdealGasPressure::Save(FESurface& surf, FEDataStream& a)
+{
+	if (m_binit == false)
+	{
+		if (!Init()) return false;
+		m_binit = true;
+	}
+	if (m_load == nullptr) return false;
+
+	if (&m_load->GetSurface() == &surf)
+	{
+		a << m_load->GetCurrentPressure();
+		return true;
+	}
+	else return false;
+}

FEBioMech/FEBioMechPlot.h

@@ -1415,3 +1415,17 @@ class FEPlotBeamCurvature : public FEPlotDomainData
 	FEPlotBeamCurvature(FEModel* pfem) : FEPlotDomainData(pfem, PLT_VEC3F, FMT_ITEM) {}
 	bool Save(FEDomain& dom, FEDataStream& a);
 };
+
+// plot the current pressure value of the FEIdealGasPressure load
+class FEIdealGasPressure;
+class FEPlotIdealGasPressure : public FEPlotSurfaceData
+{
+public:
+	FEPlotIdealGasPressure(FEModel* fem) : FEPlotSurfaceData(fem, PLT_FLOAT, FMT_REGION) {}
+	bool Save(FESurface& surf, FEDataStream& a) override;
+
+private:
+	bool Init() override;
+	bool m_binit = false;
+	FEIdealGasPressure* m_load = nullptr;
+};

FEBioMech/FEIdealGasPressure.cpp

@@ -26,7 +26,7 @@ SOFTWARE.*/
 #include "stdafx.h"
 #include "FEIdealGasPressure.h"
 #include "FEBioMech.h"
-#include <FECore/FEFacetSet.h>
+#include <FECore/log.h>
 
 BEGIN_FECORE_CLASS(FEIdealGasPressure, FESurfaceLoad)
 	ADD_PARAMETER(m_initialPressure, "P0")->setUnits(UNIT_PRESSURE);
@@ -39,6 +39,8 @@ FEIdealGasPressure::FEIdealGasPressure(FEModel* pfem) : FESurfaceLoad(pfem)
 	m_initialPressure = 0.0;
 	m_bsymm = true;
 	m_bshellb = false;
+	m_currentPressure = 0;
+	m_currentVolume = 0;
 }
 
 bool FEIdealGasPressure::Init()
@@ -64,6 +66,8 @@ void FEIdealGasPressure::Activate()
 	m_currentVolume = m_initialVolume;
 
 	m_currentPressure = m_initialPressure;
+
+	feLogDebug("initial volume: %lg\n", m_initialVolume);
 
 	FESurfaceLoad::Activate();
 }
@@ -72,6 +76,9 @@ void FEIdealGasPressure::Update()
 {
 	m_currentVolume = CalculateSurfaceVolume(GetSurface());
 	m_currentPressure = m_initialPressure * (m_initialVolume / m_currentVolume);
+
+	feLogDebug("volume: %lg, pressure: %lg\n", m_currentVolume, m_currentPressure);
+
 	FESurfaceLoad::Update();
 }
 

FEBioMech/FEIdealGasPressure.h

@@ -44,6 +44,9 @@ class FEIdealGasPressure : public FESurfaceLoad
 
 	void Update() override;
 
+public:
+	double GetCurrentPressure() const { return m_currentPressure; }
+
 public:
 	//! calculate residual
 	void LoadVector(FEGlobalVector& R) override;