sympyで「薄肉円筒(Thin cylinder)」をやってみた。技術士機械令和元年度問Ⅲ-10

技術士第一次試験－令和元年度機械部門Ⅲ-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