Example #1
0
        private XSubdomain2D_old FindSubdomainWhereNodeEnrichmentIsNotSingular(XCluster2D cluster, XNode node,
                                                                               IEnrichmentItem2D enrichment)
        {
            foreach (XSubdomain2D_old subdomain in cluster.Subdomains)
            {
                if (subdomain.BoundaryNodes.Contains(node))
                {
                    bool isSingularityNode = subdomain.DofOrderer.SingularHeavisideEnrichments.TryGetValue(node,
                                                                                                           out HashSet <IEnrichmentItem2D> singularEnrichments);
                    if (!isSingularityNode)
                    {
                        return(subdomain);
                    }
                    else if (!singularEnrichments.Contains(enrichment))
                    {
                        return(subdomain);
                    }
                }
            }

            // This point should not be reached under normal circumstances
            if (SingularHeavisideEnrichments.ContainsKey(node))
            {
                throw new IncorrectDecompositionException(
                          $"Heaviside dofs of {node} are singular in all subdomains. Investigate further");
            }
            else
            {
                throw new ArgumentException(
                          $"In this subdomain, {node} is not a boundary node with singular Heaviside enrichment");
            }
        }
Example #2
0
        /// <summary>
        /// If a fixed enrichment area is applied, all nodes inside a circle around the tip are enriched with tip
        /// functions. They can still be enriched with Heaviside functions, if they do not belong to the tip
        /// element(s).
        /// </summary>
        /// <param name="tipNodes"></param>
        /// <param name="tipElement"></param>
        private void ApplyFixedEnrichmentArea(CartesianPoint crackTip, XContinuumElement2D tipElement,
                                              HashSet <XNode> tipNodes, IEnrichmentItem2D tipEnrichments)
        {
            if (enrichmentRadiusOverElementSize > 0)
            {
                var    outline        = ConvexPolygon2D.CreateUnsafe(tipElement.Nodes);
                double elementArea    = outline.ComputeArea();
                double radius         = enrichmentRadiusOverElementSize * Math.Sqrt(elementArea);
                var    enrichmentArea = new Circle2D(crackTip, radius);

                foreach (var element in Mesh.FindElementsInsideCircle(enrichmentArea, tipElement))
                {
                    bool completelyInside = true;
                    foreach (var node in element.Nodes)
                    {
                        CirclePointPosition position = enrichmentArea.FindRelativePositionOfPoint(node);
                        if ((position == CirclePointPosition.Inside) || (position == CirclePointPosition.On))
                        {
                            tipNodes.Add(node);
                        }
                        else
                        {
                            completelyInside = false;
                        }
                    }

                    if (completelyInside)
                    {
                        element.EnrichmentItems.Add(tipEnrichments);
                    }
                }
            }
        }
Example #3
0
 public void AddEnrichment(IEnrichmentItem2D enrichment)
 {
     if (enrichments.Contains(enrichment))
     {
         throw new ArgumentException("This enrichment is already inserted");
     }
     enrichments.Add(enrichment);
 }