It is essential to understand the nonlinear free vibration behaviour of a thin-walled composite structure as they find their applications in harsher environments prone to large deformations. Further, these structures made of fibrous materials, have directional properties for different fibre orientations resulting in a diversified and a complex nonlinear behaviour. In this regard, present work provides a comprehensive mathematical formulation on nonlinear vibration of pre-twisted thin-walled anisotropic box beam. The mathematical model is derived as coupled (flap-lag-extension-torsion) model considering shear deformation and green's strain tensor for large displacements in order to study the influence of nonlinear couplings on the dynamic behaviour. The derived energy equations are presented depending on degree of nonlinearity. These nonlinear equations those derived are nondimensionalized and solved with in the frame work of energy method using classical Ritz approximation. An iterative method is used to solve the nonlinear equations pertaining displacement terms. This nonlinear behaviour is evaluated for two important types of structural tailoring techniques available for thin-walled composite beams namely Circumferentially Asymmetric Stiffness (CAS) and Circumferentially Uniform Stiffness (CUS). It is found that CAS layup experiences significantly severe nonlinear effects compared to CUS due to the presence of 3rd degree nonlinear coupling terms. The variation of nonlinear frequency ratios over different ply angles showing the influence of nonlinearity on fibres orientation is presented for the first time for the two ply lay ups. The influence of degree of nonlinearity on the accuracy of the nonlinear frequencies is evaluated and presented. Moreover, the impact of nonlinearity on the pre-twisted beam is presented for different pre-twist angles.