structuralcodes 0.5.0__py3-none-any.whl → 0.6.0__py3-none-any.whl

structuralcodes/geometry/profiles/_upn.py

@@ -0,0 +1,315 @@
+"""UPN profiles."""
+
+from shapely import (
+    Polygon,
+)
+
+from ._base_profile import BaseProfile
+from ._common_functions import (
+    _create_taper_U_section,
+)
+
+
+class UPN(BaseProfile):
+    """Simple class for representing an UPN profile.
+
+    European standard channels UPN 50 - 400.
+
+    Taper flange Channels.
+
+    14% slope in flange.
+    """
+
+    parameters = {
+        'UPN50': {
+            'h': 50.0,
+            'b': 38.0,
+            'tw': 5.0,
+            'tf': 7.0,
+            'r1': 7.0,
+            'r2': 4.0,
+            'd': 21.0,
+        },
+        'UPN65': {
+            'h': 65.0,
+            'b': 42.0,
+            'tw': 5.5,
+            'tf': 7.5,
+            'r1': 8.0,
+            'r2': 4.0,
+            'd': 34.0,
+        },
+        'UPN80': {
+            'h': 80.0,
+            'b': 45.0,
+            'tw': 6.0,
+            'tf': 8.0,
+            'r1': 8.0,
+            'r2': 4.0,
+            'd': 47.0,
+        },
+        'UPN100': {
+            'h': 100.0,
+            'b': 50.0,
+            'tw': 6.0,
+            'tf': 8.5,
+            'r1': 9.0,
+            'r2': 5.0,
+            'd': 64.0,
+        },
+        'UPN120': {
+            'h': 120.0,
+            'b': 55.0,
+            'tw': 7.0,
+            'tf': 9.0,
+            'r1': 9.0,
+            'r2': 5.0,
+            'd': 82.0,
+        },
+        'UPN140': {
+            'h': 140.0,
+            'b': 60.0,
+            'tw': 7.0,
+            'tf': 10.0,
+            'r1': 10.0,
+            'r2': 5.0,
+            'd': 98.0,
+        },
+        'UPN160': {
+            'h': 160.0,
+            'b': 65.0,
+            'tw': 7.5,
+            'tf': 10.5,
+            'r1': 11.0,
+            'r2': 6.0,
+            'd': 115.0,
+        },
+        'UPN180': {
+            'h': 180.0,
+            'b': 70.0,
+            'tw': 8.0,
+            'tf': 11.0,
+            'r1': 11.0,
+            'r2': 6.0,
+            'd': 133.0,
+        },
+        'UPN200': {
+            'h': 200.0,
+            'b': 75.0,
+            'tw': 8.5,
+            'tf': 11.5,
+            'r1': 12.0,
+            'r2': 6.0,
+            'd': 151.0,
+        },
+        'UPN220': {
+            'h': 220.0,
+            'b': 80.0,
+            'tw': 9.0,
+            'tf': 12.5,
+            'r1': 13.0,
+            'r2': 7.0,
+            'd': 167.0,
+        },
+        'UPN240': {
+            'h': 240.0,
+            'b': 85.0,
+            'tw': 9.5,
+            'tf': 13.0,
+            'r1': 13.0,
+            'r2': 7.0,
+            'd': 184.0,
+        },
+        'UPN260': {
+            'h': 260.0,
+            'b': 90.0,
+            'tw': 10.0,
+            'tf': 14.0,
+            'r1': 14.0,
+            'r2': 7.0,
+            'd': 200.0,
+        },
+        'UPN280': {
+            'h': 280.0,
+            'b': 95.0,
+            'tw': 10.0,
+            'tf': 15.0,
+            'r1': 15.0,
+            'r2': 8.0,
+            'd': 216.0,
+        },
+        'UPN300': {
+            'h': 300.0,
+            'b': 100.0,
+            'tw': 10.0,
+            'tf': 16.0,
+            'r1': 16.0,
+            'r2': 8.0,
+            'd': 232.0,
+        },
+        'UPN320': {
+            'h': 320.0,
+            'b': 100.0,
+            'tw': 14.0,
+            'tf': 17.5,
+            'r1': 18.0,
+            'r2': 9.0,
+            'd': 246.0,
+        },
+        'UPN350': {
+            'h': 350.0,
+            'b': 100.0,
+            'tw': 14.0,
+            'tf': 16.0,
+            'r1': 16.0,
+            'r2': 8.0,
+            'd': 282.0,
+        },
+        'UPN380': {
+            'h': 380.0,
+            'b': 102.0,
+            'tw': 13.5,
+            'tf': 16.0,
+            'r1': 16.0,
+            'r2': 8.0,
+            'd': 313.0,
+        },
+        'UPN400': {
+            'h': 400.0,
+            'b': 110.0,
+            'tw': 14.0,
+            'tf': 18.0,
+            'r1': 18.0,
+            'r2': 9.0,
+            'd': 324.0,
+        },
+    }
+
+    @classmethod
+    def get_polygon(cls, name: str) -> Polygon:
+        """Returns a shapely polygon representing an UPN section."""
+        if isinstance(name, (float, int)):
+            name = f'UPN{int(name):0d}'
+        parameters = cls.parameters.get(name)
+        if parameters is None:
+            raise ValueError(
+                f"Profile '{name}' not found in UPN sections. "
+                "Select a valid profile (available ones: "
+                f"{cls.profiles()})"
+            )
+        if parameters['h'] <= 300:
+            parameters['slope'] = 0.08
+            parameters['u'] = parameters['b'] / 2.0
+        else:
+            parameters['slope'] = 0.05
+            parameters['u'] = (parameters['b'] - parameters['tw']) / 2.0
+        return _create_taper_U_section(
+            **{
+                key: parameters[key]
+                for key in parameters
+                if key in ['h', 'b', 'tw', 'tf', 'r1', 'r2', 'slope', 'u']
+            }
+        )
+
+    @classmethod
+    def profiles(cls) -> list:
+        """Returns a list containing all available profiles."""
+        return list(cls.parameters.keys())
+
+    def __init__(self, name: str) -> None:
+        """Creates a new UPN object."""
+        if isinstance(name, (float, int)):
+            name = f'UPN{int(name):0d}'
+        parameters = self.parameters.get(name)
+        if parameters is None:
+            raise ValueError(
+                f"Profile '{name}' not found in UPN sections. "
+                "Select a valid profile (available ones: "
+                f"{self.profiles()})"
+            )
+        super().__init__()
+        self._h = parameters.get('h')
+        self._b = parameters.get('b')
+        self._tw = parameters.get('tw')
+        self._tf = parameters.get('tf')
+        self._r1 = parameters.get('r1')
+        self._r2 = parameters.get('r2')
+        if self._h <= 300:
+            self._flange_slope = 0.08
+            self._u = self._b / 2.0
+        else:
+            self._flange_slope = 0.05
+            self._u = (self._b - self._tw) / 2.0
+        self._polygon = _create_taper_U_section(
+            h=self._h,
+            b=self._b,
+            tw=self._tw,
+            tf=self._tf,
+            r1=self._r1,
+            r2=self._r2,
+            slope=self._flange_slope,
+            u=self._u,
+        )
+
+    @property
+    def polygon(self) -> Polygon:
+        """Returns shapely Polygon of section.
+
+        Returns:
+            Polygon: The represention of the UPN section.
+        """
+        return self._polygon
+
+    @property
+    def h(self) -> float:
+        """Returns height of UPN section.
+
+        Returns:
+            float: Height h of UPN section.
+        """
+        return self._h
+
+    @property
+    def b(self) -> float:
+        """Returns width of UPN section.
+
+        Returns:
+            float: Width b of UPN section.
+        """
+        return self._b
+
+    @property
+    def tw(self) -> float:
+        """Returns thickness of web of UPN section.
+
+        Returns:
+            float: Web thickness tw of UPN section.
+        """
+        return self._tw
+
+    @property
+    def tf(self) -> float:
+        """Returns thickness of flange of UPN section.
+
+        Returns:
+            float: Flange thickness tw of UPN section.
+        """
+        return self._tf
+
+    @property
+    def r1(self) -> float:
+        """Returns fillet radius of UPN section.
+
+        Returns:
+            float: Fillet radius r1 of UPN section.
+        """
+        return self._r1
+
+    @property
+    def r2(self) -> float:
+        """Returns fillet radius of UPN section.
+
+        Returns:
+            float: Fillet radius r2 of UPN section.
+        """
+        return self._r2

structuralcodes/materials/basic/_elastic.py

@@ -15,6 +15,9 @@ class ElasticMaterial(Material):
         self,
         E: float,
         density: float,
+        initial_strain: t.Optional[float] = None,
+        initial_stress: t.Optional[float] = None,
+        strain_compatibility: t.Optional[float] = None,
         name: t.Optional[str] = None,
     ):
         """Initialize a material with an elastic plastic constitutive law.
@@ -22,11 +25,25 @@ class ElasticMaterial(Material):
         Arguments:
             E (float): The Young's modulus.
             density (float): The density.
+            initial_strain (Optional[float]): Initial strain of the material.
+            initial_stress (Optional[float]): Initial stress of the material.
+            strain_compatibility (Optional[bool]): Only relevant if
+                initial_strain or initial_stress are different from zero. If
+                True, the material deforms with the geometry. If False, the
+                stress in the material upon loading is kept constant
+                corresponding to the initial strain.
             name (str, optional): The name of the material, default value None.
         """
-        super().__init__(density=density, name=name)
+        super().__init__(
+            density=density,
+            initial_strain=initial_strain,
+            initial_stress=initial_stress,
+            strain_compatibility=strain_compatibility,
+            name=name if name else 'ElasticMaterial',
+        )
         self._E = E
         self._constitutive_law = create_constitutive_law('elastic', self)
+        self._apply_initial_strain()
 
     @property
     def E(self) -> float:

structuralcodes/materials/basic/_elasticplastic.py

@@ -21,6 +21,9 @@ class ElasticPlasticMaterial(Material):
         density: float,
         Eh: float = 0,
         eps_su: t.Optional[float] = None,
+        initial_strain: t.Optional[float] = None,
+        initial_stress: t.Optional[float] = None,
+        strain_compatibility: t.Optional[float] = None,
         name: t.Optional[str] = None,
     ):
         """Initialize a material with an elastic plastic constitutive law.
@@ -31,9 +34,22 @@ class ElasticPlasticMaterial(Material):
             density (float): The density.
             Eh (float, optional): The hardening modulus, default value 0.
             eps_su (float, optional): The ultimate strain, default value None.
+            initial_strain (Optional[float]): Initial strain of the material.
+            initial_stress (Optional[float]): Initial stress of the material.
+            strain_compatibility (Optional[bool]): Only relevant if
+                initial_strain or initial_stress are different from zero. If
+                True, the material deforms with the geometry. If False, the
+                stress in the material upon loading is kept constant
+                corresponding to the initial strain.
             name (str, optional): The name of the material, default value None.
         """
-        super().__init__(density=density, name=name)
+        super().__init__(
+            density=density,
+            initial_strain=initial_strain,
+            initial_stress=initial_stress,
+            strain_compatibility=strain_compatibility,
+            name=name if name else 'ElasticPlasticMaterial',
+        )
         self._E = E
         self._fy = fy
         self._Eh = Eh
@@ -42,6 +58,7 @@ class ElasticPlasticMaterial(Material):
         self._constitutive_law = create_constitutive_law(
             'elasticplastic', self
         )
+        self._apply_initial_strain()
 
     @property
     def E(self) -> float:

structuralcodes/materials/basic/_generic.py

@@ -12,6 +12,9 @@ class GenericMaterial(Material):
         self,
         density: float,
         constitutive_law: ConstitutiveLaw,
+        initial_strain: t.Optional[float] = None,
+        initial_stress: t.Optional[float] = None,
+        strain_compatibility: t.Optional[bool] = None,
         name: t.Optional[str] = None,
     ):
         """Initialize a material with a constitutive law.
@@ -20,7 +23,21 @@ class GenericMaterial(Material):
             density (float): The density.
             constitutive_law (ConstitutiveLaw): The constitutive law of the
                 material.
+            initial_strain (Optional[float]): Initial strain of the material.
+            initial_stress (Optional[float]): Initial stress of the material.
+            strain_compatibility (Optional[bool]): Only relevant if
+                initial_strain or initial_stress are different from zero. If
+                True, the material deforms with the geometry. If False, the
+                stress in the material upon loading is kept constant
+                corresponding to the initial strain.
             name (str, optional): The name of the material, default value None.
         """
-        super().__init__(density=density, name=name)
+        super().__init__(
+            density=density,
+            initial_strain=initial_strain,
+            initial_stress=initial_stress,
+            strain_compatibility=strain_compatibility,
+            name=name if name else 'GenericMaterial',
+        )
         self._constitutive_law = constitutive_law
+        self._apply_initial_strain()

structuralcodes/materials/concrete/__init__.py

@@ -49,6 +49,9 @@ def create_concrete(
             desired standard. If None (default) the globally used design
             standard will be adopted. Otherwise the design standard specified
             will be used for the instance of the material.
+        **kwargs: Other valid keyword arguments that are collected and passed
+            to the specific concrete material. Please inspect the documentation
+            of the other concrete materials to see valid arguments.
 
     Raises:
         ValueError: if the design code is not valid or does not cover concrete

structuralcodes/materials/concrete/_concrete.py

@@ -21,10 +21,19 @@ class Concrete(Material):
         density: float = 2400,
         gamma_c: t.Optional[float] = None,
         existing: t.Optional[bool] = False,
+        initial_strain: t.Optional[float] = None,
+        initial_stress: t.Optional[float] = None,
+        strain_compatibility: t.Optional[bool] = None,
     ) -> None:
         """Initializes an abstract concrete material."""
         name = name if name is not None else 'Concrete'
-        super().__init__(density=density, name=name)
+        super().__init__(
+            density=density,
+            initial_strain=initial_strain,
+            initial_stress=initial_stress,
+            strain_compatibility=strain_compatibility,
+            name=name,
+        )
 
         self._fck = abs(fck)
         if existing:

structuralcodes/materials/concrete/_concreteEC2_2004.py

@@ -46,6 +46,9 @@ class ConcreteEC2_2004(Concrete):  # noqa: N801
                 ConstitutiveLaw,
             ]
         ] = 'parabolarectangle',
+        initial_strain: t.Optional[float] = None,
+        initial_stress: t.Optional[float] = None,
+        strain_compatibility: t.Optional[bool] = None,
         fcm: t.Optional[float] = None,
         fctm: t.Optional[float] = None,
         fctk_5: t.Optional[float] = None,
@@ -80,6 +83,13 @@ class ConcreteEC2_2004(Concrete):  # noqa: N801
                 law type for concrete. (valid options for string: 'elastic',
                 'parabolarectangle', 'bilinearcompression', 'sargin',
                 'popovics').
+            initial_strain (Optional[float]): Initial strain of the material.
+            initial_stress (Optional[float]): Initial stress of the material.
+            strain_compatibility (Optional[bool]): Only relevant if
+                initial_strain or initial_stress are different from zero. If
+                True, the material deforms with the geometry. If False, the
+                stress in the material upon loading is kept constant
+                corresponding to the initial strain.
             fcm (float, optional): The mean compressive strength.
             fctm (float, optional): The mean tensile strength.
             fctk_5 (float, optional): The 5% fractile for the tensile strength.
@@ -130,6 +140,9 @@ class ConcreteEC2_2004(Concrete):  # noqa: N801
             density=density,
             existing=False,
             gamma_c=gamma_c,
+            initial_strain=initial_strain,
+            initial_stress=initial_stress,
+            strain_compatibility=strain_compatibility,
         )
         self._alpha_cc = alpha_cc
         self._fcm = abs(fcm) if fcm is not None else None
@@ -165,6 +178,7 @@ class ConcreteEC2_2004(Concrete):  # noqa: N801
             raise ValueError(
                 'The provided constitutive law is not valid for concrete.'
             )
+        self._apply_initial_strain()
 
     def __post_init__(self):
         """Validator for the attributes that are set in the constructor."""

structuralcodes/materials/concrete/_concreteEC2_2023.py

@@ -48,6 +48,9 @@ class ConcreteEC2_2023(Concrete):  # noqa: N801
                 ConstitutiveLaw,
             ]
         ] = 'parabolarectangle',
+        initial_strain: t.Optional[float] = None,
+        initial_stress: t.Optional[float] = None,
+        strain_compatibility: t.Optional[bool] = None,
         fcm: t.Optional[float] = None,
         fctm: t.Optional[float] = None,
         fctk_5: t.Optional[float] = None,
@@ -80,6 +83,13 @@ class ConcreteEC2_2023(Concrete):  # noqa: N801
                 law type for concrete. (valid options for string: 'elastic',
                 'parabolarectangle', 'bilinearcompression', 'sargin',
                 'popovics').
+            initial_strain (Optional[float]): Initial strain of the material.
+            initial_stress (Optional[float]): Initial stress of the material.
+            strain_compatibility (Optional[bool]): Only relevant if
+                initial_strain or initial_stress are different from zero. If
+                True, the material deforms with the geometry. If False, the
+                stress in the material upon loading is kept constant
+                corresponding to the initial strain.
             fcm (float, optional): The mean compressive strength.
             fctm (float, optional): The mean tensile strength.
             fctk_5 (float, optional): The 5% fractile for the tensile strength.
@@ -124,6 +134,9 @@ class ConcreteEC2_2023(Concrete):  # noqa: N801
             density=density,
             existing=False,
             gamma_c=gamma_c,
+            initial_strain=initial_strain,
+            initial_stress=initial_stress,
+            strain_compatibility=strain_compatibility,
         )
         self._kE = kE
         self._strength_dev_class = strength_dev_class.strip().lower()
@@ -158,6 +171,7 @@ class ConcreteEC2_2023(Concrete):  # noqa: N801
             raise ValueError(
                 'The provided constitutive law is not valid for concrete.'
             )
+        self._apply_initial_strain()
 
     def __post_init__(self):
         """Validator for the attributes that are set in the constructor."""

structuralcodes/materials/concrete/_concreteMC2010.py

@@ -48,6 +48,9 @@ class ConcreteMC2010(Concrete):
                 ConstitutiveLaw,
             ]
         ] = 'parabolarectangle',
+        initial_strain: t.Optional[float] = None,
+        initial_stress: t.Optional[float] = None,
+        strain_compatibility: t.Optional[bool] = None,
         fcm: t.Optional[float] = None,
         fctm: t.Optional[float] = None,
         fctkmin: t.Optional[float] = None,
@@ -77,6 +80,18 @@ class ConcreteMC2010(Concrete):
             alpha_cc (float, optional): A factor for considering long-term
                 effects on the strength, and effects that arise from the way
                 the load is applied.
+            consitutive_law (ConstitutiveLaw | str): A valid ConstitutiveLaw
+                object for concrete or a string defining a valid constitutive
+                law type for concrete. (valid options for string: 'elastic',
+                'parabolarectangle', 'bilinearcompression', 'sargin',
+                'popovics').
+            initial_strain (Optional[float]): Initial strain of the material.
+            initial_stress (Optional[float]): Initial stress of the material.
+            strain_compatibility (Optional[bool]): Only relevant if
+                initial_strain or initial_stress are different from zero. If
+                True, the material deforms with the geometry. If False, the
+                stress in the material upon loading is kept constant
+                corresponding to the initial strain.
             fcm (float, optional): The mean compressive strength.
             fctm (float, optional): The mean tensile strength.
             fctkmin (float, optional): The minimum tensile strength.
@@ -126,6 +141,9 @@ class ConcreteMC2010(Concrete):
             name=name,
             density=density,
             gamma_c=gamma_c,
+            initial_strain=initial_strain,
+            initial_stress=initial_stress,
+            strain_compatibility=strain_compatibility,
         )
         self._alpha_cc = alpha_cc
         self._fcm = abs(fcm) if fcm is not None else None
@@ -162,6 +180,7 @@ class ConcreteMC2010(Concrete):
             raise ValueError(
                 'The provided constitutive law is not valid for concrete.'
             )
+        self._apply_initial_strain()
 
     def __post_init__(self):
         """Validator for the attributes that are set in the constructor."""

structuralcodes/materials/constitutive_laws/__init__.py

@@ -6,6 +6,7 @@ from ...core.base import ConstitutiveLaw, Material
 from ._bilinearcompression import BilinearCompression
 from ._elastic import Elastic
 from ._elasticplastic import ElasticPlastic
+from ._initial_strain import InitialStrain
 from ._parabolarectangle import ParabolaRectangle
 from ._popovics import Popovics
 from ._sargin import Sargin
@@ -19,6 +20,7 @@ __all__ = [
     'Popovics',
     'Sargin',
     'UserDefined',
+    'InitialStrain',
     'get_constitutive_laws_list',
     'create_constitutive_law',
 ]
@@ -31,6 +33,7 @@ CONSTITUTIVE_LAWS: t.Dict[str, ConstitutiveLaw] = {
     'parabolarectangle': ParabolaRectangle,
     'popovics': Popovics,
     'sargin': Sargin,
+    'initialstrain': InitialStrain,
 }
 
 

structuralcodes/materials/constitutive_laws/_elasticplastic.py

@@ -34,8 +34,8 @@ class ElasticPlastic(ConstitutiveLaw):
 
         Keyword Arguments:
             Eh (float): The hardening modulus.
-            eps_su (float): The ultimate strain.
-            name (str): A descriptive name for the constitutive law.
+            eps_su (float, optional): The ultimate strain.
+            name (str, optional): A descriptive name for the constitutive law.
         """
         name = name if name is not None else 'ElasticPlasticLaw'
         super().__init__(name=name)