Hello,
I want to use a mesh I created in COMSOL in FEniCSx. I exported it to a NASTRAN file, but I cannot use the mesh right away since it is missing the boundaries information. So when I open it in gmsh, and hover over the boundary lines, I don’t have anything information popping up. When I add these boundaries manually in python , it still does not appear. Using the help of ChatGPT, I managed to get the boundaries but I have these new “yellow lines” that pop up. I include some of my code below. I can’t help but think there is an easier way to do this?
Thanks. I also upload a picture of the mesh with the yellow lines. .
Code:
# Add the rectangle manually using known dimensions
** p0 = gmsh.model.geo.addPoint(0, 0, 0)**
** p1 = gmsh.model.geo.addPoint(length, 0, 0)**
** p2 = gmsh.model.geo.addPoint(length, height, 0)**
** p3 = gmsh.model.geo.addPoint(0, height, 0)**
** # Manually create lines (curves) for the boundaries**
** l0 = gmsh.model.geo.addLine(p0, p1)**
** l1 = gmsh.model.geo.addLine(p1, p2)**
** l2 = gmsh.model.geo.addLine(p2, p3)**
** l3 = gmsh.model.geo.addLine(p3, p0)**
** # Add physical groups for the boundary lines**
** gmsh.model.addPhysicalGroup(1, [l0], 1) # Bottom**
** gmsh.model.addPhysicalGroup(1, [l1], 2) # Right**
** gmsh.model.addPhysicalGroup(1, [l2], 3) # Top**
** gmsh.model.addPhysicalGroup(1, [l3], 4) # Left**
** # Create a surface for the domain**
** line_loop = gmsh.model.geo.addCurveLoop([l0, l1, l2, l3])**
** surface = gmsh.model.geo.addPlaneSurface([line_loop])**
** # Add the surface to a physical group**
** gmsh.model.addPhysicalGroup(2, [surface], 5) # Entire domain**
** print(“Manually added corner points, boundaries, and surface.”)**