The changing process of internal defects of steel fiber reinforced concrete,such as initiation and expansion,needs to be studied from a mesoscopic perspective.In order to establish a more realistic meso-scopic random distribution model of steel fiber reinforced concrete,an efficient and reliable mesoscopic modeling method of steel fiber reinforced concrete is proposed.The secondary development of the ABAQUS program is carried out in Python language,which avoids the tedious problem of the information transfer process between different software in modeling.A new interference judgment method for polygon aggregate and steel fiber is proposed-point line method,and the equation of taking points on the boundary of arbi-trary elliptical aggregate is derived,which improves the efficiency of aggregate delivery.Based on the above method,the failure process of steel fiber concrete specimen under uniaxial compression is simulated and compared with the test results.The results show that the delivery ratio of second graded polygon aggregate is 70%,the delivery ratio of second graded elliptical aggregate is 69%in 2D mesoscopic model and the intersection of steel fiber in modeling process is eliminated.The failure patterns and stress-strain curves obtained by numerical simulation are in good agreement with the experimental results.The feasibility of the proposed method is further verified.