diff --git a/src/Mod/Fem/femexamples/ccx_cantilever_base_solid.py b/src/Mod/Fem/femexamples/ccx_cantilever_base_solid.py
index f92fc3377a..0b15acc2c1 100644
--- a/src/Mod/Fem/femexamples/ccx_cantilever_base_solid.py
+++ b/src/Mod/Fem/femexamples/ccx_cantilever_base_solid.py
@@ -84,7 +84,6 @@ def setup_cantilever_base_solid(doc=None, solvertype="ccxtools"):
     mat["Name"] = "CalculiX-Steel"
     mat["YoungsModulus"] = "210000 MPa"
     mat["PoissonRatio"] = "0.30"
-    mat["Density"] = "7900 kg/m^3"
     material_obj.Material = mat
     analysis.addObject(material_obj)
 
diff --git a/src/Mod/Fem/femexamples/constraint_transform_beam_hinged.py b/src/Mod/Fem/femexamples/constraint_transform_beam_hinged.py
index 5bd5c4e158..bb6e326c06 100644
--- a/src/Mod/Fem/femexamples/constraint_transform_beam_hinged.py
+++ b/src/Mod/Fem/femexamples/constraint_transform_beam_hinged.py
@@ -141,8 +141,6 @@ def setup(doc=None, solvertype="ccxtools"):
     mat["Name"] = "CalculiX-Steel"
     mat["YoungsModulus"] = "210000 MPa"
     mat["PoissonRatio"] = "0.30"
-    mat["Density"] = "7900 kg/m^3"
-    mat["ThermalExpansionCoefficient"] = "0.012 mm/m/K"
     material_obj.Material = mat
     analysis.addObject(material_obj)
 
diff --git a/src/Mod/Fem/femexamples/material_multiple_bendingbeam_fiveboxes.py b/src/Mod/Fem/femexamples/material_multiple_bendingbeam_fiveboxes.py
index 02c7afbff9..7f90c44a4b 100644
--- a/src/Mod/Fem/femexamples/material_multiple_bendingbeam_fiveboxes.py
+++ b/src/Mod/Fem/femexamples/material_multiple_bendingbeam_fiveboxes.py
@@ -148,7 +148,6 @@ def setup(doc=None, solvertype="ccxtools"):
     mat["Name"] = "Concrete-Generic"
     mat["YoungsModulus"] = "32000 MPa"
     mat["PoissonRatio"] = "0.17"
-    mat["Density"] = "0 kg/m^3"
     material_obj1.Material = mat
     analysis.addObject(material_obj1)
 
@@ -158,7 +157,6 @@ def setup(doc=None, solvertype="ccxtools"):
     mat["Name"] = "PLA"
     mat["YoungsModulus"] = "3640 MPa"
     mat["PoissonRatio"] = "0.36"
-    mat["Density"] = "0 kg/m^3"
     material_obj2.Material = mat
     analysis.addObject(material_obj2)
 
@@ -168,7 +166,6 @@ def setup(doc=None, solvertype="ccxtools"):
     mat["Name"] = "Steel-Generic"
     mat["YoungsModulus"] = "200000 MPa"
     mat["PoissonRatio"] = "0.30"
-    mat["Density"] = "7900 kg/m^3"
     material_obj3.Material = mat
     analysis.addObject(material_obj3)
 
diff --git a/src/Mod/Fem/femexamples/material_multiple_bendingbeam_fivefaces.py b/src/Mod/Fem/femexamples/material_multiple_bendingbeam_fivefaces.py
index 0d57d93635..d5bac29b75 100644
--- a/src/Mod/Fem/femexamples/material_multiple_bendingbeam_fivefaces.py
+++ b/src/Mod/Fem/femexamples/material_multiple_bendingbeam_fivefaces.py
@@ -141,7 +141,6 @@ def setup(doc=None, solvertype="ccxtools"):
     mat["Name"] = "Concrete-Generic"
     mat["YoungsModulus"] = "32000 MPa"
     mat["PoissonRatio"] = "0.17"
-    mat["Density"] = "0 kg/m^3"
     material_obj1.Material = mat
     analysis.addObject(material_obj1)
 
@@ -154,7 +153,6 @@ def setup(doc=None, solvertype="ccxtools"):
     mat["Name"] = "PLA"
     mat["YoungsModulus"] = "3640 MPa"
     mat["PoissonRatio"] = "0.36"
-    mat["Density"] = "0 kg/m^3"
     material_obj2.Material = mat
     analysis.addObject(material_obj2)
 
@@ -164,7 +162,6 @@ def setup(doc=None, solvertype="ccxtools"):
     mat["Name"] = "Steel-Generic"
     mat["YoungsModulus"] = "200000 MPa"
     mat["PoissonRatio"] = "0.30"
-    mat["Density"] = "7900 kg/m^3"
     material_obj3.Material = mat
     analysis.addObject(material_obj3)
 
diff --git a/src/Mod/Fem/femexamples/material_multiple_tensionrod_twoboxes.py b/src/Mod/Fem/femexamples/material_multiple_tensionrod_twoboxes.py
index 89674f8db6..8ae258cf40 100644
--- a/src/Mod/Fem/femexamples/material_multiple_tensionrod_twoboxes.py
+++ b/src/Mod/Fem/femexamples/material_multiple_tensionrod_twoboxes.py
@@ -132,7 +132,6 @@ def setup(doc=None, solvertype="ccxtools"):
     mat["Name"] = "Aluminium-Generic"
     mat["YoungsModulus"] = "70000 MPa"
     mat["PoissonRatio"] = "0.35"
-    mat["Density"] = "2700  kg/m^3"
     material_obj_low.Material = mat
     material_obj_low.References = [(boxlow, "Solid1")]
     analysis.addObject(material_obj_low)
@@ -142,7 +141,6 @@ def setup(doc=None, solvertype="ccxtools"):
     mat["Name"] = "Steel-Generic"
     mat["YoungsModulus"] = "200000 MPa"
     mat["PoissonRatio"] = "0.30"
-    mat["Density"] = "7980 kg/m^3"
     material_obj_upp.Material = mat
     material_obj_upp.References = [(boxupp, "Solid1")]
     analysis.addObject(material_obj_upp)
diff --git a/src/Mod/Fem/femexamples/material_nl_platewithhole.py b/src/Mod/Fem/femexamples/material_nl_platewithhole.py
index 2941aae1f8..e81d72f71a 100644
--- a/src/Mod/Fem/femexamples/material_nl_platewithhole.py
+++ b/src/Mod/Fem/femexamples/material_nl_platewithhole.py
@@ -143,7 +143,6 @@ def setup(doc=None, solvertype="ccxtools"):
     matprop["Name"] = "CalculiX-Steel"
     matprop["YoungsModulus"] = "210000 MPa"
     matprop["PoissonRatio"] = "0.30"
-    matprop["Density"] = "7900 kg/m^3"
     material_obj.Material = matprop
     analysis.addObject(material_obj)
 
diff --git a/src/Mod/Fem/femexamples/rc_wall_2d.py b/src/Mod/Fem/femexamples/rc_wall_2d.py
index 8baa02739c..aee9eb067d 100644
--- a/src/Mod/Fem/femexamples/rc_wall_2d.py
+++ b/src/Mod/Fem/femexamples/rc_wall_2d.py
@@ -133,7 +133,6 @@ def setup(doc=None, solvertype="ccxtools"):
     matrixprop["CompressiveStrength"] = "15.75 MPa"
     # make some hint on the possible angle units in material system
     matrixprop["AngleOfFriction"] = "30 deg"
-    matrixprop["Density"] = "2500 kg/m^3"
     reinfoprop = {}
     reinfoprop["Name"] = "Reinforcement-FIB-B500"
     reinfoprop["YieldStrength"] = "315 MPa"
diff --git a/src/Mod/Fem/femexamples/square_pipe_end_twisted_edgeforces.py b/src/Mod/Fem/femexamples/square_pipe_end_twisted_edgeforces.py
index b47461d95a..e34bd0da16 100644
--- a/src/Mod/Fem/femexamples/square_pipe_end_twisted_edgeforces.py
+++ b/src/Mod/Fem/femexamples/square_pipe_end_twisted_edgeforces.py
@@ -120,7 +120,6 @@ def setup(doc=None, solvertype="ccxtools"):
     mat["Name"] = "Steel-Generic"
     mat["YoungsModulus"] = "200000 MPa"
     mat["PoissonRatio"] = "0.30"
-    mat["Density"] = "7900 kg/m^3"
     material_obj.Material = mat
     analysis.addObject(material_obj)
 
diff --git a/src/Mod/Fem/femexamples/square_pipe_end_twisted_nodeforces.py b/src/Mod/Fem/femexamples/square_pipe_end_twisted_nodeforces.py
index 8ddebcfa7a..1b70ffb811 100644
--- a/src/Mod/Fem/femexamples/square_pipe_end_twisted_nodeforces.py
+++ b/src/Mod/Fem/femexamples/square_pipe_end_twisted_nodeforces.py
@@ -230,7 +230,6 @@ def setup(doc=None, solvertype="ccxtools"):
     mat["Name"] = "Steel-Generic"
     mat["YoungsModulus"] = "200000 MPa"
     mat["PoissonRatio"] = "0.30"
-    mat["Density"] = "7900 kg/m^3"
     material_obj.Material = mat
     analysis.addObject(material_obj)
 
diff --git a/src/Mod/Fem/femexamples/thermomech_bimetall.py b/src/Mod/Fem/femexamples/thermomech_bimetall.py
index a12c95ae2d..ef4c4b2a88 100644
--- a/src/Mod/Fem/femexamples/thermomech_bimetall.py
+++ b/src/Mod/Fem/femexamples/thermomech_bimetall.py
@@ -174,7 +174,6 @@ def setup(doc=None, solvertype="ccxtools"):
     mat["SpecificHeat"] = "510 J/kg/K"
     mat["ThermalConductivity"] = "13 W/m/K"
     mat["ThermalExpansionCoefficient"] = "0.0000012 m/m/K"
-    mat["Density"] = "1.00 kg/m^3"
     material_obj_top.Material = mat
     material_obj_top.References = [(geom_obj, "Solid2")]
     analysis.addObject(material_obj_top)