$ For CST=1, the integration points are arranged circumferentially (located at a radius of 
$ 0.707 times the outside radius of the cross-section.  Thus for 4x4 integration, the IP 
$ are 22.5 degrees apart.  For an illustration, use  History > Int pt > Etype: beams > Axial stress   
$ to plot the integration point stresses in the attached example of a cantilever beam in bending.
*KEYWORD
*control_implicit_general
1,.1
$ use damping for explicit run
$*damping_global
$,30.
*title
9 vs. 16 IP in circular sections
$ beam IP output in numerical order (see p. 3.14 in 940 structured user's manual)
*interface_springback_lsdyna
5
*set_part
5
1
*CONTROL_TERMINATION
 0.9000000         0 0.0000000         0 0.0000000
*DATABASE_BINARY_D3PLOT
 0.0100000                   0
*DATABASE_EXTENT_BINARY
         0         0         0         0         0         0         0         0
         0         0        16         0         0         0
*DATABASE_NODOUT
 2.00000-3
*DATABASE_HISTORY_BEAM
1,11
*DATABASE_ELOUT
 2.00000-3
*DATABASE_SPCFORC
 2.00000-3
*DATABASE_BNDOUT
 2.00000-3
*MAT_ELASTIC
         1 3.00000-4 1.00000+7 0.2000000
$*mat_plastic_kinematic
$$1,2e-4,1e7,.2,  20,0,1
$1,2e-4,1e7,.2,  15,1e4,1
$
*SECTION_BEAM
1,1,,3,1
1,1,0,0
$1,1,.5,.5
*SECTION_BEAM
2,1,,4,1
1,1,0,0
$1,1,.5,.5
*PART
material type # 1                                                               
         1         1         1         0         0         0         0         0
*PART
material type # 1                                                               
         2         2         1         0         0         0         0         0
*NODE
       1 1.800000000E+01 0.000000000E+00 0.000000000E+00       0       0
       2 3.600000000E+01 0.000000000E+00 0.000000000E+00       0       0
       3 5.400000000E+01 0.000000000E+00 0.000000000E+00       0       0
       4 7.200000000E+01 0.000000000E+00 0.000000000E+00       0       0
       5 0.000000000E+00 1.000000000E+00 0.000000000E+00       0       0
      11 1.800000000E+01 0.000000000E+00 9.000000000E+00       0       0
      12 3.600000000E+01 0.000000000E+00 9.000000000E+00       0       0
      13 5.400000000E+01 0.000000000E+00 9.000000000E+00       0       0
      14 7.200000000E+01 0.000000000E+00 9.000000000E+00       0       0
      15 0.000000000E+00 1.000000000E+00 9.000000000E+00       0       0
*ELEMENT_BEAM
       1       1       1       2       5
       2       1       2       3       5
       3       1       3       4       5
      11       2      11      12      15
      12       2      12      13      15
      13       2      13      14      15
*BOUNDARY_SPC_NODE
         1         0         1         1         1         1         1         1
        11         0         1         1         1         1         1         1
*DEFINE_CURVE
         1         0 0.0000000  2. 
      0.00000000E+00      0.00000000E+00
      5.00000001E-01      1.0
      1.00000000E+00      1.0
*LOAD_NODE_POINT
         4         2         1 
        14         2         1 
*END