Disassembly enterprises face significant challenges in managing increasingly intensive workloads with a workforce alone. To address this, current study focuses on balancing a human-robot collaborative disassembly line. For the task operator allocation problem, a mapping constraint mechanism between the disassembly tasks and the task operators is established based on the remanufacture demand and hazard characteristics of components. A mixed-integer programming model is formulated, encompassing four optimization objectives: the number of workstations, the task operator idle balancing index, the demand index, and the hazard index. A novel modified teaching and learning optimization approach is proposed for the addressed human-robot collaborative disassembly line balancing problem. A two-layer encoding strategy is designed based on the solving characteristics of the problem and an embedded perturbation strategy based on randomly transforming the task operator is introduced to enhance the optimization performance of the proposed method. The validity of the mixed-integer programming model and the efficacy of the proposed algorithm are demonstrated through two small-scale human-robot collaborative disassembly case studies. The proposed algorithm is then applied to two real-life cases of human-robot collaborative disassembly lines with different scales. The results of the proposed optimization method are compared with other advanced algorithms, demonstrating its superior performance in both single-objective and multi-objective optimization based on multiple evaluation indicators.