from sympy import *
from sympy.physics.continuum_mechanics.beam import Beam
E, I = symbols('E, I')
b = Beam(9, E, I)
b.apply_load(-12, 9, -1)
b.apply_load(-50, 5, -2)
b.apply_load(-8, 0, 0, end=5)
b.apply_support(0, type='fixed')
b.apply_support(4.5, type='roller')
b.load
after adding the above
I suppose you are referring to R4.5000. This is a way to extract the
forces/moments from the above equation:
fs = list(b.load.free_symbols)
R45, R0, M0 = fs[1:]
M0, R0, R45
Now you can create a substitution dictionary and compute the load at a
given location.
Davide.
Il giorno lun 5 lug 20
had to rush through my project so i jerry rigged this piece of code i know
its not the most elegant way but yeah, it works.
b.applied_loads #gives all the loads applied to the beam as a
list of tuples and extracted all the unknow reaction forces and appending
them to a list
named r
free_symbols is a set. You shouldn't assume it will come out in any
particular order. The order for the exact same code might not even be
the same across different Python sessions.
Aaron Meurer
On Mon, Jul 5, 2021 at 7:23 AM Davide Sandona' wrote:
>
> I suppose you are referring to R4.5000.