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

книги / Физические основы прогнозирования долговечности конструкционных материалов

..pdf
Скачиваний:
1
Добавлен:
19.11.2023
Размер:
29.97 Mб
Скачать

что соответствует

закону

Гука е = а/Е,

где а = Го/а2— напря­

жение; E = f/a — модуль растяжения (модуль Юнга).

 

Отсюда видно, что закон Гука возникает в гармоническом

приближении межатомного потенциала

(это следует из того,

что в

выражении

(1.11а)

отсутствует параметр

ангармонично­

сти g)

и справедлив только при малых

силах (Г0

Fm). При

больших силах

возникает отклонение от закона Гука: связь

е(Го)

становится

нелинейной. Наконец, при Го ^

Fm деформа­

ция £

(1.11) и параметр f

(1.10,6), ответственный

за колеба-

Рис. 1.7. Деформация

Рис.

1.8.

Зависимость

е межатомной

связи

силы

Fy

притяжения

в зависимости

or рас­

от удлинения q меж­

тягивающей

силы

атомной связи

ния, становятся мнимыми, т. е. при Fo=Fm колебательный ре­ жим разрушается. В этой связи величина

Fт

1L

( 1. 12)

 

2g

 

является разрушающей силой. Ее физический смысл заключа­ ется в следующем. Действующая при колебаниях «возвращаю­ щая» сила

_

dtp

о

Гф —

dq

Чч + ■fq-

«имеет» график, приведенный на рис. 1.8, откуда видно, что при q ^ f / g имеется максимум силы, равный (1.10).

Таким образом, разрушение наступает, когда приложенная сила больше максимального значения возвращающей силы, осу­ ществляющей колебательный режим. Заметим, что при / = 0 потенциал приобретает вид, приведенный на рис. 1.6, линия 3: исчезает потенциальная яма, в которой могут происходить ко­ лебания. Иначе можно сказать, что при разрушении исчезает «связанное» состояние, характеризуемое наличием потенциаль­ ной ямы, осуществляющей колебательный характер движения системы.

Подчеркнем, что разрушить можно только ангармонический осциллятор: в отсутствие энгармонизма при g = 0 уравнение

31

движения осциллятора (уже гармонического) под действием внешней силы следующее:

тс/ + fq = F0

т. е. q = (Fo/f) +А cos (щ/ — Ф). Движение носит характер гар­ монических колебаний вокруг нового положения равновесия растянутой по закону Гука пружины.

Вернемся к выражению (1.106). Подстановка в него (1.11) дает

J * = fil - 6 U t),

(1.13)

где введен так называемый коэффициент Грюнайзена

G =

(U 3a)

Таким образом, учитывая (1.106), можно сказать, что под действием деформации растяжения в межатомного расстояния а (е=А а/а) частота колебаний становится равной

6 = д / - ^ - = со V l — 6Ge.

(1.14)

При малых e«Cl/6G выражение (1.14) принимает вид

— — = 3 Ое, А«о = со — 6 .

(1.14а)

Результат (1.14а), называемый соотношением Грюнайзена, устанавливает относительное уменьшение частоты —Да>/<о при деформации е. Как видно, эффект уменьшения частоты обуслов­ лен энгармонизмом межатомного потенциала и исчезает при G =0. С учетом закона Гука е = а/£ соотношение Грюнайзена принимает вид

Лео

(1.146)

10

 

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

Наконец, разделив выражение (1.12) на а2, с учетом (1.13а) получим выражение для напряжения:

oth —

Е

0,1В,

(1.15)

аг

126

 

 

 

которое выдерживает, не разрушаясь, ангармонический осцил­ лятор.

При порядковой оценке учтено, что для твердых тел G= = 1ч-3. Выражение (1.15) дает хорошую оценку предельной прочности твердых тел (так называемого теоретического пре­ дела прочности).

32

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

Уравнения движения такой системы имеют вид

tniji = —

[ff (qt q .- ,) + <f (9« + I — 9Л]. * = 1,

2, .... N,

где qt — смещение

от положения

равновесия

i

массы

(атома) т.

 

 

 

 

 

 

После алгебраических преобразований имеем

 

 

Ф= ~ - ~ Г (2<7>— <7, - . — q>4- ,)[l

(<?■: +, —

<7i)J •

(1.16)

Уравнения

(1.16)

при g = 0 совпадают

с соответствующими

уравнениями движения своего гармонического аналога. В ан­ гармоническом приближении уравнения движения являются не­ линейными и не допускают точного аналитического решения. В зависимости от величины деформаций e=dq/dx возможно два рассматриваемых ниже типа поведения системы: нелинейные колебания (при малых в) и стохастичность (при больших с).

Длинноволновые колебания в стержне. Перейдем от дискрет­ ной последовательности масс кристаллической решетки к не­ прерывной (сплошной) среде (континууму), осуществив для этого предельный переход а~> 0 — расстояние между мас­ сами в отсутствие колебаний). Применительно к кристаллу, где межатомное расстояние а фиксировано, такой переход озна­ чает, что рассматриваются только длинные волны с длиной волны л а, или ак со аГк <С 1. В этом случае дисперсионное со­ отношение (1.9) принимает вид

01=2V^rls,n'f'h aV^r и*

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

 

v = — = а л / — ^=л / — 9

 

к

\

m

V

р

где р = т / а 3— плотность.

 

(1.16)

при переходе от / к ко­

Смена

уравнений движения

ординате

x = ia превращается в

единственное конечно-разност­

ное уравнение

 

 

 

 

 

q (*) — ---- ^

[2<7(х) — q(x — а) —

 

<7(* + a)]{l ~§^[q(x + a) <7(x)]}.

3 Заказ X® 248

33

При разложении

q

(х±а) по а -> 0 с удержанием трех

членов

 

 

 

 

 

<j(x±a) = q(x)±a-% jr +

^ - ^ ! r

выражение (1.16) дает

 

 

 

 

 

=

о2/

Л

ва

дд \

д2д

q ~

т

V

f

дх /

дх*

или (с учетом выражения для скорости v длинных волн и опре­ деления коэффициента Грюнайзена G (1.13а))

< м б а >

Выражение (1.16а) при G= 0 приводит к классическому вол­ новому уравнению в стержне, деформация которого подчиня­ ется закону Гука. Считая далее величину G малым параметром, в первом приближении, сохраняющем лишь первые степени G, т. е. представляя q в виде

q = q 0 + GGq\

имеем

q*= v2 д*д° . (Эх2 *

Я = гг

эу

(Эх2 V дх дх2 •

Таким образом, нулевое приближение определяется одно­ родным волновым уравнением, а возмущение ql— решением неоднородного волнового уравнения (с «источником»).

Изотропный континуум. Обобщим приведенные выше соот­ ношения для одномерной сплошной среды на случай трехмер­ ного тела, свойства которого одинаковы во всех направлениях. Его точка с радиусом-вектором г при деформации тела смеща­ ется на величину q(r). Квадрат расстояния между двумя близ­ кими точками

d r 12= dr2+ 2tik dxi dxj,y

где Xi — декартовы координаты, i, k= 1, 2, 3.

По всем дважды повторяющимся индексам производится

суммирование

 

(aik + а.„ + ацаш), щк

(1.17)

где е*А — тензор деформации.

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

34

случае в кубическом приближении (аналогичном (1.2)) с точ­ ностью до константы (являющейся потенциалом недсформиро­ ванного тела)

2 . / к

и Л 2 ,

А

,

-

л

С

2

q = ре,* + \^~2----- з-) Ы1+

“3-

+

B*ik£u + “

 

£//»

где ft и К — модули сжатия и сдвига; А , В, С — константы (все эти коэффициенты в рассматриваемом случае изотропного тела являются скалярами).

Подставляя сюда выражение (1.17) и оставляя члены до третьего порядка включительно, получаем

ф =

4~“ a l i ) 4 “ £-------- а ц

+

4 ---- ~ ^ CLikCLiiCLik 4“

, / Д + К

BN

2 . -А

 

, В

. С з

+ ^----- О----------- з" j

4- ”j~2a ika kla li 4---9 O'ik^kiG'll

4--- 3a iu

 

 

 

 

 

(1.18)

Уравнения движения среды следующие:

 

 

 

Mi

d*Hk

*

(1.19)

 

 

 

дщ

 

где р — плотность недеформированного тела; 0,7*— тензор на­ пряжений,

(1.19а)

При равновесии смещения не зависят от времени н опреде­ ляются из равенства нулю правой части выражения (1.19). Рели смещения малы настолько, что в выражениях (1.17) и (1.18) можно ограничиваться соответственно линейным и квад­ ратичным членами п тензор деформации принимает вид

г

1

[

dqi

| dqk \

lk

2

\

dxk

dxt ) ’

то тензор напряжений (1.19а) есть

Oik= 2U£//; -J-(к ----g- *'/ 6ik,

I,

i = k;

bill = { 0,

i =5^=k

(символ Кронекера).

Это выражение называется законом Гука для изотропного тела. Введя модуль Юнга Е и коэффициент Пуассона vr] соот­

ношениями

[г= Е/2 (1 4- vn); К = Е/3 (1 — 2vn)»

з *

35

из закона Гука для компонент тензора деформации имеем

E t x x =

о х х

( р у у + o 2Z) ;

Е&УУ=

(Jyy

V l l ipXX"f* ®zz)\

E&z — ozz — vn (oxx + am/);

E tXy = (1 -f- Vnj oXffi

Егху = (1 -f- vn) OJCZI

Егтуг — {\ -f- vn) a,,z.

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

а все другие компоненты равны нулю. В приближении закона Гука отличны от нуля только диагональные компоненты сць. (1.17):

 

 

-= Ео;

 

 

 

 

 

&уу — йщ--

 

 

 

Если

коэффициент Пуассона 0 < vn < 7г

оказывается

до­

статочно

близок

к нулю, то величины aTijtj

и аГ^

малы

но

сравнению с аг .

В этой связи, переходя к выражению для 0,л

(1.19а) при условии <р (1.18) в случае одноосного растяжения стержня, будем пренебрегать всеми компонентами, кроме аХх=

=- р . В этом приближении имеем

дх

ддх

охх Е дх +(х+к+4 +в+с)(^гУ-

Сравнивая это выражение (с учетом движения (1.19)) с (1.16а), находим, что коэффициент в скобках равен 3GE.

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

F = f idV,

V

где f — объемная плотность силы.

36

Этот интеграл состоит из двух слагаемых. Первое слагаемое представляет собой сумму внутренних сил, равную нулю по третьему закону Ньютона. Второе слагаемое является суммой сил, действующих через поверхность данного объема со сто­ роны его окружения. Для того чтобы интеграл по объему пре­ образовать в интеграл по поверхности, необходимо представ­ ление

с_ duik

''о х к

Тогда

Тензор Gih называется тензором напряжений. Из его опреде­ ления видно, что произведение cnt dSh есть /-я компонента силы, действующей на элемент поверхности dS. Соответственно ощ — /-я компонента силы, действующая на единицу поверхности, пер­ пендикулярную декартовой координатной оси xtl. Как и пол­ ная сила, действующая на объем К, се момент также должен выражаться в виде интеграла по поверхности объема. Это означает, что в выражении для компонент тензора момента силы

Mik = .((f,x, -

ft,Xi)dV = J

x, -

x )d V =

Г

V

 

 

-■ Ф (ацХк — oHXi) dS[ — J (dk — On) d v

(1.20)

5

V

 

последнее слагаемое равно нулю.

Другими словами, мы установили симметрию тензора напря­ жений. Нормальное напряжение на произвольно ориентирован­

ной площадке равно

з

 

 

 

=

2 в1}П'П;.

 

i. !=\

 

Касательная составляющая определяется по формуле

 

(Тюс п= V p t + Pz + Ps — On,

где

3

 

£ Oijtij — компоненты вектора напряжений на пло-

щадне с нормалью п. Сумма

 

 

®= 4-

2

 

°

it / —I

37

называется средним, или гидростатическим, давлением, а тензор Sij = а*/ — 067

называется девнатором напряжений; — символ Кроиекера. Введем понятие главных напряжений. В прямоугольной де­ картовой системе координат охи направление, определяемое еди­ ничным вектором п, называется главным направлением тензора

Оц, если вектор напряжений Р (его

компоненты

Р,*— £ cr/j/i/ —

это так называемая формула Коши)

параллелен

/^1

п. Это озна­

чает, что

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

Р,-= ал/, где

а — некоторая

скалярная

величина, подлежащая

определению

из системы

уравнений

 

 

 

 

 

 

п

 

 

 

 

 

Е (О,7 — оЬп) t i i — O,

 

i = 1, 2, 3.

 

 

Условием существования ненулевых решений является ра­ венство нулю определителя, составленного из коэффициентов этой системы:

| оц obtj | = 0.

Это уравнение является кубическим и может быть приве­ дено к виду

гГ — /,( Г

+ /л т —

/а =

о ,

 

где

 

 

 

 

 

 

 

 

 

Л =

ап +

 

+

°зз;

а,,

а2, I

О22

032

+

 

<Ь|

+

С о Л

033

О,

Озз

Н|2

022 I

 

 

 

ои

Оо,

(УЯ1

 

 

 

/з =

 

СУ10

022

032

 

 

 

 

 

013

023

033

 

 

По свойству корней кубического уравнения 0*, /= 1, 2, 3, со­ отношения между ними и коэффициентами / г имеют следующий вид

 

 

Л =

01 ~J~ 02 + 0:ь

I о О

01

о

 

 

=010.» + 010з + 0203;

+

О

03

+

 

 

 

 

0

0

 

 

 

0

о2

0

-- 010203'

 

 

0

0

о3

 

38

Корни oi называются главными компонентами тензора на­ пряжений. Они всегда вещественны н не зависят от выбора ко­ ординатной системы. Соответственно являются инвариантными и величины U.

Три главных напряжения вводят инвариантное (для опреде­ ленности декартово) о-простраиство. В нем любые комбинации компонент <Xj, 02, 0.3 остаются неизменными (инвариантными) при преобразовании системы координат в физическом простран­ стве. Так как тензор напряжений имеет три базисных инвари­ анта, в качестве которых могут быть выбраны три любых независимых инварианта этого тензора, то можно образовать, вообще говоря, неограниченное множество инвариантных про­ странств тензора напряжений. Выбор того или иного инвари­ антного пространства неоднозначен и определяется характером задачи. Например, в ряде случаев удобно рассматривать Н-про- странство [192]:

И, = а , +02 +

03;

Пг =

0! +

о- +

05;

11з =

0? +

02 +

03-

называются главными касательными напряжениями. Комбинацию

называют интенсивностью касательных напряжений.

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

Механические критерии разрушения изотропных материалов при сложном нагружении. При механическом подходе, не рас­ сматривающем природы разрушения и опирающемся на пред­ ставление о существовании предельного состояния (предела прочности) материала, тела считаются твердыми (неразрушен­ ными), если изображающие их напряженные состояния точки в о-пространстве (главных напряжений 0/) не выходят за гра­ ницы некоторой (называемой предельной) поверхности. Эта поверхность отвечает представлению о пределе прочности. Вы­ бор ее неоднозначен и связан с «теорией прочности» (см. вве-

39

деиис) и проверкой на опыте. Однако выбор ограничен тем* что предельная поверхность должна удовлетворять ряду усло­ вий. Во-первых, для изотропных материалов уравнение пре­ дельной поверхности в су-пространетве

F (o d = 0, / = 1 , 2 , 3 ,

должно иметь действительные корни и быть инвариантом лю­ бых перестановок переменных су*. Это выполняется для пере­

менных Л ь Л ^, Л я.

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

у M i > К i.kLi d°i<'r'k

где 6* — некоторые произвольные переменные. (Коль скоро по­ нятие прочности сохраняется в кинетической теории, в ней эти условия должны быть также выполнены.)

Не ставя перед собой цели проведения детального обзора ме­ ханических критериев разрушения (они содержатся в работах [113, 192, 235]), рассмотрим основные из них главным образом для контраста с развиваемыми в монографии кинетическими подходами.

1. Критерий Галилея—Ренкина наибольших нормальных на­ пряжений. Принимается, что разрушение наступает, когда наи­ большее главное напряжение ст/ достигает предельного значе­ ния аПр. Здесь предельной поверхностью является куб со сто­ роной 2а1ф.

2. Критерий Кулона наибольших касательных напряжений. Наступление разрушения связывается с критическими значе­ ниями максимального касательного напряжения:

шах (| о, — сь |, 2 — ст., I, |оа — а,|) — а11р.

Это условие в а-простраистве задаст поверхность правиль­ ной шестигранной призмы с осью 01=02 = 03, называемой гидро­ статической (точкам на этой оси отвечают гидростатические на­ пряженные состояния). При плоском напряженном состоянии, когда a/ = ai, 02, 0, предельное условие определяется сечением призмы плотностью су3 = 0, имеющим форму шестиугольника.

3. Критерий Максвелла наибольшей интенсивности касатель­ ных напряжений. Здесь предельная поверхность есть

(<У1 — ОI f (о? — Оз)~ “Г (оз — 0\) = 2оПр.

Она соответствует поверхности цилиндра с гидростатической осью, описанной вокруг призмы Кулона (предыдущий кри­

терий).

Стохастичность. Чтобы выяснить характер движения атомов в ангармоническом кристалле при больших деформациях, пре-

40