Fem: Add radiation heat transfer (#13466)

* FEM: Add radiation heat transfer * FEM: Add radiation heat transfer to CalculiX writer * FEM: Update .inp files * Fem: Rename heat flux ui members --------- Co-authored-by: FEA-eng <59876896+FEA-eng@users.noreply.github.com>
2024-04-29 12:53:59 -03:00
parent 5d69945f6e
commit ff285b6c60
40 changed files with 278 additions and 18 deletions
--- a/src/Mod/Fem/App/FemConstraintHeatflux.cpp
+++ b/src/Mod/Fem/App/FemConstraintHeatflux.cpp
@@ -31,19 +31,20 @@ using namespace Fem;

 PROPERTY_SOURCE(Fem::ConstraintHeatflux, Fem::Constraint)

-static const char* ConstraintTypes[] = {"DFlux", "Convection", nullptr};
+static const char* ConstraintTypes[] = {"DFlux", "Convection", "Radiation", nullptr};

 ConstraintHeatflux::ConstraintHeatflux()
 {
    ADD_PROPERTY(AmbientTemp, (0.0));
    /*ADD_PROPERTY(FaceTemp,(0.0));*/
    ADD_PROPERTY(FilmCoef, (0.0));
+    ADD_PROPERTY(Emissivity, (0.0));
    ADD_PROPERTY(DFlux, (0.0));
    ADD_PROPERTY_TYPE(ConstraintType,
                      (1),
                      "ConstraintHeatflux",
                      (App::PropertyType)(App::Prop_None),
-                      "Type of constraint, surface convection or surface heat flux");
+                      "Type of constraint, surface convection, radiation or surface heat flux");
    ConstraintType.setEnums(ConstraintTypes);
 }

--- a/src/Mod/Fem/App/FemConstraintHeatflux.h
+++ b/src/Mod/Fem/App/FemConstraintHeatflux.h
@@ -42,6 +42,7 @@ public:
    App::PropertyFloat AmbientTemp;
    /*App::PropertyFloat FaceTemp;*/
    App::PropertyFloat FilmCoef;
+    App::PropertyFloat Emissivity;
    App::PropertyFloat DFlux;
    App::PropertyEnumeration ConstraintType;

--- a/src/Mod/Fem/Gui/Command.cpp
+++ b/src/Mod/Fem/Gui/Command.cpp
@@ -578,6 +578,9 @@ void CmdFemConstraintHeatflux::activated(int)
    doCommand(Doc,
              "App.activeDocument().%s.FilmCoef = 10.0",
              FeatName.c_str());  // OvG: set default not equal to 0
+    doCommand(Doc,
+              "App.activeDocument().%s.Emissivity = 1.0",
+              FeatName.c_str());  // OvG: set default not equal to 0
    doCommand(Doc,
              "App.activeDocument().%s.Scale = 1",
              FeatName.c_str());  // OvG: set initial scale to 1
--- a/src/Mod/Fem/Gui/TaskFemConstraintHeatflux.cpp
+++ b/src/Mod/Fem/Gui/TaskFemConstraintHeatflux.cpp
@@ -62,13 +62,14 @@ TaskFemConstraintHeatflux::TaskFemConstraintHeatflux(
            &TaskFemConstraintHeatflux::onReferenceDeleted);

    connect(ui->rb_convection, &QRadioButton::clicked, this, &TaskFemConstraintHeatflux::Conv);
+    connect(ui->rb_radiation, &QRadioButton::clicked, this, &TaskFemConstraintHeatflux::Rad);
    connect(ui->rb_dflux, &QRadioButton::clicked, this, &TaskFemConstraintHeatflux::Flux);

    connect(ui->if_heatflux,
            qOverload<double>(&InputField::valueChanged),
            this,
            &TaskFemConstraintHeatflux::onHeatFluxChanged);
-    connect(ui->if_ambienttemp,
+    connect(ui->if_ambienttemp_conv,
            qOverload<double>(&InputField::valueChanged),
            this,
            &TaskFemConstraintHeatflux::onAmbientTempChanged);
@@ -76,6 +77,14 @@ TaskFemConstraintHeatflux::TaskFemConstraintHeatflux(
            qOverload<double>(&InputField::valueChanged),
            this,
            &TaskFemConstraintHeatflux::onFilmCoefChanged);
+    connect(ui->if_emissivity,
+            qOverload<double>(&InputField::valueChanged),
+            this,
+            &TaskFemConstraintHeatflux::onEmissivityChanged);
+    connect(ui->if_ambienttemp_rad,
+            qOverload<double>(&InputField::valueChanged),
+            this,
+            &TaskFemConstraintHeatflux::onAmbientTempChanged);
    connect(ui->lw_references,
            &QListWidget::itemClicked,
            this,
@@ -84,9 +93,11 @@ TaskFemConstraintHeatflux::TaskFemConstraintHeatflux(
    this->groupLayout()->addWidget(proxy);

    // Temporarily prevent unnecessary feature recomputes
-    ui->if_ambienttemp->blockSignals(true);
+    ui->if_ambienttemp_conv->blockSignals(true);
    // ui->if_facetemp->blockSignals(true);
    ui->if_filmcoef->blockSignals(true);
+    ui->if_emissivity->blockSignals(true);
+    ui->if_ambienttemp_rad->blockSignals(true);
    ui->lw_references->blockSignals(true);
    ui->btnAdd->blockSignals(true);
    ui->btnRemove->blockSignals(true);
@@ -98,26 +109,41 @@ TaskFemConstraintHeatflux::TaskFemConstraintHeatflux(
    std::vector<std::string> SubElements = pcConstraint->References.getSubValues();

    // Fill data into dialog elements
-    ui->if_ambienttemp->setMinimum(0);
-    ui->if_ambienttemp->setMaximum(FLOAT_MAX);
+    ui->if_ambienttemp_conv->setMinimum(0);
+    ui->if_ambienttemp_conv->setMaximum(FLOAT_MAX);

    ui->if_filmcoef->setMinimum(0);
    ui->if_filmcoef->setMaximum(FLOAT_MAX);

+    ui->if_emissivity->setMinimum(0);
+    ui->if_emissivity->setMaximum(FLOAT_MAX);
+
+    ui->if_ambienttemp_rad->setMinimum(0);
+    ui->if_ambienttemp_rad->setMaximum(FLOAT_MAX);
+
    std::string constraint_type = pcConstraint->ConstraintType.getValueAsString();
    if (constraint_type == "Convection") {
        ui->rb_convection->setChecked(true);
        ui->sw_heatflux->setCurrentIndex(0);
        Base::Quantity t =
            Base::Quantity(pcConstraint->AmbientTemp.getValue(), Base::Unit::Temperature);
-        ui->if_ambienttemp->setValue(t);
+        ui->if_ambienttemp_conv->setValue(t);
        Base::Quantity f = Base::Quantity(pcConstraint->FilmCoef.getValue(),
                                          Base::Unit::ThermalTransferCoefficient);
        ui->if_filmcoef->setValue(f);
    }
+    else if (constraint_type == "Radiation") {
+        ui->rb_radiation->setChecked(true);
+        ui->sw_heatflux->setCurrentIndex(1);
+        Base::Quantity t =
+            Base::Quantity(pcConstraint->AmbientTemp.getValue(), Base::Unit::Temperature);
+        ui->if_ambienttemp_rad->setValue(t);
+        Base::Quantity e = Base::Quantity(pcConstraint->Emissivity.getValue(), Base::Unit());
+        ui->if_emissivity->setValue(e);
+    }
    else if (constraint_type == "DFlux") {
        ui->rb_dflux->setChecked(true);
-        ui->sw_heatflux->setCurrentIndex(1);
+        ui->sw_heatflux->setCurrentIndex(2);
        Base::Quantity c = Base::Quantity(pcConstraint->DFlux.getValue(), Base::Unit::HeatFlux);
        ui->if_heatflux->setValue(c);
    }
@@ -134,9 +160,11 @@ TaskFemConstraintHeatflux::TaskFemConstraintHeatflux(
    buttonGroup->addButton(ui->btnAdd, (int)SelectionChangeModes::refAdd);
    buttonGroup->addButton(ui->btnRemove, (int)SelectionChangeModes::refRemove);

-    ui->if_ambienttemp->blockSignals(false);
+    ui->if_ambienttemp_conv->blockSignals(false);
    // ui->if_facetemp->blockSignals(false);
    ui->if_filmcoef->blockSignals(false);
+    ui->if_emissivity->blockSignals(false);
+    ui->if_ambienttemp_rad->blockSignals(false);
    ui->lw_references->blockSignals(false);
    ui->btnAdd->blockSignals(false);
    ui->btnRemove->blockSignals(false);
@@ -169,6 +197,13 @@ void TaskFemConstraintHeatflux::onFilmCoefChanged(double val)
    pcConstraint->FilmCoef.setValue(val);  // [W]/[[m^2]/[K]]
 }

+void TaskFemConstraintHeatflux::onEmissivityChanged(double val)
+{
+    Fem::ConstraintHeatflux* pcConstraint =
+        static_cast<Fem::ConstraintHeatflux*>(ConstraintView->getObject());
+    pcConstraint->Emissivity.setValue(val);  // [-]
+}
+
 void TaskFemConstraintHeatflux::onHeatFluxChanged(double val)
 {
    Fem::ConstraintHeatflux* pcConstraint =
@@ -186,7 +221,7 @@ void TaskFemConstraintHeatflux::Conv()
                            name.c_str(),
                            get_constraint_type().c_str());
    Base::Quantity t = Base::Quantity(300, Base::Unit::Temperature);
-    ui->if_ambienttemp->setValue(t);
+    ui->if_ambienttemp_conv->setValue(t);
    pcConstraint->AmbientTemp.setValue(300);
    Base::Quantity f = Base::Quantity(10, Base::Unit::ThermalTransferCoefficient);
    ui->if_filmcoef->setValue(f);
@@ -194,6 +229,24 @@ void TaskFemConstraintHeatflux::Conv()
    ui->sw_heatflux->setCurrentIndex(0);
 }

+void TaskFemConstraintHeatflux::Rad()
+{
+    Fem::ConstraintHeatflux* pcConstraint =
+        static_cast<Fem::ConstraintHeatflux*>(ConstraintView->getObject());
+    std::string name = ConstraintView->getObject()->getNameInDocument();
+    Gui::Command::doCommand(Gui::Command::Doc,
+                            "App.ActiveDocument.%s.ConstraintType = %s",
+                            name.c_str(),
+                            get_constraint_type().c_str());
+    Base::Quantity t = Base::Quantity(300, Base::Unit::Temperature);
+    ui->if_ambienttemp_rad->setValue(t);
+    pcConstraint->AmbientTemp.setValue(300);
+    Base::Quantity e = Base::Quantity(1, Base::Unit());
+    ui->if_emissivity->setValue(e);
+    pcConstraint->Emissivity.setValue(1);
+    ui->sw_heatflux->setCurrentIndex(1);
+}
+
 void TaskFemConstraintHeatflux::Flux()
 {
    Fem::ConstraintHeatflux* pcConstraint =
@@ -206,7 +259,7 @@ void TaskFemConstraintHeatflux::Flux()
    Base::Quantity c = Base::Quantity(0, Base::Unit::HeatFlux);
    ui->if_heatflux->setValue(c);
    pcConstraint->DFlux.setValue(0);
-    ui->sw_heatflux->setCurrentIndex(1);
+    ui->sw_heatflux->setCurrentIndex(2);
 }

 void TaskFemConstraintHeatflux::addToSelection()
@@ -363,7 +416,13 @@ const std::string TaskFemConstraintHeatflux::getReferences() const

 double TaskFemConstraintHeatflux::getAmbientTemp() const
 {
-    Base::Quantity temperature = ui->if_ambienttemp->getQuantity();
+    Base::Quantity temperature;
+    if (ui->rb_convection->isChecked()) {
+        temperature = ui->if_ambienttemp_conv->getQuantity();
+    }
+    else if (ui->rb_radiation->isChecked()) {
+        temperature = ui->if_ambienttemp_rad->getQuantity();
+    }
    double temperature_in_kelvin = temperature.getValueAs(Base::Quantity::Kelvin);
    return temperature_in_kelvin;
 }
@@ -376,12 +435,22 @@ double TaskFemConstraintHeatflux::getFilmCoef() const
    return filmcoef_in_units;
 }

+double TaskFemConstraintHeatflux::getEmissivity() const
+{
+    Base::Quantity emissivity = ui->if_emissivity->getQuantity();
+    double emissivity_in_units = emissivity.getValueAs(Base::Quantity(1.0, Base::Unit()));
+    return emissivity_in_units;
+}
+
 std::string TaskFemConstraintHeatflux::get_constraint_type() const
 {
    std::string type;
    if (ui->rb_convection->isChecked()) {
        type = "\"Convection\"";
    }
+    else if (ui->rb_radiation->isChecked()) {
+        type = "\"Radiation\"";
+    }
    else if (ui->rb_dflux->isChecked()) {
        type = "\"DFlux\"";
    }
@@ -397,11 +466,15 @@ void TaskFemConstraintHeatflux::changeEvent(QEvent* e)
 {
    TaskBox::changeEvent(e);
    if (e->type() == QEvent::LanguageChange) {
-        ui->if_ambienttemp->blockSignals(true);
+        ui->if_ambienttemp_conv->blockSignals(true);
        ui->if_filmcoef->blockSignals(true);
+        ui->if_emissivity->blockSignals(true);
+        ui->if_ambienttemp_rad->blockSignals(true);
        ui->retranslateUi(proxy);
-        ui->if_ambienttemp->blockSignals(false);
+        ui->if_ambienttemp_conv->blockSignals(false);
        ui->if_filmcoef->blockSignals(false);
+        ui->if_emissivity->blockSignals(true);
+        ui->if_ambienttemp_rad->blockSignals(false);
    }
 }

@@ -465,6 +538,11 @@ bool TaskDlgFemConstraintHeatflux::accept()
                                name.c_str(),
                                parameterHeatflux->getFilmCoef());

+        Gui::Command::doCommand(Gui::Command::Doc,
+                                "App.ActiveDocument.%s.Emissivity = %f",
+                                name.c_str(),
+                                parameterHeatflux->getEmissivity());
+
        scale = parameterHeatflux->getScale();  // OvG: determine modified scale
        Gui::Command::doCommand(Gui::Command::Doc,
                                "App.ActiveDocument.%s.Scale = %s",
--- a/src/Mod/Fem/Gui/TaskFemConstraintHeatflux.h
+++ b/src/Mod/Fem/Gui/TaskFemConstraintHeatflux.h
@@ -48,6 +48,7 @@ public:
    double getAmbientTemp() const;
    /*double getFaceTemp(void) const;*/
    double getFilmCoef() const;
+    double getEmissivity() const;
    std::string get_constraint_type() const;
    const std::string getReferences() const override;

@@ -56,8 +57,10 @@ private Q_SLOTS:
    void onAmbientTempChanged(double val);
    /*void onFaceTempChanged(double val);*/
    void onFilmCoefChanged(double val);
+    void onEmissivityChanged(double val);
    void onHeatFluxChanged(double val);
    void Conv();
+    void Rad();
    void Flux();
    void addToSelection() override;
    void removeFromSelection() override;
--- a/src/Mod/Fem/Gui/TaskFemConstraintHeatflux.ui
+++ b/src/Mod/Fem/Gui/TaskFemConstraintHeatflux.ui
@@ -69,6 +69,13 @@
       </property>
      </widget>
     </item>
+     <item>
+      <widget class="QRadioButton" name="rb_radiation">
+       <property name="text">
+        <string>Surface Radiation</string>
+       </property>
+      </widget>
+     </item>
     <item>
      <widget class="QRadioButton" name="rb_dflux">
       <property name="text">
@@ -81,7 +88,7 @@
   <item>
    <widget class="QStackedWidget" name="sw_heatflux">
     <property name="currentIndex">
-      <number>1</number>
+      <number>0</number>
     </property>
     <widget class="QWidget" name="page">
      <layout class="QVBoxLayout" name="verticalLayout_4">
@@ -112,16 +119,16 @@
          </layout>
         </item>
         <item>
-          <layout class="QHBoxLayout" name="layoutAmbientTemp">
+          <layout class="QHBoxLayout" name="layoutAmbientTempCov">
           <item>
-            <widget class="QLabel" name="lbl_ambienttemp">
+            <widget class="QLabel" name="lbl_ambienttemp_conv">
             <property name="text">
              <string>Ambient Temperature</string>
             </property>
            </widget>
           </item>
           <item>
-            <widget class="Gui::InputField" name="if_ambienttemp">
+            <widget class="Gui::InputField" name="if_ambienttemp_conv">
             <property name="text">
              <string>300 K</string>
             </property>
@@ -139,6 +146,65 @@
       </item>
      </layout>
     </widget>
+     <property name="currentIndex">
+      <number>1</number>
+     </property>
+     <widget class="QWidget" name="page_1">
+      <layout class="QVBoxLayout" name="verticalLayout_5">
+       <item>
+        <layout class="QVBoxLayout" name="verticalLayout_6">
+         <item>
+          <layout class="QHBoxLayout" name="layoutEmissivity">
+           <item>
+            <widget class="QLabel" name="lbl_emissivity">
+             <property name="text">
+              <string>Emissivity</string>
+             </property>
+            </widget>
+           </item>
+           <item>
+            <widget class="Gui::InputField" name="if_emissivity">
+             <property name="text">
+              <string>1 </string>
+             </property>
+             <property name="quantity" stdset="0">
+              <double>1.000000000000000</double>
+             </property>
+            </widget>
+           </item>
+          </layout>
+         </item>
+         <item>
+          <layout class="QHBoxLayout" name="layoutAmbientTempRad">
+           <item>
+            <widget class="QLabel" name="lbl_ambienttemp_rad">
+             <property name="text">
+              <string>Ambient Temperature</string>
+             </property>
+            </widget>
+           </item>
+           <item>
+            <widget class="Gui::InputField" name="if_ambienttemp_rad">
+             <property name="text">
+              <string>300 K</string>
+             </property>
+             <property name="unit" stdset="0">
+              <string notr="true">K</string>
+             </property>
+             <property name="quantity" stdset="0">
+              <double>300.000000000000000</double>
+             </property>
+            </widget>
+           </item>
+          </layout>
+         </item>
+        </layout>
+       </item>
+      </layout>
+     </widget>
+     <property name="currentIndex">
+      <number>2</number>
+     </property>
     <widget class="QWidget" name="page_2">
      <layout class="QVBoxLayout" name="verticalLayout_3">
       <item>
--- a/src/Mod/Fem/femsolver/calculix/write_constraint_heatflux.py
+++ b/src/Mod/Fem/femsolver/calculix/write_constraint_heatflux.py
@@ -54,6 +54,15 @@ def write_meshdata_constraint(f, femobj, heatflux_obj, ccxwriter):
            heatflux_obj.AmbientTemp,
            heatflux_obj.FilmCoef * 0.001
        )
+
+    elif heatflux_obj.ConstraintType == "Radiation":
+        heatflux_key_word = "RADIATE"
+        heatflux_facetype = "R"
+        heatflux_values = "{:.13G},{:.13G}".format(
+            heatflux_obj.AmbientTemp,
+            heatflux_obj.Emissivity
+        )
+
    elif heatflux_obj.ConstraintType == "DFlux":
        heatflux_key_word = "DFLUX"
        heatflux_facetype = "S"
--- a/src/Mod/Fem/femsolver/calculix/write_femelement_material.py
+++ b/src/Mod/Fem/femsolver/calculix/write_femelement_material.py
@@ -48,6 +48,9 @@ def write_femelement_material(f, ccxwriter):
            return True
        return False

+    f.write("\n** Physical constants for SI(mm) unit system with Kelvins\n")
+    f.write("*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11\n")
+
    f.write("\n{}\n".format(59 * "*"))
    f.write("** Materials\n")
    f.write("** see information about units at file end\n")
--- a/src/Mod/Fem/femtest/data/calculix/box_frequency.inp
+++ b/src/Mod/Fem/femtest/data/calculix/box_frequency.inp
@@ -428,6 +428,9 @@ Evolumes
 *ELSET,ELSET=MechanicalMaterialSolid
 Evolumes

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/box_static.inp
+++ b/src/Mod/Fem/femtest/data/calculix/box_static.inp
@@ -474,6 +474,9 @@ Evolumes
 68,
 69,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/ccx_buckling_flexuralbuckling.inp
+++ b/src/Mod/Fem/femtest/data/calculix/ccx_buckling_flexuralbuckling.inp
@@ -730,6 +730,9 @@ Evolumes
 393,
 410,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_beam_circle.inp
+++ b/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_beam_circle.inp
@@ -45,6 +45,9 @@ Eedges
 *NSET,NSET=ConstraintFixed
 1,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_beam_pipe.inp
+++ b/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_beam_pipe.inp
@@ -45,6 +45,9 @@ Eedges
 *NSET,NSET=ConstraintFixed
 1,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_beam_rect.inp
+++ b/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_beam_rect.inp
@@ -45,6 +45,9 @@ Eedges
 *NSET,NSET=ConstraintFixed
 1,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_ele_hexa20.inp
+++ b/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_ele_hexa20.inp
@@ -370,6 +370,9 @@ Evolumes
 96,
 97,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_ele_quad4.inp
+++ b/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_ele_quad4.inp
@@ -70,6 +70,9 @@ Efaces
 2,
 5,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_ele_quad8.inp
+++ b/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_ele_quad8.inp
@@ -58,6 +58,9 @@ Efaces
 6,
 7,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_ele_seg2.inp
+++ b/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_ele_seg2.inp
@@ -187,6 +187,9 @@ Eedges
 *NSET,NSET=ConstraintFixed
 1,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_ele_seg3.inp
+++ b/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_ele_seg3.inp
@@ -45,6 +45,9 @@ Eedges
 *NSET,NSET=ConstraintFixed
 1,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_ele_tria3.inp
+++ b/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_ele_tria3.inp
@@ -1546,6 +1546,9 @@ Efaces
 9,
 10,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_ele_tria6.inp
+++ b/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_ele_tria6.inp
@@ -276,6 +276,9 @@ Efaces
 6,
 7,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_faceload.inp
+++ b/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_faceload.inp
@@ -344,6 +344,9 @@ Evolumes
 194,
 195,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_nodeload.inp
+++ b/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_nodeload.inp
@@ -344,6 +344,9 @@ Evolumes
 194,
 195,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_prescribeddisplacement.inp
+++ b/src/Mod/Fem/femtest/data/calculix/ccx_cantilever_prescribeddisplacement.inp
@@ -362,6 +362,9 @@ Evolumes
 190,
 191,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/constraint_centrif.inp
+++ b/src/Mod/Fem/femtest/data/calculix/constraint_centrif.inp
@@ -18721,6 +18721,9 @@ Evolumes
 2745,
 2746,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/constraint_contact_shell_shell.inp
+++ b/src/Mod/Fem/femtest/data/calculix/constraint_contact_shell_shell.inp
@@ -38348,6 +38348,9 @@ Efaces
 15603,S2
 15604,S2

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/constraint_contact_solid_solid.inp
+++ b/src/Mod/Fem/femtest/data/calculix/constraint_contact_solid_solid.inp
@@ -5556,6 +5556,9 @@ Evolumes
 1726,S1
 2345,S4

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/constraint_sectionprint.inp
+++ b/src/Mod/Fem/femtest/data/calculix/constraint_sectionprint.inp
@@ -3386,6 +3386,9 @@ Evolumes
 1819,S4
 1823,S4

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/constraint_selfweight_cantilever.inp
+++ b/src/Mod/Fem/femtest/data/calculix/constraint_selfweight_cantilever.inp
@@ -2136,6 +2136,9 @@ Evolumes
 12,
 13,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/constraint_tie.inp
+++ b/src/Mod/Fem/femtest/data/calculix/constraint_tie.inp
@@ -18591,6 +18591,9 @@ Evolumes
 9684,S1
 10647,S2

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/constraint_transform_beam_hinged.inp
+++ b/src/Mod/Fem/femtest/data/calculix/constraint_transform_beam_hinged.inp
@@ -3615,6 +3615,9 @@ Evolumes
 241,
 2088,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/constraint_transform_torque.inp
+++ b/src/Mod/Fem/femtest/data/calculix/constraint_transform_torque.inp
@@ -10959,6 +10959,9 @@ Evolumes
 3933,
 3934,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/frequency_beamsimple.inp
+++ b/src/Mod/Fem/femtest/data/calculix/frequency_beamsimple.inp
@@ -17021,6 +17021,9 @@ Evolumes
 17,
 18,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/material_multiple_bendingbeam_fiveboxes.inp
+++ b/src/Mod/Fem/femtest/data/calculix/material_multiple_bendingbeam_fiveboxes.inp
@@ -27607,6 +27607,9 @@ Evolumes
 4152,
 4153,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/material_multiple_bendingbeam_fivefaces.inp
+++ b/src/Mod/Fem/femtest/data/calculix/material_multiple_bendingbeam_fivefaces.inp
@@ -2518,6 +2518,9 @@ Efaces
 460,
 461,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/material_multiple_tensionrod_twoboxes.inp
+++ b/src/Mod/Fem/femtest/data/calculix/material_multiple_tensionrod_twoboxes.inp
@@ -1210,6 +1210,9 @@ Evolumes
 296,
 297,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/material_nonlinear.inp
+++ b/src/Mod/Fem/femtest/data/calculix/material_nonlinear.inp
@@ -19985,6 +19985,9 @@ Evolumes
 1328,
 1329,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/square_pipe_end_twisted_edgeforces.inp
+++ b/src/Mod/Fem/femtest/data/calculix/square_pipe_end_twisted_edgeforces.inp
@@ -2544,6 +2544,9 @@ Efaces
 1402,
 1407,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/square_pipe_end_twisted_nodeforces.inp
+++ b/src/Mod/Fem/femtest/data/calculix/square_pipe_end_twisted_nodeforces.inp
@@ -2544,6 +2544,9 @@ Efaces
 1402,
 1407,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end
--- a/src/Mod/Fem/femtest/data/calculix/thermomech_bimetall.inp
+++ b/src/Mod/Fem/femtest/data/calculix/thermomech_bimetall.inp
@@ -7040,6 +7040,9 @@ Evolumes
 1520,
 1521,

+** Physical constants for SI(mm) unit system with Kelvins
+*PHYSICAL CONSTANTS, ABSOLUTE ZERO=0, STEFAN BOLTZMANN=5.670374419e-11
+
 ***********************************************************
 ** Materials
 ** see information about units at file end