This study is aimed to analyse the laser glass bonding process numerically. Due to the viscoelastic behaviour of glass, the extremely large deformation of the frit seal is resulted continuously over the transition temperature, so that the thermal boundary condition be changed in the entire calculation process. The commercial FEM algorithm is restrictively able to remesh the large geometrical boundary shape and to adapt the boundary conditions simultaneously. According to our manual adaptation of increasing the laser line intensity to 700 mW/mm, it is possible to estimate the thermal glass bonding process under the fracture stress in principle. But it should be studied further in the case of high laser line intensity.