from model import predefined_spaces from model import fix_node_6dof from materials import typical_materials E= 2.1e6 # Steel Young modulus. nu= 0.3 # Poisson's coefficient. h= 0.1 # Espesor. dens= 1.33 # Densidad kg/m2. L= 2 p= 10 # Tipo de problema prueba= xc.ProblemaEF() mdlr= prueba.getModelador nodos= mdlr.getNodeLoader predefined_spaces.gdls_resist_materiales3D(nodos) nodos.newNodeIDXYZ(1,0,0,0) nodos.newNodeIDXYZ(2,L,0,0) nodos.newNodeIDXYZ(3,L,L,0) nodos.newNodeIDXYZ(4,0,L,0) # Definimos materiales memb1= typical_materials.defElasticMembranePlateSection(mdlr,"memb1",E,nu,dens,h) elementos= mdlr.getElementLoader elementos.defaultMaterial= "memb1" elem= elementos.newElement("shell_mitc4",xc.ID([1,2,3,4]))
from model import predefined_spaces from model import fix_node_6dof from materials import typical_materials E= 2.1e6 # Steel Young modulus. nu= 0.3 # Poisson's coefficient. h= 0.1 # Espesor. dens= 1.33 # Density kg/m2. L= 2 p= 10 # Problem type prueba= xc.ProblemaEF() prep= prueba.getPreprocessor nodos= prep.getNodeHandler predefined_spaces.gdls_resist_materiales3D(nodos) nodos.newNodeIDXYZ(1,0,0,0) nodos.newNodeIDXYZ(2,L,0,0) nodos.newNodeIDXYZ(3,L,L,0) nodos.newNodeIDXYZ(4,0,L,0) # Definimos materiales memb1= typical_materials.defElasticMembranePlateSection(prep,"memb1",E,nu,dens,h) elementos= prep.getElementHandler elementos.defaultMaterial= "memb1" elem= elementos.newElement("shell_mitc4",xc.ID([1,2,3,4]))