It must be the way I defined the element. I used a more complicated mesh from an actual meshing algorithm and got a non-singular matrix that yielded reasonable displacement values.
Thanks for recomending those books, had one of them on my wishlist already.
So, in other words, the above methods I used should work if I have my stiffness matrix properly constrained?
Right now I'm running the program on a single Q8 element. Each node has two in plane translational Dof's. I've constrained the bottom three nodes in both directions, and applied a load...
I've written a Matlab program that creates a global stiffness matrix for a group of Q8 elements.
I have a force vector, now I'm trying to apply boundary conditions and solve for the displacements. The problem I've run into, is whether I eliminate rows and columns correlated to Dof's that equal...
