Effects of Kernel Parameters of Peridynamic Theory on Heat Conduction Numerical Solution for Non-Homogeneous Material

Based on peridynamic theory (PD), heat conduction PD formulas for non-homogeneous material have been presented and numerical method to calculate non-homogeneous material properties for PD is given. Discretization process of PD heat conduction numerical formula is illustrated and the numerical equations are presented. A FORTRAN program is coded to implement PD calculations. Accuracy of this PD program is validated by results of analytical method. Effects of kernel parameters (equivalent thermal conductivity proportion factor, surface correction factor and shape factor) on PD heat conduction numerical results for non-homogeneous material are primarily studied. It is shown that equivalent thermal conductivity proportion factor and surface correction factor can affect thermal response, but effect of shape factor can be negligible. When surface correction factor is taken into consider and the equivalent thermal conductivity is specified as average of two material points within a horizon, solution of temperature response of non-homogeneous material based on PD theory approaches analytical solution.