Home health care companies provide health care services to patients in their homes. Due to increasing demand, the provision of home health care services requires effective management of operational costs while satisfying both patients and caregivers. In practice, uncertain service times might lead to considerable delays that adversely affect service quality. To this end, this paper proposes a new bi-objective optimization problem to model the routing and scheduling problems under uncertainty in home health care, considering the qualification and workload of caregivers. A mixed-integer linear programming formulation is developed. Motivated by the challenge of computational time, we propose the Adaptive Large Neighborhood Search embedded in an Enhanced Multi-Directional Local Search framework (ALNS-EMDLS). A stochastic ALNS-EMDLS is introduced to handle uncertain service times for patients. Three kinds of metrics for evaluating the Pareto fronts highlight the efficiency of our proposed method. The sensitivity analysis validates the robustness of the proposed model and method. Finally, we apply the method to a real-life case and provide managerial recommendations.