[Цветочек]

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. Какие могут быть дефекты контакта если контакт изначально открыт - хотя детали рядом, и сколько не переставляю нормали - всё так и есть?

[СергейД]

друзья! Уж извините. занят очень срочной и важной работой.аэродинамика вовсе. Форум просматривать нет времени. ЕСли действительно важные вопросы (не учебные..) пишите в личку.
Только советы общего плана дать могу.
=
Пилотными узлами не пользуюсь принципиально (когда-то что-то не понравилось в них = возможно, что тоже самое, о чем волнуется Цветочек) и разбираться в них сейчас нет времени.
с обычным контактом какие проблемы?

А всех интересующихся (кстати) приглашаю на курсы ансис (начало октября) и на конференции (конец октября и середина ноября)в МГСУ (MGSU.ru).

[aka raduntsev-nias]

Цветочек, извините, не могу найти) Пропадаю на время. Tekla(AdvanceSteel) ждет, кроме расчетов я конструктор
мало времени на ansys остается

[derik]

Сообщение от aka raduntsev-nias Цветочек, извините, не могу найти)

а что за книга? как называется, а то так каждую открывать и смотреть автора...

[aka raduntsev-nias]

Offtop: хоть убейте, не помню.
контакт - дело темное
WB = future
APDL = always

[Цветочек]

Вы уже мне сильно помогли друзья! Постоянно забываю про conta175. Моей основной заботой является контакт.
В столе лежит стопка неоконченных программ поскольку вот такие - как выше,дефекты, портят мне жизнь. Одну из головных болей я сняла. Но есть и другая!,Как заменить пилотные узлы если сама задача кажется что сформулирована через них, например, внутри обода колёс имеется контактная пара так что на этом рисунке есть по сути 3 контактные пары - в идеале нужна еще одна пара поскольку у одного из цилиндров бандаж 2 слойный. Если например принять контакт с пилотами, то получается тупо - один узел закрываем от перемещений, другой подвигаем вверх. А что тут поделать если без пилотов?
В моём примере контакт правда кривой(несколько не физичен), но думаю хорошая иллюстрация.
Толщина пластин - чисто номинально плоскости, это безразлично.

[derik]

Offtop: а как книга то называется?

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

[aka raduntsev-nias]

Сообщение от Цветочек 3 контактные пары - в идеале нужна еще одна пара поскольку у одного из цилиндров бандаж 2 слойный.

возможно контакт симметричный, если нет, то слой и target, и contact
те 2 контакта

[derik]

Код:

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

fini
/clear,nostart  
/outp,outfile,txt,C:\ansys\output   
/cwd,'C:\ansys\out' 

a=.600
h=.08
t=.02
r=.35
p=3e6

/prep7

et,1,185

mptemp,,,,,,,
mptemp,1,0  
mpdata,ex,1,,3e11
mpdata,prxy,1,,.3
mptemp,,,,,,,
mptemp,1,0
mpdata,ex,2,,6e11
mpdata,prxy,2,,.3

cylind,r,,,h,,90,
blc5,a/2,a/2,a,a,-t 
vgen,2,2,,,,,h+t,,,0
vsel,s,,,2,3
vatt,2, 
allsel,all  
!sfa,13,1,pres,p

mshkey,1
mshape,0,3d 
esize,t*2
vsweep,all  
allsel,all  

et,2,170
et,3,174
keyopt,3,2,1
keyopt,3,5,1
keyopt,3,9,1
keyopt,3,10,5
keyopt,3,12,5

mp,mu,3,.3  
mat,3   

r,3,,,.01,,,50
r,4,,,.01,,,50

type,2  
real,3  
asel,s,,,12 
nsla,s,1
esln,s,0
esurf   
type,3  
asel,s,,,2  
nsla,s,1
esln,s,0
esurf   

type,2  
real,4  
asel,s,,,7  
nsla,s,1
esln,s,0
esurf   
type,3  
asel,s,,,1  
nsla,s,1
esln,s,0
esurf   
allsel,all  
!cncheck,auto

asel,s,loc,y
asel,a,loc,x
da,all,symm 
da,6,all
da,13,uz,-1.3
allsel,all  
fini

/sol

/gst,1  
nldiag,nrre,on  
nropt,unsym,,on 
antype,0
nlgeom,1
nsubst,100,1000,50   
outres,erase
outres,all,all  
time,1
solve   
fini

/post1  

set,first   
fini

вот что получилось у меня, правда тоже не фонтан...нужно еще поиграться с опциями - но это позже чуть-чуть, счас другие дела есть немножко

[Цветочек]

В том смысле что их тяжело перевести в безпилотное состояние. Поэтому если кто-то не знает пилотов, я ему объясню - это не сложно.
За коды спасибо!

Итак у меня трудный орешек. Рисунок имеется выше. Задача такова:
два трёхмерных цилиндра из стали контактируют друг с другом, но при этом их внутренние диаметры остаются без изменений - ну как будто начинка этих цилиндров абсолютно твёрдая. 3 контактные пары:
две внутри цилиндров, и одна между ними - обычная задача Герца.
Код,

Код:

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

/prep7
ex,1,1 ! модуль юнга 1 МПа
nuxy,1,0.3 ! коэффициент пуассона 0.3
et,1,82 ! элемент код 82
r=0.1 ! внутренний радиус верхнего цилиндра
dr=0.02 ! толщина бандажа верхнего цилиндра
k,1 ! 1 точка = координаты 0,0,0
k,2,r ! 2 точка = r,0,0
k,3,r+dr ! 3 точка r+dr,0,0
k,4,0,r ! 4 точка 0,r,0
k,5,0,r+dr ! 5 точка 0,r+dr,0
k,6,-r !  6 точка -r,0,0
k,7,-r-dr ! 7 точка -r-dr,0,0 ! здесь -r-dr - наружный диаметр 
k,8,0,-r ! 8 точка 0,-r,0
k,9,0,-r-dr! 9 точка 0,-r-dr
LARC,2,4,1,r, !1 кривая по 3 точкам и радиусу
LARC,4,6,1,r, 
LARC,6,8,1,r, 
LARC,8,2,1,r, 
LARC,3,5,1,r+dr, !5
LARC,5,7,1,r+dr, 
LARC,7,9,1,r+dr, 
LARC,9,3,1,r+dr, 
lstr,2,3 !9 ! прямая линия в глобальной системе координат - линия 9
lstr,4,5
lstr,6,7
lstr,8,9
al,1,10,5,9 ! площадь - четверть окружности малого вала 1 площадь
LESIZE,9,,,3,,,,,0 ! разбиваем линии с 9 по 12 на три участка - число ! элементов по толщине
LESIZE,10,,,3,,,,,0 
LESIZE,11,,,3,,,,,0 
LESIZE,12,,,3,,,,,0 
al,2,11,6,10 ! 2 площадь
al,3,12,7,11 ! 3 - - 
al,4,9,8,12 ! 4 - -
type,1 ! применяем элемент 82
real,1 ! реальные константы номер 1
mat,1 ! материал 1
MSHAPE,0,2D ! сетка из 4-угольников!
MSHKEY,1 ! mapped meshing
esize,0.01 ! характерный размер 1 элемента - 0.01
AMAP,1,2,3,5,4 ! mapped mesh площади 1 по угловым точкам 2,3,5,4
AMAP,2,5,4,7,6 !
AMAP,3,7,6,9,8 !
AMAP,4,9,8,3,2 !


ex,2,1! материал номер 2
nuxy,2,0.3 ! коэффициент пуассона номер 2
et,2,82 ! элемент 82 номер 2
r1=0.3 ! внутренний диаметр нижнего цилиндра
dr1=0.02 ! толщина стенки нижнего цилиндра
dy=-r-dr-r1-dr1 ! величина dy - межосевое расстояние между !верхним и нижним цилиндрами
k,11,0,dy ! центр нижнего цил коорд 0,dy,0
k,12,r1,dy ! далее и ниже прорисовка нижнего цилиндра
k,13,r1+dr1,dy
k,14,0,r1+dy
k,15,0,r1+dr1+dy
k,16,-r1,dy
k,17,-r1-dr1,dy
k,18,0,-r1+dy
k,19,0,-r1-dr1+dy
LARC,12,14,11,r1, !13 дуги окружности нижнего цилиндра
LARC,14,16,11,r1, 
LARC,16,18,11,r1, 
LARC,18,12,11,r1, 
LARC,13,15,11,r1+dr1, !17
LARC,15,17,11,r1+dr1, 
LARC,17,19,11,r1+dr1, 
LARC,19,13,11,r1+dr1, 
lstr,12,13 !21 прямые линии с 21 по 24 для указания количества слоев элементов по толщине
lstr,14,15
lstr,16,17
lstr,18,19
LESIZE,21,,,3,,,,,0 ! разбивка вышеперечисленных линий на 3 подлинии
LESIZE,22,,,3,,,,,0 
LESIZE,23,,,3,,,,,0 
LESIZE,24,,,3,,,,,0 
num=12 ! количество линий в верхнем цилиндре
al,num+1,num+10,num+5,num+9 ! создание 4 подобластей для слоя нижнего цилиндра 
al,num+2,num+11,num+6,num+10
al,num+3,num+12,num+7,num+11
al,num+4,num+9,num+8,num+12
type,2 ! тип элементов теперь 2 - реальные константы 2 материал 2 
real,2
mat,2
MSHAPE,0,2D  ! сетка 4-угольников
MSHKEY,1 ! mapped meshing
esize,0.01 ! характерный размер элемента 0.01
AMAP,5,12,13,15,14 ! mapped meshing всех 4-подобластей нижнего цилиндра
AMAP,6,15,14,17,16 
AMAP,7,17,16,19,18 
AMAP,8,19,18,13,12 

mat,1
ET,10,SOLID185
TYPE,10
ESIZE,,1
VEXT,1,4,,,,0.05

mat,2
ET,12,SOLID185   
TYPE,12
ESIZE,,1
VEXT,5,8,,,,0.05

!!!!!!!!!!!!! внеконтактная геометрия, сетка и элементы с материалами !созданы
!!!!!!! создание первой контактной пары - верхний цилиндр - между жестким валом и покрытием
mat,3 ! материал 3 не создан но мы все равно его используем что отвлечь ансис
ET,13,TARGE170 ! элемент типа недеформируемая цель
KEYOPT,13,2,1 ! граничные условия для этого элемента будут задаваться пользователем
type,13
R,13,r,0, ! радиус цели r
real,13 
tshap,cyli  ! определяем форму поверхности цели
n,100000 ! узел с номером 100000 имеет координаты 0,0,0
n,120000,0,0,0.05
en,100000,100000,120000 ! задаем элемент cyli номера 100000 вида цель с узлом 100000 центре верхнего - здесь элемент 100000 это цилиндрический элемент цель.
TSHAP,PILO ! создаем мастер узел - элемент 100001
En,100001,100000 ! элемент 100001 pilot создаем с узлом в центре верхнего 
et,14,conta173 ! декларируем элемент вида контакт 174
KEYOPT,14,12,5 ! он всегда будет always bonded
type,14
real,13! используем  реальные константы номер 13  
CSYS,1
ASEL,S,LOC,x,r
NSLA,,1                    
ESURF              
csys,0
allsel,all

mat,4 ! материал номер 15 - его тоже нету
ET,15,TARGE170 ! элемент 15 - недеформируемая цель
KEYOPT,15,2,1 ! граничные условия задает пользователь
type,15
R,15,r1,0, ! радиус поверхности цели - мишени
real,15
tshap,cyli ! задаем форму поверхности цели
n,1000000,0,dy ! создать узел номером 10000 с координатами 0,dy,0
n,1200000,0,dy,0.05
en,1000000,1000000,1200000! создать элемент типа элемент-цель номером 1000000
TSHAP,PILO  ! создаем мастер-узел с номером 1000001
En,1000001,1000000 ! создание элемента 10001 с мастер-узлом
et,16,conta173 ! элемент вида conta173
KEYOPT,16,12,5  ! этот контактный элемент будет always bonded
type,16 ! создаем теперь контактные элементы
real,15 ! реальные константы номер 5 
local,11,1,0,dy,0
csys,11
asel,s,loc,x,r1
NSLA,,1                    
ESURF                
csys,0
allsel,all ! выделяем все!

! контактные пары отвечающие за внутренние контакты - между слоями и жесткими валами созданы

! Создание контактной пары для контакта двух слоев!
mat,5 ! задаем этот самый 17 материал Его тоже нет
ET,17,TARGE170 ! используем элемент типа деформируемая цель - мишень
type,17
real,17! реальные постоянные номер 17
r,17 ! реальные постоянные 17
csys,1 
esize,,1
asel,s,loc,X,R+DR
AMESH,ALL
allsel,all ! выделяем все!
et,18,conta173 ! 18 вид элемента - conta173
KEYOPT,18,4,2 ! L
KEYOPT,18,10,5 ! 
KEYOPT,18,2,3 ! 
type,18
real,17
local,11,1,0,dy,0
csys,11
asel,s,loc,x,r1+dr1
amesh,all
ESEL,S,REAL,,17 
ESURF,,REVERSE 
allsel,all
csys,0


/solu
nlgeom,on ! включить нелинейные геометрические величины - большие перемещения и вращения
nropt,unsym,,on ! Несимметричный решатель
lnsrch,on ! Включить line-search
pred,on ! Предсказывать
!CNVTOL,u,1,0.0001,2, ! определить сходимость по перемещениям, точность 0.01%, норма сходимости L2 
deltim,0.1,0.1,0.1
!deltim,0.001,0.0001,0.01 ! шаг решения 0.1 - очень быстрое решение
time,1 !
!
!d,100000,all
!ddele,10000,uy
!d,1000000,uy,0.1*dr
!d,1200000,uy,0.1*dr
!d,1000000,ux
!d,1200000,ux
 ! задать узлу 100000 перемещения по y вниз на 0.1*dr - верхний цилиндр
!d,1000000,ux ! закрепить центральный узел нижнего цилиндра по x
!d,1000000,uy,0.0 ! закрепить центральный узел нижнего цилиндра по y
!d,120000,all
!d,100000,all
!

d,all,rotx,0
d,all,roty,0
d,all,rotz,0
!solve ! решаем статическую задачу - вращения нет
d,120000,all
d,100000,all
d,1000000,uy,0.1*dr
d,1000000,ux,0
d,1000000,uz,0
d,1200000,uy,0.1*dr
d,1200000,ux,0
d,1200000,uz,0

Честно говоря вид этой проги вызывает ужос,Но - два цилиндра - вначале два блока, а потом последовательно 3 контактные пары, потом блок решения и всё
Думаю что эта задача достаточно сложна, как говориться дорогу осилит идущий, не так ли?
Всем спасибо!

[Кулик Алексей aka kpblc]

Цветочек, про тэг code не забывай

[Цветочек]

Простите меня! Я просто не знала.

[derik]

а в чем вопрос? задача решается, все хорошо!

[aka raduntsev-nias]

25 пост, дама хочет уйти от пилотных узлов
учебник в тему во вложении Александрова и Чебакова
учебник только что нашел в сети) может кому-то еще будет интересен
У Басова похожую задачу тоже наверное видели(SOLID95). Решается без пилотов.
вот еще

[derik]

ухты...не земетил
нужно попробывать, решить в 2Д, а там до трехмерки 1 шаг
будет время макрос напишу

[aka raduntsev-nias]

Offtop: ты уже вижу их штампуешь. переведи всю книжку в макросы)

[derik]

Offtop:
Девушка, а у Вас есть гайд по контакту от кадфем?

[aka raduntsev-nias]

Брачо, попробуй дать похожие закрепления как в примере. У меня нессиметричный контакт решается комп висит
Точечные закрепления убрать из кода и все
все идеально просто

в эти минуты жалеешь. что плохо понимаешь в apdl
суть не понятна
что на что действует чтобы решить в wb
скорей всего и не только в apdl
эх. ладно. поработать и спать)

[Цветочек]

1. Скажу честно Guide от CadFem это тупой перевод помощи Ансис
2. Все книги Александрова я хорошо знаю
3. Не точечные - а узловые.
4. У Басова видели, у него всё очень просто и потому грустно - у меня-то не так!

[derik]

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

[aka raduntsev-nias]

код в студию. будем изучать apdl
workbench освобождает на это время