In order to modeling seismic isolation system such as lead-rubber bearing (LRB), bilinear model is widely used by many researchers. In general, an actual force-displacement relationship for LRB has a smooth hysteretic shape. So, Bouc-Wen model with smooth hysteretic shape represents more accurately actual hysteretic shape than bilinear model. In this study, seismic responses for seismically isolated nuclear power plant (NPP) with LRB modelled by Bouc-Wen and bilinear models are compared with those of NPP without seismic isolation system. To evaluate effect of earthquake characteristics for seismic responses of NPP isolated by LRB, 5 different site class earthquakes distinguished by Geomatrix 3rd Letter Site Classification and artificially generated earthquakes corresponding to standard design spectrum by Reg. Guide 1.60 are used as input earthquakes. From the seismic response results of seismically isolated NPP, it can be observed that maximum displacements of seismic isolation modelled by Bouc-Wen model are larger than those by bilinear model. Seismic responses of NPP with LRB is significantly reduced than those without LRB. This reduction effect for seismic responses of NPP subjected to Site A (rock) earthquakes is larger than that to Site E (soft soil) earthquakes.