This work reports a numerical investigation dealing with the post-buckling behaviour, strength, and Direct Strength Method (DSM) design of cold-formed steel fixed-ended lipped channel columns experiencing coupling between distortional and global (flexural-torsional) buckling − distortional-global (D–G) interaction. Various cross-section dimensions and lengths are considered, to ensure that the columns selected undergo different levels and types (“true”, “secondary-bifurcation distortional” or “secondary-bifurcation global”) of D–G interaction. The results presented and discussed, determined by means of Abaqus shell finite element geometrically and materially non-linear analyses, consist of elastic and elastic-plastic post-buckling equilibrium paths, failure loads and collapse modes. The steel material behaviour is deemed elastic-perfectly plastic and several yield stresses are considered, thus making it possible to cover a wide column D–G slenderness range. Particular attention is devoted to identifying the most detrimental initial geometrical imperfection shapes, in the sense that they lead to the lowest column failure loads. The numerical failure load data gathered are subsequently used to assess the merits of the available DSM-based design approaches developed to handle cold-formed steel columns undergoing D–G interaction. Since these design approaches are shown to be either inefficient and/or improvable, the above failure load data are also used to propose modifications/improvements aimed at achieving an efficient failure load prediction in the specific context of fixed-ended lipped channel columns. The success of this endeavour provides encouragement to the authors in their search for a safe, accurate and reliable DSM-based design approach capable of handling cold-formed steel columns with arbitrary cross-section shapes and/or end support conditions that fail in D–G interactive modes.