In this talk, we will introduce a hybrid iterative method (Int-deep) for numerically solving the elliptic variational inequality (EVI) of the second kind, through combining the regularized semi-smooth Newton method and the deep learning method. The convergence rate analysis and numerical examples on contact problems show this algorithm converges rapidly and is efficient for solving EVIs. This work is in collaboration with Jianguo Huang and Yujian Cao.