Chizhov_lektsiya_4 (1183951), страница 2
Текст из файла (страница 2)
Изменение кинетического момента.Теорема об изменении кинетического момента частицы вещества постоянного составав элементарном объеме V также может быть представлена в форме балансногосоотношения. Пусть масса частицы m V , а ее скорость v vi ni .Проекции вектора кинетического момента частицыL r mvна орты декартовой системыFx v 02S 02sin ,Li ni L m ijk x j v j .В соответствии с теоремой динамики, скорость изменения вектора кинетическогомомента частицы определяется суммой моментов внешних сил. В рассматриваемомслучае – объемных FV и поверхностных F , действиующих на частицу:dL r FV r dF .dtЗдесь dF dFknk nk kn dSn - элемент поверхностной силы.Для проекций на орты ni получимdm ijkx j vk ijk x j f k V ijk x j kn dSn .dtПереходя к интегрированию по объему и используя теорему о среднем, получимd ijkx j vk ijk x j gk ijkx j kn .dtxnЗдесь выражение приведено для объемной силы – силы тяжести f i gi V .Воспользуемся и здесь уравнением непрерывностиd divv 0 ,dtумноженным на векторное произведение r v (записанное в тензорных обозначениях):vd ijk x j vk ijk x j vk n 0 .dtxnСкладывая оба уравнения, получимd ijk x j vk ijk x j vk ijk x j kn ijk x j gk .dtxn5Раскрывая выражение для субстанциальной производной, получим дивергентнуюформу теоремы об изменении кинетического момента: ijk x j vk ijk x j vk vn ijk x j kn ijk x j gk .txnТеперь осталось выполнить интегрирование по объему наблюдения, не зависящему отвремени, и применить теорему Гаусса к дивергенции тензора от плотности потокаимпульса ik vi vk ik , чтобы получить балансное соотношение для кинетического моментавещества Li ijk x j vk dV в объеме наблюдения:VLi ijk x j kn dSn Mi ,tгде Mi ijk x j gk dV - момент объемной силы тяжести.V4.3.3.
Симметричность тензора напряжений.Воспользуемся теоремой об изменении кинетического момента ijk x j vk ijk x j vk vn ijk x j kn ijk x j f k ,txnчтобы доказать симметричность тензора напряжений.Умножим уравнение Эйлера vk vk vn kn f ktxnна ijk x j и вычтем из уравнения момента. Выполняя дифференцирование, получим ijk vkx jx jt ijk vk vnx jxn ijkx jxn kn 0 .x j jn , поэтому ijk vk v j ijk kj 0 .txnТак как ijk vk v j 0 в силу симметрии, то ijk kj 0 для любого тензора напряжений, аНо0, аследовательно, kj jk .4.4.
Криволинейные ортогональные координаты.Уравнение Эйлера можно записать в ортогональных координатах, если заданыплотность объемных сил и поле скоростей заданных проекциями на локально-декартовыорты f f k ek , v vk ek . Дифференцирование вектора скорости даетde dvdv dvi ei vi i i ei vi ei .dt dtdtdtЗдесь вектор - скорость изменения ориентации ортов при перемещении частицы.В полярных координатах вектор скорости задается проекциями v vr e1 v e2 vz e3 , а e3 ve3 , поэтому уравнение имеет видr dv v2 vv p p p dv dv r e1 r e2 z e3 f r e1 f e2 f z e3 .r r dtr r z dt dt Субстанциальная производная также должна быть записана в криволинейныхкоординатах:d v vr v vz .dt ttrrz6Аналогично записывается и уравнение в сферических координатах, где скорость задаетсяпроекциями v vr e1 v e2 v e3 , а скорость изменения ориентации ортов cos e1 sin e2 e3 vctg e1 ve2 ve3 .rrrСубстанциальная производная в сферических координатахd . v vr v vdt ttrrr sin Обобщенные криволинейные координаты.Рассматривая сплошную среду, как систему материальных точек, можно воспользоватьсялагранжевым описанием движения в произвольных криволниейных координатах.Уравнение движения элементарной массы m постоянного состава, которая в данный моментоказалась в точке наблюдения.
Второй закон Ньютона для частицы массой m (элементарноймассы постоянного состава) имеет вид:mdv f R .dtРазбивая силы на поверхностные и объемные, запишем соответствующие выражения для этихсил:f mg gV , R pdS pV .Векторное уравнение движения в простейшем случае имеет вид: dv g p .dtПереход к эйлерову описанию проводится путем задания положения точки наблюдения, вкоторой оказывается частица в момент наблдения t. Пусть r вектор, задающий положение точкинаблюдения /выделенного элементарного объема V /. Введем обобщенные координаты q i ,задающие положение рассматриваемой точки пространства: r r qk .Для записи уравнения движения в проекциях на локальные направления r порождаемые r q k , t криволинейными координатами r q i (например, цилиндрическими илиq iсферическими), вычислим скалярные произведения 1 dv r r g r p .dt Выражения для поверхностных и объемных потенциальных сил g представимы в видеградиентов — производных по соответствующим обобщенным координатам r r g g r q k q k ,qk qk r pr p p r qk qk .q k q kДля вычисления проекций ускорения преобразуем левую часть, используя тот факт, что полескоростей отражает перемещение частиц за элементарный промежуток времени, а значит, связано r q k , t q i под знак полнойс изменением обобщенных координат.
Внося выражение r q iпроизводной d r d r r dv dv r qk q k dt q v v dt qdt k q k dt k .7 drrrРассматривая поле скоростей v как функции независимых переменныхqk dt q ktv qk , q k , t - обобщенных координат qk , скоростей q k и времени t, получим полезноеvrсоотношение,, которое позволяет преобразовать левую часть уравнения кq k q kкомпактному виду: d v v dv v v r qk dt dtqk qk d v 2 qk dt q k 2 qk v 2 . 2 Уравнения движения в криволинейных координатах получаются из соотношения:d v2 v2 1 p. dt qk 2 qk 2 qk qkДля получения уравнения Эйлера в криволинейных координатах необходимо выразитьобобщенные скорости q k как функции положения точки и проекций физической скорости навыбранные направления v( qk ) ~ r v qk , qk , t .
Для локально-декартовых координат(цилиндрических или сферических) в качестве таких величин используются проекции вектораскорости на соответствующие орты.Пример. Уравнение эйлера в сферических координатахr rer , v2 vr2 v2 v2 r 2 r 2 2 r 2 sin 2 2 .Вычислим производныеvv v2 v2 2 v2 22,, r vr , r r sin r 2 2 r 2 r sin v2r 2v2 v2222, r sin r v2 2 v2 22rsincosvctg, 0. 2 2 Учитывая выражение для субстанциальной производнойvd v , r vr dt tr tr r r sin проекции уравнения Эйлера на орты сферической системы можно записать в виде системы:v vr v2 v2dvr vrv v v1 p , vr r r dttrr r sin r r rv v vr v v2dv vv v v1 p , vr ctg dttrr r sin rr r rdvdtvt vrvrv v vr v v vpv v11 .ctg r r sin rr r sin r sin 8.