The paper presents the analysis of heavy-loaded low-speed multi-pair spiroid gears under the effect of elastic contact, bending and shearing, and elastoplastic contact interaction of spiroid gearwheel teeth and spiroid worm threads. The main stages of the developed algorithm for analysis of load distribution and plastic strain of flanks with account of elastic and elastoplastic contact conditions, multi-pair contact and macro-roughnesses on flanks are described in short. At elastoplastic contact the zones of plastic contact on flanks of spiroid gearwheel teeth are determined and the value of plastic strain is analyzed. The flanks of the spiroid gearwheel teeth and worm threads are represented in the algorithm as the set of segments (cells) with coordinates of centers determined with account of the analyzed factors that influence the load distribution in spiroid gearing. These factors are the errors of manufacturing and (or) assembly, micro-roughnesses, and deformations of spiroid gear supports. The paper gives the numerical results of research of the real spiroid gears of gearboxes for pipeline valves (PV) at elastic and elastoplastic contact represented as summary tables and diagrams of contact stress distribution along contact areas and of plastic strain, that show the workability of the algorithm. The presented numerical results of the algorithm functioning are well agreed with the experience of operation of spiroid gears and can be applied for analysis of load-carrying capacity of spiroid gears.