Example 4

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

E = sm.Symbol('E')
I = sm.Symbol('I')
b1 = Beam(5, E, I)
b2 = Beam(3, E, I)
b = b1.join(b2, via="hinge")
b.apply_support(0, type='fixed')
b.apply_support(8, type='pin')
b.apply_load(-10, 5, -1)
r0, m0, r8 = sm.symbols('R_0, M_0, R_8')
b.solve_for_reaction_loads(r0, m0, r8)

b.plot_bending_moment();