静力学と材料力学 応用例題1.32 P23「材料力学演習(20221021)」をsympyでやってみた。 9部材Truss

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

オリジナル

応用例題1.32 P23 (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')
P25,P30,P50=25,30,50  # var('P1')
L3=3  # var('L1')
H4=4  # var('L1')
t = Truss()
t.add_node     (          ('E',L3*1,H4),('F',L3*2,H4)) 
t.add_node     (('A',0,0),('B',L3*1,0 ),('C',L3*2,0 ),('D',L3*3,0)) 
t.add_member   (('AE','A','E'))
t.add_member   (('AB','A','B'))   
t.add_member   (('BC','B','C'))   
t.add_member   (('BE','B','E'))   
t.add_member   (('EC','E','C'))   
t.add_member   (('EF','E','F'))   
t.add_member   (('CF','C','F')) 
t.add_member   (('CD','C','D'))   
t.add_member   (('DF','D','F'))   
t.apply_support(('A' ,'pinned'),('D','roller'))
t.apply_load   (('F',P50,180))
t.apply_load   (('B',P25,270))
t.apply_load   (('C',P30,270))
p = t.draw()
p.show()
t.solve()
print("#",{k:round(v,1) for k,v in t.reaction_loads .items()})           
print("#",{k:round(v,1) for k,v in t.internal_forces.items()})           
# {'R_A_x': 50, 'R_A_y': 48.9, 'R_D_y': 6.1}
# {'AE': -61.1, 'AB': -13.3, 'BC': -13.3, 'BE': 25, 'EC': 29.9, 'EF': -54.6, 'CF': 6.1, 'CD': 4.6, 'DF': -7.6}

ver1.3

・作成中