| Правила | Регистрация | Пользователи | Поиск | Сообщения за день | Все разделы прочитаны |  Справка по форуму | Файлообменник |

Вернуться   Форум DWG.RU > Программное обеспечение > Расчетные программы > ANSYS > Пара вопросов по контактным задачам Ansys

Пара вопросов по контактным задачам Ansys

Ответ
Поиск в этой теме
Непрочитано 03.09.2010, 22:28
Пара вопросов по контактным задачам Ansys
Цветочек
 
Регистрация: 29.03.2009
Сообщений: 38

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. Какие могут быть дефекты контакта если контакт изначально открыт - хотя детали рядом, и сколько не переставляю нормали - всё так и есть?
Просмотров: 13222
 
Автор темы   Непрочитано 06.09.2010, 16:24
#41
Цветочек


 
Регистрация: 29.03.2009
Сообщений: 38


Простите что задерживаю - были важные дела.
Оставляю все комментарии - ибо раз хотите изучить apdl - примеры самое то, пусть и неверные

Код:
[Выделить все]
/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,3 ! материал 3 не создан но мы все равно его используем что отвлеч ансис
ET,3,TARGE169 ! элемент типа недеформируемая цель
KEYOPT,3,2,1 ! граничные условия для этого элемента будут задаваться пользователем
type,3
R,3,r,0, ! радиус цели r
real,3 
tshap,circ  ! определяем форму поверхности цели
n,100000, ! узел с номером 100000 имеет координаты 0,0,0
en,100000,100000 ! задаем элемент circ номера 100000 вида цель с узлом 100000 центре верхнего
TSHAP,PILO ! создаем мастер узел
En,100001,100000 ! элемент 100001 pilot создаем с узлом в центре верхнего 
et,4,CONTA172 ! декларируем элемент вида контакт 172
KEYOPT,4,12,5 ! он всегда будет always bonded
type,4
real,3! используем  реальные константы номер 3
lsel,s,,,1,4 ! выбираем линии с 1 по 4 
nsll,s,1! выделяем там все узлы
esurf  ! покрываем все выделенные  узлы контактными элементами
allsel,all ! выделяем все!
! создание второй контактной пары -нижний цилиндр - между жестким валом и покрытием
mat,5 ! материал номер 5
ET,5,TARGE169 ! элемент 5 - недеформируемая цель
KEYOPT,5,2,1 ! граничные условия задает пользователь
type,5
R,5,r1,0, ! радиус поверхности цели - мишени
real,5
tshap,circ ! задаем форму поверхности цели
n,100010,0,dy ! создать узел номером 100010 с координатами 0,dy,0
en,100010,100010! создать элемент типа цель номером 100010
TSHAP,PILO  ! создаем мастер-узел
En,100011,100010 ! создание элемента 100011 с мастер-узлом
et,6,CONTA172 ! элемент вида conta172
KEYOPT,6,12,5  ! этот контактный элемент будет always bonded
type,6 ! создаем теперь контактные элементы
real,5 ! реальные константы номер 5 
lsel,s,,,num+1,num+4! выделяем линии с 13 по 17
nsll,s,1 ! выделяем все узлы на выделенных линиях
esurf ! покрываем эти узлы контактными элементами
allsel,all ! выделяем все!
! контактные пары отвечающие за внутренние контакты - между слоями и жесткими валами созданы
! Создание контактной пары для контакта двух слоев!
MU,7,0.95 ! создаем седьмой материал к коэффициентом трения 0.8
mat,7 ! задаем этот самый 7 материал
ET,7,TARGE169 ! используем элемент типа деформируемая цель - мишень
type,7
real,7! реальные постоянные номер 7
r,7 ! реальные постоянные 7
RMODIF,7,13,0.0002 !меняем множество реальных постоянных 7, а именно 13 постоянную - cohe = 0.0002 
lsel,s,,,5,8! выделяем  кривые с 5-8 - внешний круг верхнего цилиндра
nsll,s,1! выделяем все узлы
esurf  ! покрываем элементами target169
et,8,CONTA172 ! 8 вид элемента - conta172
KEYOPT,8,4,2 ! Location of contact detection point : On nodal point - normal to target surface
KEYOPT,8,10,5 ! Contact stiffness update : Each iteration based on current mean stress of underlying elements ! (individual element based).
KEYOPT,8,2,3 ! Contact algorithm : Lagrange multiplier on contact normal and penalty on tangent
type,8
real,7
lsel,s,,,num+5,num+8 ! выделяем линии с 17 по 20 -  внешние линии нижнего цилиндра
nsll,s,1 ! выделяем все узлы на выделенных кривых
esurf ! покрываем их все контактными элементами
allsel,all ! выделяем все!
! блок решения
/solu
CNVTOL,U,1,0.0001,2,  ! определить сходимость по перемещениям, точность 0.01%, норма сходимости L2 
deltim,0.1,0.1,0.1 ! шаг решения 0.1 - очень быстрое решение
time,1 ! время идет до 1 
nlgeom,on ! включить нелинейные геометрические величины - большие перемещения и вращения
d,100000,ux ! запретить узлу 100000 перемещения по x
d,100000,uy,-0.1*dr ! задать узлу 100000 перемещения по y вниз на !0.1*dr - верхний цилиндр
d,100000,rotz,0 ! вращение узла 100000 невозможно
d,100010,ux ! закрепить центральный узел нижнего цилиндра по x
d,100010,uy,0.0 ! закрепить центральный узел нижнего цилиндра по y
d,100010,rotz,0 ! нижний цилиндр не вращается
solve ! решаем статическую задачу - вращения нет
Цветочек вне форума  
 
Автор темы   Непрочитано 07.09.2010, 18:04 Уф
#42
Цветочек


 
Регистрация: 29.03.2009
Сообщений: 38


Бра40 - Вы рано расстраиваетесь - Слава Богу! Здесь нету никаких сил кориолиса и так далее - у меня эти самые цилиндры не вращаются, это просто задача Герца в чистом виде! Просто она большая и всё
Цветочек вне форума  
 
Непрочитано 07.09.2010, 18:11
#43
derik

Конструякер
 
Регистрация: 27.07.2009
Подмосковье
Сообщений: 581
Отправить сообщение для derik с помощью Skype™


да не расстраиваюсь я если задача Герца тада я в теме, а если добавлять вращение...тут нада хелпа почитать :-)
Вы же написали в макросе: что решение есть, а вот вращения нет!
derik вне форума  
 
Автор темы   Непрочитано 07.09.2010, 18:13
#44
Цветочек


 
Регистрация: 29.03.2009
Сообщений: 38


Я написала для себя что запрещаю вращение по rotz - как факт чтобы не было наклонов. ТАМА НЕТУ ВРАЩЕНИЯ! Я ж писала чтобы себя успокоить!
Цветочек вне форума  
 
Непрочитано 07.09.2010, 18:23
#45
derik

Конструякер
 
Регистрация: 27.07.2009
Подмосковье
Сообщений: 581
Отправить сообщение для derik с помощью Skype™


писали для себя а не понял я хD
derik вне форума  
 
Автор темы   Непрочитано 07.09.2010, 18:30 Как многого я хотела!
#46
Цветочек


 
Регистрация: 29.03.2009
Сообщений: 38


Тут вот в чём дело! Я изначально задумывала и вращение делать, но только когда смогу решить задачу Герца в статитке, но поскольку у меня даже статика обламилась, то о какой динамике думать?
Цветочек вне форума  
 
Непрочитано 07.09.2010, 18:58
#47
derik

Конструякер
 
Регистрация: 27.07.2009
Подмосковье
Сообщений: 581
Отправить сообщение для derik с помощью Skype™


так пойдет?
(эт набросал на скрорую руку)
Код:
[Выделить все]
fini
/clear,nostart
/outp,outfile,txt,C:\ansys\output
/cwd,'C:\ansys\out'

t=.05 ! поставте значение которое нужно
t1=.05 ! поставте значение которое нужно
r=.3 ! поставте значение которое нужно
r1=.3 ! поставте значение которое нужно
a=r+r1 ! поставте значение которое нужно
p=5e2 ! поставте значение которое нужно

/prep7

et,1,183

mptemp,,,,,,,
mptemp,1,0
mpdata,ex,1,,2.06e8
mpdata,prxy,1,,.3

pcirc,r,r-t,270,450,
pcirc,r1,r1-t1,270,450,
agen,,2,,,,-a,,,,1

mshkey,1
mshape,2,2d
esize,t/4
amesh,all

kref,6,,,1,2,,
kref,1,,,1,2,,

et,2,169
et,3,172
keyopt,3,2,1
keyopt,3,5,1
keyopt,3,9,3
keyopt,3,10,5
keyopt,3,12,5
mp,mu,3,.3

mat,3
real,3

type,2
lsel,s,,,5
nsll,s,1
esln
esurf
type,3
lsel,s,,,1
nsll,s,1
esln
esurf
!cncheck,auto

lsel,s,,,3,7,4
dl,all,,all
lsel,s,,,2,8,2
dl,all,,symm
allsel,all
fini

/sol

/gst,1
antype,0
nlgeom,1
nsubst,50,1000,20
outres,erase
outres,all,all
time,1

solve
fini
/post1
set,first
/expand,2,rect,half,.00001

Последний раз редактировалось derik, 07.09.2010 в 21:21.
derik вне форума  
 
Автор темы   Непрочитано 08.09.2010, 01:35
#48
Цветочек


 
Регистрация: 29.03.2009
Сообщений: 38


Спасибо за поддержку, Бра4о! Я посмотрю ваш код
Цветочек вне форума  
 
Непрочитано 08.09.2010, 14:16
#49
derik

Конструякер
 
Регистрация: 27.07.2009
Подмосковье
Сообщений: 581
Отправить сообщение для derik с помощью Skype™


Код:
[Выделить все]
fini
/clear,nostart
/outp,outfile,txt,C:\ansys\output
/cwd,'C:\ansys\out'

r=.03 ! поставте значение которое нужно
r1=.09 ! поставте значение которое нужно

/prep7

et,1,183

mptemp,,,,,,,
mptemp,1,0
mpdata,ex,1,,7e8
mpdata,prxy,1,,.3

cyl4,,r,r,270,,450,
cyl4,,-r1,r1,270,,450

mshkey,0
mshape,2,2d
esize,r/7
amesh,all

kref,5,,,2,3,,
kref,1,,,2,3,,

et,2,169
et,3,172
keyopt,3,2,1
keyopt,3,5,1
keyopt,3,10,5
mp,mu,3,.3
r,3,0,0,11,0,0,40

mat,3
real,3
type,2
lsel,s,,,4
nsll,s,1
esln
esurf
type,3
lsel,s,,,1
nsll,s,1
esln
esurf
!cncheck,auto

ksel,s,,,3,6,3
dk,all,,,,,,all,
lsel,s,,,2,3
lsel,a,,,5,6
dl,all,,symm
allsel,all
fini

/sol

/gst,1
antype,0
nlgeom,1
nldiag,nrre,on
nropt,unsym,,on
nsubst,50,5000,10
deltim,0.1,0.1,0.1
outres,erase
outres,all,all
time,1

solve
fini
/post1
set,first
/expand,2,rect,half,.00001
а вот если два цилиндра
derik вне форума  
 
Автор темы   Непрочитано 08.09.2010, 22:31 Благодарствую за помощь!
#50
Цветочек


 
Регистрация: 29.03.2009
Сообщений: 38


Спасибо Вам бра4о большое! К сожалению то, что вы написали лишний раз подтверждает что свои ошибки каждый должен исправлять сам. К сожалению, я еще свои не исправила.
Всё же придется мне искать ошибки в переходе от 2-d примера в 3-d ситуацию - где колёса почему-то кривятся.
Цветочек вне форума  
 
Непрочитано 08.09.2010, 23:00
#51
derik

Конструякер
 
Регистрация: 27.07.2009
Подмосковье
Сообщений: 581
Отправить сообщение для derik с помощью Skype™


:-) контакт решен же...а по поводу кручения - команда omega и cnomega...большим помочь не могу :-(
derik вне форума  
 
Непрочитано 10.09.2010, 15:03
#52
derik

Конструякер
 
Регистрация: 27.07.2009
Подмосковье
Сообщений: 581
Отправить сообщение для derik с помощью Skype™


Цветочек, как там дела у Вас? что нибудь помогло?
derik вне форума  
 
Автор темы   Непрочитано 10.09.2010, 23:36 Приветик бра4о
#53
Цветочек


 
Регистрация: 29.03.2009
Сообщений: 38


У меня к сожалению начался рабочий год с 1 сентября, поэтому не могу полностью отдать себя решению этой проблемы. Нет, практически ничего не помогло
Цветочек вне форума  
 
Непрочитано 19.09.2010, 21:39
#54
derik

Конструякер
 
Регистрация: 27.07.2009
Подмосковье
Сообщений: 581
Отправить сообщение для derik с помощью Skype™


Цветочек, а Вы не пробывали добавить плотность материалу и приложить команду омега?
derik вне форума  
 
Автор темы   Непрочитано 19.09.2010, 22:13 А зачем?
#55
Цветочек


 
Регистрация: 29.03.2009
Сообщений: 38


Нет ли тут мании по поводу вращения? Уже устала уверять что вращения нету - просто сложный контакт - конечно можно добавить нулевое вращение и плотность, и ?
Вы наверное думаете что раз я так хочу пилотные узлы значит хочу вращать? Нет, вращать не хочу. Хочу задавать контакт через перемещения колёс
Цветочек вне форума  
 
Непрочитано 19.09.2010, 23:22
#56
derik

Конструякер
 
Регистрация: 27.07.2009
Подмосковье
Сообщений: 581
Отправить сообщение для derik с помощью Skype™


чтобы получить результат нужно смоделировать вращение валов и получим напряжения в точке соприкосновения.
в покое правильного результата не будет...цилиндры, даже в двумерке должны крутиться.
вращение - запрещаете перемещения по y, раздвигаете колеса и по x задаете перемещения друг к другу. либо одно колесо жестко закрепляете и приближаете второе, опять по x:
asel,s,,,nA
da,nA,ux,1(-1)
derik вне форума  
 
Автор темы   Непрочитано 19.09.2010, 23:41 Уф
#57
Цветочек


 
Регистрация: 29.03.2009
Сообщений: 38


У меня и в статике не работает эта задача - в 2-d всё крутиться , но только всё плохо крутиться нельзя в ансис смоделировать качение - то есть процессы трения при качении и соответственно нельзя получить верное НДС в месте контакта
Поэтому забудьте про качение - я его сделала в 2-d и всё неверное.
Давайте просто оставим качение, меня сейчас это уже не трогает - я сделала это другим методом. Но вот расчётная схема про пилоты и два колеса всё еще нужна - правда без вращения. Поверьте мне Бра4o. Если вы озаботитесь кодом на 3 странице наших разговоров и поможете мне понять почему пилоты неверно работают я правда буду счастлива .
А если вам просто интересно как можно задавать вращение, могу помочь - d,pilot node, rotz,1 радиан

Последний раз редактировалось Цветочек, 19.09.2010 в 23:47.
Цветочек вне форума  
 
Автор темы   Непрочитано 03.10.2010, 15:38 Определение длины кривой линии контакта
#58
Цветочек


 
Регистрация: 29.03.2009
Сообщений: 38


Кое что написала - может и неплохо - в любом случае поправите если надо
Эти строки определяют ширину площадки контакта - работает с элементами 171, 172. С 175 не работает.
Надеюсь всем пригодиться! как и мне
Код:
[Выделить все]
/post1
allsel,all
lengs=0
ereal=2
*get,nodmax,elem,,num,max
*do,i,1,nodmax,1
*GET,elemstype,ELEM,i,ATTR,ENAM
  *if,elemstype,eq,172,then
    elemstype=0
    *GET,realelem,ELEM,i,ATTR,REAL
      *if,realelem,eq,ereal,then
           realelem=0     
        *GET,flag1,NODE,nelem(i,1),CONT,STAT
        *GET,flag2,NODE,nelem(i,2),CONT,STAT
           *if,flag1+flag2,ge,4,then 
             flag1=0
             flag2=0      
             delta1=0
             *GET,delta1,ELEM,i,LENG 
              *status,i
              lengs=lengs+delta1
          *endif
       *endif
   *endif
*enddo
 *status,lengs

Последний раз редактировалось Цветочек, 11.10.2010 в 14:26.
Цветочек вне форума  
Ответ
Вернуться   Форум DWG.RU > Программное обеспечение > Расчетные программы > ANSYS > Пара вопросов по контактным задачам Ansys

Размещение рекламы
Опции темы Поиск в этой теме
Поиск в этой теме:

Расширенный поиск


Похожие темы
Тема Автор Раздел Ответов Последнее сообщение
несколько вопросов по пострпоцессору ANSYS Mechanical autocader ANSYS 17 22.06.2010 21:37
К проектировщикам по водоотводу, водопроводу, канализации есть пара вопросов. FRAER Водоснабжение и водоотведение 20 04.02.2010 15:15