ひずみエネルギとそれを使った材料力学 基本例題8.18 解説3(途中) P277 「材料力学演習(20221021)」をsympyでやってみたい。

　　　基本例題8.18 P276　　　　　　(2024-11-14)
　　　基本例題8.18 解説1,2 P276　　(2024-11-14)

オリジナル

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

sympyで

・ ver0.1 p277 基本例題8.18 解説3
・変数Eを変数myEにしました。
・(途中) 申し訳ありません。後半は、理解できませんでした

# ver0.1 p276 基本例題8.18 解説3
#　p276 基本例題8.18 解説3
from sympy import *
var('myE,G,Iz,Ip' )  # var('E,G,Iz,Ip' )
var('P,a,b,c'   )
var('x,x1,x2,x3,tA')
def myStrFactor(expr,myStr):
    return myStr+"* ("+str( (expr/sympify(myStr)).expand())+")"
str_d2vAB ="-1/(myE*Iz)*(P*c-P*a+P*x3)"                                        ;print("#",str_d2vAB)
str_tAB   =myStrFactor(integrate(sympify(str_d2vAB),x3),"-P/(myE*Iz)")+" + tA" ;print("#",str_tAB)
tAB       =sympify(str_tAB)
tA_sol    =solve  (Eq(tAB.subs({x3:0}),0),tA)[0]                               ;print("#",tA_sol)
tB        =tAB.subs({tA:tA_sol}).subs({x3:a})                                  ;print("#",factor(tB))
#
# -1/(myE*Iz)*(P*c-P*a+P*x3)
# -P/(myE*Iz)* (-a*x3 + c*x3 + x3**2/2) + tA
# 0
# P*a*(a - 2*c)/(2*Iz*myE)

・(途中)

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

(テンプレート)

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

参考文献

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