・計算結果の表示があと少しです。
・作図は未です。
・タイトルを変更予定です。
オリジナル
発展例題6.36 p189＜「材料力学演習(20221021)一括(ver.3.3)」を勉強したい。#sympy
？？？タブレット等で、pdfを開く事ができないかも。？？？
sympyで
# ver0.1 [発展例題6.36 p189]
from sympy import *
var('P,MB,θA,θB,EI1,EI2,x1,x2,l1,l2,RA')
RA =MB/l1
MAB =RA*x1
vAB =θA*x1-1/EI1**1/6*RA*x1**3
θAB =diff(vAB,x1)
sol_θA =solve(Eq(vAB.subs({x1:l1}),0),θA)[0] #;print(sol_θA)
sol_θB =θAB.subs({x1:l1,θA:sol_θA}) #;print(sol_θB)
MBC =MB+Rational(1,2)*P*x2 #;print(MBC)
vBC =θB*x2-1/EI2*(Rational(1,2)*MB*x2**2+Rational(1,12)*P*x2**3) #;print(vBC)
θBC =diff(vBC,x2).subs({θB:sol_θB}) #;print(θBC)
sol_MB =solve( Eq(θBC.subs({x2:l2}),0),MB)[0] #;print(sol_MB)
sol_θB =sol_θB.subs({MB:sol_MB}) #;print(sol_θB)
vBC =vBC.subs({x2:l2,MB:sol_MB,θB:sol_θB}).simplify()
print("#",vBC)
com=Rational(1,24)*(l2**3*P)/EI2
print("#",com," * ",vBC/com)
# P*l2**3*(3*EI1*l2 + 4*EI2*l1)/(24*EI2*(3*EI1*l2 + EI2*l1))
# P*l2**3/(24*EI2) * (3*EI1*l2 + 4*EI2*l1)/(3*EI1*l2 + EI2*l1)
いつもの? sympyの実行環境と 参考のおすすめです。
いつもと違うおすすめです。
参考文献
＞一部に等分布荷重を受ける単純支持はり
＞JSME p89
＞単純支持ばり 7)
＞構造力学公式集 p136