According to the high iterative feature of shell and tyre structure, the multi-level substructure technique is adopted to set up the multi-body contact model between tyres and long shell. Applying parametric quadratic programming method, the multi-body contact FEM analysis of the No.2 rotary kiln in Henan Branch of China Aluminum Company is done, and the conclusions are as follows： The equivalent stress of shell at supports changes 5 times, and at other parts changes 4 times; the strength of shell at supports is weak, and at other parts is affluent; the strength of shell at the supports 2 and 3 is the weakest.