<p>In this paper, we propose a modified permeable model to analyze the fracture behavior of two collinear cracks under a thermo-electro-elastic load that is described by a cubic function. By employing this modified model and Fourier transformation, the mixed boundary value problem is firstly reduced to a set of singular integral equations. Then, based on using the Cauchy kernel theory, we derive analytical solutions for these equations, with a focus on key crack-tip physical quantities and intensity factors. Numerical results demonstrate that the crack size, initial unit heat flux adjustment parameters, and initial electric displacement adjustment parameter, accounting for the effects of initial heat flux, electric displacement, and irregular crack surfaces, significantly influence the temperature gradients and intensity factors.</p>