Commit 95a1e12a authored by Andreas Zilian's avatar Andreas Zilian
Browse files

Add Rayleigh damping.

parent fba8a3b8
......@@ -56,6 +56,12 @@ classdef RheoKelvinVoigt1D < mafe.Element
% local element matrix
ele_mat = d*[ +1, -1; ...
-1, +1 ];
% provide damping matrix based on Rayleigh assumption
% get parameters
d_alpha = self.sect.d_alpha;
d_beta = self.sect.d_beta;
% D = alpha * M + beta * K
ele_mat = ele_mat + d_alpha * self.mass() + d_beta * self.stiffness();
end
%% calculation of element force vector
function ele_vec = force(self)
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment