Abstract: Based on the measured data or numerical simulate results of underground mining of coal mine,taking the "beam model hypothesis" of stope roof as the research background,using the functions of curve fitting,polynomial interpolation,differential and integral calculus and solve equation of MATLAB software,the curve function of the discrete data of stope roof beam load obtained by actual measurement or numerical calculation are obtained,the distributive curve of shear and bending moment on the beam is obtained by integral operation,the shear value and bending moment value are acquired,besides that,the maximum value of shear and bending moment are calculated by using the roof value calculation method,and the location of the maximum value of shear and bending moment are also determined.The goal of the quantitative analysis of complex mechanical problems of stope roof under the nonlinear distributed loadings is achieved,which provide a more accurate method to predict roof beam damage.

Key words: MATLAB,Curve fitting, Quantitative analysis, Stope roof, Complex stress