求助 ANSYS 两种对称加载的轴孔过盈结果不同？

malong

最近做轴孔过盈分析时，遇到一些困惑的问题，两种对称方式产生的接触对问题，后一种接触居然产生两个接触对，造成计算接触力的时候，结果相差2倍，希望高手指点一下，以下为一个例子，在两种不同加载下的结果     


期待高手能帮我解决一下，谢谢


Finish
/clear,start
/TITLE,Contact analysis with initial interference
/PREP7
!带孔圆盘的基本尺寸；
Rpin=35
Rpout=100
Hp=25
!轴的基本尺寸；
Rain=25
Raout=35
La=150
! 过盈量f;
f=0.01
!实体的单元类型为带中间节点的2阶六面体单元；
ET,1,solid186
MP,EX,1,2.1E5 !弹性模量；
MP,PRXY,1,0.3 !poisson系数；
!生成带孔圆盘的1/4实体模型；
CYL4,0,0,Raout,0,Rpout,90,Hp
!轴的1/4实体模型；
CYL4,0,0,Rain,0,Raout,90,La
!把轴的位置沿着轴向移动一段距离；
VGEN, ,2, , , , ,-10, , ,1
!
!
!
!************对实体划分网格 ****************************
!*
LESIZE,17, , ,15, , , , ,1
LESIZE,19, , ,15, , , , ,1
!*
LESIZE,18, , ,2, , , , ,1
LESIZE,20, , ,2, , , , ,1
!*
LESIZE,22, , ,20, , , , ,1
!*
LESIZE,5, , ,10, , , , ,1
LESIZE,7, , ,10, , , , ,1
!*
LESIZE,6, , ,8, , , , ,1
LESIZE,8, , ,8, , , , ,1
!*
LESIZE,10, , ,3, , , , ,1
!*
VSWEEP,ALL
!*************网格划分完毕**********************************
!
!
/COM, CONTACT PAIR CREATION - START
MP,MU,1,0.2
MAT,1
R,3
REAL,3
ET,2,170
ET,3,174
R,3,,,0.2,0.2,0.9,0 !Icont系数设置为0.2;
!设置过盈量为f;
RMORE,,,1.0E20,f,1.0,0
RMORE,0.0,0,1.0,,1.0,0
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,3,4,0
KEYOPT,3,5,0
NROPT,UNSYM
KEYOPT,3,7,0
KEYOPT,3,8,0
!KEYOPT(9)=4;
!不考虑初始几何位置造成的过盈或者间隙，
!只考虑CNOF参数设置的值，即过盈量。
KEYOPT,3,9,4
KEYOPT,3,10,2
KEYOPT,3,11,0
KEYOPT,3,12,0
KEYOPT,3,2,0
KEYOPT,2,5,0
! Generate the target surface
ASEL,S,,,4
CM,_TARGET,AREA
TYPE,2
NSLA,S,1
ESLN,S,0
ESLL,U
ESEL,U,ENAME,,188,189
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
ASEL,S,,,9
CM,_CONTACT,AREA
TYPE,3
NSLA,S,1
ESLN,S,0
ESURF
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,2
ESEL,A,TYPE,,3
ESEL,R,REAL,,3
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
/COM, CONTACT PAIR CREATION - END
FINISH
/SOL
FLST,2,4,5,ORDE,4
FITEM,2,5
FITEM,2,-6
FITEM,2,11
FITEM,2,-12
DA,P51X,SYMM !设置对称约束；
FLST,2,1,5,ORDE,1
FITEM,2,3
/GO
DA,P51X,ALL,
!*
ANTYPE,0
NLGEOM,1 !必须打开大变形效应；
NSUBST,1,0,0
AUTOTS,0
TIME,100
/STATUS,SOLU
SOLVE
FINISH
/POST1
SET,1,LAST,1,
PLNSOL,s,eqv !查看单元的Von Mises stress;
esel,s,ename,,174
PLNSOL, CONT,PRES, 0,1.0 !在contact 单元上查看contact pressure;

！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！修改对称约束加载方式

Finish
/clear,start
/TITLE,Contact analysis with initial interference
/PREP7
!带孔圆盘的基本尺寸；
Rpin=35
Rpout=100
Hp=25
!轴的基本尺寸；
Rain=25
Raout=35
La=150
! 过盈量f;
f=0.01
!实体的单元类型为带中间节点的2阶六面体单元；
ET,1,solid186
MP,EX,1,2.1E5 !弹性模量；
MP,PRXY,1,0.3 !poisson系数；
!生成带孔圆盘的1/4实体模型；
CYL4,0,0,Raout,0,Rpout,90,Hp
!轴的1/4实体模型；
CYL4,0,0,Rain,0,Raout,90,La
!把轴的位置沿着轴向移动一段距离；
VGEN, ,2, , , , ,-10, , ,1
!
!
!
!************对实体划分网格 ****************************
!*
LESIZE,17, , ,15, , , , ,1
LESIZE,19, , ,15, , , , ,1
!*
LESIZE,18, , ,2, , , , ,1
LESIZE,20, , ,2, , , , ,1
!*
LESIZE,22, , ,20, , , , ,1
!*
LESIZE,5, , ,10, , , , ,1
LESIZE,7, , ,10, , , , ,1
!*
LESIZE,6, , ,8, , , , ,1
LESIZE,8, , ,8, , , , ,1
!*
LESIZE,10, , ,3, , , , ,1
!*
VSWEEP,ALL
!*************网格划分完毕**********************************
!
!
/COM, CONTACT PAIR CREATION - START
MP,MU,1,0.2
MAT,1
R,3
REAL,3
ET,2,170
ET,3,174
R,3,,,0.2,0.2,0.9,0 !Icont系数设置为0.2;
!设置过盈量为f;
RMORE,,,1.0E20,f,1.0,0
RMORE,0.0,0,1.0,,1.0,0
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,3,4,0
KEYOPT,3,5,0
NROPT,UNSYM
KEYOPT,3,7,0
KEYOPT,3,8,0
!KEYOPT(9)=4;
!不考虑初始几何位置造成的过盈或者间隙，
!只考虑CNOF参数设置的值，即过盈量。
KEYOPT,3,9,4
KEYOPT,3,10,2
KEYOPT,3,11,0
KEYOPT,3,12,0
KEYOPT,3,2,0
KEYOPT,2,5,0
! Generate the target surface
ASEL,S,,,4
CM,_TARGET,AREA
TYPE,2
NSLA,S,1
ESLN,S,0
ESLL,U
ESEL,U,ENAME,,188,189
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
ASEL,S,,,9
CM,_CONTACT,AREA
TYPE,3
NSLA,S,1
ESLN,S,0
ESURF
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,2
ESEL,A,TYPE,,3
ESEL,R,REAL,,3
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
/COM, CONTACT PAIR CREATION - END
FINISH
/SOL
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!原来的对称约束加载
!FLST,2,4,5,ORDE,4
!FITEM,2,5
!FITEM,2,-6
!FITEM,2,11
!FITEM,2,-12
!DA,P51X,SYMM !设置对称约束；
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!更改对称约束加载
/PREP7  
csys,1
CYCLIC, , , ,'CYCLIC'


FLST,2,1,5,ORDE,1
FITEM,2,3
/GO
DA,P51X,ALL,
!*
ANTYPE,0
NLGEOM,1 !必须打开大变形效应；
!NSUBST,1,0,0
!AUTOTS,0
!TIME,100
/STATUS,SOLU
SOLVE
FINISH
!/POST1
!SET,1,LAST,1,
!PLNSOL,s,eqv !查看单元的Von Mises stress;
!esel,s,ename,,174
!PLNSOL, CONT,PRES, 0,1.0 !在contact 单元上查看contact pressure;


正确结果
*** NOTE ***                            CP =     322.266   TIME= 18:04:49
Summations based on final geometry and will not agree with solution     
reactions.                                                              
***** SUMMATION OF TOTAL    FORCES AND MOMENTS IN COORDINATE SYSTEM  1
       NOTE: THE SUM IS DONE IN COORDINATE SYSTEM  1
  FX  =  -24151.83   
  FY  = -0.7195955E-04
  FZ  = -0.6043672E-01
  MX  = -0.3046919E-03
  MY  =  -302803.7   
  MZ  = -0.9499917E-02
SUMMATION POINT=  0.0000      0.0000      0.0000   



错误结果

*** NOTE ***                            CP =     299.312   TIME= 17:56:22
Summations based on final geometry and will not agree with solution     
reactions.                                                              
***** SUMMATION OF TOTAL    FORCES AND MOMENTS IN COORDINATE SYSTEM  1
       NOTE: THE SUM IS DONE IN COORDINATE SYSTEM  1
  FX  =  -48303.98   
  FY  =  0.4374504E-03
  FZ  =  0.5176381E-01
  MX  = -0.5855726E-02
  MY  =  -605616.9   
  MZ  =  0.4383126E-03
SUMMATION POINT=  0.0000      0.0000      0.0000

[ 本帖最后由 malong 于 2009-6-9 18:12 编辑 ]

gy9991

我是新手，虽然回答不了你的问题。但是咱两人做的是同一个问题，可否帮忙看看下面这个问题啊？
http://forum.vibunion.com/forum/thread-83588-1-1.html
我先谢谢啦！

WASPCN

这个命令有点问题 R,3,,,0.2,0.2,0.9,0 !Icont系数设置为0.2;

我认为应该是
R,3,,,0.1,0.1,0.2,0 !Icont系数设置为0.2;