Beispiel #1
0
        private void _correctMergedElements()
        {
            Console.WriteLine("Correcting merged elements...");

            int numElements     = Elements.Count;
            int correctionCount = 0;

            for (int i = 0; i < numElements; ++i)
            {
                GmshMeshElement e = Elements[i];
                if (GmshMeshElement.MergedElementConversionTable.ContainsKey(e.ElementType))
                {
                    // This is an element that may be corrected
                    GmshMeshElementType correctedType = GmshMeshElement.MergedElementConversionTable[e.ElementType];
                    HashSet <int>       uniqueIds     = e.GetUniqueNodeIds();
                    if (uniqueIds.Count != GmshMeshElement.ElementTypeToNumNodes[e.ElementType])
                    {
                        // There are merged nodes
                        if (uniqueIds.Count == GmshMeshElement.ElementTypeToNumNodes[correctedType])
                        {
                            // The merging is consistent with the correction
                            GmshMeshElement correctedElement = new GmshMeshElement(correctedType, e.ID, uniqueIds.ToList(), e.ElementaryTag, e.PhysicalTag);
                            Elements[i] = correctedElement;
                            ++correctionCount;
                        }
                        else
                        {
                            Console.WriteLine(string.Format("Input element id {0} has {1} unique nodes, inconsistent with corrected element that expects {2} unique nodes. Skipping.", e.ID, uniqueIds.Count, GmshMeshElement.ElementTypeToNumNodes[correctedType]));
                        }
                    }
                }
            }

            Console.WriteLine(string.Format("Corrected {0} elements.", correctionCount));
        }
Beispiel #2
0
        private void _tetrahedralizeMergedQuads()
        {
            Console.WriteLine("Tetrahedralizing...");

            int numElements           = Elements.Count;
            int triangulatedQuads     = 0;
            int tetrahedrealizedHexes = 0;

            for (int i = 0; i < numElements; ++i)
            {
                GmshMeshElement e = Elements[i];

                switch (e.ElementType)
                {
                case GmshMeshElementType.QUAD_4NODE:
                {
                    HashSet <int> uniqueIds = e.GetUniqueNodeIds();

                    if (uniqueIds.Count != GmshMeshElement.ElementTypeToNumNodes[GmshMeshElementType.QUAD_4NODE])
                    {
                        // There are repeated nodes
                        if (uniqueIds.Count == GmshMeshElement.ElementTypeToNumNodes[GmshMeshElementType.TRIANGLE_3NODE])
                        {
                            GmshMeshTri3NodeElement eTri3Node = new GmshMeshTri3NodeElement(e.ID, uniqueIds.ToList(), e.ElementaryTag, e.PhysicalTag);
                            Elements[i] = eTri3Node;
                            ++triangulatedQuads;
                        }
                        else
                        {
                            Console.WriteLine(String.Format("Quad element id {0} has {1} unique nodes. Cannot convert to triangle.", e.ID, uniqueIds.Count));
                        }
                    }
                }
                break;

                case GmshMeshElementType.QUAD_8NODE:
                {
                    HashSet <int> uniqueIds = e.GetUniqueNodeIds();

                    if (uniqueIds.Count != GmshMeshElement.ElementTypeToNumNodes[GmshMeshElementType.QUAD_8NODE])
                    {
                        // There are repeated nodes
                        if (uniqueIds.Count == GmshMeshElement.ElementTypeToNumNodes[GmshMeshElementType.TRIANGLE_6NODE])
                        {
                            GmshMeshTri6NodeElement eTri3Node = new GmshMeshTri6NodeElement(e.ID, uniqueIds.ToList(), e.ElementaryTag, e.PhysicalTag);
                            Elements[i] = eTri3Node;
                            ++triangulatedQuads;
                        }
                        else
                        {
                            Console.WriteLine(String.Format("Quad element id {0} has {1} unique nodes. Cannot convert to triangle.", e.ID, uniqueIds.Count));
                        }
                    }
                }
                break;

                case GmshMeshElementType.HEXA_8NODE:
                {
                    HashSet <int> uniqueIds = e.GetUniqueNodeIds();

                    if (uniqueIds.Count != GmshMeshElement.ElementTypeToNumNodes[GmshMeshElementType.HEXA_8NODE])
                    {
                        // There are repeated nodes
                        if (uniqueIds.Count == GmshMeshElement.ElementTypeToNumNodes[GmshMeshElementType.TET_4NODE])
                        {
                            GmshMeshTet4NodeElement eTri3Node = new GmshMeshTet4NodeElement(e.ID, uniqueIds.ToList(), e.ElementaryTag, e.PhysicalTag);
                            Elements[i] = eTri3Node;
                            ++tetrahedrealizedHexes;
                        }
                        else
                        {
                            Console.WriteLine(String.Format("Hexa element id {0} has {1} unique nodes. Cannot convert to tetrahedron.", e.ID, uniqueIds.Count));
                        }
                    }
                }
                break;

                case GmshMeshElementType.HEXA_20NODE:
                {
                    HashSet <int> uniqueIds = e.GetUniqueNodeIds();

                    if (uniqueIds.Count != GmshMeshElement.ElementTypeToNumNodes[GmshMeshElementType.HEXA_8NODE])
                    {
                        // There are repeated nodes
                        if (uniqueIds.Count == GmshMeshElement.ElementTypeToNumNodes[GmshMeshElementType.TET_4NODE])
                        {
                            GmshMeshTet10NodeElement eTri3Node = new GmshMeshTet10NodeElement(e.ID, uniqueIds.ToList(), e.ElementaryTag, e.PhysicalTag);
                            Elements[i] = eTri3Node;
                            ++tetrahedrealizedHexes;
                        }
                        else
                        {
                            Console.WriteLine(String.Format("Hexa element id {0} has {1} unique nodes. Cannot convert to tetrahedron.", e.ID, uniqueIds.Count));
                        }
                    }
                }
                break;

                default:
                    break;
                }
            }

            Console.WriteLine("Triangulated {0} quads.", triangulatedQuads);
            Console.WriteLine("Tetrahedralized {0} hexes.", tetrahedrealizedHexes);
        }