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)
