SimWe仿真论坛's Archiver

COMSOL 2008年会

mgh_mgh 发表于 2006-10-14 13:27

总体说来,楼主做的比较好。尤其是把应力集中显现出来,这是关键的。你的单位如果是N/M的话,结果是正确的。能把完整的命令流文件传上来吗?

mgh_mgh 发表于 2006-10-14 13:31

l楼主你做的结果不对,与实际结论不符。咱们可以交流。

jwooh 发表于 2006-10-17 11:57

型号6206,我为了节省计算资源,把这个模型切掉了一半,分别试了一下加集中载荷和面载荷,这个约束不好弄了,低面加y向的位移约束,剖开部分加对称约束,集中力加在上面的中点处,不能求解,这个时候我又把侧面的一条线加了y位移约束。这样可以求解,但是我觉得不对。我把模型分完网格的命令流放上来,请楼上的哥们帮我看看。
*SET,roll_r,4.7625   
*SET,inside_r1,15
*SET,INSIDE_R2 , 18.24
*SET,outside_r1,27.765
*SET,OUTSIDE_R2 , 31
*SET,deep,8
*SET,con_r,roll_r/5
/PREP7  
!*  
ET,1,SOLID45
!*  
!*  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,2.06e5
MPDATA,PRXY,1,,0.3  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,MU,1,,0.3
SPH4, , ,roll_r
/REPLOT,RESIZE  
CYL4, ,inside_r2+roll_r,inside_r1,-90,inside_r2,-60,deep
CYL4, ,outside_r1-roll_r,outside_r1,-90,outside_r2,-60,deep
SPH4, ,roll_r,con_r
SPH4, ,-roll_r,con_r
FINISH  
/SOL
FINISH  
/PREP7  
FLST,2,2,6,ORDE,2   
FITEM,2,1   
FITEM,2,4   
VSBW,P51X   
/VIEW,1,1   
/ANG,1  
/REP,FAST   
FLST,2,2,6,ORDE,2   
FITEM,2,6   
FITEM,2,8   
VDELE,P51X, , ,1
wpro,,,30.000000
wpro,,,30.000000
wpro,,,30.000000
/VIEW,1,,,1
/ANG,1  
/REP,FAST   
FLST,2,2,6,ORDE,2   
FITEM,2,7   
FITEM,2,9   
VSBW,P51X   
FLST,2,2,6,ORDE,2   
FITEM,2,4   
FITEM,2,8   
VDELE,P51X, , ,1
VSBA,       2,      26  
VSBA,       1,      29  
VSBA,       8,      17  
!*  
LOCAL,11,0,0,-roll_r,, , , ,1,1,
WPCSYS,-1   
CSYS,11,
!*  
LOCAL,11,1,0,-roll_r,0, , , ,1,1,   
CSYS,11,
FLST,3,1,6,ORDE,1   
FITEM,3,5   
VGEN, ,P51X, , , ,90, , , ,1
VSBA,       3,      18  
VDELE,       6, , ,1
VDELE,       5, , ,1
/REPLOT,RESIZE
LPLOT
FLST,5,6,4,ORDE,6   
FITEM,5,4   
FITEM,5,-5  
FITEM,5,8   
FITEM,5,32  
FITEM,5,39  
FITEM,5,49  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,30,10, , , ,1
!*  
FLST,5,1,4,ORDE,1   
FITEM,5,39  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,30,1/10, , , ,1  
!*  
FLST,5,6,4,ORDE,6   
FITEM,5,14  
FITEM,5,37  
FITEM,5,-38
FITEM,5,45  
FITEM,5,50  
FITEM,5,-51
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,20, , , , ,1
!*  
MSHKEY,0
MSHAPE,1,3d
FLST,5,2,6,ORDE,2   
FITEM,5,2   
FITEM,5,4   
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
LPLOT   
FLST,5,5,4,ORDE,5   
FITEM,5,1   
FITEM,5,3   
FITEM,5,29  
FITEM,5,55  
FITEM,5,62  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,30,10, , , ,1
!*  
FLST,5,1,4,ORDE,1   
FITEM,5,31  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,30,1/10, , , ,1  
!*  
FLST,5,6,4,ORDE,6   
FITEM,5,36  
FITEM,5,41  
FITEM,5,54  
FITEM,5,56  
FITEM,5,-57
FITEM,5,61  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,20, , , , ,1
!*  
FLST,5,1,4,ORDE,1   
FITEM,5,31  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,30,10, , , ,1
!*
MSHKEY,0
MSHAPE,1,3d
FLST,5,2,6,ORDE,2   
FITEM,5,1   
FITEM,5,8   
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
LPLOT   
/AUTO,1
/REP,FAST   
MSHKEY,0
MSHAPE,1,3d
CM,_Y,VOLU  
VSEL, , , ,       9
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
MSHKEY,0
MSHAPE,1,3d
MSHKEY,0
MSHAPE,1,3d
FLST,5,2,6,ORDE,2   
FITEM,5,7   
FITEM,5,10  
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*

[[i] 本帖最后由 jinqinghai 于 2006-11-15 03:33 编辑 [/i]]

jwooh 发表于 2006-10-18 09:27

怎么没有做这类分析的同行么,我还想请教一下轴承除了做结构分析,还需要做模态分析么,还有别的需要分析别的么?请高手给我点思路!

liang7910 发表于 2006-10-19 18:27

我也是做轴承分析的。不过也是学Ansys没多长时间,你提到的这些问题,我也遇到了,也不知道怎么解决,等楼主帮我们解决吧。我的QQ:313433126.有空长联系。同病相怜啊!

jinqinghai 发表于 2006-10-23 13:50

No response for so much time !
anybody can help him?!

2 credits points will be awarded to anyone can help him.

jwooh 发表于 2006-10-23 16:17

版主:thank you very much!  But, i think no body can help me ,呜呜..........其实我的这个问题对高手来说应该不难了,模型的命令流绝对没有问题,就差做接触对,接触对用接触向导做也很快,偶就是不明白怎么加载合理,唉..............

davidzdb 发表于 2006-10-23 17:08

根据轴承的受力,我想内圈内表面应该加径向载荷,同样外圈外表面加径向约束,如果思路正确的话,加载和加约束前需要将坐标系改为柱坐标,然后将模型中节点的节点坐标系旋转到柱坐标系下,不知对否?

Evan 发表于 2006-10-24 00:12

我没用过ansys,不过看起来你的建模不对,除去轴承安装后的preload不管,轴承内圈的受力面是整个下半个圈, 中间载荷最大,两边逐渐减小。如果按照对称原理简化模型,至多可以减至左或右一半。其实我觉得可以用2D就好了。如果再减至一个滚珠, 光是粗略的模拟接触情况,那就只用上下两个平板夹中间一个圆圈就够了。注意限制模型左右translation freedom。

仅供参考。

Evan 发表于 2006-10-24 01:24

有点错误。。。轴承内圈受力不明,还是加上一个轴,加载荷在轴的截面上比较稳妥

jwooh 发表于 2006-10-25 08:57

楼上得朋友davidzdb 说将模型中节点的节点坐标系旋转到柱坐标系下,这样做是为了干什么啊,我不太懂;evan说得1/4模型,我发得命令流就是这样得,不过我不知道怎么加这个力,由零变大,不过这个问题应该好解决的,我想知道的是它均匀增大么?要是再加个轴,轴和内圈不是又有接触么,我觉得这样更复杂了。谢谢大家的帮忙!

crazyterry 发表于 2006-10-25 14:23

对于滚动轴承来说,由于各个滚动体的受力不一致,再建立虚拟模型仿真的时候,通常需要建立至少1/4模型(在轴向与径向分别切一刀而产生的1/4剖分体)。对于你的初始模型来说,在接触区域划分局部接触面是很明智的,但你所建模型过于简单,所以你如果想施加合理的边界条件是十分困难的;另外,我建议接触对不能全部选用Close Gap,因为在轴承的接触对中,这种设置不是很合实际;同时,这么一个简单模型,如果计算30多个小时,我不清楚你的机器配置如何,这个对于一个普通的P4-3.0,1G内存的机器来说,是不正常的。

jwooh 发表于 2006-10-25 22:25

crazyterry,这位朋友,谢谢你的提点,你说的挺有道理,刚开始的时候我也打算弄个1/4模型的,想省资源就作了个45度的模型,我再做做看,对于设置的问题,等我看看计算结果再来向大家汇报!

FreddyMusic 发表于 2006-10-31 15:02

内外圈应力情况完全不对,检验你的载荷!

Evan 发表于 2006-10-31 21:24

模型简化要考虑几何,材料属性和载荷的对称性,不能只看一个, 简化后的约束必须使模型等同于简化之前的状况

jwooh 发表于 2006-11-1 10:14

我现在正在做6300的静态分析,弄完后发上来让大家帮我瞧瞧!谢谢这么多热心的朋友!

FreddyMusic 发表于 2006-11-1 14:05

[quote]原帖由 [i]jwooh[/i] 于 2006/11/1 10:14 发表
我现在正在做6300的静态分析,弄完后发上来让大家帮我瞧瞧!谢谢这么多热心的朋友! [/quote]


把命令流注释和结构写清楚。我有空看看。

jwooh 发表于 2006-11-4 19:36

这个是6300的命令流

模型只有1/4的,上面有两个小球需要和内外圈接触,接触区域要细化网格,所以我就把接触的内外的面割出1mm深,然后用旋转拉伸,上下各分出了几块,这样就能把和小球接触得区域画得细一些。分网格的时候,接触的地方我用size0.1,别的地方就大一些.内圈割出了一个规则得园环体,是为了画出六面体网格用得,因为我要再上面加变化得表面力,力得大小有2480N变化到0,和1/4小球接触得内圈那个地方得力最大,一直到边缘为0。
下面得这段命令流没有网格,网格得命令流不好做,分割得份数不一样就变化了。

*SET,db,35                  !外圈的直径
*SET,ds,10                  !内圈的直径
*SET,b,11                   !宽度
*SET,dw,6.35                !珠子直径
*SET,dwp,23.5               !珠子中心圆直径  
*SET,z,7                    !珠子数目
*SET,RE       , 3.32        !外圈沟道曲率半径
*SET,RI       , 3.26        !内圈沟道曲率半径
*SET,di,17.15               !内圈沟道直径
*SET,de,29.85               !外圈沟道直径
*SET,he,3.7                 !外圈边缘高度
*SET,hi,4.7                 !内圈边缘高度
*set,cut1,6                !内圈分的份数
*set,cut2,8               !外圈分的份数
wpro,,,-90.000000   
/PREP7  
ET,1,SOLID45
!*  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,2.07e5
MPDATA,PRXY,1,,0.3  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,MU,1,,0.3
/VIEW,1,1
RECTNG,0,b/2,-db/2, -db/2+he,
CYL4, ,-de/2+re,re   
ASBA,       1,       2  
RECTNG,0,b/2,-ds/2-hi,-ds/2,  
CYL4, ,-di/2-ri,ri   
ASBA,       1,       2  
SPH4, ,-dwp/2,dw/2   
CYL4, ,-dwp/2,dw/2+1
FLST,2,3,5,ORDE,2   
FITEM,2,3   
FITEM,2,-5  
AOVLAP,P51X
ADELE,      10, , ,1
K,101,,,,   
K,102,,,b,
l,101,102   

K,103,,-dwp/2+dw/2+1,-b/2,  
K,103,,-dwp/2+dw/2+1,b/2,   
l,103,16
GPLOT   
ASBL,       9,      10  



FLST,2,2,5,ORDE,2   
FITEM,2,3   
FITEM,2,-4  
AGLUE,P51X



FLST,2,3,5,ORDE,3   
FITEM,2,3   
FITEM,2,-4  
FITEM,2,8   
FLST,8,2,3  
FITEM,8,101
FITEM,8,102
VROTAT,P51X, , , , , ,P51X, ,90, ,
CSYS,1  
FINISH  
/SOL
FINISH  
/PREP7  
FLST,3,1,6,ORDE,1   
FITEM,3,1   
VGEN,2,P51X, , , ,360/z, , ,0  
VSBW,       1   
FLST,2,4,3  
FITEM,2,21  
FITEM,2,7   
FITEM,2,1   
FITEM,2,25  
A,P51X  
FLST,2,2,6,ORDE,2   
FITEM,2,5   
FITEM,2,-6  
VSBA,P51X,       1  
FLST,2,3,6,ORDE,3   
FITEM,2,1   
FITEM,2,7   
FITEM,2,10  
VDELE,P51X, , ,1
FLST,2,1,5,ORDE,1   
FITEM,2,7   
FLST,8,2,3  
FITEM,8,101
FITEM,8,102
VROTAT,P51X, , , , , ,P51X, ,90,cut1,  
FLST,2,1,5,ORDE,1   
FITEM,2,6   
FLST,8,2,3  
FITEM,8,101
FITEM,8,102
VROTAT,P51X, , , , , ,P51X, ,90,cut2,
LDELE,       9, , ,1

jwooh 发表于 2006-11-4 19:40

深沟球轴承6300静态的最大接触应力

理论计算

6300的最大接触应力:(珠子的数量是7个)受力大小是额定静载荷3472N!

内圈P=4.03Gpa,平均接触应力为2.69Gpa,接触尺寸:椭圆长轴3.46mm,短轴0.34mm,弹性趋近量为22.4微米;

外圈P=3.35Gpa,平均接触应力2.24Gpa,接触尺寸:椭圆长轴为3.03mm,短轴为0.466mm,弹性趋近量为22.5微米。

jwooh 发表于 2006-11-4 19:58

我做得过程中出现了这些问题:

我做得过程中出现了这些问题:
把体割开得时候有缝隙,我就把工作平面得公差改为0.000003了,这对结果影响大不大还不清楚;
我用得是接触向导,面面接触,小球得接触面没有割开,直接用半球做接触,可能这个占用计算资源了;
接触对得设置,我只是设置了FKN=1,FTOLN=0.1,pmin和pmax我当时设置了0.0001和0.001,求解得时候不让设,说小球得整个面都用来接触了;
最后一个就是函数加载,我用得是正弦函数,这段模型在-90度到0之间,我把它转换到柱坐标下加,但是看不到表面力得情况,这个能不能显示出来?

FreddyMusic 发表于 2006-11-4 20:15

1. 理论计算基本正确!见我的计算书。

2. 几何模型建的还行,有限元模型需要再进行修订。
    下一步分网和加载负荷,参见
    Verification Manual | Part I. Verification Test Case Descriptions |
    VM63 Static Hertz Contact Problem

你先再琢磨琢磨,有问题欢迎探讨!

FreddyMusic 发表于 2006-11-4 20:24

[quote]原帖由 [i]jwooh[/i] 于 2006/11/4 19:58 发表
我做得过程中出现了这些问题:
把体割开得时候有缝隙,我就把工作平面得公差改为0.000003了,这对结果影响大不大还不清楚;
我用得是接触向导,面面接触,小球得接触面没有割开,直接用半球做接触,可能这个占用 ... [/quote]

先不要一上来就搞 3D 的有限元模型,先搞个 2 D 的 模型。
等你 2D 成功了, 3D问题 自然 迎刃而解。先退一步,才能更大的前进。

不要考虑滚动轴承的游隙,就是个简单二维的 hertz 接触。见下图
把几何尺寸计算正确,ANSYS 操作没有问题。

你把二维的搞好了,我再奖励你 3 分 上 simwe ftp。

FreddyMusic 发表于 2006-11-5 09:47

jwooh ,

GOOD NEWS. I have something for you. See picture.

目前只要分析一个 二维的 外圈 + 钢球

jwooh 发表于 2006-11-6 13:42

谢谢FreddyMusic,把小球和内外接触分别计算,这样我觉得挺好,还能省计算资源,我现在正在计算的是表面加的面载荷2480,没有用函数加载,(前2次函数加的不对),看看这次算得结果怎样,估计明天晚上能算完,不行先改成单个接触计算。

二维模型,我已经计算过了,不是用的实际工况,我当时考虑的是加个预紧位移,想分析一下预紧后的应力分布情况,预紧应该是在运动时候分析用的,在静力学分析预紧看不到什么效果,没法下结论。等我把这个三维的算完了,我再把原来的那个二维模型拿出来算算看!

FreddyMusic 发表于 2006-11-6 14:53

[quote]原帖由 [i]jwooh[/i] 于 2006/11/6 13:42 发表
谢谢FreddyMusic,把小球和内外接触分别计算,这样我觉得挺好,还能省计算资源,我现在正在计算的是表面加的面载荷2480,没有用函数加载,(前2次函数加的不对),看看这次算得结果怎样,估计明天晚上能算完,不行 ... [/quote]

你说得也有道理!两维的算不出 Hertz 接触椭圆。
我的建议是这样,先算一个三维的小球和外接触。
一般面载荷加载,外圈沟道约束。
先把简单模型和思路验证了,在往复杂的搞。
这周有空我也和你一起跑跑命令流。

jwooh 发表于 2006-11-6 21:55

晚上和liang7910讨论了一下接触的问题,对于接触的求解算法还不太清楚,在接触设置这块,这位朋友说接触刚度的设置应该是弹性模量的1/100,我记得好多例子都是设置比例因子,我设置的是1,要是真按弹性模量设置,我记得书上讲的是输入负号,这块内容有没有熟悉的朋友给指点一下!
FreddyMusic:谢谢你一直都关心着我的问题,我跑的这段6300的命令流,设置了50步,计算的相当慢,应该是设置的问题,你运行的时候也要注意一下,小心它收敛得慢啊!

FreddyMusic 发表于 2006-11-7 09:34

有问题!为啥你只有 SOLID45 的单元 ?

平面赫兹接触用这些单元。
2-D 8-Node Structural Solid Elements (PLANE82 )  
2-D 6-Node Triangular Structural Solid Elements (PLANE2)
3-D Node-to-Node Contact Elements (CONTA178)

有限元模型也不一样,外面壳体/接触处是定义 接触单元,内部为结构单元。
“ 接触刚度的设置应该是弹性模量的1/100。” 为啥 ?
只要定义 Material Properties E = 1000 N/mm2  υ = 0.3  即可。

接触的求解算法,我们没必要搞得太清楚,那是 FEA软件公司的事。
但是我们要把问题描述清楚定义清楚。

我感兴趣的!就是最近要出差了。
没有整块的时间来好好思考。你再试试!我回来帮你。
还有,你能把之前做的二维的命令流帖出来吗?
我过一下,打点基础。谢谢!

FreddyMusic 发表于 2006-11-7 10:42

这是 VM63 的 命令流

1. 理解掌握 平面 / 小球接触形式
2. 略作修改变成 轴承 的接触形式,沟道 / 小球接触形式

======================

你们谁干完了,贴出命令流,我给他加 3 分!

[code]

/COM,ANSYS MEDIA REL. 10.0 (05/31/2005) REF. VERIF. MANUAL: REL. 10.0
/VERIFY,VM63
/PREP7
SMRT,OFF
/TITLE, VM63, STATIC HERTZ CONTACT PROBLEM SOLVED USING CONTAC178 ELEMENTS
/COM  REF:  TIMOSHENKO AND GOODIER, THEORY OF ELASTICITY, 3RD ED., ART. 140.
ET,1,PLANE82,,,1        ! AXISYMMETRIC ELEMENTS
ET,2,PLANE2,,,1
ET,3,CONTA178,,4           ! NODAL CONTACT
R,1  
RMOD,1,7,1              !CONTACT NORMAL ALONG UY
MP,EX,1,1E3
MP,NUXY,1,.3
LOCAL,11,1,0,8,0        ! LOCAL CYLINDRICAL C.S. AT CENTERLINE
K,1,8,-90               ! DEFINE KEYPOINTS
K,2,8
K,3,7.5,-90
K,4,7.5
K,5
K,6,8,-82.65            ! PLACE KEYPOINT AND NODE AT EXPECTED CONTACT RADIUS
K,7,7.5,-82.65
L,1,3                   ! DEFINE LINES
L,2,4  
L,6,7  
LESIZE,ALL,,,1          ! DEFINE ELEMENT DIVISIONS ON ALL EXISTING LINES
A,1,6,7,3               ! DEFINE AREAS
A,6,2,4,7
A,3,7,4,5
LOCAL,12,0,0,8,0
ARSYM,Y,1,3,1           ! CREATE HALF-SYMMETRY MODEL
NUMMRG,KPOI
ESIZE,,4                ! DEFINE ELEMENT DIVISIONS ON REMAINING LINES
LESIZE,4,,,5
*REPEAT,2,1
LESIZE,6,,,8,8
LESIZE,7,,,8,(1/8)
LESIZE,10,,,1
*REPEAT,2,2
LESIZE,9,,,6,.2
TYPE,1                  ! CREATE NODES AND ELEMENTS
AMESH,1,2,1
AMESH,4,5,1
TYPE,2
AMESH,3,6,3
CSYS,0
N,1001,-1,1E-8          !NODE 1001 IS THE GROUND
D,1001,ALL              !X POSITION DOES NOT MATTER IN THIS CASE BECAUSE
                        !THE CONTACT NORMAL IS ONLY ALONG UY
TYPE,3
REAL,1
EN,205,1001,2           !USE THE SAME ELEMENT NUMBERS AS VM63 FOR POST-PROC
EN,201,1001,4
EN,202,1001,6
EN,203,1001,8
EN,204,1001,10
EN,206,1001,31

MODMSH,NOCHECK
TYPE,1                  ! REMOVE MIDSIDE NODES ALONG CONTACT SURFACE
EMODIF,7,7,0
*REPEAT,6,1
MODMSH,CHECK
FINISH
/SOLU   
NSEL,S,LOC,X,-.01,.01   ! BOUNDARY CONDITIONS AND LOADING
D,ALL,UX,0
NSEL,R,LOC,Y,0
D,ALL,UY,0
NSEL,ALL            
LOAD=0
*CREATE,LOADSTEP        ! MACRO TO INCREMENTALLY APPLY LOAD
FK,8,FY,ARG1
SOLVE
*END
*DO,I,1,3
LOAD=LOAD-10
*USE,LOADSTEP,LOAD*6.2831853
*ENDDO
FINISH
/POST1                  ! POSTPROCESS
/OUT,
SET,3
ESEL,,TYPE,,3
ETABLE,RFOR,SMISC,1
NSLE
PRETAB,RFOR             ! PRINT REACTION FORCE TO DETERMINE CONTACT AREA
SSUM                    ! SUM OF REACTION FORCE
NLIST                   ! LIST COORDINATES OF NODES OF CONTACT SURFACE
PRNSOL,U,COMP           ! LIST DISPLACEMENTS OF NODES
/COM    CALCULATE RATIO OF A - ACTUAL TO A - TARGET
PI=(4*ATAN(1))
LOAD=-(LOAD)*(2*PI)
ATAR=(0.88*((LOAD*0.008)**(1/3)))   ! A - TARGET
*GET,EMAX,ELEM,,NUM,MAX
*DO,ENUM,201,EMAX                   ! START SEARCH FROM ELEM 201
*GET,GRFR,ELEM,ENUM,ETAB,RFOR      ! FIND LAST ELEMENT IN CONTACT
*IF,GRFR,EQ,0.0,EXIT
*ENDDO
ESEL,,ELEM,,(ENUM-1)    ! SELECT LAST CONTACTING ELEMENT
NSLE                    ! SELECT NODES ATTACHED TO SELECTED ELEMENTS
*GET,NMIN,NODE,0,NUM,MIN
NODX=NX(NMIN)
NODY=NY(NMIN)
NUX =UX(NMIN)
NUY =UY(NMIN)
AACT=NODX+NUX           ! A - ACTUAL
YCHK=NODY+NUY
RATA=(AACT/ATAR)        ! RATIO
*status,parm
*DIM,LABEL,CHAR,1,2
*DIM,VALUE,,1,3
LABEL(1,1) = '      A,'
LABEL(1,2) = ' mm     '
*VFILL,VALUE(1,1),DATA,1.010
*VFILL,VALUE(1,2),DATA,AACT
*VFILL,VALUE(1,3),DATA,ABS(AACT/1.010)
/COM
/OUT,vm63,vrt,,append
/COM,---------------------------------------------------------------------------
/COM,----------VM63 RESULTS COMPARISON (OBTAINED USING CONTACT178 ELEMENTS)-----
/COM,
/COM,                 |   TARGET   |   ANSYS   |   RATIO
/COM,
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,'   ',F10.3,'  ',F10.3,'   ',1F5.3)
/COM,---------------------------------------------------------------------------
/OUT
FINISH
*LIST,vm63,vrt

[/code]

jwooh 发表于 2006-11-7 19:49

6300的计算结果

这个计算跑了好几天,结果却是很失望,可能是网格的问题,我记得划分完了之后,放大后能看到割开的缝隙,我把结果的等效应力云图,还有几个角度的网格划分贴上来,请朋友们帮忙琢磨琢磨。

jwooh 发表于 2006-11-7 19:57

局部网格是我放大了外圈粗细网格分界处的地方,有缝隙!

jwooh 发表于 2006-11-7 20:11

to:FreddyMusic   "有问题!为啥你只有 SOLID45 的单元 ?"

    我用的是接触向导做的,它会自动给我选择170和174,我在向导中更改设置就可以了。不需要提前定义出来,而且每个接触对还需对应不同的实常数。这个我想是没有问题!

“ 接触刚度的设置应该是弹性模量的1/100。” 为啥 ?
——这个是一个朋友告诉我的,从一篇博士论文找到的,说这个设置比较合理(对我的这个问题)


接触的求解算法,我们没必要搞得太清楚,那是 FEA软件公司的事。
————我指的是选择哪种算法计算好一些!

FreddyMusic 发表于 2006-11-8 09:11

看了,应力图还是不对。

应力最大值应该是在沟道下层,现在的是在网络过度处。

网络划分不理想,看VM的例子网络为是如何过渡的。

jwooh 发表于 2006-11-8 14:56

这个模型我用的参数建模,p1_high和p2_high需要提前知道该轴承的内外圈得边缘厚度,定义单元得时候只定义一个,接触部分得单元在接触设置得时候定,接触是用的线线,不过我定义接触的时候是通过接触部分的面定义的,因为我担心不同型号的轴承,可能割出的部分大小不一样,线线接触就会出问题,对于我的这个尺寸,我分别用了选择线接触,和定义组件选择接触面上的节点接触,结果一样。理的大小是2480,加在上面的线上,结果在附件中,但是我没法比较结果,这个结果算的过大,有高手明白的,请给指点一下!

jwooh 发表于 2006-11-9 17:55

三维的弄出来了,让大家看看对不对!

这个模型得加载是直接在内圈表面加了均布力,实际应该逐渐加力,最下面的力最大,到两边力逐渐为0,用函数加,我现在还不是很清楚,加了不知道怎么看加得对不对,有朋友给一些提示啊!

jwooh 发表于 2006-11-9 18:32

我加载得均布力2480N,额定静载荷3472N

wfk 发表于 2006-11-13 15:19

材料常数

轴承钢材料系数是:
弹性模量:E=2.07*E5 MPa
泊松比:    U=0.3

FreddyMusic 发表于 2006-11-13 22:28

不错有点像了!恭喜先!

我还没细看,过两天来挑刺!

welyou 发表于 2006-11-13 22:55

请教轴承接触设置

求教:分析滚子轴承的接触应力,先建立的是一个完整的轴承,然后我用切除命令建立1/4模型,接触对我是用接触向导做的,内外圈是目标面,滚子是接触面。边界条件是这样:用对称约束分别约束了内外圈剖切端面以及两个滚子的剖切面;外圈的外柱面约束三个方向的自由度。对接触区域也没有细化。分析的结果是只有内圈有应力云图,滚子和外圈都没有应力。 不知为什么设置的接触没有起作用?下面是我的命令流:
/BATCH  
/COM,ANSYS RELEASE  9.0    UP20041104       15:43:07    11/07/2006
RESUME,'zhoucheng5','dbb','.'   
/COM,ANSYS RELEASE  9.0    UP20041104       15:43:07    11/07/2006
!*  
/NOPR   
/PMETH,OFF,0
KEYW,PR_SET,1   
KEYW,PR_STRUC,1
KEYW,PR_THERM,0
KEYW,PR_FLUID,0
KEYW,PR_ELMAG,0
KEYW,MAGNOD,0   
KEYW,MAGEDG,0   
KEYW,MAGHFE,0   
KEYW,MAGELC,0   
KEYW,PR_MULTI,0
KEYW,PR_CFD,0   
/GO
!*  
/COM,   
/COM,Preferences for GUI filtering have been set to display:
/COM,  Structural   
!*  
/PREP7  
!*  
ET,1,SOLID186   
!*  
!*  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,7.8E05   
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,0.8E3  
MPDATA,PRXY,1,,0.3  
SMRT,6  
SMRT,7  
SMRT,7  
SMRT,8  
MSHAPE,1,3D
MSHKEY,0
!*  
FLST,5,25,6,ORDE,2  
FITEM,5,1   
FITEM,5,-25
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
FINISH  
/SOL
!*  
ANTYPE,0
/ANG,1,-30,ZS,1
/REP,FAST   
/ANG,1,-30,ZS,1
/REP,FAST   
/ANG,1,-30,ZS,1
/REP,FAST   
/ANG,1,-30,ZS,1
/REP,FAST   
/ANG,1,-30,ZS,1
/REP,FAST   
/ANG,1,-30,ZS,1
/REP,FAST   
/ANG,1,-30,ZS,1
/REP,FAST   
/ANG,1,30,YS,1  
/REP,FAST   
/ANG,1,30,YS,1  
/REP,FAST   
/ANG,1,30,YS,1  
/REP,FAST   
dA,135,ALL  
DA,27,ALL   
/ANG,1,30,YS,1  
/REP,FAST   
/ANG,1,30,YS,1  
/REP,FAST   
/ANG,1,30,YS,1  
/REP,FAST   
/ANG,1,30,YS,1  
/REP,FAST   
DLIST, ALL  
SFLIS, ALL  
/VIEW,1,1,2,3   
/ANG,1  
/REP,FAST   
FINISH  
! /EXIT,NOSAV

jwooh 发表于 2006-11-14 09:19

FreddyMusic ,谢谢啦,欢迎多挑一些啊,我打算赶快把这个静力学的弄好之后,做点动力学的东西~~

FreddyMusic 发表于 2006-11-14 09:33

不用急着做动力学,还远着呢!

1. 内圈应力过大,好像不太可能?
    帖张你施加载荷的图,描述一下。

2. 下一步搞,第10章径向积分与轴向积分,你都看了?

页: 1 [2] 3
 

Powered by Discuz! Archiver 6.1.0  © 2001-2007 Comsenz Inc.