This work presents a refinement of the finite strip method (FSM) based on the Carrera unified formulation (CUF) and is applied to free vibration analysis of variable stiffness composite laminated (VSCL) plates. VSCL plates considered here are composed of layers with curvilinear fibers in which their orientation angle changes linearly with respect to one of the in‐plane coordinates. The FSM permits the plate to be divided into some finite strips connected along the so‐called nodal lines. The conventional finite strip method involves an approximation of the solution of a plate problem by using continuously harmonic series in the longitudinal direction along the nodal lines and simple polynomial shape functions (the usual beam shape functions) in the transverse direction of the plate. However, refined finite strip model that presented in this paper allows one to express the unknown variables as a set of thickness functions that only depend on the thickness coordinate and the corresponding variable that depends on the in‐plane coordinates which involve continuously harmonic series and polynomial shape functions. Thanks to this new finite strip model, the shape functions used in the transverse direction can be chosen as 1D bar shape functions, and also, all three components of the displacement field are considered in each nodal line. Moreover, based on this approach, three‐dimensional displacement field is approximated in a compact form as a generic N‐order expansion. Therefore, explicit expressions of fundamental nuclei of finite strip matrices are obtained in a compact form and presented here. The results obtained by the present formulation for VSCL plates are compared with those obtained from published data. The results show that this kind of strip is efficient and useful. Further, several numerical examples are presented, and the effect of the order of expansion, the number of strips, curvilinear fiber angle, and various boundary conditions are examined.