In the damping design of blisks with hard coating, changing coating thickness is one of the frequently used ways for achieving better damping performance. Nevertheless, for the vibration analysis with variable coating thickness by using commercial finite element software, it is necessary to regenerate the finite element mesh, which results in inconvenience in computation. In this paper, a new laminated element is developed to solve this problem, which needs solely the node coordinates corresponding to the maximum coating thickness. The additional virtual layers to simulate the variations of coating thickness are introduced for the presented element to avoid regenerating new node coordinates. By comparing the analysis results obtained by SOLID185 in ANSYS software, experiment and the developed element, this new element is verified that it remains effective for variable coating thickness parameter input. Furthermore, the effects of the coating thickness on the natural frequency and the forced response of an actual blisk are analysed and due to using the same node coordinates, the computational efficiency is improved significantly. The analysis results indicate that with the increase of coating thickness, the resonance response is reduced significantly and even the thinner coating could also achieve better damping performance.