The purpose of this study is to simulate the early age concrete behaviors and evaluate the cracking risk with the thermal and thermal stress analysis. A new finite element method program associated with ANSYS program is developed for the computation of thermal field and thermal stress field for early age concrete considering the following characters: degree of hydration, thermal properties (such as specific heat, thermal diffusivity), thermal boundary conditions, and mechanical properties (such as shrinkage, creep) which occur at early age. The results from simulation compared with experimental values found in the literature show a good agreement. Finally, based on this user-developed subroutine, the effects of hydration heat, ambient temperature, wind velocity, shrinkage, and length-height ratio on cracking risk were analyzed for a concrete wall which is one part of the structure of Maridal culvert in Norway. By which, the measures to control the cracking were provided for the engineering application.