Abstract: The numerical implementation of pump behavior's constraint in operation opt imization for pipeline, including applicability predication and solution of pump behavior equation about pump joint-work forms in which each form contains severa l different pump types, is outlined. When the constraint conditions are preset f or every running pump type, such as minimal efficiency and maximal axial po wer and so on, by the method in this paper, the applicability predication can be quite well made for all possible pump joint-work forms which are determined by the number of pumps and pump types installed in pumping/heating stations. The un available forms will be excluded in the analysis progress if they do not satisfy operation requirements or have too great difference in their characteristic amon g join t-work pumps. The calculation cases show that the numerical implementation metho d has a high precision and provides a good reference in selection and decision o f optimal joint-work form in practice.