Kossel_studienarbeit.pdf 1.4 MB - Institut für Strömungsmechanik ...

74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
figure;hist(nF(:,6),0:12:360);
xlabel('nautische Richtung [°], N = 0°')
ylabel('Häufigkeit [-]')
set(gca,'XLim',[0 360])
set(gca,'XTick',0:30:360)
set(findobj(gca,'Type','patch'),'FaceColor','g')
set(gca,'XGrid','On')
set(gca,'YGrid','On')
figure;hist(nF(:,5),0:50:1500);
xlabel('Fxy [kN]')
ylabel('Häufigkeit [-]')
set(gca,'XLim',[0 1500])
set(gca,'XTick',0:200:1500)
set(findobj(gca,'Type','patch'),'FaceColor','b')
set(gca,'XGrid','On')
set(gca,'YGrid','On')
fid = fopen(savefile,'a');
fprintf(fid,'ZONE T= "ANSYS"\n');
for i = 1 : size(nF,1)
fprintf(fid,'%.6f\t%.6f\t%.6f\t%.6f\t%.6f\t%.6f\n',nF(i,:));
end
fclose(fid);
scatterplot.m
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
% lese if.plt der Rechnungen mit s=10, s= 25, s=75
% trage die Beträge der resultierenden Kräfte gegeneinander auf (Scatter-Plot)
function [data] = scatterplot
[filename, pathname] = uigetfile({'*.plt',' Tecplot-Datei (*.plt)';'*.*','Alle Dateien
(*.*)'},'if.plt s=10 öffnen','if.plt');
if filename == 0
data = 0; dt = 0; header = 0;
return
end
openfile{1} = strcat(pathname,filename);
[filename, pathname] = uigetfile({'*.plt','Tecplot-Datei (*.plt)';'*.*','Alle Dateien
(*.*)'},'if.plt s=25 öffnen','if.plt');
if filename == 0
data = 0; dt = 0; header = 0;
return
end
openfile{2} = strcat(pathname,filename);
[filename, pathname] = uigetfile({'*.plt', 'Tecplot-Datei (*.plt)';'*.*','Alle Dateien
(*.*)'},'if.plt s=75 öffnen','if.plt');
if filename == 0
data = 0; dt = 0; header = 0;
return
end
openfile{3} = strcat(pathname,filename);
for i = 1 : 3
fid = fopen(openfile{i});
for j = 1 : 7
header{j,1} = fgetl(fid);
end
rdata = fscanf(fid,'%f',[6 inf])';
end
data(:,i) = rdata(:,5);
fclose(fid);
dim = max(max(data));
figure;
plot(data(:,1),data(:,2),'.');
hold on
plot([0 dim],[0 dim]);
xlabel('s = 10, Fxy [kN]')
ylabel('s = 25, Fxy [kN]')
set(gca,'XLim',[0 dim])
set(gca,'YLim',[0 dim])
27