Based on the Bayesian Inference (BI) method, the Mixture-Model approach is improved to combine all kinematic data, including the coordinative position($\vec{x}$), proper motion ($\vec{\mu}$) and radial velocity(v), to separate the motion of the cluster from field stars, as well as to determine the intrinsic kinematic status and dynamical effects of the cluster, such as the mass segregation, anisotropy etc.. Meanwhile, the membership probability of individual stars are estimated as by product results. This method has been testified by simulation of toy models and also successfully used for well studied open clusters, such as M67 and NGC188. It is expected to largely help the studies of open clusters while combine the coming GAIA data.