Based on simulation of molecular models of polymeric fluids，a new micro-macro method for numerical calculation of viscoelastic flow is presented in this paper. This method coupled the mass and momentum conservation equations at the macroscopic level，with a stochastic differential equation which models the evolution of the polymer configures at the microscopic level. Accordingly，the velocity filed was determined by solving the conservation equation with a meshless method，and the stress was computed from the molecular configuration rather than from closed-form constitutive equations. Thus this method bypassed the need for a rheological constitutive equation to describe the polymeric fluid. In order to validate this method and to demonstrate its robustness，the start-up planar Couette flow was studied for Hookean，FENE and FENE-P dumbbell models and the lid-driven cavity flow was studied for Hookean dumbbell models.