sympyの 「Solving Beam Bending Problems using Singularity Functions」を調べた。

example-1

from sympy import *
x, y, z = symbols('x y z')
init_printing(use_unicode=True, wrap_line=False)
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.load
b.apply_load(50, 5, -2)
b.apply_load(8, 0, 0, end=5)
b.bc_deflection.append((0, 0))
b.bc_slope.append((0, 0))
# 
R, M = symbols('R, M')
b.apply_load(R, 0, -1)
b.apply_load(M, 0, -2)
b.load
b.solve_for_reaction_loads(R, M)
b.reaction_loads                      ;print("# b.reaction_loads =",b.reaction_loads)
b.load                                ;print("# b.load           =",b.load)

# b.shear_force()
# b.bending_moment()

# b.plot_shear_force()  
# b.plot_bending_moment()  

いつもの? sympyの実行環境と 参考のおすすめです。

(テンプレート)

いつもと違うおすすめです。

＞一部に等分布荷重を受ける単純支持はり
＞JSME p89
　
＞単純支持ばり 10)
＞構造力学公式集 p138