(16a) is inspired by Newton’s law of universal gravitation in physics, where r is the distance from the
origin, γ is positive parameter, ε is a small positive constant to prevent division by zero at the origin.
(16b) is the Gaussian shape function, where ζ is positive parameter to control the decrease.
The VFC external force is: