[step:Write the map in associated bundle charts and prove smoothness]
Let $x_0\in M$. Choose an open neighbourhood $U\subset M$ of $x_0$ and a smooth local section
\begin{align*}
s:U\to P
\end{align*}
of $\pi_P$. Choose an open neighbourhood $V\subset N$ of $f(x_0)$ with $f(U)\subset V$ and a smooth local section
\begin{align*}
t:V\to Q
\end{align*}
of $\pi_Q$. Define
\begin{align*}
\gamma:U\to G
\end{align*}
to be the unique smooth map satisfying
\begin{align*}
\Phi(s(x))=t(f(x))\gamma(x)
\end{align*}
for every $x\in U$. The uniqueness follows from the freeness and transitivity of the right $G$-action on each fiber of $Q$, and smoothness follows by reading $\gamma$ as the $G$-coordinate of $\Phi\circ s$ in the principal trivialization determined by $t$.
The associated bundle chart determined by $s$ is
\begin{align*}
\Psi_s:\rho_P^{-1}(U)\to U\times F,\quad [s(x)a,u]\mapsto (x,a\cdot u),
\end{align*}
where $x\in U$, $a\in G$, and $u\in F$. The associated bundle chart determined by $t$ is
\begin{align*}
\Psi_t:\rho_Q^{-1}(V)\to V\times F,\quad [t(y)b,u]\mapsto (y,b\cdot u),
\end{align*}
where $y\in V$, $b\in G$, and $u\in F$. For $x\in U$, $a\in G$, and $u\in F$,
\begin{align*}
\Phi(s(x)a)=\Phi(s(x))a=t(f(x))\gamma(x)a.
\end{align*}
Therefore
\begin{align*}
(\Psi_t\circ \Phi_F\circ \Psi_s^{-1})(x,a\cdot u)=(f(x),\gamma(x)a\cdot u).
\end{align*}
Equivalently, for $v=a\cdot u\in F$,
\begin{align*}
(\Psi_t\circ \Phi_F\circ \Psi_s^{-1})(x,v)=(f(x),\gamma(x)\cdot v).
\end{align*}
The map $(x,v)\mapsto (f(x),\gamma(x)\cdot v)$ is smooth because $f$ is smooth, $\gamma$ is smooth, and the left action map $G\times F\to F$ is smooth. Since smoothness of a map between smooth bundles is local in charts, $\Phi_F$ is smooth near every point of $P\times_G F$.
[/step]