More than 1 year has passed since last update.
sympyで「薄肉円筒(Thin cylinder)」をやってみた。技術士機械令和元年度問Ⅲ-10
1
Last updated at Posted at 2020-09-08
技術士第一次試験－令和元年度機械部門Ⅲ-10
＜ 公式ホームページ
（参考）技術士第一次試験－令和元年度機械部門Ⅲ-10解答例
入力
from sympy import *
var('σt σz d t p')
s = solve(Eq(p*d, 2*t*σt), σt)
print("σt=",s)
print("σt=",Float(s[0].subs({d:370,t:2.5,p:3.0}),4))
var('σz π t r p')
s = solve(Eq(π*(d/2)**2*p, σz*π*d*t), σz)
print("σz=",s)
print("σz=",Float(s[0].subs({d:370,t:2.5,p:3.0}),4))
出力
σt= [d*p/(2*t)]
σt= 222.0
σz= [d*p/(4*t)]
σz= 111.0
Register as a new user and use Qiita more conveniently
- You get articles that match your needs
- You can efficiently read back useful information
- You can use dark theme