The regions of stability of two collinear quasiclassical phases within the zerotemperature quantum phase diagram of the spin-1 2 J 1 -J 2 -J ⊥ 1 model on an AA-stacked bilayer honeycomb lattice are investigated using the coupled cluster method (CCM). The model comprises two monolayers in each of which the spins, residing on honeycomb-lattice sites, interact via both nearestneighbor (NN) and frustrating next-nearest-neighbor isotropic antiferromagnetic (AFM) Heisenberg exchange iteractions, with respective strengths J 1 > 0 and J 2 ≡ κJ 1 > 0. The two layers are coupled via a comparable Heisenberg exchange interaction between NN interlayer pairs, with a strength J ⊥ 1 ≡ δJ 1 . The complete phase boundaries of two quasiclassical collinear AFM phases, namely the Néel and Néel-II phases on each monolayer, with the two layers coupled so that NN spins between them are antiparallel, are calculated in the κδ half-plane with κ > 0. Whereas on each monolayer in the Néel state all NN pairs of spins are antiparallel, in the Néel-II state NN pairs of spins on zigzag chains along one of the three equivalent honeycomb-lattice directions are antiparallel, while NN interchain spins are parallel. We calculate directly in the thermodynamic (infinite-lattice) limit both the magnetic order parameter M and the excitation energy ∆ from the s z T = 0 ground state to the lowest-lying |s z T | = 1 excited state (where s z T is the total z component of spin for the system as a whole, and where the collinear ordering lies along the z direction), for both quasiclassical states used (separately) as the CCM model state, on top of which the multispin quantum correlations are then calculated to high orders (n ≤ 10) in a systematic series of approximations involving n-spin clusters. The sole approximation made is then to extrapolate the sequences of nth-order results for M and ∆ to the exact limit, n → ∞.