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

книги / Механика деформаций гибких тел

..pdf
Скачиваний:
4
Добавлен:
12.11.2023
Размер:
4.96 Mб
Скачать

CnM=

Cfnn Cfim’

Так как Эз — единичный вектор,

то а® = a^, а®-аз,п = Спз = 0.

Формулы вида (1.1.3) справедливы только при строчных значешшх индексов:

Cnm = \|2a^^^ (йтк.п + ^nft,m — ^nin,*).

(1.2.5)

C номощыо (1.2.3) устанавливаются следующие правила диффе])енциропанш1 произвольного двумерного вектора Av(X)=W^ai,= = WM заданного в ко- и контравариантпом базисах:

^nAV =

=

a^'дnWмf

 

OnW^ =

W^n +

CnMW^^f

(1.2.6)

дпЮм =

WM,п — CnMWL-

 

Здесь 0„ — оператор ковариаптиой производной по переменной .с’ относительно начального двумерного базиса.

^^ормулы

d ^ = jdx^dx^, 7 — (detta;,,v])‘'^

d'^ =

Idx^, Cv = Cv •а”,

dX = (Cnmdar^Jx'")

задают в окрестности точки х элемент площади d^, ограничеи- HbTii координатными линиями, элемеит длины граничного конту­ ра d*^ II расстояние dX между двумя бесконечно близкими точKaAIII X и X + dX.

2.2.Формулировка кинематических уравнений

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

точки в мгиовеппом состоянии принимает значение х = х + ц. Начальпый базис а^(х) преобразуется в мгповеипый базис Э;/)(х).

Векторы а„, = 5„х = а„ + 5„ц = а„ -4- Wn в любой момент распола­ гаются на мгновенной касательной плоскости. Вектор а,) может отклоняться от мгновенной нормали вследствие поворота мате­ риальной точки.

Из полного преобразования базиса выделяется жесткий по­ ворот, переводящий начальпый базис а^(х)в повернутый ал)(х). Деформация двумерного континуума не может изменить длину нормального вектора Эз, по может повернуть его как жесткое целое. Это условие выражается равенством

Два остальных вектора мгповеппого базиса в общем случае не совпадают с соответствующими векторами повернутого базиса. Векторы приращения

и „зап )-ап |

определяют поле деформаций двумерного континуума, согласовап-

пое C условием

(1.2.7).

 

Преобразование поворота выражается через натуральное поле

поворотов у(х)

согласоваппымя с (1.1-8) уравпеппянп

 

 

Эп! = а„ + ф,у Х а» + (psV Х (у X ап),

 

 

On = а„, - ф,у X а„) + ф2v X (V X 0«,),

(1.2.9)

Ф,з(ф)-‘ Sinф, ф*®(ф)-*(1-со8ф), ф = (VV)

II подчиняется условиям

Од'] •а.лг] = ад- •ам ^ Ядглг,

(ад-] X Пдг]) •BI,] = (ад- X Плг) •Эг, S ditML,

так что вдгм — метрический, ак»ь — дискриминантный тензоры начального и повернутого базисов одиовромеппо. Кроме того, по­ вернутый базис сохраняет присущее начальному базису свойство ортогональности

а „ ,а ,, = 0.

(1.2.10)

Согласованное с (1.1.12) равенство

^па,г, = Л '„Х а„,

(1.2.11)

вводит тензорное поло изгибаний ТЛ»(х) и определяет его через поле поворотов нелинейным уравпепием

У„ = фАУ + фгУХ5„У + (фа5„ф)У, ф* = (ф)"‘ (1-ф 1)- (1.2.12)

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

Тензорное поле деформаций (1.2.8) выражается через поля перемещений и поворотов эквивалентными пелипейпыми урав­ нениями

«» = ^nU— ф,у X а„ - фгУ X (у X а„) =*

 

 

= 9 „и -ф ,уХ ап ) + ф2УХ (у Х а „ ,).

(1.2.13)

Поля изгибаний и деформаций подчиняются выводимым по­

добно (1.1.18)

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

 

 

 

а’"’"(5 „ У „ -1 / 2 У „ Х У „ ) = 0 ,

 

 

 

’""Ч 5 « « « + У „ Х а „ ,)= 0 .

 

При введении вспомогательных векторных функций у" =

Vm,

и" =

они приобретают более простую запись: дпЧ" — 1/2У« X

X у" =

0. ^nU" -

а„, X у“ = 0. ‘

 

Для скоростей измедспия кидсматических деремеплых слраведлдвы аналогичные (1.1.26) — (1.1.28) равенства

ба*-, = боУ X ал], боУ = ф,6у + ср,у X 6у + (ср*бф)у,

б\'„=*5„боУ + боУХУ„, би„ = 5„би-боУХа„1, (1.2.14); боУ„ = 5„боУ, боН„ = 5„би — боУ X а„,.

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

 

U = н.ма” = Wjfja"*,

у = Ома" =

 

 

 

 

U,. =

н„.„а" =

 

У„ =

7я.„а" =

7„,г,

 

 

и" =

н""а.и = н""*а.и„

у" =

в""алг =

н""*азп

 

 

(по опредсленшо

 

Переход от одного из них к другому

осуществляется

с

помощью уравнений (1.2.9),

дредставлеппых

в виде

 

 

 

 

 

 

 

 

 

 

 

а.у)

*

(я.у.м "Ь ф]\’.\г) а " ,

ал/ =

(йл-,1г +

фл‘.«) а^ *

 

( 1 .2 . 1 5 )

и вводящих тсп.'юр поворотов

 

 

 

 

 

 

 

ф/» = ф|У X ад-+ ф,у X (УX а^) * ф,у X aJfJ - ф^у X (у X а^,,)

(1.2.16)

C компонентами

 

 

 

 

 

 

 

 

 

 

Ф.У.М*

ф.%- •ад, =

ф,ал.«ьУ^ + ф2(олгНд/ — aкмVl.v^‘).

(1.2.17¾

Наря,\ C

оператором

д„ ковариаптиого

дифференцирования

в начальном ба.энсе может быть

введед оператор

 

ковариаит-

ного дифференциропаипя в повернутом базисе.

При

этом для

любого вектора лу(х) будут выполняться равенства

 

 

5,лу =

л’''д„и;к =

Икд„и''^ =

 

 

 

 

 

 

люторые вместе с (1.2.15) дают

следующие связи между кова-

рпаптны.лш пропзводны.ми в разных базисах:

 

 

 

 

^njiWti =(Лья- + Фьк)5„м;^

дniW^'^ =

 

+

 

 

 

Из равенств

й„лу^ 5„(ш"*алг1) = а.и)5„,м>"*

и

(1.2.11)

выводится

с1)0рмула

 

 

 

 

 

 

 

 

 

 

(1.2.18)

 

 

 

 

 

 

 

 

 

 

 

позволяющая

вычислять

ковариаптпые производные

от

вектора

в довернутом базисе непосредственно.

 

 

 

 

 

 

Ковариаптпые производные от векторов и тензоров в началь­ ном базисе вычисляются по согласоваппым с (1.2.6) формулам:

дп^М — CnMlVLt

^nlVMK = Шдис.п — CI^K ICML CnMlVLK,

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

Из уравнения

(1.2,12)

определяются

тензора

пзгибанин

 

 

 

 

 

VnL = (ф1Лла +

Фгамх.кУ"’ )

-I- (фз^пф)

, 1 .2 .1 9 )

1'^ПЛГ 1=

(Ллг к +

фм к )

VnX,,

(р = (VL V^')'

 

из (1.2.13) — компоненты тензора деформации

 

НпЬ = М-’„Х,-ф„ь,

ИпМ] — (блГК + ф.МА-)а^'*Н.1Г., WnL^dnJtL.

(1.2.20)

Ураппешш совместности деформаци1т (условия сплошности двумерного моментного континуума) формулируются в поверну­

том базисе:.

 

 

 

5„цпл/] + а^.к (^nL -H UnL])

=

о,

= а^^^а^>'”'иттл‘

Для скоростей изменения кипематияескн.у векторов естествсн-

пымп являются согласованные с

(1.1.30)

разложения

боУ = а " ’бУлг,

би *= а''би.>г,

а„боУ =

боУ„ = а*” бУ„м1,

^„би - боУ X а„) =

боН„ =

a*'^бUn.^fl.

Компоненты вектора поворота и вектора перемещения обра­ зуют систему шести независимых кнпематпчоскпх параметров двумерного момептпого континуума. Формулы (1.2.19) и (1^2.20) выражают через эти параметры компоненты тензоров изгибаний H деформаций.

Из последовательной цепочки равенств

а„, •а„ * а*, - а„ = (а„, -I- и„) •аз, = и„ •а„ * u ,,,

следует, что мгновенный базис нс сохраняет, вообще говоря, свойство ортогональности повернутого базиса п компоненты u^si тензора деформации являются мерой отклонения мгповеппого вектора аз, и совпадающего с ппм повернутого вектора аз, от мгповенной нормали к поверхпости. В отличие от мгновенных векторов векторы 0«,, ортогональные Эз,, в таком случае пе мо­ гут находиться в плоскости, касательной к деформироваппой поверхпостп

2.3.Формулировка динамических уравнений

Ввиду отсутствия зависимости от перемеппой х* выра­ жение (1.1.41) для удельной мощности деформации двумерного

коитипуума запишется так:

 

P = Z- . (5 „ б и -б ,У Х А п ,)-Ь Y -

(б ,У Х А „ ).

Равенство пулю этой величины для группы жестких перемеще­ ний обеспечивается дополпительпой динамической связью

Z^= O,

(1.2.22)

вполне совместимой с понятием двумерного континуума, не име-

ющего сечепнп х = const, которых определено векторное поле Z^

Согласованное с условием (1.2.22) значеппе удельной мощ­

ности впутренннх сил выражается формулой

 

 

P = Z -. (а„би - боУ X а„,) + Y - . 5„боУ.

(1.2.23).

Здесь P — поверхностная

плотность

мощности

внутренних сил,

рассчитанная на единицу

начальной

площади;

Z-(х)

п У "(х) —«

поля внутренних усилий II моментов, распределенные по сече­ нию X” = COHSt II рассчитанные на единицу начальной длины. Для обозначения поля внутренних моментов оставлена пропис­ ная буква, поскольку строчная буква используется для обозна­ чения традиционного поля моментов оболочки.

Внешние механические воздействия па двумерный конти­ нуум могут HbiTIi распределены по поверхности ^ и/илп по кон­

туру

Сн.чы H аюмеиты, распределенные по поверхности, изме­

ряются

поверхностными нлотностямн Z(X) и Y(X), рассчитан­

ными на единицу начальной площади. На граничном контуре до­ пускаются силовые II кинематические воздействия. Силы и мо­ менты, раенределенные по участку контура, измеряются

контурными п.'ютпостямп Zv(X) II Yv(X), рассчитаппымп па еди­ ницу начальной длины. Кпне.матическпе воздействия па участке

задаются но.’1ямн перемещений и поворотов н(х) и у (х ). Пользуясь принципом Даламбера, условимся включать по­

верхностные плотности инерционных сил п моментов в поверх­

ностные поля

г(х)

II

Y(X).

 

 

 

Аналогично

(1.1.^^5)

и (1.1.16)

можпо

ввести определенные

на контуре

поля:

 

 

 

 

 

 

 

Vx

 

= F

Vx е ® '” ;

 

 

V

^ c ^

~ \ \

V x e ®^'^;

к-пг”

V x e « ’“;

_ к Л

V x e ® '“ ;

 

Zv

V x e ® ’•^;

 

Y v

V x e ¾ ’+

Из (1.1.47) следует уравнение баланса механической энер­ гии двумерного момептпого континуума

J (^.би + У.б„у - р ) 4 ^ + J (г.-бн'^ + YVaoV^') dW = 0.

(1.2.24)

C учетом выражения (1.2.23) оно может быть преобразовано к виду

I* ((г + ^«г-).би + (Y ^-a»)Xz" + 5,^Y'^).бov)

J (^vnZ-. (би— би) "Ь ^vnY •(б^у —

dW 4"

■уравнения (1.2.24), (1.2.25) справедливы п для внртуальны.х: кппематическпх полей. Поэтому их следствием являются ло­ кальные динамические уравнения

5„г" + г = 0, №

+ а „ )Х г" + У = 0

(1.2.26)

H граппчпые условия

 

 

 

 

U =

U,

V = V

V xe^ g'-;

(1.2.27)

CvnZ" =

Zv,

^vnY" =

Yv V x e ¾ '^

(1.2.28)

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

Z =

Z OJ

 

 

 

м у

 

°

Zv а.

Щ.

— J V

— 1

V а^^^,

(1.2.29)

Zv =

-'M]^

 

"а.и = г""'алг1,

Y' =

=

У""'амз.

 

Здесь г"“ и г " " ’ — компопепты

тензора виутропппх yciuinii, У ""

H У ""! — моментов.

 

 

 

 

 

G помощью формул (1.2.15) из равенств (1.2.29) устанавли­ ваются связи между компонентами по разным базисам, в част­ ности

^уплг _ упЧдЛГК

Динамические уравнения (1.2.26) п граничные условия (1.2.27), (1.2.28) в начальном базисе имеют скалярную формулировку:

 

5„г’“'^ +

г-'^ =

0;

 

 

+ (а^1гк +

a^^KWnL)

=

0;

Wjtf = Wjif,

Wjtf ~

W-V

Vx S

 

(1.2.30)

^vnZ

 

 

 

 

 

 

в повернутом —

 

 

 

0;

 

 

 

 

 

=

 

 

+ ( Л к +

ayiu nLi)

 

=

0;

WifJ =

Wjtf],

V u =

V u

V x е <8?-;

 

(1.2.31)

 

 

 

=

V x e i ? + .

Присутствующие здесь ковариантные производные тензоров вы­ числяются по согласованной с (1.2.18) формуле

Выражение (1.2.23) с помощью разложений (1.2.29)' преобразуется iк виду р = 2""^бЦпМ1 + У"^^бУ„лг], показывающему, что

46

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

Сформулированные кинематические и динамические уравне­ ния составляют систему механических уравнений двумерного континуума (материалыюн поверхности) Коссера.

2.4. Формулировка термодинамических и определяющих уравнений

!Сели ограинчиться изучением термомехапичеехшх про­ цессов деформирования двумерного момептиого континуума, то систему его лгеханнческих уравнений следует дополнить форму-

лпровко!! двумерных термодинамических и определяющих урав- IieiiIiii. Их можно получить пз соответствующих уравпешш трех­

мерного континуума путем формального уменьшения размерно­ сти пространства.

Для двумерного термомсхапического континуума постулиру­ ется существование новерхпостных плотностей внутренней энер­ гии з(х) II энтропии .S(X) II скалярного поля абсолютной темпе­ ратуры ^ (х )> 0 . 13 качестве параметров, характеризующих изме­

нение Tennonoii энергии, вводятся скалярное поле д(х) скорости HOBcpXHocTUoii н.лотности поступления внешнего тепла и вектор­ ное поле д(х) скорости внутреннего теплового потока через едпницу начальной длины произвольного внутреннего контура ма­ териальной поверхности.

В соответствии с (1.1.01) первый постулат термодинамики двумерного континуума формулируется как уравнение притока

тепла к произвольно!!

области Д.^, ограниченной контуром А*??;

б /

J

= J (p + Я )d^-

C v J - (1.2.32)

\Д^

 

J

 

Второй постулат выражается следующим из (1.1.62) неравенст­ вом диссипации

Ь1 J

J

J e v q r W .

(1.2.33)

/

а ЗЗ

д'®’

 

Для дифференцируемых полей д(х) и *(х) возможно преоб­ разование (1.2.32) и (1.2.33) к однородным интегральным формам

J

ЬгйЗЗ =

J

(р -f 5 — diV д)

 

39

а39

 

J

(^ O s-g +

Oiv д — д -Ь )й ^ ?> 0 .

(1.2.34)

а ЗВ

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

б2 == P + g — diVд,

(1.2.35)

—- g + diV g — g •h ^ 0.

При заппсп (1.2.34) и (1.2.35) введены операторы вычисления дивергенции и градиента в двумерном пространстве поверхности isf и коллппеарпый температурному градиенту вектор h =*

=grad t — grad In t.

Любой термомехаппческий процесс деформации двумерного момептпого континуума должен удовлетворять условиям (1.2.32) H (1.2.33) плп IIX модпфицироваппым эквивалентам. C ними должна быть согласована структура двумерных определяющих

уравпеппй.

Рассмотрим, например, двумерный континуум без внутрен­ ней диссппацип (совершенный, плп термоупругий):

г~ = tds — Q+ diV g = 0.

(1.2.30)

За основной термодинамический потенциал примем поверхност­ ную плотность свободной энергии у ^ z - st. При этом условие (1.2.36) выражается равенством 6у = р — sSt п допускает уравпспие состояния вида

У = UiUnUU PnMi, О

(1.2.37)

II определяющие уравнения

 

ду/дУпщ , S = -O yfdt.

(1.2.38)

Последние должны удовлетворять термодннамическпм ограниче­

ниям типа (1.1.82).

Благодаря (1.2.36) неравенство диссипации (1.2.35) форму­

лируется в виде

= —g •Ii >

0 п может быть

обеспечено за

счет представления

д =

- г - Ь ,

(1.2.39)

 

где г — неотрицательный тензор теплопроводности, который мо­ жет явным образом зависеть от определяющих параметров Unarj, PnMi, t и от вектора h:

г = !(UnJfJ, PnM], t, h ).

(1.2.40)

При условии (1.2.39) равенство (1.2.36) обретает смысл урав­ нения теплопроводности.

Уравнения (1 .2 .3 6 )-(1 .2 .4 0 ) замыкают систему механиче­ ских уравнений и вместе с ними формулируют нелинейную за­ дачу термоупругости для двумерного моментного континуума.

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

рая может служить математической моделью деформируемого физического тела, если эмпирически обнаружены некоторые закоиомериости, математическое выражепие которых позволяет установить определенное соответствие между трехмерным и дву­ мерным коитинуумами. Например, закономерности, которые поз­ воляют моделировать двумерным коитппуумом плоскую дефор­ мацию физического тола, отличаются от проявляюп^ихся при мо­ делировании изгиба топких физических тел. Поэтому формули­ ровка определяющих уравпепии двумерного континуума должна опираться на соответствующие уравнеппя трехмерного континуу­ ма и учитывать эмпирические связи, идентифицирующие их. Такая ирогралгма реализована в следующей главе при модели­ ровании деформации материальной оболочки. Здесь уместно за­ метить, что д.’гя специальным образодг построенной системы тер­ модинамических и определяющих уравнений двумерный конти­ нуум Koccepa обретает смысл математической модели оболочки C жесткими поперечными волокнами.

2.5. Материальная

поверхность Кирхгофа

 

и безмоментная поверхность

 

Шестикомпопептный

тензор деформаций

не мо­

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

метризации — вырождение

в четырехкомпопентный

тензор. Это

достигается наложением кинематической связи

 

Un

■Яз] = U nii ~

( 1 . 2 . 4 1 )

которая устанавливает две скалярные зависимости между век­ торами поворота и перемещения. Иначе говоря, связь (1.2.41) сокращает число кинематически независимых параметров с ше­ сти до четырех, определяя две компоненты у„ вектора поворота через четыре параметра Нл-) Уз. При этом динамические парамет­

ры

приобретают смысл силовых реакций, отвечающих

связи

(1.2.41), и находятся из динамических уравпепии (1.2.31),

тог­

да как определяющие уравнения для

них теряют смысл.

 

 

 

Равенства (1.2.41) в совокупности

с (1.2.7) указывают

на

то,

что мгновенный базисный вектор Зз в любой момент времени на­ правлен по нормали к деформированной поверхности. А именно это условие выполняется для материальной поверхности Кирх­ гофа, моделирующей деформируемую оболочку. Но наложение связи (1.2.41) не превращает еще построенную модель двумер­ ного моментпого континуума в модель Кирхгофа. От последней она отличается наличием независимого поворота относительно пормали и наличием нормальных компонент у тензора внут­ ренних моментов. Для полного совпадения с моделью Кирхгофа

надо дополнить условие (1.2.41) кинематической

связью

= О

(1.2.42)

и динамической связью

 

Первая

из

них симметризует четырехкомпопептпый тензор Un^j

и вместе

с

(1.2,41) образует систему трех уравнений, связываю­

щих компоненты вектора поворота с компонентами вектора пере­

мещения. Вторая налагает соответствующие ограничения

на

формулировку определяющих уравнений.

 

Каждая пара связей (1 .2 .4 1 )-(1 .2 .4 3 ) и каждая из них

в

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

нуум. Совокупность же

всех

трех порождает

псевдомомеитный

континуум — поверхность

Кирхгофа. Каждая

из связей (1.2.41)

и (1.2.43) сокращает число

скалярных граничных условий на

единицу, а связь (1.2.42)

не изменяет их числа. Поэтому конти-

нуул1 Кирхгофа допускает постановку па граничном контуре лишь четырех скалярных условий.

Сочетание связей (1.2.41) и (1.2.42) образует замкнутую си­ стему уравнений для определения компонент вектора поворота через компоненты вектора перемещения и приводит к формули­ ровке двумерного аналога псевдоыоментного континуума со стес­ ненным вращением, который требует задания пяти скалярных условий на граничном контуре.

Псевдомомеитный континуум с сочетанием связей (1.2.42) и (1.2.43) также требует задания пяти контурных условий и от­ личается от поверхности Кпрхгофа наличием кпиематически не­ зависимого поля угловых перемещений нормали. Такой копти-

нуу.м моделирует деформацию

оболочки,

сопровождающуюся

сдвигами

в ее нормальных сечениях (поперечными

сдвигами).

Двумерный континуум с

сочетанием

связей

(1.2.41)

и

(1.2.43)

, требующий задания

четырех контурных условий,

отли­

чается от поверхности Кирхгофа наличием независимого пово­ рота отпосптельпо нормалп.

Добавление к кинематическим связям (1.2.41) и (1.2.42) бо-

.чее жесткой, чем (1.2.43), дппампческой связи

 

Y" = О,

(1.2.44)

совместимой с заданными условиями

 

Y = O, Yv = о .

(1.2.45)

приводит к модели двумерного безмомептного континуума — безмо.ментпой материальной повер.хности.

Следствием условий (1,2.41), (1.2.44) п (1.2.45)

являются

равенства

(1.2.46)

. а=*) = 2" • = о, а„) X 2" = о,

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

Согласно (1.2,28) скалярные

равенства (1.2.46)

совместимы

C заданным условием

4 = о ,

 

Zv•a®^ =

(1.2.47)