Recently, many people have become interested in seismic stability enhancement and a chain of research and development be proceed for application of nuclear power plant according to increase the frequency and magnitude of earthquake event. Such as seismic isolation system is applied to general structure (architecture, bridge and LNG tank etc.) from ancient times. But the application results is limited for Nuclear power plant. In this paper, we proposed a stability of variable axial load from beyond design basis earthquake in Nuclear power plant. Also, the change of stiffness in isolator from the application of generally design equation is not equal to according to change in axial load compare with the experimental result in variable axial load. Therefore we proposed the empirical formula of design equation from test result of full-scale multi-lead rubber bearing for seismic analysis with real behavior (variable axial load) in the earthquake motion.