Перспективными материалами для рабочих лопаток газовых турбин являются эвтектические композиты [34], в состав которых входят нитевидные монокристаллы и дисперсно-упрочненная монокристаллическая матрица.
Помимо обычных факторов, определяющих свойства композиций, таких как природа и соотношение компонентов, распределение их и др., термическое расширение композиционного материала в значительной степени определяется наличием в нем остаточных напряжений. В композиционном материале возникают напряжения, обусловленные различием коэффициентов теплового расширения компонентов. При некоторой температуре эти напряжения достигают предела текучести матрицы [35].
Для оценки коэффициента теплового расширения композита aкм в целом необходимо знать его составляющие: вдоль (α1) и поперек (α2) волокон в зависимости от коэффициентов теплового расширения волокон и матрицы αв, aм. Для изотропных волокон и матрицы такие соотношения приведены в работах [36, 37]. Для монокристаллов анизотропия a может достигать 300 % и более [38].
Задача определения αв, aм для монокристаллов решается, например, с учетом межатомного взаимодействия в кристаллической решетке. Так, коэффициент теплового расширения связан с температурой плавления и может характеризовать силу межатомной связи [38].
Максимальная сила межатомной связи приходится на металлы с максимальными температурами плавления и минимальными значениями коэффициента термического расширения. Соответственно, эти металлы имеют и минимальные межатомные расстояния.
Для металлов с кубической решеткой коэффициент a связан с температурой плавления соотношением . Кроме того, для веществ с металлическим характером связи имеется эмпирическое соотношение [38]. Однако эти соотношения не позволяют установить зависимость коэффициента теплового расширения от эксплуатационной температуры.
Расчет a, учитывающий ангармоническую составляющую колебаний атомов, по формуле [39]
где g и b – функции потенциала U; x – расстояние; Т – температура, предполагает знание потенциала кристаллической решетки U и также не имеет методики инженерного расчета коэффициента теплового расширения.
Для оценки коэффициента теплового расширения по формуле Гилмана [38] необходимо знать силу межатомного сцепления
где k – постоянная Больцмана; V – средний атомный объем.
Поэтому актуальной является задача вывода аналитической формулы и разработки инженерной методики расчета коэффициента теплового расширения применительно к монокристаллам материалов лопаток газовых турбин.
Модель и сравнение с экспериментом
Армирующей фазой в эвтектических композитах является нитевидный монокристалл. В нем рассматривается элементарная атомная ячейка. Кристаллическая решетка, находится в равновесном состоянии (силы притяжения уравновешены силами отталкивания). Потенциальная энергия взаимодействия зарядов U0 – величина постоянная, то есть решается статическая задача.
Расчет изменения периода кристаллической решетки вследствие повышения температуры основан на учете изменения энергии связи двух разноименных зарядов в атомной ячейке кристаллической решетки, т.е. между узлом решетки и электронным облаком в межузловом пространстве.
Энергия электростатической связи двух разноименных зарядов атома в кристаллической решетке в соответствии с законом Кулона определяется по формуле
где а0/2 – расстояние между разноименными зарядами (узлом кристаллической решетки и электронным облаком).
Тепловая энергия атома при нагреве ячейки монокристалла на величину DТ определяется как
где – теплоемкость атома при постоянном объеме. В соответствии с законом Дюлонга и Пти
Дж/К,
где NA – число Авогадро.
Энергия связи двух разноименных зарядов атома в ячейке кристаллической решетки вследствие нагрева уменьшается на величину
Изменение среднего расстояния между зарядами вследствие нагрева кристаллической решетки
Изменение периода кристаллической решетки вследствие нагрева
После преобразований зависимость изменения периода кристаллической решетки от температуры примет следующий вид [22]:
(1.12)
Вследствие периодичности и симметрии бездефектного монокристалла свойства элементарной атомной ячейки передаются кристаллической решетке. Поэтому для атомной ячейки и кристаллической решетки в целом коэффициент теплового расширения определяется как
(1.13)
Результаты расчета а0 в зависимости от температуры для некоторых металлов приведены в табл. 1.5.
Таблица 1.5
T, °С |
а0 для W |
а0эксп для W |
Расхождение, % |
а0 для Mo |
а0эксп для Mo |
Расхождение, % |
а0 для Cr |
а0эксп для Cr |
Расхождение, % |
а0 для Nb |
а0 эксп для Nb |
Расхождение, % |
а0 для V |
а0эксп для V |
Расхождение, % |
20 |
3,165 |
3,165 |
0,00 |
3,147 |
3,147 |
0,00 |
2,884 |
2,884 |
0,00 |
3,3 |
3,3 |
0,00 |
3,027 |
3,027 |
0,00 |
100 |
3,174 |
3,166 |
0,25 |
3,156 |
3,149 |
0,22 |
2,891 |
2,885 |
0,21 |
3,31 |
3,303 |
0,21 |
3,035 |
3,032 |
0,10 |
200 |
3,183 |
3,167 |
0,51 |
3,165 |
3,151 |
0,44 |
2,899 |
2,887 |
0,42 |
3,32 |
3,305 |
0,45 |
3,044 |
3,034 |
0,33 |
300 |
3,192 |
3,17 |
0,69 |
3,174 |
3,153 |
0,67 |
2,907 |
2,891 |
0,55 |
3,33 |
3,307 |
0,70 |
3,052 |
3,036 |
0,53 |
400 |
3,202 |
3,171 |
0,98 |
3,183 |
3,155 |
0,89 |
2,914 |
2,894 |
0,69 |
3,34 |
3,310 |
0,91 |
3,06 |
3,039 |
0,69 |
500 |
3,211 |
3,172 |
1,23 |
3,192 |
3,157 |
1,11 |
2,922 |
2,897 |
0,86 |
3,35 |
3,313 |
1,12 |
3,069 |
3,042 |
0,89 |
600 |
3,22 |
3,174 |
1,45 |
3,201 |
3,16 |
1,30 |
2,93 |
2,9 |
1,03 |
3,36 |
3,315 |
1,36 |
3,077 |
3,045 |
1,05 |
700 |
3,229 |
3,175 |
1,70 |
3,21 |
3,162 |
1,52 |
2,937 |
2,904 |
1,14 |
3,37 |
3,318 |
1,57 |
3,086 |
3,048 |
1,25 |
800 |
3,239 |
3,177 |
1,95 |
3,22 |
3,165 |
1,74 |
2,945 |
2,907 |
1,31 |
3,38 |
3,323 |
1,72 |
3,094 |
3,053 |
1,34 |
900 |
3,248 |
3,178 |
2,20 |
3,229 |
3,167 |
1,96 |
2,953 |
2,91 |
1,48 |
3,391 |
3,325 |
1,98 |
3,103 |
3,057 |
1,50 |
1000 |
3,258 |
3,18 |
2,45 |
3,238 |
3,169 |
2,18 |
2,961 |
2,915 |
1,58 |
3,401 |
3,327 |
2,22 |
3,112 |
3,060 |
1,70 |
1100 |
3,267 |
3,181 |
2,70 |
3,248 |
3,171 |
2,43 |
2,969 |
2,92 |
1,68 |
3,411 |
3,330 |
2,43 |
3,12 |
3,065 |
1,79 |
1200 |
3,277 |
3,183 |
2,95 |
3,257 |
3,174 |
2,61 |
2,976 |
2,922 |
1,85 |
3,422 |
3,333 |
2,67 |
3,129 |
3,070 |
1,92 |
1300 |
3,287 |
3,185 |
3,20 |
3,267 |
3,176 |
2,87 |
2,984 |
2,927 |
1,95 |
3,432 |
3,336 |
2,88 |
3,138 |
3,075 |
2,05 |
1400 |
3,296 |
3,188 |
3,39 |
3,277 |
3,178 |
3,12 |
2,992 |
2,934 |
1,98 |
3,443 |
3,340 |
3,08 |
3,147 |
3,079 |
2,21 |
1500 |
3,306 |
3,19 |
3,64 |
3,286 |
3,182 |
3,27 |
3,001 |
2,94 |
2,07 |
3,454 |
3,343 |
3,32 |
3,156 |
3,084 |
2,33 |
1600 |
- |
3,192 |
- |
- |
3,185 |
- |
- |
2,948 |
- |
- |
3,347 |
- |
- |
- |
- |
2000 |
- |
3,2 |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
Примечание. Величина периода кристаллической решётки а0 приведена в ангстремах. Экспериментальные значения а0эксп взяты
из работы [38].
Графическая зависимость периода кристаллической решетки от температуры, в частности, для хрома, приведена на рис. 1.10.
Рис. 1.10. Графическая зависимость периода кристаллической решетки
(в ангстремах) от температуры для хрома:
1 – эксперимент [38]; 2 – расчет
Из табл. 1.5 и рис. 1.10 видно, что результаты расчета а0 в зависимости от температуры удовлетворительно сходятся (с расхождением до 3,3 %) с экспериментальными данными работы [38]. На рис. 1.11 приведены экспериментальные зависимости [38].
Рис. 1.11. Экспериментальные зависимости [38]:
а – W; б – Mo; в – Cr; г – Nb; д –V
Результаты расчета периода кристаллической решетки в зависимости от температуры а0 в дальнейшем используются для определения:
– коэффициента теплового расширения [22];
– модуля упругости [5];
– предела упругости и др. характеристик [40].
В частности, в работе [5] приведены результаты расчета модуля упругости углеродной нанотрубки (УНТ). Уточненные значения для УНТ приведены в табл. 1.6.
Таблица 1.6
Расчетные значения модуля упругости УНТ в зависимости от температуры
t, °C |
20 |
100 |
200 |
300 |
400 |
500 |
600 |
700 |
800 |
900 |
E, ГПа |
1750 |
1731 |
1705 |
1671 |
1631 |
1584 |
1531 |
1473 |
1409 |
1342 |
Полученные формулы (1.12)–(1.13) позволяют определять теоретическую зависимость периода кристаллической решетки от температуры и значительно сократить объем дорогостоящих экспериментов.
Результаты расчета коэффициента теплового расширения армирующей фазы и матрицы эвтектического композита являются исходными данными для оценки термических напряжений в материале лопатки турбины.