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

книги / Численные методы в теории упругости и пластичности

..pdf
Скачиваний:
5
Добавлен:
20.11.2023
Размер:
10.57 Mб
Скачать

где 5 € Их, т.е. некоторой площадке, на которой вычисляется оператор [■

Упражнение 3.7. Показать, что для вектора V (3.48) справед­

ливы соотношения

 

 

 

= 1 м ю

+ 1 Р ?\ п, уЫ уу,

(3.49)

 

«ГО») =

+ / р^ Ы уЫ у) ^ ,

(3.50)

 

 

I!

 

где величины

и о,~ определяются по формулам (3.41).

 

Решение второй краевой задачи (3.13), (3.14) теории упругости для внутренней области V будем искать в виде суммы двух решений (3.21), первое из которых «{• находится по формуле (3.3),

а второе «" — в виде потенциала простого слоя (3.4):

 

«{'(*) = / ^ ( М Ы ^ П , .

(3.51)

Е

 

Учитывая соотношения (3.50) для плотности р(^), получаем интег­ ральное уравнение Фредгольма второго рода

\п(п) + I Р } % , У)Рь(у)<Ку = 5?(Ф -

|Е.

(3.52)

Е

 

 

Упражнение 3.8. Показать, что для случая внешней области для задачи (3.13), (3.14) справедливы те же формулы (3.21), (3.3) , (3.51), где х $ V, а формулу (3.52) нужно заменить формулой

+ 1'Р$к)% у)рк(у)Жу = 5?0» - /,;«;• Ы-

(3-53)

Е

§ 4. ТЕОРИЯ МАЛЫХ УПРУГОПЛАСТИЧЕСКИХ ДЕФОРМАЦИЙ

Теория малых упругопластических деформаций разработана Л Л Ильюшиным для простых процессов. В пятимерном прост- |>111н,тне Ильюшина И5 такой процерс

еч(<) = 4 х (*)

(4.1)

описывается прямой,

причем интенсивность тензора .деформации

е и является «длиной

дуги» этой прямой в:

€„(<) = \]е%еач |х(*)| = в.

(4.2)

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

 

 

+ 1 ,

если х > 0,

 

 

ф = 31§П X =

< 0,

если х = 0,

(4.3)

 

 

1 ,

если х < 0,

 

и обозначим

 

 

 

 

 

 

 

= <р<ти.

(4.4)

I-

5*

и -----------------------

 

 

 

0

 

 

 

 

Рис. 5

Рассмотрим теперь процесс деформации, характеризующийся изменением длины дуги (рис. 5), где в,- — точки, соответствующие изменению направления процесса. Тогда величина связанная с процессом, в котором х (<) меняет знак п раз, записывается в виде

г*

-

±8 1 + 2 Ё ( - 1)'Л(8 - 8')

, (4.5)

'и(м)

=

|=1

где Н(в —в,) — единичная функция Хеви­ сайда, а знак перед выражением правой части (4.5) определяется первоначаль­ ным направлением процесса. Пусть е* — последний достигнутый предел текучес­ ти, т.е. значение е$, соответствующее последнему из 8„, когда началась разг­ рузка, а е*7- — соответствующий этому моменту девиатор тензора деформации.

Пусть е** — предел текучести, который может наступить или на- ' I Vпил при протекании процесса в одном направлении (на рис. 6 для момента 54 < в < «5 этому пределу соответствует значение в*, а для в > «5 — значение «4). Тогда процесс нагрузки (значение депиатора тензора напряжения) описывается следующим образом:

вц = 2/|[еу - ы(е*)еи] +

(4.6)

- Ц е и)е,7]й(еи - еП -

Из соотношений (4.6) следует сразу формулировка определяющих соотношений отдельно для активных и пассивных процессов, дан­ ная в § 4 гл. 1. В случае активных процессов получаем уравнения равновесия в перемещениях в виде

ц(1 - и)Ащ +

^ ( 1 ) +

0,- -

+ «,-.,-)+

 

 

+2*^4» +

0.

 

(4.7)

и для пассивных — в виде

 

 

 

р А щ +

-|- К ^0 '4 - р У * + рР{ = 0.

(4.8)

.Чдесь

 

 

 

 

у ,*= ы(е;)д « : + ^

 

 

(4.9)

1

Ли (

2 „

\

 

ш,к"

 

+ щ'>киз* ~ 1 вв>к)>

(4.ДО)

нм чдочкой помечены все величины, относящиеся к началу разгРУ ши (всего тела). Следует добавить еще граничные условия

щ кг = И?, [*уп,- + К вп{]я2 = в?.

(4.11)

И |и к, статическая задача теории малых упругопластических де­ формаций заключается в решении системы дифференциальных урнипоний (4.7) или (4.8) при удовлетворении граничным услови­ ям (4.11). Бели рассматривается несвязанная задача термоплас- 1ичп(1гти, то массовые и поверхностные силы в (4.7), (4.8), (4.11) •подует изменить согласно формулам (1.71).

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

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

Упражнение 4 .1 . Доказать теорему о разгрузке [27, с. 118]. Если всюду в теле происходит разгрузка, то перемещения точ­ ки тела в некоторый момент стадии разгрузки отличаются от их значений в момент начала разгрузки на величины упругих перемещений, которые возникали бы в теле, если бы в естест­ венном (ненапряженном и недеформированном) состоянии к нему были приложены внешние силы, равные разностям внешних сил, действующих на тело в указанные моменты.

Упражнение 4.2. Доказать теорему об остаточных напряже­ ниях, деформациях и перемещениях [27, с. 120]. Если для тела при

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

Мы видим, что задача теории малых упругопластических де­ формаций при пассивной нагрузке (4.8), (4.11), по существу, яв­ ляется задачей теории упругости. Единственность ее решения была доказана в § 1. В случае активной нагрузки решение задачи (4.7), (4.11) будет единственным, если выполняется неравенство (7.60) гл. 1:

6Ц] ^ 2рбецбеч = ^(<$ги)2 + ^(<56>)2|,

р > 0.

(4.12)

Используя (4.55) гл. 1, получим

 

 

 

> 2Л [(1 - и >)6а6п - ^

- е- ^ 6 е

к1 =

 

д е к1

 

 

 

 

=

|(1 - ы)6еу -

 

(413)

Поэтому

 

 

 

 

\д<П] Ьеы

 

л

 

(4.14)

I - ш - еи— - (6еи)2 + К (ё в )

\деы

Из сравнения (4.12) и (4.14) видно,что для единственности решения достаточно потребовать, чтобы

Л<ти

<1ы \

 

Ле»

1 ы ' Ж / > 0, К > 0.

(4.15)

Первое требование (4.15) можно записать в виде

ш + е и- — < 1,

(416)

а второе — в виде

(4.17)

Введем упругопластический потенциал деформации ЧУ и потенци­ альную энергию деформации <р:

&V

 

 

ЧУ = I <тис1еи + ^Кв2,

= ^ Ш У .

(4.18)

о

V

 

Упражнение 4.3. Сформулировать вариационный принцип Лагранжа для задач (4.7), (4.11) и (4.8), (4.11) и доказать, что стационарная точка лагранжиана является точкой минимума. ■

Соотношения (4.6) можно разрешить относительно деформа­ ций. Тогда

(419)

<ти

при активных процессах и

Ч = еЬ - ^ ( Ч - * Ъ )

(4-20)

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

еи = еи{<ти).

(4-21)

Тогда можно ввести упругопластический потенциал напряжения и> и потенциальную энергию напряжения Ф следующим образом:

<т«

(4.22)

го</<7° + у к * 2’

о

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

= 0

(4.23)

^л-г г н и «_ л У1

Ае,'у — еис,к] е^к,ы д[^,0 "Ь ^уД0]>

(4.24)

выраженные через напряжения, и граничные условия

= 5?.

(4.25)

Упражнение 4.4. Сформулировать для задачи теории ма­ лых упругопластических деформаций вариационный принцип Кастильяно и доказать, что стационарная точка Кастильяно является точкой максимума, если выполняются неравенства (4.17) и

4-^- > 0,

ц > 0.

(4.26)

и(Гу

 

 

Упражнение 4.5. Доказать теорему о простом нагружении, свидетельствующую о том, что существуют реальные среды, для которых можно выбрать входные данные таким образом, что в каждой точке среды одновременно осуществляется простой про­ цесс нагружения [27, с. 115]. Если материал несжимаем (0 = 0) и интенсивности тензоров напряжений и деформаций связаны меж­ ду собой по степенному закону

= сеС,

(4.27)

где с и п — некоторые постоянные, и, кроме того, массовые и

поверхностные силы рГ и 3° возрастают пропорционально одно­ му параметру х х , а заданное перемещение й° пропорционально другому параметру х„, причем

(* х )" =

(4.28)

то процесс деформации в каждой точке будет простым.

§5. ТЕОРИ Я ЛИНЕЙНОЙ ВЯЗКОУПРУГОСТИ

Влинейной теории вязкоупругости связь между напряжениями

идеформациями задается соотношениями (4.3), (4.6) § 4 гл.1, которые можно записать сокращенно в операторном виде:

«

 

<П)= У>Гум(*,т)еы(т)<1г = Гуыеы,

(5.1)

о

 

*

 

= /^•ТЬ,(*>Г)^«(Г)^Т =

(5.2)

о

 

или для случая ядер релаксации и ползучести разностного типа в виде соотношений (4.13) и (4.14) гл. 1:

<

(5.3)

Я цы Ц - т)<1т(т) = В.ЦЫШ,

о

 

^ П,7*((<, г)(1<Гы(г) = П1}к1<?к1-

(5.4)

О

 

При этом тензоры ядер релаксации различных типов анизотро­ пии представлены соотношениями (4.22), (4.24), (4.26) гл. 1. За­ метим, что сокращенные соотношения (5.1), (5-3) напоминают по записи соотношения Гука (2.2), причем в последних вместо тензо­ ра модулей упругости следует подставить тензор-оператор

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

РЩ = рРг + [?у-ы(2)«мЪ

(5.5)

при выполнении граничных условий

««Ь, =

Г,-7*|(2)и*,|П7-|Е2 - 5 °

(5.6)

и начальных данных:

при 1 —0 щ = V,, щ = 1$,

(5.7)

где запись Гуы(Й) отражает тот факт, что ядра релаксации могут зависеть от координат. Бели они не зависят от координат, то вязкоупругая среда называется однородной.

Упражнение 5.1. Показать, что квазистатическая задача ли­ нейной теории вязкоупругости заключается в решении трех урав­ нений равновесия

[ГУ,ы(*)«*,|Ь' + РЪ = 0

(5-8)

при удовлетворении граничным условиям 5.6. В В частности, для изотропной однородной вязкоупругой сре­

ды квазистатическая задача заключается в решении уравнений

равновесия

 

+ й)0,{ + рАщ + рР% = 0

(5-9)

« •к = И°, [А0п,- + /!(«,•,; + и;>)п;]2а = 5,?,

(5.10)

где А, Д — соответствующие скалярные интегральные операторы. Бели они соответствуют ядрам разностного типа, то будем их обозначать А и /х. Положим

2/| = Г,

2/х = А, А - ф

= Гх, А + |/х = Аь

(5.11)

т.е.

 

 

 

 

} = Геу =

[ Г(*,т)еу(т)А-,

«0- = Ае„ =

[ А(* -

т)Ае^(т),

 

«>

 

*

(5-12)

 

«

«

 

<г = Ггв =

^ Г х (*, т)0(г)</г,

«г = Дх0 = у

Дх(« - т)<*0(т). (5.13)

 

о

о

 

 

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

^Г

Ш = I

А(*>г) { / ц(т,п)Ят1)<1т1^<{т,

(5.14)

о

о

 

1 т

А/2/ = У А(< - т)Ау /х(г - п)А /(п) =

оо

<г

У - т)А ^ ц{т - Гх)с//(гх) = (хХ/. (5.15)

Единица, деленная на данный оператор, обозначает оператор,

обратный данному,

т.е. если

 

 

 

 

т> 1

ту

1

А 1

А

1

(5.16)

К — —,

К\ — т~,

П —

Пх =

,

Г

 

Гх

А

 

Дх

 

то из (5.12) и (5.13) следует соответственно

 

 

 

е»х

К 3)],

ву — Пвц,

 

(5.17)

 

0 =

Кх<г,

0 = Пх<г.

 

(5.18)

Исходя из «основных» операторов Г, Гх (Я ,Я 1) или К ,К \ (6 , 6 1 )

можно построить

другие

 

 

 

 

 

 

 

и’1= гГА'1,

Ш2 =

- К \ Т ,

Т1 = — ,

Я 2 =

_1

_

 

А >

 

6

 

О

 

Ш1

 

2

(5.19)

 

1 - -

. .

.

1

.

1

—й

 

 

О) = О»! = 0 )2 =

Г Т ? П 1 ,

1Г = ТГ1

= ТГ2

= — ,

V =

2

+ и>

 

 

О

 

 

й)

 

 

Если теперь, пользуясь алгеброй операторов, мы получим фор­ мальное решение задачи (5.9), (5.10) или (5.8), (5.6), то для по­ лучения решения задачи линейной теории вязкоупругости для однородных сред будет необходимо «расшифровать», встречаю­ щиеся в решении функции от операторов. В этом и состоит принцип Вольтерры. Следует иметь, однако ввиду, что в случае ядер релаксации и ползучести неразностного типа умножение опе­ раторов не является коммутативной операцией,, и поэтому при использовании принципа Вольтерры нужно проследить за мето­ дом получения аналитического решения соответствующей задачи теории упругости с тем, чтобы правильно записать произведение упругих постоянных, входящих в ее решение. Основная трудность при решении указанных задач возникает при «расшифровке» опе­ раторов. Для упрощения этой процедуры часто «основные» операторы выбираются в специальном виде, а экспериментально найденные ядра релаксации и ползучести аппроксимируются яд­ рами, соответствующими данному специальному виду этих опера­ торов [99]. Для случая ядер разностного типа часто применяется метод преобразования Лапласа [33]. При «расшифровке» вяз­ коупругих операторов большое значение имеет так называемый оператор А.А. Ильюшина др:

1

х

1

9 Р ~ 1 + /Й1’

* * ~ 1

+ Р ш

где (3 — некоторое действительное число. Ядро, соответству­ ющее оператору др, может быть найдено из эксперимента на релаксацию образца, к которому последовательно присоединена пружина заданной жесткости в зависимости от числа /? при нерелаксирующем объеме [33]. Если же и сдвиговые и объемные свойства материала зависят от времени, то вместо пружины сле­ дует присоединить скручиваемый образец, изготовленный из того же материала, что и основной образец [70]. В этом последнем случае тонкостенный образец « 1 » толщиной 6 и длиной скру­ чивается силой (?({), которая одновременно растягивает другой образец « 2» из того же материала, имеющий длину рабочей час­ ти /2 и площадь поперечного сечения Р^. Если теперь задать перемещение образца в виде

« = Щ ( 0 ^

(5.21)

и снимать показания динамометра об изменении нагрузки <3(2), то ядро др(1) для 0 ^ ^ | будет найдено по формуле

*/>(*) = 1 - Я Н ) ,

/? = 22/ ( ^ / 1 + 222) , ^ = *(62 - о2), 6 - о = 5. (5.22)

Для определения ядер др(1) при других значениях /? можно зак­ репить снизу образец «2», а нагрузку ф(2), которая скручивает образец «1», приложить сверху к образцу «2», заставляя послед­ ний работать на сжатие. В этом случае

/? =/2/ (222- ^ 21) .

(5.23)

Упражнение 5.2. Показать, что справедливы следующие фор­ мулы:

ЯрЛр, ~ р _ ^ (01ЯРг ~ РъЯр,),

(5.24)

йяр= ^(1 - яр)-

(5.25)

Упражнение 5.2. Используя формулу (5.20), показать, что для экспоненциального ядра

 

ы(2) = А + Ве~а%

(5.26)

ядро др(1) также является экспоненциальным:

 

дрЦ) = С + О е~ *,

(5.27)

здесь А, В , С, Б, а, у

— постоянные, причем

а О + ^ А

) _ ------- д _

В ?

1 1 + 0А + Р В ’

1 + р А ’

(1 + 0А){1 + 0А + /ЗВ)

 

 

(5.28)

Упражнение 5.4. Пользуясь результатом предыдущего упраж­ нения, показать, что ядро тг(2) (5.19) также является экспоненци­

альным:

 

тг(*) = А1 + В 1е01<,

(5.29)

 

 

причем

 

 

 

«1

а А

В

(5.30)

1 + В '

В г = -

 

А(А + В)

 

Соседние файлы в папке книги