Dear @Kevin.Chang,
Ah, I see. I would suggest defining the linear stiffness as a representative stiffness, e.g., as determined by the secant method at some representative value of y. And then I would implement the nonlinearities of the p-y curve as differences relative to this line.
Best regards,