Merge pull request #405 from Blueprints-org/403-feature-request-add-formulas-from-chapter-59-from-nen_en_1992_1_1

403 feature request add formulas from chapter 5.9 from nen-en-1992-1-1

Author: egarciamendez

blueprints/codes/eurocode/nen_en_1992_1_1_c2_2011/chapter_5_structural_analysis/formula_5_40a.py

@@ -0,0 +1,63 @@
+"""Formula 5.40a from NEN-EN 1992-1-1+C2:2011: Chapter 5 - Structural Analysis."""
+
+from blueprints.codes.eurocode.nen_en_1992_1_1_c2_2011 import NEN_EN_1992_1_1_C2_2011
+from blueprints.codes.formula import Formula
+from blueprints.codes.latex_formula import LatexFormula
+from blueprints.type_alias import M
+from blueprints.validations import raise_if_less_or_equal_to_zero, raise_if_negative
+
+
+class Form5Dot40aCheckLateralInstability(Formula):
+    r"""Class representing formula 5.40a for checking lateral instability in persistent situations."""
+
+    label = "5.40a"
+    source_document = NEN_EN_1992_1_1_C2_2011
+
+    def __init__(
+        self,
+        l_0t: M,
+        b: M,
+        h: M,
+    ) -> None:
+        r"""Check the conditions to ignore second order effects in connection with lateral instability.
+
+        NEN-EN 1992-1-1+C2:2011 art.5.9(3) - Formula (5.40a)
+
+        Parameters
+        ----------
+        l_0t : M
+            [:math:`l_{0t}`] is the distance between torsional restraints [m].
+        b : M
+            [:math:`b`] is the width of compression flange [m].
+        h : M
+            [:math:`h`] is the total depth of beam in central part of l_0t [m].
+        """
+        super().__init__()
+        self.l_0t = l_0t
+        self.b = b
+        self.h = h
+
+    @staticmethod
+    def _evaluate(
+        l_0t: float,
+        b: float,
+        h: float,
+    ) -> bool:
+        """Evaluates the formula, for more information see the __init__ method."""
+        raise_if_less_or_equal_to_zero(b=b, h=h)
+        raise_if_negative(l_0t=l_0t)
+
+        return (l_0t / b <= 50 / (h / b) ** (1 / 3)) and (h / b <= 2.5)
+
+    def latex(self) -> LatexFormula:
+        """Returns LatexFormula object for formula 5.40a."""
+        return LatexFormula(
+            return_symbol=r"CHECK",
+            result="OK" if self.__bool__() else "\\text{Not OK}",
+            equation=r"\left( \frac{l_{0t}}{b} \leq \frac{50}{\left( h/b \right)^{1/3}} \text{ and } \frac{h}{b} \leq 2.5 \right)",
+            numeric_equation=rf"\left( \frac{{{self.l_0t:.3f}}}{{{self.b:.3f}}} \leq \frac{{50}}{{\left("
+            rf" {self.h:.3f}/{self.b:.3f} \right)^{{1/3}}}} \text{{ and }} "
+            rf"\frac{{{self.h:.3f}}}{{{self.b:.3f}}} \leq 2.5 \right)",
+            comparison_operator_label="\\to",
+            unit="",
+        )

blueprints/codes/eurocode/nen_en_1992_1_1_c2_2011/chapter_5_structural_analysis/formula_5_40b.py

@@ -0,0 +1,63 @@
+"""Formula 5.40b from NEN-EN 1992-1-1+C2:2011: Chapter 5 - Structural Analysis."""
+
+from blueprints.codes.eurocode.nen_en_1992_1_1_c2_2011 import NEN_EN_1992_1_1_C2_2011
+from blueprints.codes.formula import Formula
+from blueprints.codes.latex_formula import LatexFormula
+from blueprints.type_alias import M
+from blueprints.validations import raise_if_less_or_equal_to_zero, raise_if_negative
+
+
+class Form5Dot40bCheckLateralInstability(Formula):
+    r"""Class representing formula 5.40b for checking lateral instability in transient situations."""
+
+    label = "5.40b"
+    source_document = NEN_EN_1992_1_1_C2_2011
+
+    def __init__(
+        self,
+        l_0t: M,
+        b: M,
+        h: M,
+    ) -> None:
+        r"""Check the conditions to ignore second order effects in connection with lateral instability.
+
+        NEN-EN 1992-1-1+C2:2011 art.5.9(3) - Formula (5.40b)
+
+        Parameters
+        ----------
+        l_0t : M
+            [:math:`l_{0t}`] is the distance between torsional restraints [m].
+        b : M
+            [:math:`b`] is the width of compression flange [m].
+        h : M
+            [:math:`h`] is the total depth of beam in central part of l_0t [m].
+        """
+        super().__init__()
+        self.l_0t = l_0t
+        self.b = b
+        self.h = h
+
+    @staticmethod
+    def _evaluate(
+        l_0t: float,
+        b: float,
+        h: float,
+    ) -> bool:
+        """Evaluates the formula, for more information see the __init__ method."""
+        raise_if_less_or_equal_to_zero(b=b, h=h)
+        raise_if_negative(l_0t=l_0t)
+
+        return (l_0t / b <= 70 / (h / b) ** (1 / 3)) and (h / b <= 3.5)
+
+    def latex(self) -> LatexFormula:
+        """Returns LatexFormula object for formula 5.40b."""
+        return LatexFormula(
+            return_symbol=r"CHECK",
+            result="OK" if self.__bool__() else "\\text{Not OK}",
+            equation=r"\left( \frac{l_{0t}}{b} \leq \frac{70}{\left( h/b \right)^{1/3}} \text{ and } \frac{h}{b} \leq 3.5 \right)",
+            numeric_equation=rf"\left( \frac{{{self.l_0t:.3f}}}{{{self.b:.3f}}} \leq \frac{{70}}{{\left("
+            rf" {self.h:.3f}/{self.b:.3f} \right)^{{1/3}}}} \text{{ and }} "
+            rf"\frac{{{self.h:.3f}}}{{{self.b:.3f}}} \leq 3.5 \right)",
+            comparison_operator_label="\\to",
+            unit="",
+        )

docs/source/codes/eurocode/ec2_1992_1_1_2011/formulas.md

@@ -89,8 +89,8 @@ Total of 304 formulas present.
 | 5.38a          | :heavy_check_mark: |         | Form5Dot38aCheckRelativeSlendernessRatio                  |
 | 5.38b          | :heavy_check_mark: |         | Form5Dot38bCheckRelativeEccentricityRatio                 |
 | 5.39           | :heavy_check_mark: |         | Form5Dot39SimplifiedCriterionBiaxialBending               |
-| 5.40a          |        :x:         |         |                                                           |
-| 5.40b          |        :x:         |         |                                                           |
+| 5.40a          |        :x:         |         | Form5Dot40aCheckLateralInstability                        |
+| 5.40b          |        :x:         |         | Form5Dot40bCheckLateralInstability                        |
 | 5.41           |        :x:         |         |                                                           |
 | 5.42           |        :x:         |         |                                                           |
 | 5.43           |        :x:         |         |                                                           |

tests/codes/eurocode/nen_en_1992_1_1_c2_2011/chapter_5_structural_analysis/test_formula_5_40a.py

@@ -0,0 +1,69 @@
+"""Testing formula 5.40a of NEN-EN 1992-1-1+C2:2011."""
+
+import pytest
+
+from blueprints.codes.eurocode.nen_en_1992_1_1_c2_2011.chapter_5_structural_analysis.formula_5_40a import Form5Dot40aCheckLateralInstability
+from blueprints.validations import LessOrEqualToZeroError, NegativeValueError
+
+
+class TestForm5Dot40aCheckLateralInstability:
+    """Validation for formula 5.40a from NEN-EN 1992-1-1+C2:2011."""
+
+    def test_evaluation(self) -> None:
+        """Tests the evaluation of the result."""
+        # Example values
+        l_0t = 6.0  # m
+        b = 0.3  # m
+        h = 0.5  # m
+
+        # Object to test
+        formula = Form5Dot40aCheckLateralInstability(l_0t=l_0t, b=b, h=h)
+
+        # Expected result, manually calculated
+        expected_result = True
+
+        assert formula == expected_result
+
+    @pytest.mark.parametrize(
+        ("l_0t", "b", "h"),
+        [
+            (-6.0, 0.3, 0.5),  # l_0t is negative
+            (6.0, -0.3, 0.5),  # b is negative
+            (6.0, 0.3, -0.5),  # h is negative
+            (6.0, 0.0, 0.5),  # b is zero
+            (6.0, 0.3, 0.0),  # h is zero
+        ],
+    )
+    def test_raise_error_when_invalid_values_are_given(self, l_0t: float, b: float, h: float) -> None:
+        """Test invalid values."""
+        with pytest.raises((NegativeValueError, LessOrEqualToZeroError)):
+            Form5Dot40aCheckLateralInstability(l_0t=l_0t, b=b, h=h)
+
+    @pytest.mark.parametrize(
+        ("representation", "expected"),
+        [
+            (
+                "complete",
+                r"CHECK \to \left( \frac{l_{0t}}{b} \leq \frac{50}{\left( h/b \right)^{1/3}} \text{ and } \frac{h}{b} \leq 2.5 \right) \to "
+                r"\left( \frac{6.000}{0.300} \leq \frac{50}{\left( 0.500/0.300 \right)^{1/3}} \text{ and } "
+                r"\frac{0.500}{0.300} \leq 2.5 \right) \to OK",
+            ),
+            ("short", r"CHECK \to OK"),
+        ],
+    )
+    def test_latex(self, representation: str, expected: str) -> None:
+        """Test the latex representation of the formula."""
+        # Example values
+        l_0t = 6.0  # m
+        b = 0.3  # m
+        h = 0.5  # m
+
+        # Object to test
+        latex = Form5Dot40aCheckLateralInstability(l_0t=l_0t, b=b, h=h).latex()
+
+        actual = {
+            "complete": latex.complete,
+            "short": latex.short,
+        }
+
+        assert expected == actual[representation], f"{representation} representation failed."

tests/codes/eurocode/nen_en_1992_1_1_c2_2011/chapter_5_structural_analysis/test_formula_5_40b.py

@@ -0,0 +1,69 @@
+"""Testing formula 5.40b of NEN-EN 1992-1-1+C2:2011."""
+
+import pytest
+
+from blueprints.codes.eurocode.nen_en_1992_1_1_c2_2011.chapter_5_structural_analysis.formula_5_40b import Form5Dot40bCheckLateralInstability
+from blueprints.validations import LessOrEqualToZeroError, NegativeValueError
+
+
+class TestForm5Dot40bCheckLateralInstability:
+    """Validation for formula 5.40b from NEN-EN 1992-1-1+C2:2011."""
+
+    def test_evaluation(self) -> None:
+        """Tests the evaluation of the result."""
+        # Example values
+        l_0t = 6.0  # m
+        b = 0.3  # m
+        h = 0.5  # m
+
+        # Object to test
+        formula = Form5Dot40bCheckLateralInstability(l_0t=l_0t, b=b, h=h)
+
+        # Expected result, manually calculated
+        expected_result = True
+
+        assert formula == expected_result
+
+    @pytest.mark.parametrize(
+        ("l_0t", "b", "h"),
+        [
+            (-6.0, 0.3, 0.5),  # l_0t is negative
+            (6.0, -0.3, 0.5),  # b is negative
+            (6.0, 0.3, -0.5),  # h is negative
+            (6.0, 0.0, 0.5),  # b is zero
+            (6.0, 0.3, 0.0),  # h is zero
+        ],
+    )
+    def test_raise_error_when_invalid_values_are_given(self, l_0t: float, b: float, h: float) -> None:
+        """Test invalid values."""
+        with pytest.raises((NegativeValueError, LessOrEqualToZeroError)):
+            Form5Dot40bCheckLateralInstability(l_0t=l_0t, b=b, h=h)
+
+    @pytest.mark.parametrize(
+        ("representation", "expected"),
+        [
+            (
+                "complete",
+                r"CHECK \to \left( \frac{l_{0t}}{b} \leq \frac{70}{\left( h/b \right)^{1/3}} \text{ and } \frac{h}{b} \leq 3.5 \right) \to "
+                r"\left( \frac{6.000}{0.300} \leq \frac{70}{\left( 0.500/0.300 \right)^{1/3}} \text{ and } "
+                r"\frac{0.500}{0.300} \leq 3.5 \right) \to OK",
+            ),
+            ("short", r"CHECK \to OK"),
+        ],
+    )
+    def test_latex(self, representation: str, expected: str) -> None:
+        """Test the latex representation of the formula."""
+        # Example values
+        l_0t = 6.0  # m
+        b = 0.3  # m
+        h = 0.5  # m
+
+        # Object to test
+        latex = Form5Dot40bCheckLateralInstability(l_0t=l_0t, b=b, h=h).latex()
+
+        actual = {
+            "complete": latex.complete,
+            "short": latex.short,
+        }
+
+        assert expected == actual[representation], f"{representation} representation failed."