The squaring tube process is examined by an incremental elasto-plastic finite-element method based on an updated Lagrangian formulation in which a sliding-sticking friction mode is specifically considered. The high nonlinearity of the process due to the geometric changes, the inelastic constitutive behavior, and the deformation-dependent boundary conditions are taken into account in an incremental manner. A static explicit approach to the solution is applied, the tangent stiffness matrix equation is solved without iteration and a weighting factor r min is employed to limit the step size to linear relation. The simulated geometries of squaring clearly demonstrate the processes of square tube until unloading. The formation of squaring defects both collapse and asymmetry are reported in a theoretical manner. Accordingly, the effects of various parameters of the process, such as geometric ratio R/t, strain hardening exponent n, and the friction coefficient µ, on the occurrence of collapse (collapse ratio C/t) and on the extent of asymmetry (deviation ratio C1/C2) for the squaring process are discussed and interpreted in simulation. Mainly it is expected that formation of a square tube for industrial use that does not collapse will be found during the design stage, before trials begin. The present work may be expected to improve the understanding of the formation of the square tube.