structuralcodes 0.0.1__py3-none-any.whl

structuralcodes/codes/mc2010/_concrete_torsion.py

@@ -0,0 +1,164 @@
+"""Covers torsion in Model code 2010, 7.3.4."""
+
+from math import cos, pi, sin, tan
+
+from ._concrete_shear import epsilon_x, eta_fc, v_rd_max
+
+
+def v_ed_ti(t_ed: float, a_k: float, z_i: float):
+    """Shear force due to torsion.
+
+    fib Model Code 2010, eq. (7.3-53).
+
+    Args:
+        t_ed: The acting torsion force in the cross section in Nmm.
+        z_i: Can be found in figure 7.3-18.
+        a_k: Can be found in figure 7.3-18.
+
+    Returns:
+        float: The shear force that will occur due to torsion moment.
+    """
+    return t_ed * z_i / (2 * a_k)
+
+
+def t_rd_max(
+    f_ck: float,
+    d_k: float,
+    a_k: float,
+    theta: float,
+    approx_lvl: int,
+    z: float,
+    E_s: float,
+    As: float,
+    loads: dict,
+    gamma_c: float = 1.5,
+) -> float:
+    """The maximum allowed torsion.
+
+    fib Model Code 2010, eq. (7.3-56).
+
+    Args:
+        f_ck (float): Characteristic strength in MPa.
+        gamma_c (float): Concrete safety factor.
+        d_k (float): Is the diameter in the smallest circel in the cross
+            section.
+        a_k: Can be found in figure 7.3-18.
+        theta (float): inclination of the compression stressfield in degrees.
+        approx_lvl (int): Approximation method for concrete with reinforcement.
+        z (float): distances between the centerline of the compressive chord
+            and the reinforcement in mm.
+        E_s (float): The E_modulus to the material in MPa.
+        As (float): The cross-section reinforcement in mm^2.
+        loads (dict): The given loads in a dictionary. See create_load_dict.
+        gamma_c (float): Concrete safety factor.
+
+    Returns:
+        float: The maximum allowed torsion.
+    """
+    t_ef = d_k / 8
+
+    if approx_lvl == 1:
+        k_epsilon = 0.55
+
+    elif approx_lvl == 2:
+        epsilonx = epsilon_x(E_s, As, z, loads)
+        epsilon_1 = epsilonx + (epsilonx + 0.002) * (
+            (1 / tan(theta * pi / 180)) ** 2
+        )
+        k_epsilon = min(1 / (1.2 + 55 * epsilon_1), 0.65)
+
+    elif approx_lvl == 3:
+        epsilonx = epsilon_x(E_s, As, z, loads)
+        theta_min = 20 + 10000 * epsilonx
+        epsilon_1 = epsilonx + (epsilonx + 0.002) * (
+            (1 / tan(theta_min * pi / 180)) ** 2
+        )
+        k_epsilon = min(1 / (1.2 + 55 * epsilon_1), 0.65)
+
+    k_c = eta_fc(f_ck) * k_epsilon
+    return (
+        k_c
+        * f_ck
+        * t_ef
+        * 2
+        * a_k
+        * sin(theta * pi / 180)
+        * cos(theta * pi / 180)
+        / gamma_c
+    )
+
+
+def t_rd(  # pylint: disable=r0801
+    t_ed: float,
+    approx_lvl: int,
+    fck: float,
+    bw: float,
+    theta: float,
+    z: float,
+    E_s: float,
+    As: float,
+    loads: dict,
+    d_k: float,
+    a_k: float,
+    alpha: float = 90.0,
+    gamma_c: float = 1.5,
+) -> bool:
+    """Checks if the combination of torsion and shear is ok.
+
+    fib Model Code 2010, eq. (7.3-56).
+
+    Args:
+        t_ed (float): The torsion working on the material in Nmm.
+        approx_lvl (int): Approximation level chosen for shear resistance.
+        fck (float): Characteristic strength in MPa.
+        bw: (float): Thickness of web in cross section.
+        theta (float): inclination of the compression stressfield in degrees.
+        z: (float): The length to the areasenter of cross-section in mm.
+        E_s: (float): The E_modulus to the material in MPa.
+        As: (float): The cross-section area in mm^2.
+        loads (dict): The given loads in a dictionary. See create_load_dict.
+        d_k (float): Is the diameter in the smallest circel in the cross
+            section.
+        a_k (float): Can be found in figure 7.3-18.
+        alpha (float): Inclination of the stirrups in degrees.
+        gamma_c (float): Concrete safety factor.
+
+    Returns:
+        bool: Returns a bool that is true if the criterion for torsion and
+        shear is fulfilled.
+    """
+    return bool(
+        (
+            t_ed
+            / t_rd_max(
+                fck,
+                d_k,
+                a_k,
+                theta,
+                approx_lvl,
+                z,
+                E_s,
+                As,
+                loads,
+                gamma_c,
+            )
+        )
+        ** 2
+        + (
+            loads.get('Ved')
+            / v_rd_max(
+                approx_lvl,
+                fck,
+                bw,
+                theta,
+                z,
+                E_s,
+                As,
+                loads,
+                alpha,
+                gamma_c,
+            )
+        )
+        ** 2
+        <= 1
+    )

structuralcodes/codes/mc2010/_reinforcement_material_properties.py

@@ -0,0 +1,105 @@
+"""Material properties for reinforcement steel."""
+
+import typing as t
+
+DUCTILITY_CLASSES = {
+    'A': {
+        'epsuk': 2.5e-2,
+        'k': 1.05,
+    },
+    'B': {
+        'epsuk': 5.0e-2,
+        'k': 1.08,
+    },
+    'C': {
+        'epsuk': 7.5e-2,
+        'k': 1.15,
+    },
+    'D': {
+        'epsuk': 8.0e-2,
+        'k': 1.25,
+    },
+}
+
+
+def fyd(fyk: float, gamma_s: float = 1.15) -> float:
+    """Calculate the design value of the reinforcement yield strength.
+
+    fib Model Code 2010, Sec. 4.5.2.2.3.
+
+    Args:
+        fyk (float): The characteristic yield strength in MPa.
+        gamma_s (float): The partial factor. Default value 1.15.
+
+    Returns:
+        float: The design yield strength in MPa.
+
+    Raises:
+        ValueError: If fyk is less than 0.
+        ValueError: If gamma_s is less than 1.
+    """
+    if fyk < 0:
+        raise ValueError(f'fyk={fyk} cannot be less than 0')
+    if gamma_s < 1:
+        raise ValueError(f'gamma_s={gamma_s} must be larger or equal to 1')
+    return fyk / gamma_s
+
+
+def epsud(epsuk: float, gamma_eps: float = 0.9) -> float:
+    """Calculate the design value of the reinforcement ultimate strain.
+
+    fib Model Code 2010, Sec. 7.2.3.2.
+
+    Args:
+        epsuk (float): The characteristic ultimate strain.
+        gamma_eps (float): The partial factor. Default value 0.9.
+
+    Returns:
+        float: The design ultimate strain.
+
+    Raises:
+        ValueError: If epsuk is less than 0.
+        ValueError: If gamma_eps is greater than 1.
+    """
+    if epsuk < 0:
+        raise ValueError(f'epsuk={epsuk} cannot be less than 0')
+    if gamma_eps > 1:
+        raise ValueError(
+            f'gamma_eps={gamma_eps} must be smaller or equal to 1'
+        )
+    return epsuk * gamma_eps
+
+
+def reinforcement_duct_props(
+    fyk: float,
+    ductility_class: t.Literal['A', 'B', 'C', 'D'],
+) -> t.Dict[str, float]:
+    """Return a dict with the minimum characteristic ductility properties for
+    reinforcement ductility class.
+
+    fib Model Code 2010, Sec. 5.2.5.4.
+
+    Args:
+        fyk (float): The characteristic yield strength.
+        ductility_class (Literal['A', 'B', 'C', 'D']): The reinforcement
+            ductility class designation.
+
+    Returns:
+        Dict[str, float]: A dict with the characteristik strain value at the
+        ultimate stress level (epsuk), and the characteristic ultimate stress
+        (ftk).
+
+    Raises:
+        ValueError: When the ductility_class does not define a valid ductility
+            class.
+    """
+    duct_props = DUCTILITY_CLASSES.get(ductility_class.upper(), None)
+    if duct_props is None:
+        raise ValueError(
+            'No properties were found for the provided ductility class '
+            f'({ductility_class}).'
+        )
+    return {
+        'epsuk': duct_props['epsuk'],
+        'ftk': duct_props['k'] * fyk,
+    }

structuralcodes/core/__init__.py

@@ -0,0 +1 @@
+"""Core functionality shared between other modules."""

structuralcodes/core/_section_results.py

@@ -0,0 +1,211 @@
+"""Results."""
+
+from __future__ import annotations  # To have clean hints of ArrayLike in docs
+
+from dataclasses import dataclass, field, fields
+
+from numpy.typing import ArrayLike
+
+
+@dataclass
+class GrossProperties:
+    """Simple dataclass for storing gross section properties."""
+
+    # section areas
+    area: float = field(default=0, metadata={'description': 'Total area'})
+    area_reinforcement: float = field(
+        default=0, metadata={'description': 'Reinforcement area'}
+    )
+
+    # axial rigidity
+    ea: float = field(
+        default=0, metadata={'description': 'Axial rigidity (E * A)'}
+    )
+    # ea = {'value': 0, 'description': 'Axial rigidity (E_A)'}
+
+    # section mass
+    mass: float = field(
+        default=0, metadata={'description': 'Mass per unit length'}
+    )
+
+    # section perimeter
+    perimeter: float = field(default=0, metadata={'description': 'Perimeter'})
+
+    # first moments of area
+    # global axes of section
+    sy: float = field(
+        default=0, metadata={'description': 'First moment of area (Sy)'}
+    )
+    sz: float = field(
+        default=0, metadata={'description': 'Second moment of area (Sz)'}
+    )
+
+    # first moments of area * E
+    # global axes of section
+    e_sy: float = field(default=0, metadata={'description': 'E * Sy'})
+    e_sz: float = field(default=0, metadata={'description': 'E * Sz'})
+
+    # centroids
+    cy: float = field(
+        default=0, metadata={'description': 'Centroid y Coordinate'}
+    )
+    cz: float = field(
+        default=0, metadata={'description': 'Centroid z Coordinate'}
+    )
+
+    # second moments of area
+    # global axes of section
+    iyy: float = field(
+        default=0, metadata={'description': 'Second moment (Iyy)'}
+    )
+    izz: float = field(
+        default=0, metadata={'description': 'Second moment (Izz)'}
+    )
+    iyz: float = field(
+        default=0, metadata={'description': 'Product moment (Iyz)'}
+    )
+    # centroidal axes of section
+    iyy_c: float = field(
+        default=0, metadata={'description': 'Centroidal Second moment (Iyy_c)'}
+    )
+    izz_c: float = field(
+        default=0, metadata={'description': 'Centroidal Second moment (Izz_c)'}
+    )
+    iyz_c: float = field(
+        default=0,
+        metadata={'description': 'Centroidal Product moment (Iyz_c)'},
+    )
+    # principal axes of section
+    i11: float = field(
+        default=0, metadata={'description': 'Principal Second Moment (I11)'}
+    )
+    i22: float = field(
+        default=0, metadata={'description': 'Principal Second Moment (I22)'}
+    )
+    theta: float = field(
+        default=0, metadata={'description': 'Principal axis angle (theta)'}
+    )
+
+    # section flexural rigidity
+    # global axes of section
+    e_iyy: float = field(default=0, metadata={'description': 'E * Iyy'})
+    e_izz: float = field(default=0, metadata={'description': 'E * Izz'})
+    e_iyz: float = field(default=0, metadata={'description': 'E * Iyz'})
+    # centroidal axes of section
+    e_iyy_c: float = field(default=0, metadata={'description': 'E * Iyy_c'})
+    e_izz_c: float = field(default=0, metadata={'description': 'E * Izz_c'})
+    e_iyz_c: float = field(default=0, metadata={'description': 'E * Iyz_c'})
+    # principal axes of section
+    e_i11: float = field(default=0, metadata={'description': 'E * I11'})
+    e_i22: float = field(default=0, metadata={'description': 'E * I22'})
+    e_theta: float = field(
+        default=0, metadata={'description': 'Principal axis angle (theta)'}
+    )
+
+    def __format__(self, spec: str) -> str:
+        """Defines the format for returning the string representation.
+
+        Arguments:
+            spec (str): The string specifying the format.
+        """
+        output_string = 'Gross Concrete Section Properties:\n'
+        for f in fields(self):
+            value = getattr(self, f.name)
+            description = f.metadata.get(
+                'description', 'No description available'
+            )
+            output_string += f'{description}: {value:{spec}}\n'
+
+        # etc. all other characteristics
+        return output_string
+
+    def __str__(self) -> str:
+        """Returns the informal string representation of the gross concrete
+        section properties.
+        """
+        return f'{self}'
+
+
+@dataclass
+class CrackedProperties:
+    """Simple dataclass for storing cracked section properties."""
+
+    # second moments of area
+    i_yy: float = 0
+    i_zz: float = 0
+    i_yz: float = 0
+
+    # section cracked flexural rigidity
+    ei_yy: float = 0
+    ei_zz: float = 0
+
+
+@dataclass
+class MomentCurvatureResults:
+    """Class for storing moment curvature results.
+
+    The analysis will be done in general for a given inclination of n.a.
+    """
+
+    theta: float = 0  # the inclination of n.a.
+    n: float = 0  # axial load - mantained constant during analysis
+    chi_y: ArrayLike = None  # the curvatures
+    chi_z: ArrayLike = None  # the curvatures
+    eps_axial: ArrayLike = 0  # the axial strain (at section 0,0)
+    m_y: ArrayLike = None  # the moment
+    m_z: ArrayLike = None  # the moment
+
+    # The strains can be reconstructed at each step from chi and eps_axial
+    # The stresses can be recomputed if needed on the fly? Or storing them?
+
+
+@dataclass
+class UltimateBendingMomentResults:
+    """Class for storing the ultimate bending moment computation for a given
+    inclination of n.a. and axial load.
+    """
+
+    theta: float = 0  # the inclination of n.a.
+    n: float = 0  # axial load - mantained constant during analysis
+    m_y: float = 0  # the ultimate moment for given theta and n
+    m_z: float = 0  # the ultimate moment for given theta and n
+    chi_y: float = 0  # the curvature corresponding to the ultimate moment
+    chi_z: float = 0  # the curvature corresponding to the ultimate moment
+    eps_a: float = 0  # the axial strain at 0,0 corresponding to Mult
+
+
+@dataclass
+class NMMInteractionDomain:
+    """Class for storing the NMM interaction domain results."""
+
+    num_theta: int = 0  # number of discretizations along the angle
+    num_axial: int = 0  # number of discretizations along axial load axis
+
+    strains: ArrayLike = None  # array with shape (n,3) containing strains
+    forces: ArrayLike = None  # array with shape(n,3) containing N, My, Mz
+
+
+@dataclass
+class NMInteractionDomain:
+    """Class for storing the NM interaction domain results."""
+
+    theta: float = 0  # the inclination of n.a.
+    num_axial: float = 0  # number of discretizations along axial load axis
+
+    n: ArrayLike = None  # Axial loads
+    m_y: ArrayLike = None  # Moments My
+    m_z: ArrayLike = None  # Moments Mz
+
+    strains: ArrayLike = None
+
+
+@dataclass
+class MMInteractionDomain:
+    """Class for storing the MM interaction domain results."""
+
+    num_theta: float = 0  # number of discretizations along the angle
+    n: float = 0  # axial load
+
+    theta: ArrayLike = None  # Angle theta respect axis Y
+    m_y: ArrayLike = None  # Moments My
+    m_z: ArrayLike = None  # Moments Mz