Какие CAD умеют решать задачи термоупругости?

[soklakov 1 475]

Но в стали нет полимерных цепочек.

[Fedor 1 597]

Нет. Что там с энтропией происходит не понятно при растяжении. Хотя проволочка при сгибании- разгибании нагревается за счет внутреннего  трения скорее. То есть энтропия вроде уменьшается... :)

[Борман 2 375]

Если принять какие-то гипотезы, то задача распадается на тепловую и упругую.

Нам это выводили на парах. Так оно и было.

[AlexKaz 1 822]

Связанную вроде никакие не решают... :)

Нет. Что там с энтропией происходит не понятно при растяжении. Хотя проволочка при сгибании- разгибании нагревается за счет внутреннего трения скорее. То есть энтропия вроде уменьшается... :)

Уравнение (8) - оно?

В таком виде код прикручивали к Abaqus вроде асперы в ИМСС.

[Fedor 1 597]

Извините наизусть не помню. Вроде у Новацкого в Теории упругости было. http://pnu.edu.ru/media/filer_public/2013/04/10/5-2_novacky_1975.pdfСмотрите сами - стр. 88-89 :)

Изменено 26 апреля 2016 пользователем Fedor

[AlexKaz 1 822]

В 22.04.2016 в 16:24, Fedor сказал:

Да, это примерно и имел ввиду, только с точностью до наоборот. Области растяжения становятся горячее, а не как у них холоднее. Любой может убедиться в этом на надувном шарике. С физиками такое бывает, когда в механики подаются  :)

Хитрая задачка. Если считать через энтропию при адиабатическом процессе, то температура убывает. Если считать через "классическую связанность полей деформации и температуры" - растёт. Чёт не понятно.

 

Уравнение "классической связанности" (Бьюи стр. 132): (lambda+mu) * grad (e_ii) + mu * laplacian u - 3 * K * alpha * grad T = 0 => здесь видно, что увеличение деформаций e_ii даёт прирост T.

 

Уравнение энтропии при адиабате (Бьюи стр. 138+-):           

eq1: 3*K*alpha*eii+c*dT/T0=0; =>  dT = -3*alpha*e_ii*K*T0/c

 

Cталь:

lambda=mu=0.8*10^11==8e10 Н/м^2;
c = 462 Дж/(кг*K)- теплоёмкость;
alpha = 1.5*10^-5==1.5e-5 - к-т теплового объёмного расширения;
K=(3*lambda+2*mu)/3==1.33e11 Па; точнее K=E/(3*(1-2*v))=1.75e11 Па
T0=273.15 K - начальная температура

=> приращение температуры: dT = -3547402.6 * e_ii

 

Резина:

K=E/(3*(1-2*v))), alpha=77e-6,E=5e6, v=0.46,c=1420,T0=273.15

=> приращение температуры: dT = -925.73 * e_ii

 

Изменено 12 марта 2018 пользователем AlexKaz

[Fedor 1 597]

"увеличение деформаций e_ii даёт прирост T"  это вроде похоже на сжимаемость. А у резины она близка к нулю. То и градиент шаровой части тензора напряжений маленький.  Там из - за девиатора должна меняться наверное...

Изменено 13 марта 2018 пользователем Fedor