A family of n-sub-step composite time integration methods, which employs the trapezoidal rule in the first n-1 sub-steps and a general formula in the last one, is discussed in this paper. A universal approach to optimize the parameters is provided for any cases of nge 2, and two optimal sub-families of the method are given for different purposes. From linear analysis, the first sub-family can achieve nth-order accuracy and unconditional stability with controllable algorithmic dissipation, so it is recommended for high-accuracy purposes. The second sub-family has second-order accuracy, unconditional stability with controllable algorithmic dissipation, and it is designed for heuristic energy-conserving purposes, by preserving as much low-frequency content as possible. Finally, some illustrative examples are solved to check the performance in linear and nonlinear systems.