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

基本例題8.27のつづきです。
https://qiita.com/mrrclb48z/items/d7e3fc5426d4ce132ea0

オリジナル

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

sympyで

・ver0.1 p284 基本例題8.27 解説
・??? 単位荷重法 ステキ
・NameError: name 'replace' is not defined

# ver0.1 p284 基本例題8.27 解説
from sympy import *
var('q,l')
var('x,P,X')
kou_str1="1/2*q*l*x-1/2*q*x**2"                         ;print("#",kou_str1)
kou_str2="(l-X)/l*P*x"                                  ;print("#",kou_str2)
kou1    =   "("+ kou_str1 +")**2"                       ;print("#",kou1)
kou2    = "2*("+ kou_str1 +")*" +  "("+ kou_str2 +")"   ;print("#",kou2)
kou3    =   "("+ kou_str2 +")**2"                       ;print("#",kou3)
rep0={P:0}
kou1diff=diff(kou1,P) .subs(rep0)                       ;print("#",kou1diff)
kou2diff=kou2.replace("1*P","")                         ;print("#",kou2diff)
kou3diff=diff(kou3,P) .subs(rep0)                       ;print("#",kou3diff)
# 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)**2
# 2*(1/2*q*l*x-1/2*q*x**2)*((l-X)/l*P*x)
# ((l-X)/l*P*x)**2
# 0
# 2*(1/2*q*l*x-1/2*q*x**2)*((l-X)/l*P*x)
# 0

途中です。

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

(テンプレート)

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

参考文献

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