Plot3D[(x^2+y^2)/2+Log[Cosh[Sqrt[x^2+y^2]-2]]+Log[Cosh[Sqrt[(x-1)^2+y^2]-1]],{x,-2,2},{y,-2,2}]