[Цветочек]

1. Можно ли считать контакт bonded закреплением?
Я в общем-то считаю что можно, но поясню - решаю задачу сжатия двумя жесткими плоскостями(targe170) стального пластичного диска в трехмерной постановке. Так вот нагрузку прикладываю к пилотным узлам,
d,30001,all,0
d,30005,uz,-0.5
d,30005,ux,0
d,30005,uy,0
Поскольку использую keyopt,all,12,3 кажется что дополнительные граничные условия не нужны - права ли я?
2. Какие могут быть дефекты контакта если контакт изначально открыт - хотя детали рядом, и сколько не переставляю нормали - всё так и есть?

[aka raduntsev-nias]

Для связных контакт дб MPC или Penalty с bonded
пусть поправят , если не прав

[derik]

а код мона посмотреть или это военная тайна?

http://forum.dwg.ru/showthread.php?t=52057&page=9

тут почитайте...

[Цветочек]

Вот код по которому я спрашиваю

Код:

[Выделить все]

/prep7
ET,1,SOLID185   
!*  
KEYOPT,1,2,0
KEYOPT,1,6,1
KEYOPT,1,10,0   
!*  
!*  
MPTEMP,1,0                             
MPDATA,EX,1,,2.1E11            
MPDATA,PRXY,1,,0.3  
TB,BKIN,1,1,2,0
TBTEMP,0  
TBDATA,1,2.6E8,0,,,,  

real,1
CYL4,0,0,17.5, , , ,8   
vsweep,all


! Создание контактной пары на нижней площадке
ET,2,CONTA173 ! можно делать conta173
keyopt,2,12,3 ! контакт bonded
KEYOPT,2,5,4
ET,3,TARGE170 ! 3-D TARGET ELEMENT
KEYOPT,3,2,1
! Геометрия твердой поверхности = target
n,30001,30,-30
n,30002,30,30
n,30003,-30,30
n,30004,-30,-30
! Создание площадки
alls
real,2
type,3
tshap,quad
E,30001,30002,30003,30004 ! плоскость площадки элемент
! Сама поставит нормаль
tshape,pilo
EN,6001,30001 ! Создаётся элемент 6001 pilo для контроля
allsel,all
type,2
NSEL,S,LOC,z,0 ! Выделяю плоскость c z=0;
NSEL,U,NODE,,30001,30004
ESURF,, ! GENERATE
allsel,all
! Эта пара начально закрыта - всё хорошо

real,3
! Создание контактной пары на инденторе
ET,4,CONTA173 
KEYOPT,4,5,4
keyopt,4,12,3 ! контакт bonded
ET,5,TARGE170 
KEYOPT,5,2,1
! Геометрия твердой поверхности = target
n,30005,30,-30,8
n,30006,30,30,8
n,30007,-30,30,8
n,30008,-30,-30,8
! Создание площадки
alls
type,5
tshap,quad
E,30005,30006,30007,30008 ! плоскость площадки элемент
tshape,pilo
EN,7001,30005 ! Создаётся элемент 7001 pilo для контроля

allsel,all
type,4
NSEL,S,LOC,z,8 ! Выделяю плоскость c z=8;
NSEL,U,NODE,,30005,30008
esurf,,
allsel,all
fini
/solu
nlgeom,on
lnsrch,on
d,30001,all,0
d,30005,uz,-0.5
d,30005,ux,0
d,30005,uy,0
solve

тупой Цветочек седня!

[derik]

вроде как криманала нет...кейопт,3,2,1 стоит, пилотные узлы закреплены...Сергей Иванович скажет больше, у него опыта больше в подобных делах.
Цветочек не тупой, цветочек, вероятно всего, просто устал! ;-)

[Цветочек]

Дело в том, что вот эта ошибка - вообще преследует меня - часто пишу в ансис с двумя контактными парами и очень часто эта ошибка - первая пара всегда закрыта, а вот вторая! всегда открыта - неприятно что-то недопонимать, но видимо в данном случае это так: реальные константы для разных пар разные, определение quad и pilo, как обычно; контактор можно покрывать и esurf и amesh, тоже как-бы обычно; если меняю направление - esurf,,reve - контактная пара вообще не создаётся. Так что видимо дело не в направлении нормалей, но в чём? - неужели я правильно грешу на граничные условия?

[aka raduntsev-nias]

в WB в случае открытого контакта использовал pinball
Но, насколько понимаю, открытый есть не всегда плохо

Yellow: "The contact status is open. This may be acceptable."
Workbench has detected an open contact Status condition on a nonlinear contact type, Frictionless, Rough, or Frictional, which is probably acceptable under certain conditions as stated in their descriptions.
If the Status is Far Open, the Penetration and the Gap will be set to zero even though the Resulting Pinball is non-zero.

в APDL
Use real constant PINB to specify a scaling factor (positive value for PINB when using command input) or absolute value (negative value for PINB when using command input) for the pinball region. You can specify PINB to have any value. By default, and assuming that large deflection effects apply (NLGEOM,ON), ANSYS defines the pinball region as a circle for 2-D or a sphere for 3-D of radius 4*depth (if rigid-to-flexible contact) or 2*depth (if flexible-to-flexible contact) of the underlying element. (See the discussion of element depth in Positive and Negative Real Constant Values.) If you include no large-deflection effects (NLGEOM,OFF), the default pinball region is half that of the large-deflection case. (For the no-separation (KEYOPT(12) = 4) and bonded-always (KEYOPT(12) = 5 options, the PINB default is different than described here. See Selecting Surface Interaction Models for more information.)

[derik]

а ты попробуй поставить keyopt,,5,1 - для автоматического закрытия зазора.5,4 - закрытие зазора и уменьшение проникновения.

[aka raduntsev-nias]

только что закрыл такую же ситуацию функцией нащупывания в WB
auto CNOF adjustment to close geometric gap only = keyopt,,5,1
keyo,cid,5,1 в wb

[derik]

вот а я про што! а вы мне не верите

[aka raduntsev-nias]

лучше разобраться с контактами можно по файлику, что я скинул у тебя в теме
подробно расписано

[Цветочек]

Нет, автоматическое закрытие зазора не то - не помогло - контакт открыт и контактная пара вне контакта!

[derik]

а выходной файлик мона посмотреть? (перед прер7 пропиши - /outp,outfile,txt,C:\ansys\output или в любую другую директорию

[Цветочек]

Поясню позицию иначе, когда запускается решение можно лицезреть две параллельные линии - норма и решение, вообще без пересечений

cncheck первой пары


*** NOTE *** CP = 6528.000 TIME= 19:15:20
Internal nodes from 30521 to 30691 are created.
1 internal nodes are used for handling DOF on pilot nodes of rigid
target surfaces.
170 internal nodes are used for 18x elements with mixed formulation.

*** NOTE *** CP = 6528.031 TIME= 19:15:20
Internal nodes from 30521 to 30691 are created.
1 internal nodes are used for handling DOF on pilot nodes of rigid
target surfaces.
170 internal nodes are used for 18x elements with mixed formulation.

*** NOTE *** CP = 6528.047 TIME= 19:15:20
Rigid-deformable contact pair identified by real constant set 2 and
contact element type 2 has been set up. The DOF of rigid surface is
driven by pilot node 30001.
Please verify constraints on the pilot node which may be
automatically fixed by ANSYS.
Contact algorithm: Augmented Lagrange method
Contact detection at: Gauss integration point
Default contact stiffness factor FKN 1.0000
The resulting contact stiffness 0.12571E+11
Default penetration tolerance factor FTOLN 0.10000
The resulting penetration tolerance 0.33410
Default opening contact stiffness OPSF will be used.
Default tangent stiffness factor FKT 1.0000
Default Max. friction stress TAUMAX 0.10000E+21
Average contact surface length 2.3968
Average contact pair depth 3.3410
Default pinball region factor PINB 3.0000
The resulting pinball region 10.023
Auto contact offset used to close gap 0.0000
*WARNING*: Initial penetration is included.
Bonded contact is defined.

*** NOTE *** CP = 6528.047 TIME= 19:15:20
Max. Initial penetration 0 was detected between contact element 6002
and target element 511.
****************************************


1 CONTACT PAIR IS SELECTED
CONTACT PAIR HAVING REAL ID = 2 IS INITIALLY CLOSED

cncheck второй - проблемной пары


NOTE *** CP = 6528.500 TIME= 19:16:00
Internal nodes from 30521 to 30691 are created.
1 internal nodes are used for handling DOF on pilot nodes of rigid
target surfaces.
170 internal nodes are used for 18x elements with mixed formulation.

*** NOTE *** CP = 6528.547 TIME= 19:16:00
Internal nodes from 30521 to 30691 are created.
1 internal nodes are used for handling DOF on pilot nodes of rigid
target surfaces.
170 internal nodes are used for 18x elements with mixed formulation.

*** NOTE *** CP = 6528.562 TIME= 19:16:00
Rigid-deformable contact pair identified by real constant set 3 and
contact element type 4 has been set up. The DOF of rigid surface is
driven by pilot node 30005.
Please verify constraints on the pilot node which may be
automatically fixed by ANSYS.
Contact algorithm: Augmented Lagrange method
Contact detection at: Gauss integration point
Default contact stiffness factor FKN 1.0000
The resulting contact stiffness 0.18029E+11
Default penetration tolerance factor FTOLN 0.10000
The resulting penetration tolerance 0.23295
Default opening contact stiffness OPSF will be used.
Default tangent stiffness factor FKT 1.0000
Default Max. friction stress TAUMAX 0.10000E+21
Average contact surface length 2.3968
Average contact pair depth 2.3295
Default pinball region factor PINB 3.0000
The resulting pinball region 6.9886
Auto contact offset used to close gap 0.0000
*WARNING*: Initial penetration is included.
Bonded contact is defined.

*** NOTE *** CP = 6528.562 TIME= 19:16:00
No contact was detected for this contact pair.
****************************************


*** WARNING *** CP = 6528.562 TIME= 19:16:00
All selected contact pairs are initially open. Rigid body motion can
occur. You may use auto CNOF/ICONT by setting KEYOPT(5) to close
small gaps.

1 CONTACT PAIR IS SELECTED
CONTACT PAIR HAVING REAL ID = 3 IS INITIALLY OPEN
BUT MAY BE AUTOMATICALLY CLOSED

[aka raduntsev-nias]

Contact algorithm: Augmented Lagrange method
Contact detection at: Gauss integration point
Default contact stiffness factor FKN 1.0000

для bonded penalty или MPC
для Augmented Lagrange method FKN 0,1
для bonded 1-10 можно поставить на проверку в интерациях если есть сомнение

и PINB бы не помешал, вроде у вас он по умолчанию. контакт точно закроется с ф-ией автоматического закрытия зазора

вообще. нафиг бондеды - связную сетку и все

Сергей Иванович в теме брачо писал что делать, когда нет сходимости
геометрию дадите (3d), прогоню через WB. быстрей найдем ошибки

[Цветочек]

Цилиндр 35 радиусом и 8 толщиной - блин, с двух сторон сжимается 2 плоскостями - вот и геометрия.

[aka raduntsev-nias]

точно) открыл ваш файл.
решил, но видно, где то не то задал.
у меня пока это получилось

[Цветочек]

Хе-хе!, вот видите как быстро дурь распространяется на умы!
Здесь какая-то тонкость - дело в том, что та же задача но с одной плоскостью и с равномерно распределённой нагрузкой на другой стороне нормально решается, а вот когда нагрузка прикладывается к контактной паре - через пилотный узел, что-то никак.

[aka raduntsev-nias]

Цветочек. Нужно разбираться, я корректировал только контакт в классике и все. Остальное не трогал

а ну дак конечно. у вас контакт твердое тело- плоскость

в одной из книг Огородниковой об этом написано

вот код

Код:

[Выделить все]

/prep7
ET,1,SOLID185   
!*  
KEYOPT,1,2,0
KEYOPT,1,6,1
KEYOPT,1,10,0   
!*  
!*  
MPTEMP,1,0                             
MPDATA,EX,1,,2.1E11            
MPDATA,PRXY,1,,0.3  
TB,BKIN,1,1,2,0
TBTEMP,0  
TBDATA,1,2.6E8,0,,,,  

real,1
CYL4,0,0,17.5, , , ,8   
vsweep,all


! Создание контактной пары на нижней площадке
ET,2,CONTA175 ! можно делать conta173
keyopt,2,12,3 ! контакт bonded
KEYOPT,2,5,4
ET,3,TARGE170 ! 3-D TARGET ELEMENT
KEYOPT,3,2,1
! Геометрия твердой поверхности = target
n,30001,30,-30
n,30002,30,30
n,30003,-30,30
n,30004,-30,-30
! Создание площадки
alls
real,2
type,3
tshap,quad
E,30001,30002,30003,30004 ! плоскость площадки элемент
! Сама поставит нормаль
tshape,pilo
EN,6001,30001 ! Создаётся элемент 6001 pilo для контроля
allsel,all
type,2
NSEL,S,LOC,z,0 ! Выделяю плоскость c z=0;
NSEL,U,NODE,,30001,30004
ESURF,, ! GENERATE
allsel,all
! Эта пара начально закрыта - всё хорошо

real,3
! Создание контактной пары на инденторе
ET,4,CONTA175 
KEYOPT,4,5,4
keyopt,4,12,3 ! контакт bonded
ET,5,TARGE170 
KEYOPT,5,2,1
! Геометрия твердой поверхности = target
n,30005,30,-30,8
n,30006,30,30,8
n,30007,-30,30,8
n,30008,-30,-30,8
! Создание площадки
alls
type,5
tshap,quad
E,30005,30006,30007,30008 ! плоскость площадки элемент
tshape,pilo
EN,7001,30005 ! Создаётся элемент 7001 pilo для контроля

allsel,all
type,4
NSEL,S,LOC,z,8 ! Выделяю плоскость c z=8;
NSEL,U,NODE,,30005,30008
esurf,,
allsel,all
fini
/solu
nlgeom,on
lnsrch,on
d,30001,all,0
d,30005,uz,-0.5
d,30005,ux,0
d,30005,uy,0
solve

фу. устал. WB рулит!

[Цветочек]

процитируйте пжлста эту фразу из Огородниковой, по поводу контакта с плоскостью.