・垂直方向変位,水平方向変位
オリジナル
基本例題11.02 P352<「材料力学演習(20221021)一括(ver.3.3)」を勉強したい。#sympy
???タブレット等で、pdfを開く事ができないかも。???
sympyで
・ver0.1
# ver0.1
# 基本例題11.02 解説 P351
from sympy import *
var('P,EI,l,a,θ,Q',real=True)
M =-P*a*(1-sin(θ))+Q*a*cos(θ)
δBV =(1/EI*a**2*integrate(M* cos(θ) ,(θ,0,pi/2) )) .expand() ;print("#",δBV )
Q_ans=solve(δBV,Q)[0] ;print("#",Q_ans)
δBH =1/EI*a**2*integrate(-M*(1-sin(θ)),(θ,0,pi/2) ).factor() ;print("#",δBH )
δBH = δBH.subs({Q:Q_ans}).factor() #;print("#",δBH )
δBH =(δBH*(1*pi*EI)/(P*a**3*pi)).expand()*(P*a**3*pi/(1*pi*EI)) ;print("#",δBH )
# -P*a**3/(2*EI) + pi*Q*a**3/(4*EI)
# 2*P/pi
# a**3*(-8*P + 3*pi*P - 2*Q)/(4*EI)
# P*a**3*(-2 - 1/pi + 3*pi/4)/EI
いつもの? sympyの実行環境と 参考のおすすめです。
いつもと違うおすすめです。
参考文献
>一部に等分布荷重を受ける単純支持はり
>JSME p89
>単純支持ばり 7)
>構造力学公式集 p136