This paper presents a thermo-hydro-mechanics theory and corresponding computational framework to capture the freezethaw action of frozen porous media and associated frost action under chilled gas pipelines. Based on the mixture theory, frost-susceptible soils are formulated to capture the Darcy flux and thermal actions below the pipelines. Constitutive models that combine the cryo-suction are presented to reproduce the changes in volume, strength, and thermal characteristics of solid grain, pore water, and ice crystal. A generalized hardening rule is adopted to replicate the elasto-plastic responses which strengthens the frozen porous media due to ice crystallization. Changes in permeability and thermal diffusivity are also incorporated by considering the phase transitions of pore water and ice crystal. Numerical examples for pipeline applications are designed to analyze the influence of the freezing and melting process around the pipelines.