Update.

Author: carlos-souto

FAT-SM.vcxproj

@@ -128,6 +128,7 @@
     <ClCompile Include="Source\Gui\MainWindow.cpp" />
     <ClCompile Include="Source\Gui\ParisLawTab.cpp" />
     <ClCompile Include="Source\Gui\RainflowTab.cpp" />
+    <ClCompile Include="Source\Gui\RemainingLifeTab.cpp" />
     <ClCompile Include="Source\Gui\SnCurveTab.cpp" />
     <ClCompile Include="Source\Model\Detail.cpp" />
     <ClCompile Include="Source\Model\Project.cpp" />
@@ -173,6 +174,7 @@
     <QtMoc Include="Source\Gui\DamageTab.hpp" />
     <QtMoc Include="Source\Gui\DetailSelectionTab.hpp" />
     <QtMoc Include="Source\Gui\ParisLawTab.hpp" />
+    <QtMoc Include="Source\Gui\RemainingLifeTab.hpp" />
     <ClInclude Include="Source\Model.hpp" />
     <QtMoc Include="Source\Model\Project.hpp" />
     <QtMoc Include="Source\Model\LefmDetail.hpp" />

FAT-SM.vcxproj.filters

@@ -75,6 +75,9 @@
     <ClCompile Include="Source\Gui\ParisLawTab.cpp">
       <Filter>Source\Gui</Filter>
     </ClCompile>
+    <ClCompile Include="Source\Gui\RemainingLifeTab.cpp">
+      <Filter>Source\Gui</Filter>
+    </ClCompile>
   </ItemGroup>
   <ItemGroup>
     <None Include=".gitignore" />
@@ -182,6 +185,9 @@
     <QtMoc Include="Source\Gui\ParisLawTab.hpp">
       <Filter>Source\Gui</Filter>
     </QtMoc>
+    <QtMoc Include="Source\Gui\RemainingLifeTab.hpp">
+      <Filter>Source\Gui</Filter>
+    </QtMoc>
   </ItemGroup>
   <ItemGroup>
     <QtRcc Include="Resource\Resource.qrc">

Resource/StyleSheets/NavBar.qss

@@ -16,11 +16,12 @@ QListWidget::item {
 }
 
 QListWidget::item:selected {
-    color: #FFFFFF;
-    background-color: #006CD9;
+    color: #000000;
+    border: 1px solid black;
+    background-color: #E6E6E6;
 }
 
 QListWidget::item:hover:!selected {
     color: #000000;
-    background-color: #F0F0F0;
+    background-color: #E6E6E6;
 }

Source/Core/Matrix.cpp

@@ -1,5 +1,6 @@
 #include <string>
 #include <QString>
+#include <algorithm>
 #include <stdexcept>
 #include "Core/Matrix.hpp"
 
@@ -55,6 +56,21 @@ const T& Matrix<T>::at(int i, int j) const {
     return _data[i*_columnCount + j];
 }
 
+template<typename T>
+T& Matrix<T>::max() {
+    if (_data.empty()) throw std::runtime_error("Empty matrix.");
+    return *std::max_element(_data.begin(), _data.end());
+}
+
+template<typename T>
+const T& Matrix<T>::max() const {
+    if (_data.empty()) throw std::runtime_error("Empty matrix.");
+    return *std::max_element(_data.begin(), _data.end());
+}
+
+template<typename T>
+int Matrix<T>::size() const { return _rowCount*_columnCount; }
+
 template<typename T>
 int Matrix<T>::rowCount() const { return _rowCount; }
 

Source/Core/Matrix.hpp

@@ -33,10 +33,17 @@ class Matrix {
     T& at(int k);
     const T& at(int k) const;
 
+    /* Returns the maximum value in the matrix. */
+    T& max();
+    const T& max() const;
+
     /* Access to the matrix value at the i-th row and j-th column. */
     T& at(int i, int j);
     const T& at(int i, int j) const;
 
+    /* Total number of elements. */
+    int size() const;
+
     /* The number of rows. */
     int rowCount() const;
 

Source/Core/Utility.cpp

@@ -1,6 +1,7 @@
 #include <tuple>
 #include <cmath>
 #include <limits>
+#include <numbers>
 #include <stdexcept>
 #include "Core/Utility.hpp"
 
@@ -161,3 +162,57 @@ Matrix<double> Utility::executeMinerSummation(double category, double resistance
     return results;
 
 }
+
+std::pair<double, QString> Utility::computeRemainingFatigueLife(const QString& detail,
+    const std::unordered_map<QString, double>& parameters, const Matrix<double>& history,
+    const Matrix<double>& rainflow, double C, double m, double ΔKth, double ΔKcr, double a0) {
+    double sampledTimePeriod = (history.at(history.rowCount() - 1, 0) - history.at(0, 0))/(60*60*24*365.2425); // years
+    auto c = rainflow.getColumn(0);
+    auto Δσ = rainflow.getColumn(1);
+    if (detail == "IIW - Cruciform joint K-butt weld") {
+        double w = parameters.at("w"); // Out-of-plane width         [mm]
+        double H = parameters.at("H"); // Fillet weld leg length     [mm]
+        double T = parameters.at("T"); // Attachment plate thickness [mm]
+        double t = parameters.at("t"); // Main plate thickness       [mm]
+        double W = parameters.at("W"); // Fillet weld leg length     [mm]
+        double A = parameters.at("A"); // Weld throat size           [mm]
+        double θ = parameters.at("θ"); // Weld angle                 [deg]
+        int i = 0;
+        int N = 0;
+        Matrix<double> ΔK(Δσ.size(), 1);
+        double a = a0;
+        double af = T;
+        double rfl = 0.0; // years
+        double dfl = 250.0; // years
+        double π = std::numbers::pi;
+        while (true) {
+            double limit1 = (2*a)/w;
+            if (a >= af) return { rfl, "a ≥ a_f" };
+            if (ΔK.at(i) >= ΔKcr) return {rfl, "ΔK ≥ ΔK_cr"};
+            if (limit1 >= 0.8) return { rfl, "(2*a)/w ≥ 0.8" };
+            if (rfl >= dfl) return { dfl, "RFL ≥ 250 years" };
+            double fw = std::sqrt(1.0/std::cos(π*a/w*std::sqrt(a/T)));
+            double c1 =  0.7061 - 0.4091*H/T + 0.1596*std::pow(H/T, 2) + 0.3739*W/T - 0.1329*std::pow(W/T, 2);
+            double k1 = -0.2434 - 0.3939*H/T + 0.1536*std::pow(H/T, 2) + 0.3004*W/T - 0.0995*std::pow(W/T, 2);
+            double Mk = c1*std::pow(a/T, k1);
+            double Mm = (1.04 + 0.202*std::pow(a/T, 2) - 0.106*std::pow(a/T, 4))/2.464;
+            double Mb = (1.0 - 1.34*a/T - 0.03*std::pow(a/T, 2))*Mm;
+            double Y = fw*Mk*Mm;
+            ΔK.at(i) = Y * Δσ.at(i) * std::sqrt(π * a);
+            if (ΔK.at(i) >= ΔKth) {
+                double Δa = c.at(i)*C*std::pow(ΔK.at(i), m);
+                a += Δa;
+            }
+            if (i < Δσ.size() - 1) {
+                ++i;
+            }
+            else {
+                if (ΔK.max() < ΔKth) return { dfl, "ΔK < ΔK_th" };
+                ++N;
+                i = 0;
+            }
+            rfl = sampledTimePeriod*N;
+        }
+    }
+    throw std::logic_error("Not implemented.");
+}

Source/Core/Utility.hpp

@@ -1,7 +1,9 @@
-#pragma once
+#pragma once
 
 #include <vector>
+#include <unordered_map>
 #include <utility>
+#include <QtCore>
 #include "Core/Matrix.hpp"
 
 /* Static class of utilities. */
@@ -28,4 +30,9 @@ class Utility {
     static Matrix<double> executeMinerSummation(double category, double resistanceFactor, const Matrix<double>& slopes,
         const Matrix<double>& rainflow, double stressFactor);
 
+    /* Computes the remaining fatigue life of a LEFM detail. */
+    static std::pair<double, QString> computeRemainingFatigueLife(const QString& detail,
+        const std::unordered_map<QString, double>& parameters, const Matrix<double>& history,
+        const Matrix<double>& rainflow, double C, double m, double ΔKth, double ΔKcr, double a0);
+
 };

Source/Gui/HistoryTab.cpp

@@ -160,6 +160,8 @@ void HistoryTab::buildGui() {
 
 void HistoryTab::refreshGui() {
     if (_detail == nullptr) {
+        _timeGroupBox->setVisible(true);
+        _repsGroupBox->setVisible(true);
         _timeGroupBox->setEnabled(false);
         _timeUnitsBox->setEnabled(false);
         _repsBox->setEnabled(false);
@@ -172,6 +174,8 @@ void HistoryTab::refreshGui() {
         _sampleTable->setRowCount(0);
         _sampleTable->setColumnCount(0);
     } else {
+        _timeGroupBox->setVisible(_detail->approach() != Detail::Approach::Lefm);
+        _repsGroupBox->setVisible(_detail->approach() != Detail::Approach::Lefm);
         _timeGroupBox->setEnabled(true);
         _timeUnitsBox->setEnabled(!_detail->ignoreTime());
         _repsBox->setEnabled(true);

Source/Gui/MainWindow.cpp

@@ -134,6 +134,11 @@ void MainWindow::buildGui() {
     _lefmRainflowItem->setText("Rainflow Counting");
     _lefmNavBar->addItem(_lefmRainflowItem);
 
+    // remaining life item
+    _remainingLifeItem = new QListWidgetItem(_lefmNavBar);
+    _remainingLifeItem->setText("Remaining Fatigue Life");
+    _lefmNavBar->addItem(_remainingLifeItem);
+
     // main stack
     _mainStack = new QStackedWidget(_centralWidget);
     _mainStack->setProperty("parent", "true");
@@ -168,6 +173,10 @@ void MainWindow::buildGui() {
     _parisLawTab = new ParisLawTab(_mainStack);
     _mainStack->addWidget(_parisLawTab);
 
+    // remaining life tab
+    _remainingLifeTab = new RemainingLifeTab(_mainStack);
+    _mainStack->addWidget(_remainingLifeTab);
+
     // connections
     QObject::connect(_homeTab, &HomeTab::editDetailRequestReceived, this, &MainWindow::onEditDetailRequestReceived);
     QObject::connect(_snNavBar, &QListWidget::currentRowChanged, this, &MainWindow::onSnNavBarCurrentRowChanged);
@@ -309,5 +318,9 @@ void MainWindow::onLefmNavBarCurrentRowChanged() {
         _mainStack->setCurrentWidget(_rainflowTab);
         _rainflowTab->setDetail(_detail);
     }
+    else if (_lefmNavBar->currentItem() == _remainingLifeItem) {
+        _mainStack->setCurrentWidget(_remainingLifeTab);
+        _remainingLifeTab->setDetail(dynamic_cast<LefmDetail*>(_detail));
+    }
 
 }

Source/Gui/MainWindow.hpp

@@ -11,6 +11,7 @@
 #include "Gui/DamageTab.hpp"
 #include "Gui/DetailSelectionTab.hpp"
 #include "Gui/ParisLawTab.hpp"
+#include "Gui/RemainingLifeTab.hpp"
 
 /* The main window of the application. */
 class MainWindow : public QMainWindow {
@@ -43,6 +44,7 @@ class MainWindow : public QMainWindow {
     QListWidgetItem* _parisLawItem;
     QListWidgetItem* _lefmHistoryItem;
     QListWidgetItem* _lefmRainflowItem;
+    QListWidgetItem* _remainingLifeItem;
     QStackedWidget* _mainStack;
     HomeTab* _homeTab;
     SnCurveTab* _snCurveTab;
@@ -51,6 +53,7 @@ class MainWindow : public QMainWindow {
     DamageTab* _damageTab;
     DetailSelectionTab* _detailSelectionTab;
     ParisLawTab* _parisLawTab;
+    RemainingLifeTab* _remainingLifeTab;
     Module _module;
     Project* _project;
     Detail* _detail;