The numerical model of carbon fiber reinforced polymer (CFRP) laminates with electrically modified interlayers subjected to lightning strike is constructed through finite element simulation, in which both intra-laminar and inter-laminar lightning damages are considered by means of coupled electrical-thermal-pyrolytic analysis method. Then the lightning damage extents including the damage volume and maximum damage depth are investigated. The results reveal that the simulated lightning damages could be qualitatively compared to the experimental counterparts of CFRP laminates with interlayers modified by nickel-coated multi-walled carbon nanotubes (Ni-MWCNTs). With higher electrical conductivity of modified interlayer and more amount of modified interlayers, both damage volume and maximum damage depth are reduced. This work provides an effective guidance to the anti-lightning optimization of CFRP laminates.