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

книги из ГПНТБ / Талыпов Г.Б. Сварочные деформации и напряжения

.pdf
Скачиваний:
18
Добавлен:
25.10.2023
Размер:
12.43 Mб
Скачать

и нужно рассматривать не уравнения равновесия, а уравнения

движения:

+ д*худу

+

дххг

 

д2и

.

'

дахх

 

дх

дг

= Р d t 2

'

 

дххУ

dGyy

 

дХуг

= Р

d2v

.

(3.9)

дх

ду

 

дг

dt2

'

 

&*XZдх

+ dtyzду

+ дадггг

= Р ddtw2

'

 

 

 

 

 

 

2

 

 

где p —масса единицы объема тела. Но Дюгамелем [140] было показано, что изменение температуры во времени во многих слу­ чаях происходит с достаточно малой скоростью и влиянием инер­ ционных членов можно пренебречь, рассматривая движение как последовательность состояний равновесия (гипотеза Дюгамеля). Новые исследования [33, 34, 138] показали, что влияние инер­ ционных членов оказывается существенным только в массивных телах [33, 34], а в других случаях незначительно. Поэтому будем пренебрегать влиянием ускорений. Тогда при отсутствии объемных сил уравнения равновесия будут иметь вид:

дахх

і

дхху

+

дххг

=

0;

 

дх f

 

ду

дг

 

дт:ху

I

дауу

і

dxyz

=

0;

(3.10)

дх

1

ду

1

дг

 

дххг

1

dXyZ

+

dOzz

=

0.

 

дх

1

ду

дг

 

При упругих деформациях справедливы соотношения (3.7), которые дадут:

хх

 

2G

д2и .

ц

/ д2и _ j _

d2v

 

d2w

до

 

дх2

 

 

 

 

дх ду

+ дх дг

дх

 

 

1

- 2 ц

\дх2

 

 

 

 

а (1+1*)

д

(Т- Т0)};

 

 

 

 

 

 

 

1

дх

 

 

 

 

 

 

 

 

дгху

_ п

( д ? ± ,

d2v \ .

 

 

 

 

 

 

ду

~

\ду2

"Т" дхду

) '

 

 

 

 

 

 

д^хг

_ г

( д 2

и

.

d2w

\

 

 

 

 

 

 

дг

 

\ дг2

*

дхдг

)'

 

 

Подставив последние в первое из

уравнений

(3.10), получим

д2и

, д2и . д2и

. д2и

.

d2v

,

d2w

+

дх

2

"т" ду

2

дг

"г*

дх

 

дхду

 

дх дг

 

 

 

2

1

22

 

 

 

 

 

• 2ц / д2и . дЧ

| d2w \

 

(1 + ц)

д

 

1 — 2ц \ дх*

 

дхду

+

дх дг

 

1

дх

(Т — Т0) = 0.

Аналогично можно получить еще два уравнения. Используя опе­ ратор Лапласа

дх2

дг* ( ) = Д ( ),

 

9

2

эти уравнения можно написать в виде:

Дн •

1

де

2а(1 + ц)

д

(Т-

Го) =

0;

 

 

I — 2ц

дх

1 — 2ц

 

 

 

 

 

 

 

Av-

1

де

2 « ( 1 + Ц)

a

{Т-

0

)

=

0;

(З.П)

1 — 2ц

д</

1 — 2ц

ду

Т

 

Aw-

1

да

2а(1 + ц)

д

(Т-

Т0)

=

0.

 

1 — 2ц

аг

1 — 2ц

дг

J

Последние впервые были получены почти одновременно Ней­ маном и Дюгамелем и называются уравнениями Дугамеля— Неймана. Они отличаются от обычных уравнений теории упру­ гости, например от уравнений Ляме, наличием членов, пропорцио­ нальных градиентам температуры. Таким образом, учет влияния неравномерного нагрева сводится к учету дополнительных массо­ вых сил, пропорциональных градиентам температуры.

 

 

12. ГРАНИЧНЫЕ

УСЛОВИЯ

При

наличии

внешних напряжений в точках поверхности

тела с

составляющими Xv, Yv, Z v

на этой' поверхности должны

быть выполнены

условия:

 

 

 

 

(3.12)

 

 

xxzl -j- xyzm -(- <jzzti = Zv,

где I, m, n — направляющие косинусы нормали v к поверхности тела в данной ее точке. Те же уравнения, если использовать (3.7), запишутся в виде:

( T V + *E-)'+(-S-+-S> +

( ^

ди+ i r

 

2 а ( 1 +Ц) (Т-Т0)1

 

 

 

1 — 2 ц

 

 

 

 

dv

+

 

) ' + ( b V +2

^ ) " I

(3.13)

dw

 

 

 

(~дг + дх

2а(1 + Ц) ( Т

 

 

 

. / 2ц

, п

dw \

rp\n

\ Z v

 

Сравнивая уравнения (3.11) и (3.13) с соответствующими уравне­ ниями теории упругости, видим, что температурная задача теории

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

2 а ( 1 + ц )

д

( Г - Г 0 ) ;

-

2ct (1+1-0

д

(Т-Т0);

1 — 2 ц

дх

 

 

1 — 2ц

ду

 

 

 

2 а ( 1 + ц ) д

(Т-Т0)

 

 

 

 

1

дг

 

 

 

и дополнительное поверхностное давление с интенсивностью

2а(1 +ц) , г

_ _ .

 

) .

1

(7- - Г

0

 

 

 

 

Компоненты дополнительной массовой силы, обусловленной неравномерностью нагрева, могут быть исключены из уравнений (3.11), если известно частное решение этих уравнений. Предста­ вим компоненты перемещения в виде:

и — иг + и 2, v = v1 + v2; w = wx + w2

(а)

и предположим, что имеют место равенства:

и, = дх '

dF .

 

6F

(б)

ду '

2

dz

Если подставим (а) и (б) в (3.11) и выберем функцию F так, чтобы она удовлетворяла уравнению Пуассона

A F

= a l ± i i ( T -

Т0),

(3.14)

 

 

1 - ц

 

 

 

 

то вместо системы (3.11) получим:

 

 

 

 

A " i

+

1

 

0;

 

 

1 — 2ц

 

4

.15)

Д УХ

+

1

дег

0;

(3

1

д#

Д o»i

+

1

аег

0.

 

1

дг

 

 

где

 

 

 

dwt

 

 

 

 

 

 

 

 

 

 

ду

дг

 

 

Граничные условия (3.13), если иметь в виду (3.14) и учесть, что

dF

. .

dF

dF

dF

дх

1

ду

dz

dv '

где v — нормаль к поверхности, примут вид:

(Т^.+^)'+(ТГ+ТЙГ)"1 +

\ду 1 дх

J -

1 V

2 , -

- - д у

) т +

і da>i ) п = 2 а т { ± і і ( Г - Г 0

) _ 2 - | 7 ^

 

 

 

1 - ц

 

a# av

5«x .

дал

 

 

 

(3.16)

 

V dz

 

 

 

Эх

і * 1

1

 

2ц.

 

 

1 - І * ( Г - Т 0 ) - 2 az av

 

 

Уравнения (3.15), (3.16) являются обычными в теории упругости при наличии поверхностных сил с компонентами:

9

1 + ц

- al(T —T0)-

2

a

af .

 

1 — Ц

 

 

ax

av

'

0

1 + ц

-am (T —•T0)-

2

a

af .

&

1 — ц

dy

av

'

 

 

 

о 1 + и

an (T —T0)-

2-

a dF

 

 

1 —(і

 

 

аг

dv

 

Функция F определяется по теореме Пуассона для объемного потенциала как решение уравнения (3.14)

,(3.17)

V(x-t)*

+ (y-r))* + (z-Q*

где

О - т 0 ,

г (6,'т), о = г (g, л,

интегрирование проводится по всему объему со тела и dco = db, dr\ dt,.

13. УРАВНЕНИЕ БЕЛЬТРАМИ-МИТЧЕЛЯ

Кроме уравнений равновесия и условий на поверхности должньї быть выполнены уравнения совместности деформаций (3.5). Последние, выписанные через компоненты напряжения, назы­ ваются уравнениями Бельтрами—Митчеля. Выведем эти уравне­ ния с учетом температурных членов в предположении, что дефор­ мации остаются упругими. Для этого продифференцируем первое из уравнений (3.11) по х:

л

д и Л.

 

1

 

 

д*е

2а(1 + ц.)

д* ,Т

т л _ о

/о і о\

А

— +

і _

9

 

д т

1 - 2 ц

~д^~^

°'

>

 

дх

1

дх%

 

 

 

 

и выразим все входящие сюда члены через компоненты напряже­ ния. Из соотношений (3.7) имеем:

 

 

1 + 1*

Зо(1 + ц )

(Т-Т0)

=

УУ

2G 1

1

 

20(1 +Ц)

[ в - З а

( Г - Г о ) ] ,

 

(3.19)

 

1 — 2 ц

 

 

 

 

откуда

(3.20)

Продифференцируем уравнения (3.11) по х, у, z и сложим. Это даст

Ae = ^ ± f

А(Т-Т0).

(3.21)

1 [X

Из (3.19) имеем

А о" = 2 G , ( 1 + Г }

[ А е -

 

1 — 2ц

 

или, учитывая

(3.21), получим

ACT _

2G(l + Ю Г a O + j x )

 

1 — 2 ц

1 - ц

откуда

За А ( 7 - 7„)]

За А ( Г - Г 0 ) ,

А о = - 4 в < ; + ц Р » ° А ( Г - 7 - 0 ) .

(3.22)

Далее, имея в виду (3.19), из первого соотношения (3.7) получим

да

схх

JMJ

+

а(Т-Т0).

(3.23)

~дх~

2G"

2(1 + ц ) 0

Подставив в (3.18) выражения (3.20), (3.23) с учетом (3.22), полу­ чим первое уравнение Бельтрами—Митчеля. Аналогично получим еще два уравнения, которые в совокупности можно представить в виде:

д а

л

1

д*о

2ctG(l

+(х)

, т

г , .

А с*, +

утр^-^г +

А (V -

Г0 ) +

Л

і

+ 2 a G - ^ ( 7 ' - 7 '

0 ) = 0;

т . .

1

d % .

2aG(l +

ц)

г

 

 

+ 2 а О ^ г - ( Г - Г 0

) = 0;

(3.24)

 

 

 

A a < z

+

1

д2а

 

 

 

 

 

1 + (А йг2 + ^ t " * А ( Г - Г 0 ) +

 

 

+ 2 а О - ^ ( Г - 7 0

) = 0.

 

Получим остальные три уравнения Бельтрами. Для этого сначала продифференцируем первое из уравнений (3.11) по у, а второе — по л; и сложим. При этом получим

2

д2е

4 а ( 1 + ц )

(Т-Т0)

= 0. .

1 — 2ц

дх ду

1 — 2щ дх ду

Аналогично получим еще два уравнения. Имея в виду (3.20) и последние три из соотношений (3.7), эти уравнения приведем к виду:

 

1

 

1

д2а

+

2aG дхду

(Т-Т0)

= 0;

 

 

1

1

+

дхду

 

А т , г

1

І

1

дхдг

 

 

а2

( Т - Г 0

) _ 0 ;

(3.25)

 

1

+

 

 

 

 

 

А т.

+ . 1

дудг

+

2aG

а2

( 7 - 7 0 ) _ 0 .

 

У*

'

l+ii

 

 

дудг

 

 

 

Кроме трех уравнений равновесия и условий на поверхности компоненты напряжения должны удовлетворять шести уравне­ ниям Бельтрами (3.24) и (3.25).

14. ПОТЕНЦИАЛ ТЕРМОУПРУГИХ ПЕРЕМЕЩЕНИЙ

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

 

dF

dF .

 

dF

(3.26)

и =

дх

V = • ду '

W =

дг

При этом, имея в виду, что

A F; е= AF,

систему (3.11) можно привести к одному уравнению (3.14):

(3.27)

т. е. искомая функция F должна удовлетворять уравнению Пуас­ сона и называется потенциалом термоупругих перемещений. Если функция F найдена как решение уравнения Пуассона, то компоненты перемещения определяются по формулам (3.26), а для деформации и напряжений получим:

 

дгр

 

 

 

дгр

Qzp

 

Є х х

~ дх2

 

еуу~

 

й/2

»

 

 

 

d*F

 

_

0

d*F

' 7^ = 2

d2F

 

УхУ~2'дхду

'

У

х г ~

1

dx dz

дудг

 

 

— 2G (

d2F

. JPF_

— 2G

d2F

(3.28)

 

ду2

*~ dz2

д х д у

 

 

 

 

 

 

 

d*F

 

 

 

 

 

 

 

 

dx dz

 

 

 

 

 

 

 

d*F

 

^ = - 2 G ( ^

 

+ | ^ ) '

^ 2 G d y d z -

 

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

15. МЕТОДЫ РЕШЕНИЯ УРАВНЕНИЙ ТЕРМОУПРУГОСТИ

Общий интеграл уравнений Дюгамеля—Неймана

Решение обратной задачи термоупругости не вызывает особых трудностей: когда заданы напряжения или деформации, она сво­ дится лишь к интегрированию уравнений Коши. Решение прямой задачи термоупругости приводится к интегрированию сложной системы дифференциальных уравнений в напряжениях или в пере­ мещениях. В первом из двух последних случаев необходимо найти шесть неизвестных функций ахх, . . ., ххг из трех дифференциаль­ ных уравнений равновесия (3.10) и шести дифференциальных уравнений Бельтрами (3.24), (3.25), удовлетворяющих простым алгебраическим условиям на поверхности (3.12). Во втором из этих случаев необходимо найти три неизвестные функции и, v, w решением трех дифференциальных уравнений Дюгамеля— Неймана (3.11), удовлетворяющих граничным условиям в пере­ мещениях (3.15). Общий интеграл этих трех уравнений можно

выписать [91 j в виде: и — «о +

« ! + « 2 ;

 

v =

v0 +

vt + vi;

 

(3.29)

w = w0-\-w1-\-wz,

,

 

где при наличии объемных сил с составляющими X,

Y и Z:

 

 

1

р)

dW0

 

 

2

( 1 -

дх

 

А Т 2

 

1

дУ0

(3.30)

2

( 1 -

V)

ду

 

w0 = Л

 

1

 

дУ«

 

- 2 ( 1 - -ц)

дг

 

здесь Wly Ws — любые непрерывные, включительно до своих третьих производных, решения уравнений:

Д Д У ,

=

- • G

(3.31)

Д Д Ч г з =

-

Z

G

 

дх

+ Elду

+ дУдг3

 

 

 

 

4(1—ц)

dx

(жФі + ^

+

гФз +

Фо);

 

 

 

 

 

1

- ^ ( х Ф і + ^ +

гФз +

Фо);

(3.32)

 

° і - ф « — 4 ( 1 - | і )

ф

 

 

 

 

 

причем

Ф 1 ;

Ф 2 ,

Ф 3 — три независимых

общих интеграла

урав­

нения

Лапласа

 

 

АФ, =

0;

 

 

 

(3.33)

 

 

 

 

 

 

 

 

 

 

 

 

 

1 —(X

дф

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дф

 

 

 

(3.34)

 

 

 

 

 

1 —ц

ду

 

 

 

 

 

 

 

W9 =

1+(Х

дф

 

 

 

 

 

 

 

 

1 —Ц дг

 

 

 

 

а Ф — частное

непрерывное,

включительно

до

своих вторых

производных, решение

уравнения

Пуассона

 

 

 

 

 

 

 

АФ = аТ.

 

 

 

(3.35)

Можно

показать

[91], что при ц = 0,25

функция Ф 0 в выраже­

ниях (3.32)

является

лишней.

Найденные

перемещения

(3.29)

должны удовлетворять граничным условиям (3.13). Эти условия, связывающие функции Ф, с заданными напряжениями на поверх­ ности тела, значительно сложнее алгебраических равенств (3.12). Поэтому в некоторых случаях интегрирование уравнений (3.10), (3.24), (3.25) может оказаться проще, чем разыскание трех функ­ ций Ф,-, удовлетворяющих условиям (3.13). Нахождение трех гармонических функций Ф£ , удовлетворяющих уравнениям (3.13), остается основной трудностью на пути получения полного решения общей прямой задачи теории упругости. Учет влияния объемных сил требует нахождения частных решений уравнений (3.31), а учет влияния температуры — нахождения частного решения урав­ нения (3.35).

Полуобратный метод Сен-Венана

При этом методе часть решения рассматриваемой задачи задается (угадывается), а другая часть определяется так, чтобы была удовлетворена система основных уравнений термоупругости. Та часть решения, которая задается (угадывается), должна соот­ ветствовать рассматриваемой задаче и не должна быть в проти­ воречии с основными уравнениями термоупругости [8, 91 ] .

Потенциал термоупругих перемещений

Для получения частных решений статических, квазистатиче­ ских и динамических задач термоупругости можно использовать введенный в п. 15 потенциал термоупругих перемещений. Един-

ственным ограничением, наложенным на этот потенциал F, яв­ ляется то, что эта функция F — частное решение уравнения Пуассона (3.27). Поэтому найденные через эту функцию F напря­ жения не обязательно должны удовлетворять заданным условиям на поверхности тела. Если полученное решение не удовлетворяет заданным условиям на поверхности, то приходится решать до­ полнительную задачу (см. ниже п. 16, 17). Изложение и приме­ нение метода потенциала термоупругих перемещений можно найти в работах [67, 92].

Энергетические методы

Для решения практических задач термоупругости часто исполь­ зуют следующие энергетические методы.

Принцип Даламбера. В соответствии с этим принципом работа, совершенная внешними силами и силами инерции на возможном перемещении тела из некоторого мгновенного состояния, равна изменению энергии деформации, т. е.

bA-\P(^&u

 

 

+

^8v

+ ^-6w)d<*

= W,

(3.36)

где

 

 

 

 

 

 

 

 

t/ = G J [

^ +

^

+

4 +

T ( v L + ^

+ vL) +

 

+

Т ^

е 2

-

 

1 Г = І

2 Г а ^ ] ^

 

(3.37)

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

= Ш.

(3.38)

Принцип Гамильтона. В

соответствии

с этим принципом

в случае, когда внешние силы имеют потенциал, вариация инте­ грала по времени от разности суммарного потенциала внутрен­

них и внешних сил и кинетической энергии равна

нулю, т. е.

t

 

b\(n — K)dt = Q.

(3.39)

о

 

В случае, когда имеем статические и - квазистатические задачи этот принцип переходит в принцип минимума потенциальной энергии

ЬП = 0.

(3.40)

Принцип виртуальных изменений напряженного состояния.

В соответствии с этим принципом при всяких виртуальных изме-

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

2 (и8Х + v8Y + w8Z) = 8U,

(3.41)

где X, Y, Z — означают составляющие объемных

и поверхност­

ных сил. Из этого принципа вытекают теорема Кастильяно и теорема о наименьшей работе.

Начало взаимности. По этому началу, если рассматриваются два состояния упругого тела, работа сил первого состояния на соответствующих им перемещениях второго состояния равна работе сил второго состояния на соответствующих им перемеще­ ниях первого состояния. Методы применения этого начала к за­ дачам термоупругости разработаны в работе [62].

Эти основные энергетические методы широко применяются для решения изотермических и неизотермических задач строи­ тельной механики (стержни, стержневые системы, пластины и т. д.). Изложение этих методов и их применение можно найти в работах [8, 26, 62, 67, 91, 92].

Численные методы

Из численных методов наиболее приспособленным к машин­ ному счету является метод конечных разностей [47].

Методы плоской задачи

Наиболее эффективными методами решения плоской задачи термоупругости являются метод функции Эри и метод комплекс­ ных переменных. Изложение и применение этих методов можно найти в работах [8, 68, 92, 119].

16.ПЛОСКАЯ ДЕФОРМАЦИЯ

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

мерным распределением температуры Т (х, у) = Т (х, у) — Т0. Если ось oz совместим с геометрической осью этого тела, то в этом случае получим, что

и= и(х, у);

v =

v{x,

у);

 

dw

_

dw

Q_

(3.42)

дх

~

ду

'

 

dw

=

const = с;

 

~дг

 

 

 

 

4 Г. Б. Талыпов

49

Соседние файлы в папке книги из ГПНТБ