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

ifspec.m
1 % lese integrale Lasten aus sea_substr.plt
2
3
%
%
berechne Spektrum der Lasten
lese Reaktionskräfte aus ANSYS_Knotendatei
4 % berechne Spektrum der Reaktionen
5
6
% schreibe Tecplot-Datei
% TK
7
8 function [IFspec,nFspec] = ifspec(M)
9
10 % IFspec : Spektrum der Seegangslasten [Frequenz, X-Amplitude, Y-Amplitude,
Z-Amplitude]
11 % nFspec : Spektrum der Lagerreaktionen [ Frequenz, X-Amplitude, Y-Amplitude,
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
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
Z-Amplitude]
pathname = uigetdir('','Verzeichnis wählen');
if pathname == 0
IF = 0;
return
end
[IF] = readsubstrif(pathname);
if nargin < 1
M = size(IF,1);
end
for i = 2 : 4
spec = SpectralAnalysis(IF(:,i),IF(2,1)-IF(1,1),'AMP',M);
IFspec(:,i) = spec(:,2);
end
IFspec(:,1) = spec(:,1);
[savefile, savepath] = uiputfile({'*.plt','Tecplot-Datei (*.plt)';'*.*','Alle Dateien
(*.*)'},'Speichern',strcat(pathname,'\if_spec',num2str(M),'.plt'));
if savefile == 0
return
else
savefile = strcat(savepath,savefile);
fid = fopen(savefile,'w');
header={'VARIABLES = "f [Hz]"', '"X-Amplitude [kNs]"', '"Y-Amplitude [kNs]"',
'"Z-Amplitudfprintf(fid,'%s\n',header{:}); fprintf(fid,'ZONE T= "substr"\n');
[kNs]"'};
for i = 1 : size(IFspec,1)
fprintf(fid,'%.6f\t%.6f\t%.6f\t%.6f\n',IFspec(i,:));
end
end
fclose(fid);
[filename, pathname] = uigetfile( {'*.plt','Tecplot-Datei (*.plt)';'*.*','Alle Dateien
(*.*)'},'ANSYS-Knotendaten öffnen', strcat(pathname,'\'));
if filename == 0
return
end
ndata = readplt(strcat(pathname,filename),9,13);
M = min([M,size(ndata,1)]);
for i = 2 : 4
spec = SpectralAnalysis(ndata(:,i+6)/1000,ndata(2,1)-ndata(1,1),'AMP',M);
nFspec(:,i) = spec(:, 2);
end
nFspec(:,1) = spec(:,1);
fid = fopen(savefile,'a');
fprintf(fid,'ZONE T= "ANSYS"\n');
for i = 1 : size(nFspec,1)
fprintf(fid,'%.6f\t%.6f\t%.6f\t%.6f\n',nFspec(i,:));
end
fclose(fid);
25