feProblem= xc.FEProblem()
preprocessor= feProblem.getPreprocessor
nodes= preprocessor.getNodeHandler
# Problem type
modelSpace= predefined_spaces.StructuralMechanics3D(nodes)
nodes.defaultTag= 1 #First node number.
nodes.newNodeXYZ(0,0.0,0.0)
nodes.newNodeXYZ(L,0.0,0.0)
lin= modelSpace.newLinearCrdTransf("lin",xc.Vector([0,1,0]))
# Materials
sectionProperties= xc.CrossSectionProperties3d()
sectionProperties.A= A; sectionProperties.E= E; sectionProperties.G= G;
sectionProperties.Iz= Iz; sectionProperties.Iy= Iy; sectionProperties.J= J
section= typical_materials.defElasticSectionFromMechProp3d(preprocessor, "section",sectionProperties)
# Elements definition
elements= preprocessor.getElementHandler
elements.defaultTransformation= "lin"
elements.defaultMaterial= "section"
elements.defaultTag= 1 #Tag for the next element.
beam3d= elements.newElement("ElasticBeam3d",xc.ID([1,2]));
modelSpace.fixNode000_000(1)
loadHandler= preprocessor.getLoadHandler
lPatterns= loadHandler.getLoadPatterns
#Load modulation.
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
caracMecSeccion = xc.CrossSectionProperties3d()
caracMecSeccion.A = A
caracMecSeccion.E = E
caracMecSeccion.G = G
caracMecSeccion.Iz = Iz
caracMecSeccion.Iy = Iy
caracMecSeccion.J = J
seccion = typical_materials.defElasticSectionFromMechProp3d(
preprocessor, "seccion", caracMecSeccion)
# Elements definition
elementos = preprocessor.getElementLoader
elementos.defaultTransformation = "lin"
elementos.defaultMaterial = "seccion"
elementos.defaultTag = 1 #Tag for the next element.
beam3d = elementos.newElement("elastic_beam_3d", xc.ID([1, 2]))
# Constraints
coacciones = preprocessor.getConstraintLoader
fix_node_6dof.fixNode6DOF(coacciones, 1)
# Loads definition
cargas = preprocessor.getLoadLoader
casos = cargas.getLoadPatterns
feProblem= xc.FEProblem()
preprocessor= feProblem.getPreprocessor
nodes= preprocessor.getNodeHandler
# Problem type
modelSpace= predefined_spaces.StructuralMechanics3D(nodes)
nodes.defaultTag= 1 #First node number.
nodes.newNodeXYZ(0,0.0,0.0)
nodes.newNodeXYZ(L,0.0,0.0)
lin= modelSpace.newLinearCrdTransf("lin",xc.Vector([0,0,1]))
# Materials
sectionProperties= xc.CrossSectionProperties3d()
sectionProperties.A= A; sectionProperties.E= E; sectionProperties.G= G
sectionProperties.Iz= Iz; sectionProperties.Iy= Iy; sectionProperties.J= J
section= typical_materials.defElasticSectionFromMechProp3d(preprocessor, "section",sectionProperties)
# Elements definition
elements= preprocessor.getElementHandler
elements.defaultTransformation= "lin"
elements.defaultMaterial= "section"
elements.defaultTag= 1 #Tag for the next element.
beam3d= elements.newElement("ElasticBeam3d",xc.ID([1,2]))
modelSpace.fixNode000_000(1)
elem=preprocessor.getElementHandler.getElement(1) #element 1
loadHandler= preprocessor.getLoadHandler
lPatterns= loadHandler.getLoadPatterns