This paper addresses the challenges posed by the low photon flux and limited atmospheric penetration capabilities of X-ray pulsar signals, which are further exacerbated by technical, economic, and temporal constraints during in-orbit observations. These challenges result in a lack of continuous and high-fidelity data support for ground-based X-ray pulsar navigation (XPNAV) research and validation. To overcome these obstacles, we propose an innovative method for the efficient and rapid simulation of X-ray pulsar signals as received by spacecraft. Building on the physical process that, although the arrival times of the same pulsar signal differ between the spacecraft and the Solar System Barycenter (SSB), these times can be calibrated to the same phase at the SSB using a timing model, we leverage the relationship between the pulsar signal’s arrival rate at the SSB and its arrival phase. Using the Order Statistical Method (OSM), we directly and effectively generate the arrival phase series of pulsar signals received by the spacecraft, which follow a non-homogeneous Poisson process (NHPP). Subsequently, from these phase series, we derive the arrival time series of spacecraft-received pulsar signals, incorporating multiple physical effects based on a relationship model between arrival time series and phases. This method significantly reduces the required simulation time while maintaining high accuracy, thereby greatly enhancing simulation efficiency. Validation against data from the Rossi X-ray Timing Explorer (RXTE) and the Neutron Star Interior Composition Explorer (NICER) shows that the correlation coefficients between the accumulated profiles of simulated and measured data exceed 0.99. Furthermore, our method reduces the simulation time of high-flux pulsar signals by 25.7%, high spin frequency signals by 40.8%, and long observation duration signals by 25.38% compared to the latest fast simulation techniques, highlighting its advantages in accuracy, speed, and versatility for simulating pulsar signals across varying flux rates, spin frequencies, and observation durations. This approach provides continuous and effective data support for XPNAV research, significantly advancing theoretical and practical developments in the field.