structuralcodes 0.1.0__tar.gz → 0.1.1__tar.gz

{structuralcodes-0.1.0 → structuralcodes-0.1.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: structuralcodes
-Version: 0.1.0
+Version: 0.1.1
 Summary: A Python package that contains models from structural design codes.
 Author-email: fib - International Federation for Structural Concrete <info@fib-international.org>
 Requires-Python: >=3.8

{structuralcodes-0.1.0 → structuralcodes-0.1.1}/structuralcodes/__init__.py

@@ -3,7 +3,7 @@
 from . import codes, core, geometry, materials, sections
 from .codes import get_design_codes, set_design_code, set_national_annex
 
-__version__ = '0.1.0'
+__version__ = '0.1.1'
 
 __all__ = [
     'set_design_code',

{structuralcodes-0.1.0 → structuralcodes-0.1.1}/structuralcodes/codes/ec2_2004/annex_b_shrink_and_creep.py

@@ -191,7 +191,7 @@ def eps_cs(
         ValueError: Checks if the cement class equals R, N or S.
     """
     _cement_class = cement_class.upper().strip()
-    beta_ds = (t - t_S) / (t - t_S + 0.04 * h_0 ** (1 / 3))  # (3.10)
+    beta_ds = (t - t_S) / (t - t_S + 0.04 * h_0 ** (3 / 2))  # (3.10)
     beta_as = 1 - math.exp(-0.2 * t**0.5)  # (3.13)
 
     # k_h is defined in Table 3.3 under (3.9)

{structuralcodes-0.1.0 → structuralcodes-0.1.1}/structuralcodes/codes/ec2_2004/shear.py

@@ -1,6 +1,7 @@
 """Design rules according to EN 1992-1-1 regarding shear."""
 
 import math
+import typing as t
 
 # General functions
 
@@ -171,6 +172,7 @@ def VRdc(
     fcd: float,
     k1: float = 0.15,
     gamma_c: float = 1.5,
+    CRdc: t.Optional[float] = None,
 ) -> float:
     """Compute the design strength of the shear resistance.
 
@@ -194,11 +196,13 @@ def VRdc(
             value might differ between National Annexes.
         gamma_c (float): Partial factor for concrete. Default value = 1.5,
             value might differ between National Annexes.
+        CRdc (Optional[float]): Scaling factor for the shear resistance.
+            Default value is 0.18 / gamma_c.
 
     Returns:
         float: The concrete shear resistance in MPa.
     """
-    CRdc = 0.18 / gamma_c
+    CRdc = CRdc or 0.18 / gamma_c
     return (
         max(
             CRdc * _k(d) * (100 * _rho_L(Asl, bw, d) * fck) ** (1.0 / 3.0)

{structuralcodes-0.1.0 → structuralcodes-0.1.1}/structuralcodes/codes/mc2010/_concrete_creep_and_shrinkage.py

@@ -291,7 +291,9 @@ def beta_ds(
         numpy.ndarray: Multiplication factor used for calculating the drying
         shrinkage as a function of time.
     """
-    return np.sqrt((time - ts) / (0.035 * (notional_size) ** 2 + (time - ts)))
+    time_drying = np.atleast_1d(time - ts)
+    time_drying[time_drying < 0.0] = 0.0
+    return np.sqrt(time_drying / (0.035 * (notional_size) ** 2 + time_drying))
 
 
 def beta_s1(fcm: float) -> float:

{structuralcodes-0.1.0 → structuralcodes-0.1.1}/structuralcodes/geometry/_reinforcement.py

@@ -18,6 +18,7 @@ def add_reinforcement(
     coords: t.Tuple[float, float],
     diameter: float,
     material: t.Union[Material, ConstitutiveLaw],
+    group_label: t.Optional[str] = None,
 ) -> CompoundGeometry:
     """Add a single bar given coordinate.
 
@@ -28,12 +29,16 @@ def add_reinforcement(
         diameter (float): The diameter of the reinforcement.
         material (Union(Material, ConstitutiveLaw)): A material or a
             constitutive law for the behavior of the reinforcement.
+        group_label (Optional(str)): A label for grouping several objects
+            (default is None).
 
     Returns:
         CompoundGeometry: A compound geometry with the original geometry and
         the reinforcement.
     """
-    bar = PointGeometry(Point(coords), diameter, material)
+    bar = PointGeometry(
+        Point(coords), diameter, material, group_label=group_label
+    )
     return geo + bar
 
 
@@ -47,6 +52,7 @@ def add_reinforcement_line(
     s: float = 0.0,
     first: bool = True,
     last: bool = True,
+    group_label: t.Optional[str] = None,
 ) -> CompoundGeometry:
     """Adds a set of bars distributed in a line.
 
@@ -66,6 +72,8 @@ def add_reinforcement_line(
             True).
         last (bool): Boolean indicating if placing the last bar (default =
             True).
+        group_label (Optional(str)): A label for grouping several objects
+            (default is None).
 
     Note:
         At least n or s should be greater than zero.
@@ -110,6 +118,10 @@ def add_reinforcement_line(
             continue
         coords = p1 + v * s * i
         geo = add_reinforcement(
-            geo, (coords[0], coords[1]), diameter, material
+            geo,
+            (coords[0], coords[1]),
+            diameter,
+            material,
+            group_label=group_label,
         )
     return geo

{structuralcodes-0.1.0 → structuralcodes-0.1.1}/structuralcodes/sections/_generic.py

@@ -2,9 +2,9 @@
 
 from __future__ import annotations  # To have clean hints of ArrayLike in docs
 
+import math
 import typing as t
 import warnings
-from math import cos, sin
 
 import numpy as np
 from numpy.typing import ArrayLike
@@ -208,7 +208,13 @@ class GenericSectionCalculator(SectionCalculator):
             )
             # Change sign due to moment sign convention
             izz *= -1
-            if abs(abs(izy) - abs(iyz)) > 10:
+
+            # Compute reasonable value for absolute tolerance for checking iyz
+            rel_tol = 1e-9
+            abs_tol = 0.5 * (iyy + izz) * rel_tol
+
+            # Check calculated cross moment
+            if not math.isclose(iyz, izy, rel_tol=rel_tol, abs_tol=abs_tol):
                 error_str = 'Something went wrong with computation of '
                 error_str += f'moments of area: iyz = {iyz}, izy = {izy}.\n'
                 error_str += 'They should be equal but are not!'
@@ -574,7 +580,12 @@ class GenericSectionCalculator(SectionCalculator):
         """Rotate triangulated data of angle theta."""
         rotated_triangulated_data = []
         for tr in self.triangulated_data:
-            T = np.array([[cos(theta), -sin(theta)], [sin(theta), cos(theta)]])
+            T = np.array(
+                [
+                    [np.cos(theta), -np.sin(theta)],
+                    [np.sin(theta), np.cos(theta)],
+                ]
+            )
             coords = np.vstack((tr[0], tr[1]))
             coords_r = T @ coords
             rotated_triangulated_data.append(
@@ -636,7 +647,9 @@ class GenericSectionCalculator(SectionCalculator):
         )
 
         # Rotate back to section CRS TODO Check
-        T = np.array([[cos(theta), -sin(theta)], [sin(theta), cos(theta)]])
+        T = np.array(
+            [[np.cos(theta), -np.sin(theta)], [np.sin(theta), np.cos(theta)]]
+        )
         M = T @ np.array([[My], [Mz]])
         if self.triangulated_data is not None:
             # Rotate back also triangulated data!
@@ -771,7 +784,12 @@ class GenericSectionCalculator(SectionCalculator):
                 geo=rotated_geom, strain=strain, tri=self.triangulated_data
             )
             # Rotate back to section CRS
-            T = np.array([[cos(theta), -sin(theta)], [sin(theta), cos(theta)]])
+            T = np.array(
+                [
+                    [np.cos(theta), -np.sin(theta)],
+                    [np.sin(theta), np.cos(theta)],
+                ]
+            )
             M = T @ np.array([[My], [Mz]])
             eps_a[i] = strain[0]
             my[i] = M[0, 0]
@@ -1106,7 +1124,9 @@ class GenericSectionCalculator(SectionCalculator):
         eps_a = eps_n - kappa_y * y_n
 
         # rotate back components to work in section CRS
-        T = np.array([[cos(theta), -sin(theta)], [sin(theta), cos(theta)]])
+        T = np.array(
+            [[np.cos(theta), -np.sin(theta)], [np.sin(theta), np.cos(theta)]]
+        )
         components = np.vstack((kappa_y, np.zeros_like(kappa_y)))
         rotated_components = T @ components
         return np.column_stack((eps_a, rotated_components.T)), field_num