% B matrix for CST B = zeros(3, 6); for i = 1:3 j = mod(i,3)+1; k = mod(i+1,3)+1; B(1, 2*i-1) = (y(j)-y(k)) / (2*area); B(2, 2*i) = (x(k)-x(j)) / (2*area); B(3, 2*i-1) = (x(k)-x(j)) / (2*area); B(3, 2*i) = (y(j)-y(k)) / (2*area); end
% --- Assembly --- n_dof = size(nodes,1)*2; K = zeros(n_dof); F = F_applied; matlab codes for finite element analysis m files
% Elements (triangle connectivity: node1, node2, node3) elements = [1, 2, 3; 1, 3, 4]; % B matrix for CST B = zeros(3,
% --- Assembly --- K_global = zeros(n_dof); F_global = zeros(n_dof, 1); k = mod(i+1