请高手给分析一下！！

guo721

以下是我们计算的程序，目的是通过这个程序得到体积3或5上面的反作用力，实际上是有反作用力的，但通过该程序的计算没有这个力，烦请高手给指点一下。谢谢！
/PREP7  
!*  
ET,1,SOLID95
!*  
ET,2,SOLID92
!*  
  
!*  
!*  
!*  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,2.1E11
MPDATA,PRXY,1,,0.3  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,7800
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,2,,2.1E11
MPDATA,PRXY,2,,0.3  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,2,,7800
TB,RATE,1,1,2,1
TBTEMP,0
TBDATA,,0.3419,67337.1333,,,,   
TB,MISO,1,1,5,0
TBTEMP,0
TBPT,,0.00231,4.854E8   
TBPT,,0.0159,4.92E8
TBPT,,0.02585,5.1694E8  
TBPT,,0.04059,5.4570003E8   
TBPT,,0.10063,5.875448E8
K,1,-0.020,0.020,0.7023676,
K,2,-0.0038985,0.020,0.7770086,
K,3,-0.0038985,-0.020,0.7770086,   
K,4,-0.020,-0.020,0.7023676,   
K,5,0.020,0.020,0.7023676,
K,6,0.0038985,0.020,0.7770086,
K,7,0.0038985,-0.020,0.7770086,
K,8,0.020,-0.020,0.7023676,
LSTR,       1,       2  
LSTR,       2,       3  
LSTR,       3,       4  
LSTR,       4,       1  
LSTR,       5,       6  
LSTR,       6,       7  
LSTR,       7,       8  
LSTR,       8,       5  
FLST,2,4,4  
FITEM,2,1   
FITEM,2,2   
FITEM,2,3   
FITEM,2,4   
AL,P51X
FLST,2,4,4  
FITEM,2,5   
FITEM,2,8   
FITEM,2,7   
FITEM,2,6   
AL,P51X
K,9,0,0.020,0.775964,   
LARC,2,6,9,0.004036,   
ADRAG,       2, , , , , ,       9   
CYL4,0,0,0.00725, , , ,0.785
FLST,3,3,5,ORDE,2   
FITEM,3,1   
FITEM,3,-3  
VSBA,       1,P51X  
VDELE,       3, , ,1
LDELE,       9  
NUMCMP,ALL  
CYL4,0,0,0.00725, , , ,-0.290
K,22,0.050,0.005,0.780,  
K,23,0.050,0.005,0.800,  
K,24,-0.050,0.005,0.800,
K,25,-0.050,0.005,0.780,
K,26,-0.050,-0.005,0.780,
K,27,-0.050,-0.005,0.800,
K,28,0.050,-0.005,0.800,
K,29,0.050,-0.005,0.780,
/VIEW,1,1,1,1   
V,      23,      22,      25,      24,      28,      29,      26,      27   
K,46,-0.00009,0.005,0.810,   
K,47,-0.00009,0.005,0.793,   
K,48,0.00009,0.005,0.793,
K,49,0.00009,0.005,0.810,
K,50,-0.00009,-0.005,0.810,  
K,51,-0.00009,-0.005,0.793,  
K,52,0.00009,-0.005,0.793,   
K,53,0.00009,-0.005,0.810,   
K,54,0,0.005,0.790
K,55,0,-0.005,0.790   
V,      46,      47,      48,      49,      50,      51,      52,      53   
V,      47,      54,      48,      51,      55,      52
VPLOT   
FLST,2,2,6,ORDE,2   
FITEM,2,4   
FITEM,2,-5  
VADD,P51X   
VSBV,       3,       6  
wpoff,0,0,0.8025
wprot,90,90,90  
wpoff,0.040,0,0
CYLIND,0.0025,0,-0.007,0.007,0,360,  ！体积3或5
wpoff,-0.080,0,0   
CYLIND,0.0025,0,-0.007,0.007,0,360,  ！体积3或5
NUMCMP,ALL  
TYPE,   2   
MAT,       2
REAL,   
ESYS,       0   
SECNUM,
!*  
ESIZE,0.005,0,  
MSHAPE,1,3D
MSHKEY,0
!*  
CM,_Y,VOLU  
VSEL, , , ,       1
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
VPLOT   
CM,_Y,VOLU  
VSEL, , , ,       2
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
VPLOT   
TYPE,   1   
MAT,       1
REAL,   
ESYS,       0   
SECNUM,
!*  
ESIZE,0.003,0,  
CM,_Y,VOLU  
VSEL, , , ,       4
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VSWEEP,_Y1  
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
VPLOT   
TYPE,   2   
MAT,       2
REAL,   
ESYS,       0   
SECNUM,
!*  
ESIZE,0.002,0,  
FLST,5,2,6,ORDE,2   
FITEM,5,3   
FITEM,5,5   
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*     
!*  
!*  
/COM, CONTACT PAIR CREATION - START
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,2,0   
MAT,2   
R,3
REAL,3  
ET,5,170
ET,6,174
KEYOPT,6,9,0
KEYOPT,6,10,2   
R,3,
RMORE,  
RMORE,,0
RMORE,0
! Generate the target surface   
ASEL,S,,,8  
CM,_TARGET,AREA
TYPE,5  
NSLA,S,1
ESLN,S,0
ESLL,U  
ESEL,U,ENAME,,188,189   
ESURF   
CMSEL,S,_ELEMCM
! Generate the contact surface  
ASEL,S,,,1  
CM,_CONTACT,AREA
TYPE,6  
NSLA,S,1
ESLN,S,0
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,5  
ESEL,A,TYPE,,6  
ESEL,R,REAL,,3  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,5  
ESEL,A,TYPE,,6  
ESEL,R,REAL,,3  
CMSEL,A,_NODECM
CMDEL,_NODECM   
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM   
CMSEL,S,_AREACM
CMDEL,_AREACM   
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM   
/GRES,cwz,gsav  
CMDEL,_TARGET   
CMDEL,_CONTACT  
/COM, CONTACT PAIR CREATION - END   
!*  
/COM, CONTACT PAIR CREATION - START
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,1,0   
MAT,1   
R,4
REAL,4  
ET,7,170
ET,8,174
KEYOPT,8,9,0
KEYOPT,8,10,2   
R,4,
RMORE,  
RMORE,,0
RMORE,0
! Generate the target surface   
ASEL,S,,,4  
CM,_TARGET,AREA
TYPE,7  
NSLA,S,1
ESLN,S,0
ESLL,U  
ESEL,U,ENAME,,188,189   
ESURF   
CMSEL,S,_ELEMCM
! Generate the contact surface  
ASEL,S,,,13
CM,_CONTACT,AREA
TYPE,8  
NSLA,S,1
ESLN,S,0
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,7  
ESEL,A,TYPE,,8  
ESEL,R,REAL,,4  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,7  
ESEL,A,TYPE,,8  
ESEL,R,REAL,,4  
CMSEL,A,_NODECM
CMDEL,_NODECM   
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM   
CMSEL,S,_AREACM
CMDEL,_AREACM   
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM   
/GRES,cwz,gsav  
CMDEL,_TARGET   
CMDEL,_CONTACT  
/COM, CONTACT PAIR CREATION - END  
/COM, CONTACT PAIR CREATION - START
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,1,0   
MAT,1   
R,5
REAL,5  
ET,9,170
ET,10,174   
KEYOPT,10,9,0   
KEYOPT,10,10,2  
R,5,
RMORE,  
RMORE,,0
RMORE,0
! Generate the target surface   
ASEL,S,,,19
CM,_TARGET,AREA
TYPE,9  
NSLA,S,1
ESLN,S,0
ESLL,U  
ESEL,U,ENAME,,188,189   
ESURF   
CMSEL,S,_ELEMCM
! Generate the contact surface  
ASEL,S,,,22
CM,_CONTACT,AREA
TYPE,10
NSLA,S,1
ESLN,S,0
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,9  
ESEL,A,TYPE,,10
ESEL,R,REAL,,5  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,9  
ESEL,A,TYPE,,10
ESEL,R,REAL,,5  
CMSEL,A,_NODECM
CMDEL,_NODECM   
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM   
CMSEL,S,_AREACM
CMDEL,_AREACM   
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM   
/GRES,cwz,gsav  
CMDEL,_TARGET   
CMDEL,_CONTACT  
/COM, CONTACT PAIR CREATION - END   
!*  
/COM, CONTACT PAIR CREATION - START
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,1,0   
MAT,1   
R,6
REAL,6  
ET,11,170   
ET,12,174   
KEYOPT,12,9,0   
KEYOPT,12,10,2  
R,6,
RMORE,  
RMORE,,0
RMORE,0
! Generate the target surface   
ASEL,S,,,25
CM,_TARGET,AREA
TYPE,11
NSLA,S,1
ESLN,S,0
ESLL,U  
ESEL,U,ENAME,,188,189   
ESURF   
CMSEL,S,_ELEMCM
! Generate the contact surface  
ASEL,S,,,17
CM,_CONTACT,AREA
TYPE,12
NSLA,S,1
ESLN,S,0
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,11
ESEL,A,TYPE,,12
ESEL,R,REAL,,6  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,11
ESEL,A,TYPE,,12
ESEL,R,REAL,,6  
CMSEL,A,_NODECM
CMDEL,_NODECM   
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM   
CMSEL,S,_AREACM
CMDEL,_AREACM   
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM   
/GRES,cwz,gsav  
CMDEL,_TARGET   
CMDEL,_CONTACT  
/COM, CONTACT PAIR CREATION - END  

FINISH  
/SOL
!*  
ANTYPE,4
!*  
TRNOPT,FULL
LUMPM,0
ANTYPE,4                    
NLGEOM,1
DELTIM,0.5e-6,0.5e-6,0.5e-6
!CONTROL OF OUTPUT
OUTRES,NSOL,1   
OUTRES,RSOL,1   
OUTRES,ESOL,1   
OUTRES,NLOAD,1  
OUTRES,STRS,1   
OUTRES,EPEL,1   
OUTRES,EPPL,1  
AUTOTS,1
TIME,3.6e-4  

VSEL,S, , ,       2
NSLV,S,1
FLST,2,6795,1,ORDE,2
FITEM,2,20047   
FITEM,2,-26841  
IC,P51X,UZ, ,11,

ALLSEL,ALL  
FLST,5,2,6,ORDE,2   
FITEM,5,1   
FITEM,5,-2  
VSEL,S, , ,P51X
NSLV,S,1
FLST,2,26841,1,ORDE,2   
FITEM,2,1   
FITEM,2,-26841  
!*  
/GO
D,P51X, , , , , ,UX,UY, , , ,   

ALLSEL,ALL  
FLST,5,2,6,ORDE,2   
FITEM,5,3   
FITEM,5,5   
VSEL,S, , ,P51X
NSLV,S,1
FLST,2,802,1,ORDE,2
FITEM,2,32036   
FITEM,2,-32837  
!*  
/GO
D,P51X, , , , , ,ALL, , , , ,   
ALLSEL,ALL  
/STATUS,SOLU
SOLVE