Ansys中【新手求教】如何绘制结构中某一杆件中内力和外加荷载关系曲线图？

一个张拉整体单元的ansys命令，现在需要绘制受压杆件 编号 10（或者11和12，三个杆件内力应该是一样的）的内力与外加荷载f1之间的关系曲线，荷载f1从0开始，到-40000为止
x轴为内力，y为f1
万分感谢！！！
/prep7
pi=3.14159
A1=0.25*pi*(60*60-54*54)
A2=0.25*pi*9.3*9.3
SIGY1=235
SIGY2=1770
p1=1600
p2=10000
p3=60000
f1=-5000
et,1,link8
et,2,link10
mp,ex,1,2.06e5
mp,ex,2,1.95e5
TB,BKIN,1
TBDATA,1,SIGY1
TB,BKIN,2
TBDATA,1,SIGY2
r,1,a1
r,2,a2,p2/(a2*1.95e5)
k,1,0,0,0
k,2,750,1299,0
k,3,1500,0,0
k,4,-116.1,433,1857.2
k,5,1183,1183,1857.2
k,6,1183,-317,1857.2
L,1,2
L,2,3
L,3,1
L,4,5
L,5,6
L,6,4
L,1,4
L,2,5
L,3,6
L,1,5
L,2,6
L,3,4
LSEL,s,,,1,9
/COLOR,LINE,MRED
/REPLOT
LSEL,s,,,10,12
/COLOR,LINE,BLUE
/REPLOT
Lsel,all
/RGB,INDEX,100,100,100, 0
/RGB,INDEX, 80, 80, 80,13
/RGB,INDEX, 60, 60, 60,14
/RGB,INDEX, 0, 0, 0,15
/REPLOT
LSEL,s,,,1,9
LESIZE,all,,,1
LATT,2,2,2
LMESH,all
LSEL,s,,,10,12
LESIZE,all,,,1
LATT,1,1,1
LMESH,all
allsel,all
/SOLU
TIME,1
NSUBST,1000
OUTRES,all,all
NLGEOM,1
DK,1,all
DK,2,uy
Dk,2,uz
DK,3,uz
fk,4,fz,f1
fk,5,fz,f1
fk,6,fz,f1
allsel,all
solve