For convenience write
$$ \phi_U(x) = \| x - (K-U)\|. $$
Let $N = |\mathcal{U}|$ the number of open sets in the cover.
Define
$$ f_1(x) = \frac{1}{N} \sum \phi_U(x), \qquad f_2(x) = \left( \frac1N \sum \phi_U^2(x) \right)^{1/2} $$
Note that since $\mathcal{U}$ is a cover we have
$$ 0 < f_1(x) \leq f_2(x) \leq \sqrt{N} f_1(x) $$
Let
$$ \delta_1 = \inf f_1, \quad \delta_2 = \inf f_2 $$
and we also have $0 < \delta_1 \leq \delta_2 \leq \sqrt{N} \delta_1 $. Note also that $\delta_1$ is a Lebesgue number of the covering $\mathcal{U}$.
Now, you've defined
$$ \psi_U(x) = \frac{\phi_U^2(x)}{N f_2^2(x)} $$
so
$$ N(\psi_U(x) - \psi_U(y)) = \frac{\phi_U^2(x) f_2^2(y) - \phi_U^2(y) f_2^2(x) }{f_2(x) f_2(y)} $$
$$ = \frac{(\phi_U(x) - \phi_U(y))(\phi_U(x) + \phi_U(y))}{f_2^2(x)} + \frac{\phi_U^2(y) (f_2^2(y) - f_2^2(x))}{f_2(x)^2 f_2(y)^2}. $$
So we have the following estimate
$$ N |\psi_U(x) - \psi_U(y)| \leq \left|\frac{(\phi_U(x) - \phi_U(y))(\phi_U(x) + \phi_U(y))}{f_2^2(x)}\right| + N \left| \frac{f_2^2(y) - f_2^2(x)}{f_2(x)^2} \right| $$
Now, it suffices to consider only those $x,y$ with $\|x-y\| \leq \delta \leq \delta_1$ for a fixed $\delta$. This is because
- We know that $0 \leq \psi_U(x) \leq 1$ and hence if $\|x-y\| > \delta$, then the corresponding difference quotient is $\leq \frac{1}{\delta}$.
- On the other hand, if we define
$$ \epsilon = \inf_{U\in \mathcal{U}} d( U \setminus \cup (\mathcal{U} \setminus \{U\}), K \setminus U ) $$
then essentially the same argument of Pietro Majer shows that the Lipschitz constant of the family is at least $1/\epsilon$. (Note: minimality of the covering is used to get that $U\setminus \cup(\mathcal{U}\setminus \{U\}) \neq\emptyset$.) Note that $\epsilon \geq \delta_1$. So we don't expect a much better upper bound estimate of the Lipschitz constant than $1/\delta$.
Now, $\phi_U(x)$ is well-known to be $1$-Lipschitz. So we have that $$\phi_U(x) + \phi_U(y) \leq 2 \phi_U(x) + \delta \leq (2N+1) f_1(x)$$
So
$$ \left|\frac{(\phi_U(x) - \phi_U(y))(\phi_U(x) + \phi_U(y))}{f_2^2(x)}\right| \leq \frac{2N+1}{\delta_2} \|x - y\| $$
Similarly, we have
$$ |f^2_2(x) - f^2_2(y)| \leq \frac1N \sum |\phi_U(x) - \phi_U(y)| |2\phi_U(x) + \delta| \leq \|x-y\| \cdot (\delta + 2 f_1(x)) $$
so
$$ N \left| \frac{f_2^2(y) - f_2^2(x)}{f_2(x)^2} \right| \leq \frac{3N}{\delta_2} \|x-y\|$$
So putting everything together we find that the Lipschitz constant of the family $\psi_U$ are uniformly bounded by $6 / \delta_2$.
Summary
- For the finite collection $\mathcal{U}$, each function in the corresponding partition of unity is necessarily Lipschitz.
- The uniform Lipschitz constant can be controlled by a quantity similar to the Lebesgue number of the covering.