・Function の使用は、はじめてかも。
・作図未です。
・教科書の例題3.3 p37 自重による伸びを求めよ。「図3.8 自重を受ける棒の伸び」。
オリジナル
基本例題8.07 P270<「材料力学演習(20221021)一括(ver.3.3)」を勉強したい。#sympy
???タブレット等で、pdfを開く事ができないかも。???
sympyで
オリジナル 様の方法を参考に
・ver0.1
# ver0.1
from sympy import *
var('E,P,σx,ρg,x,l,A')
u=(1/(2*E)*σx**2).subs({σx:ρg*x}) ;print("#",u)
U=integrate(u*A,(x,0,l)) ;print("#",U)
A=Function('A')
U=Integral (u*A(x),(x,0,l)).factor() ;print("#",U)
# x**2*ρg**2/(2*E)
# A*l**3*ρg**2/(6*E)
# ρg**2*Integral(x**2*A(x), (x, 0, l))/(2*E)
いつもの? sympyの実行環境と 参考のおすすめです。
いつもと違うおすすめです。
・All functions support the methods documented below, inherited from sympy.core.function.Function.
参考文献
>一部に等分布荷重を受ける単純支持はり
>JSME p89
>単純支持ばり 7)
>構造力学公式集 p136