In this work, a general method for the mathematical description of three-dimensional trimmed surface is proposed by introducing the base parametric surface and boundary curves. Since mesh density distribution for the analysis may vary by cases, a grid-based mesh generation algorithm using quadtree is proposed in the present work. For the assurance of connectivity of generated meshes among surfaces, a method for the pre-cleaning of boundary curves has been developed to be used in the automatic generation of the finite elements. In addition, mesh-smoothing algorithm is suggested which can be used in the general trimmed surface. In this algorithm nodes are moved on the original surface by the normal projection in each iterative smoothing procedure.