# #report of the section material # sectParam=sectionReport.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 # sectParam.writeReport(archTex='figures/sections/secEx73.tex', pathFigura='figures/sections/secEx73.eps') # Elements definition elementos = preprocessor.getElementLoader elementos.defaultMaterial = 'sctFibers' elementos.dimElem = 1 elementos.defaultTag = 1 elem = elementos.newElement("zero_length_section", xc.ID([1, 2])) # Constraints constCont = preprocessor.getConstraintLoader #constraints container fix_node_6dof.fixNode6DOF(constCont, 1) fix_node_6dof.Nodo6DOFMovXGiroYLibres(constCont, 2) # Loads definition cargas = preprocessor.getLoadLoader #loads container casos = cargas.getLoadPatterns #Load modulation. ts = casos.newTimeSeries("constant_ts", "ts") casos.currentTimeSeries = "ts" #Load case definition lp0 = casos.newLoadPattern("default", "0") pointLoad = xc.Vector([0, 0, 0, 0, M_y, 0]) lp0.newNodalLoad(2, pointLoad) #applies the point load on node 2 #We add the load case to domain. casos.addToDomain("0") #reads load pattern "0" and adds it to the domain
reinforcementSup.p1= geom.Pos2d(cover-width/2.0,depth/2.0-cover) # Armadura superior. reinforcementSup.p2= geom.Pos2d(width/2.0-cover,depth/2.0-cover) materiales= preprocessor.getMaterialLoader secHA= materiales.newMaterial("fiber_section_3d","secHA") fiberSectionRepr= secHA.getFiberSectionRepr() fiberSectionRepr.setGeomNamed("geomSecHA") secHA.setupFibers() banco_pruebas_scc3d.modeloSecc3d(preprocessor, "secHA") # Constraints coacciones= preprocessor.getConstraintLoader fix_node_6dof.fixNode6DOF(coacciones,1) fix_node_6dof.Nodo6DOFMovXGiroYLibres(coacciones,2) # Loads definition cargas= preprocessor.getLoadLoader casos= cargas.getLoadPatterns #Load modulation. ts= casos.newTimeSeries("constant_ts","ts") casos.currentTimeSeries= "ts" #Load case definition lp0= casos.newLoadPattern("default","0") lp0.newNodalLoad(2,xc.Vector([NDato,0,0,0,MyDato,0])) #We add the load case to domain. casos.addToDomain("0")