l = 40.8 * 0.3048 Vmax = W * l / 2.0 Av = (2 / 3.0) * As fv = Vmax / Av ratio7 = abs(fv - Fv_adj) / Fv_adj # Bearing l = 45.6 * 0.3048 bArea = 2 * 0.0254 * joistSection.b fcT = W * l / 2.0 / bArea ratio8 = abs(fcT - FcT_adj) / FcT_adj # Deflection l = 13.8 * 0.3048 W = joistSpacing * (L) beam = sb.SimpleBeam(E, Is, l) delta = beam.getDeflectionUnderUniformLoad(W, beam.l / 2.0) deltaRef = l / 360.0 ratio9 = abs(delta - deltaRef) / deltaRef # Floor vibration l = 12.5 * 0.3048 W = joistSpacing * (40 * 47.88026) beam = sb.SimpleBeam(E, Is, l) delta = beam.getDeflectionUnderUniformLoad(W, beam.l / 2.0) deltaRef = l / 360.0 ratio10 = abs(delta - deltaRef) / deltaRef ''' print('Fb= ', Fb/1e6, ' MPa') print('Fb\'= ', Fb_adj/1e6, ' MPa') print('fb= ', fb/1e6, ' MPa')
# -*- coding: utf-8 -*- from rough_calculations import ng_simple_beam as sb from rough_calculations import ng_simple_bending_reinforcement from matplotlib import pyplot from mpl_toolkits.mplot3d import Axes3D beam = sb.SimpleBeam() beam.l = 2.5 #Loads Q = 1.5 * -135e3 # N d = 1.2 q = 1.35 * (-0.25 * 2500 * 9.81 - 2.35e3 - 8.1e3) # N/m MdMax = -1e15 aMdMax = 0.0 xMdMax = 0.0 VMdMax = 0.0 VdMax = -1e15 aVdMax = 0.0 xVdMax = 0.0 MVdMax = 0.0 sequence_containing_x_vals = list() sequence_containing_y_vals = list() sequence_containing_z_vals = list() for i in range(0, 60): a = beam.l / 40.0 * i - d for j in range(0, 40):
designLoad = 77.3 * pound2Newton / foot2meter # Design load designMoment = designLoad * span**2 / 8.0 wood = dimensional_lumber.SouthernPineWood(name='SouthernPine', grade='no_2', sub_grade='') joist = mat.DimensionLumberSection(name='2x10', woodMaterial=wood) S = joist.Wzel() # Section modulus I = joist.Iz() # Moment of inertia fb = designMoment / S # Bendign stress Fb = joist.getFb() E = joist.wood.E Fv = joist.wood.Fv beam = sb.SimpleBeam(E, I, span) delta = beam.getDeflectionUnderUniformLoad(designLoad, beam.l / 2.0) ratio1 = (Fb - 5.5e6) / 5.5e6 ratio2 = (Fv - 1.2e6) / 1.2e6 ratio3 = (delta - 20.8425973321e-3) / 20.8425973321e-3 ''' print('b= ', joist.b*1000, 'mm (',joist.b/inch2meter,'in)') print('h= ', joist.h*1000, 'mm (',joist.h/inch2meter,'in)') print('S= ', S, 'm3 (',S/inch2meter**3,'in3)') print('Bending stress fb= ', fb/1e6, 'MPa (',fb/pound2Newton*inch2meter**2,'psi)') print('Fb= ', Fb/1e6, 'MPa (',Fb/pound2Newton*inch2meter**2,'psi)') print('E= ', E/1e6, 'MPa (',E/pound2Newton*inch2meter**2,'psi)') print('delta= ', delta*1000, 'mm (',delta/inch2meter,'in)') print('Fv= ', Fv/1e6, 'MPa (',Fv/pound2Newton*inch2meter**2,'psi)') print('ratio1= ',ratio1)