% Assemble global stiffness matrix K = [D11, D12, D16; D12, D22, D26; D16, D26, D66];