Thanks, Xerf.
You are right, it's strange. I compared stiffness matrix, because i want to do some manipulations in global level. So far, I've got exactly cps4, cpe8 and cps8.
Sorry, I did not make myself clear, all constitutive matrices use 3x3. Your formula are essentially the same with mine...