When studying some gauge theories approach to problems in Mechanics, I've found the following integral
$$P\exp\left[\oint A \ dt\right]=1+\dfrac{1}{2}\oint_{\partial D}\sum_{\mu,\nu}F_{\mu\nu}\gamma^{\mu}(t) \dot{\gamma}^{\nu}(t)dt,$$
where $A$ is the gauge potential, and $F$ is the field strength tensor (i.e. the pull-back of the curvature two form by a certain choice of gauge map). This integral appeared, in the articles (like this one, on page 564) I've seem, on the computation of the path-ordered exponential, but I couldn't understand where it comes from. It seems, on this formula, that we are integrating over a path, but $F$ is a $2$-form, so it should be integrated over a $2$-chain.
On the article there is one derivation, but I really didn't understand what they did, it doesn't seem very rigorous. Also, when I've studied principal fiber bundles and connections on those bundles, I didn't see this integral. I've also searched on some math books and didn't find it.
So, where this integral comes from, what it rigorously means and how it relates to the path-ordered exponential?