stan123lee,
I don't know how to define the upper limit of spring capacity either. For that, it would be needed a limit for displacement, upto which the law of variation applies and I don't know where can that be input (or if).
But, I thought of another way. It is not sophisticated, instead it may take some aditional time to model and analyze: use nonlinear hinge. This is what I thought of:
- define supports, like columns, with a separate material and input there elasticity modulus and strength of soil. Put dimension of supports such as to cover the full supported area
- go to define > frame nonlinear hinge properties and define a nonlinear hinge for AxialP (instead of using the default generated by ETABS, you can Modify/show for P the bilinear variation there)
- assign > frame line > frame nonlinear hinges > TheHinge to the supports (the columns defined as supports)
- define > static nonlinear/pushover case and add a case using your loads and afferent coefficients. Monitor on UZ a relevant point (like the location of a support)
- run the elastic analysis, then the static nonlinear analysis and Display > show static pushover curve for a graphically display of the effect of increasing load on the nonlinear behaving soil.
I advise to ask for a second opinion on this procedure. I have never verified its validity, but I hope it is useful for you. Please let me know if this works.