This paper presents a contractiveness preserving model reduction method for general Q-stable linear fractional transformations (LFTs) with guaranteed error bound. The approach is based on balancing positive definite solutions to dual linear matrix inequalities for contractiveness. A numerical example is given to demonstrate the performance of the proposed model reduction method.