# 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 app = vv.use() app.Create() app.Run()
g.setCurveMarker(curveID, 80) #MESHING: elType = 2 #Element type 2 is triangle. (3 is quad. See user manual for more element types) dofsPerNode= 2 #Degrees of freedom per node. mesher = GmshMeshGenerator(g) mesher.elSizeFactor = 0.05 mesher.elType = elType mesher.dofsPerNode = 2 #Mesh the geometry: # The first four return values are the same as those that trimesh2d() returns. # value elementmarkers is a list of markers, and is used for finding the marker of a given element (index). coords, edof, dofs, bdofs, elementmarkers = mesher.create() #VISUALISATION: #Hold left mouse button to pan. #Hold right mouse button to zoom. pcv.drawGeometry(g, labelCurves=True)#Draws the geometry. vv.figure() #New figure window pcv.drawMesh(coords=coords, edof=edof, dofsPerNode=dofsPerNode, elType=elType, filled=True, title="Example 02") #Draws the mesh. vv.gca().axis.showGrid = True #Show grid # Enter main loop: app = vv.use() app.Create() app.Run()