Loading [MathJax]/jax/output/HTML-CSS/jax.js

Страница сетевой поддержки учеников Тимура Юрьевича Альпина

07.04.2020

Задания и материалы для дистанционного занятия по векторному анализу гр. 06-822 в ср 8.04.2020, 10:10 и гр. 06-812 в пн 13.04.2020, 10:10 (Анчиков, №189, 191)

Криволинейные координаты и их орты

Положение точки в трёхмерном пространстве задаётся тремя числами. Но каким именно образом - тут возможно множество вариантов, и декартовы координаты - только один из них. Уже (хотя бы в курсе мат.анализа) использовались цилиндрические и сферические координаты, а ими возможности отнюдь не исчерпываются. Линии, на которых постоянна некоторая пара координат, и меняется третья координата - для декартовых координат являются прямыми, но в более общем случае - не обязательно, поэтому в общем случае координаты называются криволинейными. Например, если в сферических координатах зафиксировать координаты r и φ, и менять θ -- мы получим вертикальную полуокружность.

Пусть положение точки в декартовых координатах определяется координатами (x,y,z), но есть и другая система координат, в которой положение точки определяется координатами (p,q,s), причём функции {x=x(p,q,s)y=y(p,q,s)z=z(p,q,s) дифференцируемы. Радиус-вектор точки r=ix+jy+kz можно проинтегрировать по каждой криволинейной координате, например, по p: rp=ixp+jyp+kzp, Коэффициентами Ламэ называются модули полученных таким образом векторов: Hp|rp|=(xp)2+(yp)2+(zp)2, Нормируем на единицу производне типа rp: eprpHp,|ep|=1 Аналогично, Hq|rq|,Hs|rs| eqrqHq,esrsHs Векторы ep,eq,es -- орты новой системы координат. У ортогональных систем они имеют такие свойства: (epep)=1,(eqeq)=1,(eses)=1, (epeq)=0,(epes)=0,(eqes)=0, [ep×eq]=es,[eq×es]=ep,[es×ep]=eq. Например, для цилиндрических координат {x=rcosφy=rsinφz=h,r=(rcosφrsinφh) rr=(cosφsinφ0),rφ=(rsinφrcosφ0),rh=(001) Hr=cos2φ+sin2φ=1 Hφ=r2cos2φ+r2sin2φ=r Hh=0+0+1=1 тогда er=rrHr=(cosφsinφ0),eφ=rφHφ=(sinφcosφ0),eh=rhHh=(001). Причём, ereφ=(cosφsinφ0)(sinφcosφ0)=cosφsinφ+sinφcosφ=0,ereh=0,eheφ=0; er×eφ=|icosφsinφjsinφcosφk00|=k|cosφsinφsinφcosφ|=k(cos2φ(sin2φ))=k=(001)=eh, eφ×eh=eφ×[er×eφ]=er(eφeφ)eφ(eφer)=er1eφ0=er, eh×er=eh×[eφ×eh]=eφ(eheh)eh(eheφ)=eφ1eh0=eφ. Задание: найти коэффициенты Ламэ и орты сферических координат и проверить, что система сферических координат ортогональна.

Далее наша задача будет состоять в перезаписи в криволинейных координатах известных дифференциальных операторов: градиента, дивергенции и ротора. Так как те выражаются через производные по декартовым координатам, задача эта будет сводиться к замене переменных в дифференциальных выражениях (x,y,z)(p,q,s).

Градиент

Продифференцируем (1) по x, y и z: {1=xppx+xqqx+xssx0=yppx+yqqx+yssx0=zppx+zqqx+zssx{0=xppy+xqqy+xssy1=yppy+yqqy+yssy0=zppy+zqqy+zssy{0=xppz+xqqz+xssz0=yppz+yqqz+yssz1=zppz+zqqz+zssz По правилу Крамера, из первой системы px=|1xqxs0yqys0zqzs||xpxqxsypyqyszpzqzs|=(i,rq,rs)(rp,rq,rs)=HqHs(i,eq,es)HpHqHs(ep,eq,es)=1Hpi[eq×es]ep[eq×es]=1Hpiepepep=iepHp. Аналогично, из второй системы в (9) получаем py=|0xqxs1yqys0zqzs||xpxqxsypyqyszpzqzs|=jepHp,pz=|0xqxs0yqys1zqzs||xpxqxsypyqyszpzqzs|=kepHp. Так как для всякого вектора a=iax+jay+kaz=i(ia)+j(ja)+k(ka), gradp=ipx+jpy+kpz=iiepHp+jjepHp+kkepHp=epHp. Решая системы (9) относительно производных других координат, тем же путём находим, что gradq=eqHq,grads=esHs.

Градиент произвольной функции u(p,q,s) находится отсюда мгновенно: gradu=upgradp+uqgradq+usgrads=upHpep+uqHqeq+usHses Например, №189: Найти градиент в цилиндрических координатах поля (я тут дозаменил z на h) u=r2+2rcosφehsinφ ur=2r+2cosφ,uφ=2rsinφehcosφ,uh=ehsinφ gradu=urHrer+uφHφeφ+uhHheh=(2r+2cosφ)er+2rsinφehcosφreφehsinφeh. Задание: №190.

Вспомогательные формулы. Свойства производных e.

1) Возьмём свойство ортогональности ортов (извините) и продифференцируем по координате: γ|eαeβ=δαβ при этом сработает правило Лейбница eβγeα+eαγeβ=0 eβγeα=eαγeβ, где α, β и γ -- произвольные координаты из новой системы координат

2) В частности, когда α=β, eβγeβ=eβγeβ, Заменим ещё γ=α, и получим eβαeβ=0, 3) Смешанные производные не зависят от порядка дифференцирования, поэтому rαβ=rβα, подставим rβ и rα: β(Hαeα)=α(Hβeβ) Применим правило Лейбница и умножим обе части βHαeα+Hαβeα=αHβeβ+Hβαeβ|eβ В силу (20) Hαeββeα=αHβ, eββeα=1HααHβ, где αβ.

4)

Теперь умножим βHαeα+Hαβeα=αHβeβ+Hβαeβ|eγ (γα, γβ) Hαeγβeα=Hβeγαeβ eγβeα=HβHαeγαeβ Эти четыре формулы пригодятся ниже.

Дивергенция

Для поля F=iP+jQ+kR дивергенция примет вид divF=Px+Qy+Rz= =Pppx+Pqqx+Pssx+Qppy+Qqqy+Qssy+Rppz+Rqqz+Rssz= =(Pppx+Qppy+Rppz)+(Pqqx+Qqqy+Rqqz)+(Pssx+Qssy+Rssz)= =Fpgradp+Fqgradq+Fsgrads=FpepHp+FqeqHq+FsesHs

В локальном криволинейном базисе F=Fpep+Fqeq+Fses Fpep=p(Fpep+Fqeq+Fses)ep= =(pFpep+pFqeq+pFses+Fppep+Fqpeq+Fspes)ep= по (25) =pFp+Fqeppeq+Fseppes=pFp+Fq1HqqHp+Fs1HssHp= =1HqHs(HqHspFp+FqHsqHp+FsHqsHp). Так же получим, что Fqeq=1HpHs(HpHsqFq+FpHspHq+FsHpsHq), Fses=1HpHq(HpHqsFs+FpHqpHs+FqHpqHs). Тогда divF=FpepHp+FqeqHq+FsesHs= =1HqHsHp(HqHspFp+FqHsqHp+FsHqsHp)+ 1HpHsHq(HpHsqFq+FpHspHq+FsHpsHq)+ 1HpHqHs(HpHqsFs+FpHqpHs+FqHpqHs)= =1HpHqHs(HqHspFp+FpHspHq+FpHqpHs+ +HpHsqFq+FqHsqHp+FqHpqHs+ HpHqsFs+FsHqsHp+FsHpsHq)= =1HpHqHs[p(HqHsFp)+q(HpHsFq)+s(HpHqFs)] Например, в цилиндрических координатах №191 (я опять заменил z на h во избежание путаницы): F=φarctgrer+2eφh2eheh divF=1HrHφHh[r(HφHhFr)+φ(HrHhFφ)+h(HrHφFh)]= =1r[r(rFr)+φ(Fφ)+h(rFh)]= =1r[r(rφarctgr)+φ(2)+h(h2ehr)]= =1r[φ(arctgr+r11+r2)r(2heh+h2eh)]= =φrarctgr+φ1+r2(2+h)heh. Задание: №192.

Комментариев нет »

No comments yet.

RSS feed for comments on this post.

Leave a comment

Для отправки комментария вам необходимо авторизоваться.

Хостингом угостил Вадим "Moose" Калинников