## Longitudinal reinforcement
rebarDiam = 16e-3
rebarArea = EHE_materials.Fi16
numOfRebars = int(math.ceil(math.pi * (pierDiameter - 2 * nomCover) / .05))
mainReinf = def_simple_RC_section.LongReinfLayers([
def_simple_RC_section.ReinfRow(
rebarsDiam=rebarDiam,
nRebars=numOfRebars,
width=math.pi * (pierDiameter - 2 * (nomCover + shearReinfDiam)),
nominalCover=nomCover)
])
# Section geometry
section = def_column_RC_section.RCCircularSection(
name='test',
sectionDescr='circular section',
Rext=pierDiameter / 2.0,
concrType=concrete,
reinfSteelType=reinfSteel)
section.mainReinf = mainReinf
section.shReinf = shearReinf
# All this is made only to allow the
# creation of the RCMaterialDistribution
feProblem = xc.FEProblem()
preprocessor = feProblem.getPreprocessor
nodeHandler = preprocessor.getNodeHandler
nodeHandler.dimSpace = 1 # One coordinate for each node.
nodeHandler.numDOFs = 1 # One degree of freedom for each node.
n1 = nodeHandler.newNodeX(1)
n2 = nodeHandler.newNodeX(2)
elementHandler = preprocessor.getElementHandler
__author__ = "Luis C. Pérez Tato (LCPT) and Ana Ortega (A_OO)" __copyright__ = "Copyright 2015, LCPT and AO_O" __license__ = "GPL" __version__ = "3.0" __email__ = "[email protected] [email protected]" # Materials definition concr = EHE_materials.HA40 steel = EHE_materials.B500S # Section geometry diameter = 350e-3 # Cross-section diameter [m] cover = 0.025 # Cover [m] section = def_column_RC_section.RCCircularSection(name='test', Rext=diameter / 2.0, concrType=concr, reinfSteelType=steel) section.nDivIJ = 10 # Longitudinal reinforcement rebarDiam = 16e-3 rebarArea = EHE_materials.Fi16 numOfRebars = 12 # # Shear reinforcement # shearReinfArea= EHE_materials.Fi6 # shearReinfDiam= 6e-3 # numRamas= 2 # shearReinf= def_simple_RC_section.ShearReinforcement(familyName= "sh",nShReinfBranches= numRamas, areaShReinfBranch= shearReinfArea, shReinfSpacing= 0.2, angAlphaShReinf= math.pi/2.0,angThetaConcrStruts= math.pi/4.0) nCover = cover
radius = 0.75 / 2.0 # Cross-section radius expressed in meters.
rebarDiam = 20e-3 # Bar diameter expressed in meters.
nCover = 0.06 + 12e-3
cover = nCover + rebarDiam / 2.0 # Concrete cover expressed in meters.
rebarArea = math.pi * (rebarDiam /
2.0)**2 # Rebar area expressed in square meters.
feProblem = xc.FEProblem()
preprocessor = feProblem.getPreprocessor
# Materials definition
concr = EHE_materials.HA30
concr.alfacc = 0.85 #f_maxd= 0.85*fcd concrete long term compressive strength factor (normally alfacc=1)
steel = EHE_materials.B500S
section = def_column_RC_section.RCCircularSection(name='test',
Rext=radius,
concrType=concr,
reinfSteelType=steel)
section.mainReinf = def_simple_RC_section.LongReinfLayers([
def_simple_RC_section.ReinfRow(rebarsDiam=rebarDiam,
nRebars=14,
width=2 * math.pi * (radius - cover),
nominalCover=nCover)
])
section.defRCSection(preprocessor, matDiagType='d')
diagIntsecHA = section.defInteractionDiagram(preprocessor)
# Compute capacity factors.
fc1 = diagIntsecHA.getCapacityFactor(geom.Pos3d(1850e3, 0, 0))