matlab有限元轴孔过盈接触分析如何进行迭代？

如何利用matlab编程进行轴孔过盈有限元分析，利用罚函数法。有会的大佬请联系我，万分感谢