作了个简单的面面接触

freddieaid

装配问题,两个体最初在z方向有5个单位的距离,在z方向给小圆柱-20单位的强制位移,使其压入配件.见下图,采用对称结构

freddieaid

最终压力

freddieaid

最终摩擦力

freddieaid

发现一个有趣的问题,当初始接触面和目标面分离较大时,采用拉格朗日算法,节点搜索,能够很快收敛;但是罚函数却不能收敛,奇怪!

freddieaid

搞定,原来是罚函数的罚参数取得过大

freddieaid

附上命令
/prep7
et,1,solid45
mp,ex,1,1e10
mp,nuxy,1,.2
et,2,conta173
et,3,targe170
mp,mu,2,0.2
et,4,mesh200
keyopt,4,1,6
keyopt,2,9,0
!keyopt,2,2,4
keyopt,2,4,2
r,1,,,.1,0.1     !定义法向刚度因子为1

!建模
/view,,1,1,1
k,1,0,0,0
k,2,0,10,0
k,3,0,50,0
k,4,10,0
k,5,50,0,0
larc,2,4,1,10
larc,3,5,1,50
l,2,3
l,4,5
al,1,2,3,4
!用mesh200划分面
type,4
lesize,all,,,8
mshape,0,2d
mshkey,1
amesh,all
!拉伸生成体,共拉伸80
esize,,12
type,1
mat,1
vext,all,,,0,0,80
!生成小柱
*get,kmax,kp,,num,max
*get,lmax,line,,num,max
lsel,none
k,kmax+1,0,0,200
k,kmax+2,11,0,200
k,kmax+3,0,11,200
l,kmax+1,kmax+2
l,kmax+1,kmax+3
larc,kmax+2,kmax+3,kmax+1,11
asel,none
al,lmax+1,lmax+2,lmax+3
!用mesh200划分面
type,4
lesize,all,,,6
mshape,0,2d
mshkey,1
amesh,7
!拉伸生成体,共拉伸100
esize,,16
type,1
mat,1
asel,s,,,7
vext,all,,,0,0,-120
allsel,all

!生成目标面
type,3
mat,2
real,1
tshap,quad
asel,s,,,3
cm,_target,area
nsla,s,1
!esln,s,0
esurf
!生成接触面
type,2
mat,2
real,1
tshap,quad
asel,s,,,10
cm,_contact,area
nsla,s,1
!esln,s,0
esurf
allsel,all
/expand,4,polar,half,,90
/view,,1,1,1
eplot,all
fini
/solu
allsel,all
antype,static,new
outpr,all,1
outres,all,1
!NROPT,SYMm
nlgeom,on
solcontrol,on
rescontrol,define,all,-1
autots,off
nsubst,20

!施加对称边界条件
asel,s,,,4,6,2
da,all,symm
asel,s,,,9,11,2
da,all,symm
!nsel,s,loc,z,0
!d,all,all
!allsel,all
da,5,all,0

da,7,ux,0
da,7,uy,0
da,7,uz,-10
allsel,all
eplot,all
solve

zhangfl

我在一本书上看到一模一样的,但不知道如何装配的?能不能调用外部模型来装配?

freddieaid

命令一样吗,纯属原创,见笑

zhangfl

命令看不懂,但是内容相似,我一直想把它弄明白,无奈新手上路,道路坎坷,只能感觉一样---一轴一柱,如果不是对称实体,能否按你所说的方法分析.

freddieaid

下面引用由zhangfl在 2005/05/11 03:47pm 发表的内容：
命令看不懂,但是内容相似,我一直想把它弄明白,无奈新手上路,道路坎坷,只能感觉一样---一轴一柱,如果不是对称实体,能否按你所说的方法分析.

不是对称实体,解题思路还是一样的

zhourunfa66

{:sleepy:}{:sleepy:}
楼主问题解决了么？？
解决了的话教程能否分享下？？？
谢谢
381848450@qq.com