pystran 0.2.0__tar.gz → 0.2.2__tar.gz

{pystran-0.2.0 → pystran-0.2.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pystran
-Version: 0.2.0
+Version: 0.2.2
 Summary: Python Structural Analysis package
 Project-URL: Homepage, https://github.com/PetrKryslUCSD/pystran
 Project-URL: Issues, https://github.com/PetrKryslUCSD/pystran/issues

{pystran-0.2.0 → pystran-0.2.2}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "hatchling.build"
 
 [project]
 name = "pystran"
-version = "0.2.0"
+version = "0.2.2"
 authors = [
   { name="Petr Krysl", email="pkrysl@ucsd.edu" },
 ]

{pystran-0.2.0 → pystran-0.2.2}/pystran/plots.py

@@ -157,6 +157,12 @@ def setup(m, set_limits=False, fontsize=0):
     m['fig'] = fig # plotting objects saved
     if m["dim"] == 3:
         ax = fig.add_subplot(projection="3d")
+        if set_limits:
+            box = bounding_box(m)
+            cd = characteristic_dimension(m)
+            ax.set_xlim3d([box[0] - cd / 10, box[3] + cd / 10])
+            ax.set_ylim3d([box[1] - cd / 10, box[4] + cd / 10])
+            ax.set_zlim3d([box[2] - cd / 10, box[5] + cd / 10])
     else:
         ax = fig.add_subplot()
         if set_limits:
@@ -164,9 +170,25 @@ def setup(m, set_limits=False, fontsize=0):
             cd = characteristic_dimension(m)
             ax.set_xlim([box[0] - cd / 10, box[2] + cd / 10])
             ax.set_ylim([box[1] - cd / 10, box[3] + cd / 10])
+    ax.set_aspect("equal")
     m['ax'] = ax # plotting objects saved
     return ax
 
+def show(m):
+    """
+    Show the plot.
+
+    Parameters
+    ----------
+    m
+        Model dictionary.
+    """
+    if not ('fig' in m):
+        raise RuntimeError('Please first call plots.setup(m)')
+    fig = m['fig']
+    ax = m['ax'] 
+    ax.set_aspect("equal")
+    plt.show()
 
 def _area_extrema(all_members):
     min_area = numpy.inf
@@ -410,6 +432,7 @@ def _plot_member_ids_2d(m):
             i, j = m["joints"][connectivity[0]], m["joints"][connectivity[1]]
             ci, cj = i["coordinates"], j["coordinates"]
             xm = (ci + cj) / 2.0
+            ax.scatter(xm[0], xm[1], alpha=0)
             ax.text(xm[0], xm[1], str(jid),
                     bbox=dict(facecolor='white', boxstyle='square'))
     if "beam_members" in m:
@@ -419,6 +442,7 @@ def _plot_member_ids_2d(m):
             i, j = m["joints"][connectivity[0]], m["joints"][connectivity[1]]
             ci, cj = i["coordinates"], j["coordinates"]
             xm = (ci + cj) / 2.0
+            ax.scatter(xm[0], xm[1], alpha=0)
             ax.text(xm[0], xm[1], str(jid),
                     bbox=dict(facecolor='white', boxstyle='square'))
     return ax
@@ -433,6 +457,7 @@ def _plot_member_ids_3d(m):
             i, j = m["joints"][connectivity[0]], m["joints"][connectivity[1]]
             ci, cj = i["coordinates"], j["coordinates"]
             xm = (ci + cj) / 2.0
+            ax.scatter(xm[0], xm[1], xm[2], alpha=0)
             ax.text(xm[0], xm[1], xm[2], str(jid), 
                     bbox=dict(facecolor='white', boxstyle='square'))
     if "beam_members" in m:
@@ -442,6 +467,7 @@ def _plot_member_ids_3d(m):
             i, j = m["joints"][connectivity[0]], m["joints"][connectivity[1]]
             ci, cj = i["coordinates"], j["coordinates"]
             xm = (ci + cj) / 2.0
+            ax.scatter(xm[0], xm[1], xm[2], alpha=0)
             ax.text(xm[0], xm[1], xm[2], str(jid), 
                     bbox=dict(facecolor='white', boxstyle='square'))
     return ax
@@ -506,6 +532,10 @@ def plot_joint_ids(m, offsets = None):
     ax = m['ax']
     for j in m["joints"].values():
         if m["dim"] == 3:
+            ax.scatter(
+                j["coordinates"][0] + offsets[0],
+                j["coordinates"][1] + offsets[1],
+                j["coordinates"][2] + offsets[2], alpha=0)
             ax.text(
                 j["coordinates"][0] + offsets[0],
                 j["coordinates"][1] + offsets[1],
@@ -514,6 +544,9 @@ def plot_joint_ids(m, offsets = None):
                 bbox=dict(facecolor='white', boxstyle='circle')
             )
         else:
+            ax.scatter(
+                j["coordinates"][0] + offsets[0],
+                j["coordinates"][1] + offsets[1], alpha=0)
             ax.text(j["coordinates"][0] + offsets[0], 
                     j["coordinates"][1] + offsets[1], 
                     str(j["jid"]),
@@ -1338,58 +1371,6 @@ def plot_rotation_supports(m, scale=0.0, radius=0.0, shortest_arrow=1.0e-6):
                             )
     return ax
 
-
-def show(m):
-    """
-    Show the plot.
-
-    Parameters
-    ----------
-    m
-        Model dictionary.
-
-    """
-    if not ('fig' in m):
-        raise RuntimeError('Please first call plots.setup(m)')
-    fig = m['fig']
-    ax = m['ax']
-    ax.set_aspect("equal")
-    box = bounding_box(m)
-    cd = characteristic_dimension(m)
-    if m["dim"] == 3:
-        ax.set_xlim([box[0] - cd / 10, box[3] + cd / 10])
-        ax.set_ylim([box[1] - cd / 10, box[4] + cd / 10])
-        ax.set_zlim([box[2] - cd / 10, box[5] + cd / 10])
-        ax.set_xlabel("X")
-        ax.set_ylabel("Y")
-        ax.set_zlabel("Z")
-    else:
-        ax.set_xlim([box[0] - cd / 10, box[2] + cd / 10])
-        ax.set_ylim([box[1] - cd / 10, box[3] + cd / 10])
-        ax.set_xlabel("X")
-        ax.set_ylabel("Z")
-    # zoom in 3d will change ticks
-    def update_ticks(ax, n_max=6):
-        dim = m['dim']
-        # set an appropriate number of major ticks for each axis based on span
-        if dim == 2:
-            for axis in (ax.xaxis, ax.yaxis,):
-                axis.set_major_locator(MaxNLocator(n_max))
-        else:
-            for axis in (ax.xaxis, ax.yaxis, ax.zaxis):
-                axis.set_major_locator(MaxNLocator(n_max))
-        fig.canvas.draw_idle()
-    
-    def on_limits_change(event_ax):
-        # event_ax is the axis object passed by the callback
-        update_ticks(ax, n_max=6)
-    # Connect per-axis limit-change callbacks
-    ax.callbacks.connect('xlim_changed', on_limits_change)
-    ax.callbacks.connect('ylim_changed', on_limits_change)
-    ax.callbacks.connect('zlim_changed', on_limits_change)
-    plt.show()
-
-
 def plot_reaction_forces(m, scale=0.0):
     """
     Plot the reaction forces at the joints.

{pystran-0.2.0 → pystran-0.2.2}/tests/unittests_planar_frame.py

@@ -293,7 +293,7 @@ class UnitTestsPlanarFrames(unittest.TestCase):
         """
         pystran - Python package for structural analysis with trusses and beams
 
-        (C) 2025, Petr Krysl, pkrysl@ucsd.edu
+        (C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
         # Two story planar frame vibration
 
@@ -382,7 +382,7 @@ class UnitTestsPlanarFrames(unittest.TestCase):
         """
         pystran - Python package for structural analysis with trusses and beams
 
-        (C) 2025, Petr Krysl, pkrysl@ucsd.edu
+        (C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
         # Natural Frequency of Mass supported by a Beam on Springs
 

{pystran-0.2.0 → pystran-0.2.2}/tests/unittests_planar_truss.py

@@ -213,7 +213,7 @@ class UnitTestsPlanarTrusses(unittest.TestCase):
         """
         pystran - Python package for structural analysis with trusses and beams
 
-        (C) 2025, Petr Krysl, pkrysl@ucsd.edu
+        (C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
         # Truss with a combination of loads
 

{pystran-0.2.0 → pystran-0.2.2}/tests/unittests_space_frame.py

@@ -808,7 +808,7 @@ class UnitTestsSpaceFrames(unittest.TestCase):
         """
         pystran - Python package for structural analysis with trusses and beams
 
-        (C) 2025, Petr Krysl, pkrysl@ucsd.edu
+        (C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
         # Two story 3D Frame vibration
 
@@ -949,7 +949,7 @@ class UnitTestsSpaceFrames(unittest.TestCase):
         """
         pystran - Python package for structural analysis with trusses and beams
 
-        (C) 2025, Petr Krysl, pkrysl@ucsd.edu
+        (C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
         # Spatial frame with a hinge and spring supports
 

{pystran-0.2.0 → pystran-0.2.2}/tutorials/00_matrices_tut.py

@@ -1,7 +1,7 @@
 """
 pystran - Python package for structural analysis with trusses and beams 
 
-(C) 2025, Petr Krysl, pkrysl@ucsd.edu
+(C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
 # A gentle introduction to matrices in Python
 

{pystran-0.2.0 → pystran-0.2.2}/tutorials/01_alt_three_bars_tut.py

@@ -1,15 +1,15 @@
 """
 pystran - Python package for structural analysis with trusses and beams.
 
-(C) 2025, Petr Krysl, pkrysl@ucsd.edu
+(C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
 # Three-bar truss example
 
 ## Problem description:
 
 Structure consisting of three truss members. This tutorial mirrors the tutorial
-01_three_bars_tut.py, but uses the alternative (mixed) identifiers for
-creating the model.
+01_three_bars_tut.py, but uses the alternative (mixed, integers and 
+strings) identifiers for creating the model.
 
 Displacements and internal forces are provided in the verification manual.
 
@@ -31,7 +31,7 @@ from pystran import plots
 
 # Now we create the model `m`, which is a dictionary that contains all the
 # constituent parts of the system. The model is created for a planar structure
-# (2). The deformation is assumed to occur in the plane `x-z`.
+# (the 2). The deformation is assumed to occur in the plane `x-z`.
 m = model.create(2)
 freedoms = m["freedoms"]
 
@@ -46,8 +46,7 @@ model.add_joint(m, "B", [0.0, 20.0])
 model.add_joint(m, "C", [0.0, 10.0])
 model.add_joint(m, "D", (+10.0, 0.0))
 
-# For instance, the material parameters and the cross sectional area of the
-# bars.
+# The material parameters and the cross sectional area of the bars.
 E = 30000000.0
 A = 0.65700000
 
@@ -58,6 +57,8 @@ s1 = section.truss_section("steel", E, A)
 # lists of [start, end] joints. The section is also provided. Note that one
 # member is called "m1" and the other two are numbered. Also, the connectivity
 # is supplied in one case as tuple, in the other two cases as lists.
+# This is purely to show off the flexibility, there are no actual such requirements.
+
 model.add_truss_member(m, "m1", ["A", "B"], s1)
 model.add_truss_member(m, 2, ("A", "C"), s1)
 model.add_truss_member(m, 3, ["A", "D"], s1)
@@ -74,12 +75,12 @@ ax = plots.setup(m)
 plots.plot_members(m)
 plots.plot_member_orientation(m)
 plots.plot_joint_ids(m)
-plots.plot_translation_supports(m)
 ax.set_title("Structure")
 plots.show(m)
 
 # At this point we can visualize the supports. The translation supports are
-# shown with arrow heads.
+# shown with arrow heads. Only when the displacements prescribed are non zero,
+# the arrows will have visible shafts.
 ax = plots.setup(m)
 plots.plot_joint_ids(m)
 plots.plot_translation_supports(m)
@@ -90,8 +91,8 @@ plots.show(m)
 model.add_load(jA, 0, -10000.0)
 model.add_load(jA, 1, -10000.0 / 2.0)
 
-# The model is shown graphically, members are displayed with the member numbers attached.
-# The applied forces are also shown.
+# The model is shown graphically, members are displayed with the member 
+# numbers attached. The applied forces are also shown.
 plots.setup(m)
 plots.plot_members(m)
 plots.plot_applied_forces(m, 1 / 3000.0)
@@ -119,7 +120,9 @@ if norm(jA["displacements"] - [-0.0033325938, -0.001591621]) > 1.0e-3:
 print("Displacement calculation OK")
 
 # The forces in the bars can be calculated using the strain-displacement matrix
-# to compute the strain, from which we can evaluate the force.
+# to compute the strain, from which we can evaluate the force. (Note there 
+# is a function to do this, here we show how it could be done from 
+# individual components.)
 for b in m["truss_members"].values():
     connectivity = b["connectivity"]
     i, j = m["joints"][connectivity[0]], m["joints"][connectivity[1]]
@@ -129,7 +132,7 @@ for b in m["truss_members"].values():
     eps = dot(B, u)
     print("Bar " + str(connectivity) + " force = ", E * A * eps[0])
 
-# The forces in the bars can be simply calculated using the `truss_axial_force`
+# The same calculation, made easier with the `truss_axial_force`
 # function. That function does what was described in the loop above.
 for b in m["truss_members"].values():
     connectivity = b["connectivity"]
@@ -140,7 +143,6 @@ for b in m["truss_members"].values():
 # These are the reference values of the forces from the book.
 print("Reference forces: ", -0.656854250e4, -0.48528137e4, -0.15685425e4)
 
-
 # The solution is visualized with deformed shape.
 plots.setup(m)
 plots.plot_members(m)

{pystran-0.2.0 → pystran-0.2.2}/tutorials/01_three_bars_tut.py

@@ -1,7 +1,7 @@
 """
 pystran - Python package for structural analysis with trusses and beams
 
-(C) 2025, Petr Krysl, pkrysl@ucsd.edu
+(C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
 # Three-bar truss example
 

{pystran-0.2.0 → pystran-0.2.2}/tutorials/02_continuous_beam_2_spans_consist_tut.py

@@ -1,16 +1,16 @@
 """
 pystran - Python package for structural analysis with trusses and beams
 
-(C) 2025, Petr Krysl, pkrysl@ucsd.edu
+(C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
 # Continuous beam with two spans
 
 ## Problem description:
 
 Continuous simply supported beam with two spans. The beam is loaded with
-moments at the joint.
+moments at the joints.
 
-Displacements and internal forces are provided in the verification manual.
+Displacements and internal forces are provided in the literature.
 
 ## References
 

{pystran-0.2.0 → pystran-0.2.2}/tutorials/03_weaver_1_tut.py

@@ -1,7 +1,7 @@
 """
 pystran - Python package for structural analysis with trusses and beams
 
-(C) 2025, Petr Krysl, pkrysl@ucsd.edu
+(C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
 # Weaver 3D frame benchmark
 
@@ -122,7 +122,8 @@ model.add_beam_member(m, 3, [2, 4], sect_2)
 
 
 # Now we can plot the geometry of the structure. We show the members, the
-# member numbers, and the orientations of the local coordinate systems.
+# member numbers, and the orientations of the local coordinate systems. 
+# Also the joint numbers are displayed.
 
 ax = plots.setup(m)
 plots.plot_members(m)
@@ -163,7 +164,7 @@ plots.show(m)
 # Now we can solve the static equilibrium of the frame. First we number the
 # degrees of freedom, and then we call the solver that will construct the
 # stiffness matrix, the right-hand side vector of applied generalized forces
-# (forces plus moments), and solve a system of equations. The solution is then
+# (forces plus moments), and solve the system of equations. The solution is then
 # distributed to the joints.
 model.number_dofs(m)
 model.solve_statics(m)
@@ -177,12 +178,12 @@ for jid in [1, 2]:
 # The displacements of the joints can be compared to the reference values.
 # These are the displacements of joint 1:
 ref1 = [0.22267, 0.00016, -0.17182, -0.00255, 0.00217, -0.00213]
-if norm(m["joints"][1]["displacements"] - ref1) > 1.0e-1 * norm(ref1):
+if norm(m["joints"][1]["displacements"] - ref1) > 1.0e-2 * norm(ref1):
     raise ValueError("Displacement calculation error")
 print("Displacement calculation OK")
 # These are the displacements of joint 2:
 ref2 = [0.22202, -0.48119, -0.70161, -0.00802, 0.00101, -0.00435]
-if norm(m["joints"][2]["displacements"] - ref2) > 1.0e-1 * norm(ref2):
+if norm(m["joints"][2]["displacements"] - ref2) > 1.0e-2 * norm(ref2):
     raise ValueError("Displacement calculation error")
 print("Displacement calculation OK")
 

{pystran-0.2.0 → pystran-0.2.2}/tutorials/04_supp_settle_tut.py

@@ -1,7 +1,7 @@
 """
 pystran - Python package for structural analysis with trusses and beams
 
-(C) 2025, Petr Krysl, pkrysl@ucsd.edu
+(C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
 # Example of a support-settlement problem
 
@@ -33,7 +33,7 @@ from pystran import plots
 E = 2.9e6
 I = 1.0
 A = 1.0  # cross-sectional area does not influence the results
-L = 10 * 12  # span in inchesc
+L = 10 * 12  # span in inches
 
 # The model is created as two dimensional.
 m = model.create(2)
@@ -46,7 +46,7 @@ model.add_joint(m, 3, [2 * L, 0.0])
 # The left hand side is clamped (all degrees of freedom set to zero), the other
 # joints are simply supported.
 model.add_support(m["joints"][1], freedoms.ALL_DOFS)
-# The middle support moves down by 0.25 inches (notice the non zero value of the
+# The middle support moves down by 0.25 inches (notice the non-zero value of the
 # enforced displacement).
 model.add_support(m["joints"][2], freedoms.U2, -0.25)
 model.add_support(m["joints"][3], freedoms.U2)
@@ -61,7 +61,7 @@ model.add_beam_member(m, 2, [2, 3], s1)
 # scale doesn't work well in this case).
 ax = plots.setup(m, set_limits=True)
 plots.plot_members(m)
-plots.plot_translation_supports(m, 50)
+plots.plot_translation_supports(m, scale=50)
 ax.set_title("Supports (magnified 50x)")
 plots.show(m)
 

{pystran-0.2.0 → pystran-0.2.2}/tutorials/05_hinge_frame_tut.py

@@ -1,7 +1,7 @@
 """
 pystran - Python package for structural analysis with trusses and beams
 
-(C) 2025, Petr Krysl, pkrysl@ucsd.edu
+(C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
 # Example of a frame with a hinge
 
@@ -28,10 +28,10 @@ from pystran import plots
 # US customary units, inches, pounds, seconds are assumed.
 
 # The book gives the product of the modulus of elasticity and the moment of
-# inertia as 2.9e6.
+# inertia as 2.9e8.
 E = 29e6
 I = 100.0
-A = 10.0  # cross-sectional area does not influence the results
+A = 10.0  # cross-sectional area does not influence the results (much)
 L = 10 * 12  # span in inches
 
 m = model.create(2)
@@ -45,9 +45,11 @@ model.add_joint(m, 4, [L, 0.0])
 
 # There needs to be a hinge (not transferring bending moments, but ensuring
 # continuity of displacements) at the joint 2. We add another joint, 5, at the
-# same location, and we link the degrees of freedom that needs to be the same.
+# same location.
 model.add_joint(m, 5, [0, L])
 
+# Now we link the degrees of freedom that need to be the same --
+# the translations are the same, rotations can be different.
 # We supply the list of joints that need to be linked (2 and 5), and the
 # degrees of freedom that are to be the same.
 model.add_dof_links(m, [2, 5], freedoms.U1)

{pystran-0.2.0 → pystran-0.2.2}/tutorials/06_sennett_tut.py

@@ -1,7 +1,7 @@
 """
 pystran - Python package for structural analysis with trusses and beams
 
-(C) 2025, Petr Krysl, pkrysl@ucsd.edu
+(C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
 # Example of a three dimensional frame problem
 

{pystran-0.2.0 → pystran-0.2.2}/tutorials/07_frame_thermal_tut.py

@@ -1,7 +1,7 @@
 """
 pystran - Python package for structural analysis with trusses and beams
 
-(C) 2025, Petr Krysl, pkrysl@ucsd.edu
+(C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
 # Example of a two-dimensional frame problem with thermal gradient loading
 
@@ -103,7 +103,7 @@ M_T = CTE * (Ttop - Tbot) / depth * E * I
 model.add_load(i, freedoms.UR3, -M_T)
 model.add_load(j, freedoms.UR3, +M_T)
 
-# The nodal moments can be visualized with the following plot.
+# The applied nodal moments can be visualized with the following plot.
 ax = plots.setup(m)
 plots.plot_members(m)
 ax = plots.plot_applied_moments(m, 0.0, 50)
@@ -158,7 +158,7 @@ if abs((f["Qzj"] + 1710) / 1710) > 1e-2:
 if abs((f["Myj"] + M_T) / 40171 - 1) > 1e-2:
     raise ValueError("Member 2, joint i, bending moment error")
 
-# There are diagrams below are for the discrete model that is loaded with the
+# The diagrams below are for the discrete model that is loaded with the
 # negative reactions at the fixed joints. These resultants are not combined with
 # the thermal forces, they represent the response of the frame with the nodal loads.
 

{pystran-0.2.0 → pystran-0.2.2}/tutorials/08_alt_refine_member_tut.py

@@ -1,7 +1,7 @@
 """
 pystran - Python package for structural analysis with trusses and beams
 
-(C) 2025, Petr Krysl, pkrysl@ucsd.edu
+(C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
 # Continuous beam with two spans: refine a member into multiple elements
 

{pystran-0.2.0 → pystran-0.2.2}/tutorials/08_refine_member_tut.py

@@ -1,7 +1,7 @@
 """
 pystran - Python package for structural analysis with trusses and beams
 
-(C) 2025, Petr Krysl, pkrysl@ucsd.edu
+(C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
 # Continuous beam with two spans: refine a member into multiple elements
 

{pystran-0.2.0 → pystran-0.2.2}/tutorials/09_SDLX_01_89_vibration_tut.py

@@ -1,7 +1,7 @@
 """
 pystran - Python package for structural analysis with trusses and beams
 
-(C) 2025, Petr Krysl, pkrysl@ucsd.edu
+(C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
 # Two story planar frame vibration
 
@@ -32,6 +32,7 @@ B = 4.8e-3
 A = H * B
 I = H * B**3 / 12
 
+# Create the model in usual fashion.
 m = model.create(2)
 freedoms = m["freedoms"]
 
@@ -65,7 +66,8 @@ plots.show(m)
 # All members will now be refined into eight finite elements. Without the
 # refinement, the reference solutions cannot be reproduced: there simply
 # wouldn't be enough degrees of freedom. Unfortunately the reference publication
-# does not mention the numbers of finite elements used per member.
+# does not mention with how many finite elements used per member 
+# the reference solution was calculated.
 nref = 8
 for i in range(6):
     model.refine_member(m, i + 1, nref)

{pystran-0.2.0 → pystran-0.2.2}/tutorials/10_truss_5_vibration_tut.py

@@ -1,7 +1,7 @@
 """
 pystran - Python package for structural analysis with trusses and beams
 
-(C) 2025, Petr Krysl, pkrysl@ucsd.edu
+(C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
 # Statically-determinate truss vibration
 

{pystran-0.2.0 → pystran-0.2.2}/tutorials/11_samcef_frame_3d_vibration_tut.py

@@ -1,7 +1,7 @@
 """
 pystran - Python package for structural analysis with trusses and beams
 
-(C) 2025, Petr Krysl, pkrysl@ucsd.edu
+(C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
 # Two story 3D Frame vibration
 

{pystran-0.2.0 → pystran-0.2.2}/tutorials/12_timo_beam_on_springs_tut.py

@@ -1,7 +1,7 @@
 """
 pystran - Python package for structural analysis with trusses and beams
 
-(C) 2025, Petr Krysl, pkrysl@ucsd.edu
+(C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
 # Natural Frequency of Mass supported by a Beam on Springs
 

{pystran-0.2.0 → pystran-0.2.2}/tutorials/13_hinged_3d_frame_tut.py

@@ -1,7 +1,7 @@
 """
 pystran - Python package for structural analysis with trusses and beams
 
-(C) 2025, Petr Krysl, pkrysl@ucsd.edu
+(C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
 # Spatial frame with a hinge and spring supports
 

{pystran-0.2.0 → pystran-0.2.2}/tutorials/14_truss_skew_support_tut.py

@@ -1,7 +1,7 @@
 """
 pystran - Python package for structural analysis with trusses and beams
 
-(C) 2025, Petr Krysl, pkrysl@ucsd.edu
+(C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
 # Truss with a combination of loads
 

{pystran-0.2.0 → pystran-0.2.2}/tutorials/15_penalty_supports_tut.py

@@ -1,7 +1,7 @@
 """
 pystran - Python package for structural analysis with trusses and beams
 
-(C) 2025, Petr Krysl, pkrysl@ucsd.edu
+(C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
 # Penalty enforcement of supports
 

{pystran-0.2.0 → pystran-0.2.2}/tutorials/16_distributed_load_continuous_beam_tut.py

@@ -1,7 +1,7 @@
 """
 pystran - Python package for structural analysis with trusses and beams
 
-(C) 2025, Petr Krysl, pkrysl@ucsd.edu
+(C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
 # Continuous beam with distributed loading
 

{pystran-0.2.0 → pystran-0.2.2}/tutorials/17_arch_out_of_plane_tut.py

@@ -1,7 +1,7 @@
 """
 pystran - Python package for structural analysis with trusses and beams
 
-(C) 2025, Petr Krysl, pkrysl@ucsd.edu
+(C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
 # Bending of a quarter circle under out-of-plane force
 

{pystran-0.2.0 → pystran-0.2.2}/tutorials/18_SSLL05_89_tut.py

@@ -1,7 +1,7 @@
 """
 pystran - Python package for structural analysis with trusses and beams
 
-(C) 2025, Petr Krysl, pkrysl@ucsd.edu
+(C) 2025-2026, Petr Krysl, pkrysl@ucsd.edu
 
 # Bending of a two incompressible cantilevers rigidly linked at the free end
 

{pystran-0.2.0 → pystran-0.2.2}/tutorials/19_garteur_air_frame_vibration_ref_tut.py

@@ -26,8 +26,9 @@ The test-bed was designed and manufactured by ONERA, France.
 ![](IMAC97photo.png)
 
 This tutorial includes also the constraining layer of the wing, which 
-was not included in the "19_garteur_air_frame_statics_ref_tut.py" tutorial.
-The effect of the constraining layer may be overestimated when computing 
+was not included in the "19_garteur_air_frame_statics_tut.py" tutorial.
+The effect of the constraining layer may be overestimated 
+(the frequencies may end up a little bit too high) when computing 
 the natural frequencies: the flexibility of the viscoelastic layer is omitted.
 
 References

{pystran-0.2.0 → pystran-0.2.2}/tutorials/19_garteur_air_frame_vibration_simpl_tut.py

@@ -26,7 +26,7 @@ The test-bed was designed and manufactured by ONERA, France.
 ![](IMAC97photo.png)
 
 This tutorial is simplified: no damping, no constraining layer, no 
-added masses (not on the drums, no masses for connecting metal).
+added masses (not on the drums, no masses for connecting fixtures).
 
 References
 

pystran-0.2.0/Job-1.inp

@@ -1,111 +0,0 @@
-*Heading
-** Job name: Job-1 Model name: Model-1
-** Generated by: Abaqus/CAE Learning Edition 2025
-*Preprint, echo=NO, model=NO, history=NO, contact=NO
-**
-** PARTS
-**
-*Part, name=Part-frame
-*Node
-      1,           0.,           0.,         180.
-      2,         180.,           0.,         180.
-      3,         180.,         180.,         180.
-      4,         180.,           0.,           0.
-*Element, type=B33
-1, 1, 2
-2, 2, 3
-3, 2, 4
-*Nset, nset=Set-1
- 1, 2
-*Elset, elset=Set-1
- 1,
-*Nset, nset=Set-2
- 1, 2
-*Elset, elset=Set-2
- 1,
-*Nset, nset=Set-3
- 1, 2
-*Elset, elset=Set-3
- 1,
-*Nset, nset=Set-4
- 2, 3
-*Elset, elset=Set-4
- 2,
-*Nset, nset=Set-5
- 2, 3
-*Elset, elset=Set-5
- 2,
-*Nset, nset=Set-6
- 2, 4
-*Elset, elset=Set-6
- 3,
-*Nset, nset=Set-7
- 2, 4
-*Elset, elset=Set-7
- 3,
-** Section: Section-member  Profile: Profile-general
-*Beam General Section, elset=Set-1, material=Material-Steel, section=GENERAL
-7.08, 82.8, 0., 18.3, 0.35
-0.,1.,0.
-** Section: Section-member  Profile: Profile-general
-*Beam General Section, elset=Set-4, material=Material-Steel, section=GENERAL
-7.08, 82.8, 0., 18.3, 0.35
-1.,0.,0.
-** Section: Section-member  Profile: Profile-general
-*Beam General Section, elset=Set-7, material=Material-Steel, section=GENERAL
-7.08, 82.8, 0., 18.3, 0.35
-0.,1.,0.
-*End Part
-**  
-**
-** ASSEMBLY
-**
-*Assembly, name=Assembly
-**  
-*Instance, name=Part-frame-1, part=Part-frame
-*End Instance
-**  
-*Nset, nset=Set-fixed, instance=Part-frame-1
- 1, 3, 4
-*Nset, nset=Set-loaded, instance=Part-frame-1
- 2,
-*End Assembly
-** 
-** MATERIALS
-** 
-*Material, name=Material-Steel
-*Elastic
- 2.9e+07, 0.3
-** ----------------------------------------------------------------
-** 
-** STEP: Step-statics
-** 
-*Step, name=Step-statics, nlgeom=NO, perturbation
-*Static
-** 
-** BOUNDARY CONDITIONS
-** 
-** Name: BC-fixed Type: Symmetry/Antisymmetry/Encastre
-*Boundary
-Set-fixed, ENCASTRE
-** 
-** LOADS
-** 
-** Name: Load-moment   Type: Moment
-*Cload
-Set-loaded, 5, 492000.
-** 
-** OUTPUT REQUESTS
-** 
-** 
-** FIELD OUTPUT: F-Output-1
-** 
-*Output, field
-*Node Output, variable=PRESELECT
-*Element Output, directions=YES
-LE, S, SF
-** 
-** HISTORY OUTPUT: H-Output-1
-** 
-*Output, history, variable=PRESELECT
-*End Step