guo721 发表于 2008-9-28 20:25

请高手给分析一下!!

以下是我们计算的程序,目的是通过这个程序得到体积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
页: [1]
查看完整版本: 请高手给分析一下!!