Here's a solution that ends up checking fewer than $30$ cases. I don't know if it would count as an analytic solution, but it's at least partly analytic.
Note that you don't need to approach this as a search through all 9000 four-digit numbers. Instead view it as a search through all sets of four digits $\{a,b,c,d\}$, seeing if $(a+b+c+d)(a^2+b^2+c^2+d^2)^2$ is some arrangement of those digits.
You established that no digit can be $7$ or higher. So not only is the number at most $9999$, it is at most $6666$. Suppose there is a $6$. Even $(6+0+0+0)(6^2+0^2+0^2+0^2)^2>6666$, so there is no $6$ either. So not only is the number at most $6666$, it is at most $5555$.
Suppose there is a $5$. But $(5+2+0+0)(5^2+2^2+0^2+0^2)^2>5555$ and $(5+1+1+1)(5^2+1^2+1^2+1^2)^2>5555$, and so the only possibilities with a $5$ are $\{5,1,1,0\}$ with $(5+1+1+0)(5^2+1^2+1^2+0^2)^2=5103$; $\{5,1,0,0\}$ with $(5+1+0+0)(5^2+1^2+0^2+0^2)^2=4056$; and $\{5,0,0,0\}$ with $(5+0+0+0)(5^2+0^2+0^2+0^2)^2=3125$. So there is no $5$.
So not only is the number at most $5555$, it is at most $4444$.
Meanwhile, consider if the sum of digits is $10$. The minimal value of $(a+b+c+d)(a^2+b^2+c^2+d^2)^2$ would come from making the digits as close together as possible: $\{3,3,2,2\}$. But $(3+3+2+2)(3^2+3^2+2^2+2^2)^2>4444$. So the sum of the digits is at most $9$.
Now that the digits are from $\{0,1,2,3,4\}$ and the sum must be at most $9$, there are not that many combinations left to check:
- $\{4,4,*,*\}$ even with 0s, the result is more than $4444$.
- $\{4,3,0,0\}$ gives $4375$
- $\{4,3,*,*\}$ even with 1,0, the result is more than $4444$.
- $\{4,2,0,0\}$ gives $2400$, a winner!
- $\{4,2,1,0\}$ gives $3087$
- $\{4,2,1,1\}$ gives $3872$
- $\{4,2,2,*\}$ even with 0, the result is more than $4444$.
- $\{4,1,0,0\}$ gives $1445$
- $\{4,1,1,0\}$ gives $1944$
- $\{4,1,1,1\}$ gives $2527$
Any further solutions have no $4$, so are at most $3333$. The sum of the digits now can be at most $8$.
- $\{3,3,2,0\}$ gives $3827$
- $\{3,3,1,1\}$ gives $3200$
- $\{3,3,1,0\}$ gives $2527$
- $\{3,3,0,0\}$ gives $1944$
- $\{3,2,2,1\}$ gives $2592$
- $\{3,2,2,0\}$ gives $2023$, a winner!
- $\{3,2,1,1\}$ gives $1575$
- $\{3,2,1,0\}$ gives $1176$
- $\{3,2,0,0\}$ gives $845$
- $\{3,1,*,*\}$ even with 1,1 is too small
- $\{2,2,2,2\}$ gives $2048$
- $\{2,2,2,1\}$ gives $1183$
- $\{2,2,1,1\}$ gives $600$ and everything else to check is too small.