Removed the (erroneous) calculation of strain energy density in 'viscoelastic' and 'uncoupled viscoelastic' materials.

Author: gateshian

FEBioMech/FEUncoupledViscoElasticMaterial.cpp

@@ -153,7 +153,7 @@ tens4ds FEUncoupledViscoElasticMaterial::DevTangent(FEMaterialPoint& pt)
 
 //-----------------------------------------------------------------------------
 //! Strain energy density function
-double FEUncoupledViscoElasticMaterial::DevStrainEnergyDensity(FEMaterialPoint& mp)
+/*double FEUncoupledViscoElasticMaterial::DevStrainEnergyDensity(FEMaterialPoint& mp)
 {
 	// get the viscoelastic point data
 	FEViscoElasticMaterialPoint& pt = *mp.ExtractData<FEViscoElasticMaterialPoint>();
@@ -256,4 +256,4 @@ bool FEUncoupledViscoElasticMaterial::SeriesStretchExponent(FEMaterialPoint& mp)
     et.m_F = Fsafe; et.m_J = Jsafe;
     return true;
 }
-
+ */

FEBioMech/FEUncoupledViscoElasticMaterial.h

@@ -61,10 +61,10 @@ class FEUncoupledViscoElasticMaterial :	public FEUncoupledMaterial
 	tens4ds DevTangent(FEMaterialPoint& pt) override;
 
 	//! deviatoric strain energy density function
-	double DevStrainEnergyDensity(FEMaterialPoint& pt) override;
+    double DevStrainEnergyDensity(FEMaterialPoint& pt) override { return 0; }
 
     //! calculate exponent of right-stretch tensor in series spring
-    bool SeriesStretchExponent(FEMaterialPoint& pt);
+//    bool SeriesStretchExponent(FEMaterialPoint& pt);
 
 	//! returns a pointer to a new material point object
 	FEMaterialPointData* CreateMaterialPointData() override;

FEBioMech/FEViscoElasticMaterial.cpp

@@ -222,6 +222,7 @@ tens4ds FEViscoElasticMaterial::Tangent(FEMaterialPoint& pt)
 
 //-----------------------------------------------------------------------------
 //! Strain energy density function
+/*
 double FEViscoElasticMaterial::StrainEnergyDensity(FEMaterialPoint& mp)
 {
 	// get the viscoelastic point data
@@ -325,4 +326,4 @@ bool FEViscoElasticMaterial::SeriesStretchExponent(FEMaterialPoint& mp)
     et.m_F = Fsafe; et.m_J = Jsafe;
     return true;
 }
-
+*/

FEBioMech/FEViscoElasticMaterial.h

@@ -95,10 +95,10 @@ class FEViscoElasticMaterial :	public FEElasticMaterial
 	tens4ds Tangent(FEMaterialPoint& pt) override;
 
 	//! strain energy density
-	double StrainEnergyDensity(FEMaterialPoint& pt) override;
+    double StrainEnergyDensity(FEMaterialPoint& pt) override { return 0; }
 
     //! calculate exponent of right-stretch tensor in series spring
-    bool SeriesStretchExponent(FEMaterialPoint& pt);
+//    bool SeriesStretchExponent(FEMaterialPoint& pt);
 
     // returns a pointer to a new material point object
 	FEMaterialPointData* CreateMaterialPointData() override;