basic contact problem for a beam model

[dovil]

Dear all,

I have a problem modeling a contact beam models.
I am a beginning of the contact models.
I attached a very simple model made by myself.
Can anybody help me to check what is wrong with my model?

Thank you very much for your time.

Best Regards,
Dovil




FINISH
/CLEAR,START
*AFUN,RAD
/PREP7

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! Define Parameters for Geometry !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
L = 15
Ri = 10e-3/2
Ro = 15e-3/2
h = 2.1353
pi = 3.1415926
NumberOfSegments = 10
F = 5E-11
NN = 1 ! number of waves

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! DEFINE MATERIAL PROPERTIES !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
MP,EX,1,100e-3
MP,NUXY,1,0.3

!!!!!!!!!!!!!!!!!!!!!!
!DEFINE ELEMENT TYPES!
!!!!!!!!!!!!!!!!!!!!!!
ET,2,beam188

!!!!!!!!!!!!!!!!
!Build Geometry!
!!!!!!!!!!!!!!!!
*do,i,0,NumberOfSegments,1
*set,x,i/NumberOfSegments*L
*set,y,h*sin(2*pi*x/L)
k,i+1,x,y
*enddo

*do,j,1,NumberOfSegments,1
l,j,j+1
*enddo

LGEN,NN,all, , , L,, , ,0
LSEL,ALL
lglue,all

*get,max_kp_num,kp,,num,max
*do,i,0,NumberOfSegments,1
*set,x,i/NumberOfSegments*L
*set,y,-h*sin(2*pi*x/L) - 2.2*h
k,i+max_kp_num+1,x,y
*enddo

*do,j,1,NumberOfSegments,1
l,j+max_kp_num,j+1+max_kp_num
*enddo

lsel,s,line,,max_kp_num,2*(max_kp_num-1)
LGEN,NN,all, , , L,, , ,0
lglue,all

!l,11,22

lsel,all
lglue,all
LESIZE,all,,,10,,,,,

SECTYPE, 1, BEAM, CTUBE, , 0
SECOFFSET, CENT
SECDATA,Ri,Ro,0,0,0,0,0,0,0,0,0,0

LATT,1, ,2, , , ,1

lsel,all
lmesh,all



!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!R,1000,,,Fkn,Ftol,Icon,,
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
ET,11,targe170
TSHAP,line
KEYOPT,11,1,0
!KEYOPT,11,5,0


ET,12,conta176
KEYOPT,12,5,4
R,10,Ro,Ro

nsel,s,loc,y,-h - 0.0001, h + 0.0001
nsel,u,loc,x,L
cm,node_1,node

nsel,all
cmsel,u,node_1
nsel,u,loc,x,L
cm,node_2,node

allsel,
cmsel,s,node_2
real,10
type,12
esurf,,top,

allsel,
cmsel,s,node_1
!lsel,all
!lsel,u,loc,y,-h - 0.0001, h + 0.0001
real,10
type,11
esurf,,,

!!!!!!!!!!!!!!!!!!!!
!Solve the problem !
!!!!!!!!!!!!!!!!!!!!
ALLSEL
FINISH

/solu

!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!Apply boundary conditions !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!
ksel,s,kp,,1
ksel,a,kp,,(NumberOfSegments+1)*NN + 1
DK,all,all,0

!!!!!!!!!!!!!
!Apply load !
!!!!!!!!!!!!!
!ksel,all
!FK,(NumberOfSegments+1)*NN,FX,-F
!FK,(NumberOfSegments+1)*NN*2,FX,-F
ksel,all
DK,(NumberOfSegments+1)*NN,UX,-3
DK,(NumberOfSegments+1)*NN*2,UX,-3

allsel
ANTYPE,0

!NLGEOM,on
NSUBST,100,1000,1
SOLVE
FINISH
!!!!!!!!!!!!!!
!displace DX !
!!!!!!!!!!!!!!

/POST1
/DSCALE,ALL,1.0
/EFACET,1
PLNSOL, U,X, 1,1.0