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

книги из ГПНТБ / Кристенсен, Р. Введение в теорию вязкоупругости

.pdf
Скачиваний:
38
Добавлен:
22.10.2023
Размер:
10.86 Mб
Скачать

§ 2.9. Цилиндр под действием внутреннего давления

77

m

°zz(r,f) = — р о ' £

[p ln r (cij'fB (tZj) -f E (fl/)],

Hm [Д (s)/(s — а/)]

/= 1

 

 

(2.87)

где комплексные корни а/ уравнения F ( s ) = 0 можно оп­ ределить на ЭВМ с помощью стандартной программы, когда в представлении для р(^) заданы характеристики какого-либо конкретного материала.

Квазистатическое поведение сжимаемого цилиндра

Для сжимаемого вязкоупругого цилиндра при осевой симметрии единственное нетривиальное уравнение квазистатического равновесия в цилиндрических координа­ тах г, Ѳи z имеет вид

+

_ L

^

_ J L = o,

(2.88)

д г 2

г

д г

г

где u = u(r, s ) — преобразование Лапласа от радиально­

го перемещения. Решение уравнения (2.88)

будет

и = Cr D/r,

(2.89)

где C(s) и D(s) — функции переменной преобразования, подлежащие определению.

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

(2.90)

(2.91)

(2.92)

где

е = (ди/дг) -f- (и/г).

(2.93)

78

Гл. 2. Изотермические краевые задачи

Граничные условия для этой задачи те же, что и для предыдущей, и задаются зависимостями (2.78) с уче­

том (2.79). Постоянные С и D определяются из преобра­ зования Лапласа от (2.78) и (2.79). Оказывается, что

 

 

o rr =

C i (b2 г2) С2,

 

 

 

оѳѳ -

C ,+ (b */r*)C 2,

 

 

 

azz =

2svCu

 

где

 

 

 

 

п

 

Р(К sp)

 

^ 1

(К-— sp) +

(b2/a2)[K (1 — 2sv) +

sp]

 

 

г

p [/C (1 — 2 s v )+ s p ]

 

и 2

— sp) +

(b2/a2)[K (1 2sv) +

sp]

 

 

причем

 

 

 

 

 

 

~

2 * ( 1 — V *). ‘

 

(2.94)

(2.95)

(2.96)

(2.97)

(2.98)

(2.99)

Если ц(7) и k(t) принять по-прежнему в виде (2.68),

то р. и V будут отношениями полиномов относительно па­ раметра преобразования s. Примем, например, p(t) в ви­ де ступенчатой функции времени (2.85) и подставим ее

преобразование и преобразования p,(s) и v(s) в (2.94) — (2.99). Это приводит к преобразованиям напряжений, которые содержат отношения полиномов по s и могут быть обращены с помощью метода, изложенного в при­ ложении Б. Выполнив эту процедуру, считая вязкоупру­ гий коэффициент Пуассона действительной константой

и принимая sp(s) в форме (2.80), получим для агг вы­ ражение вида

Orr = Po[C(s) — (b2/r2)D (s )]/F (s ),

(2.100)

где

С (s) = — КВ (s) + А (s), D (s) = (1 — 2v) KB (s) + А (s),

F (s) = s {KB (s) — А (s) + (bW )[( 1 — 2v) KB (s) + А (s)]}.

Функцию F (s) можно разложить на множители соглас-

§ 2.10. Действие давления на границу сферической полости 79

но (2.86). Обратное преобразование Лапласа зависимо­ сти (2.100) дает

 

fC (о/) - (ЬѴг») D (aj) 1 еV

Огг 0

p°S

lim [Д (s)/(s-— fly)]

(2.101а)

 

/=1

 

 

 

Другие две компоненты

напряжений

в

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

с (2.95) и (2.96)

оказываются равными

 

fl;/

 

 

 

 

°ѳѳ 'Г>0

~ Po У '

[С (g /)+ (fca/r*) D (fl/)) e ;

lim [F (s)/(s — a.j)]

 

 

/=1

s->ay

 

(2.1016)

 

 

а

<

 

 

 

 

 

2vC (a,j e

1

 

<Дг (/■. 0 = Po 5 ] lim [F (s)/(s — aj)]

 

l—‘ s->a-

Примеры вычисления напряжений по формулам (2.101) дали Кристенсен и Шрейнер [2.5]. В этой работе функ­ ция релаксации р(^) представлялась рядом (2.68) убы­ вающих экспонент и находились соответствующие корни F (s ). В работу [2.5] включен также один частный спо­ соб учета влияния границы, которая изменяется вслед­ ствие сгорания или коррозии материала.

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

§ 2.10. Действие давления на границу сферической полости

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

80

Гл. 2. Изотермические краевые задачи

среде. В отличие от предыдущего примера эта задача используется для иллюстрации использования преобра­ зования Фурье. При этом учитываются инерционные чле­ ны; таким образом, рассматривается задача о динами­ ческом поведении, включающем эффект распространения волн. Этот эффект не будет, однако, рассматриваться здесь детально; основное внимание будет уделено мето­ ду решения задачи. Специальное исследование вопроса о распространении волн будет проведено в гл. 4. Реше­ ние, приведенное здесь, взято из работы Локкета [2.23].

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

д2и

. 2 ди

2и _

р

дЧі

 

 

dr2

г дг

г2 _ А + 2ц ä 2 ’

'

'

где и = и (г, /) — радиальная

компонента перемещения,

а X и ц — постоянные Ламе.

Следуя методу, изложенно­

му в § 2.5, превратим это уравнение движения в преобра­ зование Фурье уравнения движения вязкоупругой среды

заменой и (г, t)

на и (г, со), К на К* (іа) и ц на ц*(гсо). Это

дает уравнение

 

 

 

 

 

 

а* 5

,

2

d“

2ц = —Й2со2«,

 

(2.103)

дг2

^

г

дг

г2

 

 

где

 

 

 

 

 

 

Q2 =

р/[Я*(г<о) + 2ц* (гео)].

 

(2.104)

Решение уравнения

(2.103)

имеет вид

 

 

 

д

'

С (со) еШаг + D (со) e - iQ(ür

'

(2.105)

и = ~я~

 

 

Г

 

 

дг

 

 

 

 

где функции С (со) и D (a) подлежат определению. По­ скольку. внешней границы не существует, не может быть и отражения, вызывающего приходящие волны; в силу этого могут существовать только исходящие волны и С((о) должно обращаться в нуль. Это приводит (2.105) к виду

« = — (Цт2)(1 ! tQcor) е~ іаслг-

(2.106)

Ненулевыми являются только компоненты напряжения огг и сгѳѳ= а фф. Используя преобразование соотношения

§ 2.10. Действие давления на границу сферической полости 81

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

(2.106) в виде

 

 

 

Ъ„ - (Dir3) [4fx* (1 + iQ w ) — pcoV2J е~ійш

(2.107)

и

 

 

 

■ Стѳѳ = — (Dir3) [2ц* (1 + a im )

+

Q2coV21*\ е~ аыг.

(2.108)

На границе полости г = а

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

огг(г, а ) = р(ы ),

г = а,

(2.109)

где р — преобразование Фурье

заданного давления на

поверхности полости. Используя (2.107) и (2.109), имеем

дЗ П

 

D = 4ц* (1 + iQ m ) рш2 а2

(2.110)

Получающиеся после подстановки (2.110) в (2.106) — (2.108) формулы удобны для обращения. Формулу обращения преобразования Фурье (1.73) можно перепи­ сать в иной форме

оо

f (t) = —

Re Г / (м) еш dar,

(2.111)

ЗХ

,}

 

 

о

 

где через Re обозначена действительная часть. Напри­

мер, применяя (2.111)

к обращению перемещения грани­

цы полости и (a, t),

находим его в виде

 

и (a, t) =

R i'

(1 і®ша) реш da

(2. 112)

 

я

^ е .)

4ц* ( 1 -f-ißcoa)— pco2a2

 

Интегрирование в (2.112) можно выполнить, если из­ вестно давление p(t) и величины ц* (ш) =ц '(ш ) -Нц'^®) и К* (too) = А/(<о) -[-ІА," (со) заданы как функции частоты ш для конкретного вязкоупругого материала. В общем случае провести это интегрирование аналитически невоз­ можно, однако практически можно проделать его чис­ ленно. В качестве примера зададим такой закон измене­ ния давления, для которого численное определение (2.112) упрощается.

6—851

82 Гл. 2. Изотермические краевые задачи

Примем давление на границу сферической полости

равным

 

 

ап

 

 

p ( t ) = - PЛ± J

(2.113)

и,

 

 

Используя определение интегрального синуса

 

U

 

 

Si (и) = I*

dx,

(2.114)

о

 

 

можно переписать (2.113) в виде

 

p {t) = (Ро/п) [Si((o20

— Si (cöi^)].

(2.115)

Преобразование Фурье от этого выражения, определяе­ мое с помощью (1.72), дается формулами

 

О

при

0 < ; со <

од,

 

р(ю)

Po/t©

при

ОД ДС О) <

ОД,

(2.116)

 

О

при

о д ^ ® ^ 00-

 

Получив p ( t ) из (2.115), можно найти интегральные си­ нусы по любым стандартным таблицам, например таб­ лицам Янке и Эмде [2.20]. Схематический график p (t) в зависимости от t дан на рис. 2.3, где од = 5, ©2= 1000.

Р и с . 2.3. Давление на границе сферической полости.

§ 2.10. Действие давления на границу сферической полости 83

Из него мы видим, что давление p (t), вычисленное по (2.115), является приближением к разрывному измене­ нию давления. Преимущество использования (2.115) вместо ступенчатой функции изменения давления состо­ ит в том, что характеристики, получаемые из (2.116), позволяют нам использовать интеграл обращения на ко­ нечном интервале частот вместо обычного бесконечного.

Теперь, используя выражение для р из (2.116), прида­ дим (2.112) вид

 

 

ш2

 

и (а, t) =

Re

(1 -|- t'Qcoa) ieim da»

(2.117)

4м|х* (1 4- ійсоа) — раз2 а2

 

и,

Интегрирование в (2.117) производится в конечных пре­ делах и может быть осуществлено непосредственно, ког­ да p,*(t'co) и Х*(ш) заданы как функции частоты.

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

крешению, выраженному в конечной форме, требующей простого численного интегрирования.

Вобщем случае, когда мы ищем аналитические ре­ шения в переходных динамических вязкоупругих зада­ чах, используя методы интегральных преобразований, процесс обращения связан с точками ветвления. В дан­ ном примере такие точки ветвления вносит П в (1.117). В силу этой трудности аналитические решения удается получить редко. Более полно эта ситуация будет рас­ смотрена в гл. 4.

6*

84

Гл. 2. Изотермические краевые задачи

§2.11. Свободные колебания

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

характер поведения материала, является

класс

задач

о свободных колебаниях. Определением

этого

класса

служит то, что работа всех сил на границе там, где они существуют, тождественно равна нулю; таким образом, граничные условия имеют вид

Ui==n

На

а«'

(2.118)

а ц П/ = 0

на

Во.

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

Преобразование Лапласа уравнений движения дает­ ся зависимостью

Щ,П + (1 — 2sv) 1 Ufji = - J — 2

д 2 щ

(2.119)

dt2

Sfl (s)

 

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

Частный случай постоянного коэффициента Пуассона (ѵ = const)

Сначала напомним анализ свободных колебаний для упругих тел. Уравнения движения имеют вид

§ 2.11. Свободные колебания

85

и по-прежнему применяются граничные условия (2.118). Решение ищется в форме

ul {xh t) = U f{x i) e M .

(2.121)

Функции координат Uf(x.) определяются таким обра­ зом, чтобы удовлетворялись уравнения движения (2.120), а удовлетворение граничным условиям (2.118) приводит к задаче о собственных значениях. Соответст­ вующие собственные значения и собственные векторы определяют собственные частоты и отвечающие им фор­ мы колебаний, которые обозначаются через

С0„ И U f ’n(Х;).

Спектр частот со» может быть дискретным или непре­ рывным в зависимости от того, конечны размеры упру­ гого тела или нет. Здесь и далее в этом параграфе мы будем рассматривать только случай дискретных собст­ венных значений; однако результаты легко можно обоб­ щить и на непрерывный случай. Полное решение при этом представляется суперпозицией собственных векто­ ров, и, согласно (2.121), это дает

=

(2-122)

П

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

Рассмотрим теперь соответствующую вязкоупругую задачу. Предположим, что вязкоупругое решение име­ ет вид

’"(*,) П 0>

(2.123)

где fn(t) — неизвестные функции времени, а Uf-n (х{) ~

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

86

Гл. 2. Изотермические краевые задачи

 

и для вязкоупругой задачи. Подстановка (2.123)

в урав­

нения движения (2.119) дает

 

 

№ + (1

■ 2»Г'ир\ Г =

 

 

 

= (р; sp) [s2 г - f n(0) -

s f (0)j £/f

(2.124)

где

 

 

 

/«(0)

=/»(*)|<=о и /"(0) =

гі/"(ОА«І*=о,

(2.125)

и величины /™(0) и /” (0) определяются из начальных условий. В то же время из соответствующего уравнения (2.120) для упругой задачи следует, что

U f f , + (1 - 2ѵ)—1U f:ü = { - со2 p/pj u f ' n. (2.126)

Подставим (2.126) в (2.124) и найдем соответствующие значения fn (s). Имеем

-п /"(0) + s/” (0)

(2.127)

sa + (cü^sp/p) '

Поскольку (о2п из соответствующей упругой задачи про­

порционально р, решение, определяемое формулой (2.127) , в действительности не зависит от принятого уп­ ругого модуля сдвига р. Принимая вязкоупругую функ­ цию релаксации при сдвиге p(f) в виде ряда из убываю­ щих экспонент, как в (2.68), видим, что (2.127) имеет только простые полюсы, соответствующие которым корни легко найти, исполызуя стандартные программы для

ЭВМ, и обращение /(s) можно провести непосредствен­ но. В точности та же ситуация возникала и в задачах о квазистатическом поведении, где мы сочли полезным оставить в представлении (2.68) для р(^) столько чле­ нов, сколько требуется для достижения точности. Про­ стой пример, включающий обращение зависимости (2.127) , дан в конце этой главы в качестве задачи.

Полное решение задачи о свободных колебаниях име­ ет вид

М * ц* ) = 5 Х ,л(*і) П 0

(2.128)

п

 

при начальных условиях