Clustered Energy Based Modal Assurance Criterion (CEMAC)
Clustered Energy Based Modal Assurance Criterion (CEMAC)
Clustered Energy Based Modal Assurance Criterion (CEMAC)
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.
load('K.mat','K')
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
M6=V'*M*V;
figure()
imagesc(M6)
M7=V'*K*V;
figure()
imagesc(M7)
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
cl={1:2:20,2:2:20};
Define eigenmode ID
j=1;
Define cluster ID
k=2;
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
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