We quantitatively investigated the statistical behaviors of semiflexible polymer chains, which are simultaneously subjected to force stretching and rectangular tube confinement. Based on the wormlike chain model and Odijk deflection theory, we derived a new deflection length, by which new compact formulas are obtained for the confinement free energy and force-confinement-extension relation. These newly derived formulas have been justified by numerical solutions of an eigenvalue problem associated with the Fokker-Planck governing equation and extensive Brownian dynamics simulations based on the so-called Generalized Bead-Rod (GBR) model. We found that, comparing to the classical deflection theory, these new formulas are valid for a much extended range of the confinement-size /persistence-length ratio, and have no adjustable fitting parameters for sufficient long semiflexible chains in the whole deflection regime.