Spektrale Beschleunigungsdichte

WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

cali

Forum-Newbie


	Beiträge: 7

	Anmeldedatum: 03.01.15

	Wohnort: ---

	Version: ---

Verfasst am: 04.02.2015, 11:27 Titel: Spektrale Beschleunigungsdichte

Hallo zusammen,

ich bin momentan dabei Beschleunigungsmesswerte auszuwerten. Dazu habe ich die spektrale Beschleunigungsdichte folgendermaßen berechnet und geplottet:

Code:

fid = fopen('Dateipfad\Messreihe_data.bin','r');
[data,count]=fread(fid,[13,inf],'double'); % Anzahl Kanäle + Zeitstempel

fid2 = fopen('Dateipfad\Messreihe_obd.bin.bin','r');
[data2]=fread(fid2,[5,inf],'double'); % Anzahl Kanäle + Zeitstempel

sensitivity = 0.33;
% x - y - z
Blau3412 = [data(8,:).*(1/sensitivity); data(6,:).*(1/sensitivity); data(12,:).*(1/sensitivity); data(13,:).*(1/sensitivity)];

iFftsize = 2^12;
iFrameshift = iFftsize/4;
iReadSamplesFromFile = 2^9;
h_window = hann(iFftsize,'periodic');
% h_window = ones(iFftsize,1);
fs = 5000;

h = figure('Name','Simulation Plot Window');
k= subplot(3,1,1);
plot(data(1,:),Blau3412(1,:))
legend('Messreihe');
ylabel('Acceleration in g');
grid on;
axis tight;

k= subplot(3,1,2);
[~, F, T, P] = spectrogram(Blau3412(1,:), h_window, iFftsize-iFrameshift, iFftsize, fs);
surf(T, F(1:200), 10*log(P(1:200,:)),'edgecolor','none');
set(gca,'clim',[-200 -30])
box('off')
axis tight;
view(0, 90);
ylabel('Hz');
pos=get(gca,'pos');
set(gca,'pos',[pos(1) pos(2) pos(3) pos(4)]);
pos=get(gca,'pos');
hc=colorbar('location','northoutside','position',[pos(1) pos(2)-0.38 pos(3) 0.03]);
h = hc;

k= subplot(3,1,3);
plot(data2(1,:),data2(2,:))
axis tight;
xlabel('Zeit in Sekunden');
ylabel('RPM')
grid on;

Funktion ohne Link?

Nun wollte ich für den Wert auf dem Plot im Anhang die mittlere Beschleunigung berechnen. Sehe ich das richtig, dass ich das folgendermaßen anstelle:

Wert=10^(Wert_dB/10)

Accel=sqrt(Wert*Frequenz)

Dann komme ich für Wert_dB=-68.01 und der Frequenz=53.71 nämlich nur auf eine Beschleunigung von 0.0029. Laut dem obersten Plot müsste sie aber nahe 0.5 sein. Mache ich was falsch oder habe ich einfach nur einen Denkfehler?

Gruß Cali

test.png

Beschreibung:

Plot

Download

Dateiname:

test.png

Dateigröße:

367.33 KB

Heruntergeladen:

298 mal

Sigggi

Forum-Anfänger


	Beiträge: 46

	Anmeldedatum: 27.01.11

	Wohnort: ---

	Version: ---

Verfasst am: 11.02.2015, 17:59 Titel:

Hallo,

wir machen Schwingungsmessungen an Maschinen. Als Mittelwert für eine Messreihe nehmen wir immer den quadratischen Mittelwert aus den (gefilterten) Messdaten im Zeitbereich:

Code:

rms = sqrt(data.^2/length(data));

Funktion ohne Link?

Oder geht's Dir um den mittleren Peak bei Drehfrequenz?

Gruß, Siggi.

cali

Themenstarter

Forum-Newbie


	Beiträge: 7

	Anmeldedatum: 03.01.15

	Wohnort: ---

	Version: ---

Verfasst am: 12.02.2015, 00:25 Titel: Beschleunigungsanteil

Hallo Siggi,

vielen Dank für deine Antwort. Ich versuche momentan die Beschleunigung für die stärker ausgeprägten Frequenzen zu berechnen. Sprich ich möchte den Anteil der Gesamtbeschleunigung wissen, der in meinem angehängten Plot von den gelben Frequenzen zustande kommt.

Gruß Cali

untitled.png

Beschreibung:

Download

Dateiname:

untitled.png

Dateigröße:

113.82 KB

Heruntergeladen:

292 mal

Sigggi

Forum-Anfänger


	Beiträge: 46

	Anmeldedatum: 27.01.11

	Wohnort: ---

	Version: ---

Verfasst am: 12.02.2015, 10:23 Titel:

Hallo,

Du könntest in P (leider hab' ich die Signal Processing Toolbox und damit spectrogram nicht) - das sind wahrscheinlich die absoluten Pegel - in einem Frequenzband nach dem höchsten Wert suchen. In den Bildern siehts so aus, als wäre die maximale Anregung bei doppelter Drehfrequenz.

Code:

for k=1:numel(T)
% doppelte Drehfrequenz
fx = 2*data(2,:)/60;

% Index
[~,ind] = min(abs(F-fx));

% Pegel im Frequenzband (mit der Breite muss man etwas experimentieren
a_max(T) = sum(P(ind-10:ind+10,k));

% bezogen auf das gesamte Spektrum
a_max_rel(T) = sum(P(1:200,k))/a_max(T);
end

Funktion ohne Link?

Das sind allerdings nur die Pegelverhältnisse. Ein genaueres Bild ergibt sich möglicherweise mit dem Effektivwert aus dem Zeitbereich ...


Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:	Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen

Impressum | Nutzungsbedingungen | Datenschutz | FAQ |

RSS

Hosted by:

Copyright © 2007 - 2025 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.