We study the heat conduction of two nonlinear lattices joined by a weak harmonic link. When the system reaches a steady state, the heat conduction of the system is decided by the tunneling heat flow through the weak link. We present an analytical analysis by the combination of the self-consistent phonon theory and the heat tunneling transport formalism, and then the tunneling heat flow can be obtained. Moreover, the nonequilibrium molecular dynamics simulations are performed and the simulations results are consistent with the analytical predictions.