A formula is developed to estimate the total mass loss of projectile, based on the assumptions that the peeling of molten surface layer in projectile nose is the primary cause of mass loss, and the frictional heat is totally absorbed by the projectile. Extrapolating this formula to predict the mass loss of local area of projectile, the receding displacement on projectile surface is obtained, which is vertical to the symmetry axis of projectile. Thereby, a finite difference method model is constructed to simulate the variation of projectile shape. The shape of residual projectile, depth of penetration of projectile and its mass loss obtained by calculation are found in good consistency with respective experimental data.The high-speed projectile penetrating into concrete has a significant mass loss, 1-4 the cause of which is commonly thought as the peeling of molten surface layer in the projectile nose. The heat is totally supplied by the frictional heat between target and projectile. 5,6 However, the dynamic friction between target and projectile is a much complex problem. 7 Despite the efforts made, the study of dynamic friction is still in a pre-mature stage.On the other hand, the mass loss of projectile occurs mainly on projectile nose surface and thus induces projectile nose blunting. The blunter the projectile nose, the worse the performance of projectile. Therefore, it is necessary to obtain the global or even local shape variation of projectile nose. Generally, the variation of projectile shape is simulated by numerical methods, such as finite element method (FEM), 8-10 finite difference method (FDM), 11-13 etc. However, the FEM model for shape variation of projectile usually demands insertion of algorithm of mass abrasion into certain commercial FEM software, 8-10 which is a much complicated procedure. The FDM model is simpler, but it still has its own shortcomings, e.g., the algorithm developed in Ref. 11 is too ambiguous and complicated to be implemented, and Wen et al. confronted singularity in the projectile tip when the impact velocity was high. 13 In the present manuscript, based on Refs. 5 and 10, the frictional work is derived from the energy conservation law, and it is found proportional to the work done by pressure. Since the projectile mass loss occurs only within a thin layer on the projectile nose surface, 2,3,5,6 it is reasonable to assume that the projectile is rigid except its lost mass portion. The displacement on projectile surface is thus obtained during a small time interval, and an FDM model is developed accordingly to a) Corresponding author. Email: chenxiaoweintu@yahoo.com. simulate the variation of projectile shape.Silling and Forrestal found that the mass loss of projectile was proportional to its initial kinetic energy by graphical analysis. 10 He and Chen applied this relationship at higher impact velocities. 14 Therefore, the increment of mass could be expressed aswhere K, v, m are the instant kinetic energy, penetration velocity and mass of projectile, respectively. Since the n...