Thermosonic wire bonding is a common fabrication process for connecting devices in electronic packaging. However, when the free air ball (FAB) is compressed onto the I/O pad of the chip during bonding procedure, chip cracking may occur if the contact pressure is too large. This study proposes an effective simulation technique that can predict the wire ball geometry after bonding in an accurate range. The contact force obtained in the simulation can be used for possible die cracking behavior evaluation. The simulation in this study used the explicit time integration scheme to deal with the time marching problem, and the second-order precision arbitrary Lagrangian-Eulerian (ALE) algorithm was used to deal with the large deformation of the wire ball during the bonding process. In addition, the equilibrium smoothing algorithm in LS-DYNA can make the contact behavior and geometry of the bonding wire almost the same as the experiment, which can also significantly reduce the distortion of the mesh geometry after remeshing.
