This paper presents a method of calculating the stress intensity factor (K) and crack propagation direction (${ heta}_0$) at the crack-tip that is associated with delamination in the large scale integration(LSI) package. To establish a reasonable strength evaluation method and life prediction, it is necessary to assess fracture parameters under various fracture conditions. Therefore, we conducted quantitative stress singularity analysis considering thermal stress simulating the changes of crack length (a), (h) and (v) in delamination using the 2-dimensional elastic boundary element method (BEM), and from these results predicted crack propagation direction and path.