Example #1
0
 def defDiagK(self,preprocessor):
   ''' Returns XC uniaxial material (characteristic values). '''
   self.materialDiagramK= typical_materials.defSteel01(preprocessor,self.nmbDiagK,self.Es,self.fyk,self.bsh())
   self.matTagK= self.materialDiagramK.tag
   return self.materialDiagramK #30160925 was 'return self.matTagK'
Example #2
0
# Problem type
feProblem= xc.FEProblem()
preprocessor=  feProblem.getPreprocessor
nodes= preprocessor.getNodeLoader
modelSpace= predefined_spaces.StructuralMechanics2D(nodes)

nodes.defaultTag= 1 #First node number.
nod= nodes.newNodeXY(0,0.0)
nod= nodes.newNodeXY(L,0.0)

# Geometric transformation(s)
lin= modelSpace.newLinearCrdTransf("lin")

# Materials definition
fy= 275e6 # Yield stress of the steel.
steel= typical_materials.defSteel01(preprocessor, "steel",E,fy,0.001)

respT= typical_materials.defElasticMaterial(preprocessor, "respT",G*J) # Torsion response.
respVy= typical_materials.defElasticMaterial(preprocessor, "respVy",1e9) # Shear response in y direction.
respVz= typical_materials.defElasticMaterial(preprocessor, "respVz",1e9) # Shear response in z direction.
# Secciones
import os
pth= os.path.dirname(__file__)
#print "pth= ", pth
if(not pth):
  pth= "."
execfile(pth+"/../aux/testQuadRegion.py")

materiales= preprocessor.getMaterialLoader
quadFibers= materiales.newMaterial("fiber_section_3d","quadFibers")
fiberSectionRepr= quadFibers.getFiberSectionRepr()
Example #3
0
nDivJK= 32 # number of cells in JK direction
nFibTeor= nDivIJ*nDivJK
areaTeor= width*depth
iyTeor= 1/12.0*width*depth**3
izTeor= 1/12.0*depth*width**3
y0= 5
z0= 3.2

feProblem= xc.FEProblem()
preprocessor=  feProblem.getPreprocessor

#  Materials definition: uniaxial bilinear steel
fy= 50.0   # yield strength
E= 30000.0 # initial elastic tangent
b=0.001    # strain-hardening ratio: ratio between post-yield tangent and initial elastic tangent
steel= typical_materials.defSteel01(preprocessor=preprocessor,name="steel",E=E,fy=fy,b=b)

# Section geometry
# creation
geomRectang= preprocessor.getMaterialHandler.newSectionGeometry("geomRectang")
y1= width/2.0
z1= depth/2.0
#filling with regions
regions= geomRectang.getRegions
#generation of a quadrilateral region of the specified sizes and number of
#divisions for the cells (fibers) generation
steelRegion= regions.newQuadRegion("steel")   
steelRegion.nDivIJ= nDivIJ
steelRegion.nDivJK= nDivJK
steelRegion.pMin= geom.Pos2d(y0-y1,z0-z1)
steelRegion.pMax= geom.Pos2d(y0+y1,z0+z1)
Example #4
0
 def defDiagD(self,preprocessor):
   ''' Returns XC uniaxial material (design values). '''
   self.materialDiagramD= typical_materials.defSteel01(preprocessor,self.nmbDiagD,self.Es,self.fyd(),self.bsh())
   self.matTagD= self.materialDiagramD.tag
   return self.materialDiagramD #30160925 was 'return self.matTagD'
Example #5
0
preprocessor = prueba.getPreprocessor
nodes = preprocessor.getNodeLoader
# Problem type
modelSpace = predefined_spaces.StructuralMechanics3D(nodes)
nodes.defaultTag = 1  #First node number.
nod = nodes.newNodeXYZ(0, 0.0, 0.0)
nod = nodes.newNodeXYZ(L, 0.0, 0.0)

trfs = preprocessor.getTransfCooLoader
lin = trfs.newLinearCrdTransf3d("lin")
lin.xzVector = xc.Vector([0, 1, 0])

# Materials definition
fy = 275e6  # Tensión de cedencia of the steel
E = 210e9  # Young modulus of the steel.
acero = typical_materials.defSteel01(preprocessor, "acero", E, fy, 0.001)

import os
pth = os.path.dirname(__file__)
#print "pth= ", pth
if (not pth):
    pth = "."
execfile(pth + "/geomCuadFibrasTN.py"
         )  #Definition of section geometry (regions and rebars)
# Definition of a new empty fiber section named 'cuadFibrasTN' and stored in a
# Python variable of the same name (surprisingly enough).
cuadFibrasTN = preprocessor.getMaterialLoader.newMaterial(
    "fiber_section_3d", "cuadFibrasTN")
fiberSectionRepr = cuadFibrasTN.getFiberSectionRepr(
)  #Fiber section representation
# of 'cuadFibrasTN'
Example #6
0
feProblem= xc.FEProblem()
preprocessor=  feProblem.getPreprocessor   
nodes= preprocessor.getNodeHandler
# Problem type
modelSpace= predefined_spaces.StructuralMechanics3D(nodes)
nodes.defaultTag= 1 #First node number.
nod= nodes.newNodeXYZ(0,0.0,0.0)
nod= nodes.newNodeXYZ(L,0.0,0.0)


lin= modelSpace.newLinearCrdTransf("lin",xc.Vector([0,1,0]))

# Materials definition
fy= 275e6 # Yield stress of the steel.
E= 210e9 # Young modulus of the steel.
steel= typical_materials.defSteel01(preprocessor, "steel",E,fy,0.001)

respT= typical_materials.defElasticMaterial(preprocessor, "respT",1e10) # Torsion response.
respVy= typical_materials.defElasticMaterial(preprocessor, "respVy",1e9) # Shear response in y direction.
respVz= typical_materials.defElasticMaterial(preprocessor, "respVz",1e9) # Shear response in z direction.
# Sections
import os
pth= os.path.dirname(__file__)
#print "pth= ", pth
if(not pth):
  pth= "."
execfile(pth+"/../../aux/testQuadRegion.py")

materiales= preprocessor.getMaterialHandler
quadFibers= materiales.newMaterial("fiber_section_3d","quadFibers")
fiberSectionRepr= quadFibers.getFiberSectionRepr()
Example #7
0
# Problem type
feProblem= xc.FEProblem()
preprocessor=  feProblem.getPreprocessor
nodes= preprocessor.getNodeHandler
modelSpace= predefined_spaces.StructuralMechanics2D(nodes)

nodes.defaultTag= 1 #First node number.
nod= nodes.newNodeXY(1,0)
nod= nodes.newNodeXY(1,0)

# Materials definition: uniaxial bilinear steel
fy= 2600   # yield strength
E= 2.1e6   # initial elastic tangent
b=0.001    # strain-hardening ratio: ratio between post-yield tangent and initial elastic tangent
steel= typical_materials.defSteel01(preprocessor=preprocessor,name="steel",E=E,fy=fy,b=b)


quadFibersGeom= preprocessor.getMaterialHandler.newSectionGeometry("quadFibersGeom")
y1= width/2.0
z1= depth/2.0

reinforcement= quadFibersGeom.getReinfLayers
reinforcementA= reinforcement.newStraightReinfLayer("steel")
reinforcementA.numReinfBars= 2
reinforcementA.barArea= As
reinforcementA.p1= geom.Pos2d(y0-depth/2.0,z0-width/2.0)
reinforcementA.p2= geom.Pos2d(y0+depth/2.0,z0-width/2.0)

reinforcementB= reinforcement.newStraightReinfLayer("steel")
reinforcementB.numReinfBars= 2
Example #8
0
 def defDiagD(self,preprocessor):
   ''' Returns XC uniaxial material (design values). '''
   self.materialDiagramD= typical_materials.defSteel01(preprocessor,self.nmbDiagD,self.Es,self.fyd(),self.bsh())
   self.matTagD= self.materialDiagramD.tag
   return self.materialDiagramD #30160925 was 'return self.matTagD'
Example #9
0
 def defDiagK(self,preprocessor):
   ''' Returns XC uniaxial material (characteristic values). '''
   self.materialDiagramK= typical_materials.defSteel01(preprocessor,self.nmbDiagK,self.Es,self.fyk,self.bsh())
   self.matTagK= self.materialDiagramK.tag
   return self.materialDiagramK #30160925 was 'return self.matTagK'