静力学と材料力学 基本例題1.06 P3「材料力学演習(20221021)」をsympyでやってみた。

２通りで

オリジナル

基本例題1.06 P3
http://zairikiweb.starfree.jp/zai_enshuh/zai_enshuh_ver.3.3.pdf#page=9

「材料力学演習(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で(pdfの手順で実行)

from sympy import *
PA, PB, PC=symbols('PA  PB  PC',real=True)
NAB,NBC   =symbols('NAB NBC   ',real=True)
#
myAF=Eq(-PA +NAB   ,0)
myBF=Eq(-NAB-PB+NBC,0)
myCF=Eq(-NBC+PC    ,0)
myABCF=Eq( myAF.lhs+myBF.lhs+myCF.lhs,myAF.rhs+myBF.rhs+myCF.rhs)
print("#",myABCF)
PA_ans=solve(myABCF,PA)[0]
print("#",PA_ans)
NAB_ans=solve(myAF,NAB)[0].subs({PA:PA_ans})
NBC_ans=solve(myCF,NBC)[0].subs({PA:PA_ans})
print("#",NAB_ans,NBC_ans)
mySubs={PB:7,PC:10}
print("#",NAB_ans.subs(mySubs),NBC_ans.subs(mySubs))
# Eq(-PA - PB + PC, 0)
# -PB + PC
# -PB + PC PC
# 3 10

sympyで(solveを1回で実行)

from sympy import *
PA, PB, PC=symbols('PA  PB  PC',real=True)
NAB,NBC   =symbols('NAB NBC   ',real=True)
myAF=Eq(-PA+NAB    ,0)
myBF=Eq(-NAB-PB+NBC,0)
myCF=Eq(-NBC+PC    ,0)
ans=solve([myAF,myBF,myCF],[PA,NAB,NBC])
print("#",ans)
mySubs={PB:7,PC:10}
print("#",ans[NAB].subs(mySubs))
print("#",ans[NBC].subs(mySubs))
# {PA: -PB + PC, NAB: -PB + PC, NBC: PC}
# 3
# 10

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

(テンプレート)

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

Mathematicaで

勉強中

作図

勉強中

Qitta内

参考文献

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