Structure-preserving approach to multi-port dynamic compact models of nonlinear heat conduction