PyCBA 0.7.0__tar.gz → 0.8.0__tar.gz

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

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

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

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

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

@@ -13,12 +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_moving_load.py
+tests/test_nonlinear.py
+tests/test_nonprismatic.py
 tests/test_results.py

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

@@ -2,10 +2,11 @@
 PyCBA - Continuous Beam Analysis in Python
 """
 
-__version__ = "0.7.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.7.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:
@@ -100,6 +104,9 @@ class BeamAnalysis:
             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:
@@ -243,6 +250,29 @@ 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,