オリジナル
基本例題1.22 P16
http://zairikiweb.starfree.jp/zai_enshuh/zai_enshuh_ver.3.3.pdf#page=21
「材料力学演習(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 *
a ,b =symbols('a b' ,real=True)
t =symbols('t' ,real=True)
P =symbols('P' ,real=True)
T =symbols('T' ,real=True)
RAH,RAV=symbols('RAH,RAV',real=True)
myFx=Eq( RAH-T*cos(t) ,0)
myFy=Eq( RAV+T*sin(t)-P ,0)
myM =Eq((a+b)*T*sin(t)-a*P,0) # (Aまわり)
ans=solve([myFx,myFy,myM],[T,RAH,RAV])
print("#T ",ans[T ].simplify())
print("#RAH",ans[RAH].simplify())
print("#RAv",ans[RAV].simplify())
#T P*a/((a + b)*sin(t))
#RAH P*a/((a + b)*tan(t))
#RAv P*b/(a + b)
いつもの? sympyの実行環境と 参考のおすすめです。
いつもと違うおすすめです。
参考文献
>一部に等分布荷重を受ける単純支持はり
>JSME p89
>単純ばり 2)
>構造力学公式集 表5.2 単純ばりの公式 p134