28 lines
986 B
Mathematica
28 lines
986 B
Mathematica
|
function y = BeamAssemble(K,k,i,j)
|
||
|
|
%BeamAssemble This function assembles the element stiffness
|
||
|
|
% matrix k of the beam element with nodes
|
||
|
|
% i and j into the global stiffness matrix K.
|
||
|
|
% This function returns the global stiffness
|
||
|
|
% matrix K after the element stiffness matrix
|
||
|
|
% k is assembled.
|
||
|
|
K(2*i-1,2*i-1) = K(2*i-1,2*i-1) + k(1,1);
|
||
|
|
K(2*i-1,2*i) = K(2*i-1,2*i) + k(1,2);
|
||
|
|
K(2*i-1,2*j-1) = K(2*i-1,2*j-1) + k(1,3);
|
||
|
|
K(2*i-1,2*j) = K(2*i-1,2*j) + k(1,4);
|
||
|
|
K(2*i,2*i-1) = K(2*i,2*i-1) + k(2,1);
|
||
|
|
K(2*i,2*i) = K(2*i,2*i) + k(2,2);
|
||
|
|
K(2*i,2*j-1) = K(2*i,2*j-1) + k(2,3);
|
||
|
|
K(2*i,2*j) = K(2*i,2*j) + k(2,4);
|
||
|
|
K(2*j-1,2*i-1) = K(2*j-1,2*i-1) + k(3,1);
|
||
|
|
K(2*j-1,2*i) = K(2*j-1,2*i) + k(3,2);
|
||
|
|
K(2*j-1,2*j-1) = K(2*j-1,2*j-1) + k(3,3);
|
||
|
|
K(2*j-1,2*j) = K(2*j-1,2*j) + k(3,4);
|
||
|
|
K(2*j,2*i-1) = K(2*j,2*i-1) + k(4,1);
|
||
|
|
K(2*j,2*i) = K(2*j,2*i) + k(4,2);
|
||
|
|
K(2*j,2*j-1) = K(2*j,2*j-1) + k(4,3);
|
||
|
|
K(2*j,2*j) = K(2*j,2*j) + k(4,4);
|
||
|
|
y = K;
|
||
|
|
|
||
|
|
|
||
|
|
|