Check and warn in case of excitations at eigenfrequency (Inf solutions).

+mafe/FeAnalysisDynamicFD.m

@@ -119,6 +119,12 @@ classdef FeAnalysisDynamicFD < mafe.FeAnalysis
                 ff = f + fa + fv + fu;
                 % calculate system response associated with time function tf
                 self.Z(act,ii) = ( -tf.Omega^2 * self.M(act,act) + 1i * tf.Omega * self.D(act,act) + self.K(act,act) ) \ ff(act);
+                % check resonance case
+                if length( find( self.Z(act,ii) == Inf ) ) > 0
+                  warning('FeAnalysisDynamicFD: Given excitation frequency of TimeFunction seems to match system eigenfrequency!');
+                  tf
+                  self.Z(act,ii) = 0; % neglect this solution
+                end
                 % put kinematic constraint values to solution
                 self.Z(dbc,ii) = fu(dbc);
                 % put local force vector