An improved model for predicting stiffness of single-lap composites bolted joints using Matlab/Simulink