Example 11

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

E = sm.Symbol('E')
I = sm.Symbol('I')
b = Beam(20, E, I)
r0 = b.apply_support(0, type="pin")
r20, m20 = b.apply_support(20, type="fixed")
b.apply_rotation_hinge(5)
b.solve_for_ild_reactions(-1, r0, r20, m20)

b.plot_ild_reactions();