はりに一つ以上の集中荷重。エネルギーで「youtube【東京理科大学】第11回 材料力学基礎授業[例3]」をsympyでやってみたい。

[例3]はりに一つ以上の集中荷重。エネルギーで

youtubeの板書です。

途中24:00から (0:00〜1:22:49)

sympyで

・ver0.1

# ver0.1
#  [例3] 
from sympy import *
var('P,l,R,a,x,EI',real=True)
M1=-P*(a-x)+R*(l-x)                                                      #;print('#',M1)
M2= R*(l-x)                                                              #;print('#',M2)
U =1/(2*EI)*integrate(M1**2,(x,0,a))+1/(2*EI)*integrate(M2**2,(x,a,l))   #;print('#',U )
δB=diff(U,R)                                                            #;print('#',δB )                                              
R_ans=solve(δB.subs({x:0}),R)[0]                                         ;print('#', R_ans)
# P*a**2*(-a + 3*l)/(2*l**3)

(勉強中) Mathematicaで

(勉強中) MATLAB 追加のオプションSymbolic Math Toolbox で

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

(テンプレート)

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

参考文献

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