静力学と材料力学 演習問題6 p28「材料力学演習(20221021)」をsympyでやってみた。 7部材Truss

・sympyのTrussで作図しました。縦横比が1:1でありません。

オリジナル

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

sympyで

・勉強中

sympyのTRussで

# ver0.1 
from sympy.physics.continuum_mechanics.truss import Truss
import math
from sympy import *
var('P')
P1=1  # var('P1')
L1=1  # var('L1')
t = Truss()
t.add_node     (('E',0,2*L1)) 
t.add_node     (('C',0,1*L1),('D',Rational(3,4)*L1,1*L1)) 
t.add_node     (('A',0,0   ),('B',Rational(3,4)*L1,0  )) 
t.add_member   (('AB','A','B'))   
t.add_member   (('AC','A','C'),('BC','B','C'),('BD','B','D'))   
t.add_member   (('CE','C','E'),('CD','C','D'),('DE','D','E'))   
t.add_member   ()   
t.apply_support(('A' ,'pinned'),('B','roller'))
t.apply_load   (('E',P1,0.0))
t.apply_load   (('C',P1,0.0))
p = t.draw()
p.show()
t.solve()
print("#",{k:P*v for k,v in t.reaction_loads .items()})           
print("#",{k:P*v for k,v in t.internal_forces.items()})           
# {'R_A_x': -2*P, 'R_A_y': -4*P, 'R_B_y': 4*P}
# {'AB': 2*P, 'AC': 4*P, 'BC': -10*P/3, 'BD': -4*P/3, 'CE': 4*P/3, 'CD': P, 'DE': -5*P/3}

ver1.0

・Trussの 「property member_lengths」がありました。
　作成中です。