|
1 | 1 | close all; clear all;
|
2 | 2 |
|
3 |
| -coords1 = csvread('01-refe.csv'); |
| 3 | +coords1 = csvread('clamped-1-refe.csv'); |
4 | 4 | N1 = length(coords1);
|
5 | 5 | n1 = N1/3;
|
6 | 6 | x1 = coords1(:,1);
|
7 | 7 | y1 = coords1(:,2);
|
8 | 8 | z1 = coords1(:,3);
|
9 | 9 | ind1 = find(x1 == min(x1) & y1 == min(y1));
|
10 |
| -coords2 = csvread('02-refe.csv'); |
| 10 | +coords2 = csvread('clamped-2-refe.csv'); |
11 | 11 | N2 = length(coords2);
|
12 | 12 | n2 = N2/3;
|
13 | 13 | x2 = coords2(:,1);
|
|
26 | 26 | ctr = 1;
|
27 | 27 | for i=0:100:10000
|
28 | 28 | if (i < 10)
|
29 |
| - disp1_file_name = strcat('01-disp-000', num2str(i), '.csv'); |
30 |
| - velo1_file_name = strcat('01-velo-000', num2str(i), '.csv'); |
31 |
| - acce1_file_name = strcat('01-acce-000', num2str(i), '.csv'); |
32 |
| - time1_file_name = strcat('01-time-000', num2str(i), '.csv'); |
33 |
| - disp2_file_name = strcat('02-disp-000', num2str(i), '.csv'); |
34 |
| - velo2_file_name = strcat('02-velo-000', num2str(i), '.csv'); |
35 |
| - acce2_file_name = strcat('02-acce-000', num2str(i), '.csv'); |
36 |
| - time2_file_name = strcat('02-time-000', num2str(i), '.csv'); |
37 |
| - dispRef1_file_name = strcat('../../dynamic-opinf-fom-neohookean/01-disp-000', num2str(i), '.csv'); |
38 |
| - veloRef1_file_name = strcat('../../dynamic-opinf-fom-neohookean/01-velo-000', num2str(i), '.csv'); |
39 |
| - acceRef1_file_name = strcat('../../dynamic-opinf-fom-neohookean/01-acce-000', num2str(i), '.csv'); |
40 |
| - dispRef2_file_name = strcat('../../dynamic-opinf-fom-neohookean/02-disp-000', num2str(i), '.csv'); |
41 |
| - veloRef2_file_name = strcat('../../dynamic-opinf-fom-neohookean/02-velo-000', num2str(i), '.csv'); |
42 |
| - acceRef2_file_name = strcat('../../dynamic-opinf-fom-neohookean/02-acce-000', num2str(i), '.csv'); |
| 29 | + disp1_file_name = strcat('clamped-1-disp-000', num2str(i), '.csv'); |
| 30 | + velo1_file_name = strcat('clamped-1-velo-000', num2str(i), '.csv'); |
| 31 | + acce1_file_name = strcat('clamped-1-acce-000', num2str(i), '.csv'); |
| 32 | + time1_file_name = strcat('clamped-1-time-000', num2str(i), '.csv'); |
| 33 | + disp2_file_name = strcat('clamped-2-disp-000', num2str(i), '.csv'); |
| 34 | + velo2_file_name = strcat('clamped-2-velo-000', num2str(i), '.csv'); |
| 35 | + acce2_file_name = strcat('clamped-2-acce-000', num2str(i), '.csv'); |
| 36 | + time2_file_name = strcat('clamped-2-time-000', num2str(i), '.csv'); |
| 37 | + dispRef1_file_name = strcat('../../dynamic-opinf-fom-neohookean/clamped-1-disp-000', num2str(i), '.csv'); |
| 38 | + veloRef1_file_name = strcat('../../dynamic-opinf-fom-neohookean/clamped-1-velo-000', num2str(i), '.csv'); |
| 39 | + acceRef1_file_name = strcat('../../dynamic-opinf-fom-neohookean/clamped-1-acce-000', num2str(i), '.csv'); |
| 40 | + dispRef2_file_name = strcat('../../dynamic-opinf-fom-neohookean/clamped-2-disp-000', num2str(i), '.csv'); |
| 41 | + veloRef2_file_name = strcat('../../dynamic-opinf-fom-neohookean/clamped-2-velo-000', num2str(i), '.csv'); |
| 42 | + acceRef2_file_name = strcat('../../dynamic-opinf-fom-neohookean/clamped-2-acce-000', num2str(i), '.csv'); |
43 | 43 | elseif (i < 100)
|
44 |
| - disp1_file_name = strcat('01-disp-00', num2str(i), '.csv'); |
45 |
| - velo1_file_name = strcat('01-velo-00', num2str(i), '.csv'); |
46 |
| - acce1_file_name = strcat('01-acce-00', num2str(i), '.csv'); |
47 |
| - time1_file_name = strcat('01-time-00', num2str(i), '.csv'); |
48 |
| - disp2_file_name = strcat('02-disp-00', num2str(i), '.csv'); |
49 |
| - velo2_file_name = strcat('02-velo-00', num2str(i), '.csv'); |
50 |
| - acce2_file_name = strcat('02-acce-00', num2str(i), '.csv'); |
51 |
| - time2_file_name = strcat('02-time-00', num2str(i), '.csv'); |
52 |
| - dispRef1_file_name = strcat('../../dynamic-opinf-fom-neohookean/01-disp-00', num2str(i), '.csv'); |
53 |
| - veloRef1_file_name = strcat('../../dynamic-opinf-fom-neohookean/01-velo-00', num2str(i), '.csv'); |
54 |
| - acceRef1_file_name = strcat('../../dynamic-opinf-fom-neohookean/01-acce-00', num2str(i), '.csv'); |
55 |
| - dispRef2_file_name = strcat('../../dynamic-opinf-fom-neohookean/02-disp-00', num2str(i), '.csv'); |
56 |
| - veloRef2_file_name = strcat('../../dynamic-opinf-fom-neohookean/02-velo-00', num2str(i), '.csv'); |
57 |
| - acceRef2_file_name = strcat('../../dynamic-opinf-fom-neohookean/02-acce-00', num2str(i), '.csv'); |
| 44 | + disp1_file_name = strcat('clamped-1-disp-00', num2str(i), '.csv'); |
| 45 | + velo1_file_name = strcat('clamped-1-velo-00', num2str(i), '.csv'); |
| 46 | + acce1_file_name = strcat('clamped-1-acce-00', num2str(i), '.csv'); |
| 47 | + time1_file_name = strcat('clamped-1-time-00', num2str(i), '.csv'); |
| 48 | + disp2_file_name = strcat('clamped-2-disp-00', num2str(i), '.csv'); |
| 49 | + velo2_file_name = strcat('clamped-2-velo-00', num2str(i), '.csv'); |
| 50 | + acce2_file_name = strcat('clamped-2-acce-00', num2str(i), '.csv'); |
| 51 | + time2_file_name = strcat('clamped-2-time-00', num2str(i), '.csv'); |
| 52 | + dispRef1_file_name = strcat('../../dynamic-opinf-fom-neohookean/clamped-1-disp-00', num2str(i), '.csv'); |
| 53 | + veloRef1_file_name = strcat('../../dynamic-opinf-fom-neohookean/clamped-1-velo-00', num2str(i), '.csv'); |
| 54 | + acceRef1_file_name = strcat('../../dynamic-opinf-fom-neohookean/clamped-1-acce-00', num2str(i), '.csv'); |
| 55 | + dispRef2_file_name = strcat('../../dynamic-opinf-fom-neohookean/clamped-2-disp-00', num2str(i), '.csv'); |
| 56 | + veloRef2_file_name = strcat('../../dynamic-opinf-fom-neohookean/clamped-2-velo-00', num2str(i), '.csv'); |
| 57 | + acceRef2_file_name = strcat('../../dynamic-opinf-fom-neohookean/clamped-2-acce-00', num2str(i), '.csv'); |
58 | 58 | elseif (i < 1000)
|
59 |
| - disp1_file_name = strcat('01-disp-0', num2str(i), '.csv'); |
60 |
| - velo1_file_name = strcat('01-velo-0', num2str(i), '.csv'); |
61 |
| - acce1_file_name = strcat('01-acce-0', num2str(i), '.csv'); |
62 |
| - time1_file_name = strcat('01-time-0', num2str(i), '.csv'); |
63 |
| - disp2_file_name = strcat('02-disp-0', num2str(i), '.csv'); |
64 |
| - velo2_file_name = strcat('02-velo-0', num2str(i), '.csv'); |
65 |
| - acce2_file_name = strcat('02-acce-0', num2str(i), '.csv'); |
66 |
| - time2_file_name = strcat('02-time-0', num2str(i), '.csv'); |
67 |
| - dispRef1_file_name = strcat('../../dynamic-opinf-fom-neohookean/01-disp-0', num2str(i), '.csv'); |
68 |
| - veloRef1_file_name = strcat('../../dynamic-opinf-fom-neohookean/01-velo-0', num2str(i), '.csv'); |
69 |
| - acceRef1_file_name = strcat('../../dynamic-opinf-fom-neohookean/01-acce-0', num2str(i), '.csv'); |
70 |
| - dispRef2_file_name = strcat('../../dynamic-opinf-fom-neohookean/02-disp-0', num2str(i), '.csv'); |
71 |
| - veloRef2_file_name = strcat('../../dynamic-opinf-fom-neohookean/02-velo-0', num2str(i), '.csv'); |
72 |
| - acceRef2_file_name = strcat('../../dynamic-opinf-fom-neohookean/02-acce-0', num2str(i), '.csv'); |
| 59 | + disp1_file_name = strcat('clamped-1-disp-0', num2str(i), '.csv'); |
| 60 | + velo1_file_name = strcat('clamped-1-velo-0', num2str(i), '.csv'); |
| 61 | + acce1_file_name = strcat('clamped-1-acce-0', num2str(i), '.csv'); |
| 62 | + time1_file_name = strcat('clamped-1-time-0', num2str(i), '.csv'); |
| 63 | + disp2_file_name = strcat('clamped-2-disp-0', num2str(i), '.csv'); |
| 64 | + velo2_file_name = strcat('clamped-2-velo-0', num2str(i), '.csv'); |
| 65 | + acce2_file_name = strcat('clamped-2-acce-0', num2str(i), '.csv'); |
| 66 | + time2_file_name = strcat('clamped-2-time-0', num2str(i), '.csv'); |
| 67 | + dispRef1_file_name = strcat('../../dynamic-opinf-fom-neohookean/clamped-1-disp-0', num2str(i), '.csv'); |
| 68 | + veloRef1_file_name = strcat('../../dynamic-opinf-fom-neohookean/clamped-1-velo-0', num2str(i), '.csv'); |
| 69 | + acceRef1_file_name = strcat('../../dynamic-opinf-fom-neohookean/clamped-1-acce-0', num2str(i), '.csv'); |
| 70 | + dispRef2_file_name = strcat('../../dynamic-opinf-fom-neohookean/clamped-2-disp-0', num2str(i), '.csv'); |
| 71 | + veloRef2_file_name = strcat('../../dynamic-opinf-fom-neohookean/clamped-2-velo-0', num2str(i), '.csv'); |
| 72 | + acceRef2_file_name = strcat('../../dynamic-opinf-fom-neohookean/clamped-2-acce-0', num2str(i), '.csv'); |
73 | 73 | else
|
74 |
| - disp1_file_name = strcat('01-disp-', num2str(i), '.csv'); |
75 |
| - velo1_file_name = strcat('01-velo-', num2str(i), '.csv'); |
76 |
| - acce1_file_name = strcat('01-acce-', num2str(i), '.csv'); |
77 |
| - time1_file_name = strcat('01-time-', num2str(i), '.csv'); |
78 |
| - disp2_file_name = strcat('02-disp-', num2str(i), '.csv'); |
79 |
| - velo2_file_name = strcat('02-velo-', num2str(i), '.csv'); |
80 |
| - acce2_file_name = strcat('02-acce-', num2str(i), '.csv'); |
81 |
| - time2_file_name = strcat('02-time-', num2str(i), '.csv'); |
82 |
| - dispRef1_file_name = strcat('../../dynamic-opinf-fom-neohookean/01-disp-', num2str(i), '.csv'); |
83 |
| - veloRef1_file_name = strcat('../../dynamic-opinf-fom-neohookean/01-velo-', num2str(i), '.csv'); |
84 |
| - acceRef1_file_name = strcat('../../dynamic-opinf-fom-neohookean/01-acce-', num2str(i), '.csv'); |
85 |
| - dispRef2_file_name = strcat('../../dynamic-opinf-fom-neohookean/02-disp-', num2str(i), '.csv'); |
86 |
| - veloRef2_file_name = strcat('../../dynamic-opinf-fom-neohookean/02-velo-', num2str(i), '.csv'); |
87 |
| - acceRef2_file_name = strcat('../../dynamic-opinf-fom-neohookean/02-acce-', num2str(i), '.csv'); |
| 74 | + disp1_file_name = strcat('clamped-1-disp-', num2str(i), '.csv'); |
| 75 | + velo1_file_name = strcat('clamped-1-velo-', num2str(i), '.csv'); |
| 76 | + acce1_file_name = strcat('clamped-1-acce-', num2str(i), '.csv'); |
| 77 | + time1_file_name = strcat('clamped-1-time-', num2str(i), '.csv'); |
| 78 | + disp2_file_name = strcat('clamped-2-disp-', num2str(i), '.csv'); |
| 79 | + velo2_file_name = strcat('clamped-2-velo-', num2str(i), '.csv'); |
| 80 | + acce2_file_name = strcat('clamped-2-acce-', num2str(i), '.csv'); |
| 81 | + time2_file_name = strcat('clamped-2-time-', num2str(i), '.csv'); |
| 82 | + dispRef1_file_name = strcat('../../dynamic-opinf-fom-neohookean/clamped-1-disp-', num2str(i), '.csv'); |
| 83 | + veloRef1_file_name = strcat('../../dynamic-opinf-fom-neohookean/clamped-1-velo-', num2str(i), '.csv'); |
| 84 | + acceRef1_file_name = strcat('../../dynamic-opinf-fom-neohookean/clamped-1-acce-', num2str(i), '.csv'); |
| 85 | + dispRef2_file_name = strcat('../../dynamic-opinf-fom-neohookean/clamped-2-disp-', num2str(i), '.csv'); |
| 86 | + veloRef2_file_name = strcat('../../dynamic-opinf-fom-neohookean/clamped-2-velo-', num2str(i), '.csv'); |
| 87 | + acceRef2_file_name = strcat('../../dynamic-opinf-fom-neohookean/clamped-2-acce-', num2str(i), '.csv'); |
88 | 88 | end
|
89 | 89 | d1 = dlmread(disp1_file_name);
|
90 | 90 | v1 = dlmread(velo1_file_name);
|
|
0 commit comments