自己做的轴承静力学分析，请大家帮忙找找问题在哪?(有命令流）
[attach]164877[/attach]还不会贴图，先把命令流发上来。（我把图压缩了，在附件里面）命令流：
*SET,db,19 !外圈的直径
*SET,ds,5 !内圈的直径
*SET,b,6 !宽度
*SET,dw,2.5 !珠子直径
*SET,dwp,12.5 !珠子中心圆直径
*SET,z,11 !珠子数目
*SET,de,15 !外圈沟道直径
*SET,dde,13.8 !外圈挡边直径
*SET,re,1.28 !外圈沟道曲率半径
*SET,di,10 !内圈沟道直径
*SET,ddi,11.3 !内圈高挡边直径
*SET,ddii,10.4 !内圈低挡边直径
*SET,ri,1.28 !内圈沟道曲率半径
*SET,l,0 !游系
*SET,size1,0.1 !钢球接触处单元长度
*SET,size2,0.1 !游系
/PREP7
MP,EX,1,2.8e7 !定义材料
MP,PRXY,1,0.4
MP,MU,1,0.07
MP,EX,2,3.2e11
MP,PRXY,2,0.25
ET,1,SOLID185 !定义单元类型
ET,2,TARGE170
ET,3,CONTA174
KEYOPT,3,9,0
KEYOPT,3,10,1
wpro,,,-90.000000 !建立轴承模型
/VIEW,-1,-1
RECTNG,-b/2,b/2,ds/2,ddi/2,
CYL4, ,di/2+ri,ri
ASBA,1,2
RECTNG,-b/2,0,ddii/2,ddi/2
ASBA,3,1
RECTNG,-b/2,b/2,dde/2,db/2,
CYL4, ,de/2-re,re
ASBA, 1, 3
RECTNG,0,b/2,dde/2,de/2
ASBA, 4, 1
k,17,,,
k,18,,,b
K,19,0,dwp/2,0
CYL4,0,dwp/2,3*dw/5
FLST,3,2,4,ORDE,2
FITEM,3,3
FITEM,3,-4
ASBL, 3,P51X
FLST,3,2,4,ORDE,2
FITEM,3,7
FITEM,3,10
ASBL, 2,P51X
ADELE, 1, , ,1
FLST,2,2,5,ORDE,2
FITEM,2,4
FITEM,2,-5
AGLUE,P51X
FLST,2,2,5,ORDE,2
FITEM,2,3
FITEM,2,6
AGLUE,P51X
FLST,2,4,5,ORDE,2
FITEM,2,3
FITEM,2,-6
FLST,8,2,3
FITEM,8,17
FITEM,8,18
VROTAT,P51X, , , , , ,P51X, ,360/22, ,
FLST,2,4,5,ORDE,2
FITEM,2,3
FITEM,2,-6
FLST,8,2,3
FITEM,8,17
FITEM,8,18
VROTAT,P51X, , , , , ,P51X, ,-360/22, ,
VADD,3,7
VADD,4,8
VADD,2,6
VADD,1,5
SPH4, ,dwp/2,dw/2
NUMCMP,KP
NUMCMP,LINE
NUMCMP,AREA
NUMCMP,VOLU
wpoff,0,6.25,0
VSEL,S, , , 1
VPLOT
VSBW, 1
wprot,0,0,-90
FLST,2,2,6,ORDE,2
FITEM,2,6
FITEM,2,-7
VSBW,P51X
wprot,0,90,0
FLST,2,4,6,ORDE,3
FITEM,2,1
FITEM,2,8
FITEM,2,-10
VSBW,P51X
ALLSEL,ALL
FLST,2,2,6,ORDE,2
FITEM,2,4
FITEM,2,-5
VGLUE,P51X
FLST,2,2,6,ORDE,2
FITEM,2,2
FITEM,2,-3
VGLUE,P51X
FLST,5,8,6,ORDE,4 !划分网格
FITEM,5,6
FITEM,5,-7
FITEM,5,11
FITEM,5,-16
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 2, , 1, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
ESIZE,0.3,0,
MSHAPE,0,3D
MSHKEY,1
!*
FLST,5,8,6,ORDE,4
FITEM,5,6
FITEM,5,-7
FITEM,5,11
FITEM,5,-16
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
ALLSEL,ALL
VSEL,S, , , 4
VPLOT
CM,_Y,VOLU
VSEL, , , , 4
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 1, , 1, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
ESIZE,0.2,0,
CM,_Y,VOLU
VSEL, , , , 4
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VSWEEP,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
ALLSEL,ALL
VSEL,S, , , 2
VPLOT
CM,_Y,VOLU
VSEL, , , , 2
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VSWEEP,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
ALLSEL,ALL !!!!!
VPLOT
CM,_Y,VOLU
VSEL, , , , 5
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 1, , 1, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
ESIZE,0.4,0,
CM,_Y,VOLU
VSEL, , , , 5
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VSWEEP,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!* !!!!!!!!!
ALLSEL,ALL
VPLOT
CM,_Y,VOLU
VSEL, , , , 3
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VSWEEP,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
FLST,5,4,5,ORDE,4 !建立解除对
FITEM,5,3
FITEM,5,61
FITEM,5,64
FITEM,5,67
ASEL,S, , ,P51X
NSLA,S,1
CM,no1,NODE
FLST,5,6,5,ORDE,4
FITEM,5,12
FITEM,5,-14
FITEM,5,37
FITEM,5,-39
ASEL,S, , ,P51X
APLOT
NSLA,S,1
NPLOT
CM,no2,NODE
ALLSEL,ALL
FLST,5,4,5,ORDE,4
FITEM,5,4
FITEM,5,62
FITEM,5,65
FITEM,5,68
ASEL,S, , ,P51X
NSLA,S,1
NPLOT
CM,no3,NODE
ALLSEL,ALL
FLST,5,8,5,ORDE,5
FITEM,5,2
FITEM,5,6
FITEM,5,-8
FITEM,5,30
FITEM,5,-33
ASEL,S, , ,P51X
NSLA,S,1
NPLOT
CM,no4,NODE
EPLOT
/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.07
MAT,1
R,3
REAL,3
ET,4,170
ET,5,174
KEYOPT,5,9,0
KEYOPT,5,10,1
R,3,
RMORE,
RMORE,,0
RMORE,0
! Generate the target surface
NSEL,S,,,NO1
CM,_TARGET,NODE
TYPE,4
ESLN,S,0
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
NSEL,S,,,NO2
CM,_CONTACT,NODE
TYPE,5
ESLN,S,0
ESURF
!* Create Companion Pair - Start
R,4
REAL,4
ET,6,170
ET,7,174
KEYOPT,7,9,0
KEYOPT,7,10,1
R,4,
RMORE,
RMORE,,0
RMORE,0
TYPE,6
ESEL,S,TYPE,,5
NSLE,S
ESLN,S,0
ESURF
TYPE,7
ESEL,S,TYPE,,4
NSLE,S
ESLN,S,0
ESURF
!* Create Companion Pair - End
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,4
ESEL,A,TYPE,,5
ESEL,R,REAL,,3
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,4
ESEL,A,TYPE,,5
ESEL,R,REAL,,3
ESEL,A,TYPE,,6
ESEL,A,TYPE,,7
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.07
MAT,1
R,5
REAL,5
ET,8,170
ET,9,174
KEYOPT,9,9,0
KEYOPT,9,10,1
R,5,
RMORE,
RMORE,,0
RMORE,0
! Generate the target surface
NSEL,S,,,NO3
CM,_TARGET,NODE
TYPE,8
ESLN,S,0
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
NSEL,S,,,NO4
CM,_CONTACT,NODE
TYPE,9
ESLN,S,0
ESURF
!* Create Companion Pair - Start
R,6
REAL,6
ET,10,170
ET,11,174
KEYOPT,11,9,0
KEYOPT,11,10,1
R,6,
RMORE,
RMORE,,0
RMORE,0
TYPE,10
ESEL,S,TYPE,,9
NSLE,S
ESLN,S,0
ESURF
TYPE,11
ESEL,S,TYPE,,8
NSLE,S
ESLN,S,0
ESURF
!* Create Companion Pair - End
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,8
ESEL,A,TYPE,,9
ESEL,R,REAL,,5
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,8
ESEL,A,TYPE,,9
ESEL,R,REAL,,5
ESEL,A,TYPE,,10
ESEL,A,TYPE,,11
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
DELTIM,0.2,0.05,0.5
PRED,ON,,ON
TIME,2
ALLSEL,ALL
CSYS,1
FLST,2,2,5,ORDE,2
FITEM,2,20
FITEM,2,45
/GO
DA,P51X,ALL,
FLST,2,2,5,ORDE,2
FITEM,2,21
FITEM,2,46
!*
/GO
DA,P51X,UZ,
FLST,2,2,5,ORDE,2
FITEM,2,27
FITEM,2,52
!*
/GO
DA,P51X,UZ,
FLST,2,2,5,ORDE,2
FITEM,2,24
FITEM,2,49
!*
/GO
DA,P51X,UZ,
FLST,2,2,5,ORDE,2
FITEM,2,22
FITEM,2,28
!*
/GO
DA,P51X,UY,
FLST,2,2,5,ORDE,2
FITEM,2,47
FITEM,2,53
!*
/GO
DA,P51X,UY,
FLST,2,2,5,ORDE,2
FITEM,2,23
FITEM,2,48
/GO
!*
SFA,P51X,1,PRES,0.5
[[i] 本帖最后由 muouou 于 2008-6-18 10:36 编辑 [/i]] 出现什么错误？ 我不清楚问题出在哪，告诉我怎么贴图吧
我截了几个图想发上来 plotctrl—>hardcopy—>to file 谢谢shv lee75，又学了一招。
接触对我是用接触向导直接创建的，结果不知为什么就不收敛。 你的滚子是如何约束的？ 滚珠没施加约束，怎么给它施加约束我还不懂。
这种是不是就是不收敛啊？
可能的原因是什么？ 1图表示你的计算是收敛的 [size=4]为什么不用半个钢球？用半个钢球可以通过对称模拟整个钢球，可以减小计算量。[/size]
[size=4]在初始调整中选择CLOSE GAP试一下。[/size] 我选了close gap，结果变成这个样子了
[[i] 本帖最后由 muouou 于 2008-6-25 09:45 编辑 [/i]]
回复 1# 的帖子
你做的珠子怎么被分为8份，而且相互之间没有接触对另外你作的轴承套圈与珠子之间的部分应该叫镶块吧，它与套圈之间也应建立接触对
还有多余约束存在
总之，你建立的模型问题挺多 珠子分八份是为了划分六面体网格，以前做过这种例子，应该没什么问题
“另外你作的轴承套圈与珠子之间的部分应该叫镶块吧，它与套圈之间也应建立接触对”
这句话不太明白，能再解释详细一点吗？ 自己感觉结果有点像了，但出现了新问题。加不同的载荷，应力变化非常小
加1N的载荷和40N的载荷，只是从87045变成87049。
另外球的最大应力和内外圈的最大应力不一样。外圈是从22103变成22104
可能的原因是什么？
[[i] 本帖最后由 muouou 于 2008-6-28 17:14 编辑 [/i]] [quote]原帖由 [i]muouou[/i] 于 2008-6-27 10:05 发表 [url=http://www.simwe.com/forum/redirect.php?goto=findpost&pid=1368783&ptid=837394][img]http://www.simwe.com/forum/images/common/back.gif[/img][/url]
珠子分八份是为了划分六面体网格，以前做过这种例子，应该没什么问题
“另外你作的轴承套圈与珠子之间的部分应该叫镶块吧，它与套圈之间也应建立接触对”
这句话不太明白，能再解释详细一点吗？ [/quote]
我执行了你的命令流发现在套圈和珠子之间还有一个与滚道形状相似的体，我还以为是陶瓷或其他特殊材料的镶块呢，难道不是吗？ [quote]原帖由 [i]muouou[/i] 于 2008-6-28 17:12 发表 [url=http://www.simwe.com/forum/redirect.php?goto=findpost&pid=1369815&ptid=837394][img]http://www.simwe.com/forum/images/common/back.gif[/img][/url]
自己感觉结果有点像了，但出现了新问题。加不同的载荷，应力变化非常小
加1N的载荷和40N的载荷，只是从87045变成87049。
另外球的最大应力和内外圈的最大应力不一样。外圈是从22103变成22104
可能的原因是什么？ [/quote]
这个结果应该是不对的，检查一下你所加的约束，可能存在过定位的情况
页:
[1]