Problem with 3D Contact Analysis 1

[jujumoi]

Hello,
I'am having some troubles trying to set up a 3D contact analysis of a cylinder and an arc. The contact is not occuring, the contact volume is not interacting with target elements. I think Iam probably not constraining the model correctly. The input file is listed below, if anyone can give me a hint I would really appreciate.

/prep7
d=1
r1=0.8
t=0.1
mf=20
et,1,185
mp,ex,1,30e6
mp,prxy,1,.3
Cyl4,0,0,r1,90,,,1
wpof,,-r1*9
Cyl4,0,0,r1*10,90,r1*10+t,,1
lsel,s,,,11,13,2
lsel,a,,,14,16,2
lesize,all,,,mf,1/mf
lsel,s,,,12,15,3
lesize,all,,,mf/5
lsel,s,,,7
lsel,a,,,18,20,2
lesize,all,,,mf/5*(d/t)
allsel,all
MSHAPE,0,3d
MSHKEY,1
VMESH,2
lsel,s,,,3,4,
lesize,all,,,mf,1/mf
lsel,a,,,2,5,
lesize,all,,,mf,mf
lplot
MSHAPE,0,3d
VMESH,1
allsel,all
et,2,170
et,3,173
keyopt,3,9,1
keyopt,3,10,1
r,3,0.0,0.0,0.5,0.05,0.0,0.0
rmore,0.0,0.0,0.1E+21,0.0,1.0,0.0
rmore,0.0,0.0,1.0,0.0,1.0,0.50,
type,3
real,3
asel,s,,,3,
nsla,s,1
esln,s
eplot
esurf
allsel,all
asel,s,,,9,
nsla,s,1
esln,s
eplot
esurf
allsel,all

/solu
nlgeom,on
nsubst,1
asel,s,loc,x,0
asel,a,loc,z,0
asel,a,loc,z,d
da,all,symm
allsel
asel,s,loc,y,0
asel,a,,,10,
da,all,uy,0
asel,s,,,8,
da,all,uy,-.02
allsel
/STATUS,SOLU
SOLVE
/STATUS,SOLU
FINISH
/POST1
!*
/EFACET,1
PLNSOL, S,EQV, 0,1.0

Thanks.

 

[doatis]

Just a quick look, and run on my machine, you don't have any target elements defined. You defined the element type, but never created them in the active session. Here's the first half of your script with a single line added:

/prep7
d=1
r1=0.8
t=0.1
mf=20
et,1,185
mp,ex,1,30e6
mp,prxy,1,.3
Cyl4,0,0,r1,90,,,1
wpof,,-r1*9
Cyl4,0,0,r1*10,90,r1*10+t,,1
lsel,s,,,11,13,2
lsel,a,,,14,16,2
lesize,all,,,mf,1/mf
lsel,s,,,12,15,3
lesize,all,,,mf/5
lsel,s,,,7
lsel,a,,,18,20,2
lesize,all,,,mf/5*(d/t)
allsel,all
MSHAPE,0,3d
MSHKEY,1
VMESH,2
lsel,s,,,3,4,
lesize,all,,,mf,1/mf
lsel,a,,,2,5,
lesize,all,,,mf,mf
lplot
MSHAPE,0,3d
VMESH,1
allsel,all
et,2,170
et,3,173
keyopt,3,9,1
keyopt,3,10,1
r,3,0.0,0.0,0.5,0.05,0.0,0.0
rmore,0.0,0.0,0.1E+21,0.0,1.0,0.0
rmore,0.0,0.0,1.0,0.0,1.0,0.50,
type,3
real,3
asel,s,,,3,
nsla,s,1
esln,s
eplot
esurf
type,2 !Switch to target elements
allsel,all
asel,s,,,9,
nsla,s,1
esln,s
eplot
esurf
allsel,all

This at least defines the contact surface for me, and it's initial condition is closed.

Hope this helps,
Doug

 

[doatis]

Sorry, just another follow-up, your problem solved on my machine after ~10 minutes.

Good luck,
Doug

 

[jujumoi]

Thak you very much Doug,
I was killing myself here, I though I`ve already defined the target elements on the active CS, and I couldn`t find this little simple detail.
Thanks again!