To support the airworthiness certification of aircraft engines regarding CCAR 33.63 provisions, a three-dimensional contact dynamics-based computational model for fan blade-flexible casing rubbing vibration was established. Modal reduction was employed to reduce the size of the computational model, and a cubic B-spline surface fitting technique was used to enhance the discretization accuracy of the casing inner surface. A bilinear elastic-plastic model was adopted to describe the rubbing mechanical properties of the abradable coating, and contact forces between blades and casings were calculated based on coating wear depth, blade thickness, and element shape functions. Using this model, the rubbing response characteristics between fan blades and flexible casings were investigated. The results showed that blade damping is a key factor affecting rubbing vibration response: for blades with lower damping, in addition to coupling vibrations between bladecasing single nodal diameter modes during rubbing process, there also exist blade-casing coupled vibrations caused by traveling waves associated with multi-nodal diameter modes of casings. The variations in contact strength during rubbing process leading to changes in dynamic characteristics of components caused the above interactions. The dynamic model established in this study provides a new approach for identifying and evaluating blade-casing rubbing coupled vibrations.