Плоский контакте в Ansys

[Цветочек]

Здравствуйте!
Недавно, взялась за ансис и кажется есть небольшие успехи - но все омрачает построение контактных пар.
Способов для их создания очень много и самые простые похоже применяют реже всего. У меня такой случай, написала простенькую программу для расчета контакта резиновой детали с абсолютно жесткой целью. К сожалению, контакта так и не наблюдается - даже не знаю - что и делать.
Если Вам эта тема знакома, посмотрите пожалуйста где я ошиблась.

/prep7
*SET,R,0.15
*SET,ar,0.1480
smrt,off
ET,1,PLANE182, , ,2
KEYOPT,1,6,1
keyopt,2,5,3
ET,2,CONTA171 !
R,2,,,-2000
ET,3,TARGE169
MP,EX,1,2.82
C10 = 0.293
C01 = 0.177
NU1 = 0.49967
DD = 2*(1-2*NU1)/(C10+C01)
TB,HYPER,1,1,2,MOONEY
TBDATA,1,C10,C01,DD
CSYS,1
k,1
K,2,R,-80
K,3,R,90
K,4,R,180
K,5,R,260
K,6,ar,-80
K,7,ar,90
K,8,ar,180
K,9,ar,260
k,10,R+ar,270
L,2,5
L,6,9
L,2,6
L,5,9
AL,1,3,2,4
type,1
ESIZE,,20
AMESH,ALL
local,11,1,0,-R-ar,0
n,1003,ar,80
n,1002,ar,100
n,1001,0,0,0
real,2
type,3
tshap,arc
E,1002,1003,1001
alls
real,2
type,2
lsel,s,line,,1
nsll,s,1
esurf
ALLSEL,ALL
FINISH

/SOLU
lsel,s,line,,2
nsll,s,1
D,all,UX,0
D,all,UY,-0.0005
alls
NSEL,S,NODE,,1001,1003
D,all,all
ANTYPE,STATIC
CNVTOL,F,,,,-1
NLGEOM,ON
NSUBST,600
OUTRES,,1
SOLVE
FINISH

Буду очень признательна

[Zai]

Посмотрите сообщение
http://dxdy.ru/topic2066-15.html

[Цветочек]

Спасибо вам большое за внимание!
Только все конечно в том коде по ссылке верно - единственное, с чем у меня проблема, код который здесь представлен не работает если за целевой элемент выбрать дугу или окружность, а вот если целевой элемент - линия - все работает прекрасно!
Если в том коде, на который вы указали - поверхность target будет не плоской - окружность, например, как будет выглядеть код? С линией у меня все работает - а я хотела бы чтобы работало и с дугами, окружностями!
Я приведу пример своего кода, если целевая поверхность - линия:
/prep7
smrt,off
ET,1,PLANE182, , ,2
KEYOPT,1,6,1
ET,2,CONTA175
R,2, , ,-2000
ET,3,TARGE169
MP,EX,1,2.82
MP,NUXY,1,0.49967
C10 = 0.293
C01 = 0.177
NU1 = 0.49967
DD = 2*(1-2*NU1)/(C10+C01)
TB,HYPER,1,1,2,MOONEY
TBDATA,1,C10,C01,DD
CSYS,1
*SET,R,0.15
*SET,ar,0.1480
k,1
K,2,R,-80
K,3,R,90
K,4,R,180
K,5,R,260
K,6,ar,-80
K,7,ar,90
K,8,ar,180
K,9,ar,260
k,10,R+ar,270
L,2,5
L,6,9
L,2,6
L,5,9
AL,1,3,2,4
DL,2,1,UX,0
DL,2,1,UY,-0.00025
ESIZE,,100
AMESH,ALL !
SAVE
csys,0
N,9001,(-2*R),-R
N,9002,(2*R) ,-R
NSEL,S,NODE,,9001,9002
TYPE,3
REAL,2
TSHAP,LINE
E,9002,9001
D,ALL,ALL,0
lsel,s,line,,1
nsll,s,1
MAT,2
TYPE,2
REAL,2
esurf
ALLSEL,ALL
FINISH

/SOLU
ANTYPE,STATIC
CNVTOL,F,,,,-1
NLGEOM,ON
NSUBST,600
OUTRES,,1
SOLVE
FINISH

[Zai]

В Вашем первом примере
Замените conta171 на conta172.
При создании дуги нужно проверять ее нормаль. Элементы созданные без линий лучше задавать следующей командой(чтобы знать его номер)
...
tshap,arc
En,10000,1002,1003,1001
ESEL,S,,,10000
ESURF,,REVERSE ! Очень важно - нормаль меняет направление

alls
...

В Вашем втором примере
Замените conta175 на conta172. Ваше решение что-то плохо сходится. Для начала проверьте решение на меньшем смещении (DL,2,1,UY,-0.0002)

[Цветочек]

Спасибо за помощь, все поняла и наконец-то итерации сошлись!
Благодарю Вас от чистого сердца!
Цветочек