Skip to content

GitLab

Commits (2)
Hide whitespace changes
Inline Side-by-side
+mafe/Member2D.m
View file @ 9f1d9cb9
...@@ -52,15 +52,15 @@ classdef Member2D < mafe.Li2Element2D ...@@ -52,15 +52,15 @@ classdef Member2D < mafe.Li2Element2D
EI = self.sect.E * self.sect.Iy; EI = self.sect.E * self.sect.Iy;
% local element matrix % local element matrix
h1 = 2*EI/L; h1 = 2*EI/L;
h2 = 6*EI/L^2; h2 = 6*EI/L^2;
h3 = 12*EI/L^3; h3 = 12*EI/L^3;
h4 = EA/L; h4 = EA/L;
ele_mat_local = [ +h4, 0, 0, -h4, 0, 0 ; ele_mat_local = [ +h4, 0, 0, -h4, 0, 0 ;
0, +h3, +h2, 0, -h3, +h2 ; 0, +h3, +h2, 0, -h3, +h2 ;
0, +h2, 2*h1, 0, -h2, +h1 ; 0, +h2, 2*h1, 0, -h2, +h1 ;
-h4, 0, 0, +h4, 0, 0 ; -h4, 0, 0, +h4, 0, 0 ;
0, -h3, -h2, 0, +h3, -h2 ; 0, -h3, -h2, 0, +h3, -h2 ;
0, +h2, +h1, 0, -h2, 2*h1 ]; 0, +h2, +h1, 0, -h2, 2*h1 ];
% global element matrix % global element matrix
ele_mat = T' * (ele_mat_local) * T; ele_mat = T' * (ele_mat_local) * T;
end end
......
+mafe/RheoKelvinVoigt1D.m
View file @ 9f1d9cb9
...@@ -56,6 +56,12 @@ classdef RheoKelvinVoigt1D < mafe.Element ...@@ -56,6 +56,12 @@ classdef RheoKelvinVoigt1D < mafe.Element
% local element matrix % local element matrix
ele_mat = d*[ +1, -1; ... ele_mat = d*[ +1, -1; ...
-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 end
%% calculation of element force vector %% calculation of element force vector
function ele_vec = force(self) function ele_vec = force(self)
......