It is of practical significance to analyze the multi-physics process of nuclear system, which includes neutronics, heat transfer and thermoelasticity. Fission reaction is the heat source in system, the heat source will affect the distribution of temperature field, which will lead to the change of strain. Strain in turn will affect the distribution of neutron field. Therefore, it is necessary to analyze the distribution of neutron flux, temperature and strain in system. Three aspects of work have been carried out: 1) Based on Galerkin finite element theory, the governing equations of neutronics, heat transfer and thermoelasticity are established; 2) The multi-physics analysis code is developed; 3) The calculation results are analyzed and discussed.