Merge pull request #524 from Blueprints-org/523-formula-516-from-nen-en-1993-52008

523 formula 516 from nen en 1993 52008

Author: GerjanDorgelo

blueprints/codes/eurocode/nen_en_1993_5_2008/chapter_5_ultimate_limit_states/formula_5_16.py

@@ -0,0 +1,61 @@
+"""Formula 5.16 from NEN-EN 1993-5:2008 Chapter 5 - Ultimate limit state."""
+
+from blueprints.codes.eurocode.nen_en_1993_5_2008 import NEN_EN_1993_5_2008
+from blueprints.codes.formula import Formula
+from blueprints.codes.latex_formula import LatexFormula
+from blueprints.type_alias import DIMENSIONLESS, MM2, MPA, N
+from blueprints.validations import raise_if_less_or_equal_to_zero, raise_if_negative
+
+
+class Form5Dot16PlasticDesignResistance(Formula):
+    """Class representing formula 5.16 for the calculation of the plastic design resistance of the cross-section, [$N_{pl,Rd}$]."""
+
+    label = "5.16"
+    source_document = NEN_EN_1993_5_2008
+
+    def __init__(
+        self,
+        a: MM2,
+        f_y: MPA,
+        gamma_m0: DIMENSIONLESS,
+    ) -> None:
+        r"""[$N_{pl,Rd}$] Plastic design resistance of the cross-section [$N$].
+
+        NEN-EN 1993-5:2008 art.5.2.3 (9) - Formula (5.16)
+
+        Parameters
+        ----------
+        a : MM2
+            [$A$] Area of the cross-section [$mm^2$].
+        f_y : MPA
+            [$f_y$] Yield strength [$MPa$].
+        gamma_m0 : DIMENSIONLESS
+            [$\gamma_{M0}$] Partial factor according to 5.1.1 (4) [$-$].
+        """
+        super().__init__()
+        self.a = a
+        self.f_y = f_y
+        self.gamma_m0 = gamma_m0
+
+    @staticmethod
+    def _evaluate(
+        a: MM2,
+        f_y: MPA,
+        gamma_m0: DIMENSIONLESS,
+    ) -> N:
+        """Evaluates the formula, for more information see the __init__ method."""
+        raise_if_less_or_equal_to_zero(gamma_m0=gamma_m0)
+        raise_if_negative(f_y=f_y, a=a)
+
+        return (a * f_y) / gamma_m0
+
+    def latex(self) -> LatexFormula:
+        """Returns LatexFormula object for formula 5.16."""
+        return LatexFormula(
+            return_symbol=r"N_{pl,Rd}",
+            result=f"{self:.3f}",
+            equation=r"\frac{A \cdot f_y}{\gamma_{M0}}",
+            numeric_equation=rf"\frac{{{self.a:.3f} \cdot {self.f_y:.3f}}}{{{self.gamma_m0:.3f}}}",
+            comparison_operator_label="=",
+            unit="N",
+        )

docs/objects_overview/eurocode/nen_en_1993_5_2008/formulas.md

@@ -24,7 +24,7 @@ Total of 63 formulas present.
 | 5.13           | :heavy_check_mark: |         | Form5Dot13SimplifiedBucklingCheck        |
 | 5.14           |        :x:         |         |                                          |
 | 5.15           |        :x:         |         |                                          |
-| 5.16           |        :x:         |         |                                          |
+| 5.16           | :heavy_check_mark: |         | Form5Dot16PlasticDesignResistance        |
 | 5.17           |        :x:         |         |                                          |
 | 5.18           |        :x:         |         |                                          |
 | 5.19           |        :x:         |         |                                          |

tests/codes/eurocode/nen_en_1993_5_2008/chapter_5_ultimate_limit_states/test_formula_5_16.py

@@ -0,0 +1,66 @@
+"""Testing formula 5.16 of NEN-EN 1993-5:2008."""
+
+import pytest
+
+from blueprints.codes.eurocode.nen_en_1993_5_2008.chapter_5_ultimate_limit_states.formula_5_16 import Form5Dot16PlasticDesignResistance
+from blueprints.validations import LessOrEqualToZeroError, NegativeValueError
+
+
+class TestForm5Dot16PlasticDesignResistance:
+    """Validation for formula 5.16 from NEN-EN 1993-5:2008."""
+
+    def test_evaluation(self) -> None:
+        """Test the evaluation of the result."""
+        # Example values
+        a = 5000  # MM2
+        f_y = 355  # MPA
+        gamma_m0 = 1.0  # DIMENSIONLESS
+
+        # Object to test
+        formula = Form5Dot16PlasticDesignResistance(a=a, f_y=f_y, gamma_m0=gamma_m0)
+
+        # Expected result, manually calculated
+        manually_calculated_result = 1775000.0  # N
+
+        assert formula == pytest.approx(manually_calculated_result, rel=1e-4)
+
+    @pytest.mark.parametrize(
+        ("a", "f_y", "gamma_m0"),
+        [
+            (-5000, 355, 1.0),  # a is negative
+            (5000, -355, 1.0),  # f_y is negative
+            (5000, 355, -1.0),  # gamma_m0 is negative
+            (5000, 355, 0),  # gamma_m0 is zero
+        ],
+    )
+    def test_raise_error_when_invalid_values_are_given(self, a: float, f_y: float, gamma_m0: float) -> None:
+        """Test invalid values."""
+        with pytest.raises((LessOrEqualToZeroError, NegativeValueError)):
+            Form5Dot16PlasticDesignResistance(a=a, f_y=f_y, gamma_m0=gamma_m0)
+
+    @pytest.mark.parametrize(
+        ("representation", "expected"),
+        [
+            (
+                "complete",
+                r"N_{pl,Rd} = \frac{A \cdot f_y}{\gamma_{M0}} = \frac{5000.000 \cdot 355.000}{1.000} = 1775000.000 N",
+            ),
+            ("short", r"N_{pl,Rd} = 1775000.000 N"),
+        ],
+    )
+    def test_latex(self, representation: str, expected: str) -> None:
+        """Test the latex representation of the formula."""
+        # Example values
+        a = 5000  # MM2
+        f_y = 355  # MPA
+        gamma_m0 = 1.0  # DIMENSIONLESS
+
+        # Object to test
+        latex = Form5Dot16PlasticDesignResistance(a=a, f_y=f_y, gamma_m0=gamma_m0).latex()
+
+        actual = {
+            "complete": latex.complete,
+            "short": latex.short,
+        }
+
+        assert expected == actual[representation], f"{representation} representation failed."