The chapter deals with numerical solution of 2D contact shape optimization problems with Tresca friction and a solution dependent coefficient of friction. The discretized problem is transformed into a nonsmooth minimization problem and handled by the bundle trust method.
Necessary subgradient information is computed using the generalized differential calculus of Mordukchovich. Numerical results of model examples are presented.