DVD-Inhalt Dieser Studienarbeit ist eine DVD mit den durchgeführten Berechnungen und den erstellten und genutzten Matlab-Tools beigefügt. Die Verzeichnis-Struktur gliedert sich wie folgt: - Pfahl Berechnungen zu Kapitel 3 - - t1000 1D-/ 2D-Berechnung, Kapitel 3.1 und 3.2 t300 versch. Spread-Parameter, Kapitel 3.3 - t60 kurze Simulatonsdauer, nicht <strong>für</strong> Auswertung - GWK Berechnungen zu Kapitel 4, Modelle siehe Abb. 4.3 - - Feder - trocken - Cm1 - Cm2 Masse - trocken - Cm1 - Cm2 Modelle mit Biegefedern statt Einspannung Modelle mit Massenpunkt - Matlab-Tools erstellte und genutzte Matlab-Tools Weiterhin steht auf der DVD dieses Dokument im PDF-Format zur Verfügung. 23
Anhang A: Matlab-Tools spec.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 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 % Lese Wasserspiegelauslenkung aus surface1D.plt % berechne 1D-Spektrum % Lese 1D-Spektrum aus spec1d.plt bzw. Spectrum1D_in.plt % berechne m0, Hm0 <strong>für</strong> beide Spektren % schreibe Tecplot-Datei % TK function [spec,m0_1,Hm0_1] = spec(M,erw) % spe c : Spektrum aus Wasserspiegelauslenkung % m0_1, Hm0_1 : aus Spektrum aus Wasserspiegelauslenkung % % M : Glättungsparameter <strong>für</strong> Fourier-Transformation % erw : Zeitreihe der Wasserspiegelauslenkungen erweitern % => bei 0 beginnen, aufhören % erw = Anzahl der Stützstellen [filenam e, pathname] = uigetfile({'*.plt','Tecplot-Datei (*.plt)';'*.*','Alle Dateien (*.*)'},'Surface-Datei öffnen','..\surface1d.plt'); if filename == 0 spec = 0; m0_1 = 0; Hm0_1 = 0; return end [data,dt] = readsurface1D(strcat(pathname,filename)); if nargin > 1 data(erw+1:length(data)+erw,1) = data; for i = 1 : erw data(i,1) = data(erw+1,1) * (i-1) /erw; data(length(data)+1,1) = data(length(data)+1-i,1) * (erw-i) /erw; end erw = strca t('_mod',num2str(erw)); else erw = ''; end if nargin < 1 M = size(data,1); end psdt = SpectralAnalysis(data,dt,'PSD',M); m0_1 = sum(psdt(2:length(psdt),2)) * psdt(2,1); Hm0_1 = 4.0 * sqrt(m0_1); spec(:,1) = psdt(:,1) * (2.*pi); spec(:,2) = psdt(:,2) / (2.*pi); data = readspec1d(pathname); m0_0 = trapz(data(:,1),data(:,2)); Hm0_0 = 4.0 * sqrt(m0_0); [filename, newpath] = uiputfile({'*.dat','Datei (*.dat)';'*.*','Alle Dateien(*.*)'}, 'm0, Hm0 speichern',strcat(pathname,'m0_Hm0_',num2str(M),erw,'.dat')); if filename ~= 0 fid = fopen(strcat(newpath,filename),'w'); fprintf(fid,'m0_0 = %.8f\nHm0_0 = %.8f\nm0_1 = %.8f\nHm0_1 = %.8f', m0_0,Hm0_0,m0_1,Hm0_1); fclose(fid); pathname = newpath; end [filename, pathname] = uiputfile({'*.plt','Tecplot-Datei (*.plt)';'*.*','Alle Dateien(*.*)'},'Spektrum speichern',strcat(pathname,'spec_',num2str(M),erw,'.plt')); if filename ~= 0 filename = strcat(pathname,filename); writespec(spec,'w','surface1d',filename); if strcmp(questdlg('Input-Spektrum anhängen?','Frage','Ja','Nein','Ja'),'Ja') writespec(data,'a','spectrum1d_in',filename); end end 24