In this PR <https://github.com/sympy/sympy/pull/11374> Sampad Saha 
mentioned to check if the beam bending works with units or not.

Currently it doesn't work and I think it requires some changes in the units 
module as well as in continuum_mechanics package.

For the code snippet

from sympy import *
from sympy.physics.units import meter, kilogram, newton, second, kilo
from sympy.physics.continuum_mechanics.beam import Beam
E, I = symbols('E, I')
R1, R2 = symbols('R1, R2')
b = Beam(3*meter, E*newton/meter**2, I*meter**4)
b.apply_load(-8*kilo*newton, 0, -1)
b.apply_load(R1*kilo*newton, 1, -1)
b.apply_load(R2*kilo*newton, 3, -1)
b.apply_load(12*kilo*newton, 3, -2)
b.bc_deflection = [(1, 0), (3, 0)]
b.solve_for_reaction_loads(R1, R2)
print(b.reaction_loads)

 It currently gives ValueError 

-- 
You received this message because you are subscribed to the Google Groups 
"sympy" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to sympy+unsubscr...@googlegroups.com.
To post to this group, send email to sympy@googlegroups.com.
Visit this group at https://groups.google.com/group/sympy.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/sympy/3b72dbfd-43e6-4ea7-82f0-c44819189a55%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

Reply via email to