Abstract: The main four types of damage mode (matrix-cracking, fiber-matrix shear, delamination and fiber breakage) were comprehensively considered. A parametric modeling program was developed based on the finite element analysis software ANSYS. The progressive damage process of two kinds of layer composite laminates with a circular hole under tensile loading were analyzed by this 3D progressive damage model and finite element method. The results indicate that the progressive damage process and ultimate tensile strength of the composite laminates with arbitrary material system and layer parameters at tensile load could be predicted by this program. The prediction accuracy of this method was satisfactory by comparison with experimental data in literatures, which was proved the analysis method and program were valid.