「(2018)平成30年二級建築士試験学科Ⅲ（建築構造）〔Ｎｏ．３〕」単純梁に集中荷重と一部等分布荷重をsympyのBeamでやってみた。

Last updated at 2025-08-04Posted at 2025-08-02

Qiita_二級建築士目次

オリジナル

・過去問.com 二級建築士様

sympyで(オリジナル様の方法を参考に)

ver0.1

# ver0.1
from sympy import *
var('x,Vb,Ma')             
Vb=solve(Eq(-2*10+Vb*8-12*2,0),Vb)[0] #;print(Vb)
Ma=solve(Eq(-2* 6+Vb*4+Ma  ,0),Ma)[0] #;print(Ma)
print("#",abs(Ma))
# 10

ver0.2

# ver0.2
from sympy import *
var('VI')             
VI=solve(Eq(-8*VI+2*10+12*2,0),VI)[0]
print("#",4*VI-6*2)
# 10

sympyで(Beamで)

ver0.7

・シンボリックで、していません。

# ver0.7
from sympy.physics.continuum_mechanics.beam import Beam
from sympy import *
var('x')             
var('R1,M1,R2,M2')         
def my_Overhange():
    P,w,l2,l4,E,I=2,3,2,4,1,1
    b=Beam(l2+l4+l4,E,I)
    b.apply_load(P ,0       ,-1)
    b.apply_load(R1,l2      ,-1)
    b.apply_load(w ,l2+l4   , 0,l2+l4+l4)
    b.apply_load(R2,l2+l4+l4,-1)
    b.bc_deflection=[(l2      ,0), \
                     (l2+l4+l4,0)]
    b.solve_for_reaction_loads(R1,R2)
    print("#",b.bending_moment().subs({x:l2+l4}))
    b.plot_shear_force()
    b.plot_bending_moment()
    b.remove_load  (R1,l2      ,-1)
    b.remove_load  (R2,l2+l4+l4,-1)
    b.apply_support(l2         ,"pin")
    b.apply_support(l2+l4+l4   ,"roller")
    b.draw().show()
    return 
my_Overhange()
# 10