SUMMARYComputer-aided mesh generation (CAMG) dictated solely by the minimal key set of requirements of geometry, material, loading and support condition can produce 'mega-sized', arbitrary-shaped distorted elements. However, this may result in substantial cost saving and reduced bookkeeping for the subsequent ÿnite element analysis (FEA) and reduced engineering manpower requirement for ÿnal quality assurance. A method, denoted as c-type, has been proposed by constructively deÿning a ÿnite element space whereby the above hurdles may be overcome with a minimal number of hyper-sized elements. Bezier (and de Boor) control vectors are used as the generalized displacements and the Bernstein polynomials (and B-splines) as the elemental basis functions. A concomitant idea of coerced parametry and inter-element continuity on demand uniÿes modelling and ÿnite element method. The c-type method may introduce additional control, namely, an inter-element continuity condition to the existing h-type and p-type methods. Adaptation of the c-type method to existing commercial and general-purpose computer programs based on a conventional displacement-based ÿnite element method is straightforward. The c-type method with associated subdivision technique can be easily made into a hierarchic adaptive computer method with a suitable a posteriori error analysis. In this context, a summary of a geometrically exact non-linear formulation for the two-dimensional curved beams=arches is presented. Several beam problems ranging from truly three-dimensional tortuous linear curved beams to geometrically extremely non-linear two-dimensional arches are solved to establish numerical e ciency of the method. Incremental Lagrangian curvilinear formulation may be extended to overcome rotational singularity in 3D geometric non-linearity and to treat general material non-linearity.