The vapor-liquid equilibrium(VLE)for the system acetic acid,acetaldehyde and crotonaldehyde was investigated by combining the Gibbs ensemble Monte Carlo(GEMC)method and activity coefficient models prediction.The VLE of above pure components and their binary systems were firstly studied using GEMC method by Towhee software.Based on the strong association characteristics of acetic acid studied by molecular dynamic method using Material Studio software,the revised TraPPE-UA-D force field was developed with quantum chemistry methods for the VLE calculation of acetic acid system.The original and mixed TraPPE-UA force field was introduced for the VLE calculation of acetaldehyde and crotonaldehyde system,respectively.The results showed that the predictions agreed fairly well with the average relative deviations of the saturated densities smaller than 4.60%.Moreover,binary VLE diagram was calculated using above force fields.The simulation results of the mole composition were in good agreement with the experimental values,which verified the high simulation accuracy of the adopted force fields.On this basis,the VLE of the ternary system were predicted with the NRTL(Non-random Two-liquid)and NRTL-HOC(Hayden-O'Connell)models,and binary interaction parameters of the models were obtained by fitting the VLE simulation results of above binary systems.The simulation results were in accordance with the predictions obtained by the parameters from the Aspen data source.This work also provided a new and feasible way to obtain the VLE values of mixtures.