Ordering of Gauss points

[K.Li]

Dear FEAP users:

It seems the ordering (numbering) of Gauss points in a 8-node brick element doesn't follow element node numbering in FEAP. The first Gauss point seems to be located near the 5th node of the element. This is different from other FE software. Could anyone please tell me how this ordering is determined or point me to the routine where this was defined? Can this ordering be changed for a specific element? 


Thanks a lot!

[Prof. R.L. Taylor]

The quadrature order is defined in int3d.f (/program directory).

For the 2x2x2 quadrature the points would follow node numbering for the 8-node brick by setting:
        s(3,i    ) = -g
        s(3,i+4) =  g

[K.Li]

Thank you for your help, Prof Taylor!