g.addStructuredSurface([0,1,2,3]) g.addSurface([4,1]) # ---- Create mesh ---------------------------------------------------------- # Element type 3 is quad. (2 is triangle. See user manual for more element types) elType = 3 # Degrees of freedom per node. dofsPerNode = 1 coords, edof, dofs, bdofs, elementmarkers = cfm.createGmshMesh(g, elSizeFactor = 0.05, elType = elType, dofsPerNode = dofsPerNode) # ---- Visualise mesh ------------------------------------------------------- # Draw geometry cfv.drawGeometry(g) # Draw mesh vv.figure() cfv.drawMesh(coords=coords, edof=edof, dofsPerNode=dofsPerNode, elType=elType, filled=True) # Enter main loop
g.addStructuredVolume([0,1,2,3,4,5], 0, marker=90) # ---- Create mesh ---------------------------------------------------------- # Element type 5 is hexahedron. (See user manual for more element types) elType = 5 # Degrees of freedom per node. dofsPerNode= 1 # Create mesh coords, edof, dofs, bdofs, _ = cfm.createGmshMesh(geometry = g, elType = elType, dofsPerNode = dofsPerNode) # ---- Visualise mesh ------------------------------------------------------- # Hold Left Mouse button to rotate. # Hold right mouse button to zoom. # Hold SHIFT and left mouse button to pan. # Hold SHIFT and right mouse button to change the field of view. # Hold Ctrl and left mouse button to roll the camera. # Draw geometry cfv.drawGeometry(g, drawPoints=False) # Draw mesh