# Beam.get_reactions would be wrong.
reactions = [PointLoad(1, 0), PointLoad(1, a + b + c)]
# Initialize the beam object.
# For now we don't need the actual value of Jz since we are
# going to compute the actual dimensions based on the values
# of reaction forces.
beam = Beam(
loads,
reactions,
l=a + b + c,
E=E,
Jz=1.,
)
Rs = beam.get_reactions()
print('reactions:\n R_A = {:.2e}N\n R_B = {:.2e}N'.format(Rs[0], Rs[1]))
### computation and plots of bending force and moment ###
n = 12 * 5 + 1
x = np.linspace(0, a + b + c, n)
plt.plot(x, [beam.force(xi) for xi in x], '.-', label='force')
Mo = np.array([beam.moment(xi) for xi in x])
plt.plot(x, Mo, '.-', label='moment')
plt.legend(loc='best')
plt.grid()
### actual dimensions of the cross-section ###
# reactions - list of unknown loads
# All must have magnitude = 1, otherwise the results of
# Beam.get_reactions would be wrong.
reactions = [PointLoad(1, 0), PointLoad(1, a+b+c)]
# Initialize the beam object.
# For now we don't need the actual value of Jz since we are
# going to compute the actual dimensions based on the values
# of reaction forces.
beam = Beam(loads, reactions,
l = a+b+c,
E = E,
Jz = 1.,
)
Rs = beam.get_reactions()
print('reactions:\n R_A = {:.2e}N\n R_B = {:.2e}N'.format(Rs[0], Rs[1]))
### computation and plots of bending force and moment ###
n = 12*5+1
x = np.linspace(0, a+b+c, n)
plt.plot(x, [beam.force(xi) for xi in x], '.-', label='force')
Mo = np.array([beam.moment(xi) for xi in x])
plt.plot(x, Mo, '.-', label='moment')
plt.legend(loc='best')
plt.grid()
### actual dimensions of the cross-section ###
