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

オリジナル

基本例題2.11 P38
http://zairikiweb.starfree.jp/zai_enshuh/zai_enshuh_ver.3.3.pdf#page=43

「材料力学演習(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 *
N1,N2  =symbols('N1  N2' ,real=True)
P      =symbols('P'      ,real=True)
l      =symbols('l'      ,real=True) 
A1,A2  =symbols('A1  A2' ,real=True)
EE1,EE2=symbols('EE1 EE2',real=True)
myTikara  =Eq(2*N1+N2      ,2*P)
myHenkaryo=Eq(N1*l/(A1*EE1),N2*l/(A2*EE2))
ans=solve([myTikara,myHenkaryo],(N1,N2))
print("#",ans)
#
N1=P*(A1*EE1)/(2*A1*EE1+A2*EE2)
N2=P*(A2*EE2)/(2*A1*EE1+A2*EE2)
print("#?????",(2*N1+N2).simplify())
print("#?????",(N1*l)/(A1*EE1)-(N2*l)/(A2*EE2) )
# {N1: 2*A1*EE1*P/(2*A1*EE1 + A2*EE2), N2: 2*A2*EE2*P/(2*A1*EE1 + A2*EE2)}
#????? P
#????? 0

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

Mathematicaで

勉強中

作図

勉強中