Abstract: High intensity laser ablation is a coplex physical and chemical process involving many complicated factors, such as mass removal, phase transition and dynamical boundary conditions. If the ablating target is a composite, the mechanism of laser ablation will be further complicated and it is hard to find a general solution. Only a numerical method may simulate this process perfectly. Based on the theory and the model of high intensity laser ablation, and the method of smoothed particle hydrodynamics (SPH), a meshfree particle method is used in the numerical simulation of high intensity laser ablation of typical materials. The ablation pits and temperature fields of one layer of aluminum and two layers of epoxy/aluminum irradiated by high intensity laser are obtained in this paper. The numerical results show that the ablation will expand in a radial way at the interface of different materials.