Mechanical integrity of the dielectric stack is challenged by the trend towards porous, lower dielectric constant interlayer dielectric (ILD) materials. As a result, fracture in the ILD stacks caused either by assembly process or by the dicing process is an important reliability consideration. In general, there is a need to either assess the propensity of the structure to fracture under assembly conditions, or to design crackarrest features that prevent propagation of cracks into active areas.In the case of wire bonded packages, the reliability concern associated with the fracture of Ultra Low-k (ULK) dielectrics while bonding over the active circuits (BOAC) is a significant challenge due to the impact load and the high ultrasonic energy transmitted to the ILD stack. In this paper, a multilevel modeling procedure is presented to assess the risk of fracture in ILD stacks during wire bonding process. First, a nonlinear, dynamic finite element model is developed to simulate the process steps -impact stage and last cycle of ultrasonic vibration and study the mechanical response of the ball, pad, and the underlying ULK under pad during copper wire bonding. Further, a simulation framework based on enriched isogeometric approximations is presented to compute damage accumulation in the ULK stacks using a cohesive damage description. The simulation framework is employed to develop insights on the potential crack initiation sites within the ILD stack and to evaluate the risk of fracture during each process step. KEY WORDS: wire bonding, ULK, damage, fracture, reliability NOMENCLATURE C bi-variate parametric geometry D damage E elastic modulus, N/m 2 G energy release rate, J/m 2 N NURBS basis functions V tri-variate parametric geometry d distance, m f behavioral field, unit w weight field s parameter of bi-variate geometry, C t parameter of bi-variate geometry, C û field unknowns at the control points, unit Greek symbols Ω domain Γ boundary ρ density, kg/m 3 ν Poisson's ratio σ stress, N/m 2 ϒ fracture toughness, J/m 2 ψ enriching function δ separation distance, m ζparameter of tri-variate geometry, V η parameter of tri-variate geometry, V ξ parameter of tri-variate geometry, V
