def getDiagIntSection1(mdlr):
ancho = 1
canto = 0.3
recpos = 0.076 + 0.008
recneg = 0.076 + 0.008
geomSecHA2 = mdlr.getMaterialHandler.newSectionGeometry('geomSecHA2')
regions = geomSecHA2.getRegions
concrete= regions.newQuadRegion(EHE_materials.HA30.nmbDiagD)
concrete.nDivIJ = 10
concrete.nDivJK = 10
concrete.pMin = geom.Pos2d(-ancho / 2, -canto / 2)
concrete.pMax = geom.Pos2d(ancho / 2, canto / 2)
reinforcement = geomSecHA2.getReinfLayers
bottomReinforcement = reinforcement.newStraightReinfLayer(EHE_materials.B500S.nmbDiagD)
bottomReinforcement.numReinfBars = 7
bottomReinforcement.barArea = areaFi20
bottomReinforcement.p1 = geom.Pos2d(-ancho / 2 + recneg, -canto / 2 + recneg)
bottomReinforcement.p2 = geom.Pos2d(ancho / 2 - recneg, -canto / 2 + recneg)
topReinforcement = reinforcement.newStraightReinfLayer(EHE_materials.B500S.nmbDiagD)
topReinforcement.numReinfBars = 7
topReinforcement.barArea = areaFi20
topReinforcement.p1 = geom.Pos2d(-ancho / 2 + recpos, canto / 2 - recpos)
topReinforcement.p2 = geom.Pos2d(ancho / 2 - recpos, canto / 2 - recpos)
materialHandler = mdlr.getMaterialHandler
secHA2 = materialHandler.newMaterial('fiber_section_3d', 'secHA2')
fiberSectionRepr = secHA2.getFiberSectionRepr()
fiberSectionRepr.setGeomNamed('geomSecHA2')
secHA2.setupFibers()
param = xc.InteractionDiagramParameters()
param.concreteTag = EHE_materials.HA30.tagDiagD
param.reinforcementTag = EHE_materials.B500S.tagDiagD
diagIntSecHA2 = materialHandler.calcInteractionDiagram('secHA2', param)
return diagIntSecHA2
def getDiagIntSection1(mdlr):
ancho = 1
canto = 0.3
recpos = 0.076 + 0.008
recneg = 0.076 + 0.008
geomSecHA2 = mdlr.getMaterialHandler.newSectionGeometry('geomSecHA2')
regiones = geomSecHA2.getRegions
concrete = regiones.newQuadRegion(EHE_materials.HA30.nmbDiagD)
concrete.nDivIJ = 10
concrete.nDivJK = 10
concrete.pMin = geom.Pos2d(-ancho / 2, -canto / 2)
concrete.pMax = geom.Pos2d(ancho / 2, canto / 2)
armaduras = geomSecHA2.getReinfLayers
armaduraInf = armaduras.newStraightReinfLayer(EHE_materials.B500S.nmbDiagD)
armaduraInf.numReinfBars = 7
armaduraInf.barArea = areaFi20
armaduraInf.p1 = geom.Pos2d(-ancho / 2 + recneg, -canto / 2 + recneg)
armaduraInf.p2 = geom.Pos2d(ancho / 2 - recneg, -canto / 2 + recneg)
armaduraSup = armaduras.newStraightReinfLayer(EHE_materials.B500S.nmbDiagD)
armaduraSup.numReinfBars = 7
armaduraSup.barArea = areaFi20
armaduraSup.p1 = geom.Pos2d(-ancho / 2 + recpos, canto / 2 - recpos)
armaduraSup.p2 = geom.Pos2d(ancho / 2 - recpos, canto / 2 - recpos)
materiales = mdlr.getMaterialHandler
secHA2 = materiales.newMaterial('fiber_section_3d', 'secHA2')
fiberSectionRepr = secHA2.getFiberSectionRepr()
fiberSectionRepr.setGeomNamed('geomSecHA2')
secHA2.setupFibers()
param = xc.InteractionDiagramParameters()
param.concreteTag = EHE_materials.HA30.tagDiagD
param.tagArmadura = EHE_materials.B500S.tagDiagD
diagIntSecHA2 = materiales.calcInteractionDiagram('secHA2', param)
return diagIntSecHA2
def getDiagIntSection2(mdlr):
ancho = 1
canto = 0.3
recpos = 0.06 + 0.008
recneg = 0.06 + 0.008
geomSecHA1 = mdlr.getMaterialLoader.newSectionGeometry('geomSecHA1')
regiones = geomSecHA1.getRegions
hormigon = regiones.newQuadRegion(hormigonesEHE.HA30.nmbDiagD)
hormigon.nDivIJ = 10
hormigon.nDivJK = 10
hormigon.pMin = geom.Pos2d(-ancho / 2, -canto / 2)
hormigon.pMax = geom.Pos2d(ancho / 2, canto / 2)
armaduras = geomSecHA1.getReinfLayers
armaduraInf = armaduras.newStraightReinfLayer(acerosEHE.B500S.nmbDiagD)
armaduraInf.numReinfBars = 5
armaduraInf.barArea = areaFi12
armaduraInf.p1 = geom.Pos2d(-ancho / 2 + recneg, -canto / 2 + recneg)
armaduraInf.p2 = geom.Pos2d(ancho / 2 - recneg, -canto / 2 + recneg)
armaduraSup = armaduras.newStraightReinfLayer(acerosEHE.B500S.nmbDiagD)
armaduraSup.numReinfBars = 5
armaduraSup.barArea = areaFi12
armaduraSup.p1 = geom.Pos2d(-ancho / 2 + recpos, canto / 2 - recpos)
armaduraSup.p2 = geom.Pos2d(ancho / 2 - recpos, canto / 2 - recpos)
materiales = mdlr.getMaterialLoader
secHA1 = materiales.newMaterial('fiber_section_3d', 'secHA1')
fiberSectionRepr = secHA1.getFiberSectionRepr()
fiberSectionRepr.setGeomNamed('geomSecHA1')
secHA1.setupFibers()
param = xc.InteractionDiagramParameters()
param.tagHormigon = hormigonesEHE.HA30.tagDiagD
param.tagArmadura = acerosEHE.B500S.tagDiagD
diagIntSecHA1 = materiales.calcInteractionDiagram('secHA1', param)
return diagIntSecHA1
def defInteractionDiagramParameters(self, preprocessor):
''' parameters for interaction diagrams. '''
self.param = xc.InteractionDiagramParameters()
if (self.diagType == "d"):
self.param.concreteTag = self.concrType.matTagD
self.param.reinforcementTag = self.reinfSteelType.matTagD
elif (self.diagType == "k"):
self.param.concreteTag = self.concrType.matTagK
self.param.reinforcementTag = self.reinfSteelType.matTagK
return self.param
def defInteractionDiagramParameters(self, preprocessor):
''' parameters for interaction diagrams. '''
self.param= xc.InteractionDiagramParameters()
if(self.diagType=="d"):
self.param.tagHormigon= self.concrType.matTagD
self.param.tagArmadura= self.reinfSteelType.matTagD
elif(self.diagType=="k"):
self.param.tagHormigon= self.concrType.matTagK
self.param.tagArmadura= self.reinfSteelType.matTagK
return self.param
reinforcement = pileGeometry.getReinfLayers
reinforcement = reinforcement.newCircReinfLayer(EHE_materials.B500S.nmbDiagD)
reinforcement.numReinfBars = 14
reinforcement.barArea = areaFi20
reinforcement.initAngle = 0.0
reinforcement.finalAngle = 2 * math.pi
reinforcement.radius = concrete.extRad - cover
materiales = preprocessor.getMaterialLoader
secHA = materiales.newMaterial("fiber_section_3d", "secHA")
fiberSectionRepr = secHA.getFiberSectionRepr()
fiberSectionRepr.setGeomNamed("pileGeometry")
secHA.setupFibers()
fibras = secHA.getFibers()
param = xc.InteractionDiagramParameters()
param.concreteTag = EHE_materials.HA30.matTagD
param.tagArmadura = EHE_materials.B500S.matTagD
diagIntsecHA = materiales.calcInteractionDiagram("secHA", param)
fc1 = diagIntsecHA.getCapacityFactor(geom.Pos3d(1850e3, 0, 0))
fc2 = diagIntsecHA.getCapacityFactor(geom.Pos3d(-152e3, 530e3, 0))
fc3 = diagIntsecHA.getCapacityFactor(geom.Pos3d(-152e3, 0, 590e3))
fc4 = diagIntsecHA.getCapacityFactor(geom.Pos3d(-9097e3, 0.0, 0.0))
ratio1 = abs(fc1 - 1)
ratio2 = abs(fc2 - 1)
ratio3 = abs(fc3 - 1.0)
ratio4 = abs(fc4 - 1.0)
'''
print "fc1= ",fc1
