Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Температурные напряжения и малоцикловая усталость

..pdf
Скачиваний:
11
Добавлен:
12.11.2023
Размер:
20.99 Mб
Скачать

равновесия, и определить из них необходимый ряд, удовлетво­ ряющий уравнениям совместности. Однако правильность выбора такого ряда может быть определена расчетом с помощью дополни­ тельной энергии тела, которая выражена через напряжения. Принцип минимума дополнительной энергии заключается в том, что из всей серии функций, удовлетворяющих уравнениям равно­ весия и согласующихся с л юбыми заданными напряжениями на поверхности, правильной является такая, которая сводит допол­ нительную энергию тела к минимуму. Можно показать, что при использовании вариационных методов выполнение условия ми­ нимума дополнительной энергии эквивалентно удовлетворению уравнений совместности.

1.7.4. Применение принципа минимума дополнительной энер­ гии к задачам о температурных напряжениях. Преимущество ис­ пользования принципа дополнительной энергии заключается в том, что он требует выбора системы функций, автоматически удовлетворяющих уравнению равновесия. Такой функцией яв­ ляется, например, функция напряжений Эри. Если соотношения Эри

_д2(р

__ д2ср

_

д2ф

(1.64)

ду2 ’

дх2 9

 

дхду

 

 

подставить в выражение (1.62), то результат будет выражен ис­ ключительно через ф и полученное уравнение может быть раз­ решено относительно ф, если найти минимум дополнительной энергии, следуя обычному вариационному методу.

К сожалению, если ф остается в общем виде неизвестной функ­ цией х и у, то вариационный метод дает уже известное бигармоническое уравнение У 4ф = —Е а у 2Т, что соответствует эквивалент­ ности принципа дополнительной энергии и уравнений совмест­ ности.

Достоинство энергетического метода заключено в получении приближенного решения по методам, подобным методам Релея и Ритца, используемым в других случаях численного анализа. Для решения принимается такая форма функции напряжений ф, в которую входят как известные, так и неизвестные функции. Применение теоремы о дополнительной энергии позволяет соста­ вить уравнения для определения неизвестных функций в форме, сходной с принятыми известными функциями. Если выбранные известные функции соответствуют точному решению, то неизвест­ ные также дадут точное решение. Если функции выбраны не­ точно, то решение окажется приближенным, причем степень приближения зависит от того, насколько правильно сделан выбор.

Подставляя уравнения (1.64) в уравнение (1.62), получаем дополнительную энергию

11=1] + Ш Т + 2 Е а Г (&+Ш

('-в5»

Примем, что ф выражается как сумма произведений функций от х на функции от у, причем функции от у известны, а функции от х неизвестны:

тогда

Ф =

Л (У) Фх (х) +

Р2 (У) ф2 (*) +

рз (У) Фз (*);

О -66)

 

 

 

 

 

0

=

*>х(У) % (х) +

Р2 (у) ф" (х) +

Р3 (У) Фз (х);

(1 -67)

0

=

Р\ (У) Фх(х) +

Р; (у) ф2 (х) +

Р'' (у) фз (х).

(1.68)

Так кад Р х (у), Р 2 (у), Ра )—известные функции от у, то продиф­ ференцировать можно точно. После подстановки уравнений (1.67) и (1.68) в уравнение (1.65) дополнительная энергия будет выра­ жена через различные комбинации ф и <р", и никаких других про­ изводных не будет. Проинтегрировать можно сначала по у, так как все выражения через у известны, при этом дополнительная энер­ гия выражается в виде интеграла от различных функций ф и их вторых производных:

(1.69)

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

(1.70)

причем уравнений столько, сколько имеется различных функций ф. Таким образом, вариационный процесс приводит к числу диффе­ ренциальных уравнений, равному числу неизвестных. Совместное решение этих уравнений относительно фь ф2. соответствует оптимальному выбору этих функций для принятых значений Р х (у), Р 2 (у) и т. д. Таким образом, нужно решить, сколько следует взять функций ф и какой должен быть вид функций Р (у).

1.8. МЕТОД ХЕЛДЕНФЕЛСА И РОБЕРТСА

1.8.1.Основное уравнение. Р. Хелденфелс и В. Робертс [1.14]

применили этот метод к расчету

температурных

напряжений

в длинной плоской

пластине, ограниченной х = ± а ,

у = ± Ь =

= ±1 и имеющей

распределение

температур

 

Т = Т0 + X (х) V (у).

(1.71)

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

Ф = Н х)ё(у).

(1.72)

Тогда по уравнению (1.64)

<Ъ = ё"!\ °у = ёГ\

= — !'ё',

(1-73)

где штрихами показано дифференцирование функций (/ или §) по соответствующим переменным или у).

В задаче, решенной аналитически и экспериментально Хелденфельсом и Робертсом, рассматривалась длинная плоская пластина,

вкоторой температура изменялась только по меньшему размеру у

ибыла постоянна по большему размеру х. Как уже отмечалось (см. раздел 1.6), распределение напряжений в такой балке (в дан­ ном случае в пластине) может быть легко определено по аналогии Дюамеля. Эти напряжения определяются в областях, достаточно удаленных от краев, где нормальные и касательные напряжения должны снижаться до нуля. Следовательно, только вблизи краев требуется строгое решение, основанное на использовании функции напряжений. Так как на большей части пластины напряжения из­ вестны, форма функции напряжений выбирается такой, чтобы

напряжения, полученные в этой области, совпадали с известным решением. По уравнению (1.45)

-[-1

-|-1

§" = оХа> = - а Е [ Т - 1/й \

Т й у - Ц - \Т у й у ), (1.74)

—1

-1

где Т = Т (у) — функция только от у.

Можно ожидать, что функция / (х) оказывается равной единице в большей части пластины (тогда сА = [§" — §") и только на краях будет отличаться от единицы. Для получения выражения для § через значения двух произвольных констант, определяемых из граничных условий вдоль длинных краев, уравнение (1.74) можно проинтегрировать дважды. Тргда

и

 

 

йу Лу +

Вгу В2,

(1-75)

где В г и В 2 определяются в каждом

отдельном случае так, что

§ = §' = 0 при у = 1. Хотя выбор формы для § (у) основывается на рассмотрении длинной балки с изменением температуры только в направлении меньшего размера у, метод может быть обобщен для случая, когда температура изменяется в обоих направлениях по уравнению (1.71), если вместо Т в уравнение (1.74) подставить компонент 7 (у) уравнения (1.71). Подставив выражение для §

в уравнение для дополнительной энергии и интегрируя по у, так как все выражения через у являются явными, получаем

Ц = т

[А^ +

А Г +

2Лз/Г + .2Я а/ (Л6 +

Л „*) +

где

+

2ЕаГ(А1 +

АгХ)]4х,

(1.76)

 

+ ь

 

 

^

 

 

 

 

 

 

 

> » 1 = [

§"2ау,

 

 

 

—ь

 

 

 

 

 

 

-\-Ь

 

 

 

 

 

А2 = \ ё2Лу,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аз = ё ё '

— А ;

 

 

 

'

 

 

 

 

 

У

 

ё п &У\

 

 

 

—ь

 

 

 

 

 

+6

 

(1.77)

 

 

 

 

 

 

Л

= П

г ' |

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

—Ь

УдГйу,

 

 

 

 

 

 

 

 

А7== Т'о

У| ё &У\

 

 

 

 

 

 

 

 

 

 

 

 

 

 

АВ = \

Уе йу-

 

 

 

ь

 

 

 

 

1.8.2. Решения уравнения. Применяя вариационные уравне­ ния (1.70) к выражению для дополнительной энергии (1.76) и.зная, что при у = ±Ь, тху — ['д' = 0, т. е. §' = 0, а также полагая Л3 = —Л4 и Лв = 0, получаем результирующее уравнение

Л*Г" - 2Л4Г + АЛ= - Еа (АаХ + А8Х").

(1.78)

Уравнение (1.78) — это обычное дифференциальное уравнение с постоянными коэффициентами, решение которого может быть получено в окончательном виде. Решение однородного уравнения дает

= Сх соз Нрх соз дх + С2 соз Нрх з т

дх +

+ С3 з т Нрх соз дх + С4 з т Нрх з т дх,

(1.79)

где р и д — действительная и мнимая части корней ±

(р + 1д)

и ± (р — 1д) характеристического уравнения

 

Л 2т 4 — 2Л4т 2 + А х = 0.

(1.80)

Частное решение уравнения (1.78) может быть найдено из до­ пущения, что решение будет иметь такую же форму, что и X,

если функция X

выражена в простом виде — полиномиальной,

экспоненциальной или гармонической функцией. Для обеспечения

равенства правой и левой частей уравнения вводится достаточное

число неопределенных коэффициентов. Методика дана ниже для

случая, когда температура является функцией только У (X =

1).

Она существенно не меняется, если X задана такой функцией,

для которой частное решение легко может быть найдено. Частное

решение уравнения (1.78) для X =

1 можно получить при условии

/ = соп$1

=

тогда /р = —(ЕаАв/Ах). Однако подстановкой

уравнения (1.74) в выражения для А-х и Л„ в уравнениях (1.77)

можно легко показать, что /р =

1. Следовательно, полным реше­

нием является 1

+ / с, где дано соотношением (1.79). Константы

от Схдо С4 определяются из условия, что / = /' = () при х =

± а .

Тогда

^

__ р соз Н р а

д а - { - д Н р а соз д а

^

 

 

 

 

1

р соз д а

зш д а +

Ц з т Н р а со5 Н ра 9

\

' /

 

 

С2 — Сз — 0*,

(1.82)

 

 

р Нра соз да д соз Нра з т да

(1.83)

 

 

р соз да з т да +

д з т Нра соз Нра

 

 

 

 

и напряжения

 

 

 

 

 

 

ах — ах>т (1 •+- Сг соз крхсоз дх + С4 з т к рх зШ дх);

 

 

 

у

 

 

 

 

 

 

^ = ( 1о1 а*°^У йУ + В1У+ В») {[Сх (Рг— <72) + 2С4рд] X

 

 

X созкрх соз дх + [С4 (р2 — д2) — 2С&д] з!п крх з!п дх};

(1.85)

 

у

 

 

 

 

 

 

хХу=

—( | Охойу + 5 Х)[(— Схд + с 4р)созкрх зт дх +

 

 

 

 

+ (СхР +

С,д) з1пкрх соз дх].

(1.86)

1.8.3.

Другие случаи применения. Метод Хелденфелса и Ро­

бертса применен

Сингером

и

другими исследователями

[1.321

для определения напряжений в плоской пластине при различных видах распределения температур в поперечном сечении. Сингер с сотрудниками распространил этот метод также на пластины с пе­ ременными по хорде толщиной и температурой, подобные крыльям самолета и лопаткам турбин. В принципе метод расчета анало­ гичен методу для плоской пластины, за исключением того, что в выражении для дополнительной энергии учитывается толщина. Дополнительные подробности приведены в работе [1.32].

1.9. МЕТОД ОРТОГОНАЛЬНЫХ ПОЛИНОМОВ

Одно из наиболее существенных ограничений метода Хелденфелса и Робертса заключается в том, что решение принимается в форме только произведения функций от X и от У Если необ­ ходимо включить большее количество членов, то возникает во­ прос, в какой форме следует брать § (у).

Метод, основанный на использовании ортогональных полиномов самоуравновешенных параметров, предложен Хорви, который при­ менил его при решении задач о температурных напряжениях в пря­ моугольных полосах с различным распределением температур. Описание этого метода приведено в работах [1.16, 1.19, 1.2, 1.17].

1.10.ОБОБЩЕННЫЙ ПОДХОД С ИСПОЛЬЗОВАНИЕМ ФУНКЦИЙ ПОЛОЖЕНИЯ

Оба рассмотренных вариационных метода имеют серьезные ограничения. Метод Хелденфелса и Робертса по существу огра­ ничен выбором функции напряжений, состоящей из произведения функций от х и от у, и тем фактом, что если температура изменяется в направлении х и у, вариация должна быть выражена как произ­ ведение функций от х и от у. Метод Хорви ограничен случаями, когда температура изменяется только в одном направлении и ши­ рина в несколько раз больше высоты. Так как встречаются случаи, когда ограничения обоих методов являются чрезмерно жесткими, необходим подход, который не связан с этими ограничениями. Единственный недостаток описываемого ниже метода заключается в том, что он связан со значительным объемом вычислений, за­ висящим от требуемой точности. Поскольку, однако, большинству организаций, решающих подобные задачи, доступны быстродей­ ствующие вычислительные машины, это ограничение несуще­ ственно.

Основы этого подхода подробно описаны в разделе 1.6, где функции положения применительно к методу коллокации опре­ делены в более общем виде. Здесь можно указать только, что используются функции напряжений в форме уравнения (1.66), где коэффициенты Р (у), Р 2 (у). . имеют форму полинома. Подстав­ ляя выражение (1.66) для ср в уравнение (1.65) для дополнительной энергии и затем применяя условия минимизации (1.70), получают систему совместных дифференциальных уравнений для срх (х), Ф2 (х) ит. д., при этом может быть использовано любое количество функций ф„ (х), хотя для многих задач достаточно двух или трех функций. Уравнения могут быть проинтегрированы классиче­ скими методами, однако может также быть использован метод коллокации, значительно облегчающий труд. В этом случае диф­ ференциальные уравнения удовлетворяются только при некоторых значениях х, а не при всех его значениях (см. раздел 1.16).

На рис. 1.14 представлены результаты некоторых расчетов пла­

стины с соотношением размеров 6 x 2

и распределением темпера­

тур

в

поперечном

направлении

Т =

Т0 (1

у)9, проведенных

тремя

методами. Расчеты мето­

 

 

 

 

 

 

 

дами

Хорви,

Хелденфелса и

 

 

 

 

 

 

 

Робертса

взяты

из

 

работы

 

 

 

 

 

 

 

[1.321. Для энергетического ме­

 

 

 

 

 

 

 

тода,

использующего

функции

 

 

 

 

 

 

 

положения,

эти

функции

вы­

 

 

 

 

 

 

 

браны для значений у, равных

 

 

 

 

 

 

 

*/4 и 3/4, и получена система

 

 

 

 

 

 

 

обыкновенных дифференциаль­

 

 

 

 

 

 

 

ных уравнений

для

функций

 

 

 

 

 

 

 

напряжений

подстановкой

в

 

 

 

 

 

 

 

уравнение для энергии (1.65) и

 

 

 

 

 

 

 

на

основе

применения

вариа­

 

 

 

 

 

 

 

ционных уравнений (1.70). Эти

 

 

 

 

 

 

 

дифференциальные

уравнения

 

 

 

 

 

 

 

были решены двумя способами:

 

 

 

 

 

 

 

1) точным, который дал распре­

 

 

 

 

 

 

 

деление напряжений

ах (при

 

 

 

 

 

 

 

у =

—0,9

очень

близкое

к

 

 

 

 

 

 

 

распределению,

 

полученному

 

 

 

 

 

 

 

методом

Хорви

и Хелденфел­

 

 

 

 

 

 

 

с а — Робертса

(см. рис.

1.14);

 

 

 

 

 

 

 

2) методом

коллокации,

когда

 

 

 

 

 

 

 

эти

дифференциальные уравне­

Рис. 1.14. Распределение

осевых

на­

ния

удовлетворяются

 

только

в

отдельных

точках.

Метод

пряжений в пластине размерами 2 X 6

с постоянной температурой в осевом

двойной

коллокации

 

описан

направлении

и переменной — в попе­

в разделе

1.16.6;

цель

его —

речном Г =

Г# (1 — уУ при Е а Т 0 =

обеспечить решение таких урав­

 

92,5 кге/мм-

п у =

— 0,9

 

0 — 2 X 3 ,

коллокацня;

А — данные

нений.

 

 

 

 

 

 

 

 

рис.

1.14

видно,

что

Хелденфелса

и Робертса;

□ — данные

На

Хорвея;

■ — 2 X 4,

коллокацня;

V —

всего для шести точек (две в на­

3 х

4,

коллокацня;

■ н V

энергети­

ческий

метод с использованием функции

правлении у при трех в направ­

 

 

 

положения

 

 

лении

х)

полученные

значения

 

 

 

 

 

 

 

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

Хотя для этой частной задачи хорошие результаты получены всеми методами расчета, можно полагать, что в ряде случаев различие окажется существенным. Например, если температура изменяется в обоих направлениях, то метод Хорви применим, а метод Хелденфелса — Робертса может быть использован только

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

1.12ПРИМЕНЕНИЕ ТЕОРЕМЫ ВЗАИМНОСТИ

Д.Гудьер [1.91 применил теорему взаимности Бетти к выводу некоторых формул в задачах термоупругости. Этим способом в ряде случаев можно непосредственно получить конечный результат, который было бы сложно получить при детальном анализе напря­ жений и деформаций.

Теорема взаимности дает соотношение между перемещениями и силами двух независимых, действующих на тело систем внешних нагрузок, и формулируется следующим образом:

«Работа, произведенная силами первой системы на соответ­ ствующих перемещениях второй системы, равна работе, произ­ веденной силами второй системы на перемещениях первой».

Одна из выбранных систем сил обычно предельно проста, так что перемещения легко могут быть определены из предваритель­ ного анализа. В этом случае из теоремы вытекает соотношение между силами и перемещениями второй системы. Хотя это соот­ ношение является по существу следствием обычных уравнений теории упругости и граничных условий, оно часто получается в фор­ ме, позволяющей получить результат в конечном виде.

Возможность рассмотрения системы температурных напряже­ ний как эквивалентной системы объемных и поверхностных сил уже отмечалась. Теорема о взаимности перемещений может быть использована при рассмотрении этой системы сил в сочетании с не­ которой простой системой, для которой перемещения известны. Эффективность решения зависит от изобретательности исследова­ теля в выборе второй системы.

Гудьер получил некоторые результаты, представляющие не только академический интерес, но имеющие, практическое зна­

чение. Например, он показал, что изменение объема тела

 

Д Г = | З а 7 ^ .

(1.87)

Этот интеграл представляет собой изменение объема тела только за счет теплового расширения без учета температурных напряже­ ний, которые не оказывают влияния на изменение объема. Этот результат качественно очевиден, так как система температурных напряжений — это полностью внутренне уравновешенная система и в условиях упругости можно ожидать, что изменения объема от температурных напряжений компенсируются.

Этим методом можно также показать, что удлинение стержня постоянного сечения

М = -~ ^ аТ й У ,

(1.88)

где А — площадь поперечного сечения (сплошного стержня).

Изменение угла между двумя сечениями стержня (поворот при изгибе)

 

ш =

- Ы < *ТхЛ',

(1.89)

где

— момент инерции сечения относительно главной централь­

ной оси Оу; х — расстояние от центра в направлении

прогиба.

 

Прогиб консольной балки относительно заделки

 

 

б = — ±

| аТх (I — 8) йУ,

(1.90)

где

— момент инерции относительно главной центральной оси;

х — расстояние от центра в направлении прогиба; I — длина стержня; з — расстояние от заделки.

1.13.РЕШЕНИЕ МЕТОДОМ КОНЕЧНЫХ РАЗНОСТЕЙ

Дифференциальные уравнения равновесия и совместности обеспечивают соотношения между напряжениями и деформациями в точках, бесконечно близких друг к другу. Если эти уравнения разрешимы, то решение, естественно, дает точные значения ве­ личин. В большинстве случаев, однако, невозможно получить точное решение, поэтому необходимо прибегать к приближенным методам. Одним из таких подходов, который часто оказывается весьма плодотворным, является метод конечных разностей. Диф­ ференциальные уравнения заменяют системой таких конечно-раз­ ностных уравнений, которые связывают напряжения, деформации или другие величины, например, функции напряжений в точках, удаленных друг от друга на конечные расстояния. Эти уравнения, хотя и не точные, дают решения, которые близки к точным и удов­ летворяют инженерным целям. Решение таких уравнений обычно состоит из ряда последовательных расчетов, в которых значения определяемых функций получаются из известных значений функ­ ций на соседних участках. Обычно этот процесс дает численное решение задачи, так как решения в общем виде приводят к гро­ моздким формулам.

Уравнения в конечных разностях могут быть получены одним из двух способов:

1) на основе той же методики, которая используется при вы­ воде дифференциальных уравнений, уравнение в конечных раз­ ностях выводится относительно напряжений или деформаций при конечных расстояниях между точками. При составлении диф­ ференциального уравнения вначале рассматривают точки, удален­ ные друг от друга на конечные расстояния, а затем расстояния между точками устремляются к нулю. При составлении разност­ ных уравнений этот этап опускается;

2) дифференциальное уравнение формально превращают в раз­ ностное на основе известных определений для обычной и частной производных.

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

1.14. РЕШЕНИЕ ОБЫКНОВЕННЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ В КОНЕЧНЫХ РАЗНОСТЯХ

1.14.1.Пример. Вращающийся диск с радиальным температур­ ным перепадом. Метод конечных разностей иногда очень полезен при решении обыкновенных дифференциальных уравнений, осо­ бенно когда коэффициентами их являются сложные функции не­ зависимых переменных. Такое положение возникает в случае вращающегося диска с температурным градиентом. Если толщина диска изменяется по радиусу и свойства материала существенно меняются с изменением температуры, окончательные дифферен­ циальные уравнения имеют переменные коэффициенты и это де­ лает решение в замкнутой форме невозможным. Однако эта труд­ ность может быть легко преодолена при использовании метода конечных разностей.

Расчетную область делят на дискретные участки, и дифферен­ циальное уравнение превращается в разностное. Так, если напря­ жения на одном участке известны, напряжения на других уча­ стках могут быть подсчитаны с помощью последовательного при­ менения уравнений в конечных разностях. Когда напряжения на рассматриваемом участке не известны, напряжения на всех осталь­ ных участках могут быть выражены через значения неизвестных напряжений, которые затем определяются из заданных граничных условий.

Задача рассмотрена в работе [1.20]. В этом случае вывод уравнений проведен по методу, описанному в разделе 1.3.1. Уравнение совместности, превращенное в уравнение для напря­ жений на основе закона Гука, вместе с уравнением равновесия образует систему для определения радиального и кольцевого на­ пряжений в каждой точке радиуса.Уравнение равновесия имеет вид

(гНог) — /ш0 + р(д2г2Н=

О,

 

(1.91)

а видоизмененное уравнение совместности

 

 

)=°,

йг

Ег

сг0)

(1 +|х) (оу

 

(1.92)