We study the self‐entanglement of a single linear polymer chain with N monomers confined to a cubic box (L × L × L) using the bond‐fluctuation lattice model and primitive path analysis. We probe chains with N between 30 and 750 and vary the degree of confinement L/Rg0 between 0.4 and 12, where Rg0 is the radius of gyration of an unconfined polymer. We find that the conformational properties Rg/Rg0 and Lp/Rg0, where Lp is the average primitive path length, collapse onto a single master curve as a function of the degree of confinement. In the strongly confined regime, L/Rg0 < 1, we find that Rg/Rg0 ∼ (L/Rg0)0.8 and (Lp/Rg0) ∼ (L/Rg0)−2. We verify that the simulation methodology used is quantitatively consistent with experimental data, and the Colby‐Rubinstein entanglement model for unconfined concentrated polymer solutions. The most significant difference between unconfined and confined systems is the variation of Lp with monomer density ϕ; Lp ∼ ϕ5/9, in the former, and Lp ∼ ϕ2/3, in the latter. © 2014 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2014, 52, 1283–1290