A method proposed earlier, relying on the use of harmonic Cartesian polynomial and rational functions, is extended here to find a semi-analytical solution to the uncoupled, two-dimensional problem of thermo-magnetoelasticity for a system of long parallel, non-intersecting, transversely isotropic elastic cylindrical electrical conductors. Results are presented for two conductors of equal circular normal cross-sections carrying currents of equal densities flowing along the same direction, subjected to Robin-type thermal boundary conditions. Quantities of practical interest are represented graphically and discussed. Consideration of a system of electrical conductors is of practical importance in power plants and in various technological instruments, where it is required to assess the interaction between conductors. The obtained formulas for the magnetic vector potential may be of importance for the determination of the coefficients of self- and mutual inductance of long electric conductors, otherwise difficult to calculate by standard methods. Comparing the results with those of a single conductor allows us to assess the interaction between conductors.