Understanding and characterizing semi-crystalline models with crystalline and amorphous segments is crucial for industrial applications. A coarse-grained molecular dynamics (CGMD) simulations study probed the crystal network formation in high-density polyethylene (HDPE) from melt, and shed light on tensile properties for microstructure analysis. Modified Paul–Yoon–Smith (PYS/R) forcefield parameters are used to compute the interatomic forces among the PE chains. The isothermal crystallization at 300 K and 1 atm predicts the multi-nucleus crystal growth; moreover, the lamellar crystal stems and amorphous region are alternatively oriented. A one-dimensional density distribution along the alternative lamellar stems further confirms the ordering of the lamellar-stack orientation. Using this plastic model preparation approach, the semi-crystalline model density (ρcr) of ca. 0.913 g·cm−3 and amorphous model density (ρam) of ca. 0.856 g·cm−3 are obtained. Furthermore, the ratio of ρcr/ρam ≈ 1.06 is in good agreement with computational (≈1.096) and experimental (≈1.14) data, ensuring the reliability of the simulations. The degree of crystallinity (χc) of the model is ca. 52% at 300 K. Nevertheless, there is a gradual increase in crystallinity over the specified time, indicating the alignment of the lamellar stems during crystallization. The characteristic stress–strain curve mimicking tensile tests along the z-axis orientation exhibits a reversible sharp elastic regime, tensile strength at yield ca. 100 MPa, and a non-reversible tensile strength at break of 350%. The cavitation mechanism embraces the alignment of lamellar stems along the deformation axis. The study highlights an explanatory model of crystal network formation for the PE model using a PYS/R forcefield, and it produces a microstructure with ordered lamellar and amorphous segments with robust mechanical properties, which aids in predicting the microstructure–mechanical property relationships in plastics under applied forces.