sgio 0.2.13__py3-none-any.whl → 0.2.15__py3-none-any.whl

sgio/model/beam.py

@@ -1,111 +1,118 @@
 from __future__ import annotations
 
 import math
+from typing import Optional, List
+from pydantic import BaseModel, Field, field_validator, computed_field
 
 # from dataclasses import dataclass
 
 # @dataclass
-class EulerBernoulliBeamModel:
+class EulerBernoulliBeamModel(BaseModel):
     """Euler-Bernoulli Beam Model
     """
 
-    dim = 1
-    label = 'bm1'
-    model_name = 'Euler-Bernoulli beam model'
-
-    def __init__(self):
-        self.name = ''
-        self.id = None
-
-        # #: float: Geometric center location in x2 direction
-        # self.xg2 = 0.
-        # #: float: Geometric center location in x3 direction
-        # self.xg3 = 0.
-
-        # #: float: Area of the cross-section
-        # self.area = 0.
-
-
-        # Inertial properties
-        # -------------------
-
-        #: list of list of floats:
-        #: The 6x6 mass matrix
-        self.mass = None
-
-        #: list of lists of floats:
-        #: The 6x6 mass matrix at the mass center
-        self.mass_mc = None
-
-        #: float: Mass center location in x2 direction
-        self.xm2 = None
-        #: float: Mass center location in x3 direction
-        self.xm3 = None
-
-        #: float: Mass per unit span
-        self.mu = None
-        #: float: Mass moments of inertia i11
-        self.i11 = None
-        #: float: Principal mass moments of inertia i22
-        self.i22 = None
-        #: float: Principle mass moments of inertia i33
-        self.i33 = None
-        #: float: Principal inertial axes rotation angle in degree
-        self.phi_pia = 0
-        #: float: mass-weighted radius of gyration
-        self.rg = None
-
-
-        # Structural properties
-        # ---------------------
-
-        #: list of lists of floats:
-        #: Classical stiffness matrix (1-extension; 2-twist; 3,4-bending)
-        self.stff = None
-        #: list of lists of floats:
-        #: Classical compliance matrix (1-extension; 2-twist; 3,4-bending)
-        self.cmpl = None
-
-        #: float: Tension center location in x2 direction
-        self.xt2 = None
-        #: float: Tension center location in x3 direction
-        self.xt3 = None
-
-        #: float: Extension stiffness EA
-        self.ea = None
-        #: float: Torsional stiffness GJ
-        self.gj = None
-        #: float: Principal bending stiffness EI22
-        self.ei22 = None
-        #: float: Principal bending stiffness EI33
-        self.ei33 = None
-        #: float: Principle bending axes rotation angle in degree
-        self.phi_pba = 0
-
-        # #: list of lists of floats:
-        # #: Timoshenko stiffness matrix (1-extension; 2,3-shear, 4-twist; 5,6-bending)
-        # self.stff_t = []
-        # #: list of lists of floats:
-        # #: Timoshenko compliance matrix (1-extension; 2,3-shear, 4-twist; 5,6-bending)
-        # self.cmpl_t = []
-
-        # #: float: Generalized shear center location in x2 direction
-        # self.xs2 = None
-        # #: float: Generalized shear center location in x3 direction
-        # self.xs3 = None
-        # #: float: Principal shear stiffness GA22
-        # self.ga22 = None
-        # #: float: Principal shear stiffness GA33
-        # self.ga33 = None
-        # #: float: Principal shear axes rotation angle in degree
-        # self.phi_psa = None
-
+    # Class attributes (not Pydantic fields)
+    dim: int = 1
+    label: str = 'bm1'
+    model_name: str = 'Euler-Bernoulli beam model'
+
+    # Basic properties
+    name: str = Field(default='', description="Beam name")
+    id: Optional[int] = Field(default=None, description="Beam ID")
+
+    # Inertial properties
+    mass: Optional[List[List[float]]] = Field(
+        default=None,
+        description="The 6x6 mass matrix"
+    )
+    mass_mc: Optional[List[List[float]]] = Field(
+        default=None,
+        description="The 6x6 mass matrix at the mass center"
+    )
+
+    # Mass center locations
+    xm2: Optional[float] = Field(default=None, description="Mass center location in x2 direction")
+    xm3: Optional[float] = Field(default=None, description="Mass center location in x3 direction")
+
+    # Mass properties
+    mu: Optional[float] = Field(default=None, ge=0, description="Mass per unit span")
+    i11: Optional[float] = Field(default=None, ge=0, description="Mass moments of inertia i11")
+    i22: Optional[float] = Field(default=None, ge=0, description="Principal mass moments of inertia i22")
+    i33: Optional[float] = Field(default=None, ge=0, description="Principal mass moments of inertia i33")
+    phi_pia: float = Field(default=0, description="Principal inertial axes rotation angle in degree")
+    rg: Optional[float] = Field(default=None, ge=0, description="Mass-weighted radius of gyration")
+
+    # Structural properties
+    stff: Optional[List[List[float]]] = Field(
+        default=None,
+        description="Classical stiffness matrix (1-extension; 2-twist; 3,4-bending)"
+    )
+    cmpl: Optional[List[List[float]]] = Field(
+        default=None,
+        description="Classical compliance matrix (1-extension; 2-twist; 3,4-bending)"
+    )
+
+    # Tension center locations
+    xt2: Optional[float] = Field(default=None, description="Tension center location in x2 direction")
+    xt3: Optional[float] = Field(default=None, description="Tension center location in x3 direction")
+
+    # Stiffness properties
+    ea: Optional[float] = Field(default=None, ge=0, description="Extension stiffness EA")
+    gj: Optional[float] = Field(default=None, ge=0, description="Torsional stiffness GJ")
+    ei22: Optional[float] = Field(default=None, ge=0, description="Principal bending stiffness EI22")
+    ei33: Optional[float] = Field(default=None, ge=0, description="Principal bending stiffness EI33")
+    phi_pba: float = Field(default=0, description="Principal bending axes rotation angle in degree")
+
+    # Pydantic configuration
+    model_config = {"arbitrary_types_allowed": True}
+
+    # Field validators
+    @field_validator('mass', 'mass_mc')
+    @classmethod
+    def validate_6x6_matrix(cls, v):
+        """Validate that mass matrices are 6x6"""
+        if v is not None:
+            if not isinstance(v, list) or len(v) != 6:
+                raise ValueError('Matrix must be 6x6 (6 rows)')
+            for i, row in enumerate(v):
+                if not isinstance(row, list) or len(row) != 6:
+                    raise ValueError(f'Row {i} must have 6 columns')
+        return v
+
+    @field_validator('stff', 'cmpl')
+    @classmethod
+    def validate_4x4_matrix(cls, v):
+        """Validate that stiffness/compliance matrices are 4x4"""
+        if v is not None:
+            if not isinstance(v, list) or len(v) != 4:
+                raise ValueError('Matrix must be 4x4 (4 rows)')
+            for i, row in enumerate(v):
+                if not isinstance(row, list) or len(row) != 4:
+                    raise ValueError(f'Row {i} must have 4 columns')
+        return v
+
+    # Computed properties
+    @computed_field
     @property
-    def gyr1(self): return self.rg
+    def gyr1(self) -> Optional[float]:
+        """Mass-weighted radius of gyration (same as rg)"""
+        return self.rg
+
+    @computed_field
     @property
-    def gyr2(self): return math.sqrt(self.i22/self.mu)
+    def gyr2(self) -> Optional[float]:
+        """Radius of gyration about x2 axis"""
+        if self.i22 is not None and self.mu is not None and self.mu > 0:
+            return math.sqrt(self.i22 / self.mu)
+        return None
+
+    @computed_field
     @property
-    def gyr3(self): return math.sqrt(self.i33/self.mu)
+    def gyr3(self) -> Optional[float]:
+        """Radius of gyration about x3 axis"""
+        if self.i33 is not None and self.mu is not None and self.mu > 0:
+            return math.sqrt(self.i33 / self.mu)
+        return None
 
 
     def __repr__(self):
@@ -282,7 +289,9 @@ class EulerBernoulliBeamModel:
 
             # Mass
             if name.startswith('ms'):
-                return self.mass[int(name[2])-1][int(name[3])-1]
+                if self.mass is not None:
+                    return self.mass[int(name[2])-1][int(name[3])-1]
+                return None
             if name == 'mu':
                 return self.mu
             if name == 'mmoi1':
@@ -300,11 +309,15 @@ class EulerBernoulliBeamModel:
 
             # Stiffness
             if name.startswith('stf'):
-                return self.stff[int(name[3])-1][int(name[4])-1]
+                if self.stff is not None:
+                    return self.stff[int(name[3])-1][int(name[4])-1]
+                return None
 
             # Compliance
             if name.startswith('cmp'):
-                return self.cmpl[int(name[3])-1][int(name[4])-1]
+                if self.cmpl is not None:
+                    return self.cmpl[int(name[3])-1][int(name[4])-1]
+                return None
 
             if name == 'ea':
                 return self.ea
@@ -331,12 +344,17 @@ class EulerBernoulliBeamModel:
             if name == 'phi_pba':
                 return self.phi_pba
 
+            # Return None for unrecognized properties (e.g., ga22, ga33 which are Timoshenko-only)
+            return None
+
         elif isinstance(name, list) or isinstance(name, tuple):
             props = []
             for n in name:
                 props.append(self.get(n))
             return props
 
+        return None
+
 
     def getAll(self):
         """Get all beam properties.
@@ -815,884 +833,3 @@ class TimoshenkoBeamModel:
 
         return dict_prop
 
-
-
-
-
-
-
-
-
-# # Legacy
-
-# class BeamModel():
-#     """A beam property class (smdim = 1)
-
-#     """
-
-#     def __init__(self):
-#         # MaterialSection.__init__(self, smdim=1)
-
-#         #: float: Geometric center location in x2 direction
-#         self.xg2 = 0.
-#         #: float: Geometric center location in x3 direction
-#         self.xg3 = 0.
-
-#         #: float: Area of the cross-section
-#         self.area = 0.
-
-
-#         # Inertial properties
-#         # -------------------
-
-#         #: list of list of floats:
-#         #: The 6x6 mass matrix
-#         self.mass = []
-
-#         #: list of lists of floats:
-#         #: The 6x6 mass matrix at the mass center
-#         self.mass_cs = []
-
-#         #: float: Mass center location in x2 direction
-#         self.xm2 = 0.
-#         #: float: Mass center location in x3 direction
-#         self.xm3 = 0.
-
-#         #: float: Mass per unit span
-#         self.mu = 0.
-#         #: float: Mass moments of inertia i11
-#         self.i11 = 0.
-#         #: float: Principal mass moments of inertia i22
-#         self.i22 = 0.
-#         #: float: Principle mass moments of inertia i33
-#         self.i33 = 0.
-#         #: float: Principal inertial axes rotation angle in degree
-#         self.phi_pia = 0.
-#         #: float: mass-weighted radius of gyration
-#         self.rg = 0.
-
-
-#         # Structural properties
-#         # ---------------------
-
-#         #: list of lists of floats:
-#         #: Classical stiffness matrix (1-extension; 2-twist; 3,4-bending)
-#         self.stff = []
-#         #: list of lists of floats:
-#         #: Classical compliance matrix (1-extension; 2-twist; 3,4-bending)
-#         self.cmpl = []
-
-#         #: float: Tension center location in x2 direction
-#         self.xt2 = 0.
-#         #: float: Tension center location in x3 direction
-#         self.xt3 = 0.
-
-#         #: float: Extension stiffness EA
-#         self.ea = 0.
-#         #: float: Torsional stiffness GJ
-#         self.gj = 0.
-#         #: float: Principal bending stiffness EI22
-#         self.ei22 = 0.
-#         #: float: Principal bending stiffness EI33
-#         self.ei33 = 0.
-#         #: float: Principle bending axes rotation angle in degree
-#         self.phi_pba = 0.
-
-#         #: list of lists of floats:
-#         #: Timoshenko stiffness matrix (1-extension; 2,3-shear, 4-twist; 5,6-bending)
-#         self.stff_t = []
-#         #: list of lists of floats:
-#         #: Timoshenko compliance matrix (1-extension; 2,3-shear, 4-twist; 5,6-bending)
-#         self.cmpl_t = []
-
-#         #: float: Generalized shear center location in x2 direction
-#         self.xs2 = 0.
-#         #: float: Generalized shear center location in x3 direction
-#         self.xs3 = 0.
-#         #: float: Principal shear stiffness GA22
-#         self.ga22 = 0.
-#         #: float: Principal shear stiffness GA33
-#         self.ga33 = 0.
-#         #: float: Principal shear axes rotation angle in degree
-#         self.phi_psa = 0.
-
-
-#     def printData(self):
-#         fmt = '20.10E'
-#         print('Summary')
-#         print()
-#         print('The 6x6 mass matrix')
-#         print(utio.matrixToString(self.mass, fmt=fmt))
-#         print()
-#         print('The Mass Center of the Cross-Section')
-#         # print(f'({self.xm2:{fmt}}, {self.xm3:{fmt}})')
-#         print('({:{fmt}}, {:{fmt}})'.format(self.xm2, self.xm3, fmt=fmt))
-#         print()
-#         print('The Mass Properties with respect to Principal Inertial Axes')
-#         # print(f'Mass per unit span = {self.mu:{fmt}}')
-#         print('Mass per unit span = {:{fmt}}'.format(self.mu, fmt=fmt))
-#         # print(f'Mass moments of innertia i11 = {self.i11:{fmt}}')
-#         print('Mass moments of innertia i11 = {:{fmt}}'.format(self.i11, fmt=fmt))
-#         # print(f'Principle mass moments of innertia i22 = {self.i22:{fmt}}')
-#         print('Principle mass moments of innertia i22 = {:{fmt}}'.format(self.i22, fmt=fmt))
-#         # print(f'Principle mass moments of innertia i33 = {self.i33:{fmt}}')
-#         print('Principle mass moments of innertia i33 = {:{fmt}}'.format(self.i33, fmt=fmt))
-#         # print(f'The principal inertial axes rotation = {self.phi_pia:{fmt}}')
-#         print('The principal inertial axes rotation = {:{fmt}}'.format(self.phi_pia, fmt=fmt))
-#         # print(f'The mass-weighted radius of gyration = {self.rg:{fmt}}')
-#         print('The mass-weighted radius of gyration = {:{fmt}}'.format(self.rg, fmt=fmt))
-#         print()
-#         print('Classical Stiffness Matrix (1-extension; 2-twist; 3,4-bending)')
-#         print(utio.matrixToString(self.stff, fmt=fmt))
-#         print('Classical Flexibility Matrix (1-extension; 2-twist; 3,4-bending)')
-#         print(utio.matrixToString(self.cmpl, fmt=fmt))
-#         print()
-#         print('The Tension Center of the Cross-Section')
-#         # print(f'({self.xt2:{fmt}}, {self.xt3:{fmt}})')
-#         print('({:{fmt}}, {:{fmt}})'.format(self.xt2, self.xt3, fmt=fmt))
-#         print()
-#         # print(f'The extension stiffness EA = {self.ea:{fmt}}')
-#         print('The extension stiffness EA = {:{fmt}}'.format(self.ea, fmt=fmt))
-#         # print(f'The torsional stiffness GJ = {self.gj:{fmt}}')
-#         print('The torsional stiffness GJ = {:{fmt}}'.format(self.gj, fmt=fmt))
-#         # print(f'Principal bending stiffness EI22 = {self.ei22:{fmt}}')
-#         print('Principal bending stiffness EI22 = {:{fmt}}'.format(self.ei22, fmt=fmt))
-#         # print(f'Principal bending stiffness EI33 = {self.ei33:{fmt}}')
-#         print('Principal bending stiffness EI33 = {:{fmt}}'.format(self.ei33, fmt=fmt))
-#         # print(f'The principal bending axes rotation = {self.phi_pba:{fmt}}')
-#         print('The principal bending axes rotation = {:{fmt}}'.format(self.phi_pba, fmt=fmt))
-#         print()
-#         print('Timoshenko Stiffness Matrix (1-extension; 2,3-shear, 4-twist; 5,6-bending)')
-#         print(utio.matrixToString(self.stff_t, fmt=fmt))
-#         print('Timoshenko Flexibility Matrix (1-extension; 2,3-shear, 4-twist; 5,6-bending)')
-#         print(utio.matrixToString(self.cmpl_t, fmt=fmt))
-#         print()
-#         print('The Generalized Shear Center of the Cross-Section')
-#         # print(f'({self.xs2:{fmt}}, {self.xs3:{fmt}})')
-#         print('({:{fmt}}, {:{fmt}})'.format(self.xs2, self.xs3, fmt=fmt))
-#         print()
-#         # print(f'Principal shear stiffness GA22 = {self.ga22:{fmt}}')
-#         print('Principal shear stiffness GA22 = {:{fmt}}'.format(self.ga22, fmt=fmt))
-#         # print(f'Principal shear stiffness GA33 = {self.ga33:{fmt}}')
-#         print('Principal shear stiffness GA33 = {:{fmt}}'.format(self.ga33, fmt=fmt))
-#         # print(f'The principal shear axes rotation = {self.phi_psa:{fmt}}')
-#         print('The principal shear axes rotation = {:{fmt}}'.format(self.phi_psa, fmt=fmt))
-
-
-#     def get(self, name):
-#         """Get beam properties using specific names.
-
-#         Parameters
-#         ----------
-#         name : str
-#             Name of the property that will be returned.
-
-#         Returns
-#         -------
-#         float:
-#             Value of the specified beam property.
-
-#         Notes
-#         -----
-
-#         ..  list-table:: Inertial properties
-#             :header-rows: 1
-
-#             * - Name
-#               - Description
-#             * - ``msijo`` (``i``, ``j`` are numbers 1 to 6)
-#               - Entry (i, j) of the 6x6 mass matrix at the origin
-#             * - ``msijc`` (``i``, ``j`` are numbers 1 to 6)
-#               - Entry (i, j) of the 6x6 mass matrix at the mass center
-#             * - ``mu``
-#               - Mass per unit length
-#             * - ``mmoi1`` | ``mmoi2`` | ``mmoi3``
-#               - Mass moment of inertia about x1/x2/x3 axis
-
-#         ..  list-table:: Stiffness properties
-#             :header-rows: 1
-
-#             * - Name
-#               - Description
-#             * - ``stfijc`` (``i``, ``j`` are numbers 1 to 6)
-#               - Entry (i, j) of the 4x4 classical stiffness matrix
-#             * - ``stfijr`` (``i``, ``j`` are numbers 1 to 6)
-#               - Entry (i, j) of the 6x6 refined stiffness matrix
-#             * - ``eac`` | ``ear``
-#               - Axial stiffness of the classical/refined model
-#             * - ``gjc`` | ``gjr``
-#               - Torsional stiffness of the classical/refined model
-#             * - ``ei2c`` | ``eifc`` | ``ei2r`` | ``eifr``
-#               - Bending stiffness around x2 (flapwise) of the classical/refined model
-#             * - ``ei3c`` | ``eicc`` | ``ei3r`` | ``eicr``
-#               - Bending stiffness around x3 (chordwise or lead-lag) of the classical/refined model
-#             * - ``cmpijc`` (``i``, ``j`` are numbers 1 to 6)
-#               - Entry (i, j) of the 4x4 classical compliance matrix
-#             * - ``cmpijr`` (``i``, ``j`` are numbers 1 to 6)
-#               - Entry (i, j) of the 6x6 refined compliance matrix
-
-#         ..  list-table:: Center offsets
-#             :header-rows: 1
-
-#             * - Name
-#               - Description
-#             * - ``mcy`` | ``mc2``
-#               - y (or x2) component of the mass center
-#             * - ``mcz`` | ``mc3``
-#               - z (or x3) component of the mass center
-#             * - ``tcy`` | ``tc2``
-#               - y (or x2) component of the tension center
-#             * - ``tcz`` | ``tc3``
-#               - z (or x3) component of the tension center
-#             * - ``scy`` | ``sc2``
-#               - y (or x2) component of the shear center
-#             * - ``scz`` | ``sc3``
-#               - z (or x3) component of the shear center
-
-#         .
-
-#         """
-
-#         if isinstance(name, str):
-#             name = name.lower()
-
-#             # Mass
-#             if name.startswith('ms'):
-#                 return self.mass[int(name[2])-1][int(name[3])-1]
-#             if name == 'mu':
-#                 return self.mu
-#             if name == 'mmoi1':
-#                 return self.i11
-#             if name == 'mmoi2':
-#                 return self.i22
-#             if name == 'mmoi3':
-#                 return self.i33
-
-#             # Stiffness
-#             if name.startswith('stf'):
-#                 if name[-1] == 'c':
-#                     return self.stff[int(name[3])-1][int(name[4])-1]
-#                 elif name[-1] == 'r':
-#                     return self.stff_t[int(name[3])-1][int(name[4])-1]
-
-#             # Compliance
-#             if name.startswith('cmp'):
-#                 if name[-1] == 'c':
-#                     return self.cmpl[int(name[3])-1][int(name[4])-1]
-#                 elif name[-1] == 'r':
-#                     return self.cmpl_t[int(name[3])-1][int(name[4])-1]
-
-#             if name == 'ea':
-#                 return self.ea
-#             if name in ['ga22', 'gayy', 'ga2', 'gay']:
-#                 return self.ga22
-#             if name in ['ga33', 'gazz', 'ga3', 'gaz']:
-#                 return self.ga33
-#             if name == 'gj':
-#                 return self.gj
-#             if name in ['ei22', 'eiyy', 'ei2', 'eiy']:
-#                 return self.ei22
-#             if name in ['ei33', 'eizz', 'ei3', 'eiz']:
-#                 return self.ei33
-
-#             # Various centers
-#             if name == 'mcy' or name == 'mc2':
-#                 return self.xm2
-#             if name == 'mcz' or name == 'mc3':
-#                 return self.xm3
-#             if name == 'tcy' or name == 'tc2':
-#                 return self.xt2
-#             if name == 'tcz' or name == 'tc3':
-#                 return self.xt3
-#             if name == 'scy' or name == 'sc2':
-#                 return self.xs2
-#             if name == 'scz' or name == 'sc3':
-#                 return self.xs3
-
-#         elif isinstance(name, list) or isinstance(name, tuple):
-#             props = []
-#             for n in name:
-#                 props.append(self.get(n))
-#             return props
-
-
-#     def getAll(self):
-#         """Get all beam properties.
-
-#         Returns
-#         -------
-#         dict:
-#             A Dictionary of all beam properties.
-
-#         Notes
-#         -----
-
-#         Names are
-
-#         - mu, mmoi1, mmoi2, mmoi3
-#         - ea, ga22, ga33, gj, ei22, ei33
-#         - mc2, mc3, tc2, tc3, sc2, sc3
-#         - msij, stfijc, cmpijc, stfijr, cmpijr
-
-#         """
-#         names = [
-#             'mu', 'mmoi1', 'mmoi2', 'mmoi3',
-#             'ea', 'ga22', 'ga33', 'gj', 'ei22', 'ei33',
-#             'mc2', 'mc3', 'tc2', 'tc3', 'sc2', 'sc3'
-#         ]
-#         for i in range(4):
-#             for j in range(4):
-#                 names.append('stf{}{}c'.format(i+1, j+1))
-#                 names.append('cmp{}{}c'.format(i+1, j+1))
-#         for i in range(6):
-#             for j in range(6):
-#                 names.append('ms{}{}'.format(i+1, j+1))
-#                 names.append('stf{}{}r'.format(i+1, j+1))
-#                 names.append('cmp{}{}r'.format(i+1, j+1))
-
-#         dict_prop = {}
-#         for n in names:
-#             dict_prop[n] = self.get(n)
-
-#         return dict_prop
-
-
-
-
-
-
-
-
-
-# import copy
-# import numpy as np
-# import sgio.utils.io as utio
-
-# from .general import MaterialSection
-
-
-# class BeamProperty(MaterialSection):
-#     """A beam property class (smdim = 1)
-
-#     """
-
-#     def __init__(self):
-#         MaterialSection.__init__(self, smdim=1)
-
-#         #: float: Geometric center location in x2 direction
-#         self.xg2 = 0.
-#         #: float: Geometric center location in x3 direction
-#         self.xg3 = 0.
-
-#         #: float: Area of the cross-section
-#         self.area = 0.
-
-
-#         # Inertial properties
-#         # -------------------
-
-#         #: list of list of floats:
-#         #: The 6x6 mass matrix
-#         self.mass = []
-
-#         #: list of lists of floats:
-#         #: The 6x6 mass matrix at the mass center
-#         self.mass_cs = []
-
-#         #: float: Mass center location in x2 direction
-#         self.xm2 = 0.
-#         #: float: Mass center location in x3 direction
-#         self.xm3 = 0.
-
-#         #: float: Mass per unit span
-#         self.mu = 0.
-#         #: float: Mass moments of inertia i11
-#         self.i11 = 0.
-#         #: float: Principal mass moments of inertia i22
-#         self.i22 = 0.
-#         #: float: Principle mass moments of inertia i33
-#         self.i33 = 0.
-#         #: float: Principal inertial axes rotation angle in degree
-#         self.phi_pia = 0.
-#         #: float: mass-weighted radius of gyration
-#         self.rg = 0.
-
-
-#         # Structural properties
-#         # ---------------------
-
-#         #: list of lists of floats:
-#         #: Classical stiffness matrix (1-extension; 2-twist; 3,4-bending)
-#         self.stff = []
-#         #: list of lists of floats:
-#         #: Classical compliance matrix (1-extension; 2-twist; 3,4-bending)
-#         self.cmpl = []
-
-#         #: float: Tension center location in x2 direction
-#         self.xt2 = 0.
-#         #: float: Tension center location in x3 direction
-#         self.xt3 = 0.
-
-#         #: float: Extension stiffness EA
-#         self.ea = 0.
-#         #: float: Torsional stiffness GJ
-#         self.gj = 0.
-#         #: float: Principal bending stiffness EI22
-#         self.ei22 = 0.
-#         #: float: Principal bending stiffness EI33
-#         self.ei33 = 0.
-#         #: float: Principle bending axes rotation angle in degree
-#         self.phi_pba = 0.
-
-#         #: list of lists of floats:
-#         #: Timoshenko stiffness matrix (1-extension; 2,3-shear, 4-twist; 5,6-bending)
-#         self.stff_t = []
-#         #: list of lists of floats:
-#         #: Timoshenko compliance matrix (1-extension; 2,3-shear, 4-twist; 5,6-bending)
-#         self.cmpl_t = []
-
-#         #: float: Generalized shear center location in x2 direction
-#         self.xs2 = 0.
-#         #: float: Generalized shear center location in x3 direction
-#         self.xs3 = 0.
-#         #: float: Principal shear stiffness GA22
-#         self.ga22 = 0.
-#         #: float: Principal shear stiffness GA33
-#         self.ga33 = 0.
-#         #: float: Principal shear axes rotation angle in degree
-#         self.phi_psa = 0.
-
-
-#     def printData(self):
-#         fmt = '20.10E'
-#         print('Summary')
-#         print()
-#         print('The 6x6 mass matrix')
-#         print(utio.matrixToString(self.mass, fmt=fmt))
-#         print()
-#         print('The Mass Center of the Cross-Section')
-#         # print(f'({self.xm2:{fmt}}, {self.xm3:{fmt}})')
-#         print('({:{fmt}}, {:{fmt}})'.format(self.xm2, self.xm3, fmt=fmt))
-#         print()
-#         print('The Mass Properties with respect to Principal Inertial Axes')
-#         # print(f'Mass per unit span = {self.mu:{fmt}}')
-#         print('Mass per unit span = {:{fmt}}'.format(self.mu, fmt=fmt))
-#         # print(f'Mass moments of innertia i11 = {self.i11:{fmt}}')
-#         print('Mass moments of innertia i11 = {:{fmt}}'.format(self.i11, fmt=fmt))
-#         # print(f'Principle mass moments of innertia i22 = {self.i22:{fmt}}')
-#         print('Principle mass moments of innertia i22 = {:{fmt}}'.format(self.i22, fmt=fmt))
-#         # print(f'Principle mass moments of innertia i33 = {self.i33:{fmt}}')
-#         print('Principle mass moments of innertia i33 = {:{fmt}}'.format(self.i33, fmt=fmt))
-#         # print(f'The principal inertial axes rotation = {self.phi_pia:{fmt}}')
-#         print('The principal inertial axes rotation = {:{fmt}}'.format(self.phi_pia, fmt=fmt))
-#         # print(f'The mass-weighted radius of gyration = {self.rg:{fmt}}')
-#         print('The mass-weighted radius of gyration = {:{fmt}}'.format(self.rg, fmt=fmt))
-#         print()
-#         print('Classical Stiffness Matrix (1-extension; 2-twist; 3,4-bending)')
-#         print(utio.matrixToString(self.stff, fmt=fmt))
-#         print('Classical Flexibility Matrix (1-extension; 2-twist; 3,4-bending)')
-#         print(utio.matrixToString(self.cmpl, fmt=fmt))
-#         print()
-#         print('The Tension Center of the Cross-Section')
-#         # print(f'({self.xt2:{fmt}}, {self.xt3:{fmt}})')
-#         print('({:{fmt}}, {:{fmt}})'.format(self.xt2, self.xt3, fmt=fmt))
-#         print()
-#         # print(f'The extension stiffness EA = {self.ea:{fmt}}')
-#         print('The extension stiffness EA = {:{fmt}}'.format(self.ea, fmt=fmt))
-#         # print(f'The torsional stiffness GJ = {self.gj:{fmt}}')
-#         print('The torsional stiffness GJ = {:{fmt}}'.format(self.gj, fmt=fmt))
-#         # print(f'Principal bending stiffness EI22 = {self.ei22:{fmt}}')
-#         print('Principal bending stiffness EI22 = {:{fmt}}'.format(self.ei22, fmt=fmt))
-#         # print(f'Principal bending stiffness EI33 = {self.ei33:{fmt}}')
-#         print('Principal bending stiffness EI33 = {:{fmt}}'.format(self.ei33, fmt=fmt))
-#         # print(f'The principal bending axes rotation = {self.phi_pba:{fmt}}')
-#         print('The principal bending axes rotation = {:{fmt}}'.format(self.phi_pba, fmt=fmt))
-#         print()
-#         print('Timoshenko Stiffness Matrix (1-extension; 2,3-shear, 4-twist; 5,6-bending)')
-#         print(utio.matrixToString(self.stff_t, fmt=fmt))
-#         print('Timoshenko Flexibility Matrix (1-extension; 2,3-shear, 4-twist; 5,6-bending)')
-#         print(utio.matrixToString(self.cmpl_t, fmt=fmt))
-#         print()
-#         print('The Generalized Shear Center of the Cross-Section')
-#         # print(f'({self.xs2:{fmt}}, {self.xs3:{fmt}})')
-#         print('({:{fmt}}, {:{fmt}})'.format(self.xs2, self.xs3, fmt=fmt))
-#         print()
-#         # print(f'Principal shear stiffness GA22 = {self.ga22:{fmt}}')
-#         print('Principal shear stiffness GA22 = {:{fmt}}'.format(self.ga22, fmt=fmt))
-#         # print(f'Principal shear stiffness GA33 = {self.ga33:{fmt}}')
-#         print('Principal shear stiffness GA33 = {:{fmt}}'.format(self.ga33, fmt=fmt))
-#         # print(f'The principal shear axes rotation = {self.phi_psa:{fmt}}')
-#         print('The principal shear axes rotation = {:{fmt}}'.format(self.phi_psa, fmt=fmt))
-
-
-#     def get(self, name):
-#         """Get beam properties using specific names.
-
-#         Parameters
-#         ----------
-#         name : str
-#             Name of the property that will be returned.
-
-#         Returns
-#         -------
-#         float:
-#             Value of the specified beam property.
-
-#         Notes
-#         -----
-
-#         ..  list-table:: Inertial properties
-#             :header-rows: 1
-
-#             * - Name
-#               - Description
-#             * - ``msijo`` (``i``, ``j`` are numbers 1 to 6)
-#               - Entry (i, j) of the 6x6 mass matrix at the origin
-#             * - ``msijc`` (``i``, ``j`` are numbers 1 to 6)
-#               - Entry (i, j) of the 6x6 mass matrix at the mass center
-#             * - ``mu``
-#               - Mass per unit length
-#             * - ``mmoi1`` | ``mmoi2`` | ``mmoi3``
-#               - Mass moment of inertia about x1/x2/x3 axis
-
-#         ..  list-table:: Stiffness properties
-#             :header-rows: 1
-
-#             * - Name
-#               - Description
-#             * - ``stfijc`` (``i``, ``j`` are numbers 1 to 6)
-#               - Entry (i, j) of the 4x4 classical stiffness matrix
-#             * - ``stfijr`` (``i``, ``j`` are numbers 1 to 6)
-#               - Entry (i, j) of the 6x6 refined stiffness matrix
-#             * - ``eac`` | ``ear``
-#               - Axial stiffness of the classical/refined model
-#             * - ``gjc`` | ``gjr``
-#               - Torsional stiffness of the classical/refined model
-#             * - ``ei2c`` | ``eifc`` | ``ei2r`` | ``eifr``
-#               - Bending stiffness around x2 (flapwise) of the classical/refined model
-#             * - ``ei3c`` | ``eicc`` | ``ei3r`` | ``eicr``
-#               - Bending stiffness around x3 (chordwise or lead-lag) of the classical/refined model
-#             * - ``cmpijc`` (``i``, ``j`` are numbers 1 to 6)
-#               - Entry (i, j) of the 4x4 classical compliance matrix
-#             * - ``cmpijr`` (``i``, ``j`` are numbers 1 to 6)
-#               - Entry (i, j) of the 6x6 refined compliance matrix
-
-#         ..  list-table:: Center offsets
-#             :header-rows: 1
-
-#             * - Name
-#               - Description
-#             * - ``mcy`` | ``mc2``
-#               - y (or x2) component of the mass center
-#             * - ``mcz`` | ``mc3``
-#               - z (or x3) component of the mass center
-#             * - ``tcy`` | ``tc2``
-#               - y (or x2) component of the tension center
-#             * - ``tcz`` | ``tc3``
-#               - z (or x3) component of the tension center
-#             * - ``scy`` | ``sc2``
-#               - y (or x2) component of the shear center
-#             * - ``scz`` | ``sc3``
-#               - z (or x3) component of the shear center
-
-#         .
-
-#         """
-
-#         # mm = self.mass_origin
-#         # stf_c = self.stiffness
-#         # cmp_c = self.compliance
-#         # stf_r = self.stiffness_refined
-#         # cmp_r = self.compliance_refined
-
-#         # if type(center).__name__ == 'str':
-#         #     if center == 'tc':
-#         #         mm, stf_c, cmp_c, stf_r, cmp_r = calcOffsetBeamProperty(self.tension_center[1], self.tension_center[2])
-#         #     elif center == 'sc':
-#         #         mm, stf_c, cmp_c, stf_r, cmp_r = calcOffsetBeamProperty(self.shear_center[1], self.shear_center[2])
-#         #     elif center == 'mc':
-#         #         mm, stf_c, cmp_c, stf_r, cmp_r = calcOffsetBeamProperty(self.mass_center[1], self.mass_center[2])
-
-#         if isinstance(name, str):
-#             name = name.lower()
-
-#             # Mass
-#             if name.startswith('ms'):
-#                 return self.mass[int(name[2])-1][int(name[3])-1]
-#             if name == 'mu':
-#                 return self.mu
-#             if name == 'mmoi1':
-#                 return self.i11
-#             if name == 'mmoi2':
-#                 return self.i22
-#             if name == 'mmoi3':
-#                 return self.i33
-
-#             # Stiffness
-#             if name.startswith('stf'):
-#                 if name[-1] == 'c':
-#                     return self.stff[int(name[3])-1][int(name[4])-1]
-#                 elif name[-1] == 'r':
-#                     return self.stff_t[int(name[3])-1][int(name[4])-1]
-
-#             # Compliance
-#             if name.startswith('cmp'):
-#                 if name[-1] == 'c':
-#                     return self.cmpl[int(name[3])-1][int(name[4])-1]
-#                 elif name[-1] == 'r':
-#                     return self.cmpl_t[int(name[3])-1][int(name[4])-1]
-
-#             if name == 'ea':
-#                 return self.ea
-#             if name in ['ga22', 'gayy', 'ga2', 'gay']:
-#                 return self.ga22
-#             if name in ['ga33', 'gazz', 'ga3', 'gaz']:
-#                 return self.ga33
-#             if name == 'gj':
-#                 return self.gj
-#             if name in ['ei22', 'eiyy', 'ei2', 'eiy']:
-#                 return self.ei22
-#             if name in ['ei33', 'eizz', 'ei3', 'eiz']:
-#                 return self.ei33
-
-#             # Various centers
-#             if name == 'mcy' or name == 'mc2':
-#                 return self.xm2
-#             if name == 'mcz' or name == 'mc3':
-#                 return self.xm3
-#             if name == 'tcy' or name == 'tc2':
-#                 return self.xt2
-#             if name == 'tcz' or name == 'tc3':
-#                 return self.xt3
-#             if name == 'scy' or name == 'sc2':
-#                 return self.xs2
-#             if name == 'scz' or name == 'sc3':
-#                 return self.xs3
-
-#         elif isinstance(name, list) or isinstance(name, tuple):
-#             props = []
-#             for n in name:
-#                 props.append(self.get(n))
-#             return props
-
-
-#     def getAll(self):
-#         """Get all beam properties.
-
-#         Returns
-#         -------
-#         dict:
-#             A Dictionary of all beam properties.
-
-#         Notes
-#         -----
-
-#         Names are
-
-#         - mu, mmoi1, mmoi2, mmoi3
-#         - ea, ga22, ga33, gj, ei22, ei33
-#         - mc2, mc3, tc2, tc3, sc2, sc3
-#         - msij, stfijc, cmpijc, stfijr, cmpijr
-
-#         """
-#         names = [
-#             'mu', 'mmoi1', 'mmoi2', 'mmoi3',
-#             'ea', 'ga22', 'ga33', 'gj', 'ei22', 'ei33',
-#             'mc2', 'mc3', 'tc2', 'tc3', 'sc2', 'sc3'
-#         ]
-#         for i in range(4):
-#             for j in range(4):
-#                 names.append('stf{}{}c'.format(i+1, j+1))
-#                 names.append('cmp{}{}c'.format(i+1, j+1))
-#         for i in range(6):
-#             for j in range(6):
-#                 names.append('ms{}{}'.format(i+1, j+1))
-#                 names.append('stf{}{}r'.format(i+1, j+1))
-#                 names.append('cmp{}{}r'.format(i+1, j+1))
-
-#         dict_prop = {}
-#         for n in names:
-#             dict_prop[n] = self.get(n)
-
-#         return dict_prop
-
-
-#     def calcPropertyAt(self, new_origin):
-#         """Offset the beam reference center and recalculate beam properties.
-
-#         Parameters
-#         ----------
-#         offset_x2 : float
-#             x2 of the offset of the new center with respect to the current one.
-#         offset_x3 : float
-#             x3 of the offset of the new center with respect to the current one.
-
-#         """
-
-#         bp_new = copy.copy(self)
-
-#         # offset = [0.0, offset[0], offset[1]]
-#         if isinstance(new_origin, str):
-#             if new_origin == 'gc':
-#                 offset_x2 = self.xg2
-#                 offset_x3 = self.xg3
-#             elif new_origin == 'mc':
-#                 offset_x2 = self.xm2
-#                 offset_x3 = self.xm3
-#             elif new_origin == 'tc':
-#                 offset_x2 = self.xt2
-#                 offset_x3 = self.xt3
-#             elif new_origin == 'sc':
-#                 offset_x2 = self.xs2
-#                 offset_x3 = self.xs3
-#         else:
-#             offset_x2 = new_origin[0]
-#             offset_x3 = new_origin[1]
-
-#         bp_new.xg2 -= offset_x2
-#         bp_new.xg3 -= offset_x3
-
-#         # Offset mass matrix
-#         mm_o = np.asarray(self.mass_cs)
-#         if (offset_x2 != self.xm2) or (offset_x3 != self.xm3):
-#             # mm_c = np.asarray(self.mass_mc)
-#             mu = mm_o[0, 0]
-#             mi_c = mm_o[3:, 3:]
-
-#             x2 = self.xm2 - offset_x2
-#             x3 = self.xm3 - offset_x3
-#             r_tilde = np.array([
-#                 [0, -x3, x2],
-#                 [x3, 0, 0],
-#                 [-x2, 0, 0]
-#             ])
-
-#             mm_o[:3, 3:] = mu * r_tilde.T
-#             mm_o[3:, :3] = mu * r_tilde
-
-#             # I_o = I_c + m * r_tilde.r_tilde^T
-#             mm_o[3:, 3:] = mm_o[3:, 3:] + mu * np.dot(r_tilde, r_tilde.T)
-#         bp_new.mass = mm_o
-#         bp_new.xm2 -= offset_x2
-#         bp_new.xm3 -= offset_x3
-
-
-#         # Offset stiffness and compliance
-#         trfm_4 = np.eye(4)
-#         trfm_6 = np.eye(6)
-
-#         trfm_4[2, 0] = offset_x3
-#         trfm_4[3, 0] = -offset_x2
-
-#         trfm_6[4, 0] = offset_x3
-#         trfm_6[5, 0] = -offset_x2
-#         trfm_6[3, 1] = -offset_x3
-#         trfm_6[3, 2] = offset_x2
-
-#         cmp_4 = np.asarray(self.cmpl)
-#         cmp_6 = np.asarray(self.cmpl_t)
-
-#         bp_new.cmpl = np.dot(trfm_4.T, np.dot(cmp_4, trfm_4))
-#         bp_new.cmpl_t = np.dot(trfm_6.T, np.dot(cmp_6, trfm_6))
-
-#         bp_new.stff = np.linalg.inv(bp_new.cmpl)
-#         bp_new.stff_t = np.linalg.inv(bp_new.cmpl_t)
-
-#         bp_new.xt2 -= offset_x2
-#         bp_new.xt3 -= offset_x3
-
-#         bp_new.xs2 -= offset_x2
-#         bp_new.xs3 -= offset_x3
-
-#         return bp_new
-
-
-#     def writeToFile(self, fn, fmt='vabs'):
-#         with open(fn, 'w') as fo:
-#             if fmt.startswith('v'):
-#                 self.writeToFileVABS(fo)
-
-#         return
-
-
-#     def writeToFileVABS(self, fo):
-#         fmt_float = '20.10E'
-
-#         fo.write('\n The 6X6 Mass Matrix\n')
-#         fo.write(' '+'='*56)
-#         fo.write('\n\n')
-#         utio.writeFormatFloatsMatrix(fo, self.mass, fmt=fmt_float, indent=1)
-
-#         fo.write('\n The Mass Center of the Cross Section\n')
-#         fo.write(' '+'='*56)
-#         fo.write('\n\n')
-#         utio.writeFormatFloats(fo, (self.xm2, self.xm3), fmt=fmt_float, indent=1)
-
-#         fo.write('\n The 6X6 Mass Matrix at the Mass Center\n')
-#         fo.write(' '+'='*56)
-#         fo.write('\n\n')
-#         utio.writeFormatFloatsMatrix(fo, self.mass_cs, fmt='20.12E', indent=1)
-
-#         fo.write('\n The Mass Properties with respect to Principal Inertial Axes\n')
-#         fo.write(' '+'='*56)
-#         fo.write('\n\n')
-#         fo.write(' {:39s}={:20.10E}\n'.format('Mass per unit span', self.mu))
-#         fo.write(' {:39s}={:20.10E}\n'.format('Mass moment of inertia i11', self.i11))
-#         fo.write(' {:39s}={:20.10E}\n'.format('Principal mass moments of inertia i22', self.i22))
-#         fo.write(' {:39s}={:20.10E}\n'.format('Principal mass moments of inertia i33', self.i33))
-#         fo.write(' The principal inertial axes rotated from user coordinate system by {} degrees about the positive direction of x1 axis.\n'.format(self.phi_pia))
-#         fo.write(' {:39s}={:20.10E}\n'.format('The mass-weighted radius of gyration', self.rg))
-
-#         fo.write('\n The Geometric Center of the Cross Section\n')
-#         fo.write(' '+'='*56)
-#         fo.write('\n\n')
-#         utio.writeFormatFloats(fo, (self.xg2, self.xg3), fmt=fmt_float, indent=1)
-
-#         fo.write('\n The Area of the Cross Section\n')
-#         fo.write(' '+'='*56)
-#         fo.write('\n\n')
-#         fo.write(' Area ={:20.10E}\n'.format(self.area))
-
-#         fo.write('\n Classical Stiffness Matrix (1-extension; 2-twist; 3,4-bending)\n')
-#         fo.write(' '+'='*56)
-#         fo.write('\n\n')
-#         utio.writeFormatFloatsMatrix(fo, self.stff, fmt='20.10E', indent=1)
-
-#         fo.write('\n Classical Compliance Matrix (1-extension; 2-twist; 3,4-bending)\n')
-#         fo.write(' '+'='*56)
-#         fo.write('\n\n')
-#         utio.writeFormatFloatsMatrix(fo, self.cmpl, fmt='20.10E', indent=1)
-
-#         fo.write('\n The Tension Center of the Cross Section\n')
-#         fo.write(' '+'='*56)
-#         fo.write('\n\n')
-#         utio.writeFormatFloats(fo, (self.xt2, self.xt3), fmt=fmt_float, indent=1)
-
-#         fo.write('\n')
-#         fo.write(' {:34s}={:20.10E}\n'.format('The extension stiffness EA', self.ea))
-#         fo.write(' {:34s}={:20.10E}\n'.format('The torsional stiffness GJ', self.gj))
-#         fo.write(' {:34s}={:20.10E}\n'.format('Principal bending stiffness EI22', self.ei22))
-#         fo.write(' {:34s}={:20.10E}\n'.format('Principal bending stiffness EI33', self.ei33))
-#         fo.write(' The principal bending axes rotated from the user coordinate system by {} degrees about the positive direction of x1 axis.\n'.format(self.phi_pba))
-
-#         fo.write('\n Timoshenko Stiffness Matrix (1-extension; 2,3-shear, 4-twist; 5,6-bending)\n')
-#         fo.write(' '+'='*56)
-#         fo.write('\n\n')
-#         utio.writeFormatFloatsMatrix(fo, self.stff_t, fmt=fmt_float, indent=1)
-
-#         fo.write('\n Timoshenko Compliance Matrix (1-extension; 2,3-shear, 4-twist; 5,6-bending)\n')
-#         fo.write(' '+'='*56)
-#         fo.write('\n\n')
-#         utio.writeFormatFloatsMatrix(fo, self.cmpl_t, fmt=fmt_float, indent=1)
-
-#         fo.write('\n The Shear Center of the Cross Section in the User Coordinate System\n')
-#         fo.write(' '+'='*56)
-#         fo.write('\n\n')
-#         utio.writeFormatFloats(fo, (self.xs2, self.xs3), fmt=fmt_float, indent=1)
-
-#         fo.write('\n')
-#         fo.write(' {:31s}={:20.10E}\n'.format('Principal shear stiffness GA22', self.ga22))
-#         fo.write(' {:31s}={:20.10E}\n'.format('Principal shear stiffness GA33', self.ga33))
-#         fo.write(' The principal shear axes rotated from user coordinate system by {} degrees about the positive direction of x1 axis.\n'.format(self.phi_psa))
-
-#         return