静力学と材料力学 発展例題1.23 P17「材料力学演習(20221021)」をsympyでやってみた。

オリジナル

発展例題1.23 P17
http://zairikiweb.starfree.jp/zai_enshuh/zai_enshuh_ver.3.3.pdf#page=22

「材料力学演習(20221021)一括(ver.3.3)」を勉強したい。＞sympy
https://qiita.com/mrrclb48z/items/c3274daf5f7a99cfe027#%E4%BD%9C%E6%A5%AD%E4%B8%AD%E7%A7%81%E3%81%AE%E8%A7%A3%E7%AD%94sympy

sympyで

from sympy import *
l1,l2,l3=symbols('l1,l2,l3',real=True)
P      =symbols('P'        ,real=True)
T      =symbols('T'        ,real=True)
RAH,RAV=symbols('RAH,RAV'  ,real=True)
RBH,RBV=symbols('RBH,RBV'  ,real=True)
myFx=Eq( RAH-RBH                ,0)
myFy=Eq( RAV+RBV-P              ,0)
myM =Eq((l1+l2)*RBV-(l1+l2+l3)*P,0)          # (Aまわり)
ans=solve([myFx,myFy,myM],[RAH,RAV,RBH,RBV])
# print(ans)
mySubs={P:500,l1:4,l2:2,l3:1}
RBV_ans=ans[RBV].subs(mySubs)
print("#",round(RBV_ans), round(ans[RAV].subs(mySubs)),ans[RAH].simplify())
# 
myCM=Eq(2*RBV-4*RBH,0) 
RBH_ans=solve(myCM,RBH)[0].subs({RBV:RBV_ans})
RAH_ans=RBH_ans.subs({RBV:RBV_ans})
RCH_ans=RBH_ans.subs({RBV:RBV_ans})
RCV_ans=RBV_ans
print("#",round(RBH_ans))
print("#",round(RCH_ans),round(RCV_ans))
# 583 -83 RBH
# 292
# 292 583

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

(テンプレート)

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

参考文献

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