「建築士H27構No.3【構造力学】#1 単純ばりの反力」をsympyのBeamでやってみた。

・タイトルを変更予定です。
・同じ？答えがでました。
・専門家様のアドバイスをいただければ幸いです。
・どなたかAIで実行してほしいです。

オリジナル

【構造力学】#1 単純ばりの反力
スチのり 様(0:00〜12:18)

普通に？

勉強中

sympyのBeamで

# H27構No.3
# ver0.1
from sympy.physics.continuum_mechanics.beam import Beam
from sympy import symbols, Piecewise
P , E, I = symbols('P,E,I')
VA, VB = symbols('VA, VB')
b = Beam(6, E, I)
b.apply_load(VA,0,-1)
b.apply_load(-P,2,-1)
b.apply_load(VB,6,-1)
b.bc_deflection = [(0, 0), (6, 0)]
b.boundary_conditions
b.load
b.solve_for_reaction_loads(VA, VB)
b.load
b.shear_force()
b.bending_moment()
#b.slope()
#b.deflection()
#b.deflection().rewrite(Piecewise)
#p = b.draw()  
#p
print("#VA=",b.reaction_loads[VA])
print("#VB=",b.reaction_loads[VB])
#VA= 2*P/3
#VB= P/3

　
https://live.sympy.org/

参考

sympyのdoc

Qiita