function Ke = element_stiffness(xy, ABD, As) % xy: 4x2 matrix of element nodal coordinates (x,y) % ABD: 6x6 matrix [A B; B D] % As: 2x2 shear stiffness matrix % Returns 20x20 element stiffness matrix
A = zeros( ); B = zeros( ); D = zeros(
