concr = EHE_materials.HA25 concr.alfacc = 0.85 # f_maxd= 0.85*fcd concrete long term compressive strength factor (normally alfacc=1) steel = EHE_materials.B500S # Section geometry section = def_simple_RC_section.RCRectangularSection(name='test', width=width, depth=depth, concrType=concr, reinfSteelType=steel) section.positvRebarRows = def_simple_RC_section.LongReinfLayers([lowerRow]) section.negatvRebarRows = def_simple_RC_section.LongReinfLayers([upperRow]) # section.shReinfY= shearReinf section.defRCSection(preprocessor, matDiagType='d') scc3d_testing_bench.sectionModel(preprocessor, section.sectionName) # Constraints modelSpace = predefined_spaces.getStructuralMechanics3DSpace(preprocessor) modelSpace.fixNode000_000(1) # Loads definition lp0 = modelSpace.newLoadPattern(name='0') lp0.newNodalLoad(2, xc.Vector([NDato, 0, 0, 0, MyDato, MzDato])) lp1 = modelSpace.newLoadPattern(name='1') lp1.newNodalLoad(2, xc.Vector([NDato, 0, 0, 0, MzDato / 10.0, MyDato / 10.0])) lp2 = modelSpace.newLoadPattern(name='2') lp2.newNodalLoad(2, xc.Vector([NDato, 0, 0, 0, 0, 0])) # We add the load case to domain. modelSpace.addLoadCaseToDomain(lp0.name)
dgDB500S = EHE_materials.B500S.getDiagD(preprocessor) import os pth = os.path.dirname(__file__) if (not pth): pth = "." #print "pth= ", pth execfile(pth + "/geomSeccBarras.py") secBarras = preprocessor.getMaterialHandler.newMaterial( "fiber_section_3d", "secBarras") fiberSectionRepr = secBarras.getFiberSectionRepr() fiberSectionRepr.setGeomNamed("geomSecBarras") secBarras.setupFibers() scc3d_testing_bench.sectionModel(preprocessor, "secBarras") # Constraints modelSpace = predefined_spaces.getStructuralMechanics3DSpace(preprocessor) modelSpace.fixNode000_000(1) modelSpace.fixNodeF00_00F(2) # Loads definition loadHandler = preprocessor.getLoadHandler lPatterns = loadHandler.getLoadPatterns #Load modulation. ts = lPatterns.newTimeSeries("constant_ts", "ts") lPatterns.currentTimeSeries = "ts" #Load case definition
#SHS50x50x2_5.A= 4.51e-4 #SHS50x50x2_5.Iz= 17.9e-8 #SHS50x50x2_5.Iy= 17.9e-8 SHS50x50x2_5.Wzel = 7.16e-6 SHS50x50x2_5.Wzpl = 8.47e-6 SHS50x50x2_5.iz = 1.99e-2 SHS50x50x2_5.iy = 1.99e-2 #SHS50x50x2_5.J= 26.8e-8 SHS50x50x2_5.E = 210000e6 SHS50x50x2_5.nu = 0.3e6 #SHS50x50x2_5.G= SHS50x50x2_5.E/(2*(1+SHS50x50x2_5.nu)) SHS50x50x2_5.alpha = 5 / 6 # XXX Corregir # Materials definition mat = SHS50x50x2_5.defElasticSection3d(preprocessor, matSHS50x50x2_5) elemZLS = scc3d_testing_bench.sectionModel(preprocessor, SHS50x50x2_5.sectionName) # Constraints modelSpace = predefined_spaces.getStructuralMechanics3DSpace(preprocessor) modelSpace.fixNode000_000(1) modelSpace.fixNodeF00_00F(2) # Loads definition loadHandler = preprocessor.getLoadHandler lPatterns = loadHandler.getLoadPatterns #Load modulation. ts = lPatterns.newTimeSeries("constant_ts", "ts") lPatterns.currentTimeSeries = "ts" #Load case definition
#datosScc1LosC.negatvRebars.setUp(0,0.0,0.0,datosScc1LosC.b,0.1) #datosScc1LosC.positvRebars.setUp(6,26e-3,(areaFi22+areaFi26)/2.0,datosScc1LosC.b,0.05+0.026/2.0) NDato = 0 # Axial force for crack control checking. MyDato = 117e3 # Bending moment force for crack control checking. feProblem = xc.FEProblem() preprocessor = feProblem.getPreprocessor # Materials definition # concreteMatTag= EHE_materials.HA30.defDiagK(preprocessor) # reinfSteelMaterialTag= EHE_materials.B500S.defDiagK(preprocessor) datosScc1LosC.defRCSimpleSection(preprocessor, "k") scc3d_testing_bench.sectionModel(preprocessor, datosScc1LosC.sectionName) # Constraints modelSpace = predefined_spaces.getStructuralMechanics3DSpace(preprocessor) modelSpace.fixNode000_000(1) modelSpace.fixNodeF00_0F0(2) # Loads definition cargas = preprocessor.getLoadLoader casos = cargas.getLoadPatterns #Load modulation. ts = casos.newTimeSeries("constant_ts", "ts") casos.currentTimeSeries = "ts" #Load case definition
EHE_materials.Y1860S7.tInic()) tag = EHE_materials.HP45.defDiagD(preprocessor) import os pth = os.path.dirname(__file__) if (not pth): pth = "." #print "pth= ", pth execfile(pth + "/prestressed_concrete_section_01.py") materialHandler = preprocessor.getMaterialHandler secHP = materialHandler.newMaterial("fiber_section_3d", "secHP") fiberSectionRepr = secHP.getFiberSectionRepr() fiberSectionRepr.setGeomNamed("prestressedConcretSectionGeom01") secHP.setupFibers() elem = scc3d_testing_bench.sectionModel(preprocessor, "secHP") # Constraints modelSpace = predefined_spaces.getStructuralMechanics3DSpace(preprocessor) modelSpace.fixNode000_000(1) modelSpace.fixNodeF00_00F(2) # Loads definition loadHandler = preprocessor.getLoadHandler lPatterns = loadHandler.getLoadPatterns #Load modulation. ts = lPatterns.newTimeSeries("constant_ts", "ts") lPatterns.currentTimeSeries = "ts" #Load case definition
respVz= typical_materials.defElasticMaterial(preprocessor, "respVz",1e6) # Shear response in y direction. # Section geometry #creation geomRectang= preprocessor.getMaterialHandler.newSectionGeometry("geomRectang") reg= scc1x1.getRegion(geomRectang,"elast") sa= preprocessor.getMaterialHandler.newMaterial("fiberSectionShear3d","sa") fiberSectionRepr= sa.getFiberSectionRepr() fiberSectionRepr.setGeomNamed("geomRectang") sa.setupFibers() extractFiberSectionProperties(sa,scc1x1) sa.setRespVyByName("respVy") sa.setRespVzByName("respVz") sa.setRespTByName("respT") scc3d_testing_bench.sectionModel(preprocessor, "sa") # Constraints modelSpace= predefined_spaces.getStructuralMechanics3DSpace(preprocessor) modelSpace.fixNode000_000(1) # Loads definition loadHandler= preprocessor.getLoadHandler lPatterns= loadHandler.getLoadPatterns #Load modulation. ts= lPatterns.newTimeSeries("constant_ts","ts") lPatterns.currentTimeSeries= "ts" #Load case definition lp0= lPatterns.newLoadPattern("default","0") #lPatterns.currentLoadPattern= "0"
__email__= "*****@*****.**" feProblem= xc.FEProblem() feProblem.logFileName= "/tmp/erase.log" # Don't pring warnings preprocessor= feProblem.getPreprocessor # Rectangular cross-section definition scc10x20= section_properties.RectangularSection(name="rectang",b=.10,h=.20) # Section geometry. matscc10x20=typical_materials.MaterialData(name='mtrectang',E=2.1e6,nu=0.3,rho=2500) # Section material. feProblem= xc.FEProblem() preprocessor= feProblem.getPreprocessor # Materials definition matPoteau= scc10x20.defElasticShearSection3d(preprocessor,matscc10x20) elemZLS= scc3d_testing_bench.sectionModel(preprocessor, scc10x20.sectionName) # Constraints modelSpace= predefined_spaces.getStructuralMechanics3DSpace(preprocessor) modelSpace.fixNode000_000(1) # Load definition. lp0= modelSpace.newLoadPattern(name= '0') loadMz= 1e3 lp0.newNodalLoad(2,xc.Vector([0,0,0,0,0,loadMz])) # We add the load case to domain. modelSpace.addLoadCaseToDomain(lp0.name)
__version__= "3.0" __email__= "*****@*****.**" feProblem= xc.FEProblem() feProblem.logFileName= "/tmp/erase.log" # Ignore warning messages preprocessor= feProblem.getPreprocessor # Section geometry scc10x20= section_properties.RectangularSection(name="rectang",b=.10,h=.20) # Section geometry. # Section material matscc10x20=typical_materials.MaterialData(name='mtrectang',E=2.1e6,nu=0.3,rho=2500) # Section material. # Materials definition matPoteau= scc10x20.defElasticSection3d(preprocessor,matscc10x20) elemZLS= scc3d_testing_bench.sectionModel(preprocessor, scc10x20.sectionName) # Constraints modelSpace= predefined_spaces.getStructuralMechanics3DSpace(preprocessor) modelSpace.fixNode000_000(1) modelSpace.fixNodeF00_00F(2) # Loads definition loadHandler= preprocessor.getLoadHandler lPatterns= loadHandler.getLoadPatterns #Load modulation. ts= lPatterns.newTimeSeries("constant_ts","ts") lPatterns.currentTimeSeries= "ts" #Load case definition lp0= lPatterns.newLoadPattern("default","0") #lPatterns.currentLoadPattern= "0"
# Materials definition epp= typical_materials.defElasticPPMaterial(preprocessor, "epp",E,fy,-fy) # Section geometry # creation geomRectang= preprocessor.getMaterialHandler.newSectionGeometry("geomRectang") #generation of a quadrilateral region of the scc10x20 sizes and number of #divisions made of material nmbMat reg= scc10x20.getRegion(gm=geomRectang,nmbMat="epp") rectang= preprocessor.getMaterialHandler.newMaterial("fiber_section_3d","rectang") fiberSectionRepr= rectang.getFiberSectionRepr() fiberSectionRepr.setGeomNamed("geomRectang") rectang.setupFibers() extractFiberSectionProperties(rectang,scc10x20) scc3d_testing_bench.sectionModel(preprocessor, "rectang") # Constraints modelSpace= predefined_spaces.getStructuralMechanics3DSpace(preprocessor) modelSpace.fixNode000_000(1) modelSpace.fixNodeF00_00F(2) # Loads definition loadHandler= preprocessor.getLoadHandler lPatterns= loadHandler.getLoadPatterns #Load modulation. ts= lPatterns.newTimeSeries("constant_ts","ts") lPatterns.currentTimeSeries= "ts" #Load case definition lp0= lPatterns.newLoadPattern("default","0")
respT = typical_materials.defElasticMaterial(preprocessor, "respT", 1e10) # Torsion response. respVy = typical_materials.defElasticMaterial( preprocessor, "respVy", 1e6) # Shear response in y direction. respVz = typical_materials.defElasticMaterial( preprocessor, "respVz", 1e3) # Shear response in y direction. materialHandler = preprocessor.getMaterialHandler sctFibers = materialHandler.newMaterial("fiberSectionShear3d", "sctFibers") fiberSectionRepr = sctFibers.getFiberSectionRepr() fiberSectionRepr.setGeomNamed("geomSectFibers") sctFibers.setupFibers() sctFibers.setRespVyByName("respVy") sctFibers.setRespVzByName("respVz") sctFibers.setRespTByName("respT") scc3d_testing_bench.sectionModel(preprocessor, "sctFibers") # #report of the section material # sectParam=section_report.SectionInfo(preprocessor=preprocessor,sectName='example_7.3_EC2W',sectDescr='Test example 7.3 EC2 Worked examples. Section definition',concrete=concrete,rfSteel=rfSteel,concrDiag=concrDiagram,rfStDiag=steelDiagram,geomSection=geomSectFibers,width=width,depth=depth) #Obtains section parameters for report # Elements definition ''' elements= preprocessor.getElementHandler elements.defaultMaterial='sctFibers' elements.dimElem= 3 # Dimension of element space elements.defaultTag= 1 elem= elements.newElement("ZeroLengthSection",xc.ID([1,2])) ''' # Constraints modelSpace = predefined_spaces.getStructuralMechanics3DSpace( preprocessor
feProblem= xc.FEProblem() preprocessor= feProblem.getPreprocessor # Materials definition tag= EHE_materials.B500S.defDiagD(preprocessor) import os pth= os.path.dirname(__file__) if(not pth): pth= "." # print("pth= ", pth) exec(open(pth+"/barsSectionGeometry.py").read()) barsSection= preprocessor.getMaterialHandler.newMaterial("fiber_section_GJ","barsSection") fiberSectionRepr= barsSection.getFiberSectionRepr() fiberSectionRepr.setGeomNamed("barsSectionGeometry") barsSection.setupFibers() scc3d_testing_bench.sectionModel(preprocessor, "barsSection") # Constraints modelSpace= predefined_spaces.getStructuralMechanics3DSpace(preprocessor) modelSpace.fixNode000_000(1) modelSpace.fixNodeF00_00F(2) # Loads definition lp0= modelSpace.newLoadPattern(name= '0') lp0.newNodalLoad(2,xc.Vector([0,0,0,0,0,MzDato])) # We add the load case to domain. modelSpace.addLoadCaseToDomain(lp0.name) # Solution procedure
NDato= 0 # Axial force for crack control checking. MyDato= 1000 # Bending moment force for crack control checking. feProblem= xc.FEProblem() preprocessor= feProblem.getPreprocessor # Materials definition concreteMatTag= EHE_materials.HA25.defDiagK(preprocessor) reinfSteelMaterialTag= EHE_materials.B500S.defDiagK(preprocessor) datosScc1LosC.defRCSimpleSection(preprocessor, "k") scc3d_testing_bench.sectionModel(preprocessor, datosScc1LosC.sectionName) # Constraints modelSpace= predefined_spaces.getStructuralMechanics3DSpace(preprocessor) modelSpace.fixNode000_000(1) modelSpace.fixNodeF00_0F0(2) # Loads definition loadHandler= preprocessor.getLoadHandler lPatterns= loadHandler.getLoadPatterns #Load modulation. ts= lPatterns.newTimeSeries("constant_ts","ts") lPatterns.currentTimeSeries= "ts" #Load case definition
# Materials definition tag= EHE_materials.Y1860S7.defDiagD(preprocessor, EHE_materials.Y1860S7.tInic()) tag= EHE_materials.HP45.defDiagD(preprocessor) import os pth= os.path.dirname(__file__) if(not pth): pth= "." #print "pth= ", pth execfile(pth+"/prestressed_concrete_section_01.py") materiales= preprocessor.getMaterialHandler secHP= materiales.newMaterial("fiber_section_3d","secHP") fiberSectionRepr= secHP.getFiberSectionRepr() fiberSectionRepr.setGeomNamed("prestressedConcretSectionGeom01") secHP.setupFibers() elem= scc3d_testing_bench.sectionModel(preprocessor, "secHP") # Constraints modelSpace= predefined_spaces.getStructuralMechanics3DSpace(preprocessor) modelSpace.fixNode000_000(1) modelSpace.fixNodeF00_00F(2) # Loads definition loadHandler= preprocessor.getLoadHandler lPatterns= loadHandler.getLoadPatterns #Load modulation. ts= lPatterns.newTimeSeries("constant_ts","ts") lPatterns.currentTimeSeries= "ts" #Load case definition
#SHS50x50x2_5.Iz= 17.9e-8 #SHS50x50x2_5.Iy= 17.9e-8 SHS50x50x2_5.Wzel= 7.16e-6 SHS50x50x2_5.Wzpl= 8.47e-6 SHS50x50x2_5.iz= 1.99e-2 SHS50x50x2_5.iy= 1.99e-2 #SHS50x50x2_5.J= 26.8e-8 SHS50x50x2_5.E= 210000e6 SHS50x50x2_5.nu= 0.3e6 #SHS50x50x2_5.G= SHS50x50x2_5.E/(2*(1+SHS50x50x2_5.nu)) SHS50x50x2_5.alpha= 5/6 # XXX Corregir # Materials definition mat= SHS50x50x2_5.defElasticSection3d(preprocessor,matSHS50x50x2_5) elemZLS= scc3d_testing_bench.sectionModel(preprocessor, SHS50x50x2_5.sectionName) # Constraints modelSpace= predefined_spaces.getStructuralMechanics3DSpace(preprocessor) modelSpace.fixNode000_000(1) modelSpace.fixNodeF00_00F(2) # Loads definition loadHandler= preprocessor.getLoadHandler lPatterns= loadHandler.getLoadPatterns #Load modulation. ts= lPatterns.newTimeSeries("constant_ts","ts") lPatterns.currentTimeSeries= "ts" #Load case definition