The evaluation of post-earthquake functional recovery of buildings is pivotal for seismic resilience assessment and underpins resilience-oriented seismic design. This study proposes a generalized method that can predict the potential functional recovery trajectory of any occupancy type of building. This method distinctly models delay time and repair time. Delay time accounts for the duration required for building mobilization activities and the recovery of utility systems to full functionality. The triggering conditions for mobilization activities are clarified by component damage severity and the repair cost ratio of buildings. The repair time is modeled by embedding the critical path method into an Activity-on-Edge network, incorporating the procedure importance coefficient under limited floor space to optimize labor allocation, thus achieving dynamic coupling of repair processes and labor resource scheduling. Utilizing Monte Carlo simulation as the underlying mechanism, the method enables probabilistic modeling of the functional recovery trajectory and subsequent calculation of the seismic resilience index. A case study of a teaching building demonstrates the method’s efficacy in pinpointing critical recovery bottlenecks, optimizing labor allocation, and delivering recovery trajectory simulations. It exhibits robust adaptability and scalability, providing an effective analytical tool for seismic resilience assessment, enhancement, and design across diverse occupancy types of buildings.