ANSYS单桩非线性求解疑问

我是一个ANSYS小菜鸟，朋友给我一个单桩非线性例子，但是我该如何进行后处理得到Q-S曲线等图形呢？非常感谢~附命令流：
finish
/clear
/filname,foundation
/title，pile
/config,nres,2000000
/units,si !国际单位制
fe=0.3 !摩擦系数
fk1=10 !法向接触刚度
fk2=1 !法向接触刚度
fd=3 !切向接触刚度
ftoln=1 !初始渗透因子
A=3.1415926*0.6*0.6 !桩的横截面积
/prep7
et,1,solid45
mp,mu,1,0.2 !桩体材料属性
mp,ex,1,2.5e10
mp,nuxy,1,0.2
mp,dens,1,2500
mp,ex,2,2.0e8 !土层2材料属性
mp,nuxy,2,0.4
mp,dens,2,2000
mp,ex,3,2.5e8 !土层3材料属性
mp,nuxy,3,0.4
mp,dens,3,2000
mp,ex,4,3.0e8 !土层4材料属性
mp,nuxy,4,0.4
mp,dens,4,2000
mp,ex,5,1.8e9 !岩层5材料属性
mp,nuxy,5,0.29
mp,dens,5,2600
tb,concr,1 !桩体的参数
tbdata,,0.7,0.9,2.6,-1
tb,dp,2 !土层2的D-P参数
tbdata,1,19,32,30
tb,dp,3 !土层3的D-P参数
tbdata,1,19,32,30
tb,dp,4 !土层4的D-P参数
tbdata,1,19,32,30
tb,dp,5 !土层5的D-P参数
tbdata,1,0,35,0
cylind,0,0.6,0,21.0,0,90
cylind,0.6,3.0,0,3.0,0,90
cylind,0.6,3.0,3.0,8.0,0,90
cylind,0.6,3.0,8.0,14.0,0,90
cylind,0.6,3.0,14.0,21.0,0,90
cylind,0,0.6,21.0,24.0,0,90 !建立桩尖岩层几何体
cylind,0.6,3.0,21.0,24.0,0,90 !布尔运算粘结体所有的除桩之外的部分vsel,s,,,2,7,1
vglue,2,3,4,5,6,7
allsel,all
numcmp,all
/pnum,volu,1
vplot
vlist
save
/pnum,line,1
allsel,all
vsel,s,,,1 !桩体附属性，划分网格
aslv,s
lsla,s
lplot
lesize,1,0.1
lesize,2,0.1
lesize,3,,,6
lesize,4,,,6
lesize,5,0.1
lesize,6,0.1
lesize,7,0.5
lesize,8,0.5
lesize,9,0.5
type,1
mat,1
mshape,0,3
mshkey,1
vmesh,all
vsel,s,,,2 !土层2附属性，划分网格
vplot
aslv,s
lsla,s
lplot
lesize,10,,,6
lesize,11,,,6
lesize,12,,,6
lesize,13,,,6
lesize,14,,,6
lesize,15,,,6
lesize,16,,,6
lesize,17,,,6
lesize,18,0.5
lesize,19,0.5
lesize,20,0.5
lesize,21,0.5
type,1
mat,2
mshape,0,3
mshkey,1
vmesh,all
vsel,s,,,4 !土层3附属性，并划分网格
vplot
aslv,s
lsla,s
lplot
lesize,22,,,6
lesize,23,,,6
lesize,24,,,6
lesize,25,,,6
lesize,43,0.5
lesize,44,0.5
lesize,45,0.5
lesize,46,0.5
type,1
mat,3
mshape,0,3
mshkey,1
vmesh,all
vsel,s,,,5 !土层4附属性，并划分网格
vplot
aslv,s
lsla,s
lplot
lesize,26,,,6
lesize,27,,,6
lesize,28,,,6
lesize,29,,,6
lesize,47,0.5
lesize,48,0.5
lesize,49,0.5
lesize,50,0.5
type,1
mat,4
mshape,0,3
mshkey,1
vmesh,all
vsel,s,,,6 !岩层5附属性，并划分网格
vplot
aslv,s
lsla,s
lplot
lesize,30,,,6
lesize,31,,,6
lesize,32,,,6
lesize,33,,,6
lesize,51,0.5
lesize,52,0.5
lesize,53,0.5
lesize,54,0.5
type,1
mat,5
mshape,0,3
mshkey,1
vmesh,all
vsel,s,,,3 !桩底岩层6附属性，并划分网格
vplot
aslv,s
lsla,s
lplot
lesize,34,,,6
lesize,35,,,6
lesize,36,,,6
lesize,37,,,6
lesize,41,,,6
lesize,42,,,6
type,1
mat,5
mshape,0,3d
mshkey,1
vmesh,all
vsel,s,,,7 !桩底岩层6附属性，并划分网格
vplot
aslv,s
lsla,s
lplot
lesize,38,,,6
lesize,57,,,6
lesize,58,,,6
lesize,41,0.5
lesize,42,0.5
lesize,55,0.5
lesize,56,0.5
type,1
mat,5
mshape,0,3
mshkey,1
vmesh,all!用接触向导建立接触对
et,2,targe170
et,3,conta173
keyopt,3,4,1
keyopt,3,9,0 !消除初始渗透
keyopt,3,12,2 !不分开接触，允许滑动，法向不分开!桩侧，建立目标面r,1,,,fk1,ftoln, , ,rmore,,,,,,5
r,2,,,fk2,ftoln, , ,rmore,,,,,,5
mp,mu,6,0.2
allsel,all
vsel,s,,,1
aslv,s
asel,s,,,3
/pnum,area,1
aplot
cm,target1,area
type,2
real,1
mat,6
nsla,s,1
esurf,all! 桩侧，建立接触面
allsel
vsel,s,,,2
vsel,a,,,4,6,1
aslv,s
asel,s,,,9
asel,a,,,22
asel,a,,,26
asel,a,,,30
aplot
cm,contact1,area
type,3
real,1
mat,6
nsla,s,1
esurf,all !桩尖，建立目标面
allsel,all
vsel,s,,,3
aslv,s
asel,r,,,15
aplot
nsla,s,1
type,2
real,2
mat,6
cm,target2,elem
esurf,all
esel,s,type,,2 eplot !桩尖，建立接触面
allsel,all
vsel,s,,,1
aslv,s
asel,s,,,2
aplot
nsla,s,1
type,3
real,2
mat,6
cm,contact1,area
esurf,all
esel,s,type,,3
eplot
/psymb,esys,1
allsel
gplot
save
finish