Обоснование принципа Сен-Венана для естественно-закрученного стержня Текст научной статьи по специальности «Математика»

Владикавказский математический журнал 2010, Том 12, Выпуск 1, С. 53-67

УДК 539.3

ОБОСНОВАНИЕ ПРИНЦИПА СЕН-ВЕНАНА ДЛЯ ЕСТЕСТВЕННО-ЗАКРУЧЕННОГО СТЕРЖНЯ1

Ю. А. Устинов

В работе дается математическое обоснование принципа Сен-Венана для естественно-закрученного

стержня.

Ключевые слова: принцип Сен-Венана, спектральная задача, собственные значения, элементарные решения, обобщенная ортогональность.

Введение

В 1856 г. были опубликованы два знаменитых мемуара Сен-Венана [1, 2], заложившие основы прикладной теории упругости (широкому кругу читателей на русском языке они стали доступны благодаря переводу Г. Ю. Джанелидзе [3]). Результаты исследований, изложенные в этих трудах, под названием «задачи Сен-Венана» вошли практически во все учебники по теории упругости, начиная с изданной в 1862 г. книги Клебша «Теория упругости твердых тел» и переизданной в 1883 г. в переводе Сен-Венана [4] с громадным количеством комментариев, едва не превышающих основной объем. Клебшем же и был введен термин «задачи Сен-Венана», который включает в себя задачи о растяжении, кручении, чистом изгибе и изгибе перерезывающей силой цилиндра (призмы) с произвольным основанием, боковая поверхность которого свободна от напряжений. Следует отметить, что первая попытка исследовать задачу кручения стержня с прямоугольным поперечным сечением на основе уравнений теории упругости принадлежит Коши. Из построенного им приближенного решения вытекает, что поперечные сечения не остаются плоскими. Этот результат был использован Сен-Венаном для построения более полной теории кручения и изгиба призматических стержней. В основе этой теории лежит приближенный полуобратный метод построения решений трехмерных уравнений теории упругости и «принцип упругой равнозначности статически эквивалентных систем сил». Более общая формулировка этого принципа, названного «принципом Сен-Венана», была дана учеником Сен-Венана Буссинеском. Как показывает дальнейшая история развития теории упругости, полуобратный метод и принцип Сен-Венана стали важным инструментом исследования напряженно-деформированного состояния упругих тел.

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

© 2010 Устинов Ю. А.

1 Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований, проект №09-01-00065-а.

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

Математическое обоснование принципа Сен-Венана для призмы была дано в работах [5, 6, 7], в [8] на конкретном примере показан условный характер этого принципа.

В настоящей работе излагается метод построения решений Сен-Венана для естественно закрученного стержня и приводится строгое математическое обоснование принципа Сен-Венана.

1. Геометрическое описание естественно закрученного стержня

Под естественно закрученным стержнем (ЕЗС) принято понимать трехмерное тело, которое получается в результате винтового движения вдоль оси некоторой плоской фигуры. Геометрию такого тела можно описать следующим образом.

Пусть Ж1,Ж2,Жз = x — декартова система координат с ортами ¿1, ¿2, ¿3. Объем, занимаемый ЕЗС, получим в результате винтового движения плоской фигуры S параллельно плоскости x1 x2 вдоль оси x.

В качестве параметра винта возьмем относительный угол закручивания т, который принято называть «круткой». Таким образом, угол поворота сечения, находящегося на расстоянии x по оси от начала координат, р = тх.

Введем также сопутствующую (винтовую) систему координат £ь£2,£3 = £ так, что оси £1 и £2 жестко связаны с фигурой S в процессе ее движения. Дополнительно предположим, что центр тяжести S расположен на расстоянии Го от оси £, а оси £а (а = 1, 2) направлены по главным осям инерции. Связь между координатами обеих систем определяется соотношениями

£1 = x1 cos р + х2 sin р, £2 = —x1 sin р + х2 cos р, £ = х3. (1.1)

Орты, определяющие направление осей £а, обозначим ea (а = 1, 2). Таким образом, R = £1e1 + £2e2 + £¿3 — радиус-вектор произвольной точки области V.

Обозначим через Г боковую поверхность ЕЗС, через N — ее внешнюю нормаль. Определим проекции N на оси сопутствующей системы координат. Уравнение контура dS области S зададим соотношениями £ао = £«o(s), где £ао — координаты точки контура; s — расстояние по дуге контура от некоторой фиксированной точки.

Следующие соотношения определяют касательный t и нормальный n орты контура dS:

t = £'a0ea, n = na ea = £2oe1 — £10e2, (1.2)

где £10, £20 — производная £10, £20 по s.

Рассмотрим два сечения ЕЗС, расположенные на расстоянии dx. Точки сечений A и B будем называть эквивалентными, если им соответствуют одни и те же координаты £1, £2. Очевидно, что радиусы-векторы Ra, Rb двух бесконечно близких эквивалентных точек, лежащих на боковой поверхности Г, связаны соотношением

Rb = Ra + [¿з + т(¿3 х Г0)] dx, (1 3)

Г0 = £10 e1 + £20e2.

Очевидно также, что вектор

N = г х (Яв — Яа )

(1.4)

ортогонален поверхности Г.

Подставляя (1.2), (1.3) в (1.4), после преобразований получаем следующие выражения для компонент N в базисе винтовой системы координат:

N1 = П1 = £2о, N = П2 = —£10, N3 = тЬп,

й /г2"

6™ = - П1 £2 = у) .

(1.5)

2. Операторные формы уравнений равновесия теории упругости

Будем считать, что материал ЕЗС — изотропный, упругие свойства которого определяются двумя упругими постоянными V — модулем сдвига и коэффициентом Пуассона соответственно; 0 ^ х ^ I, где I — длина ЕЗС.

Обозначим через д расстояние между двумя произвольными точками Р1, Р2 £ дБ; Н = 8ир(д) — диаметр области; через и = (и1,и2,из)т — безразмерный (отнесенный к Н) вектор смещений. Введенные выше координаты , £к (к = 1, 2, 3) и «крутку» т, не меняя обозначений, будем считать безразмерными (и& ~ и&/Н, т ~ Нт). В дальнейшем компоненты ик рассматриваются как функции £1, £2, £.

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

(и, V) = J (и1 «1 + и2^2 + изЬз) й£1 й£2.

(2.1)

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

Ьт (д)и = д2С и + дВт и + Ат и = 0,

(2.2)

где

С=

00 10 0 2(1 + к)

Вт = Взт — В'3

3т,

тВ -т д1

В3т = т тВ д2

2кд1 2кд2 2(1 + к)тВ,

тВ -т 2кд1

В3 т = - т тВ 2кд2

д1 д2 2(1 + к)тВ,

Ат--11 |

(2.3)

Ац = 2(1 + к)д2 + д2 + т2(В2 - 1), А12 = (1 + 2к)д^ - 2т2В,

А13 = т (1 + 2к)д1Д

А22 = д2 + 2(1 + к)д2 + т2(В2 - 1),

А31 = т (1 + 2к)Вд1,

А33 = д2 + д22 + 2(1 + к)т 2В2.

А21 = (1 + 2к)д1д2 + 2т2 Д А23 = т (1 + 2к)д2А А32 = т (1 + 2к)Вд2,

Здесь

д д V

да = , д = —, В = £2д1 - &д2, к = --—.

д£а д£ 1 - 2v

Будем считать, что боковая поверхность ЕЗС свободна от напряжений. Это условие можно представить в виде

Мт (д )и = (дОт + Ет )и

От =

тЬп 0

П1

0 2кщ тЬп 2кП2 П2 2(1 + к)тЬ п

0,

Ет = || Ег^-

(2.4)

Е12 = 2кщд2 + П2д1 - т2Ьп

Ец = 2(1 + к)щд1 + П2 д2 + П3т 2ЬпВ, Е13 = 2кп1тВ + тЬпд1, Е22 = П1д1 + 2(1 + к)п2 д2 + т 2ЬпВ, Е31 = П1тВ + П2т + 2ктЬп д1, Е33 = П1д1 + П2д2 + 2(1 + к)т 2ЬпВ. Объединяя уравнение равновесия (2.2) и граничное условие (2.4), получаем:

Е21 = П1д2 + 2кп2д1 + т 2Ьп,

Е23 = 2кп2тВ + тЬпд2,

Е32 = -П1т + N2 тВ + 2ктЬпд2,

¿1т (д)и = { А (д)и,Мт (д)и} = 0.

(2.5)

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

Введем расширенный шестикомпонентный вектор

w

= (и, а)Т = («1, П2, «3, а03, а03, ст3э)Т ,

(2.6)

где

а = Сди + В3т и. (2.7)

Заметим, что /а = (лст03,/033 )Т является вектором напряжений на площадках с нормалями, параллельными оси стержня (13). Из (2.7) выразим ди через и и а

ди = -С-1В3т и + С-1а.

1

(2.8)

Второе уравнение получается на основе векторного вида уравнения равновесия в напряжениях [1, 2]

да = -В3т С-1В3т и + Ат и + В3т С-1а. (2.9)

г

Объединив (2.8) и (2.9), получаем операторное дифференциальное уравнение первого порядка относительно т

дт — Тот = 0, То = (Та/3), (2.10)

Т11 = —С-1В3т, Т12 = —С-1, Т21 = —В3 тС—1 В3т + Ат, Т22 = В3 т С-1. Граничное условие (2.4) запишется в виде

2т = 0, 2 = ) (а, в = 1, 2), (2.11)

2ц = 212 = 0, 221 = —СтС 1В3т + Ет, 222 = СтС 1.

В гильбертовом пространстве Н1 = Н ф Н вектор-функций V = (а, 6)т со скалярным произведением

(V!, V2)l = (а1, а2) + (Ьь 62) (2.12)

на основании дифференциальной формы Т0 и граничного условия (2.11) определим оператор Т и запишем уравнения равновесия и граничные условия на боковой поверхности ЕЗС в следующей форме:

дт — Тт = 0, Тт = {Т0т, 2т} = 0. (2.13)

3. Однородные элементарные решения и их свойства

3.1. В теории упругости широко используется термин однородные решения.

Определение. Однородным решением будем называть любую вектор-функцию и, удовлетворяющую уравнению (2.5), а также любую вектор-функцию т(х), удовлетворяющую уравнению (2.13).

Иными словами в данном случае однородным является любое решение однородных уравнений теории упругости, удовлетворяющее условию отсутствия напряжений на боковой поверхности ЕЗС.

Будем искать решение уравнения (2.5) в виде

и(£) = е* а(£1,£2). (3.1)

Подставляя (3.1) в уравнение (2.5), получаем спектральную задачу на поперечном сечении ЕЗС

Т1(7,ш)а = 0. (3.2)

Некоторые особенности распределения множества собственных значений {7к} и общие свойства собственных и присоединенных векторов описаны в [9—11]. Ниже они приведены в соответствии с введенной системой обозначений.

Если 7к — простое собственное значение (СЗ) задачи (3.2), ак — соответствующий ему собственный вектор (СВ), то элементарным решением уравнений (2.5), (2.13), отвечающим СЗ 7к, называется вектор-функция вида

ик (£) = е7к« ак (£1,£2). (3.3)

Соответственно, элементарным решением уравнения (2.13) называется вектор-функция

тк (£) = е7к«Vk (£1 ,£2), (3.4)

а собственные значения и собственные векторы определяются при решении спектральной

задачи

Ти - 7" = 0.

Задачи (3.2) и (3.5) эквивалентны, при этом

"к = (ак, Ък)т, Ък = (7кС + Взг )ак.

(3.5)

(3.6)

Если 7к — N-кратное собственное значение задачи, а° к (^ = 1,..., п < N) — множество линейно независимых собственных векторов, отвечающих 7к, при этом каждому СВ соответствует своя система присоединенных векторов а^ к (£ = 1,... ), то множество элементарных решений, отвечающих кратному СЗ 7к, состоит их следующих вектор-функций:

п{ (ж) = е7к х

а-1

¿1 а° >к

+

(* - 1)1

а

1 к

+-----+ а:

г , к

* = 0,...,Рк.

(3.7)

Общее количество таких элементарных решений будет

рк + рк + ••• + рП = N

где N — алгебраическая кратность, а п — собственная кратность СЗ 7к.

Множество элементарных решений уравнения (2.13), отвечающих кратному собственному значению, имеют вид

w3k (ж) = е

НкХ

г г— 1

Ж О ^ О А

1 ж -¿1 и°,к +-

(* -1)1

<к +

+ V

г,к

(3.8)

Собственный и присоединенный векторы определяются решением цепочки задач

Ти°,к - 7ки°,к = 0, Т<к -

1,к,

= 1,

1

,Рк,

(3.9)

где

и!к = (а!к, Ъ!к)Т.

Компоненты векторов и° к выражаются через векторы а° к по-прежнему соотношени-

.7

ем (3.6), а присоединенные векторы к выражаются через систему векторов {а^}р=к 1 следующим формулами:

V, к = К к, (7кС + вз)а^, к + Са''_ 1 ,к)Т, £ = 1,... ^.

(3.10)

3.2. Система собственных и присоединенных векторов оператора Т не обладает обычной ортогональностью, поскольку оператор не является самосопряженным. Однако для этой системы справедливы некоторые соотношения обобщенной ортогональности. Для построения таких соотношений определим сопряженный оператор Т*. Имеем

(Т"1, "2)1 = ("1,Т*"2)1.

Учитывая, что операторы С-1, А — самосопряженные, получаем:

Т= {ТЯ'и = 0}, (3.11)

То = (Тав), Я' = (^) (а, в = 1, 2),

г

ж

1

1

и

ТУ ГТ1 гтл1 ГТ1 гтл1 ГТ1 гтл1 ГТ1

11 — -Т22, Т 12 — Т 21; Т 21 — Т 12, Т 22 — -Т11>

1 = 2 — 0, ^21 — —^22, ^22 — ^21;

где Тав, определены соотношениями (2.10), (2.11). Введем оператор (симплектическую единицу)

3 —

0 -I

I 0

— - /* — - /-1,

где I — единичная матрица 3 х 3.

Утверждение 3.1. Оператор Т является гамильтоновым, т. е.

/Т — -Т */ (3.12)

< С учетом конкретного вида оператора Т и самосопряженности операторов С-1, А, имеем:

(7Т«1, «2)1 — -(Т21в1, а2) - (Т22&1, 02) + (Тпа1, Ь2) + (Т12&1, Ь2)

— ((ВС-1Вз - А)01, а^ - (ВС-1Ьь02) - (С-1Взаь 62) + (С-1Ьь 62)

— (01, (В3С-1Вз - А)а^ - (61, С-1Вз02) - (01, ВС+ (61, С-1 62). С другой стороны,

/Т«1, «2) — («1, (/Т)*«2) — («1,Т— -(«1 ,Т*/«2). Следовательно, оператор

_Т * / —

В3С-1Вз - -В3С-1 -С-1Вз С-1

т. е. /Т — -Т*/, что и требовалось доказать. >

Так как оператор Т является гамильтоновым, то его собственные значения расположены симметрично в комплексной плоскости 7, т. е. каждому СЗ , расположенному в первом квадранте, соответствуют еще три СЗ — -7к, 7к, -7к — 7-к. Следовательно, оператор Т* будет иметь тот же самый спектр (прямое доказательство этого свойства спектра приведено в [11]).

Утверждение 3.2. Пусть 7к, 7т — произвольные СЗ оператора Т, «к, «т — соответствующие им собственные векторы. Тогда имеют место следующие соотношения обобщенной ортогональности:

(«к, /«т)1 — 0, 7к — -7т. (3.13)

< Рассмотрим уравнения

Т«к — Тк «к, Т«т — 7т«т.

С учетом (3.13) получаем

(Т«к, /«т)1 — («к, Т*/«т)1 — -(«к,/Т«т)1 — -(«к,7т/«т)1

— -7т(«к, /«т)1 — 7к(«к, /«т)1. Из последнего равенства вытекает соотношение (3.13). >

Утверждение 3.3. Пусть Yk , Ym — кратные СЗ, {vs,k}, {vi>m} (s = 0,...,pk, t = 0,..., pm) соответствующие им жордановы цепочки; тогда для любой пары корневых векторов имеют место следующие соотношения обобщенной ортогональности:

(Vs,k, J)l =0, Yk = -7m. (3.14)

< Рассмотрим различные скалярные произведения (Tvs,k,Jvt>m)i, (vs,k,JTvt>m)i. Преобразуем их, используя соотношения (3.9), (3.10), следующим образом:

(T, JVt m) 1 = Yk(Vs,k, JVt,m)l + (Vs-l,k, JVijm)l,

(vs,k, JTvt|m)l = Ym(Vs,k, JVt,m)l + (Vs,k, JVt-1,m)l = (-TVs,k, JVi;m)b Складывая первое равенство со вторым, получаем

(Yk + Ym)(Vs,k, J Vi;m)l = -(vs-l,k, J Vi;m )l - (Vs,k, J Vt-l,m)l • (3.15)

Пусть s = 1, t = 0; тогда из соотношений (3.13), (3.19) вытекает

(Vl)fe, JVojm)l = 0.

Далее по индукции получаем доказательство соотношения (3.14). > Утверждение 3.4. У гамильтонова оператора при условии, что

j ,Jv0,k )l =0, (3-16)

все комплексные и не равные нулю вещественные собственные значения простые.

< Обратимся к соотношениям (3.14), (3.19) и положим в них Ym = Yk, t = 0. После элементарных преобразований для элементов конкретной жордановой цепочки получаем:

(Yk+7k )(vm+l,k, j v0,fc )l- J, j j )l=0. (3.17)

Из (3.17) вытекает

Vl,k = (Yk + 7k)-lVo,k, V2,k = (Yk + 7k)-2Vo,k,...,

т. е. все присоединенные векторы пропорциональны собственному и, следовательно, корневое подпространство состоит из одного собственного вектора. >

Замечание. Условие (3.16) продиктовано тем, что метрика, порождаемая оператором J — индифинитная [11]. Во всех известных примерах из теории упругости и акустики это условие выполняется.

Рассмотрим отдельно частный случай, когда Yk — чисто мнимое или нулевое кратное

j

собственное значение и ему отвечают одна или несколько жордановых цепочек {vm k}m=o (j = 1,...,n). Как будет показано ниже, множество соответствующих ЭР формируют решения Сен-Венана естественно-закрученного стержня.

Из условий (3.17) вытекают следующие соотношения ортогональности элементов жордановой цепочки.

Утверждение 3.5. Пусть ys — кратное мнимое (или нулевое) собственное значение

j

оператора T и пусть { V s }p=o — одна из соответствующих ему жордановых цепочек; тогда справедливы следующие соотношения:

Ks,J<s)l =0, t + q < pS, (3.18)

(<т>в, /<-т,з)1 — (-1«, t + 9 — / (3.19)

Здесь — (« з , 5)1 — - (а5з , Ь° 5) + (Ьз , а° 5) — инвариант жордановой цепочки.

Ре • р8' р8'

4. Общее представление решения трехмерной задачи и обоснование принципа Сен-Венана

4.1. Рассмотрим вектор твердого смещения ЕЗС

и0 — а0 + ш х г, (4.1)

где а0 — (а^, а°, а°)т — вектор поступательного смещения; ш — (ш1,ш2,шз)т — вектор малого поворота; а°, ш/ — проекции на оси основной системы координат ж/. Проектируя (4.1) на оси подвижной системы координат, получаем:

и°1 — С^ + С2е-гт« + Сз£егт« + С4£в-гт« - Себ, °2 — ¿С^ - ¿С2егт« + ¿Сз£егт« - ¿С4£е-гт« + С6£1, (4.2)

где

и° — -Сз(е^ - С4Се-г^ + С5,

С — £1 + ¿£2, С — £1 - ¿£2,

С1 — 2 (а? - ¿а°), С2 — С1, Сз — 2(ш2 + ¿Ш1), С4 — Сз, С5 — а°, С6 — шз.

Поскольку любое твердое смещение является однородным решением, а С^, £ — 1,..., 6, — произвольные постоянные, то на основании (4.2) можно сделать следующие выводы:

Утверждение 4.1. 7 — 7° — 0, 7 — 71 — ¿т, 7 — 7-1 — -¿т являются собственными значениями спектральной задачи (3.2).

Утверждение 4.2. СЗ 71 соответствуют собственный и присоединенный векторы

а° — 01 — (М, 0)т, а1 — аз — (0, 0, -£1 - ¿£2)т. (4.3)

СЗ 7-1 соответствуют собственный и присоединенный векторы

а° — 02 — (1, -¿, 0)т, а1 — 04 — (0, 0, -£1 + ¿£2)Т. (4.4)

Утверждение 4.3. СЗ 7° соответствуют два собственных вектора

а° — 05 — (0, 0,1)т, а4 — а6 — (-£2,£1, 0)т. (4.5)

Однако приведенная система собственных и присоединенных векторов не исчерпывает корневые подпространства соответствующих СЗ. В [12] показано:

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

Ьт(¿т)а] — ¥], Мг(¿т— /] — 2, 3), (4.6)

¥2 — (2к, 2к, 0)т, /1 — (2кщ, 2кп2, 2т(1 + к)6„()Т, — -(2гтС + Вт)а2 - Са1, /3 — -¿^та2 .

дБ

Утверждение 4.5. Собственный вектор а2, помимо присоединенного вектора а2, ет еще два присоеди дующих краевых задач:

22

имеет еще два присоединенных вектора а2 и 03, которые определяются решением сле■

¿т(-¿т)а2 — ¥2, Мт(-¿т)а2 „/2 — 2, 3), (4.7)

дБ

¥2 — (2к, -2кг, 0)т, f2 — (2кщ, -2гкп2, 2т(1 + к)6„()Т,

¥2 — - (2«тСт + Вт)а2 - Ста2, f 2 — -Ста2

3 х 13 2 дБ

Утверждение 4.6. Собственный вектор аЦ — (а° 1, аЦ 2, з)т имеет только один

3/333 \Т

присоединенный вектор а1 — (а^ 1, а3 2, а3 з)т , который определяется решением следующей краевой задачи:

Ат а3 — 0, Ет а1

— (2кщ, 2кп2, 2(1 + к)6„т)Т. (4.8)

дБ

Утверждение 4.7. Собственный вектор а4 имеет только один присоединенный вектор а14, который определяется решением следующей краевой задачи:

Ата4 — 0, Ета1 — ( - т6„£2, т6„£ь -6„)Т. (4.9)

дБ

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

и1(£) — егт« 01, и2(£)— «1 (£), «з(£) — егт« (£01 + аз), «4(£) — «з(£),

£2 ч 2

£3 £2

«9(£) — егтМ ^01 + Vаз + £09 + 07 ) , «1°(£) — йд(£),

«5(£) — 05, и6(£) — а6, «7(£) — егтЧ £-01 + £03 + 09) , «в(£) — «7(£),

(4.10)

62

«11 (£) — £05 + 011, «12(£) — £а6 + 012.

Здесь

а1 — а1 — (1,2, 0)т, а2 — а^ — а1, аз — а1 — (0, 0, -()т, а4 — а2 — аз, 05 — а° — (0, 0,1)т, а6 — а4 — (-£2,6, 0)т, 07 — а1, ав — а2 — 07, (4.11) 09 — а2, 01° — а2 — а9, 011 — а1, 012 — а1.

Поскольку любая линейная комбинация векторов «г — 0 (1 — 1,..., 6) относится к классу «твердых смещений», отвечающие им напряженно-деформированные состояния тождественно равны нулю и, соответственно, векторы

Л = 0 (1 — 1,..., 6). (4.12)

Для векторов ( (1 — 7,..., 12) имеем:

Л7(£)— егт« (£ Ь9 + 67), (8 (£)— Л7(£),

Л9(£) — е^ 69, Л1°(£) — Л9(£), (4.13)

Лц(£) — 611, Л12(£) — 612,

0 = (6i3,b 623,í , 633,í) ,

- 1 bg = 03

O7 = 01 = (B3t + ÍTC )07 + C Og, 0g = 02 = 67,

1 2 - (4.14)

Ьд = &2 = (Взт - ¿тС)ад + СОз, Ью = = 6д,

ЬП = Ь3 = Взт О11 + Са5, Ь12 = Ь1 = Взт 012 + С а6.

4.2. Как известно [8], решение любой трехмерной краевой задачи для призмы можно представить в виде суммы решения Сен-Венана и решения типа «погранслоя». При этом главный вектор и главный момент напряжений в поперечном сечении призмы, отвечающих решению Сен-Венана, в общем случае отличны от нуля (это свойство и определяет понятие «решение Сен-Венана»), а главный вектор и главный момент напряжений, отвечающий «погранслою», равны нулю.

Решение Сен-Венана для ЕЗС будем называть линейные комбинации вида

6 12

US (£) = J] + J] - L)

¿=7

12

^S(£) = /X£Q- L). (4.15)

¿=7

Покажем, что постоянные (£ = 7,..., 12) можно выразить через компоненты главного вектора Q = (Q^i, Q^2,Q3)T = (Qxi, QX2,Q3)T и главного момента M = (M^i,M^2, M3)T = (Mxi, Mx2, M3)T напряжений ^3 = (^13,^23,^33)^, действующих в поперечном сечении S. Действительно,

Q^fc = J o^dS, M^i = J £2^33dS,

SS

M^2 = -J £i ^dS, M3 = J (^23 - 6 ^)dS, SS

где — проекции вектора напряжений в сечении £ = const на оси сопутствующей системы координат;

Qxi = Qii cos (т£) - Q^2 sin (т£),

Qx2 = Qíi sin (т£) + Q^2 cos (t£).

Введем расширенные векторы = (o¿, 0¿)T, = (6i3,¿, 623,^, 633,¿)T и рассмотрим скалярные произведения

= (vq+6, J v¿)i, = 1,..., 6. (4.16)

Учитывая, что = (o¿, 0)T, имеем

= (0q+6, 0¿). (4.17)

будем называть элементами матрицы жесткости . В случае, когда т = 0, опираясь на утверждения 3.3, 3.5, можно показать, что матрица жесткости имеет следующую

структуру:

¿11 0 0 0 0 0

0 ¿22 0 0 0 0

¿31 0 ¿33 0 0 0

0 ¿42 0 ¿44 0 0

0 0 0 0 ¿55 ¿56

0 0 0 0 ¿56 ¿66

где

¿11 — ¿22 — -¿33 — -¿44 — - ¿ — J £633,9

Б

¿31 — ¿42 — - J С&33,7 ¿5, Б

¿55 ^ У Ь33,11 ¿5, ¿56 — ^ 633,12 ¿5, ББ

¿66 — ^ (£1^23,12 - £2613,12) ¿5.

Б

Рассмотрим скалярные произведения ((, ад) (д — 1,..., 6). С учетом конкретного вида векторов ад, ( и структуры матрицы жесткости, получаем:

-егтп^¿С7 — 3*1 - ¿3*2, е^тп ^ [(¿31 + ^¿)С7 + ¿С9] — М^2 + ¿М^1, ^(¿55С11 + ¿56 С12 ) — 33,

^(¿5бСц + ¿66С12) — М3.

(4.18)

Здесь п — £ - ¿. Поскольку

- ¿3*2 — егтп(3x1 - ¿3x2), М^(£) + ¿М*1 — егтп (Мх2 + ¿Мх1),

Мх1 — Мг*1 - £3x2, Мх2 — Мг*2 + £3x1, где М*1, М*2 — изгибающие моменты в сечении £ — 1, получаем

^¿С7 — -3x1 + ¿3x2,

^¿С9 — М*2 + ¿М*1 + ¿31 (3x1 - ¿3x2).

(4.19)

Формулы (4.18), (4.19) доказывают, что постоянные С (^ — 7,..., 12) определяются только компонентами главного вектора и главного момента напряжений, действующих в поперечном сечении ЕЗС.

4.3. Основываясь на теореме о полноте системы элементарных решений [11], решение трехмерной задачи при отсутствии напряжений на боковой поверхности ЕЗС можно представить в виде

«(£) — «б (£) + «Р (£), ((£) — (Б (£) + (р (£).

(4.20)

где (£), (£) — решения Сен-Венана (4.15),

ир(£) = Е [С+м+(£) + С"«"(£ - О ,

к

(£) = Е [сМ(£)+(£ - о.

к

Здесь -ш±(£) = (м±(£), ^±(£))т — элементарные решения, отвечающие собственным значениям, и отличным от 70 =0 и 7±1 = ±гт; С±,Ск± — произвольные постоянные, которые определяются при удовлетворении граничным условиям на торцах £ = 0, £ = 1. При удалении от сечений £ = 0 и £ = 1 во внутреннюю часть ЕЗС (£) и ст±(£) экспоненциально убывают:

м+(£) = е^а+, м"(£ - 1) = еЪ«"'>а", *■+(£) = е7+«Ь+, о-"(£ - 1) = е^

Для обоснования принципа Сен-Венана необходимо доказать следующие утверждения:

1) главные векторы и главные моменты напряжений, отвечающие элементарным решениям ^±(£), равны нулю;

2) среди 7к± нет чисто мнимых.

Из первого утверждения вытекает, что напряженное состояние, отвечающее однородному решению мр(£) в любом поперечном сечении стержня, является самоуравновешенным (главный вектор и главный момент напряжений равны нулю), из второго утверждения следует, что оно экспоненциально убывает по мере удаления от торцов £ = 0, £ = 1 (имеет характер «погранслоев», локализованных у торцов).

Доказательство первого утверждения вытекает из утверждения 3.3 и следующей цепочки равенств:

(<, ^)1 = (а±, Ь,) - (Ь±, а,) = -(Ь±, а,) = 0, I = 1,..., 6.

Здесь учтено, что Ь, = 0 (^ = 1,..., 6).

Для доказательства второго утверждения умножим уравнение (3.2) скалярно в Н на вектор а. Получаем относительно 7 квадратное уравнение

5о72 + 517 - 52 = 0. (4.21)

Здесь 50 = (Са, а), 51 = (Вта, а), = (Ата, а). Преобразуем 51:

51 = (Вта, а) = (В3та, а) - (В3та, а) = (В3та, а) - (В3та, а) = 2г 1т(В3та, а) = 2гд.

Так как 51 — чисто мнимая величина, то для доказательства второго утверждения достаточно показать, что дискриминант уравнения (4.21) Д = +4до52 =4(-52+5о52) > 0. Представим до, 5, 52 в следующем виде:

5о = с2 + 2с2 + 2с2,

5 = 1т

2 V -ао + г^а3 , + а3 + та^ + 2т(Ва3, Й3)

52 = 2/2 + /22 + 2/32 + 2^2 + 4

Здесь

ci = ||ао||, С2 = ||аз||, сз = ||аз||к, ао = ааеа, а * = (Dai - a2)ei + (Da2 + ai)e2,

fi =

V -ао + rDa3

/2 =

V аз + та*

/з = т У£аз|

¿2 — Цда^Ц2 + ||д2 а2^2, ¿2 — Цд1а2 + д2й1|2, V— е^ + е2$2, (а, Б)к — а ■ Б ¿5, ||а||К —(а, а)к.

Б

Используя неравенство Коши — Буняковского, получаем следующую оценку для д:

о о

|д| ^ 2(V -ао + т^аз,аз)к + (V аз + та *, ао) + 2т(£аз,аз)

< 2/1 ■ Сз + /2 ■ С1 + 2/з ■ С2.

Следовательно,

j ^ - (2/i ■ сз + /2 ■ ci + 2/3 ■ С2)2 + (ci + 2с2 + 2с3)(2/2 + /22 + 2/32 + 2d2 + ¿2)

2 , J2)(„2 , о „2

— 2(/1С1 - /2 сз) + 4(/1С2 - сз/зУ + 2(/2С2 - /3с^ + (2¿2 + ¿2) (с2 + 2с2 + 2с2) ^ 0.

Равенство А — 0 возможно лишь при условии, что /1 — /2 — /з — 0, ¿1 — ¿2 — 0, которые выполняются тогда и только тогда, когда вектор ао является линейной комбинацией векторов а1, а0, аЦ, а4. Следовательно, дискриминант обращается в ноль только на собственных векторах, соответствующих собственным значениям 70 — 0, 7±1 — ±ят, и, следовательно, остальные собственные значения не могут быть только чисто мнимыми. Из выражений (4.20) вытекает, что

Лр — Л(£) - ЛБ(£)— О(е-7*0), 7* — шт(Ие7+), в — шш(£,1 - £),

к *

т. е. имеет характер погранслоев, что обосновывает принцип Сен-Венана.

Литература

1. Saint-Venant B. Memoire sur la torsin des prisms, avec des conciderations sur leur flexion, ainsi que sur l'equilibre interieur des solides elastiques en general, et des formules pratiques pour le calcul de leur resistance a divers efforts s'exercant simultanement // Mem. Savants Etrang.—1856.—Vol. 14.—P. 233560.

2. Saint-Venant B. Memoire sur la flaxion des prismes, etc. // J. Math.—Liouville, 1856.—Vol. 1.—P. 89189.

3. Сен-Венан Б. Мемуар о кручении призм. Мемуар об изгибе призм. «Классики естествознания».— М.: Физматгиз, 1961.—518 с.

4. Mathieu E. L. Theorie de l'elasticite des corps solides.—Paris, 1883.

5. Clebsch R. Saint-Venant's principle // Arch. Ration. Mech. and Analysis.—1965.—Vol. 18, № 2.—P. 8396.

6. Олейник О. А., Иосифьян Г. А. Об условиях затухания и предельном поведении на бесконечности решений системы уравнений теории упругости // Докл. АН СССР.—1981.—Т. 258, № 3.—С. 550553.

7. Устинов Ю. А. К обоснованию принципа Сен-Венана // Изв. вузов Сев.-Кавк. рег.—1994.—С. 9192.

8. Устинов Ю. А. Задачи Сен-Венана для псевдоцилиндров.—М.: Физматлит, 2003.—125 c.

9. Ворович И. И., Бабешко В. А. Динамические смешанные задачи теории упругости для неклассических областей.—М., 1973.—320 с.

10. Костюченко А. Г., Оразов М. Б. Задача о колебаниях упругого полуцилиндра и связанные с ней самосопряженные квадратичные пучки // Тр. семинара им. И. Г. Петровского.—1981.—Вып. 6.— С. 97-146.

11. Гетман И. П., Устинов Ю. А. Математическая теория нерегулярных твердых волноводов.—Ростов-на-Дону: Изд-во РГУ, 1993.—144 с.

12. Друзь А. Н., Поляков Н. А., Устинов Ю. А. Однородные решения и задачи Сен-Венана для естественно закрученного стержня // ПММ.—1996.—Т. 60, вып. 4.—С. 660-668.

Статья поступила 12 апреля 2009 г.

УСТИНОВ ЮРИЙ АНАТОЛЬЕВИЧ Южный федеральный университет, проф. каф. теории упругости

РОССИЯ, 344090, Ростов-на-Дону, ул. Мильчакова, 8-а; Южный математический институт ВНЦ РАН и РСО-А, гл. науч. сотр. лаб. мат. методов механики сплошной среды РОССИЯ, 362027, Владикавказ, ул. Маркуса, 22 E-mail: [email protected]

JUSTIFICATION OF THE PRINCIPLE OF SAINT VENANT FOR A NATURALLY TWISTED ROD

Ustinov Yu. A.

The paper gives the mathematical evidence of the principle of Saint-Venant for a naturally twisted rod.

Key words: principle of Saint-Venant, spectral problem, eigenvalues, elementary solutions, generalized orthogonality.