马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。
您需要 登录 才可以下载或查看，没有账号？我要加入
x
最近做轴孔过盈分析时，遇到一些困惑的问题，两种对称方式产生的接触对问题，后一种接触居然产生两个接触对，造成计算接触力的时候，结果相差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 编辑 ] |
发表于 2009-6-9 18:06
提升卡
变色卡
显身卡