ねじり 基本例題3.07 P63「材料力学演習(20221021)」をsympyでやってみた。

オリジナル

基本例題3.07 P63
http://zairikiweb.starfree.jp/zai_enshuh/zai_enshuh_ver.3.3.pdf#page=68

「材料力学演習(20221021)一括(ver.3.3)」を勉強したい。＞sympy
https://qiita.com/mrrclb48z/items/c3274daf5f7a99cfe027#%E4%BD%9C%E6%A5%AD%E4%B8%AD%E7%A7%81%E3%81%AE%E8%A7%A3%E7%AD%94sympy

sympyで

・ver0.1

# ver0.1 
from sympy import *
TC,  TD       =symbols('TC  TD '       ,real=True)
TAC, TCB      =symbols('TAC TCB'       ,real=True)
phiA,phiB,phiC=symbols('phiA,phiB,phiC',real=True)
G,Ip,x        =symbols('G Ip x'        ,real=True)
l1,l2         =symbols('l1 l2'         ,real=True)
myEqA =Eq(TC+TCB-TAC,0)
myEqB1=Eq(phiC-phiA,integrate(TAC/(G*Ip),(x,0 ,l1    )))
myEqB2=Eq(phiB-phiC,integrate(TCB/(G*Ip),(x,l1,l1+l2 )))
myEqB3=Eq(phiB-phiA,0)
myEqB4=Eq(myEqB2.lhs+myEqB1.lhs,myEqB2.rhs+myEqB1.rhs  )
#
myEqB3mB4=Eq(myEqB3.lhs-myEqB4.lhs,((myEqB3.rhs-myEqB4.rhs)*(G*Ip)).simplify())
myEqB    =Eq(myEqB3mB4.rhs*(-1),myEqB3mB4.lhs*(-1))
ans=solve([myEqA,myEqB],[TAC,TCB])
print("#",                     ans[TAC],    ans[TCB]   )
print("#",myEqB1.rhs.subs({TAC:ans[TAC],TCB:ans[TCB]}) )
# TC*l2/(l1 + l2) -TC*l1/(l1 + l2)
# TC*l1*l2/(G*Ip*(l1 + l2))

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

(テンプレート)

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

参考文献

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