各位高人,我做了一个螺栓受力分析的文章,但总是出错,不知什么原因,烦大家帮忙看看,在线等。
!*螺栓预紧力分析
*SET,le,500
*SET,he,100
*SET,DN,89
*SET,T,4.5
*SET,D1,196
*SET,D0,80
*SET,DC,160
*SET,N,20
*SET,hod,20
*SET,NUM,4
*SET,screwh,12
*SET,screwr,18
*SET,bolth,3*n
*set,boltd,18
/PREP7
k,1,0,DN/2,0
K,2,LE,DN/2,0
K,3,LE,DN/2+HE,0
LSTR, 1, 2
LSTR, 2, 3
LFILLT,1,2,44.5, ,
wprot,0,0,90
CYL4,0,0,DN/2
CYL4,0,0,D0/2
ASBA, 1, 2
FLST,8,3,4
FITEM,8,1
FITEM,8,3
FITEM,8,2
VDRAG, 3, , , , , ,P51X
aplot
wprot,0,0,-90
wpoff,le+dn/2,he+dn/2,0
wprot,0,-90,0
!*第一个法兰
CYL4,0,0,D1/2
CYL4,0,0,Dn/2
ASBA, 29, 30
VOFFST,31,n, ,
wpoff,0,0,n
!*第二个法兰
CYL4,0,0,D1/2
CYL4,0,0,Dn/2
ASBA, 39, 40
VOFFST,41,n, ,
wpoff,0,0,n
!*第二截短管
CYL4,0,0,DN/2
CYL4,0,0,D0/2
ASBA, 49, 50
VOFFST,51,3*n, ,
!*在此做标记
!*画出第一个螺栓孔
wpoff,dc/2,0,-2*n
CYL4,0,0,hod/2
VOFFST,59,2*n, ,
!*画出第二个螺栓孔
wpoff,-dc,0,0
CYL4,0,0,hod/2
VOFFST,65,2*n, ,
!*画出第三个螺栓孔
wpoff,dc/2,dc/2,0
CYL4,0,0,hod/2
VOFFST,71,2*n, ,
!*画出第四个螺栓孔
wpoff,0,-dc,0
CYL4,0,0,hod/2
VOFFST,77,2*n, ,
!*切出四个螺栓孔
FLST,2,2,6,ORDE,2
FITEM,2,4
FITEM,2,-5
FLST,3,4,6,ORDE,2
FITEM,3,7
FITEM,3,-10
VSBV,P51X,P51X
vplot
!*画出第一个螺栓
RPR4,6,0,0,screwr,120,-screwh
CYL4,0,0,boltd/2
VOFFST,63,bolth, ,
!*画出第一个螺母
wpoff,0,0,2*n
RPR4,6,0,0,screwr,120,screwh
!*画出第二个螺栓
wpoff,0,dc,-2*n
RPR4,6,0,0,screwr,120,-screwh
CYL4,0,0,boltd/2
VOFFST,121,bolth, ,
!*画出第二个螺母
wpoff,0,0,2*n
RPR4,6,0,0,screwr,120,screwh
!*画出第三个螺栓
wpoff,dc/2,-dc/2,-2*n
RPR4,6,0,0,screwr,120,-screwh
CYL4,0,0,boltd/2
VOFFST,143,bolth, ,
!*画出第三个螺母
wpoff,0,0,2*n
RPR4,6,0,0,screwr,120,screwh
!*画出第四个螺栓
wpoff,-dc,0,-2*n
RPR4,6,0,0,screwr,120,-screwh
CYL4,0,0,boltd/2
VOFFST,165,bolth, ,
!*画出第四个螺母
wpoff,0,0,2*n
RPR4,6,0,0,screwr,120,screwh
!*合并与压缩各组件
NUMMRG,ALL, , , ,LOW
NUMCMP,ALL
!*合并各组件
FLST,2,4,6,ORDE,3
FITEM,2,1
FITEM,2,-3
FITEM,2,11
VGLUE,P51X
FLST,2,2,6,ORDE,2
FITEM,2,6
FITEM,2,12
VGLUE,P51X
FLST,2,2,6,ORDE,2
FITEM,2,4
FITEM,2,-5
VGLUE,P51X
FLST,2,2,6,ORDE,2
FITEM,2,13
FITEM,2,-14
VGLUE,P51X
FLST,2,2,6,ORDE,2
FITEM,2,8
FITEM,2,-9
VGLUE,P51X
FLST,2,2,6,ORDE,2
FITEM,2,16
FITEM,2,-17
VGLUE,P51X
!*给单元指定单元格及相应的属性
ET,1,SOLID45
!*
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,2.01e+005
MPDATA,PRXY,1,,0.28
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,7.8e-006
!*给单元格划分网格
FLST,5,18,6,ORDE,4
FITEM,5,1
FITEM,5,-15
FITEM,5,17
FITEM,5,-19
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VSWEEP,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*自由划分网格
MSHAPE,1,3D
MSHKEY,0
!*
FLST,5,4,6,ORDE,4
FITEM,5,4
FITEM,5,8
FITEM,5,13
FITEM,5,19
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*指定第一个螺栓的预紧力
PSMESH, ,m1, ,VOLU,17, 4,Z,-20, , , ,elems,
!*指定第二个螺栓的预紧力
wpoff,dc,0,0
PSMESH, ,m2, ,VOLU,14, 4,Z,-20, , , ,elems,
!*指定第三个螺栓的预紧力
wpoff,-dc/2,-dc/2,0
PSMESH, ,m3, ,VOLU,5, 4,Z,-20, , , ,elems,
!*指定第四个螺栓的预紧力
wpoff,0,dc,0
PSMESH, ,m4, ,VOLU,9, 4,Z,-20, , , ,elems,
!* 指定上法兰与下法兰的接触单元
/PREP7
CM,_TARGET,AREA
!*
CM,_CONTACT,AREA
!*
/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.2
MAT,1
MP,EMIS,1,7.88860905221e-031
R,3
REAL,3
ET,3,170
ET,4,174
R,3,,,1,0.1,0,
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,4,4,0
KEYOPT,4,5,0
KEYOPT,4,7,0
KEYOPT,4,8,0
KEYOPT,4,9,0
KEYOPT,4,10,1
KEYOPT,4,11,0
KEYOPT,4,12,0
KEYOPT,4,2,0
KEYOPT,3,5,0
! Generate the target surface
ASEL,S,,,30
ASEL,A,,,32
ASEL,A,,,33
ASEL,A,,,34
ASEL,A,,,35
ASEL,A,,,36
ASEL,A,,,37
ASEL,A,,,38
ASEL,A,,,83
ASEL,A,,,84
ASEL,A,,,85
ASEL,A,,,86
ASEL,A,,,87
ASEL,A,,,88
ASEL,A,,,89
ASEL,A,,,90
ASEL,A,,,91
ASEL,A,,,92
ASEL,A,,,93
ASEL,A,,,94
ASEL,A,,,95
ASEL,A,,,96
ASEL,A,,,97
ASEL,A,,,98
ASEL,A,,,99
ASEL,A,,,100
CM,_TARGET,AREA
TYPE,3
NSLA,S,1
ESLN,S,0
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
ASEL,S,,,40
ASEL,A,,,42
ASEL,A,,,43
ASEL,A,,,44
ASEL,A,,,45
ASEL,A,,,46
ASEL,A,,,47
ASEL,A,,,48
ASEL,A,,,100
ASEL,A,,,101
ASEL,A,,,102
ASEL,A,,,103
ASEL,A,,,104
ASEL,A,,,105
ASEL,A,,,106
ASEL,A,,,107
ASEL,A,,,108
ASEL,A,,,109
ASEL,A,,,110
ASEL,A,,,111
ASEL,A,,,112
ASEL,A,,,113
ASEL,A,,,114
ASEL,A,,,115
ASEL,A,,,116
ASEL,A,,,117
CM,_CONTACT,AREA
TYPE,4
NSLA,S,1
ESLN,S,0
ESURF
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,3
ESEL,A,TYPE,,4
ESEL,R,REAL,,3
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,3
ESEL,A,TYPE,,4
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.2
MAT,1
MP,EMIS,1,7.88860905221e-031
R,5
REAL,5
ET,7,170
ET,8,174
R,5,,,1,0.1,0,
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,8,4,0
KEYOPT,8,5,0
KEYOPT,8,7,0
KEYOPT,8,8,0
KEYOPT,8,9,0
KEYOPT,8,10,1
KEYOPT,8,11,0
KEYOPT,8,12,0
KEYOPT,8,2,0
KEYOPT,7,5,0
! Generate the target surface
ASEL,S,,,29
ASEL,A,,,77
ASEL,A,,,134
ASEL,A,,,178
CM,_TARGET,AREA
TYPE,7
NSLA,S,1
ESLN,S,0
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
ASEL,S,,,99
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,,5
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,7
ESEL,A,TYPE,,8
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.2
MAT,1
MP,EMIS,1,7.88860905221e-031
R,4
REAL,4
ET,5,170
ET,6,174
R,4,,,1,0.1,0,
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,6,4,0
KEYOPT,6,5,0
KEYOPT,6,7,0
KEYOPT,6,8,0
KEYOPT,6,9,0
KEYOPT,6,10,1
KEYOPT,6,11,0
KEYOPT,6,12,0
KEYOPT,6,2,0
KEYOPT,5,5,0
! Generate the target surface
ASEL,S,,,70
ASEL,A,,,127
ASEL,A,,,149
ASEL,A,,,171
CM,_TARGET,AREA
TYPE,5
NSLA,S,1
ESLN,S,0
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
ASEL,S,,,117
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,,4
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,5
ESEL,A,TYPE,,6
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
!*指定四个螺栓的预紧力
*DEL,_z
*DIM,_z,ARRAY,1,1,1
*SET,_z(1) ,1
*DO,_z1,1,1
SLOAD,_z(_z1),9,LOCK,FORC, 100,1,2
*ENDDO
*DEL,_z
*DIM,_z,ARRAY,1,1,1
*SET,_z(1) ,2
*DO,_z1,1,1
SLOAD,_z(_z1),9,LOCK,FORC, 100,1,2
*ENDDO
*DEL,_z
*DIM,_z,ARRAY,1,1,1
*SET,_z(1) ,3
*DO,_z1,1,1
SLOAD,_z(_z1),9,LOCK,FORC, 100,1,2
*ENDDO
*DEL,_z
*DIM,_z,ARRAY,1,1,1
*SET,_z(1) ,4
*DO,_z1,1,1
SLOAD,_z(_z1),9,LOCK,FORC, 100,1,2
*ENDDO
*DEL,_z
*DIM,_z,ARRAY,1,1,1
*SET,_z(1) ,4
*DO,_z1,1,1
SLOAD,_z(_z1),9,LOCK,FORC, 100,1,2
*ENDDO
!*施加载荷及计算应力
FLST,2,1,5,ORDE,1
FITEM,2,49
DA,P51X,ALL,0
ACEL,0,-9.8,0,
FLST,2,1,5,ORDE,1
FITEM,2,3
/GO
SFA,P51X,1,PRES,10
!*求解运算结果
/SOL
/STATUS,SOLU
SOLVE
[ 本帖最后由 jinqinghai 于 2006-5-22 08:20 编辑 ] |