The mechanical response results for variable stiffness composite laminates were obtained by employing the automated modeling and finite element analysis, facilitating the investigation into the compression properties of variable stiffness composite laminates with a hole. The findings demonstrate the feasibility of utilizing programming languages to automatically model curvilinear fiber placement within these laminates, which enables independent design of the fiber orientation angles and the thickness of each element in the finite element model. Additionally, finite element analysis provides displacement result under compressive loads for variable stiffness composite laminate with a hole, in which the uneven Y-direction deformation is attributed to the variations in effective Poisson's ratio caused by the different fiber orientations. Compared to the stress result of the constant stiffness laminate, the results for both force resultants distribution and stress distribution of the variable stiffness laminate reveal that curvilinear fibers effectively redistribute loads and reduce the shear stress concentration resulting from the hole.