We demonstrate a complete procedure for simulations of electric field poled polymeric nonlinear optical systems with the purpose to evaluate the macroscopic electro-optic coefficients. The simulations cover the electric field poling effects on the chromophore order at the liquid state, the cooling procedure from the liquid to the solid state in the presence of the poling field, and the back-relaxation of the system after the removal of the field. We use Disperse Red chromophore molecules doped in a poly(methyl methacrylate) matrix for a numerical demonstration of the total procedure. On the basis of the simulation results, the polymer mobility and the static properties of the dopant chromophores are derived. In the liquid state, the chromophore molecules are closer to the side chains than to the backbones of the polymer matrix, and after the simulated annealing, the polymer matrix tends to be closely packed, leading to a significant change in the polymer structure around the chromophore molecules. Besides predicting the absolute macroscopic electro-optic coefficient values, the results are used to derive the microscopic origin of these values in terms of geometric and electronic structure, loading, poling, and back-relaxation effects, thereby aiding to establish design principles for optimum guest-host configurations.