In this research, the nonlinear static, buckling and vibration analysis of viscoelastic micro-composite beam reinforced by various distributions of boron nitrid nanotube (BNNT) with initial geometrical imperfection by modified strain gradient theory (MSGT) using finite element method (FEM) are presented. The various distributions of BNNT are considered as UD, FG-V and FG-X and also, the extended rule of mixture is used to estimate the properties of micro-composite beam. The components of stress are dependent to mechanical, electrical and thermal terms and calculated using piezoelasticity theory. Then, the kinematic equations of micro-composite beam using the displacement fields are obtained. The governing equations of motion are derived using energy method and Hamilton\'s principle based on MSGT. Then, using FEM, these equations are solved. Finally the effects of different parameters such as initial geometrical imperfection, various distributions of nanotube, damping coefficient, piezoelectric constant, slenderness ratio, Winkler spring constant, Pasternak shear constant, various boundary conditions and three material length scale parameters on the behavior of nonlinear static, buckling and vibration of micro-composite beam are investigated. The results indicate that with an increase in the geometrical imperfection parameter, the stiffness of micro-composite beam increases and thus the non-dimensional nonlinear frequency of the micro structure reduces gradually.