Hello everyone,
From the position of the integration point on the midplane, I guess Abaqus use the newton-cotes rather than the Gauss integration scheme , but the output stiffness matrix has coupling terms between the degrees of freedom.
What integral scheme does abaqus use to numerically calculate the stiffness matrix of COH3D8 elements.