Снова помощь по результатам.

[Happybulkin 0]

Добрый день уважаемые посетители. Снова вынужден обратиться за помощью- старая тема открытая мной называлась "Помощь по результатам". Если вкратце- тонкостенная трубка нагружена внутренним давлением, зафиксирован только один торец от осевых перемещений, интересует поведение трубки за пределами упругости. Вначале задача решалась, НО результаты по радиальным перемещениям получались в 30 раз меньше чем на практике, в процессе разборок что то пытался менять. Теперь же задача просто разваливается в самом начале перехода в пластику, уменьшать сетку и шаги пробовал. Моих знаний недостаточно чтобы определить причину несходимости. Прошу у Вас помощи!! задачу решить мне необходимо, а возможности нет. По необходимости прикреплю файл. P.S. Особенно поможет если кто то вызовется помочь мне досканально разобраться в причинах, не за спасибо конечно же. Заранее благодарю.

[piden 2 874]

Выложи. Или хоть условие задачи опиши. Только полностью, а не так, как в предыдущих темах..

[Happybulkin 0]

Прикрепил фото экрана(файл не прикрепляется- 5мб весит). Задача осесимметричная. На внутренней стенке давление(из рисунка видно). Давление увеличиваю до максимального за 500 шагов по времени. Размеры трубки:  длина 0,15, диаметр 0.03; толщина стенки;0,002- все в метрах. Число элементов 2000. Кривая деформирования полилинейная по точкам для 09г2с. Материал с изотропным упрочнением по мизесу. В солюшен контролс ставлю large displacement static, солвер ставлю спарс директ, также подключаю полный ньютон рафсон(full n- r), также настраиваю критерии сходимости по силе и перемещению(взял из хелпа). Задача не сходится как только начинается переход в пластику(файл err прикрепил). А на границе пластики перемещения получаются в 30 раз меньше реальных.

obraztsy truba polylin osnova.err

[piden 2 874]

 *** ERROR ***                           CP =       9.688   TIME= 21:21:29
 For material 1 at temperature 0, the initial slope of the stress-strain 
 curve based on the TB,MISO table ( 2.102564103E+11 ) is greater than    
 EX (2.1E+11).

 

......

 

 *** ERROR ***     SUPPRESSED MESSAGE    CP =    2416.547   TIME= 21:54:31
 One or more elements have become highly distorted.   

 

Что-то с "полилинейной кривой деформирования для 09г2с" не то...

И еще сетка сильно деформируется, вырождаются элементы.

 

@Happybulkin, скажите, вам когда задачу ставили, тоже так сумбурно все излагали? Говорили, мол: смотри, давление как на картинке, свойства как у ... 9г2с, что, зачем, для чего - не важно. Считай, дорогой. Так было?

 

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

Ну а когда "полный ньютон рафсон - поцелуй кошку в голову", то разбираться не очч охота..

 

Изменено 18 января 2017 пользователем piden

[Борман 2 375]

Кольцевые напряжения получаются 225МПа.

Маловато для пластики в 09Г2С. А раз анасису кажется что не маловато, то ВЫВОД: Кривая кривая деформирования.

И да... осесимметричность то включена ???

Какая у вас последняя точка на диаграмме деформирования ?

[Happybulkin 0]

22 часа назад, piden сказал:

 *** ERROR ***                           CP =       9.688   TIME= 21:21:29
 For material 1 at temperature 0, the initial slope of the stress-strain 
 curve based on the TB,MISO table ( 2.102564103E+11 ) is greater than    
 EX (2.1E+11).

 

......

 

 *** ERROR ***     SUPPRESSED MESSAGE    CP =    2416.547   TIME= 21:54:31
 One or more elements have become highly distorted.   

 

Что-то с "полилинейной кривой деформирования для 09г2с" не то...

И еще сетка сильно деформируется, вырождаются элементы.

 

@Happybulkin, скажите, вам когда задачу ставили, тоже так сумбурно все излагали? Говорили, мол: смотри, давление как на картинке, свойства как у ... 9г2с, что, зачем, для чего - не важно. Считай, дорогой. Так было?

 

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

Ну а когда "полный ньютон рафсон - поцелуй кошку в голову", то разбираться не очч охота..

 

Piden, спасибо что ответили. А чтобы элементы не вырождались их надо больше сделать?

ERROR связанный с Ех я исправил- просто не знал что в файле можно исправлять чтобы информация полной была.

 

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

Задача мне ставилась без учета того что решаться она будет в ансисе и звучала следующим образом: смоделировать эксперимент по нагружению тонкостенной трубки внутренним давлением(прикладываю рисунок установки для похожего эксперимента- образец на ней номер 3) и детально определить ее поведение на всем протяжении эксперимента вплоть до разрушения трубки. Для начала я попробовал упростить и смоделировать поведение трубки с незакрепленными концами- собственно та задача которую я и решал.

 

Прикладываю диаграмму деформирования которую я закладывал, модуль юнга 2.1Е11. Элементы plane 183. Давление в трубке задавал 50 МПа. 

22 часа назад, Борман сказал:

Кольцевые напряжения получаются 225МПа.

Маловато для пластики в 09Г2С. А раз анасису кажется что не маловато, то ВЫВОД: Кривая кривая деформирования.

И да... осесимметричность то включена ???

Какая у вас последняя точка на диаграмме деформирования ?

Спасибо. Осесимметричность я включил в элементе(обведено красным на прилагаемом рисунке). Кривую деформирования после упругости- в полном объеме я выложил в прошлом сообщении. Давление записывал в расчете что это Паскали а все размеры в метрах.

[Борман 2 375]

Без модели не раскурить.

[Happybulkin 0]

В 20.01.2017 в 07:04, Борман сказал:

Без модели не раскурить.

Может можно модель вам прислать по почте или в личных сообщениях? Просто модель у меня больше 1 мб весит а к сообщению на форуме максимум один мегабайт можно прикреплять.

[AlexKaz 1 821]

16 часов назад, Happybulkin сказал:

модель у меня больше 1 мб весит

Кидайте ADPL-код, подчищенный log-файл или проект по методике Борман'а (заархивируйте потом):

Цитата

ALLSEL,ALL

CDWRITE,DB,TO_BORMAN,cdb,,,

ALLSEL,ALL

А вообще, самый короткий анекдот

Цитата

UNZIP.ZIP

lol

 

[kolo666 253]

В 19.01.2017 в 21:36, Happybulkin сказал:

ERROR связанный с Ех я исправил

Кривая stress-strain должна идти из точки (0,0), точно исправили? Не вижу на скриншоте..

Вторая точка лежит на прямой упругого участка, ну а далее остальное, если речь конечно идет о кривой именно stress-strein, а не stress-plastiс strain.

 

P.S. Советую для начала на кубике посчитать мультилинейную модуль материала, и по результатам расчета на растяжение построить эту кривую, и сравнить с тем что бы вы хотели там увидеть..а уж потом решать осесимитричные задачи

Изменено 27 января 2017 пользователем kolo666

[kolo666 253]

@Happybulkin  Не подскажите заодно, откуда взяли точки для кривой деформирования для 09г2с? Мне бы пригодились =)

[piden 2 874]

12 hours ago, kolo666 said:

Кривая stress-strain должна идти из точки (0,0)

По-моему, этого как раз и не нужно делать..

[kolo666 253]

19 минут назад, piden сказал:

По-моему, этого как раз и не нужно делать..

Это почему же? 

[piden 2 874]

Я не могу сейчас найти цитату из документации, но суть в том, что прога и сама понимает, что кривая выходит из точки (0, 0).

И первой точкой будет (отн. деформации при напряжении текучести, напряжение текучести)

 

!==============================================================
! Multilinear Kinematic Hardening with Stress vs. Total Strain
!==============================================================

/prep7
TB,KINH,1,2,3                  ! Activate a data table
TBTEMP,20.0                    ! Temperature = 20.0
TBPT,,0.001,1.0                ! Strain = 0.001, Stress = 1.0
TBPT,,0.1012,1.2               ! Strain = 0.1012, Stress = 1.2
TBPT,,0.2013,1.3               ! Strain = 0.2013, Stress = 1.3
TBTEMP,40.0                    ! Temperature = 40.0
TBPT,,0.008,0.9                ! Strain = 0.008, Stress = 0.9
TBPT,,0.09088,1.0              ! Strain = 0.09088, Stress = 1.0
TBPT,,0.12926,1.05             ! Strain = 0.12926, Stress = 1.05

 

[kolo666 253]

1 минуту назад, piden сказал:

Я не могу сейчас найти цитату из документации, но суть в том, что прога и сама понимает, что кривая выходит из точки (0, 0).

И первой точкой будет (отн. деформации при напряжении текучести, напряжение текучести)

 

 

!==============================================================
! Multilinear Kinematic Hardening with Stress vs. Total Strain
!==============================================================

/prep7
TB,KINH,1,2,3                  ! Activate a data table
TBTEMP,20.0                    ! Temperature = 20.0
TBPT,,0.001,1.0                ! Strain = 0.001, Stress = 1.0
TBPT,,0.1012,1.2               ! Strain = 0.1012, Stress = 1.2
TBPT,,0.2013,1.3               ! Strain = 0.2013, Stress = 1.3
TBTEMP,40.0                    ! Temperature = 40.0
TBPT,,0.008,0.9                ! Strain = 0.008, Stress = 0.9
TBPT,,0.09088,1.0              ! Strain = 0.09088, Stress = 1.0
TBPT,,0.12926,1.05             ! Strain = 0.12926, Stress = 1.05

 

 

Все понял, действительно так..хотя Nastran например требует ее..

[Борман 2 375]

8 часов назад, piden сказал:

 

Я не могу сейчас найти цитату из документации, но суть в том, что прога и сама понимает, что кривая выходит из точки (0, 0).

И первой точкой будет (отн. деформации при напряжении текучести, напряжение текучести)

 

Поскучал...

Цитата

The multilinear hardening behavior is described by a piece-wise linear stress-total strain curve, starting at the origin and defined by sets of positive stress and strain values, as shown in this figure:

...

The first stress-strain point corresponds to the yield stress. 

Как жить ?

[Happybulkin 0]

В 23.01.2017 в 14:14, AlexKaz сказал:

Кидайте ADPL-код, подчищенный log-файл или проект по методике Борман'а (заархивируйте потом):

А вообще, самый короткий анекдот

lol

 

Прикрепляю файл текстовый с кодом. 

Obraztsy truba.log

[Борман 2 375]

43 минуты назад, Happybulkin сказал:

Obraztsy truba.log

Не вижу в модели никаких проблем.

На ваши настройки решателя смотреть не стал.

Чтобы решить конкретно это задачу, после /sol

1) Закрепить

2) Нагрузить

3) nsubst,10

4) autots,on

5) Решить.

Эквивалентные напряжения равны кольцевым и совпадают с теоретическими.

Пластики нет.

 

В чем проблема ?

 

[Happybulkin 0]

Добрый день всем! Снова возникли проблемы- При нагружении этой же задачи за предел текучести задача в скором времени(подшагов через 10 после начала пластики) перестает сходится и соответственно останавливает решение. В комментариях пишет что часть элементов выродилась и соответственно счет не сходится. Пробовал увеличивать количество шагов, увеличивать количество элементов(до нескольких десятков тысяч), менял на треугольные элементы- счет все равно не сходится в самом начале пластики. Может кто то подсказать как быть в этой ситуации? Скидываю код. 

 

/BATCH  
/COM,ANSYS RELEASE 14.0    UP20111024       20:50:00    01/29/2017
/input,menust,tmp,'',,,,,,,,,,,,,,,,1   
/GRA,POWER
/GST,ON
/PLO,INFO,3
/GRO,CURL,ON
/CPLANE,1   
/REPLOT,RESIZE  
WPSTYLE,,,,,,,,0
!*  
/NOPR   
KEYW,PR_SET,1   
KEYW,PR_STRUC,1 
KEYW,PR_THERM,0 
KEYW,PR_FLUID,0 
KEYW,PR_ELMAG,0 
KEYW,MAGNOD,0   
KEYW,MAGEDG,0   
KEYW,MAGHFE,0   
KEYW,MAGELC,0   
KEYW,PR_MULTI,0 
KEYW,PR_CFD,0   
/GO 
!*  
/COM,   
/COM,Preferences for GUI filtering have been set to display:
/COM,  Structural   
!*  
/PREP7  
!*  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,2.1E11 
MPDATA,PRXY,1,,0.3  
TB,MISO,1,1,7,0 
TBTEMP,0
TBPT,,0.00195,4.095E8   
TBPT,,0.022,4.1E8   
TBPT,,0.042,4.75E8  
TBPT,,0.0875,5.75E8 
TBPT,,0.175,7.25E8  
TBPT,,0.25,7.75E8   
TBPT,,1,9.75E8  
!*  
ET,1,PLANE183   
!*  
KEYOPT,1,1,0
KEYOPT,1,3,1
KEYOPT,1,6,0
!*  
K,1,0.015,0,0,  
K,2,0.017,0,0,  
K,3,0.017,0.3,0,
K,3,0.015,0.3,0,
KDELE,       3  
K,3,0.017,0.3,0,
K,4,0.015,0.3,0,
LSTR,       1,       2  
LSTR,       2,       3  
LSTR,       3,       4  
LSTR,       4,       1  
ESIZE,0.001,0,  
FLST,2,4,4  
FITEM,2,2   
FITEM,2,1   
FITEM,2,4   
FITEM,2,3   
AL,P51X 
MSHAPE,0,2D 
MSHKEY,0
!*  
CM,_Y,AREA  
ASEL, , , ,       1 
CM,_Y1,AREA 
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
FINISH  
/SOL
!*  
CNVTOL,U, ,0.001,2,0.0001,  
!*  
CNVTOL,F, ,0.001,2,1,   
!*  
ANTYPE,0
NLGEOM,1
NSUBST,500,500,300  
EQSLV,SPAR  
TIME,1  
!*  
!*  
NLGEOM,1
NROPT,FULL, ,   
STAOPT,DEFA 
LUMPM,0 
EQSLV,SPAR, ,0, ,DELE   
MSAVE,0 
PCGOPT,0, ,AUTO, , ,AUTO
PIVCHECK,1  
PSTRESS,0   
TOFFST,0,   
!*  
FLST,2,1,4,ORDE,1   
FITEM,2,1   
!*  
/GO 
DL,P51X, ,UY,   
FLST,2,1,4,ORDE,1   
FITEM,2,4   
/GO 
FLST,2,1,4,ORDE,1   
FITEM,2,4   
/GO 
!*  
SFL,P51X,PRES,60000000, 
/STATUS,SOLU

Изменено 13 апреля 2017 пользователем Happybulkin

[AlexKaz 1 821]

@Happybulkin , удалите косячную вторую точку из кривой sigma-epsilon.

Цитата

/CLEAR,START
/PREP7
MPTEMP,,,,,,,,
MPTEMP,1,0  
MPDATA,EX,1,,2.1E11
MPDATA,PRXY,1,,0.3  
TB,MISO,1,1,7,0
TBTEMP,0
TBPT,,0.00195,4.095E8
TBPT,,0.042,4.75E8  
TBPT,,0.0875,5.75E8
TBPT,,0.175,7.25E8  
TBPT,,0.25,7.75E8   
TBPT,,1,9.75E8  
!*  
ET,1,PLANE183   
!*  
KEYOPT,1,1,0
KEYOPT,1,3,1
KEYOPT,1,6,0
!*  
K,1,0.015,0,0,  
K,2,0.017,0,0,  
K,3,0.017,0.3,0,
K,3,0.015,0.3,0,
KDELE,       3  
K,3,0.017,0.3,0,
K,4,0.015,0.3,0,
LSTR,       1,       2  
LSTR,       2,       3  
LSTR,       3,       4  
LSTR,       4,       1  
ESIZE,0.001/4,0,  
lsel,all
AL,all
MSHAPE,0,2D
MSHKEY,0
ASEL, , , , 1
AMESH,1

FINISH  
/SOL

 

ANTYPE,0
NLGEOM,1
NSUBST,10  !500,500,300
autots,on

DL,1, ,UY,
SFL,4,PRES,60000000,
/solu
solve