Drag Load variation along the pile and how to stimulate this on Plaxis 2D

Hi all, 

   Dotted line shows the capacity if there is no any NSF

My question is how to stimulate these things in Plaxis 2D and if any ways to find the neutral plane through analysis?

Thanks.