[例4]L字はり(ラーメン構造)
youtubeの板書です。
途中37:00から (0:00〜1:22:49)
sympyで
・ver0.1
# ver0.1
# [例4] L字はり(ラーメン構造)
from sympy import *
var('P,l,R,a,b,x,EA,EI',real=True)
U1=1/(2*EA)*integrate( P **2,(x,0,a))+1/(2*EI)*integrate((P*b)**2,(x,0,a)) #;print('#',U1)
U2=1/(2*EI)*integrate((-P*(b-x))**2,(x,0,b)) #;print('#',U2)
U1U2=U1+U2 #;print('#',U1U2)
δ =diff(U1U2,P) ;print('#',δ)
# P*a*b**2/EI + P*b**3/(3*EI) + P*a/EA
(勉強中) Mathematicaで
(勉強中) MATLAB 追加のオプションSymbolic Math Toolbox で
いつもの? sympyの実行環境と 参考のおすすめです。
いつもと違うおすすめです。
参考文献
>等分布荷重を受ける単純支持はり
>JSME p85
>単純支持ばり 7)
>構造力学公式集 p136