L为桩长,r为计算点与桩轴线的水平距离,z为计算点距离地面深度。都取正值 %f1,f2,f3依次是点荷载、均布荷载以及三角形荷载下的Geddes应力解函数 %Ip,Ir,It依次是对应三个荷载下的应力系数,负表示拉,正表示压