Перейти к публикации

Модель полимера


Рекомендованные сообщения

Всем доброго дня!

На разрывной машине получены экспериментальные характеристики "сила-перемещение" нескольких цилиндрических образцов из полипропилена. Эксперимент проводился по ГОСТ 11262-80 при темп. 20 С.

 

5a743a9d54b89_--001.jpg.8d3534b9c3880a0e28c23425193c7f05.jpg

Затем получена диаграмма в координатах "напряжение-деформация" 

Stress-strain-001.jpg.d9a67c4f1fdbb2c1bfadfab27aa8d6bf.jpg

Нашел полученные характеристики похожие с (файл весь не загружается) 

Скрытый текст

Miller_Axel_ANSYS_2014_Detroit.jpg.c2e4152fdc4c57062db7a9f4ed4c3dbf.jpg

Далее пытаюсь подобрать свойства материала и смоделировать этот разрыв в ANSYS MAPDL.

Для этого на отрезке от 0 до σт свойства материла задаю упругими модуль упругости (MP,EX,1,600е6 МПа ) и коеф. Пуассона   (MP,PRXY,1,0.3)

Далее на отрезке от σт до σв задаю мультилинейное изотропное упрочнение (TB,PLASTIC,1,,,MISO) с вычетом упругих деформаций (значения здесь MISO.txt)

А вот ка быть далее? После σв  напряжения релаксируют, т.к. наблюдается снижение усилия на кривой. Этот процесс решил описать свойством ползучести (но...не полностью уверен).

В хелпе есть 13 моделей ползучести. Остановился на 2-ой модели с включенной опцией RATE,ON, да и примеры кое-какие есть здесь https://www.sharcnet.ca/Software/Fluent14/help/ans_str/Hlp_G_STR8_3.html#strmiscre6801 

Пробую рассчитать/подобрать коэффициенты ползучести при помощи Curve Fitting Creep Models так, как описано здесь https://www.sharcnet.ca/Software/Fluent14/help/ans_mat/Hlp_AM_CMCF.html#ncmcftable3 

но в итоге выдает ошибку error during solution. 

Вопрос: правильно ли я выбрал модель ползучести? и Как рассчитать/подобрать коэффициенты ползучести?

В планах провести такие же эксперименты только с учетом температуры.

Прошу вашего дельного совета в данном вопросе.

Всем откликнувшимся огромная благодарность!

 

Изменено пользователем Maximal
Ссылка на сообщение
Поделиться на других сайтах


7 минут назад, Maximal сказал:

Для этого на отрезке от 0 до σт свойства материла задаю упругими

так будет лучше: на отрезке от 0 до σпр 

7 минут назад, Maximal сказал:

После σв  напряжения релаксируют, т.к. наблюдается снижение усилия на кривой.

Смотрите на тру стресс. Не падают они. И не могут падать.

9 минут назад, Maximal сказал:

Прошу вашего дельного совета в данном вопросе.

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

Ссылка на сообщение
Поделиться на других сайтах
3 минуты назад, soklakov сказал:

так будет лучше: на отрезке от 0 до σпр 

да, я итак пробовал, но чтобы меньше точек было решил от 0 до σт

4 минуты назад, soklakov сказал:

Смотрите на тру стресс. Не падают они. И не могут падать.

это я вижу. анализировал экспериментальную диаграмму. усилие снижается...напряжения - нет.

что скажите на счет вот этого?

16 минут назад, Maximal сказал:

Вопрос: правильно ли я выбрал модель ползучести? и Как рассчитать/подобрать коэффициенты ползучести?

 

21 минуту назад, soklakov сказал:

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

этого я не заметил, сори

Ссылка на сообщение
Поделиться на других сайтах
23 минуты назад, soklakov сказал:

Но это сильно позже.

насколько позже?..)

24 минуты назад, soklakov сказал:

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

чтобы Вы посоветовали? 

Ссылка на сообщение
Поделиться на других сайтах
24 минуты назад, Maximal сказал:

усилие снижается...

потому что площадь уменьшается.

23 минуты назад, Maximal сказал:

чтобы Вы посоветовали?

Моделировать пока на просто пластичности. А вязкость вводить по мере несовпадения с экспериментом.

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

Ссылка на сообщение
Поделиться на других сайтах

Спасибо, это все я уже делал и вот что получил. это зависимость реакции (Н) в узле в заделке от перемещения (м)

polimer.jpg.5e47983ac13fec9e3c60a3fccb5c2deb.jpg

 

Ссылка на сообщение
Поделиться на других сайтах
3 minutes ago, Maximal said:

угу...как-то не очень)

поэтому нужно вязкость вводить

"У меня на лыжах с маленькой горки плохо получается, падаю.. Съеду-ка я лучше с большой, вдруг лучше получится!"

 

Остаётся только поддержать аплодисментами)

:clap_1:

Ссылка на сообщение
Поделиться на других сайтах

глаза боятся руки делают  - лыжи едут

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

Остаётся только поддержать аплодисментами)

и на этом спасибо 

Ссылка на сообщение
Поделиться на других сайтах
В 02.02.2018 в 16:45, Maximal сказал:

поэтому нужно вязкость вводить

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

Ссылка на сообщение
Поделиться на других сайтах
  • 2 недели спустя...
В 05.02.2018 в 11:29, soklakov сказал:

Попробуйте сначала параметры пластичности подобрать хотя бы чтобы совпадать по пиковому значению нагрузки и по соответствующему ему перемещению.

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

5a8821056f13c_-.thumb.jpg.67d982154ab09a22e1a4f0f52bc4e536.jpg

эксперименты с разными скоростями нагружения провел

вот, что получилось

5a8821b1a97c7_-.thumb.jpg.7d50470404041217eb650bc5fcafd850.jpg

 

5a8821b8ce184_-.thumb.jpg.c660c766d2a0583f5c7f1553293e8d5f.jpg

затем вязкость ввожу посредством команды PRONY https://www.sharcnet.ca/Software/Ansys/17.2/en-us/help/ans_str/Hlp_G_STR8_3.html#strviselastlm61199512

вот так

Скрытый текст

MU=0.3

Ro=910

 

E_1=24.1e9        ! рассчитано с эксперимента 
E_2=12.75E8       !это произвольное значение

NU_1=1.495E9    !это произвольное значение
NU_2=5.390E9    !это произвольное значение

 

G_1=E_1/(2*(1+MU)) 
G_2=E_2/(2*(1+MU))
K_1=E_1/(3*(1-2*MU))
K_2=E_2/(3*(1-2*MU))

E_0=E_1+E_2
G_0=G_1+G_2
K_0=K_1+K_2

A1=G_1/G_0        ! Relative modulus
A2=G_2/G_0

T1=NU_1/E_1        ! Relative time
T2=NU_2/E_2

MP,EX,1,E_0
MP,PRXY,1,MU
MP,DENS,1,RO

 

TB,PRONY,1,,II,SHEAR    ! Define Prony constants
*IF,II,EQ,1,THEN
    TBDATA,1,A1,T1
*ELSE
    TBDATA,1,0,T1,0,T2
*ENDIF

TB,PRONY,1,,II,BULK    ! Define Prony constants
*IF,II,EQ,1,THEN
    TBDATA,1,A1,T1
*ELSE
    TBDATA,1,0,T1,0,T2
*ENDIF

в итоге получается какая-то ерунда

polimer.jpg.7860ea7f59fcc5990468e1a899c1a99e.jpg

прошу вашей помощи направить мои действия в нужном направлении. 

всем спасибо!

Ссылка на сообщение
Поделиться на других сайтах
В 18.02.2018 в 00:12, Maximal сказал:

так-то получше будет

но вопросики остались....

win-win-win

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

А где вопросики?

Ссылка на сообщение
Поделиться на других сайтах

Вам нужно истинную диаграмму задавать. На участке после σв чтоб построить кривую нужно относительное сужение после разрыва знать. Рекомендую для лучшей сходимости применять все же закон кинематического упрочнения, с изотропным были задачи, где решение останавливалось, причем статика и были так же остановки с формулировкой типа "выходит за пределы теории расчета" (хотя сами по себе напряжения эквивалентные были на порядки меньше текучки). Судя по приведенной кривой вам подойдет. По приведенной зависимости можно достаточно точно построить.

 

 

Sigm.png

Изменено пользователем Chingachguk
Ссылка на сообщение
Поделиться на других сайтах

@Chingachguk , выложили бы уже всю статью тогда. А то еще думать, причем тут N1 и N2 и что такое отвлеченные единицы.

Ссылка на сообщение
Поделиться на других сайтах
10 часов назад, Chingachguk сказал:

Вам нужно истинную диаграмму задавать

согласен...я так и делал (см. пост №1). угол наклона кривой зависимости "сигма-епсилон" должен быть >= 0.

 

10 часов назад, Chingachguk сказал:

с изотропным были задачи, где решение останавливалось

да... пока такого не наблюдал.

спасибо за рекомендации.... и за ссылку :smile:

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

выложили бы уже всю статью тогда. А то еще думать, причем тут N1 и N2 и что такое отвлеченные единицы.

это дисер

23 часа назад, soklakov сказал:

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

как сделаю обязательно поделюсь 

23 часа назад, soklakov сказал:

А где вопросики?

потом и вопросики полезут.. как червячки:smile:

Ссылка на сообщение
Поделиться на других сайтах

@Maximal Можно лукавый вопрос: Ваше изделие из термопласта будет работать под нагрузкой менее 0,1 часа?

Ссылка на сообщение
Поделиться на других сайтах

Присоединяйтесь к обсуждению

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

Гость
Ответить в тему...

×   Вставлено в виде отформатированного текста.   Вставить в виде обычного текста

  Разрешено не более 75 эмодзи.

×   Ваша ссылка была автоматически встроена.   Отобразить как ссылку

×   Ваш предыдущий контент был восстановлен.   Очистить редактор

×   Вы не можете вставить изображения напрямую. Загрузите или вставьте изображения по ссылке.

  • Сейчас на странице   0 пользователей

    Нет пользователей, просматривающих эту страницу.




  • Сообщения

    • Eng.Destructor
      Сделал так: Добавил новый формат fs2 25 0 4 0 3t #No decimal, absolute, 3 поменял  result = newfs (25, thddepth$) # изменил 16 на 25 result = newfs (25, thdfirst$) # изменил 16 на 25 result = newfs (25, thdlast$) # изменил 16 на 25 result = newfs (25, thdfinish$) # изменил 16 на 25 Спасибо за наводку
    • arsenev
    • teqn1c
      Есть динозавр stratos plus на io mb -  В нем всего 256 КБ памяти  есть ли возможность через режим rtm подгружать программу с ПК через rs 232 ? В шкафу ещё есть разъем под карту но с нее получается только считать и опять таки сохранить на память станка. Соединение ПК-станок настроил ,есть возможность отправлять на него через альфакам программы ,но они сохраняются в память станка. В альфакам есть непрерывная отправка - она отправляет - но в режиме rtm ничего не происходит 
    • Metal_Cutter
      @arsenev В FC такого нет.
    • Bot
      HCL CAMWorks 2024, новая версия программного обеспечения CAM, полностью интегрированного с SOLIDWORKS и Solid Edge включает в себя передовые технологии, позволяющие сократить время программирования, уменьшить время цикла обработки и увеличить срок службы инструмента. Смотреть полностью
    • Александр1979
      Руководство для дусь взял здесь: https://www.cnczone.com/forums/daewoo-doosan/413426-every-doosan-service-manual-currently-auto-updated.html
    • ss_162
      Может есть у кого оригинальные резцы на продажу? Нужен резец SFCC 50 - это маркировка D`Andrea, или Искаровский, он такой же, только маркировка у него другая.
    • Aveovar
      Подскажи, как ты сделал несколько экранов с разных ракурсов
    • Kaktusito
      Всем доброго времени суток! Только начал заниматься редактированием постпроцессоров и сразу же столкнулся с проблемой. Суть в том что я нашел по форумам постпроцессоры с SolidCAM 2015 года, ибо у меня в версии 2022 года их почти нет, а мне нужен был более менее адекватный пост под токарку.  Закинул все файлы Gpptool 2015 года в свою папку, выбрал подходящий мне станок и начал генерировать УП.  Как вы поняли ничего мне не сгенерировало. Пишет ошибку "please change the name "Bore" in the VMID. This name is reserved for use in the GPP"  Сделал то что просит ошибка и был уверен что это не поможет, так оно и вышло. По итогу начало долбить другие ошибки. Очень еще не понятно с файлами которые я кидал в папку Gpptool потому что вместо стандартных двух у меня их аж 4 штуки.: gpp с циклами, gpp обычный, PRP и vmid. Хотя я знаю что в версиях вроде после 2017 уже идет по 2 файла. Что я только не пробовал кидал только 2 файла, кидал все 4, кидал содержимое файла с циклами в основной gpp, шёл по пути ошибок и исправлял что меня просили. В общем ничего не помогает  От безисходности пишу сюда, вы моя последняя надежда понять что происходит и почему ничего не работает. Заранее благодарю за ответы!
    • Krusnik
×
×
  • Создать...