オリジナル
基本例題3.01 P59 (pdf)<「材料力学演習(20221021)一括(ver.3.3)」を勉強したい。#sympy
???タブレット等で、pdfを開く事ができないかも。???
sympyで
# ver0.1
from sympy import *
var('T,D,l,r,G,Ip')
rep1={l:1,T:10,D:20*10**(-3),G:25}
rep2={Ip:(pi/32*D**4).subs(rep1)} ;print("#",round(rep2[Ip]*10**9,1))
τ =round((T/Ip).subs(rep1).subs(rep2)*10**(-6))
print("#",τ,"*r")
print("#",round((τ*r/G).subs({r:D/2,G:25*10**9}).subs(rep1)*10**12),"*10**(-6)")
φ=round((T*l/(G*Ip)).subs(rep1).subs(rep2)*10**(-6),1)
print("#",φ,"*10**(-3)")
print("#",round((φ*10**(-4)*180/pi),3))
# 15.7
# 637 *r
# 255 *10**(-6)
# 25.5 *10**(-3)
# 0.146
いつもの? sympyの実行環境と 参考のおすすめです。
いつもと違うおすすめです。
参考文献
>一部に等分布荷重を受ける単純支持はり
>JSME p89
>単純ばり 2)
>構造力学公式集 表5.2 単純ばりの公式 p134