Dear Omid and Sepehr
I think a non-zero entry in the stiffness matrix doesn't necessarily imply dependence of ith DOF on the jth DOF. This is the case also for standard FEM (without any enrichments). For example for a quadrilateral 2D element the x component of displacement has interaction with y component at each node, but as we know, two independent values of translation in x and y directions can