・計算結果の表示があと少しです。
・作図は未です。
・タイトルを変更予定です。

オリジナル

発展例題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~~