XFEM-A.Posteriori.Error.Estimation/enrElem.m at main · NavidCOMSC/XFEM-A.Posteriori.Error.Estimation · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
% Written By: Matthew Jon Pais, University of Florida (2009)
% Website: http://sites.google.com/site/matthewjpais/Home
% Email: mpais@ufl.edu, matthewjpais@gmail.com

function [enrElem,IElem] = enrElem(i,pHDOF)
% This function finds the enriched elements.

global CONNEC NODES

% Find the elements consisting of the crack tip domain
I = find(NODES(:,4) > 0);                                                   % Nodes containing crack tip enrichment
[a,b,c1]    = intersect(I,CONNEC(:,2)');                                    % Find elements with crack tip enrichment
[a,b,c2]    = intersect(I,CONNEC(:,3)');                                    % Find elements with crack tip enrichment
[a,b,c3]    = intersect(I,CONNEC(:,4)');                                    % Find elements with crack tip enrichment
[a,b,c4]    = intersect(I,CONNEC(:,5)');                                    % Find elements with crack tip enrichment
CTElem      = unique([c1 c2 c3 c4]);                                        % Elements containing crack tip enrichment

if i == 1                                                                   % First iteration
    % Find the elements consiting of the crack body domain
    I = find(NODES(:,2) > 0);                                               % Nodes containing crack tip enrichment
    [a,b,c1]    = intersect(I,CONNEC(:,2)');                                % Find elements with crack tip enrichment
    [a,b,c2]    = intersect(I,CONNEC(:,3)');                                % Find elements with crack tip enrichment
    [a,b,c3]    = intersect(I,CONNEC(:,4)');                                % Find elements with crack tip enrichment
    [a,b,c4]    = intersect(I,CONNEC(:,5)');                                % Find elements with crack tip enrichment
    HElem       = unique([c1 c2 c3 c4]);                                    % Elements containing crack tip enrichment

    % Find the inclusion elements
    I = find(NODES(:,30) > 0);                                              % Nodes containing crack tip enrichment
    [a,b,c1]    = intersect(I,CONNEC(:,2)');                                % Find elements with inclusion enrichment
    [a,b,c2]    = intersect(I,CONNEC(:,3)');                                % Find elements with inclusion enrichment
    [a,b,c3]    = intersect(I,CONNEC(:,4)');                                % Find elements with inclusion enrichment
    [a,b,c4]    = intersect(I,CONNEC(:,5)');                                % Find elements with inclusion enrichment
    IElem       = unique([c1 c2 c3 c4]);                                    % Elements containing inclusion enrichment
else
    % Find the new elements containing crack body domain
    I = find(NODES(:,2) > pHDOF/2);                                         % Nodes containing crack tip enrichment
    [a,b,c1]    = intersect(I,CONNEC(:,2)');                                % Find elements with crack tip enrichment
    [a,b,c2]    = intersect(I,CONNEC(:,3)');                                % Find elements with crack tip enrichment
    [a,b,c3]    = intersect(I,CONNEC(:,4)');                                % Find elements with crack tip enrichment
    [a,b,c4]    = intersect(I,CONNEC(:,5)');                                % Find elements with crack tip enrichment
    HElem       = unique([c1 c2 c3 c4]);                                    % Elements containing crack tip enrichment

    % Return empty inclusion elements
    IElem = [];
end

% Remove duplicate entries
enrElem     = unique([CTElem HElem]);