This paper investigates the seismic performance of thin-walled stiffened steel square box columns, modeling bridge piers subjected to unidirectional cyclic lateral loading with a constant axial load, focusing on local, global, and local-global interactive buckling phenomena. Initially, the finite element model was validated against existing experimental results. The study further explored the degradation in strength and ductility of both thin-walled and compact columns under cyclic loading. Thin-walled, stiffened steel square box columns exhibited buckling near the base, forming a half-sine wave shape. The research also addresses discrepancies from different material models used to analyze steel tubular bridge piers. Analysis using a modified two-surface plasticity model (2SM) yielded results closer to experimental data than a multi-linear kinematic hardening model, particularly for compact sections. The 2SM, which accounts for cycling within the yield plateau and strain hardening regime, demonstrated enhanced accuracy over the multi-linear kinematic hardening model. Additionally, a parametric study was conducted to assess the impact of key design parameters—such as width-to-thickness ratio (Rf), column slenderness ratio (λ), and magnitude of axial load (P/Py)—on the performance of thin-walled stiffened steel square box columns. Design equations were then developed to predict the strength and ductility of bridge piers. These equations closely matched experimental results, achieving an accuracy of 95% for ultimate strength and 97% for ductility.