In this paper an analytical approach for determining load distribution in single-column multi-bolt composite joints by considering elastic nonlinear behavior of the composite materials is presented. Load distribution was calculated by writing the governing equations of the motion. This closed form solution is an integration of spring-based models with nonlinear behavior of composite plate materials. Developed method is capable to calculate taken load by each bolt and its displacement by simultaneous solving of governing equilibrium equations of the system. This manuscript specifically focused on the influence of composite material nonlinearity on the load distribution of single bolted composite joints. For this purpose, load changes versus displacement are plotted by taking into account both the linear and nonlinear material behavior. The achieved results via suggested solution revealed that displacements were increased upto 2.5-5 percent in comparison with the results of linear method available in the literature. In addition, due to the manufacturing tolerances, bolt–hole clearances can vary within allowable limits and fits. Therefore the effect of bolt hole-clearance on the composite joints with linear and nonlinear material properties was also investigated.