三モーメント法 第12章 もう分類するのが面倒くさい 発展例題12.06 p377「材料力学演習(20221021)」をsympyでやってみたい。

・作図は未です。

オリジナル

発展例題12.06 p377＜「材料力学演習(20221021)一括(ver.3.3)」を勉強したい。#sympy
　　？？？タブレット等で、pdfを開く事ができないかも。？？？　

sympyで

# ver0.1
from sympy import *
var('w,a,EI')
var('MA,MB,MC,MD')
var('Mnp1,Mn,ln')
var('rn2,rnp11')
var('Rn2,Rnp11')
# θAn1    =w*a**3/(24*EI); θAn2=θAn1
MA=0;MD=0
rep_MBMC=solve([Eq(a*MA+4*a*MB+  a*MC     ,-Rational(1,2)*w*a**3),  \
                Eq(       a*MB+4*a*MC+a*MD,-Rational(1,2)*w*a**3)], \
          [MB,MC] \
         )   #;print(rep_MBMC)
MB=rep_MBMC[MB];MC=rep_MBMC[MC]
equ_1205a1_1= (Mnp1-Mn)/ln
equ_1205a1_2=-(Mnp1-Mn)/ln
RA1         =equ_1205a1_1.subs({Mnp1:MB,Mn:MA,ln:a})   #;print("#",RA1)
RB1         =equ_1205a1_1.subs({Mnp1:MC,Mn:MB,ln:a})   #;print("#",RB1)
RC1         =equ_1205a1_1.subs({Mnp1:MD,Mn:MC,ln:a})   #;print("#",RC1)
RA2         =-RA1
RB2         =-RB1
RC2         =-RC1
equ_1205d   =rn2+rnp11+Rn2+Rnp11
rep         ={rn2:Rational(1,2)*w*a,rnp11:Rational(1,2)*w*a}
RA          =(rn2+RA1).subs(rep)                        
RB          =equ_1205d.subs(rep).subs({Rn2:RA2,Rnp11:RB1})     ;print("#",RA,",",RB)
RC          =equ_1205d.subs(rep).subs({Rn2:RB2,Rnp11:RC1})    
RD          =(rn2+RC2).subs(rep)                               ;print("#",RC,",",RD)
# 2*a*w/5 , 11*a*w/10
# 11*a*w/10 , 2*a*w/5

いつもの? sympyの実行環境と 参考のおすすめです。

(テンプレート)

いつもと違うおすすめです。

参考文献

＞一部に等分布荷重を受ける単純支持はり
＞JSME p89
　
＞単純支持ばり 7)
＞構造力学公式集 p136