Allow length greater than 10m

beambending/beam.py

@@ -304,7 +304,7 @@ def _plot_analytical(self, ax: plt.axes, sym_func, title: str = "", maxmin_hline
         :return: a matplotlib.Axes object representing the plotted data.
 
         """
-        x_vec = np.linspace(self._x0, self._x1, min(int((self.length) * 1000 + 1), 1e4))
+        x_vec = np.linspace(self._x0, self._x1, int(self.length * 1000 + 1))
         y_lam = lambdify(x, sym_func, "numpy")
         y_vec = np.array([y_lam(t) for t in x_vec])
 

examples/example_4.py

@@ -0,0 +1,28 @@
+# -*- coding: utf-8 -*-
+"""Example 4: Similar to example 1.
+"""
+import os
+from beambending import Beam, DistributedLoadV, PointLoadH, PointLoadV, x
+
+def example_4():
+    """Run example 4"""
+    beam = Beam(20)  # Initialize a Beam object of length 9 m
+    beam.pinned_support = 2    # x-coordinate of the pinned support
+    beam.rolling_support = 18  # x-coordinate of the rolling support
+    beam.add_loads((
+                    PointLoadH(10, 3),  # 10 kN pointing right, at x=3 m
+                    PointLoadV(-20, 15),  # 20 kN downwards, at x=3 m
+                    DistributedLoadV(-10, (0, 20)),  # 10 kN/m, downwards, for 0 m <= x <= 20 m
+                    DistributedLoadV(-20 + x**2, (0, 11)),  # variable load, for 0 <= x <= 11 m
+                ))
+    fig = beam.plot()
+
+    # save the png and add it to the documentation
+    mod_path = os.path.dirname(os.path.abspath(__file__))  # current module
+    save_name = os.path.basename(__file__).replace('.py', '.png')  # file name
+    save_path = os.path.join(mod_path, save_name)
+    fig.savefig(save_path, transparent=True)
+
+
+if __name__ == '__main__':  # call function when run as script
+    example_4()