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

オリジナル

基本例題2.01 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で

from sympy import *
NAB,NBC=symbols('NAB NBC',real=True)
PB     =symbols('PB'     ,real=True)
A  ,EE =symbols('A   EE' ,real=True)
l1 ,l2 =symbols('l1  l2' ,real=True)
myBsyuhen=Eq(-NAB-PB+NBC,0)
dAB=NAB*l1/(A*EE)
dBC=NBC*l2/(A*EE)
myHenkei=Eq(dAB+dBC,0)
ans=solve([myBsyuhen,myHenkei],[NAB,NBC])
print("#",ans)
mySubs={PB:7,l1:0.7,l2:0.3}
print("#",round(ans[NAB].subs(mySubs),1) )
print("#",round(ans[NBC].subs(mySubs),1) )
# {NAB: -PB*l2/(l1 + l2), NBC: PB*l1/(l1 + l2)}
# -2.1
# 4.9

①私の環境は,pycharmです。
②よく聞くのは、Jupyterです。
③web上で、上記のソースを「SymPy Live shell」に、コピー貼り付けでもできました。
黒背景の右上に、マウスを移動すると、コピーマークが発生します。
???　タブレット環境で、コピー貼り付けが実行できませんでした。???

Mathematicaで

勉強中

作図

勉強中