弾性棒に係る不静定系の力学 基本例題2.02 P31 「材料力学演習(20221021)」をsympyでやってみた。

オリジナル

基本例題2.02 P31
http://zairikiweb.starfree.jp/zai_enshuh/zai_enshuh_ver.3.3.pdf#page=36

「材料力学演習(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で

オリジナル 様の方法を参考に

・ver0.1

# ver0.1
from sympy import *
var('P,N1,N2,A1,E1,A2,E2,l')
equ_a=Eq(N1+N2,P)
δ1   =N1*l/(A1*E1)
δ2   =N2*l/(A2*E2)
equ_b=Eq(δ1,δ2)
sol  =solve([equ_a,equ_b],[N1,N2])
print("#",sol[N1]   ,sol[N2])
print("#",sol[N1]/A1,sol[N2]/A2)
# A1*E1*P/(A1*E1 + A2*E2) A2*E2*P/(A1*E1 + A2*E2)
# E1*P/(A1*E1 + A2*E2) E2*P/(A1*E1 + A2*E2)

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

(テンプレート) 　

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

参考文献

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