This paper presents a thin plate formulation with nodal integration for bending analysis using three-node triangular cells and linear interpolation functions. The formulation was based on the classic thin plate theory, in which only deflection field was required and dealt with as the field variables. They were assumed to be piecewisely linear and expressed using a set of three-node triangular cells. Based on each node, the integration domain has been further derived, where the curvature in the domain was computed using a gradient smoothing technique (GST). As a result, the curvature in each integration domain is a constant whereby the deflection is compatible in the whole problem domain. The generalized smoothed Galerkin weak form is then used to create the discretized system equations where the system stiffness is obtained using simple summation operation. The essential rotational boundary conditions are imposed in the process of constructing the curvature field in conjunction with imposing the translational boundary conditions in the same way as undertaken in the standard FEM. A number of numerical examples were studied using the present formulation, including both static and free vibration analyses. The numerical results were compared with the reference ones together with those shown in the state-of-art literatures published. Very good accuracy has been achieved using the present method.