Importing A Quadratic Mesh (GeometryFromMesh)
4 次查看(过去 30 天)
显示 更早的评论
I'm using the PDE Toolbox and the geometryFromMesh function to import meshes made externally (gmsh). The PDE Toolbox supports both linear and quadratic meshes of tetrahedra. What I'm struggling with is successfully importing an external quadratic mesh.
So I have a geometry defined inside gmsh and can produce meshes of both orders on said geometry. If I create a linear mesh and import it via geometry from mesh as:
Single_coil_3v_2split
model = createpde(1)
nodes = msh.POS'; elements = msh.TETS; groupsID = elements(:,5);
elements = elements(:,1:4); elements =elements';
geometryFromMesh(model,nodes,elements,groupsID);
pdegplot(model,'CellLabels','on','FaceAlpha',0.5)
The geometry is detected correctly. However for a quadratic mesh on the same geometry, I end up with holes & distorted elements:
Single_coil_3v_splitquad
model = createpde(1)
nodes = msh.POS'; elements = msh.TETS; groupsID = elements(:,11);
elements = msh.TETS10; elements = elements(:,1:10); elements =elements';
geometryFromMesh(model,nodes,elements,groupsID);
pdegplot(model,'CellLabels','on','FaceAlpha',0.5)
This seems to be constistently a problem regardless of how I mesh (higher/lower resolutions still have the same fundamental issue). Ideally I'd import the geometry for the linear mesh, then import the quadratic mesh however this isn't supported by the toolbox from what I can tell.
Does anyone have any ideas, solutions or suggestions?
0 个评论
回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Geometry and Mesh 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!