Recently, growing interests in frozen ground have stimulated us to advance fundamental theories and systematic researches on soil behavior under freezing conditions. Unlike the well-established soil mechanics theory, temperature variation and phase change of pore-water cause water migration to cold side, ground heaving, sharp increase in earth pressure, etc., which bring about serious problems in frozen geotechnical structures. Elasto-plastic mechanical constitutive model for frozen/unfrozen soil subjected to fully coupled THM phenomena is formulated based on a new stress variable that is continuous in frozen-unfrozen transitional regions. Numerical simulations are conducted to discuss numerical reliability and applicability of the developed constitutive model: one-dimensional heaving pressure, tri-axial compression test, and one-side freezing tests. The numerical results show that developed model can efficiently describe complex THM phenomena of frozen soil, and they can be utilized to analyze and design the geotechnical structures under freezing conditions, and predict their long-term behavior.