Analytic function expression:
F*cos(atan2(py,abs(px)))

Boundary Load expression:
load(-P0,Y-PinHoleY,Z)*(sign(X)*Z>0)