Bivariate interval censored data arises in many applications. However, both theoreticaland computational investigations for this type of data are limited because of thecomplicated structure of bivariate censoring. In this paper, we propose a two-stage splinebasedsieve estimator for the association between two event times with bivariate case 2interval censored data. A smooth and explicit estimator for the joint distribution functionis also available. The proposed estimators are shown to be asymptotically consistent andcomputationally efficient. We demonstrate the finite sample performances of the splinebasedsieve estimators using both simulation studies and real data analysis from an AIDSclinical trial study.