Цистерна с гидростатическим давлением

Deshewle · 15.02.2013, 23:43

Всем привет.
Считаю напряжения, возникающие в опорах цистерны. Напряжения возникают из-за гидростатического давления. (в цистерне вода). Посчитал с помощью градиаента давления. Задача решена линейно (контакт не захотел использовать, т.к. он мне пока не поддается), совместил узлы на опорах (solid45) и узлы на оболочке (shell63) и получил решение. (фото ниже)

[img]http://img694.**********.us/img694/5319/tank1w.jpg[/img]
[img]http://img96.**********.us/img96/1006/tank2mf.jpg[/img]

Мне нужно проделать тоже самое, но без gradP, а, применив fluid80.

[img]http://img812.**********.us/img812/9001/tank4z.jpg[/img]
[img]http://img41.**********.us/img41/5167/tank3l.jpg[/img]

Так же делаю одинаковую сетку на shell и на опорах (solid45), чтобы узлы были совмещены. Через Inertia--Gravity прикладываю ускорение свободно падения. В итоге ничего не получается, делаю solve, после чего отображаются только 4 опоры, а shell+fluid пропадают. В чем может быть проблема?

Solution выдает следующее:

[img]http://img163.**********.us/img163/6662/tank5o.jpg[/img]

Что я делаю не так? Проблема с привязкой сеток?

derik · 16.02.2013, 11:36

перед тем как входить в раздел /sol выберите все: alls. а еще лучше вывесите код, отсюда не видать где ошибка закралась

Deshewle · 16.02.2013, 11:45

Alls выбрано. Все делал через GUI, в связи с чем в макросе куча лишнего. Может быть решение не идет из-за того, что узлы сетки fluid80 и узлы сетки shell63 не стыкуются должным образом? (не узел к узлу)

derik · 16.02.2013, 11:58

а вы не пробывали екструдов выдавить цистерну? тогда и сетка бы была из кирпичей а не треугольников.

Deshewle · 16.02.2013, 12:14

Если вы имеете в виду, сделать цистерну из того же solid45, задать gravity и получить решение, то это выход, но мне нужно использовать именно shell+fluid+solid45(опоры). Может быть есть идеи как можно "провернуть" данную задачу? (главное без контакта)

Deshewle · 21.02.2013, 00:35

Ребят, с проблемой разобрался. Как задать собственный вес для fluid80 в цистерне? (плотность, вязкость, модуль юнга - в material properties задал). Далее (через select) выделяю объем fluid80 и в через acel присваиваю по оси Y 9.8. Нуль реакции.

derik · 21.02.2013, 00:41

без кода не видно что Вы там прикладываете

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

Вам нужн овыделить все - alls а потом уже назначать гравитицию. а по поводу как конкретно к элементу - смотрите в хелпе.

Deshewle · 21.02.2013, 01:04

У меня узлы shell63 связаны с узлами fluid80, соответственно (пренебрегаем весом самой цистерны) ускорение свободного падения требуется только для жидкости (fluid). Если бы я нашел ответ в хелпе, то здесь бы ни в коем случае не появился. Кто-нибудь может подсказать как задать гравитацию для отдельного элемента?

derik · 21.02.2013, 01:11

выбираете элемент и присваиваете ему гравитацию, как Вы и делаете. а может единицы измерений пляшут?

Deshewle · 21.02.2013, 01:17

Хм, странно. Система SI. Когда на свободную поверхность оказываю давление, то все отлично решается, как только использую acel, напряжений нет вообще.
Ладно, буду разбираться, где собака зарыта, раз вы говорите, что я правильно делаю.

Deshewle · 14.03.2013, 11:05

Ребят, снова прошу помощи. Почему цистерну выворачивает?

[img]http://img152.**********.us/img152/9117/123dd.jpg[/img]
[img]http://img853.**********.us/img853/2483/1234th.jpg[/img]
[img]http://img818.**********.us/img818/6298/12345d.png[/img]

Прикладываю макрос:

Код:

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

/units, si 
/PREP7 
K, ,0,0,0, 
K, ,0,0,12, 
/REPLO 
L, 1, 2 
FLST,2,2,8 
FITEM,2,0,0,0 
FITEM,2,2,0,0 
CIRCLE,P51X, , , ,360, ,
/REPLO 
K, ,1.414213562,0,1.2, 
K, ,-1.414213562,0,1.2, 
K, ,-1.414213562,-2.5,1.2, 
K, ,1.414213562,-2.5,1.2, 
K, ,1.414213562,0,1.7, 
K, ,-1.414213562,0,1.7, 
K, ,-1.414213562,-2.5,1.7, 
K, ,1.414213562,-2.5,1.7, 
K, ,1.414213562,0,10.8, 
K, ,-1.414213562,0,10.8,
K, ,-1.414213562,-2.5,10.8, 
K, ,1.414213562,-2.5,10.8, 
K, ,1.414213562,0,10.3, 
K, ,-1.414213562,0,10.3,
K, ,-1.414213562,-2.5,10.3, 
K, ,1.414213562,-2.5,10.3, 
L, 8, 12 
L, 12, 11 
L, 11, 7 
L, 7, 8 
L, 12, 13 
L, 13, 9 
L, 9, 8 
L, 13, 14 
L, 14, 10 
L, 10, 9 
L, 14, 11 
L, 7, 10 
L, 16, 15 
L, 15, 19 
L, 19, 20 
L, 20, 16 
L, 16, 17 
L, 17, 21 
L, 21, 20 
L, 21, 18 
LDELE, 25 
L, 17, 18 
L, 18, 22 
L, 22, 21 
L, 18, 15 
L, 19, 22 
FLST,2,4,4 
FITEM,2,8 
FITEM,2,7 
FITEM,2,9 
FITEM,2,6 
AL,P51X 
FLST,2,4,4 
FITEM,2,6 
FITEM,2,10 
FITEM,2,12 
FITEM,2,11 
AL,P51X 
FLST,2,4,4 
FITEM,2,11 
FITEM,2,15 
FITEM,2,13 
FITEM,2,14 
AL,P51X 
FLST,2,4,4 
FITEM,2,8 
FITEM,2,16 
FITEM,2,17 
FITEM,2,14 
AL,P51X 
FLST,2,4,4 
FITEM,2,9 
FITEM,2,12 
FITEM,2,17 
FITEM,2,15 
AL,P51X 
FLST,2,4,4 
FITEM,2,7 
FITEM,2,10 
FITEM,2,16 
FITEM,2,13 
AL,P51X 
FLST,2,6,5,ORDE,2 
FITEM,2,1 
FITEM,2,-6 
VA,P51X 
FLST,2,4,4 
FITEM,2,21 
FITEM,2,24 
FITEM,2,22 
FITEM,2,23 
AL,P51X 
FLST,2,4,4 
FITEM,2,21 
FITEM,2,20 
FITEM,2,19 
FITEM,2,18 
AL,P51X 
FLST,2,4,4 
FITEM,2,23 
FITEM,2,25 
FITEM,2,27 
FITEM,2,26 
AL,P51X 
FLST,2,4,4 
FITEM,2,19 
FITEM,2,28 
FITEM,2,29 
FITEM,2,26 
AL,P51X 
FLST,2,4,4 
FITEM,2,24 
FITEM,2,27 
FITEM,2,20 
FITEM,2,29 
AL,P51X 
FLST,2,4,4 
FITEM,2,22 
FITEM,2,18 
FITEM,2,28 
FITEM,2,25 
AL,P51X 
FLST,2,6,5,ORDE,2 
FITEM,2,7 
FITEM,2,-12 
VA,P51X 
FLST,2,4,4,ORDE,2 
FITEM,2,2 
FITEM,2,-5 
ADRAG,P51X, , , , , , 1 
FLST,2,2,6,ORDE,2 
FITEM,2,1 
FITEM,2,-2 
FLST,2,2,6,ORDE,2 
FITEM,2,1 
FITEM,2,-2 
FLST,3,2,5,ORDE,2 
FITEM,3,15 
FITEM,3,-16 
VSBA,P51X,P51X, ,KEEP,KEEP 
FLST,2,2,6,ORDE,2 
FITEM,2,3 
FITEM,2,5 
VDELE,P51X, , ,1
FLST,2,2,6,ORDE,2 
FITEM,2,1 
FITEM,2,-2 
VDELE,P51X, , ,1
/REPLOT 
APLOT 
K, ,0,-3.5,0, 
K, ,0,-3.5,12, 
L, 7, 8 
L, 8, 2 
L, 2, 1 
L, 1, 7 
FLST,2,4,4 
FITEM,2,7 
FITEM,2,1 
FITEM,2,6 
FITEM,2,8 
AL,P51X 
FLST,2,2,6,ORDE,2 
FITEM,2,4 
FITEM,2,6 
FLST,3,1,5,ORDE,1 
FITEM,3,1 
VSBA,P51X,P51X, ,KEEP,KEEP 
ADELE, 1, , ,1
FLST,2,2,6,ORDE,2 
FITEM,2,1 
FITEM,2,3 
VDELE,P51X, , ,1
FLST,2,2,6,ORDE,2 
FITEM,2,4 
FITEM,2,6 
VDELE,P51X, , ,1
FLST,2,2,5,ORDE,2 
FITEM,2,14 
FITEM,2,-15 
ADELE,P51X, , ,1 
/REPLO 
!* 
ET,1,SHELL63
!* 
ET,2,SOLID45
!* 
R,1,0.003,0.003,0.003,0.003, , ,
RMORE, , , ,
RMORE 
RMORE, ,
!* 
!* 
MPTEMP,,,,,,,, 
MPTEMP,1,0 
MPDATA,EX,1,,2e11 
MPDATA,PRXY,1,,0.3 
FLST,5,2,6,ORDE,2 
FITEM,5,2 
FITEM,5,5 
CM,_Y,VOLU 
VSEL, , , ,P51X 
CM,_Y1,VOLU 
CMSEL,S,_Y 
!* 
CMSEL,S,_Y1 
VATT, 1, 1, 2, 0
CMSEL,S,_Y 
CMDELE,_Y 
CMDELE,_Y1 
!* 
FLST,5,2,5,ORDE,2 
FITEM,5,13 
FITEM,5,16 
CM,_Y,AREA 
ASEL, , , ,P51X 
CM,_Y1,AREA 
CMSEL,S,_Y 
!* 
CMSEL,S,_Y1 
AATT, 1, 1, 1, 0, 
CMSEL,S,_Y 
CMDELE,_Y 
CMDELE,_Y1 
!* 
SMRT,6 
SMRT,1 
TYPE, 2 
MAT, 1
REAL, 1 
ESYS, 0 
SECNUM, 
!* 
ESIZE,0.1,0,
FLST,2,14,5,ORDE,12 
FITEM,2,6 
FITEM,2,-8 
FITEM,2,10 
FITEM,2,13 
FITEM,2,16 
FITEM,2,18 
FITEM,2,20 
FITEM,2,-22 
FITEM,2,24 
FITEM,2,27 
FITEM,2,30 
FITEM,2,34 
AESIZE,P51X,0.1,
MSHAPE,0,3D 
MSHKEY,1
!* 
FLST,5,2,6,ORDE,2 
FITEM,5,2 
FITEM,5,5 
CM,_Y,VOLU 
VSEL, , , ,P51X 
CM,_Y1,VOLU 
CHKMSH,'VOLU' 
CMSEL,S,_Y 
!* 
VMESH,_Y1 
!* 
CMDELE,_Y 
CMDELE,_Y1 
CMDELE,_Y2 
!* 
TYPE, 1 
MAT, 1
REAL, 1 
ESYS, 0 
SECNUM, 
!* 
MSHAPE,0,2D 
!* 
FLST,5,2,5,ORDE,2 
FITEM,5,13 
FITEM,5,16 
CM,_Y,AREA 
ASEL, , , ,P51X 
CM,_Y1,AREA 
CHKMSH,'AREA' 
CMSEL,S,_Y 
!* 
AMESH,_Y1 
!* 
CMDELE,_Y 
CMDELE,_Y1 
CMDELE,_Y2 
!* 
/UI,MESH,OFF
FLST,5,3,5,ORDE,3 
FITEM,5,16 
FITEM,5,20 
FITEM,5,30 
ASEL,S, , ,P51X 
APLOT 
NSLA,S,1
NPLOT 
NUMMRG,NODE, , , ,LOW 
ALLSEL,ALL 
FLST,5,2,5,ORDE,2 
FITEM,5,13 
FITEM,5,16 
ASEL,S, , ,P51X 
NSLA,S,1
NPLOT 
SFGRAD,PRES,1,Y,90,-217.7777, 
FLST,2,7865,1,ORDE,86 
FITEM,2,7 
FITEM,2,18 
FITEM,2,-22 
FITEM,2,74 
FITEM,2,100 
FITEM,2,-114
FITEM,2,265 
FITEM,2,-344
FITEM,2,1231
FITEM,2,1242
FITEM,2,-1246 
FITEM,2,1313
FITEM,2,1324
FITEM,2,-1338 
FITEM,2,1489
FITEM,2,-1568 
FITEM,2,2449
FITEM,2,-6516 
FITEM,2,6523
FITEM,2,-6607 
FITEM,2,6614
FITEM,2,-6635 
FITEM,2,6642
FITEM,2,-6726 
FITEM,2,6733
FITEM,2,-6754 
FITEM,2,6761
FITEM,2,-6845 
FITEM,2,6852
FITEM,2,-6873 
FITEM,2,6880
FITEM,2,-6964 
FITEM,2,6971
FITEM,2,-6992 
FITEM,2,6999
FITEM,2,-7083 
FITEM,2,7090
FITEM,2,-7111 
FITEM,2,7118
FITEM,2,-7202 
FITEM,2,7209
FITEM,2,-7230 
FITEM,2,7237
FITEM,2,-7321 
FITEM,2,7328
FITEM,2,-7349 
FITEM,2,7356
FITEM,2,-7440 
FITEM,2,7447
FITEM,2,-7468 
FITEM,2,7475
FITEM,2,-7559 
FITEM,2,7566
FITEM,2,-7587 
FITEM,2,7594
FITEM,2,-7678 
FITEM,2,7685
FITEM,2,-7706 
FITEM,2,7713
FITEM,2,-7797 
FITEM,2,7804
FITEM,2,-7825 
FITEM,2,7832
FITEM,2,-7916 
FITEM,2,7923
FITEM,2,-7944 
FITEM,2,7951
FITEM,2,-8035 
FITEM,2,8042
FITEM,2,-8063 
FITEM,2,8070
FITEM,2,-8154 
FITEM,2,8161
FITEM,2,-8182 
FITEM,2,8189
FITEM,2,-8273 
FITEM,2,8280
FITEM,2,-8301 
FITEM,2,8308
FITEM,2,-8392 
FITEM,2,8399
FITEM,2,-8420 
FITEM,2,8427
FITEM,2,-8511 
FITEM,2,8518
FITEM,2,-10313 
/GO 
!* 
SF,P51X,PRES,0 
!* 
/PSF,PRES,NORM,2,1,1
/PBF,DEFA, ,1 
/PIC,DEFA, ,1 
/PSYMB,CS,0 
/PSYMB,NDIR,0 
/PSYMB,ESYS,0 
/PSYMB,LDIV,0 
/PSYMB,LDIR,0 
/PSYMB,ADIR,0 
/PSYMB,ECON,0 
/PSYMB,XNODE,0 
/PSYMB,DOT,1
/PSYMB,PCONV, 
/PSYMB,LAYR,0 
/PSYMB,FBCS,0 
!* 
/PBC,ALL, ,1
/REP
!* 
ALLSEL,ALL 
FLST,2,2,4,ORDE,2 
FITEM,2,32 
FITEM,2,36 
DL,P51X, ,SYMM 
FLST,2,2,5,ORDE,2 
FITEM,2,6 
FITEM,2,18 
DA,P51X,SYMM
FLST,2,2,5,ORDE,2 
FITEM,2,7 
FITEM,2,21 
!* 
/GO 
DA,P51X,ALL,
FINISH 
/SOL
/STATUS,SOLU
/STATUS,SOLU
SOLVE 
FINISH 
/POST1 
!* 
/EFACET,1 
PLNSOL, S,EQV, 0,1.0

AlexLPN · 14.03.2013, 11:46

Попробовал Ваш код запустить с командой NLGEOM,ON. Т.к. считаю, что тонкостенные конструкции без нее вообще смысла нет запускать (разве что в при совсем малых деформациях). Счет не пошел, разошелся. Но если считать последний результат, то там цистерна не выворачивается. Попробуйте может тут разобраться, что бы счет заработал.

Deshewle · 15.03.2013, 00:30

Непонятно, почему такие большие расхождения? Давление вроде правильно посчитано и приложено (p=pgh=1000*9.8*4=39200 Па)
Причем, если в цистерну помещаем жидкость с плотностью 1000 кг/(м^3) и считаем, то в ответе макисмальные напряжения всего 6 МПа, а не 30500 МПа, как в случае с гидростатическим давлением.

[IMG]http://img9.**********.us/img9/7800/123456ki.jpg[/IMG]

Кулик Алексей aka kpblc · 15.03.2013, 10:36

Deshewle, раз уж тебе в пост картинки не приложить, хотя бы скрывай их тэгом [more][/more]

AlexLPN · 15.03.2013, 14:25

Вы уверены, что давление правильно прикладываете? Там его нужно к площади относить, чтобы к весу перейти.

Deshewle · 15.03.2013, 15:03

Немного не понял о чем Вы =) В случае с градиентом давления, вроде как все логично приложено - ко всем узлам оболочки по линейному закону, в случае с жидкостью (fluid80), нагрузка от узлов жидкости полностью передается на узлы оболочки (shell63).

Deshewle · 16.03.2013, 23:15

Ребят, еще один вопрос. Решил заполнить цистерну наполовину жидкостью (водой). Связываю узлы сетки fluid и узлы сетки shell с помощью NUMMRG,NODE, , , ,LOW, получаю следующее (картинка ниже), выходит, что связываются лишь узлы линий shell и fluid, а дно - свободно.
Подскажите как еще можно грамотно связать (слить) 2 сетки?
CP? CPINTF?

[IMG]http://img708.**********.us/img708/2868/1222ql.jpg[/IMG]

MrWhite · 17.03.2013, 14:28

Цитата:

Сообщение от Deshewle

Подскажите как еще можно грамотно связать (слить) 2 сетки?

Уравнения связывания или контакт. Объединение узлов не подходит для разнотипных сеток.

Deshewle · 17.03.2013, 14:31

Контакт - отбрасываем (не нужна нелинейность). А можно поподробнее про уравнения связывания?

fasol · 17.03.2013, 20:45

Контакт у вас и не получится, контактные элементы не дружат с FLUID79-80.
Уравнение связи - команда CP, нужно объединить перемещения по нормали.

15.02.2013, 23:43		#1
Цистерна с гидростатическим давлением Deshewle Регистрация: 05.03.2012 Сообщений: 20 Всем привет. Считаю напряжения, возникающие в опорах цистерны. Напряжения возникают из-за гидростатического давления. (в цистерне вода). Посчитал с помощью градиаента давления. Задача решена линейно (контакт не захотел использовать, т.к. он мне пока не поддается), совместил узлы на опорах (solid45) и узлы на оболочке (shell63) и получил решение. (фото ниже) [img]http://img694.********.us/img694/5319/tank1w.jpg[/img] [img]http://img96.******.us/img96/1006/tank2mf.jpg[/img] Мне нужно проделать тоже самое, но без gradP, а, применив fluid80. [img]http://img812.******.us/img812/9001/tank4z.jpg[/img] [img]http://img41.******.us/img41/5167/tank3l.jpg[/img] Так же делаю одинаковую сетку на shell и на опорах (solid45), чтобы узлы были совмещены. Через Inertia--Gravity прикладываю ускорение свободно падения. В итоге ничего не получается, делаю solve, после чего отображаются только 4 опоры, а shell+fluid пропадают. В чем может быть проблема? Solution выдает следующее: [img]http://img163.*******.us/img163/6662/tank5o.jpg[/img] Что я делаю не так? Проблема с привязкой сеток? Последний раз редактировалось Кулик Алексей aka kpblc, 15.03.2013 в 10:36.*
Просмотров: 5996

Похожие темы
Тема	Автор	Раздел	Ответов	Последнее сообщение
Ищу Ост 24.023-29-84 Изготовление деталей из пластмасс прессованием и литьем под давлением. Типовой технологический процесс	elena_si	Поиск литературы, чертежей, моделей и прочих материалов	2	27.03.2012 18:33
Перечень документов для получения разрешения сдачи в эксплуатацию сосуда работающего под давлением.	Дмитрий DV	Инженерные сети	1	17.05.2011 11:53
Как правильно рассчитать толщину стенок закрытой емкости под давлением	burch	Конструкции зданий и сооружений	1	20.12.2009 23:05
Самый лёгкий материал для ёмкости работающей под высоким давлением?	Pavel Samofalov	Прочее. Отраслевые разделы	42	17.04.2009 14:01