Clustered Energy Based Modal Assurance Criterion (CEMAC)

Download as pdf or txt
Download as pdf or txt
You are on page 1of 4

Clustered Energy Based Modal Assurance Criterion (MAC)

Clustered Energy Based Modal Assurance Criterion (MAC)

Contents
Reference
Description
Structural eigenvalue analysis
Modal assurance criterion
Copyright

Reference
Brehm, M., Zabel, V., & Bucher, C. (2010). An automatic mode pairing strategy using an enhanced modal assurance criterion based on
modal strain energies. Journal of Sound and Vibration, 329(25), 5375-5392.

https://doi.org/10.1016/j.jsv.2010.07.006

Description
The above reference proposes a mode pairing strategy using an enhanced modal assurance criterion based on modal strain energies. The
procedure analyzed in section 3 (Mode assignment using energy-based modal assurance criterion) of this reference is programmed in this
submission. This verification script verifies a specific result of section 4 (Benchmark study: cantilever truss), i.e. the fraction of the total
modal strain energy of numerical eigenmode 1 that corresponds to vertical degrees of freedom, shown in Figure 7.

Structural eigenvalue analysis


Load stiffness matrix

load('K.mat','K')

Load stiffness matrix

load('M.mat','M')

Eigenvalue analysis

[V,D]=eig(K,M);

Eigenfrequencies (cycles/time)

f=sqrt(diag(D))/(2*pi);
f(1)

ans =

31.4975

file:///C/Users/pc/Desktop/CEMAC/html/verification.html[2/27/2023 12:24:41 AM]


Clustered Energy Based Modal Assurance Criterion (MAC)

Check equation (6)

M6=V'*M*V;
figure()
imagesc(M6)

Check equation (7)

M7=V'*K*V;
figure()
imagesc(M7)

file:///C/Users/pc/Desktop/CEMAC/html/verification.html[2/27/2023 12:24:41 AM]


Clustered Energy Based Modal Assurance Criterion (MAC)

Modal assurance criterion


Apply equation (1)

mac=zeros(size(V,2));
for i=1:size(V,2)
for j=1:size(V,2)
mac(i,j)=(abs(V(:,i)'*V(:,j)))^2/((V(:,i)'*V(:,i))*(V(:,j)'*V(:,j)));
end
end

Define clustering (horizontal and vertical DOFs)

cl={1:2:20,2:2:20};

Define eigenmode ID

j=1;

Define cluster ID

k=2;

file:///C/Users/pc/Desktop/CEMAC/html/verification.html[2/27/2023 12:24:41 AM]


Clustered Energy Based Modal Assurance Criterion (MAC)

Apply equation (10)

MSEjk=0;
for l=1:numel(cl)
MSEjk=MSEjk+1/2*V(cl{k},j)'*K(cl{k},cl{l})*V(cl{l},j);
end

Apply equation (11)

MSEj=0;
for k=1:numel(cl)
for l=1:numel(cl)
MSEj=MSEj+1/2*V(cl{k},j)'*K(cl{k},cl{l})*V(cl{l},j);
end
end

Apply equation (12). Compare the result with the relative modal strain energy for vertical degrees of freedom of the 1st numerical mode of
Figure 7 of the above reference.

PI_jk=MSEjk/MSEj

PI_jk =

0.9720

Copyright
Copyright (c) 2023 by George Papazafeiropoulos

Major, Infrastructure Engineer, Hellenic Air Force


Civil Engineer, M.Sc., Ph.D.
Email: [email protected]

Published with MATLAB® R2017b

file:///C/Users/pc/Desktop/CEMAC/html/verification.html[2/27/2023 12:24:41 AM]

You might also like