Рефетека.ру / Промышленность и пр-во

Реферат: Теория оболочек

Введение. Основные определения


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

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

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

Снижение материалоемкости конструкции - важный фактор для многих машин и агрегатов. Выгодно это и в строительных сооружениях. Оболочки позволяют эффективно решать проблему минимизации массы.

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

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

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

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

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

Основная проблема теории тонких упругих оболочек состоит в сведении трехмерной задачи теории упругости к двумерной задачи. Таким образом, развитие общей теории тонких упругих пластин и оболочек идет по пути сведения трехмерных уравнений теории упругости к двумерным. Для решения этой проблемы предложено большое число методов, которые по классификации С.А. Амбарцумяна могут быть объединены в три группы: метод гипотез, метод разложения общих уравнений теории упругости по толщине оболочки и асимптотический метод. Все эти методы интенсивно развиваются, дополняя друг друга.

Список обозначений


a1, a2 - криволинейные ортогональные координаты срединной поверхности So оболочки на линиях главных кривизн; для оболочки вращения a1 ─ продольная, a2-окружная координаты; z ─ координата по нормали Теория оболочек к S;

А1, А2 -коэффициенты Лямэ; k1, k2 -главные кривизны;

U, V, W - компоненты вектора перемещений произвольной точки оболочки;

u, v, w - компоненты вектора перемещений точек поверхности So;

q 1, q2 - углы поворота нормали Теория оболочек;

ejk - компоненты тензора деформаций;

E11, E22, E12 - компоненты тангенциальной деформации на S: растяжения-сжатия по направлениям координат a1 и a2 и сдвиг;

K11, K22, K12 - компоненты изгибной деформации: изменения главных кривизн и кручение;

T11, T22, S - тангенциальные внутренние усилия, приведенные к So: усилия растяжения-сжатия и сдвига;

M11, M22, H - изгибающие и крутящий моменты;

Q11, Q22 - перерезывающие силы;

q1, q2, q3 - компоненты внешней поверхностной нагрузки, приведенные к S;

E, n - модуль Юнга и коэффициенты Пуассона материала оболочки;

yj -унифицированные обозначения основных независимых переменных в разрешающих системах обыкновенных дифференциальных уравнений (ОДУ);

fj - операторы правых частей канонических систем ОДУ;

Рассмотрим элемент произвольной тонкой оболочки, пусть в дальнейшем

h - толщина оболочки, принимаемая в дальнейшем постоянной.

Обозначим через R1, R2- главные радиусы кривизны срединной поверхности оболочки S. R=min {R1, R2}.

Основным геометрическим параметром оболочки является параметр тонкостенности или относительная толщина, определяемый отношением e=h/R.

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

Будем считать оболочку тонкой, если ее относительная толщина значительно меньше единицы. Обычно оболочки считают тонкими при значении e<1/20. Значения 1/20 < e < 1/10 соответствуют оболочке средней толщины, а e > 1/10 - толстой оболочке.

Для незамкнутых оболочек можно задать характерный размер размер a. Тогда параметр тонкостенности можно определить как e = min (h/a, h/R).

Поверхность оболочки S, равноотстоящая от лицевых поверхностей S+ и S - называется ее срединной поверхностью.


Криволинейные, ортогональные системы координат


Правило дифференцирования базисных векторов криволинейной ортогональной системы координат определяется следующим образом:


es,t = - (Ht,s /Hs) et - dstСHt

С = em (…),m / Hm


Здесь Hm - параметры Ляме координатной системы, имеющие вид


Теория оболочек= (r, i) 2; Hi = Ѕr, iЅ.


Здесь r, I - радиус - вектор произвольной точки тела оболочки. В частности:


e1,1 = (H1,1/H1) e1 - (H1,1/H1) e1 - (H1,2/H2) e2 - (H1,3/H3) e3

e1,2 = (H2,1/H1) e2; e3,2 = (H2,3/H3) e2; Hi (a1, a2, a3)


Запишем условие совместности, которое в принятых обозначениях имеет вид:


(e1,1),2 = (e1,2),1

(e1,2),1 = ( (H2,1/H1) e2),1 = (H2,1/H1),1 e2 + (H2,1/H1) (H1,2/H2) e1;

(e1,1),2 = - [ (H1,2/H2) e2 + (H1,3/H3) e3],2 =

= - (H1,2/H2),2 e2 + (H1,2/H2) ( (H2,1/H1) e1+ (H2,3/H3) e3) -

(H1,3/H3),2 e3 - (H1,3/H3) (H2,3/H3) e2


Тогда, приравнивая коэффициенты при базисных векторах, получим:


e1: (H2,1H1,2) / (H1H2) - (H2,1H1,2) / (H1H2) є 0 - тождество

e2: (H2,1/H1),1 + (H1,2/H2),2 + (H1,3 Ч H2,3) /Теория оболочек = 0

e3: (H1,2 Ч H2,3) / (H2H3) - (H1,3/H3),2 = 0


Круговая перестановка индексов приводит к шести уравнениям совместности параметров Ляме.


Некоторые сведения из теории поверхностей


Рассмотрим произвольную гладкую поверхность и систему декартовых координат x, y, z.

Пусть r = r (a1, a2) - радиус-вектор произвольной точки срединной поверхности оболочки. Рассмотрим производные r по переменным a1 и a2


r,1 = r1; r,2 = r2


Введем в рассмотрение базис


r1 /Ѕr1Ѕ= e1 r2 /Ѕr2Ѕ= e2


и обозначим ЅraЅ = Aa на срединной поверхности S ( a3 =0). В этом случае ri = Ai ei

Составим скалярные произведения:


ra Ч rb =Gab; G11 =Теория оболочек; G22 = Теория оболочек


G12 = G21 = 0 для ортогональной системы координат

При этом образуется тензор второго ранга Теория оболочек= Gab ra rb, который называется первым фундаментальным тензором поверхности.


ds2 = (dr) 2 = (r,1 da1 + r,2 da2) 2 =

= (r1da1 + r2 da2) 2 = G11dТеория оболочек + 2G12 da1da2 +C22 dТеория оболочек =

= Теория оболочекdТеория оболочек + Теория оболочекdТеория оболочек;


Коэффициенты А1 и А2 являются коэффициентами первой квадратичной формы и называются параметрами Ляме. Первая квадратичная форма определяет так называемую внутреннюю геометрию поверхности и определяет метрику поверхности. Введем в рассмотрение единичный вектор внешней нормали к поверхности N. Запишем очевидное соотношение N Ч N =1 и продифференцируем его по a1, a2:

2N Ч N, i = 0; очевидно, вектор N, I лежит в касательной плоскости к поверхности S и может быть представлен в виде разложения N, i = Bij rj.

При этом вводится в рассмотрение тензор второго ранга


Теория оболочек= Bab Ч ra Ч rb,


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

В главных осях тензор Теория оболочекможет быть записан в виде:


Теория оболочек= Теория оболочек= k1e1e1 + k2e2e2

k1 = 1/R1; k2 = 1/R2 –


главные кривизны

В дальнейшем координатные линии выбираются вдоль главных осей кривизны. Пусть в дальнейшем


I1 = k1 + k2 - первый инвариант (средняя кривизна) Теория оболочек

I2 = k1 Ч k2 - второй инвариант (гауссова кривизна) Теория оболочек


Специальная система координат в теории оболочек


N = e1 ґ e2


Для любой точки тела оболочки:


r (a1,a2,a3) = r (a1,a2) + a3N

Теория оболочек = (r, i) 2 = (r, i + (a3N), i) 2 = (ri + a3Bij rj) 2 (B12 = B21 =0)

Теория оболочек= (r1 + a3 N,1) 2 = (r1 + a3 Ч B11 r1) 2 = Теория оболочек (1 + a3k1) 2

H1 = A1 (1 + a3k1); H2 = A2 (1 + a3k2); (ЅriЅ= Ai)

Теория оболочек = N Ч N = 1 ® H3 = 1 –


параметры Ляме в специальной системе координат


Соотношения Гаусса и Кодацци


Уравнения совместности параметров Ляме:


(H2,1/H1),1 + (H1,2/H2),2 + (H1,3Ч H2,3) /Теория оболочек= 0

(H1,2 Ч H2,3) / (H2H3) - (H1,3/H3),2 = 0


В специальной системе координат


Hb = Ab (1 + a3kb); H3 = 1 (b = 1,2)


Рассмотрим срединную поверхность a3 = 0


(A2,1/A1),1 + (A1,2/A2),2 + k1A1 k2A2 = 0 –


соотношение Гаусса.


A1,2 k2 - (A1k1),2 = 0, (A1k1),2 = A1,2 k2


при замене индексов получаем два соотношения Кодацци


(A2k2),1 = A2,1 k1


Вектор перемещений


u = R - r = u1e1 + u2e2 + u3e3


R - текущая конфигурация

r - отсчетная конфигурация


u, i = (ukek), i = (uk), i ek + uk (ek), i


Дифференцирование ортов в специальной системе координат


e1,1 = - e2 (H1,2/H2) - e3 (H1,3/H3) = - e2 1/ (A2 (1 + a3k2)) Ч [A1 (1 + a3k1)],2 -

e3 Ч [A1 (1 + a3k1)],3 = - e2 1/ (A2 (1 + a3k2)) [A1,2 + a3 (A1k1),2] - e3 A1k1 =

= - e2 (A1,2 (1 + a3k2)) / (A2 (1 + a3k2)) - e3 A1k1 =

= - e2 (A1,2/A2) - e3 A1k1;

e1,2 = e2 (H2,1/H1) = e2 1/ (A1 (1 + a3k1)) [A2,1 + a3 (A2k2),1] =

= e2 (A2,1 (1 + a3k1)) / (A1 (1 + a3k1)) = e2 (A2,1/A1);

e1,3 = e3 (H3,1/H1) = 0 (т.к H3 = 1)

e2,1 = e1 (A1,2/A2) - получаем из e1,2 заменой (1«2)

e2,3 = e3 (H3,2/H2) = 0 e3,2 = e2 (H2,3/H3) = e2 A2k2

e3,1 = e1 (H1,3/H3) = e1 A1k1 e3,3 = 0 (H3 = 1)


Удлинения, сдвиги и повороты элемента сплошной среды


а) Рассмотрим удлинения

dr - в отсчетной конфигурации, dR - в текущей конфигурации


dR = dr Ч Теория оболочек; Теория оболочекR (Теория оболочек= ek (…),k / Hk)

R = r + u

Теория оболочек (r + u) = Теория оболочекr +Теория оболочекu = Теория оболочекu


Рассмотрим относительное удлинение


(ЅdRЅ-ЅdrЅ) /ЅdrЅ = e; ЅdRЅ = dS; ЅdrЅ = ds;

dS2 - ds2 = dR Ч dR - dr Ч dr = dr ЧТеория оболочекЧ dr Ч Теория оболочек- dr Ч dr =

= dr ЧТеория оболочекЧТеория оболочекЧ dr - dr ЧТеория оболочекЧ dr = dr (Теория оболочекЧТеория оболочек-ЧТеория оболочек) Ч dr =

= 2dr Ч Теория оболочекЧ dr; Теория оболочек= 0,5 (Теория оболочекЧТеория оболочек - Теория оболочек) - тензор деформаций Грина

Теория оболочек = 0,5 [ (Теория оболочек +Теория оболочекu) (Теория оболочек +Теория оболочекuT) - Теория оболочек] = 0,5 (Теория оболочекu +Теория оболочекuT +Теория оболочекuЧТеория оболочекuT)

dr = e ds ® e = dr /ЅdrЅ - единичный вектор

dS2 - ds2 = 2ds2 eЧeG Ч e

(dS2 - ds2) /ds2= (dS/ds) 2 - 1 = 2eЧeGЧe

dS/ds = (1 + 2eЧeGЧe) 1/2;

ee = (dS - ds) /ds = (1 + 2eЧeGЧe) 1/2 - 1 - удлинение

Пусть e = e1; Теория оболочек= (1+2Теория оболочек) 1/2 - 1 = 1 + Теория оболочек+ … - 1 = Теория оболочек » e11

e = 0,5 (СuuT) - линейный тензор деформаций Коши.


Деформации сдвига


Выделим два прямолинейных волокна, направление которых определяется единичными векторами m1 и m2


dr1 = m1ds1; dr2 = m2ds2;


dsi = ЅdriЅ - длины элементов волокон до деформаций

Деформации сдвига характеризуется изменением угла q12


cos q12 - cos Q12 = (dr1Ч dr2) / (ds1Ч ds2) - (dR1Ч dR2) / (dS1Ч dS2) =

= m1 Ч m2 - [ (drТеория оболочекЧ dr2) / ds1 (1+em1) ds2 (1+em2)] =

= m1 Ч m2 - mТеория оболочекЧ m2 = m1Ч (Теория оболочек - Теория оболочек) Ч m2 = - 2mТеория оболочекЧ m2;

Пусть m1 = e1; m2 = e2; m1 Ч m2 = 0

cos q12 = 0 0 - cos Q12 = -2Теория оболочек

cos Q12 = cos (p/2 - g12) = 2Теория оболочек = sin g12 = g12


g12 - угол сдвига; g12 » e12, если g12 - небольшой


Повороты


Рассмотрим материальное волокно dr = e ds


w = (dr ґ dR) / (ЅdrЅЧЅdRЅ) - вектор поворота материального волокна

ЅwЅ= sin j

w - нормаль, относительно которой происходит поворот

w = (dr ґ (drЧТеория оболочек)) / (dsЧds (1 + ee)) = e ґ (eЧТеория оболочек) =

= e ґ [e Ч (Теория оболочекu)] = e ґ e + e ґ (e Ч Теория оболочекu) = e ґ (e ЧТеория оболочекu)


Пусть e = et - базисные вектора t = 1,2,3

wt - вектор поворота материального волокна t


wt = et ґ (et Ч Теория оболочекu) = et ґ (et Ч ukj ek ej) =

Теория оболочекu = ukj ek ej Теория оболочек= et ґ (ukj dtk ej) = et ґ utj ej =

= utj 'tjk ek = wtk ek = wt, где wtk = utj 'tjk


'tjk - символы Леви-Чивита, которые определяются:

СЛЧ = 0, если среди r,s,t есть одинаковые

=+1, если индексы r,s,t - различные ® 123, 231, 312

= -1, если этот порядок нарушается

'rst = er Ч (es ґ et)

wtk характеризует поворот орта t относительно орта k.

Введем тензор второго ранга Теория оболочек= wtk et ek - тензор поворота


w11 = 0; w12 = u1j '1j2 = - u13; w13 = u1j '1j3 = u12

w21 = u2j '2j1 = u23; w22 = 0; w23 = - u21

w31 = u3j '3j1 = - u32; w32 = u3j '3j2 = u31;

w33 = 0;


Определим компоненты градиента вектора перемещений в специальной системе координат:


(Теория оболочек= es (…),s / Hs; Hi = Ai (1 + a3ki) i = 1,2 H3 = 1


“o” - в дальнейшем опускаем


Сu = es (u1e1 + u2e2 + u3e3),s / Hs = e1 (u1,1/H1) e1 + e1 (e1,1 /H1) u1 + e2 (u1,2/H2) e1 +

+ e2 (e1,2 /H2 ) u1 + e3 (u1,3/H3) e1 + e3 (e1,3 /H3) u1 + e1 (u2,1/H1) e2 + e1 (e2,1 /H1) u2 +

+ e2 (u2,2/H2) e2 + e2 (e2,2 /H2) u2 + e3 (u2,3/H3) e2 + e3 (e2,3 /H3) u2 + e1 (u3,1/H1) e3 +

+ e1 (e3,1 /H1) u3 + e2 (u3,2/H2) e3 + e2 (e3,2 /H2) u3 + e3 (u3,3/H3) e3


После подстановки выражений ek,j (j,k = 1,2,3)


Hi = Ai (1 + a3ki) i = 1,2; H3 = 1


h/2 Ј a3 Ј h/2; учитывая, что h/Ri " 1, т.е. оболочка тонкая, получим:


u11 = u1,1/A1 + (A1,2 / (A1A2)) u2 + u3k1

u12 = u2,1/A1 - (A1,2 / (A1A2)) u1

u13 = u3,1/A1 - u1k1

u21 = u1,2/A2 - (A2,1 / (A1A2)) u2

u22 = u2,2/A2 + (A2,1 / (A1A2)) u1 + u3k2

u23 = u3,2/A2 - u2k2

u31 = u1,3 u32 = u2,3 u33 = u3,3

Теория оболочек = 0,5 (Сu + СuT) Ю eii = uii


Для удлинений имеем:


e11 = u11; e22 = u22; e33 = u33;


Для деформаций сдвига соответственно:


e12 = 0,5 (u12 + u21) = 0,5 [ (A2/A1) (u2/A2),1 + (A1/A2) (u1/A1),2]

e13 = 0,5 (u13 + u31) = 0,5 (u3,1/A1 + u1,3 - u1k1)

e23 = 0,5 (u23 + u32) = 0,5 (u3,2/A2 + u2,3 - u2k2)


Углы поворота определяются через перемещения следующим образом: wii = 0


w12 = - u13 = - u3,1/A1 + u1k1

w21 = u23 = u3,2/A2 - u2k2; w13 = u12

w23 = - u21 w31 = - u32 = - u2,3; w32 = u31 = u1,3.


Теория малых удлинений и малых квадратов углов поворота


Рассмотрим тензор нелинейных деформаций Грина:


Теория оболочек = 0,5 (Сu + СuT + Сu Ч СuT) = Теория оболочек + 0,5Сu Ч СuT


Его нелинейная часть определяется следующим образом:


Сu Ч СuT = umn em Ч en Ч uij ej ei = umn uij dnj eei =

= umn uin em ei; Ю emi = emi + 0,5umn uin


Таким образом, компоненты тензора деформаций можно записать в виде:


e11 = e11 + 0,5 (Теория оболочек)

e12 = e12 + 0,5 (u11u21 + u12u22 + u13u23)

e13 = e13 + 0,5 (u11u31 + u12u32 + u13u33)

e21 = e21 + 0,5 (u21u11 + u22u12 + u23u13)

e22 = e22 + 0,5 (Теория оболочек)

e31 = e31 + 0,5 (u31u11 + u32u12 + u33u13)

e32 = e32 + 0,5 (u31u21 + u32u22 + u33u23)

e33 = e33 + 0,5 (Теория оболочек)

e23 = e23 + 0,5 (u21u31 + u22u32 + u23u33)


или, подставляя выражения для углов поворота:


e11 = e11 + 0,5 (Теория оболочек)

e22 = e22 + 0,5 (Теория оболочек)

e33 = e33 + 0,5 (Теория оболочек)

e12= e12 + 0,5 (-e11w23 + e22w13 - w12w21)

e13= e13 + 0,5 (e11w32 - w13w31 - w12w33)

e23= e23 + 0,5 (- w32w23 -e22w31 + e33w21)

(u21 = -w23; u23 = w21; u31 = w32; u12 = w13; u32 = -w31;

u31 = w32; u11 = e11; u22 = e22; u33 = e33)


Введем следующие предположения:

eii << 1 - деформации растяжения -сжатия малы

предполагаем, что величины поворотов w13 << 1; w23 << 1, а в отношении остальных величин можно принять, что Теория оболочек<< 1

wij - угол поворота i-го орта относительно j-го орта. Таким образом из соотношений (…) следует:


e11 = e11 + 0,5Теория оболочек

e22 = e22 + 0,5Теория оболочек

e33 = e33 + 0,5 (Теория оболочек)

e12= e12 - 0,5w12 w21

e13= e13 + 0,5 (e11w32 - w13w31 - w12w33)

e23= e23 + 0,5 (- w32w23 -e22w31 + e33w21)


Гипотезы Кирхгофа-Лява


Результаты, полученные в предыдущих параграфах, основаны на геометрических и статических соображениях. Однако их недостаточно для полного построения теории оболочек. При выборе соотношений, связывающих компоненты деформаций с перемещениями и усилия и моменты с компонентами деформаций приходится принимать некоторые упрощающие подходы. Первый заключается в том, что оболочку рассматривают как трехмерное упругое тело. Решение соответствующих уравнений теории упругости разыскиваются путем разложения всех величин в ряды по степеням точки оболочки от срединной поверхности. Этот подход, предложенный в теории пластин А. Коши, позволяет при удержании достаточного числа членов (при условии сходимости рядов) получить решение близкое к точному. Этот метод весьма громоздок, поэтому в большинстве случаев идут по другому пути. Во втором подходе предложенном также при построении теории пластин Г. Кирхгоффом принимаются гипотезы, аналогичные тем, которые используются в теории балок:

прямолинейные и нормальные к срединной поверхности волокна недеформированной оболочки остаются прямолинейными и нормальными к деформированной срединной поверхности и не меняют своей длины;

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

Первая гипотеза имеет геометрический характер, вторая - статический. Теория оболочек, основанная на гипотезах Кирхгоффа, была построена, в основном А. Лявом, поэтому в теории оболочек гипотезы 1 и 2 принято называть гипотезами Кирхгоффа-Лява. Иногда их называют гипотезой жесткой (недеформированной) нормали или гипотезой сохранения нормали.

Гипотеза 1 используется только для записи зависимостей деформации оболочки от перемещений, гипотеза 2 - для записи зависимостей деформаций от напряжений. В первом случае предполагается, что в нормальных сечениях отсутствуют сдвиги e13 = e23 = 0, и поперечные деформации e33 = 0. Во втором случае допускается, что нормальное напряжение s33 незначительно влияет на деформации e11 и e22, так что эти деформации выражаются через нормальные напряжения s11, s22 и s33 << {s11, s12, s22}. Таким образом гипотезу 1 нельзя понимать в буквальном смысле, поскольку в действительности в оболочках имеют место поперечные сдвиги - поперечные или как их иногда называют перерезывающие силы не равны нулю.

Гипотезы Кирхгоффа-Лява просты и физичны. Они позволяют свести трехмерную задачу определения напряженно-деформированного состояния оболочки к двумерной. Исследование поведения элемента оболочки в рамках этих гипотез сводится к исследованию поведения ее срединной поверхности. Следует отметить, что теория, построенная на гипотезах Кирхгоффа-Лява, является существенно приближенной. Принятие этих гипотез вносит погрешность порядка h/R, где h - толщина оболочки, R - минимальный линейный размер срединной поверхности.

Рассмотрим элемент тонкой оболочки со срединной поверхностью S. До деформирования в исходной конфигурации радиус-вектор произвольной точки оболочки, не лежащей на срединной поверхности может быть представлен в виде:


r (a1, a2, a3) = r (a1, a2) + a3n


r (a1, a2) - радиус-вектор проекции точки на S до деформации.

После деформирования (в актуальной конфигурации)


R (a1, a2, a3) = P (a1, a2) + a3N


P (a1, a2) - радиус-вектор проекции точки на S после деформации

Тогда вектор перемещений запишется в виде:


u = R - r = P - r + a3 (N - n) = = u° (a1, a2) + a3u1 (a1, a2)


u° - вектор перемещений точек, лежащих на S

В координатной форме соответственно:


u1 = Теория оболочек (a1, a2) + a3Теория оболочек (a1, a2)

u2 = Теория оболочек (a1, a2) + a3Теория оболочек (a1, a2)

u3 = Теория оболочек (a1, a2)


При этом компоненты вектора перемещений u1 и u2 линейным образом зависят от координаты a3, а функция поперечного прогиба постоянна по толщине в силу недеформируемости нормали.

Рассмотрим детальнее геометрическую гипотезу Кирхгоффа-Лява. Тот факт, что нормаль к срединной поверхности S в процессе деформирования остается нормалью приводит к соотношениям:


e13 = 0, e23 = 0


Таким образом


u3,1/A1 + u1,3 - u1k1 = 0

Теория оболочек/A1 + Теория оболочек - (Теория оболочек + a3Теория оболочек) k1 = 0

Теория оболочек (1 - a3k1) - Теория оболочекk1 + Теория оболочек/A1 = 0


считая оболочку достаточно тонкой, пренебрегаем членом a3k1 << 1


Теория оболочек (a1, a2) = -Теория оболочек/A1 + Теория оболочекk1

Теория оболочек Теория оболочек (a1, a2) = -Теория оболочек/A2 + Теория оболочекk2

w12 = - u3,1/A1 + u1k1 = -Теория оболочек/A1 + Теория оболочекk1 + a3k1Теория оболочек=

= - Теория оболочек/A1 + Теория оболочекk1 + a3k1 (-Теория оболочек/A1 + Теория оболочекk1) =

= (-Теория оболочек/A1 + Теория оболочекk1) (1 + a3k1) = Теория оболочек (1 + a3k1) » Теория оболочек=

= -Теория оболочек/A1 + Теория оболочекk1 = q1 (a1, a2) - угол поворота на поверхности S.


Аналогично:


w21 = u3,2/A2 - u2k2 » Теория оболочек= Теория оболочек/A2 - Теория оболочекk2 = - Теория оболочек = - q2 (a1, a2)


Обозначим: Теория оболочек = u; Теория оболочек = v; Теория оболочек= w тогда можно записать:


u1 = u + a3q1, u2 = v + a3q2, u3 = w


Введем в рассмотрение плоский вектор перемещений и поворотов:


u = ue1 + ve2

q = q1e1 + q2e2


Тензор кривизны в главных осях можно представить в виде:


Теория оболочек= k1e1e1 + k2e2e2; ki = 1/Ri


Окончательно кинематические соотношения, соответствующие теории Кирхгоффа-Лява запишутся в виде:


q = -Сw + Теория оболочекЧ u С = es (1/As) (¶/¶s) (s = 1,2)

u (a1, a2, a3) = u1e1 + u2e2 = u (a1, a2) + a3q (a1, a2)

u3 = w (a1, a2)


С учетом проведенных выкладок для компонентов тензора деформаций имеем:


e11 = e11 + 0,5Теория оболочек = Теория оболочек + a3Теория оболочек

Теория оболочек= Теория оболочек/A1 + A1,1 / (A1A2) Теория оболочек + Теория оболочекk1 + 0,5Теория оболочек

Теория оболочек= q1,1/A1 + A1,1 / (A1A2) q2 (w12 = q1)

e22 = e22 + 0,5Теория оболочек

Теория оболочек= Теория оболочек/A2 + A2,1 / (A1A2) Теория оболочек + Теория оболочекk2 + 0,5Теория оболочек

Теория оболочек= q2,2/A2 + A2,1 / (A1A2) q1

e12= e12 - 0,5w12w21 = e12 + 0,5q1q2 = Теория оболочек (w21 = -q2)

Теория оболочек = 0,5 [ (A2/A1) (Теория оболочек/A2),1 + (A1/A2) (Теория оболочек/A1),2] + 0,5q1q2

Теория оболочек = 0,5 [ (A2/A1) (q2 /A2),1 + (A1/A2) (q1/A1),2]


Введем в рассмотрение плоский тензор деформаций


Теория оболочек= e11e1e1 + e12 (e1e2 + e2e1) + e22e2e2


Он может быть записан в другой форме:


Теория оболочек = Теория оболочек+ a3Теория оболочек, где

Теория оболочек = Теория оболочекe1e1 + Теория оболочек (e1e2 + e2e1) + Теория оболочекe2e2 (b = 0,1)


Таким образом использование геометрической гипотезы Кирхгоффа-Лява приводит к линейному распределению перемещений и деформаций по толщине оболочки. В компактной форме можно записать:


Теория оболочек= 0,5 (Сu+ СuT) + Теория оболочекw + 0,5qq

Теория оболочек= 0,5 (Сq + СqT)


Теория оболочек характеризует деформации растяжения-сжатия срединной поверхности S, Теория оболочек - изменение кривизн и кручение срединной поверхности.

Определение напряженного состояния оболочки


dS = H2 da2da3 = A2 (1 + a3k2) da2da3

S11 = Теория оболочек A2 (1 + a3k2) da3da2= A2da2Ч

Ч Теория оболочек (1 + a3k2) da3, A2da2 –


длина средней линии.

Введем усилия на единицу длины:


T11 = S11/ (A2da2) Теория оболочек (1 + a3k2) da3


Аналогично:


T12 = Теория оболочек (1 + a3k2) da3;

T22 = Теория оболочек (1 + a3k1) da3;

T21 = Теория оболочек (1 + a3k1) da3


Tab - усилия растяжения-сжатия в срединной поверхности оболочки. В дальнейшем:


Tab = Теория оболочекda3 (a, b) = 1,2, Теория оболочек= Tabeaeb = T11e1e1 + T22e2e2 + T12 (e1e2 + e2e1)


Введем изгибающие моменты


Мab = Теория оболочекa3da3

Теория оболочек= M11e1e1 + M22e2e2 + M12 (e1e2 + e2e1)

Qb = Теория оболочекda3, b = 1,2 sb - перерезывающие силы

Q = Q1e1 + Q2e2


Свяжем напряженное состояние с ее деформир.:

Замечание о возможности использования линейных физич. соотношений

Материал: однородный, изотропный

Обобщенный закон Гука:


e11 = (1/E) [s11 - n (s22 + s33)] e12 = (1/2m) s12

e22 = (1/E) [s22 - n (s11 + s33)] e23 = (1/2m) s23

e33 = (1/E) [s33 - n (s11 + s22)] e13 = (1/2m) s13

E = 2m (1 + n)


Теория - геометр. нелинейная, но физич. - линейная.

По 2-ой гипотезе Кирхгофа-Лява:


s33 = 0

s11 - ns22 = Ee11 Ю s11 = E / (1-n2) (e11 + ne22)

s22 - ns11 = Ee22 s22 = E / (1-n2) (e22 + ne11)

s12 = 2m e12 E / (1 + n) e12

T11 = Теория оболочекda3 = Теория оболочек/ (1-n2) (e11 + ne22) da3 =

= E / (1-n2) Теория оболочек

= (Eh) / (1-n2) Теория оболочек - т.к материал - однород.

B = (Eh) / (1-n2) - жесткость на растяж. - cжатие

T11 = BТеория оболочек; T22 = BТеория оболочек;

T12 = (Eh) (1+n) Теория оболочек = B (1-n) Теория оболочек

M11 = Теория оболочек=

= E / (1-n2) Теория оболочек

= (E) / (1-n2) Теория оболочекТеория оболочек=

= (Eh3) / (12 (1-n2)) Теория оболочек

M11 = DТеория оболочек, D = (Eh3) / (12 (1-n2)) - жесткость на изгиб (цилиндрическая жесткость)

M22 = DТеория оболочек; M12 = D (1-n) Теория оболочек.


Рефетека ру refoteka@gmail.com