滚动轴承接触问题

pp1352

做的是轧机的轴承,简化了保持架,现在做不下去了,在接触对的参数开始感觉问题很大,温度加载就是在滚动体表面加了70度,轴承内圈加了40度,这个分析就是为后面的动力分析做耦合分析,还请高手的指点修改下



!*  
/COM,   
/COM,Preferences for GUI filtering have been set to display:
/COM,  Structural   
/COM,  Thermal  
!*  
/filname,CAD-3-0,1
!设置单元和材料属性
/PREP7   
ET,1,SOLID70
mp,kxx,1,66.6
!建模
/Prep7
k,1,40.00000,410.0000,0.000000
k,2,40.00000,384.6791,0.000000
k,3,135.0000,362.7319,0.000000
k,4,230.0000,384.6791,0.000000
k,5,230.0000,410.0000,0.000000
k,6,0.000000,335.0000,0.000000
k,7,0.000000,280.0000,0.000000
k,8,134.7500,280.0000,0.000000
k,9,134.7500,309.8622,0.000000
k,10,109.7500,309.8622,0.000000
k,11,31.81165,325.0000,0.000000
k,12,33.93722,335.0000,0.000000
k,13,135.2500,280.0000,0.000000
k,14,135.2500,309.8622,0.000000
k,15,160.2500,309.8622,0.000000
k,16,238.1883,325.0000,0.000000
k,17,236.0628,335.0000,0.000000
k,18,270.0000,335.0000,0.000000
k,19,270.0000,280.0000,0.000000
k,20,38.04900,354.3444,0.000000
k,21,111.4101,338.7511,0.000000
k,22,105.4472,310.6980,0.000000
k,23,231.9510,354.3444,0.000000
k,24,158.5899,338.7511,0.000000
k,25,164.5528,310.6980,0.000000
k,26,0,0,0
k,27,270,0,0

!点生成线
l,1,2
l,2,3
l,3,4
l,4,5
l,5,1
l,6,7
l,7,8
l,8,9
l,9,10
l,10,11
l,11,12
l,12,6
l,13,14
l,14,15
l,15,16
l,16,17
l,17,18
l,18,19
l,19,13

!线生成面
al,1,2,3,4,5
al,6,7,8,9,10,11,12,
al,13,14,15,16,17,18,19
!划分单元
/Prep7
!外圈
lsel,s,,,2,3,1
LESIZE,all, , ,10, , , , ,1
lsel,s,,,5
LESIZE,all, , ,20, , , , ,1
lsel,s,,,1
lsel,a,,,4
LESIZE,all, , ,3, , , , ,1
lsel,s,,,10
lsel,a,,,15
LESIZE,all, , ,10, , , , ,1
lsel,s,,,6
lsel,a,,,8
LESIZE,all, , ,4, , , , ,1
lsel,s,,,12
lsel,a,,,17
LESIZE,all, , ,3, , , , ,1

lsel,s,,,11
lsel,a,,,16
LESIZE,all, , ,1, , , , ,1
lsel,s,,,7
lsel,a,,,19
LESIZE,all, , ,15, , , , ,1
lsel,s,,,9
lsel,a,,,14
LESIZE,all, , ,2, , , , ,1

lsel,s,,,8
lsel,a,,,13
LESIZE,all, , ,3, , , , ,1

!绕26－27转生面体
vrotat,1,2,3,,,,26,27,,
lsel,s,,,25,29,1
LESIZE,all, , ,30, , , , ,1
lsel,s,,,63,67,1
LESIZE,all, , ,30, , , , ,1
lsel,s,,,101,105,1
LESIZE,all, , ,30, , , , ,1
lsel,s,,,134,138,1
LESIZE,all, , ,30, , , , ,1
lsel,s,,,134,138,1
LESIZE,all, , ,30, , , , ,1
lsel,s,,,75,81,1
LESIZE,all, , ,30, , , , ,1
lsel,s,,,89,95,1
LESIZE,all, , ,30, , , , ,1
lsel,s,,,113,119,1
LESIZE,all, , ,30, , , , ,1
lsel,s,,,127,133,1
LESIZE,all, , ,30, , , , ,1
lsel,s,,,139,152,1
LESIZE,all, , ,30, , , , ,1
lsel,s,,,37,43,1
LESIZE,all, , ,30, , , , ,1
lsel,s,,,51,57,1
LESIZE,all, , ,30, , , , ,1
vsweep,1
vsweep,2
vsweep,3
vsweep,4
vsweep,5
vsweep,6
vsweep,7
vsweep,8
vsweep,9
vsweep,10
vsweep,11
vsweep,12
!定义工作平面
kwplan,1,20,21,11
!旋转工作平面
wprota,,,90
!定义局部坐标系
cswpla,11,1
!在工作面上生成锥体
con4,,,30,28.6798,75

定义工作平面
kwplan,1,23,24,11
!旋转工作面
wprota,,,90
!定义局部坐标系
cswpla,12,1
!在工作面上生成锥体
con4,,,30,28.6798,75
!压缩定义的项目编码
numcmp,all

!滚子
lsel,s,,,153,160,1
lsel,a,,,163,170,1
LESIZE,all, , ,5, , , , ,1
lsel,s,,,161,162,1
lsel,a,,,171,172,1
LESIZE,all, , ,20, , , , ,1
vsweep,13
vsweep,14

!坐标系
csys,0
k,109,0,0,100
!三个关键点定义局部坐标系的
cskp,13,1,26,7,109,1,1

!激活先前定义的坐标系
csys,13

!滚子
vgen,30,13,14,1,0,12,0,,0,0

!生成接触对
!滚子与滚道
ALLSEL,ALL  
!*  
/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.1
MAT,1   
R,3
REAL,3  
ET,3,170
ET,4,174
KEYOPT,4,9,1
KEYOPT,4,10,1   
R,3,
RMORE,  
RMORE,,0
RMORE,0
! Generate the target surface   
!选择面
ASEL,S,,,5,71,22
CM,_TARGET,AREA
TYPE,3  
NSLA,S,1
ESLN,S,0
ESURF   
CMSEL,S,_ELEMCM
!Generate the contact surface  
NSEL,NONE   
csys,13
nsel,s,loc,y,-2,2
nsel,a,loc,y,10,14
nsel,a,loc,y,22,26
nsel,a,loc,y,34,38
nsel,a,loc,y,46,50
nsel,a,loc,y,58,62
nsel,a,loc,y,70,74
nsel,a,loc,y,82,86
nsel,a,loc,y,94,98
nsel,a,loc,y,106,110
nsel,a,loc,y,118,122
nsel,a,loc,y,130,134
nsel,a,loc,y,142,146
nsel,a,loc,y,154,158
nsel,a,loc,y,166,170
nsel,a,loc,y,178,182
nsel,a,loc,y,190,194
nsel,a,loc,y,202,206
nsel,a,loc,y,214,218
nsel,a,loc,y,226,230
nsel,a,loc,y,238,242
nsel,a,loc,y,250,254
nsel,a,loc,y,262,266
nsel,a,loc,y,274,278
nsel,a,loc,y,286,290
nsel,a,loc,y,298,302
nsel,a,loc,y,310,314
nsel,a,loc,y,322,326
nsel,a,loc,y,334,338
nsel,a,loc,y,346,350

ASEL,s,,,92,324,8
NSLA,R  
CM,_CONTACT,NODE
TYPE,4  
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.1
MAT,1   
R,4
REAL,4  
ET,5,170
ET,6,174
KEYOPT,6,9,1
KEYOPT,6,10,1   
R,4,
RMORE,  
RMORE,,0
RMORE,0
! Generate the target surface   
ASEL,S,,,6,72,22
CM,_TARGET,AREA
TYPE,5  
NSLA,S,1
ESLN,S,0
ESURF   
CMSEL,S,_ELEMCM
! Generate the contact surface  
NSEL,NONE   
csys,13
nsel,s,loc,y,-2,2
nsel,a,loc,y,10,14
nsel,a,loc,y,22,26
nsel,a,loc,y,34,38
nsel,a,loc,y,46,50
nsel,a,loc,y,58,62
nsel,a,loc,y,70,74
nsel,a,loc,y,82,86
nsel,a,loc,y,94,98
nsel,a,loc,y,106,110
nsel,a,loc,y,118,122
nsel,a,loc,y,130,134
nsel,a,loc,y,142,146
nsel,a,loc,y,154,158
nsel,a,loc,y,166,170
nsel,a,loc,y,178,182
nsel,a,loc,y,190,194
nsel,a,loc,y,202,206
nsel,a,loc,y,214,218
nsel,a,loc,y,226,230
nsel,a,loc,y,238,242
nsel,a,loc,y,250,254
nsel,a,loc,y,262,266
nsel,a,loc,y,274,278
nsel,a,loc,y,286,290
nsel,a,loc,y,298,302
nsel,a,loc,y,310,314
nsel,a,loc,y,322,326
nsel,a,loc,y,334,338
nsel,a,loc,y,346,350
ASEL,s,,,95,327,8
NSLA,R  
CM,_CONTACT,NODE
TYPE,6  
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   
!*  
!*  
/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.1
MAT,1   
R,5
REAL,5  
ET,7,170
ET,8,174
KEYOPT,8,9,1
KEYOPT,8,10,1   
R,5,
RMORE,  
RMORE,,0
RMORE,0
! Generate the target surface   
ASEL,S,,,14,58,22
ASEL,a,,,79
CM,_TARGET,AREA
TYPE,7  
NSLA,S,1
ESLN,S,0
ESURF   
CMSEL,S,_ELEMCM
! Generate the contact surface  
NSEL,NONE   
csys,13
nsel,s,loc,y,-2,2
nsel,a,loc,y,10,14
nsel,a,loc,y,22,26
nsel,a,loc,y,34,38
nsel,a,loc,y,46,50
nsel,a,loc,y,58,62
nsel,a,loc,y,70,74
nsel,a,loc,y,82,86
nsel,a,loc,y,94,98
nsel,a,loc,y,106,110
nsel,a,loc,y,118,122
nsel,a,loc,y,130,134
nsel,a,loc,y,142,146
nsel,a,loc,y,154,158
nsel,a,loc,y,166,170
nsel,a,loc,y,178,182
nsel,a,loc,y,190,194
nsel,a,loc,y,202,206
nsel,a,loc,y,214,218
nsel,a,loc,y,226,230
nsel,a,loc,y,238,242
nsel,a,loc,y,250,254
nsel,a,loc,y,262,266
nsel,a,loc,y,274,278
nsel,a,loc,y,286,290
nsel,a,loc,y,298,302
nsel,a,loc,y,310,314
nsel,a,loc,y,322,326
nsel,a,loc,y,334,338
nsel,a,loc,y,346,350
ASEL,s,,,91,323,8
NSLA,R  
CM,_CONTACT,NODE
TYPE,8  
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.1
MAT,1   
R,6
REAL,6  
ET,9,170
ET,10,174   
KEYOPT,10,9,1   
KEYOPT,10,10,1  
R,6,
RMORE,  
RMORE,,0
RMORE,0
! Generate the target surface   
ASEL,S,,,20,64,22
ASEL,a,,,84
CM,_TARGET,AREA
TYPE,9  
NSLA,S,1
ESLN,S,0
ESURF   
CMSEL,S,_ELEMCM
! Generate the contact surface  
NSEL,NONE   
csys,13
nsel,s,loc,y,-2,2
nsel,a,loc,y,10,14
nsel,a,loc,y,22,26
nsel,a,loc,y,34,38
nsel,a,loc,y,46,50
nsel,a,loc,y,58,62
nsel,a,loc,y,70,74
nsel,a,loc,y,82,86
nsel,a,loc,y,94,98
nsel,a,loc,y,106,110
nsel,a,loc,y,118,122
nsel,a,loc,y,130,134
nsel,a,loc,y,142,146
nsel,a,loc,y,154,158
nsel,a,loc,y,166,170
nsel,a,loc,y,178,182
nsel,a,loc,y,190,194
nsel,a,loc,y,202,206
nsel,a,loc,y,214,218
nsel,a,loc,y,226,230
nsel,a,loc,y,238,242
nsel,a,loc,y,250,254
nsel,a,loc,y,262,266
nsel,a,loc,y,274,278
nsel,a,loc,y,286,290
nsel,a,loc,y,298,302
nsel,a,loc,y,310,314
nsel,a,loc,y,322,326
nsel,a,loc,y,334,338
nsel,a,loc,y,346,350
ASEL,S,,,96,328,8
NSLA,R  
CM,_CONTACT,NODE
TYPE,10
ESLN,S,0
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,9  
ESEL,A,TYPE,,10
ESEL,R,REAL,,6  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,9  
ESEL,A,TYPE,,10
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   
!挡边
/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.1
MAT,1   
R,13
REAL,13
ET,23,170   
ET,24,174   
KEYOPT,24,9,1   
KEYOPT,24,10,1  
R,13,   
RMORE,  
RMORE,,0
RMORE,0
! Generate the target surface   
ASEL,S,,,15,59,22
ASEL,a,,,80
CM,_TARGET,AREA
TYPE,23
NSLA,S,1
ESLN,S,0
ESURF   
CMSEL,S,_ELEMCM
! Generate the contact surface  
asel,s,,,89,321,8
CM,_CONTACT,AREA
TYPE,24
NSLA,S,1
ESLN,S,0
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,23
ESEL,A,TYPE,,24
ESEL,R,REAL,,13
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,23
ESEL,A,TYPE,,24
ESEL,R,REAL,,13
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.1
MAT,1   
R,14
REAL,14
ET,25,170   
ET,26,174   
KEYOPT,26,9,1   
KEYOPT,26,10,1  
R,14,   
RMORE,  
RMORE,,0
RMORE,0
! Generate the target surface   
ASEL,S,,,85
ASEL,a,,,21,65,22
CM,_TARGET,AREA
TYPE,25
NSLA,S,1
ESLN,S,0
ESURF   
CMSEL,S,_ELEMCM
! Generate the contact surface  
asel,s,,,93,325,8
CM,_CONTACT,AREA
TYPE,26
NSLA,S,1
ESLN,S,0
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,25
ESEL,A,TYPE,,26
ESEL,R,REAL,,14
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,25
ESEL,A,TYPE,,26
ESEL,R,REAL,,14
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 PROPERTIES - START
KEYOPT,4,1,2
KEYOPT,6,1,2
KEYOPT,8,1,2
KEYOPT,10,1,2  
KEYOPT,24,1,2
KEYOPT,26,1,2
KEYOPT,34,1,2
KEYOPT,36,1,2
KEYOPT,38,1,2
KEYOPT,40,1,2
KEYOPT,42,1,2
KEYOPT,44,1,2
KEYOPT,46,1,2
KEYOPT,48,1,2
KEYOPT,66,1,2
KEYOPT,68,1,2
KEYOPT,72,1,2
/COM, CONTACT PAIR PROPERTIES - END
/COM, CONTACT PAIR PROPERTIES - START   
RMODIF,3,10,0.203
RMODIF,4,10,0.404
RMODIF,5,10,0.203
RMODIF,6,10,0.404
RMODIF,13,10,-0.19421414   
RMODIF,14,10,-0.19422324     
/COM, CONTACT PAIR PROPERTIES - END  
/COM, CONTACT PAIR PROPERTIES - START   
RMODIF,3,14,200
RMODIF,4,14,200
RMODIF,5,14,200
RMODIF,6,14,200
RMODIF,13,14,200
RMODIF,14,14,200
RMODIF,18,14,200
RMODIF,19,14,200
RMODIF,20,14,200
RMODIF,21,14,200
RMODIF,22,14,200
RMODIF,23,14,200
RMODIF,24,14,200
RMODIF,25,14,200
RMODIF,34,14,0
RMODIF,35,14,0   
RMODIF,37,14,0   
/COM, CONTACT PAIR PROPERTIES - END
!加载
!温度载荷
ALLSEL,ALL
!tref,30
!温度
ASEL,s,,,95,327,8
ASEL,a,,,96,328,8
ASEL,a,,,92,324,8
ASEL,a,,,91,323,8
nsla,s
nplot
d,all,temp,70
ASEL,S,,,11,55,22
asel,a,,,76
ASEL,a,,,24,68,22
asel,a,,,88
nsla,s
asel,s,,,93,325,8
asel,a,,,89,321,8
nslv,a
d,all,temp,40
nsel,s,loc,x,410
d,all,temp,30

!求解
/SOL
solcontrol,on
cnvtol,temp,0.1
toffst,273
ANTYPE,0
NLGEOM,1
NSUBST,50,100,35
TIME,1  
ALLSEL,ALL