The contact relations between rockfill and concrete can be simulated truly by Clough hyperbola constitutive model
and the contact state can also be simulated by augmented Lagrangian multiplier algorithm accurately. Clough hyperbola constitutive model has been applied on jointed element with no thickness and Desai interface element with thickness successfully. Those merits and success experiences are absorbed and combined to be applied on solution of contact friction problem. Sequential quadratic programming method (SQP) is introduced based on augmented Lagrangian multiplier contact algorithm in the same time. Numerical algorithm and implementing steps that Clough hyperbola constitutive model is inducted into augmented Lagrangian multiplier contact algorithm are derived in detail. The computation results of numerical example and the project example show that the stress and deformation state between contact surfaces can be simulated truly with the method.