With the development of factional calculus theory and applications in different fields in recent years, the rational approximation problem of fractional calculus operator has become a hot spot of research. In the early 1950s and 1960s, Carlson and Halijak proposed regular Newton iterating method to implement rational approximation of the one-nth calculus operator. Carlson regular Newton iterating method has a great sense of innovation for the rational approximation of fractional calculus operator, however, it has been used only for certain calculus operators. The aim of this paper is to achieve rational approximation of arbitrary order fractional calculus operator. The realization is achieved via the generalization of Carlson regular Newton iterating method. To construct a rational function sequence which is convergent to irrational fractional calculus operator function, the rational approximation problem of fractional calculus operator is transformed into the algebra iterating solution of arithmetic root of binomial equation. To speed up the convergence, the pre-distortion function is introduced. And the Newton iterating formula is used to solve arithmetic root. Then the approximated rational impedance function of arbitrary order fractional calculus operator is obtained. For nine different operational orders with n changing from 2 to 5, the impedance functions are calculated respectively through choosing eight different initial impedances for a certain operational order. Considering fractional order operation characteristics of the impedance function and the physical realization of network synthesis, the impedance function should satisfy these basic properties simultaneously: computational rationality, positive reality principle and operational validity. In other words, there exists only rational computation of operational variable s in the expression of impedance function. All the zeros and poles of impedance function are located on the negative real axis of s complex plane or the left-half plane of s complex plane in conjugate pairs. The frequency-domain characteristics of impedance function approximate to those of ideal fractional calculus operator over a certain frequency range. Given suitable initial impedance and for an arbitrary operational order, it is proved that the impedance function could meet all properties above by studying the zero-pole distribution and analyzing frequency-domain characteristics of the impedance function. Therefore, the impedance function could take on operational performance of the ideal fractional calculus operator and achieve the physical realization. It is of great effectiveness in the generalization of this kind of method in both theory and experiment. The results educed in this paper are the basis for further theoretic research and engineering application in constructing the arbitrary order fractional circuits and systems.