静力学と材料力学 演習問題1,2,3 P27 「材料力学演習(20221021)」をsympyでやってみた。

オリジナル

演習問題1,2,3 P27
http://zairikiweb.starfree.jp/zai_enshuh/zai_enshuh_ver.3.3.pdf#page=32

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

？？？？？四捨五入　#問題2. dAB= 44.0000000000000

# p54 ensyu mondai1
from sympy import *
NAB,NBC,NCD =symbols('NAB NBC NCD',real=True)
P1,P2       =symbols('P1  P2'     ,real=True)
RA,RB       =symbols('RA  RB'     ,real=True)
A  ,EE      =symbols('A   EE'     ,real=True)
lAB,lBC,lCD =symbols('lAB lBC lCD',real=True)
RA ,RD      =symbols('RA  RD'     ,real=True)
myAsyuhen=Eq( RA +NAB   ,0)
myBsyuhen=Eq(-NAB+P1+NBC,0)
myCsyuhen=Eq(-NBC-P2+NCD,0)
myDsyuhen=Eq( NCD+RD    ,0)
dAB=NAB*lAB/(A*EE)
dBC=NBC*lBC/(A*EE)
dCD=NCD*lCD/(A*EE)
myHenkei=Eq(dAB+dBC+dCD,0)
ans=solve([myAsyuhen,
           myBsyuhen,
           myCsyuhen,
           myDsyuhen,
           myHenkei],[RA,NAB,NBC,NCD,RD])
mySubs1={P1 :15 ,P2 :10,
         lAB:0.6,lBC:1.0,lCD:1.4,
         A  :500*10**(-6),EE:200*10**(+6)}
RA__ans=ans[RA ].subs(mySubs1)
NAB_ans=ans[NAB].subs(mySubs1)
NBC_ans=ans[NBC].subs(mySubs1)
NCD_ans=ans[NCD].subs(mySubs1)
RD__ans=ans[RD ].subs(mySubs1)
print("#問題1. RA =",round(RA__ans,2))
print("#問題1. RD =",round(RD__ans,2))
mySubs2={RA :RA__ans,
         NAB:NAB_ans,NBC:NBC_ans,NCD:NCD_ans,
         RD :RD__ans}
print("#問題2. NAB=",round(NAB_ans,2))
print("#問題2. NBC=",round(NBC_ans,2))
print("#問題2. NCD=",round(NCD_ans,2))
dAB__ans=dAB.subs(mySubs1).subs(mySubs2)
dCD__ans=dCD.subs(mySubs1).subs(mySubs2)
print("#問題2. dAB=",round(dAB__ans*10**6,2))
print("#問題2. dCD=",round(dCD__ans*10**6,1))
print('#############################################')
from sympy import *
NAB,NBC,NCD =symbols('NAB NBC NCD',real=True)
P1, P2      =symbols('P1  P2'     ,real=True)
RA, RB      =symbols('RA  RB'     ,real=True)
A , EE      =symbols('A   EE'     ,real=True)
A50         =symbols('A50'        ,real=True)
lAB,lBC,lCD =symbols('lAB lBC lCD',real=True)
RA ,RD      =symbols('RA  RD'     ,real=True)
myAsyuhen=Eq( RA +NAB   ,0)
myBsyuhen=Eq(-NAB+P1+NBC,0)
myCsyuhen=Eq(-NBC-P2+NCD,0)
myDsyuhen=Eq( NCD+RD    ,0)
dAB=NAB*lAB/(A*EE)
dBC=NBC*lBC/(A50*EE)
dCD=NCD*lCD/(A*EE)
myHenkei=Eq(dAB+dBC+dCD,0)
ans=solve([myAsyuhen,
           myBsyuhen,
           myCsyuhen,
           myDsyuhen,
           myHenkei],[RA,NAB,NBC,NCD,RD])
mySubs1={P1 :15 ,P2 :10,
         lAB:0.6,lBC:1.0,lCD:1.4,
         A  :500*10**(-6),EE:200*10**(+6),
         A50:500*10**(-6)*0.5}
RA__ans=ans[RA ].subs(mySubs1)
NAB_ans=ans[NAB].subs(mySubs1)
NBC_ans=ans[NBC].subs(mySubs1)
NCD_ans=ans[NCD].subs(mySubs1)
RD__ans=ans[RD ].subs(mySubs1)
print("#問題3. RA =",round(RA__ans,2))
print("#問題3. RD =",round(RD__ans,2))
#問題1. RA = -7.33
#問題1. RD = -2.33
#問題2. NAB= 7.33
#問題2. NBC= -7.67
#問題2. NCD= 2.33
#問題2. dAB= 44.0000000000000
#問題2. dCD= 32.7
#############################################
#問題3. RA = -9.25
#問題3. RD = -4.25

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

Mathematicaで

勉強中

作図

勉強中