PyCBA 0.6.0__tar.gz → 0.8.0__tar.gz

{pycba-0.6.0/src/PyCBA.egg-info → pycba-0.8.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: PyCBA
-Version: 0.6.0
+Version: 0.8.0
 Summary: Python Continuous Beam Analysis
 Author-email: Colin Caprani <colin.caprani@monash.edu>
 License: Apache 2.0

{pycba-0.6.0 → pycba-0.8.0/src/PyCBA.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: PyCBA
-Version: 0.6.0
+Version: 0.8.0
 Summary: Python Continuous Beam Analysis
 Author-email: Colin Caprani <colin.caprani@monash.edu>
 License: Apache 2.0

{pycba-0.6.0 → pycba-0.8.0}/src/PyCBA.egg-info/SOURCES.txt

@@ -13,11 +13,17 @@ src/pycba/beam.py
 src/pycba/bridge.py
 src/pycba/inf_lines.py
 src/pycba/load.py
+src/pycba/nonlinear.py
 src/pycba/pattern.py
 src/pycba/results.py
+src/pycba/section.py
 src/pycba/types.py
 src/pycba/utils.py
 src/pycba/vehicle.py
 tests/test_basic.py
 tests/test_bridge.py
-tests/test_inf_lines.py
+tests/test_inf_lines.py
+tests/test_moving_load.py
+tests/test_nonlinear.py
+tests/test_nonprismatic.py
+tests/test_results.py

{pycba-0.6.0 → pycba-0.8.0}/src/pycba/__init__.py

@@ -2,10 +2,11 @@
 PyCBA - Continuous Beam Analysis in Python
 """
 
-__version__ = "0.6.0"
+__version__ = "0.8.0"
 
 from .analysis import BeamAnalysis
 from .beam import Beam
+from .section import SectionEI
 from .load import (
     LoadCNL,
     MemberResults,
@@ -21,3 +22,4 @@ from .utils import parse_beam_string
 from .bridge import BridgeAnalysis
 from .vehicle import Vehicle, make_train, VehicleLibrary
 from .pattern import LoadPattern
+from .nonlinear import NonlinearBeamAnalysis, NonlinearResult, HingeEvent

{pycba-0.6.0 → pycba-0.8.0}/src/pycba/analysis.py

@@ -78,9 +78,13 @@ class BeamAnalysis:
         ----------
         L : array_like of float
             Span lengths.  Length ``N`` for an ``N``-span beam.
-        EI : float or array_like of float
-            Flexural rigidity of each span.  A scalar value is applied to all
-            spans; otherwise one value per span is required.
+        EI : float, pycba.section.SectionEI, or array_like
+            Flexural rigidity of each span.  A single scalar (or a single
+            :class:`~pycba.section.SectionEI`) is applied to all spans;
+            otherwise one entry per span is required.  A span whose rigidity is
+            given as a :class:`~pycba.section.SectionEI` is treated as
+            **non-prismatic** (variable ``EI``) and analysed by flexibility
+            integration; scalar entries use the closed-form prismatic element.
         R : array_like of int or float
             Nodal restraint vector, length ``2(N+1)``.  Two entries per node
             (vertical DOF then rotational DOF), ordered left to right:
@@ -91,13 +95,18 @@ class BeamAnalysis:
             * ``+k`` — elastic spring stiffness in consistent units.
 
         LM : list of list, optional
-            Load matrix: a list of load descriptors, each of the form
-            ``[span, load_type, value, a, c]``.  Load types:
-
-            1. UDL — ``value`` is load intensity; set ``a = c = 0``.
-            2. Point load — ``value`` at distance ``a`` from the span start.
-            3. Partial UDL — ``value`` intensity from ``a`` for length ``c``.
-            4. Moment load — ``value`` at distance ``a`` from the span start.
+            Load matrix: a list of load descriptors.  The number of columns
+            per entry depends on the load type:
+
+            1. UDL — ``[span, 1, w]``.
+            2. Point load — ``[span, 2, P, a]``.
+            3. Partial UDL — ``[span, 3, w, a, c]``.
+            4. Moment load — ``[span, 4, M, a]``.
+            5. Trapezoidal — ``[span, 5, w1, w2]`` (full span) or
+               ``[span, 5, w1, w2, a, c]`` (partial).
+            6. Imposed curvature — ``[span, 6, k0, k1, ...]`` where the free
+               curvature field is ``κ(x) = k0 + k1·x + …`` (e.g. creep,
+               shrinkage or thermal curvature).
 
         eletype : array_like of int, optional
             Element type for each span, controlling which end(s) carry moment:
@@ -241,6 +250,67 @@ class BeamAnalysis:
         load = [i_span, 4, m, a]
         self._beam.add_load(load)
 
+    def add_ic(self, i_span: int, kappa):
+        r"""
+        Append an imposed-curvature (initial-strain) member load.
+
+        The free curvature field ``κ(x) = k0 + k1·x + k2·x² + …`` is imposed
+        over the member.  On a simply-supported span it produces no internal
+        forces (only a free deflected shape); on a restrained or continuous
+        structure its restraint generates real moments and reactions.  This is
+        the mechanism for applying creep, shrinkage and thermal curvatures to a
+        continuous beam (see :class:`pycba.load.LoadIC`).
+
+        Parameters
+        ----------
+        i_span : int
+            1-based span index.
+        kappa : float or array_like of float
+            Imposed-curvature polynomial coefficients in increasing powers of
+            ``x``: ``[k0, k1, k2, ...]``.  A scalar is a uniform curvature.
+        """
+        coeffs = np.atleast_1d(np.asarray(kappa, dtype=float)).tolist()
+        load = [i_span, 6] + coeffs
+        self._beam.add_load(load)
+
+    def add_trap(
+        self,
+        i_span: int,
+        w1: float,
+        w2: float,
+        a: Optional[float] = None,
+        c: Optional[float] = None,
+    ):
+        """
+        Append a trapezoidal (linearly varying) distributed load.
+
+        When *a* and *c* are omitted the load covers the full span, varying
+        from *w1* at the left end to *w2* at the right end.  When *a* and *c*
+        are given the load covers the region from *a* to *a + c*, varying from
+        *w1* to *w2* over that length.
+
+        Parameters
+        ----------
+        i_span : int
+            1-based span index.
+        w1 : float
+            Load intensity at the start of the load.  Positive values act downward.
+        w2 : float
+            Load intensity at the end of the load.  Positive values act downward.
+        a : float, optional
+            Distance from the left end of the span to the start of the load.
+            If given, *c* must also be provided.
+        c : float, optional
+            Length (cover) of the load.  Required when *a* is provided.
+        """
+        if a is not None and c is None:
+            raise ValueError("If 'a' is specified, 'c' must also be provided")
+        if a is not None:
+            load = [i_span, 5, w1, w2, a, c]
+        else:
+            load = [i_span, 5, w1, w2]
+        self._beam.add_load(load)
+
     def analyze(self, npts: Optional[int] = None) -> int:
         """
         Execute the direct-stiffness analysis.