Health and convenience are two indispensable indicators of the society promotion. Nowadays, to improve community health levels, the comfort of patients and those in need of health services has received much attention. Providing Home Health Care (HHC) services is one of the critical issues of health care to improve the patient convenience. However, manual nurse planning which is still performed in many HHC institutes results in a waste of time, cost, and ultimately lower efficiency. In this research, a multi-objective mixed-integer model is presented for home health care planning so that in addition to focusing on the financial goals of the institution, other objectives that can help increase productivity and quality of services are highlighted. Therefore, four different objectives of the total cost, environmental emission, workload balance, and service quality are addressed. Taking into account medical staff with different service levels, and the preferences of patients in selecting a service level, as well as different vehicle types, are other aspects discussed in this model. The epsilon-constraint method is implemented in CPLEX to solve small-size instances. Moreover, a Multi-Objective Variable Neighborhood Search (MOVNS) composed of nine local neighborhood moves is developed to solve the practical-size instances. The results of the MOVNS are compared with the epsilon-constraint method, and the strengths and weaknesses of the proposed algorithm are demonstrated by a comprehensive sensitivity analysis. To show the applicability of the algorithm, a real example is designed based on a case study, and the results of the algorithm over real data are evaluated.