__license__ = "GPL"
__version__ = "3.0"
__email__ = "*****@*****.**"
inch2meter = 0.0254
pound2Newton = 4.4482216282509
foot2meter = 0.3048
span = 16 * foot2meter
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
'''
l.setProp('webOrientation', -kVector)
l.setElemSize(elementLength)
## Template coordinate transformation.
linCT = preprocessor.getTransfCooHandler.newLinearCrdTransf3d('linCT')
linCT.xzVector = xc.Vector([1.0, 0, 0])
# Materials
shortBeamsSection = structural_panels.LVLHeaderSections['1.75x11-7/8']
beamsSection = structural_panels.LVLHeaderSections['5.25x11-7/8']
longBeamsSection = structural_panels.LVLHeaderSections['7x11-7/8'] #'3.5x18']
longJoistsSection = pr_400_i_joists.pr400_i_joists['PRI-20_241']
joistsWood = dl.SouthernPineWood(name='SouthernPine',
grade='no_2',
sub_grade='')
joistsSection = AWCNDS_materials.DimensionLumberSection(
name='2x8', woodMaterial=joistsWood)
# Mesh generation.
seedElemHandler = preprocessor.getElementHandler.seedElemHandler
seedElemHandler.defaultTransformation = "linCT"
def createMesh(xcSet, section):
xcSection = section.defElasticShearSection3d(preprocessor)
seedElemHandler.defaultMaterial = section.name
for l in xcSet.getLines:
vDir = l.getProp('webOrientation')
l.setProp('crossSection', section)
linCT.xzVector = xc.Vector([vDir.x, vDir.y, vDir.z])
elem = seedElemHandler.newElement("ElasticBeam3d", xc.ID([0, 0]))