ひずみエネルギとそれを使った材料力学 基本例題8.27 P284 「材料力学演習(20221021)」をsympyでやってみたい。

基本例題8.27解説へつづきます。
https://qiita.com/mrrclb48z/items/20834c88192e10a315d8

オリジナル

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

sympyで

・ ver0.1 p284 基本例題8.27

# ver0.1 p284 基本例題8.27
from sympy import *
var('q,l')
var('x,P,X')
MAX_str ="1/2*q*l*x-1/2*q*x**2+(l-X)/l*P*x"                                             ;print("#",MAX_str)
MXB_str = MAX_str                         + "-P*(x-X)"                                  ;print("#",MXB_str)
U       = "1/(2*EI)*integrate((" + MAX_str+")**2,(x,0,X))+" \
         +"1/(2*EI)*integrate((" + MXB_str+")**2,(x,X,l))"                             #;print("#",U)    
σX      =     "1/EI*integrate("  + str(diff(sympify(MAX_str)**2,P)/2) + ",(x,0,X))+" \
              "1/EI*integrate("  + str(diff(sympify(MXB_str)**2,P)/2) + ",(x,X,l))"    #;print("#",σX)    
rep={P:0}
σX      =factor(sympify(σX).subs(rep))                                                  ;print("#",σX)    
# 1/2*q*l*x-1/2*q*x**2+(l-X)/l*P*x
# 1/2*q*l*x-1/2*q*x**2+(l-X)/l*P*x-P*(x-X)
# X*q*(-X + l)*(-X**2 + X*l + l**2)/(24*EI)

sympyのBeamで

勉強中

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

(テンプレート)

いつもと違うおすすめです。

参考文献

＞一部に等分布荷重を受ける単純支持はり
＞JSME p89
　
＞単純支持ばり 7)
＞構造力学公式集 p136