Diffeomorphic approximation of piecewise affine homeomorphisms

Given any $f$ a locally finitely piecewise affine homeomorphism of $\Omega \subset \mathbb{R}^d$ onto $\Delta \subset \mathbb{R}^d$ (for $d=3, 4$) such that $f\in W^{1,p}(\Omega, \mathbb{R}^d)$ and $f^{-1}\in W^{1,q}(\Delta, \mathbb{R}^d)$, $1\leq p ,q < \infty$ and any $\epsilon >0$ we construct a diffeomorphism $\tilde{f}$ such that $$\|f-\tilde{f}\|_{W^{1,p}(\Omega,\mathbb{R}^d)} + \|f^{-1}-\tilde{f}^{-1}\|_{W^{1,q}(\Delta,\mathbb{R}^d)} < \epsilon.$$