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

structuralcodes/geometry/profiles/_uc.py

@@ -0,0 +1,276 @@
+"""UC profiles."""
+
+from shapely import (
+    Polygon,
+)
+
+from ._base_profile import BaseProfile
+from ._common_functions import (
+    _create_I_section,
+)
+
+
+class UC(BaseProfile):
+    """Simple class for representing a UC profile.
+
+    Universal Columns.
+    """
+
+    parameters = {
+        'UC152x152x23': {
+            'h': 152.0,
+            'b': 152.0,
+            'tw': 5.8,
+            'tf': 6.8,
+            'r': 8.0,
+        },
+        'UC152x152x30': {
+            'h': 158.0,
+            'b': 153.0,
+            'tw': 6.5,
+            'tf': 9.4,
+            'r': 8.0,
+        },
+        'UC152x152x37': {
+            'h': 162.0,
+            'b': 154.0,
+            'tw': 8.0,
+            'tf': 11.5,
+            'r': 8.0,
+        },
+        'UC152x152x44': {
+            'h': 166.0,
+            'b': 156.0,
+            'tw': 9.5,
+            'tf': 13.6,
+            'r': 8.0,
+        },
+        'UC152x152x51': {
+            'h': 170.0,
+            'b': 157.0,
+            'tw': 11.0,
+            'tf': 15.7,
+            'r': 8.0,
+        },
+        'UC203x203x46': {
+            'h': 203.0,
+            'b': 204.0,
+            'tw': 7.2,
+            'tf': 11.0,
+            'r': 13.0,
+        },
+        'UC203x203x52': {
+            'h': 206.0,
+            'b': 204.0,
+            'tw': 7.9,
+            'tf': 12.5,
+            'r': 13.0,
+        },
+        'UC203x203x60': {
+            'h': 210.0,
+            'b': 206.0,
+            'tw': 9.4,
+            'tf': 14.2,
+            'r': 13.0,
+        },
+        'UC203x203x71': {
+            'h': 216.0,
+            'b': 206.0,
+            'tw': 10.0,
+            'tf': 17.3,
+            'r': 13.0,
+        },
+        'UC203x203x86': {
+            'h': 222.0,
+            'b': 209.0,
+            'tw': 12.7,
+            'tf': 20.5,
+            'r': 13.0,
+        },
+        'UC203x203x100': {
+            'h': 229.0,
+            'b': 210.0,
+            'tw': 14.5,
+            'tf': 23.7,
+            'r': 13.0,
+        },
+        'UC203x203x113': {
+            'h': 235.0,
+            'b': 212.0,
+            'tw': 16.3,
+            'tf': 26.9,
+            'r': 13.0,
+        },
+        'UC203x203x127': {
+            'h': 241.0,
+            'b': 214.0,
+            'tw': 18.1,
+            'tf': 30.1,
+            'r': 13.0,
+        },
+        'UC254x254x73': {
+            'h': 254.0,
+            'b': 255.0,
+            'tw': 8.6,
+            'tf': 14.2,
+            'r': 20.0,
+        },
+        'UC254x254x81': {
+            'h': 256.0,
+            'b': 255.0,
+            'tw': 9.4,
+            'tf': 15.6,
+            'r': 20.0,
+        },
+        'UC254x254x89': {
+            'h': 260.0,
+            'b': 256.0,
+            'tw': 10.3,
+            'tf': 17.3,
+            'r': 20.0,
+        },
+        'UC254x254x101': {
+            'h': 264.0,
+            'b': 257.0,
+            'tw': 11.9,
+            'tf': 19.6,
+            'r': 20.0,
+        },
+        'UC254x254x107': {
+            'h': 267.0,
+            'b': 259.0,
+            'tw': 12.8,
+            'tf': 20.5,
+            'r': 20.0,
+        },
+        'UC254x254x115': {
+            'h': 269.0,
+            'b': 259.0,
+            'tw': 13.5,
+            'tf': 22.1,
+            'r': 20.0,
+        },
+        'UC254x254x132': {
+            'h': 276.0,
+            'b': 261.0,
+            'tw': 15.3,
+            'tf': 25.3,
+            'r': 20.0,
+        },
+        'UC254x254x149': {
+            'h': 282.0,
+            'b': 263.0,
+            'tw': 17.3,
+            'tf': 28.4,
+            'r': 20.0,
+        },
+        'UC254x254x167': {
+            'h': 289.0,
+            'b': 265.0,
+            'tw': 19.2,
+            'tf': 31.7,
+            'r': 20.0,
+        },
+        'UC305x305x97': {
+            'h': 308.0,
+            'b': 305.0,
+            'tw': 9.9,
+            'tf': 15.4,
+            'r': 20.0,
+        },
+        'UC305x305x107': {
+            'h': 311.0,
+            'b': 306.0,
+            'tw': 10.9,
+            'tf': 17.0,
+            'r': 20.0,
+        },
+    }
+
+    @classmethod
+    def get_polygon(cls, name: str) -> Polygon:
+        """Returns a shapely polygon representing a UC section."""
+        parameters = cls.parameters.get(name)
+        if parameters is None:
+            raise ValueError(
+                f"Profile '{name}' not found in UC sections. "
+                "Select a valid profile (available ones: "
+                f"{cls.profiles()})"
+            )
+        return _create_I_section(**parameters)
+
+    @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 UC object."""
+        parameters = self.parameters.get(name)
+        if parameters is None:
+            raise ValueError(
+                f"Profile '{name}' not found in UC 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._r = parameters.get('r')
+        self._polygon = _create_I_section(**parameters)
+
+    @property
+    def polygon(self) -> Polygon:
+        """Returns shapely Polygon of section.
+
+        Returns:
+            Polygon: The represention of the UC section.
+        """
+        return self._polygon
+
+    @property
+    def h(self) -> float:
+        """Returns height of UC section.
+
+        Returns:
+            float: Height h of UC section.
+        """
+        return self._h
+
+    @property
+    def b(self) -> float:
+        """Returns width of UC section.
+
+        Returns:
+            float: Width b of UC section.
+        """
+        return self._b
+
+    @property
+    def tw(self) -> float:
+        """Returns thickness of web of UC section.
+
+        Returns:
+            float: Web thickness tw of UC section.
+        """
+        return self._tw
+
+    @property
+    def tf(self) -> float:
+        """Returns thickness of flange of UC section.
+
+        Returns:
+            float: Flange thickness tw of UC section.
+        """
+        return self._tf
+
+    @property
+    def r(self) -> float:
+        """Returns fillet radius of UC section.
+
+        Returns:
+            float: Fillet radius r of UC section.
+        """
+        return self._r

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/__init__.py

@@ -1,9 +1,10 @@
 """Main entry point for materials."""
 
-from . import concrete, constitutive_laws, reinforcement
+from . import basic, concrete, constitutive_laws, reinforcement
 
 __all__ = [
     'concrete',
     'constitutive_laws',
     'reinforcement',
+    'basic',
 ]

structuralcodes/materials/basic/__init__.py

@@ -0,0 +1,11 @@
+"""A collection of basic material classes."""
+
+from ._elastic import ElasticMaterial
+from ._elasticplastic import ElasticPlasticMaterial
+from ._generic import GenericMaterial
+
+__all__ = [
+    'ElasticMaterial',
+    'ElasticPlasticMaterial',
+    'GenericMaterial',
+]

structuralcodes/materials/basic/_elastic.py

@@ -0,0 +1,69 @@
+"""A material class with elastic properties."""
+
+import typing as t
+
+from ...core.base import Material
+from ..constitutive_laws import create_constitutive_law
+
+
+class ElasticMaterial(Material):
+    """A material class with elastic properties."""
+
+    _E: float
+
+    def __init__(
+        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.
+
+        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,
+            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:
+        """Returns the Young's modulus."""
+        return self._E
+
+    def __elastic__(self) -> dict:
+        """Returns kwargs for creating an elastic constitutive law."""
+        return {'E': self.E}
+
+    @classmethod
+    def from_material(cls, other_material: Material):
+        """Create an elastic material based on another material."""
+        # Create name of elastic material
+        name = other_material.name
+        if name is not None:
+            name += '_elastic'
+
+        return cls(
+            E=other_material.constitutive_law.get_tangent(eps=0),
+            density=other_material.density,
+            name=name,
+        )