作了一个接触的例子，有点问题，帮忙看看
我做了一个接触分析，但是在顶端施加载荷时，结果感觉球与压板好像独立似的，各位师兄帮忙看看Re:作了一个接触的例子，有点问题，帮忙看看
附上命令流，希望大家多帮忙：finish
/clear
fact=1000
rad=fact*2.5e-6
th=fact*0.125e-6
X1=fact*13.5e-6
X2=fact*20e-6
Y1=fact*20e-6
Y2=fact*40e-6
/prep7
et,1,plane182
et,2,contac48 !设置二维点面接触单元
KEYOPT,1,3,1 !设置为轴对称问题
mp,ex,1,128000
mp,nuxy,1,0.3
MP,ALPX,1,2.7e-6
mp,ex,2,70000
mp,nuxy,2,0.2
MP,ALPX,2,3.5e-6
MPTEMP,,,,,,,,
MPTEMP,,20,100,190
MPDATA,EX,3,,227500,221300,214400
MPDATA,PRXY,3,,0.3,0.32,0.38
MPTEMP,,,,,,,,
MPTEMP,,-40,100,190
MPDATA,ALPX,3,,12.5e-6,14.5e-6,15e-6
MP,REFT,3,190
MPTEMP,,,,,,,,
MPTEMP,,20,100,150,190
MPDATA,EX,4,,1600,800,200,50
MPDATA,PRXY,4,,0.4,0.45,0.47,0.49
MPTEMP,,,,,,,,
MPTEMP,,-40,130,137,190
MPDATA,ALPX,4,,70e-6,70e-6,145e-6,150e-6,
MP,REFT,4,190
MPTEMP,,,,,,,,
MPTEMP,,20,100,150,190
MPDATA,EX,5,,2000,1200,400,100
MPDATA,PRXY,5,,0.4,0.45,0.47,0.485
MPTEMP,,,,,,,,
MPTEMP,,-40,130,137,190
MPDATA,ALPX,5,,70e-6,70e-6,145e-6,150e-6,
MP,REFT,5,190
MPTEMP,,,,,,,,
MPTEMP,,20,100,190
MPDATA,EX,6,,227500,221300,214400
MPDATA,PRXY,6,,0.3,0.32,0.38
MPTEMP,,,,,,,,
MPTEMP,,-40,100,190
MPDATA,ALPX,6,,12.5e-6,14.5e-6,15e-6
MP,REFT,6,190
r,1,227500,,,,0.001
r,2,227500,,,,0.001
!建立模型
csys,1
k,1,rad-th,-90
k,2,rad,-90
kgen,19,1,2,1,,10
*do,j,1,36,2
a,j,j+1,j+3,j+2
*enddo
rectng,,X1,rad,rad+Y1
rectng,,X2,rad+Y1,rad+Y2
rectng,,X2,-rad,-rad-Y1
k, !该关键点编号为51
csys,1
a,37,19,51
a,1,19,51
asel,s,loc,x,0,rad
aglue,all
asel,s,area,,19,20,1
aglue,all
!划分网格
csys,1
lsel,s,loc,x,rad-th,rad
lesize,all,,,4
mat,6
type,1
asel,s,area,,1,18,1
amesh,all
lsel,s,loc,x,rad-th
lsel,r,loc,y,0,90
lccat,all
lsel,s,loc,x,rad-th
lsel,r,loc,y,0,-90
lccat,all
mat,5
type,1
mshape,0,2d
asel,s,area,,24,25,1
mshkey,1
amesh,all
lsel,s,line,,56
lesize,all,,,30,8
lsel,s,line,,58
lesize,all,,,30,0.125
lsel,s,line,,57,59,2
lesize,all,,,20
mat,3
type,1
asel,s,area,,19
mshape,0,2d
amesh,all
csys,0
lsel,s,line,,61,71,10
lesize,all,,,10
lsel,s,loc,y,rad+Y1
lccat,all
mat,1
type,1
asel,s,area,,22
mshape,0,2d
amesh,all
lsel,s,line,,64
lesize,all,,,30,8
lsel,s,line,,66
lesize,all,,,30,0.125
lsel,s,line,,65,67,2
lesize,all,,,30
mat,2
type,1
asel,s,area,,21
mshape,0,2d
amesh,all
save
!创建接触单元
lsel,s,line,,56
nsll,,1
nplot
cm,target1,node !定义目标节点
csys,1
nsel,s,loc,y,60,90
nsel,r,loc,x,rad
nsel,u,loc,y,90
nsel,u,node,,349
nsel,u,node,,2312
nplot
cm,contac1,node !定义接触节点
type,2
real,1
nsel,all
gcgen,contac1,target1
gcgen,target1,contac1
lsel,s,line,,66
nsll,,1
nplot
cm,target2,node !定义目标节点
csys,1
nsel,s,loc,y,-60,-90
nsel,r,loc,x,rad
nsel,u,loc,y,-90
nsel,u,node,,7
nsel,u,node,,3394
nplot
cm,contac2,node !定义接触节点
type,2
real,2
nsel,all
gcgen,contac2,target2
gcgen,target2,contac2
save
finish
/solu
antype,static
!施加约束
csys,0
nsel,s,loc,x,0
DSYM,SYMM,X
nsel,s,loc,y,-rad-Y1
d,all,all
csys,1
nsel,s,loc,y,-90
nsel,r,loc,x,rad
d,all,all
!设定非线性选项
outres,all,all
time,1
nsubst,200
!nlgeom,on
lnsrch,on
!nropt,full
lsel,s,line,,62
nsll,,1
SF,all,PRES,0.15
allsel
solve :)
页:
[1]