Example 8

%matplotlib inline
import sympy as sm
from sympy.physics.continuum_mechanics.beam import Beam

E = sm.Symbol('E')
I = sm.Symbol('I')
b = Beam(7, E, I)
r0, m0 = b.apply_support(0, type='fixed')
r4 = b.apply_support(4, type='pin')
r7 = b.apply_support(7, type='pin')
b.apply_rotation_hinge(2)
b.apply_rotation_hinge(5)
b.apply_load(-5, 0, 0, 2)
b.apply_load(-10, 6, -1)
b.solve_for_reaction_loads(r0, m0, r4, r7)

b.plot_bending_moment();