In this paper, we constrain the Cardassian expansion models from the latest observations including the updated Gamma-ray bursts (GRBs), which calibrated cosmology-independently from the Union2 compilation of type Ia supernovae (SNe Ia). By combining the GRB data to the joint observations with the Union2 SNe Ia set, along with the Cosmic Microwave Background radiation observation from the seven-year Wilkinson Microwave Anisotropy Probe result, the baryonic acoustic oscillation observation from the spectroscopic Sloan Digital Sky Survey Data Release galaxy sample, we find significant constraints on model parameters of the original Cardassian model Ω M0 = 0.282 +0.015 −0.014 , n = 0.03 +0.05 −0.05 ; and n = −0.16 +0.25 −3.26 , β = 0.76 +0.34 −0.58 of the modified polytropic Cardassian model, which are consistent with the ΛCDM model in 1-σ confidence region. From the reconstruction of the deceleration parameter q(z) in Cardassian models, we obtain the transition redshift z T = 0.73 ± 0.04 for the original Cardassian model, and z T = 0.68 ± 0.04 for the modified polytropic Cardassian model.