Poly(ϵ-caprolactone) (PCL) diol, with good biodegradation and biocompatibility, is one of the widely used soft segments (SSs) in composing bio-polyester-urethanes (Bio-PUs), which show great potential in both biomedical and tissue engineering applications. Properties of Bio-PUs are tunable by combining SS monomers with different molecular weights, structures, modifications, and ratio of components. Although numbers of research have reported many Bio-PUs properties, few studies have been done at the molecular scale. In this study, we use molecular dynamic (MD) simulation to construct atomistic models for two commonly used PCL diol SSs with different molecular weights 1247.58 Da and 1932.42 Da. We compare the simulation results by using two widely used classical force fields for organic molecules: Consistent Valence Force Field (CVFF) and CHARMM General Force Field (CGenFF), and discuss the validity and accuracy. Melt density, volume, polymer conformations, transition temperature, and mechanical properties of PCL diols are calculated and compared with experiments. Our results show that both force fields provide accurate predictions on the properties of PCL diol system at the molecular scale and could help the design of future Bio-PUs.