Newer
Older
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
close all;
clear all;
clc;
hold on;
% equi-dimensional MFD
mfd = csvread('mfd/mfd_geiger_blocking_0_0.1_to_0.9_1.csv', 1, 0);
plot(mfd(:, 6), mfd(:, 1), 'Color', 'k', 'LineWidth', 1);
% Box-DFM
dfm = csvread('boxdfm/boxdfm_geiger_blocking_0_0.1_to_0.9_1.csv', 1, 0);
plot(dfm(:, 14), dfm(:, 4), 'Color', 'g', 'LineWidth', 1);
% CC-DFM
ccdfm = csvread('ccdfm/ccdfm_geiger_blocking_0_0.1_to_0.9_1.csv', 1, 0);
plot( sqrt(ccdfm(:, 4).^2 + ccdfm(:,5).^2)-0.1, ccdfm(:, 1), 'Color', 'c', 'LineWidth', 1);
% EDFM
edfm = csvread('edfm/edfm_geiger_blocking_0_0.1_to_0.9_1.csv', 1, 0);
plot( sqrt(edfm(:, 4).^2 + edfm(:,5).^2)-0.1, edfm(:, 1), 'Color', 'r', 'LineWidth', 1);
% Mortar-DFM
mortar = csvread('mortardfm/mortardfm_geiger_blocking_0_0.1_to_0.9_1.csv', 1, 0);
plot( sqrt(mortar(:, 4).^2 + mortar(:,5).^2)-0.1, mortar(:, 1), 'Color', [1 0.7 0], 'LineWidth', 1);
% P-XFEM
p_xfem = csvread('pxfem/pxfem_geiger_blocking_0_0.1_to_0.9_1.csv', 1, 0);
plot(p_xfem(:, 3), p_xfem(:, 1), 'Color', 'b', 'LineWidth', 1);
% D-XFEM
d_xfem_fine = csvread('dxfem/dxfem_geiger_blocking_0_0.1_to_0.9_1.csv', 1, 0);
plot(d_xfem_fine(:, 3), d_xfem_fine(:, 1), 'Color', 'm', 'LineWidth', 1);
% format the plot
xlabel('arc length [m]', 'FontSize', 12, 'FontWeight', 'bold', 'Color', 'k');
ylabel('h [m]', 'FontSize', 12, 'FontWeight', 'bold', 'Color', 'k');
set(gca, 'FontSize', 12);
legend('reference', 'Box-DFM', 'CC-DFM', 'EDFM', 'Mortar-DFM', 'P-XFEM', 'D-XFEM', 'location', 'northeastoutside');
%plot2svg('hydrocoin_paper.svg');