structuralcodes 0.0.1__py3-none-any.whl

structuralcodes/__init__.py

@@ -0,0 +1,17 @@
+"""A Python package that contains models from structural design codes."""
+
+from . import codes, core, geometry, materials, sections
+from .codes import get_design_codes, set_design_code, set_national_annex
+
+__version__ = '0.0.1'
+
+__all__ = [
+    'set_design_code',
+    'get_design_codes',
+    'set_national_annex',
+    'codes',
+    'core',
+    'materials',
+    'geometry',
+    'sections',
+]

structuralcodes/codes/__init__.py

@@ -0,0 +1,79 @@
+"""Collection of functions related to design codes."""
+
+import types
+import typing as t
+
+from . import ec2_2004, ec2_2023, mc2010
+
+__all__ = [
+    'mc2010',
+    'ec2_2023',
+    'ec2_2004',
+    'set_design_code',
+    'get_design_codes',
+    'set_national_annex',
+]
+
+# Global code object used by material classes
+_CODE: t.Optional[types.ModuleType] = None
+
+# Global national annex object
+_NATIONAL_ANNEX: t.Optional[str] = None
+
+# Design code registry
+_DESIGN_CODES = {
+    'mc2010': mc2010,
+    'ec2_2004': ec2_2004,
+    'ec2_2023': ec2_2023,
+}
+
+
+def set_design_code(design_code: t.Optional[str] = None) -> None:
+    """Set the current design code globally.
+
+    Args:
+        design_code (str): The abbreviation of the code.
+
+    Note:
+        Call get_design_codes() to get a list of the available codes.
+    """
+    global _CODE  # pylint: disable=W0603
+    if design_code is not None:
+        _CODE = _DESIGN_CODES.get(design_code.lower())
+    else:
+        _CODE = None
+
+
+def get_design_codes() -> t.List[str]:
+    """Get a list of the available design codes."""
+    return list(_DESIGN_CODES.keys())
+
+
+def set_national_annex(national_annex: str) -> None:
+    """Set the current national annex globally.
+
+    Args:
+        national_annex (str): The abbreviation of the national annex.
+
+    Note:
+        Call get_national_annexes() on the relevant design code to see a list
+        of available national annexes.
+    """
+    global _NATIONAL_ANNEX  # pylint: disable=W0603
+    _NATIONAL_ANNEX = national_annex.lower()
+
+
+def _use_design_code(
+    design_code: t.Optional[str] = None,
+) -> t.Optional[types.ModuleType]:
+    """Use a design code in a class.
+
+    Kwargs:
+        design_code (str): The abbreviation of the code.
+
+    Note:
+        Call get_design_codes() to get a list of the available codes.
+    """
+    if design_code is None:
+        return _CODE
+    return _DESIGN_CODES.get(design_code.lower())

structuralcodes/codes/ec2_2004/__init__.py

@@ -0,0 +1,133 @@
+"""EUROCODE 2 1992-1-1:2004."""
+
+import typing as t
+
+from ._concrete_material_properties import (
+    Ecm,
+    eps_c1,
+    eps_c2,
+    eps_c3,
+    eps_cu1,
+    eps_cu2,
+    eps_cu3,
+    fcd,
+    fcm,
+    fctk_5,
+    fctk_95,
+    fctm,
+    k_sargin,
+    n_parabolic_rectangular,
+)
+from ._reinforcement_material_properties import (
+    epsud,
+    fyd,
+    reinforcement_duct_props,
+)
+from ._section_7_3_crack_control import (
+    As_min,
+    As_min_2,
+    As_min_p,
+    alpha_e,
+    eps_sm_eps_cm,
+    hc_eff,
+    k,
+    k1,
+    k2,
+    k3,
+    k4,
+    kc_flanges_area,
+    kc_rect_area,
+    kc_tension,
+    kt,
+    phi_eq,
+    rho_p_eff,
+    sr_max_close,
+    sr_max_far,
+    sr_max_theta,
+    w_max,
+    w_spacing,
+    wk,
+    xi1,
+)
+from .annex_b_shrink_and_creep import (
+    beta_c,
+    beta_fcm,
+    beta_H,
+    beta_RH,
+    eps_cd_0,
+    eps_cs,
+    phi,
+    phi_RH,
+    t0_adj,
+)
+from .shear import (
+    Asw_max,
+    VEdmax_unreinf,
+    VRdc,
+    VRdc_prin_stress,
+    VRdmax,
+    VRds,
+)
+
+__all__ = [
+    'As_min',
+    'As_min_2',
+    'As_min_p',
+    'Asw_max',
+    'alpha_e',
+    'eps_sm_eps_cm',
+    'hc_eff',
+    'k',
+    'k1',
+    'k2',
+    'k3',
+    'k4',
+    'kc_flanges_area',
+    'kc_rect_area',
+    'kc_tension',
+    'kt',
+    'phi_eq',
+    'reinforcement_duct_props',
+    'rho_p_eff',
+    'sr_max_close',
+    'sr_max_far',
+    'sr_max_theta',
+    'w_max',
+    'w_spacing',
+    'wk',
+    'xi1',
+    'fcd',
+    'fcm',
+    'fctm',
+    'fctk_5',
+    'fctk_95',
+    'Ecm',
+    'eps_c1',
+    'eps_cu1',
+    'k_sargin',
+    'eps_c2',
+    'eps_cu2',
+    'n_parabolic_rectangular',
+    'eps_c3',
+    'eps_cu3',
+    'fyd',
+    'epsud',
+    'VEdmax_unreinf',
+    'VRdc',
+    'VRdc_prin_stress',
+    'VRdmax',
+    'VRds',
+    'beta_c',
+    'beta_fcm',
+    'beta_H',
+    'beta_RH',
+    'eps_cd_0',
+    'eps_cs',
+    'phi',
+    'phi_RH',
+    't0_adj',
+]
+
+__title__: str = 'EUROCODE 2 1992-1-1:2004'
+__year__: str = '2004'
+__materials__: t.Tuple[str] = ('concrete', 'reinforcement')

structuralcodes/codes/ec2_2004/_concrete_material_properties.py

@@ -0,0 +1,239 @@
+"""Concrete material properties according to Tab. 3.1."""
+
+import math
+
+from structuralcodes.codes import mc2010
+
+
+def fcm(fck: float, delta_f: float = 8) -> float:
+    """The mean compressive strength of concrete.
+
+    EN 1992-1-1:2004, Table 3.1.
+
+    Args:
+        fck (float): The characteristic compressive strength of concrete in
+            MPa.
+
+    Keyword Args:
+        delta_f (float): The difference between the mean and the
+            characteristic strength.
+
+    Returns:
+        float: The mean compressive strength in MPa.
+    """
+    return mc2010.fcm(fck=abs(fck), delta_f=abs(delta_f))
+
+
+def fctm(fck: float) -> float:
+    """The mean tensile strength of concrete.
+
+    EN 1992-1-1: 2004, Table 3.1.
+
+    Args:
+        fck (float): The characteristic compressive strength of concrete in
+            MPa.
+
+    Returns:
+        float: The mean tensile strength in MPa.
+    """
+    return mc2010.fctm(fck=abs(fck))
+
+
+def fctk_5(fctm: float) -> float:
+    """The 5% fractile of the tensile strength of concrete.
+
+    EN 1992-1-1: 2004, Table 3.1.
+
+    Args:
+        fctm (float): The mean tensile strength of concrete in MPa.
+
+    Returns:
+        float: The 5% fractile of the tensile strength in MPa.
+    """
+    return mc2010.fctkmin(fctm=abs(fctm))
+
+
+def fctk_95(fctm: float) -> float:
+    """The 95% fractile of the tensile strength of concrete.
+
+    EN 1992-1-1: 2004, Table 3.1.
+
+    Args:
+        fctm (float): The mean tensile strength of concrete in MPa.
+
+    Returns:
+        float: The 95% fractile of the tensile strength in MPa.
+    """
+    return mc2010.fctkmax(fctm=abs(fctm))
+
+
+def Ecm(fcm: float) -> float:
+    """The secant modulus of concrete.
+
+    EN 1992-1-1:2004, Table 3.1.
+
+    Args:
+        fcm (float): The mean compressive strength of concrete in MPa.
+
+    Returns:
+        float: The secant modulus of concrete in MPa.
+    """
+    return 22000.0 * math.pow(abs(fcm) / 10, 0.3)
+
+
+def eps_c1(fcm: float) -> float:
+    """The strain at maximum compressive stress of concrete (fcm) for the
+    Sargin constitutive law.
+
+    EN 1992-1-1:2004, Table 3.1.
+
+    Args:
+        fcm (float): The mean compressive strength of concrete in MPa.
+
+    Returns:
+        float: The strain at maximum compressive stress, absolute value, no
+        unit.
+    """
+    return min(0.7 * math.pow(abs(fcm), 0.31), 2.8) / 1000
+
+
+def eps_cu1(fck: float) -> float:
+    """The ultimate strain for the Sargin constitutive law.
+
+    EN 1992-1-1:2004, Table 3.1.
+
+    Args:
+        fck (float): The characteristic compressive strength of concrete in
+            MPa.
+
+    Returns:
+        float: The ultimate strain, absolute value, no unit.
+    """
+    fck = abs(fck)
+    return (
+        3.5 / 1000
+        if fck < 50
+        else (2.8 + 27 * ((98 - fcm(fck)) / 100) ** 4) / 1000
+    )
+
+
+def k_sargin(
+    Ecm: float,
+    fcm: float,
+    eps_c1: float,
+) -> float:
+    """Computation of k parameter for Sargin constitutive Law.
+
+    EN 1992-1-1:2004, Eq. (3.14)
+
+    Args:
+        Ecm (float): the mean elastic modulus of concrete in MPa.
+        fcm (float): the mean compressive strength in MPa.
+        eps_c1 (float): the strain corresponding to peak stress.
+    """
+    return 1.05 * Ecm * abs(eps_c1) / fcm
+
+
+def eps_c2(fck: float) -> float:
+    """The strain at maximum compressive stress of concrete for the
+    parabolic-rectangular law.
+
+    EN 1992-1-1:2004, Table 3.1.
+
+    Args:
+        fck (float): The characteristic compressive strength of concrete in
+            MPa.
+
+    Returns:
+        float: The strain at maximum compressive stress, absolute value, no
+        unit.
+    """
+    fck = abs(fck)
+    return (
+        2.0 / 1000 if fck <= 50 else (2.0 + 0.085 * (fck - 50) ** 0.53) / 1000
+    )
+
+
+def eps_cu2(fck: float) -> float:
+    """The ultimate strain of the parabolic-rectangular law.
+
+    EN 1992-1-1:2004, Table 3.1.
+
+    Args:
+        fck (float): The characteristic compressive strength of concrete in
+            MPa.
+
+    Returns:
+        float: The ultimate strain, absolute value, no unit.
+    """
+    fck = abs(fck)
+    return (
+        3.5 / 1000
+        if fck <= 50
+        else (2.6 + 35 * ((90 - fck) / 100) ** 4) / 1000
+    )
+
+
+def n_parabolic_rectangular(fck: float) -> float:
+    """The exponent in the parabolic-rectangular law.
+
+    EN 1992-1-1:2004, Table 3.1.
+
+    Args:
+        fck (float): The characteristic compressive strength of concrete in
+            MPa.
+
+    Returns:
+        float: The exponent n, absolute value, no unit.
+    """
+    fck = abs(fck)
+    return 2.0 if fck <= 50 else (1.4 + 23.4 * ((90 - fck) / 100) ** 4)
+
+
+def eps_c3(fck: float) -> float:
+    """The strain at maximum compressive stress of the bi-linear law.
+
+    EN 1992-1-1:2004, Table 3.1.
+
+    Args:
+        fck (float): The characteristic compressive strength of concrete in
+            MPa.
+
+    Returns:
+        float: The strain at maximum compressive stress, absolute value, no
+        unit.
+    """
+    fck = abs(fck)
+    return 1.75 / 1000 if fck <= 50 else (1.75 + 0.55 * (fck - 50) / 40) / 1000
+
+
+def eps_cu3(fck: float) -> float:
+    """The ultimate strain of the bi-linear law.
+
+    EN 1992-1-1:2004, Table 3.1.
+
+    Args:
+        fck (float): The characteristic compressive strength of concrete in
+            MPa.
+
+    Returns:
+        float: The ultimate strain, absolute value, no unit.
+    """
+    return eps_cu2(fck)
+
+
+def fcd(fck: float, alpha_cc: float, gamma_c: float) -> float:
+    """The design compressive strength of concrete.
+
+    EN 1992-1-1:2004, Eq. (3.15)
+
+    Args:
+        fck (float): The characteristic compressive strength in MPa.
+        alpha_cc (float): A factor for considering long-term effects on the
+            strength, and effects that arise from the way the load is applied.
+        gamma_c (float): The partial factor of concrete.
+
+    Returns:
+        float: The design compressive strength of concrete in MPa
+    """
+    return abs(alpha_cc) * abs(fck) / abs(gamma_c)

structuralcodes/codes/ec2_2004/_reinforcement_material_properties.py

@@ -0,0 +1,104 @@
+"""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,
+    },
+}
+
+
+def fyd(fyk: float, gamma_s: float) -> float:
+    """Calculate the design value of the reinforcement yield strength.
+
+    EUROCODE 2 1992-1-1:2004, Fig. 3.8
+
+    Args:
+        fyk (float): The characteristic yield strength in MPa.
+        gamma_s (float): The partial factor.
+
+    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.
+
+    EUROCDE 2 1992-1-1:2004, Fig 3.8
+
+    Args:
+        epsuk (float): The characteristic ultimate strain
+
+    Keyword Args:
+        gamma_eps (float): The partial factor specified in NA.
+            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'],
+) -> t.Dict[str, float]:
+    """Return a dict with the minimum characteristic ductility properties for
+    reinforcement ductility class.
+
+    EUROCODE 2 1992-1-1:2004, Tab. C.1
+
+    Args:
+        fyk (float): The characteristic yield strength.
+        ductility_class (Literal['A', 'B', 'C']): 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,
+    }