Научный форум dxdy

Математика, Физика, Computer Science, LaTeX, Механика и Техника, Химия, Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки
Текущее время: Чт мар 11, 2010 07:57:50
Для набора любых формул следует использовать тег [math]. В противном случае сообщение будет отправлено в карантин.
Видите оффтопик? Жмите Пожаловаться на это сообщение
С Правилами Научного форума можно ознакомиться здесь.
Халявы здесь нет. На нашем форуме не решают задачи за вас.
Нужна подсветка синтаксиса? Есть такая возможность!
Попробуйте новый поиск по математическим формулам.


Часовой пояс: UTC + 3 часа




Начать новую тему Ответить на тему  [ Сообщений: 49 ]  На страницу Пред.  1, 2, 3, 4
Автор Сообщение
 Не в сети
 Re: Контактные задачи в Ansys
СообщениеПт июл 03, 2009 19:41:05 

Появился: 29/03/09
Сообщения: 35
Откуда: СРАНЬ
Уважаемые друзья!
Уже говорилось, что учет трения в Ansys представляет собой определенную сложность и что разные опции контакта существенно меняют поведение контактирующих поверхностей. Вопрос, который хочу задать связан с трением. Уже публиковала исходный код своих опытов, но он был неформатирован и не верен. Теперь вроде как все отладила и осталось только решить вопрос с поворотом и трением. Всегда появляется ошибка - negative jacobian 2d, Я понимаю что это значит: что-то неверно в постановке задачи. Привожу скрипт с описанием. Очень буду счастлива если поможите!

/prep7
*SET,R,1
*SET,ar,0.98
*SET,delt,0.004
smrt,off
ET,1,PLANE182, , ,2
R,1,,,,,,
KEYOPT,1,6,1 !mixed formulation - p-u
keyopt,1,3,2 ! плоская деформация
MP,EX,1,2.82 ! YOUNG'S MODULUS [MPA]
MP,NUXY,1,0.49967 ! POISSON'S RATIO
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,0,0
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
! Разбиение на элементы тела 1
LESIZE,3,,,7
LESIZE,4,,,7
AL,1,3,2,4
mat,1
type,1
real,1
ESIZE,,125
AMESH,all
! Сдвигаем начало координат на -R-ar-delt по y
local,11,1,0,-R-ar-delt,0
! Геометрия второго тела
k,20000,ar+delt,80
k,20001,ar,80
k,20002,ar+delt,100
k,20003,ar,100
L,20000,20002 ! 5 линия
L,20001,20003 ! 6 линия
L,20000,20001! 7 линия
L,20002,20003! 8 линия
AL,5,6,7,8
et,6,82
KEYOPT,6,6,0 ! pure displacement
keyopt,6,3,2 ! плоская деформация
R,6,,,,,,,,
MP,EX,6,1.54 ! YOUNG'S MODULUS [MPA]
MP,NUXY,6,0.2
TB,MELAS,6,1,4
TBTEMP,0
TBPT,defi,0.1,0.154 ! логарифмические деформации и истинные напряжения!
tbpt,defi,0.3,0.3*1.54
tbpt,defi,0.7,0.7*1.54
tbpt,defi,0.95,0.95*1.54
real,6
mat,6
type,6
ESIZE,,165
AMESH,2
! Контактная пара 1
ET,2,CONTA172 ! можно делать conta171
ET,3,TARGE169 ! 2-D TARGET ELEMENT
keyopt,2,12,5
R,2,ar , , -2000
! Геометрия твердой поверхности = target
n,30003,ar,80! узел для твердой поверхности
k,30003,ar,80
n,30002,ar,100! узел твердой поверхности
k,30002,ar,100
n,30001,0,0,0! пилотный узел
k,30001,0,0,0
! Создание цели
alls
real,2
type,3
tshap,arc
E,30002,30003,30001
tshape,pilo
E,30002,30003,30001
NSEL,S,NODE,,30001,30003
ESLN,s,0
ESURF,,REVERSE
! Создание активной контактной поверхности
alls
type,2
real,2
lsel,s,line,,6
nsll,s,1
esurf
! Контактная пара 2
alls
ET,4,CONTA172
ET,5,TARGE169
keyopt,4,5,4 ! редукция проникновения и зазора
R,3,ar+delt,,-200000
TB,FRIC,5,,,ISO
TBDATA,1,0.2
mat,5
! Создание цели
real,3
type,5
lsel,s,line,,1
nsll,s,1
esurf
! Создание активной контактной поверхности
alls
type,4
lsel,s,line,,5
nsll,s,1
esurf
! конец подготовки
alls
finish

! Блок решения
/solu
NROPT,UNSYM
alls
antype,static
nlgeom,on
deltime,0.01,0.01,0.01
lsel,s,line,,2
nsll,s,1
D,all,all,0
alls
D,30002,UY,0.0033
D,30002,ux,0
time,1
alls
lsel,s,line,,6
nsll,s,1
D,all,ux,0
alls
solve
time,2
lsel,s,line,,6
nsll,s,1
DDele,all,all
alls
D,30002,rotz,0.03
D,30002,UY,0
D,30002,ux,0
solve
/post1
/DIST,1, 0.1
set,1
PLNS,s,eqv
ANTIME,20,0.5, ,1,1,1,2

 Профиль  
                  
 Не в сети
 Re: Контактные задачи в Ansys
СообщениеСр июл 08, 2009 07:45:24 
Заслуженный участник
Аватара пользователя
Годы на форумеГоды на форуме
Появился: 11/04/07
Сообщения: 713
Откуда: Москва
Попробуйте позапускать следующий макрос. Нижний цилиндр поворачивается на 1 радиан, верхний чуть больше.
/prep7
ex,1,1
nuxy,1,0.3
et,1,82
r=0.1
dr=0.02
k,1
k,2,r
k,3,r+dr
k,4,0,r
k,5,0,r+dr
k,6,-r
k,7,-r-dr
k,8,0,-r
k,9,0,-r-dr
LARC,2,4,1,r, !1
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
lstr,4,5
lstr,6,7
lstr,8,9
al,1,10,5,9
LESIZE,9,,,3,,,,,0
LESIZE,10,,,3,,,,,0
LESIZE,11,,,3,,,,,0
LESIZE,12,,,3,,,,,0
al,2,11,6,10
al,3,12,7,11
al,4,9,8,12
type,1
real,1
mat,1
MSHAPE,0,2D
MSHKEY,1
elsize,0.01
AMAP,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
nuxy,2,0.3
et,2,82
r1=0.3
dr1=0.02
dy=-r-dr-r1-dr1
k,11,0,dy
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
lstr,14,15
lstr,16,17
lstr,18,19
LESIZE,21,,,3,,,,,0
LESIZE,22,,,3,,,,,0
LESIZE,23,,,3,,,,,0
LESIZE,24,,,3,,,,,0

num=12
al,num+1,num+10,num+5,num+9
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
real,2
mat,2
MSHAPE,0,2D
MSHKEY,1
elsize,0.01
AMAP,5,12,13,15,14
AMAP,6,15,14,17,16
AMAP,7,17,16,19,18
AMAP,8,19,18,13,12
mat,3
ET,3,TARGE169
KEYOPT,3,2,1
type,3
R,3,r,0,
real,3
tshap,circ
n,100000,
en,100000,100000
TSHAP,PILO
En,100001,100000
et,4,CONTA172
KEYOPT,4,12,5
type,4
real,3
lsel,s,,,1,4
nsll,s,1
esurf
allsel,all
mat,5
ET,5,TARGE169
KEYOPT,5,2,1
type,5
R,5,r1,0,
real,5
tshap,circ
n,100010,0,dy
en,100010,100010
TSHAP,PILO
En,100011,100010
et,6,CONTA172
KEYOPT,6,12,5
type,6
real,5
lsel,s,,,num+1,num+4
nsll,s,1
esurf
allsel,all

MU,7,0.8
mat,7
ET,7,TARGE169
type,7
real,7
r,7
RMODIF,7,13,0.0002
lsel,s,,,5,8
nsll,s,1
esurf
et,8,CONTA172
KEYOPT,8,4,2
KEYOPT,8,10,5
KEYOPT,8,2,3

type,8
real,7
lsel,s,,,num+5,num+8
nsll,s,1
esurf
allsel,all
/solu
CNVTOL,U,1,0.0001,2,
deltim,0.1,0.1,0.1
time,1
nlgeom,on
d,100000,ux
d,100000,uy,-0.1*dr
d,100000,rotz,0
d,100010,ux
d,100010,uy,0.0
d,100010,rotz,0
solve
deltim,0.005,0.005,0.005
time,2
d,100010,rotz,-1.0
DDELE,100000,ROTZ
outres,all,all
solve

 Профиль  
                  
 Не в сети
 Re: Контактные задачи в Ansys
СообщениеСр июл 08, 2009 23:27:09 

Появился: 29/03/09
Сообщения: 35
Откуда: СРАНЬ
Уважаемые друзья!
Спасибо за помощь, которую вы оказываете. Я поработала с предложенным скриптом, и обнаружила, конечно, что работает он очень быстро.
Когда я недавно писала о том, что в скрипте, опубликованном мной на форуме,возникает ошибка negative jacobian как раз в начале проворота, думала что это связано с моим стилем написания и образом мыслей. Я использовала следующие значения диаметров и толщин r=0.3, dr=0.02 верхний и r1=0.3, dr1=0.02 нижний цилиндры в предложенном скрипте. Также как и ранее возникает negative jacobian - опять же в начале проворота. С чем это связано пока не знаю. Появляется сообщение о непозволительно больших деформациях элемента и далее - ошибка negative jacobian. Буду признательна за советы и помощь. Спасибо большое!

 Профиль  
                  
 Не в сети
 Re: Контактные задачи в Ansys
СообщениеСр янв 20, 2010 11:59:44 

Появился: 20/01/10
Сообщения: 1
!

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 49 ]  На страницу Пред.  1, 2, 3, 4

Часовой пояс: UTC + 3 часа



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и гости: 0


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения

Найти:

Темы с похожим названием

 Темы   Автор   Ответы 
Электрuчество, задачи с подвохами!

в форуме Помогите решить / разобраться (Ф)

invisible1

13

Куда можно посылать олимпиадные задачи?

в форуме Олимпиадные задачи (М)

Legioner93

7

Задачи на кинематику, динамику.

в форуме Карантин

Makino

3

срочно!!!! помогите пожалуйста решить задачи!!!

в форуме Карантин

anna-kas

1

помогите решить задачи по финансам

в форуме Карантин

anna-kas

1

Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group