貯水槽の壁の厚さt 第10章 薄肉圧力容器 基本例題10.02 p335「材料力学演習(20221021)」をsympyでやってみたい。(英:thin wall structure）

パイソニスタの方へ
・solveで、不等式を教えて下さい。
・Inequality Solversの使い方を教えて下さい。
https://docs.sympy.org/latest/modules/solvers/inequalities.html

オリジナル

基本例題10.02 p335＜「材料力学演習(20221021)一括(ver.3.3)」を勉強したい。#sympy
　　　　　　　　？？？タブレット等で、pdfを開く事ができないかも。？？？　

sympyで

from sympy import *

var('σa,σt'  ,real=True)
var('Q,H'    ,real=True)
var('D,t,p,g',real=True)
var('x'      ,real=True)
eq=Eq(σt,D/(2*t)*p)                 ; print("#",eq)
eq=Eq(σt,D/(2*t)*p).subs({p:Q*g*H}) ; print("#",eq)
lt=Lt(σt,σa)                        ; print("#",lt)
lt=Lt(σt,σa).subs({σt:eq.rhs})      #; print("#",lt)
# print( solve(lt,t) )
lt=Lt(σt,σa).subs({σt:eq.rhs})      #; print("#",lt)
lt=Lt(lt.lhs*t/σa,lt.rhs*t/σa)      #; print("#",lt)
gt=Gt(lt.rhs,lt.lhs)                ; print("#",gt)
# Eq(σt, D*p/(2*t))
# Eq(σt, D*H*Q*g/(2*t))
# σt < σa
# t > D*H*Q*g/(2*σa)

いつもの? sympyの実行環境と 参考のおすすめです。

(テンプレート)