Two kinds of multidimensional IFR distribution are defined by using a partial order in Rn+, which is derived from a non-negative, strictly increasing function in Rn+. Some closure properties under operations and an application to a shock model are discussed.