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