книги / Оболочки и пластины
..pdfд ф = 0 |
при х = О, х = L, |
дхду |
|
W (х, у, 0 ) = w0 (х, у),
ф(*. у, 0) = ф0 (х,у).
Будем считать, что кривые релаксации напряжений в материале обо лочки могут быть аппроксимированы функциями вида
<*ц( 0 = М 1 Уu— ( 1 —
Опыты на релаксацию растянутых образцов стеклопластиков показы вают, что некоторые параметры Yij и atj для образцов, по-разному ориентированных к основе ткани, могут быть с достаточной для инже нерной практики точностью приняты одинаковыми:
Yu = Y22*» а1Х= а12 — а22 — а -
Так, для полиэфирного стеклопластика (ПН-1, Т-1) имеем Оц = а12 = а22 = а = 0,5,
Yu = Y*a=-g-’ Y12 = M'l’Yii> Y= 3Yu.
(3,22,25)
*ц = 1,49-10»-^-, 622 = 1,14-Ю5- ^ - ; bn — bn = 0,18-10®
6 = 0,13-105
Примем
Rn = R22 = y№~at — |
Ri2= H i*Я н |
R = 3Ru, |
(3,22,26) |
а в преобразовании Лапласа |
|
|
|
= |
= |
Зуа |
(3,22,27) |
где \i\ есть неизвестный параметр. Представим прогиб w(x, у, /) в виде
w (х, у, 0 = ш0 (л:, у) о»!(0, |
(3,22,28) |
где w0(x, у) есть упругое решение, а W\(t) выбирается в форме
w1(t) = 'K— \ie-at. |
(3,22,29) |
Параметры X и р, подлежат определению. Очевидно, что X—ц=1. Функ цию напряжений возьмем в виде
ф(*, У, 0 = Фо(*, y)-q>i(0,
где ф0 (JC, у) есть упругое решение задачи, а <pi(0 при заданной a>i(f) найдем из уравнений (2,24,16) и (2,24,12). Таким образом, задача сво дится к построению упругого решения w0, <р0 и определению параметров jxi, X и ji. Полагая в уравнениях (2,24,12) и (2,24,16) &i= 0 и <7=0, полу
|
|
|
р _ |
— |
|
еал |
|
|
|
|
|
|
|
|
40аТ1*п |
|
1 |
1 |
к |
(3,22,34) |
|
|
|
C l |
|
Схб» |
\ |
еатА |
||||
|
|
|
|
|||||||
|
[ - ( • i66?,e2 |
2^11^2 (в) / |
|
gi (в)+ |
g2(в) |
|
||||
|
|
_ М |
й р |
+ Л .С2 __ 1 _ Сг==Со> |
|
|
(3,22,35) |
|||
|
|
|
Ci |
|
8 |
2 |
|
|
|
|
Из |
уравнений |
(3,22,33) —(3,22,35) |
можно построить |
зависимость Р— |
||||||
для |
различных |
значений |
параметров |
волнообразования т) и 0. |
|
Результаты вычислений для материала (3,22,25), (3,22,26) пред ставлены на графике (рис. 3.76) в случае 0= 0,5 и -л=0,22. Определяя критические нагрузки, получим для верхнего значения
Рв = 0,558,
а для нижнего
РИ= 0,171.
Из анализа значений вторых производных от да0 и <р0 по х и у в точке (x=L/2, у) легко видеть, что для m-нечетного решение задачи о влия нии времени на прогиб получается без затруднений. Рассмотрим случай /n-четного. Преобразуя по Лапласу уравнение (3,22,28) с учетом форму лы (3,22,29), имеем
w'(p) = w0( - -------- |
(3,22,36) |
Ч р |
р + о / |
Внесем (3,22,36) в уравнение (2,24,12). Решая полученное уравнение относительно <pi, получаем после инверсии
<Pi = Ах-|- A1e~at + A3e~'>ai + A4e~2al + Abte~at, |
(3,22,37) |
где
Л = А ( * 1 » . “'о, Фо. К ц. Hi).
Внося ср* и |
w\ в уравнение |
(2,24,16) |
и |
учитывая (3,22,27), получаем |
|
после инверсии преобразования Лапласа |
|
|
|||
Я 1 + К |
2е - “ +' K3e~‘/,at + |
Кле-'Ч^ + |
Kse~2ai + K6e~™ + |
K7te~^ + |
|
|
|
+ K e t e - * * |
= |
Q , |
(3,22,38) |
где Ki = Ki(bn, w0, фо, К ц, |xi), причем в К\, /С2, Кз, Кб и К7 входит ли нейным образом нагрузка Р. Используя метод коллокации и полагая
в уравнении |
(3,22,38) |
Р = 0,171, £i = 5,396, £2=2, получаем для |
материа |
|
ла |
(3,22,25), |
(3,22,26) |
значения Ц1=0,12, К= 1,242, (х=0,242. |
|
на |
Таким образом, величина рм имеет порядок коэффициента Пуассо |
|||
[155] для |
рассматриваемого материала (для полиэфирного стекло |
|||
пластика ПН-1, Т-1, |
v = 0,11; vi2 = 0,12 [153]). В дальнейшем |
значение |
jii = 0,12 сохраним для всех Р. Рассматривая прогиб для различных зна
чений ^=const на участке 0<Р<Рь (£i = 0, £2= 0), получаем методом коллокации (полагая /= 0, 1, 2, 3, 4, с») систему уравнений, из которой
с помощью формулы Крамера найдем Х=1,2. На участке 0<Р<Рь из менение прогиба во времени при постоянных нагрузках будет выражать ся формулой
w (x,y t) = w0(x,y)(l,2 — 0,2е-°.5'),
т. е. прогиб может изменяться по сравнению с упругим на 20%, если /—>-оо. Предполагая нагрузку изменяющейся во времени
|
P = PoPi(t), |
|
|
|
|
||
где P\{t) есть безразмерная |
функция времени, |
такая, |
|
что |
P j(0 )= l, |
||
имеем |
|
|
|
|
|
|
|
д2<Ро(*. У) |
= Фх(0 = — P*Pi(t), |
|
|
||||
д у* x = 0 ,L |
ду* |
|
|
|
|
|
|
откуда 9i(f) —Pi(t), или, учитывая (3,22,37), |
|
|
|
|
|||
- ± t |
- |
± |
t |
- |
L |
t |
(3,22,39) |
Р — Ро {А -}- Л2е 2 |
А$е |
12 |
-)- |
1-f- A$te |
2 |
). |
Время критического состояния можем найти из условия dP/dw = 0, учтя, что прогиб задан в форме
|
W (х, у, t) = |
w0 {х, у) (К — р,е |
4 * )• |
|
|
|
Полагая |
оо, получаем уравнение для tкр: |
|
|
|
||
|
• А4е~‘ + Аье 2 * |
1 |
—Y |
t |
||
2 |
|
2 |
= 0 . |
|||
12 |
|
|
|
|
|
|
Это уравнение можно решить |
приближенно, |
представив |
eat в виде |
|||
ряда |
|
|
|
|
|
|
|
l + a t + |
Ш |
+ H W |
|
|
|
|
|
2! |
3! |
|
|
|
Для Х= 1,242, ц = 0,242, Ро = 0,171, ш0=1,76 (получим кубическое уравне
ние, из которого по формулам Кардана найдем |
fKp = 7,058. Далее най |
дем |
|
|
7,058 |
рн = 0,1374, wH= w0(1,242 — 0,2426 |
2 ) = 2,17. |
Таким образом, учет линейной наследственности снижает значение кри тических нагрузок для оболочек из стеклопластика.
Упругое решение РВу = 0,558; с наследственностью Рвн=0,470. Упру гое решение Рну=0,171; с наследственностью Лш=0,137. На рис. 3.76 дан график зависимости нагрузка—прогиб: а) упругое решение; б) из менение прогиба во времени при t->оо и Р = const; в) с учетом наслед ственности при нагрузке, изменяющейся во времени. Заметим, что кри тические нагрузки оболочек, материал которых обладает линейной на следственностью, существенно зависят от режимов нагружения, увели чиваясь с увеличением скорости нагружения.
3. Применение табулированных сингулярных ядер в задачах
устойчивости упруго-вязких оболочек
В предыдущих пунктах при решении задач использовалась простей шая функция влияния — ядро релаксации в виде экспоненты. Однако, как уже отмечалось в § 23 гл. II, в начале релаксационных процессов наблюдаются весьма большие скорости, которые не могут быть описа ны экспоненциальной функцией или их суммой и довольно хорошо ап проксимируются сингулярными функциями влияния типа (2,23,17) или (2,23,18). Однако их применение в расчетах ползучести элементов кон струкций ограничивалось трудностью вычисления интегралов от них, появляющихся в решениях. Эти трудности в значительной мере снижа ются при использовании таблиц сингулярных функций влияния и их интегралов [156]. Покажем, что решение задачи об устойчивости гиб ких упруго-вязких оболочек может быть выражено через табулирован-
t |
t |
ные ядра R, К, ^ R(u)du, |
j*K(u)du, где R и К даны формулами |
6 |
о |
(2.23.12) и (2,23,17). Задача сводится к решению системы уравнений (2.24.12) и (2,14,16) при определенных граничных условиях. Примем их такими же, как и в п. 1 для шарнирно закрепленных кромок обо лочки.
Приближенное представление прогиба во времени примем в виде
w(x, у, t) = w0(x, |
w±(t) = X — \ie-<x”t — yQ(t — tk), |
(3,22,40) |
где |
|
|
|
9 < < - у = { ° |
|
Считая возможным |
представить функцию напряжений |
в виде |
ф(*>У, t) =фо(*>#)<Pi(0 |
и внося w*(p) в ур.авнение неразрывности |
(2.24.12) , получим уравнение, разрешая которое относительно ф1*, по
лучим зависимость, |
аналогичную |
(3,22,13), которую здесь |
запишем |
в виде |
|
|
|
ф; = Qi (х, У) Ю * - |
Qi (х, у) R* (w*y - |
Q2 (*, У) vf[ + Q2 (х, у) R‘w\, |
(3,22,41) |
где
Q2 у) — |
|
(Ь цЬ п - Ь М | \ |
Hr h ^ |
(3,22,42) |
|
д*Фо |
1 |
|
^Фо |
||
|
|
^Фо |
|||
|
bil ~д* |
* 2Ь |
{bllh* ~ 612621“ |
466ls) д*д? + |
622~ду* |
Q l ( x , у ) = |
|
(6ll622~ |
6 12 6 2 0 |
[ ( l a y ) |
~ ( дх2 |
ду2 |
)] |
||
|
|
|
|
|
|
|
d2w0 |
*d2w0 |
|
, |
^Фо |
, |
1 |
„ |
, |
t , |
ч ^фо ' , , |
^Фо |
|
611 |
д х * |
+ |
2b |
( и |
22~ |
12 21 “ 6 |
l2) дх* ду> ф |
22 ^ |
|
величины, определяемые по упругим решениям. Используя теорему умножения операционного исчисления [157], согласно которой, если интегралы
Ш = ] м )е п * < и и r 2 ( p ) = |
] f 2 ( t ) e - p ‘ d i |
6 |
о |
сходятся абсолютно в области Rep>yi, то функция
П (Р ) = П ( Р ) П ( Р )
есть преобразование Лапласа функции
|
|
/ з |
( 0 = j f i ( i |
— s)ft (s)ds, |
|
|
|
||
|
|
|
О |
|
|
|
|
|
|
а сходимость интеграла |
|
|
|
|
|
|
|
|
|
|
|
n(p) = ] fs(t)e -p‘dt |
|
|
|
||||
в области Rep>yi абсолютная, — запишем (3,22,42) |
в виде |
|
|||||||
где |
ф; ср) = Qi к г - |
Q i^; - |
Q2«>;+ |
Q2^ |
; |
(3.22.43) |
|||
t |
|
|
|
|
< |
|
|
|
|
|
s) ш2,(s) ds, |
J£2 = |
|
|
|
||||
|
= j R (i _ |
j R (/ — s) tj»j.(s) ds. |
(3.22.44) |
||||||
Так как |
c j k (f)j* = |
£ |
c„/‘ (p), |
то, обозначая |
|
|
|
||
|
A=1 |
A=1 |
|
|
|
|
|
|
|
|
ф*(p) = ф* (p) — Qi (cei?)* H- Qi-%1 + |
Q2W1 — QH$2 = 0, |
(3,22,45) |
||||||
и пользуясь теоремой [157] о том, что «если интеграл |
|
|
|||||||
|
|
ф* (р) = |
j e-pt?i(0 dt |
|
|
|
|||
сходится в точке ро и во всех точках pk=Po + kl |
(k—0, 1, 2,...), |
функция |
|||||||
Тр, (рА) =0, то <pi(0 равна |
нулю в интервале (0, |
оо) |
всюду, за |
исключе |
нием точек разрыва функции ф!* (р)"», — получим из (3,22,45) уравнение
Фх (0 = <Pi(0 — Qi (*. у) wi (0 + |
Qi (*. у) |
(0 + |
Qz (•*. у) щ (0 — |
||||
откуда |
- - -0.2 (,Х > |
у ) |
( 0 = |
О , |
|
|
|
|
|
|
|
|
|
|
|
Фх (0 = Qx(*. lJ)w\(t) — Q2 (х , у) wl {t) — Qi (л:, у) |
(0 + Q2 |
(х, у) |
(t). |
||||
|
|
|
|
|
|
|
(3,22,46) |
Внося сюда выражения |
для,#] |
из |
(3,22,44) |
и учитывая (3,22,41),. |
|||
получим решение для <pi (t) в виде |
|
|
|
|
|
|
|
Фх (0 = Qi (Ь - |
- Y0 (t - |
У )2 - |
(Q^2 -'Л ) j Я (i - |
s) ds |
f |
||
|
|
|
|
|
о |
|
|
+ (2pX — p) J |
|
|
f |
|
|
— 2[ry I e~a\sQ(s — sk) R(t — s) ds. |
(3,22,47) |
|
о |
|
|
В это выражение входят интегралы |
от произведения ядра |
R на неко- |
|
t |
|
торые функции. Поскольку функции |
R и . | Rds протабулированы, по* |
|
|
о |
|
лезно эти интегралы свести к табулированным функциям, чтобы избе жать трудностей вычисления интегралов типа
t
^ R (t — s) Ф (s) ds,
о
где в качестве* R (t) |
выбрана сингулярная функция типа |
|||
|
R(t) = Ае-Щ«-' |
(0 < |
а < 1). |
|
В самом деле, встречающиеся в |
(3,22,47) интегралы, если учесть, чтск |
|||
Q2(s—sk) = 0(5—5ь), будут равны |
|
|
|
|
t |
|
t |
R(t — s)ds— j* R (t — s) ds, |
|
I*02(s — sk) R (t—*s)ds= |
j* |
|||
6 |
|
0 |
|
0 |
t |
t |
|
|
*k |
I*g-ais 0 (s — sk)R(t — s)ds = f e~a'sR (/,— s)ds — j e~aiSR (t — s) ds. |
||||
6 |
o |
|
o |
t |
|
|
|
|
|
Последние интегралы выражаются через/? и |
J /?^следующим образом: |
|||
|
|
|
|
о |
= j e-a'sR (t — s)ds |
|
<) |
|
|
|
о |
|
|
|
|
|
0 |
|