--
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+unsubscribe@googlegroups.com.
To post to this group, send email to sy...@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/ad2c8be6-7bee-4bff-9c60-6aed1a8a5e50%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
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)
--
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+unsubscribe@googlegroups.com.
To post to this group, send email to sy...@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.
>>> from sympy.physics.continuum_mechanics.beam import Beam
>>> from sympy import symbols, Piecewise
>>> E, I = symbols('E, I')
>>> R1, R2, R3, R4 = symbols('R1, R2, R3, R4')
>>> b = Beam(4, E, I)
>>> b.apply_load(R1, 0, -1)
>>> b.apply_load(R2, 0, -2)
>>> b.apply_load(R3, 4, -1)
>>> b.apply_load(R4, 4, -2)
>>> b.apply_load(6, 2, -1)
>>> b.bc_deflection = [(0, 0), (4, 0)]
>>> b.bc_slope = [(0,0), (4,0)]
>>> b.solve_for_reaction_loads(R1,R2,R3,R4)
>>>b.shear_force()
R3*SingularityFunction(x, 4, 0) + R4*SingularityFunction(x, 4, -1) + (-R3 - 6)*SingularityFunction(x, 0, 0) + (4*R3 - R4 + 12)*SingularityFunction(x, 0, -1) + 6*SingularityFunction(x, 2, 0)
>>> b.bending_moment()
R3*SingularityFunction(x, 4, 1) + R4*SingularityFunction(x, 4, 0) + (-R3 - 6)*SingularityFunction(x, 0, 1) + (4*R3 - R4 + 12)*SingularityFunction(x, 0, 0) + 6*SingularityFunction(x, 2, 1)
--
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+unsubscribe@googlegroups.com.
To post to this group, send email to sy...@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/67bd6cfd-eeaa-4463-adbf-3def6b427ac9%40googlegroups.com.
--
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+unsubscribe@googlegroups.com.
To post to this group, send email to sy...@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/f815cf7c-ea89-4c4c-adec-3ee911395548%40googlegroups.com.