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.