import xc_base
import geom
import xc
from materials.sections import section_properties as sp
from materials.ec3 import EC3_materials
# Problem type
mainBeam = xc.FEProblem()
mainBeam.title = 'Falsework support beams'
preprocessor = mainBeam.getPreprocessor
#Materials
S355JR = EC3_materials.S355JR
S355JR.gammaM = 1.00
chsSection = EC3_materials.CHSShape(S355JR, 'CHS_48.3_5.0')
# ^ Y
# |
#
# o o o o
#
# -> Z
#
# o o o o
#
positions = [[-0.5, -0.5], [-0.5, -0.1], [-0.5, 0.1], [-0.5, 0.5], [0.5, -0.5],
[0.5, -0.1], [0.5, 0.1], [0.5, 0.5]]
sectionList = []
for p in positions:
alpha_HE300_comp = 2
beta_HE300_comp = 5 * Nd / 2643750.0
SHS175 = EC3_materials.SHSShape(steel=S235JR, name='SHS175x175x8')
alpha_SHS175, beta_SHS175 = SHS175.getBiaxBendCoeffs(Nd,
SHS175.getNcRd(sctClass))
alpha_SHS175_comp = 1.66 / (1 - 1.13 * (Nd / 1246205.0)**2)
beta_SHS175_comp = 1.66 / (1 - 1.13 * (Nd / 1246205.0)**2)
RHS250 = EC3_materials.RHSShape(steel=S235JR, name='RHS250x150x16')
alpha_RHS250, beta_RHS250 = RHS250.getBiaxBendCoeffs(Nd / 2.,
RHS250.getNcRd(sctClass))
alpha_RHS250_comp = 1.66 / (1 - 1.13 * (0.5 * Nd / 2702500.0)**2)
beta_RHS250_comp = 1.66 / (1 - 1.13 * (0.5 * Nd / 2702500.0)**2)
CHS193 = EC3_materials.CHSShape(steel=S235JR, name='CHS_193.7_11.0')
alpha_CHS193, beta_CHS193 = CHS193.getBiaxBendCoeffs(Nd,
CHS193.getNcRd(sctClass))
alpha_CHS193_comp = 2
beta_CHS193_comp = 2
ratioLst = [
alpha_IPE300 - alpha_IPE300_comp, beta_IPE300 - beta_IPE300_comp,
alpha_IPN300 - alpha_IPN300_comp, beta_IPN300 - beta_IPN300_comp,
alpha_HE300 - alpha_HE300_comp, beta_HE300 - beta_HE300_comp,
alpha_SHS175 - alpha_SHS175_comp, beta_SHS175 - beta_SHS175_comp,
alpha_RHS250 - alpha_RHS250_comp, beta_RHS250 - beta_RHS250_comp,
alpha_CHS193 - alpha_CHS193_comp, beta_CHS193 - beta_CHS193_comp
]
import os
