We need to simulate the clinch joining of aluminum and Carbon Fiber Reinforced Polymers. We would like to predict the damaged zone and perform the simulation using Hou criteria. We have already seen that a number of authors that developed VUMAT Fortran subroutine and that perfomred simulation of drilling process using such a model within ABAQUS environment. Thus we are looking for possible collaboration in this field. All contributes are welcome.