<object id="bg47o"></object>

    
    

  1. 您好,歡迎來到山東合運電氣有限公司網(wǎng)站!

    關(guān)于合運 | 聯(lián)系我們 | 用戶須知 | sitemap

    400-088-6921155-8888-6921

    電源問答

    首頁 > 電源問答

    電勢

    時間:2022-10-16 人氣: 來源:山東合運電氣有限公司

      在這篇文章內(nèi),向量與標(biāo)量分別用粗體與斜體顯示。例如,位置向量通常用{\displaystyle\mathbf{r}\,\!}\mathbf{r}\,\!表示;而其大小則用{\displaystyle r\,\!}r\,\!來表示。Treatise_on_Electricity_and_Magnetism_Fig_01.jpg


      在靜電學(xué)里,電勢(electric potential)又稱電位[1],是描述電場中某一點之能量高低性質(zhì)的物理標(biāo)量,操作型定義為“電場中某處的電勢”等于“處于電場中該位置的單位電荷所具有的電勢能”[2],單位用伏特。


      電勢的數(shù)值不具有絕對意義,只具有相對意義,因此為了便于分析問題,必須設(shè)定一個參考位置,并把它設(shè)為零,稱為零勢能點。通常,會把無窮遠(yuǎn)處的電勢設(shè)定為零。那么,電勢可以定義如下:假設(shè)檢驗電荷從無窮遠(yuǎn)位置,經(jīng)過任意路徑,克服電場力,以緩慢、沒有產(chǎn)生加速度的方式移動到某位置,則在這位置的電勢,等于因移動檢驗電荷所做的功與檢驗電荷的電荷量的比值。在國際單位制里,電勢的單位為伏特({\displaystyle\scriptstyle{{\text{V}}={\text{J}}/{\text{C}}}}{\displaystyle\scriptstyle{{\text{V}}={\text{J}}/{\text{C}}}})(Volt),它是為了紀(jì)念意大利物理學(xué)家亞歷山德羅·伏打(Alessandro Volta)而命名。


      電勢必需滿足泊松方程,同時符合相關(guān)邊界條件;假設(shè)在某區(qū)域內(nèi)的電荷密度為零,則泊松方程約化為拉普拉斯方程,電勢必需滿足拉普拉斯方程。


      在電動力學(xué)里,當(dāng)含時電磁場存在的時候,電勢可以延伸為“廣義電勢”。特別注意,廣義電勢不能被視為電勢能每單位電荷。

    簡介


      處于外電場的帶電粒子會受到外電場施加的作用力,稱為電場力,促使帶電粒子加速運動。對于帶正電粒子,電場力與電場同方向;對于帶負(fù)電粒子,電場力與電場反方向。電場力的數(shù)值大小與電荷量、電場數(shù)值大小成正比。


      作用力與勢能之間有非常直接的關(guān)系。隨著物體朝著作用力的方向的加速運動,物體的動能變大,勢能變小。例如,一個石頭在山頂?shù)闹亓菽艽笥谠谏侥_的重力勢能。隨著物體的滾落,重力勢能變小,動能變大。


      對于某種特別作用力,科學(xué)家可以定義其向量場和其位勢,使得物體因為這向量場而具有的勢能,只與物體位置、參考位置之間的距離有關(guān)。稱這種作用力為保守力,這種向量場為保守場。


      例如,重力、靜電場的電場力,都是保守力。靜電場的標(biāo)勢稱為電勢,或稱為靜電勢。


      電勢和磁矢勢共同形成一個四維向量,稱為四維勢。從某一個慣性參考系觀察到的四維勢,應(yīng)用洛倫茲變換,可以計算出另外一個慣性參考系所觀察到的四維勢。


    靜電學(xué)里的電勢


      在靜電學(xué)里,電場{\displaystyle\mathbf{E}}\mathbf{E}內(nèi)某位置{\displaystyle\mathbf{r}}\mathbf{r}的電勢{\displaystyle\phi}\phi,以方程定義為[2]


      {\displaystyle\phi(\mathbf{r})\{\stackrel{def}{=}}\U_{\mathrm{E}}(\mathbf{r})/q}\phi(\mathbf{r})\\stackrel{def}{=}\U_\mathrm{E}(\mathbf{r})/q;


      其中,{\displaystyle U_{\mathrm{E}}}U_\mathrm{E}是在位置{\displaystyle\mathbf{r}}\mathbf{r}的檢驗電荷{\displaystyle q}q所具有的電勢能。


      電勢能的數(shù)值是人為設(shè)定的,沒有絕對意義,只有相對于某參考位置的已設(shè)定參考值時才有物理意義。假若要設(shè)定電勢能在空間任意位置的數(shù)值,必須先設(shè)定其在某參考位置{\displaystyle\mathbf{r}_{0}}\mathbf{r}_0的數(shù)值。為了方便運算,假設(shè)其參考數(shù)值為0。然后,就可以將在位置{\displaystyle\mathbf{r}}\mathbf{r}的電勢能{\displaystyle U_{\mathrm{E}}(\mathbf{r})}U_\mathrm{E}(\mathbf{r})定義為從參考位置{\displaystyle\mathbf{r}_{0}}\mathbf{r}_0緩慢地將檢驗電荷{\displaystyle q}q移動至{\displaystyle\mathbf{r}}\mathbf{r}所需做的機(jī)械功{\displaystyle W}W:


      {\displaystyle U_{\mathrm{E}}(\mathbf{r})\{\stackrel{def}{=}}\W}U_\mathrm{E}(\mathbf{r})\\stackrel{def}{=}\W。


      移動檢驗電荷時所施加的外力{\displaystyle\mathbf{F}}\mathbf{F},必須恰巧抵消處于電場{\displaystyle\mathbf{E}}\mathbf{E}的檢驗電荷{\displaystyle q}q所感受到的電場力{\displaystyle q\mathbf{E}}q\mathbf{E},即{\displaystyle\mathbf{F}=-q\mathbf{E}}\mathbf{F}=-q\mathbf{E}。其所做機(jī)械功等于外力{\displaystyle\mathbf{F}}\mathbf{F}的路徑積分:


      {\displaystyle W=\int _{\mathbb{L}}\mathbf{F}\cdot\mathrmhyxjpf1ac{\boldsymbol{\ell}}=-q\int _{\mathbb{L}}\mathbf{E}\cdot\mathrmhyxjpf1ac{\boldsymbol{\ell}}}W=\int_\mathbb{L}\mathbf{F}\cdot\mathrmhyxjpf1ac\boldsymbol{\ell}=-q\int_\mathbb{L}\mathbf{E}\cdot\mathrmhyxjpf1ac\boldsymbol{\ell};


      其中,{\displaystyle\mathbb{L}}\mathbb{L}是從參考位置{\displaystyle\mathbf{r}_{0}}\mathbf{r}_0到位置{\displaystyle\mathbf{r}}\mathbf{r}的一條任意路徑,{\displaystyle\mathrmhyxjpf1ac{\boldsymbol{\ell}}}\mathrmhyxjpf1ac{\boldsymbol{\ell}}是微小線元素。


      在靜電學(xué)里,{\displaystyle\mathbf{\nabla}\times\mathbf{E}=0}\mathbf{\nabla}\times\mathbf{E}=0,電場是保守場,所以,在積分時,可以選擇任意路徑{\displaystyle\mathbb{L}}\mathbb{L},計算出來的結(jié)果都一樣。欲知更詳盡細(xì)節(jié),請參閱條目保守力。由于這方程右邊的路徑積分跟路徑{\displaystyle\mathbb{L}}\mathbb{L}無關(guān),只跟路徑的初始位置{\displaystyle\mathbf{r}_{0}}\mathbf{r}_0、終止位置{\displaystyle\mathbf{r}}\mathbf{r}有關(guān),因此若能夠假設(shè)無窮遠(yuǎn)位置{\displaystyle\infty}\infty的電勢能為0,則可以設(shè)定參考位置{\displaystyle\mathbf{r}_{0}}\mathbf{r}_0在無窮遠(yuǎn)位置{\displaystyle\infty}\infty:


      {\displaystyle U_{\mathrm{E}}(\mathbf{r})=-q\int _{\infty}^{\mathbf{r}}\mathbf{E}\cdot\mathrmhyxjpf1ac{\boldsymbol{\ell}}}U_\mathrm{E}(\mathbf{r})=-q\int_\infty^\mathbf{r}\mathbf{E}\cdot\mathrmhyxjpf1ac\boldsymbol{\ell}。


      所以,電勢就是從無窮遠(yuǎn)位置到檢驗位置對于電場做路徑積分所得結(jié)果的負(fù)值:


      {\displaystyle\phi(\mathbf{r})=-\int _{\infty}^{\mathbf{r}}\mathbf{E}\cdot\mathrmhyxjpf1ac{\boldsymbol{\ell}}}\phi(\mathbf{r})=-\int_\infty^\mathbf{r}\mathbf{E}\cdot\mathrmhyxjpf1ac\boldsymbol{\ell}。


      在任意兩個位置{\displaystyle\mathbf{r}_{1}}\mathbf{r}_1、{\displaystyle\mathbf{r}_{2}}\mathbf{r}_2之間的“電勢差”{\displaystyle\Delta\phi}\Delta\phi為


      {\displaystyle\Delta\phi=\phi(\mathbf{r}_{2})-\phi(\mathbf{r}_{1})=-\int _{\mathbf{r}_{1}}^{\mathbf{r}_{2}}\mathbf{E}\cdot\mathrmhyxjpf1ac{\boldsymbol{\ell}}}\Delta\phi=\phi(\mathbf{r}_2)-\phi(\mathbf{r}_1)=-\int_{\mathbf{r}_1}^{\mathbf{r}_2}\mathbf{E}\cdot\mathrmhyxjpf1ac\boldsymbol{\ell}。


      由于電場{\displaystyle\mathbf{E}}\mathbf{E}是保守場,電勢差也與積分路徑無關(guān),只跟積分路徑的初始位置與終止位置有關(guān)。


    點電荷


      由點電荷Q所產(chǎn)生的電勢,在距離r時,可表示為


      {\displaystyle V={\frac{1}{4\pi\varepsilon _{0}}}{\frac{Q}{r}}}V=\frac{1}{4\pi\varepsilon_0}\frac{Q}{r}


      其中,ε0是真空電容率。


      在無限遠(yuǎn)處,電勢為零。由多個點電荷產(chǎn)生的電勢,相等于各點電荷所產(chǎn)生的電勢之和。此外,電勢場是標(biāo)量場,電場則是向量場。


    疊加原理


      電場遵守疊加原理:假設(shè)在三維空間里,由兩組完全不相交的電荷分布所產(chǎn)生的電場分別為{\displaystyle\mathbf{E}_{1}}\mathbf{E}_1、{\displaystyle\mathbf{E}_{2}}\mathbf{E}_2,則總電場為{\displaystyle\mathbf{E}_{t}=\mathbf{E}_{1}+\mathbf{E}_{2}}\mathbf{E}_t=\mathbf{E}_1+\mathbf{E}_2。


      總電勢為每單位電荷克服電場力所做的機(jī)械功之和:


      {\displaystyle\phi _{t}(\mathbf{r})=-\int _{\infty}^{\mathbf{r}}\mathbf{E}_{t}\cdot\mathrmhyxjpf1ac{\boldsymbol{\ell}}=-\int _{\infty}^{\mathbf{r}}(\mathbf{E}_{1}+\mathbf{E}_{2})\cdot\mathrmhyxjpf1ac{\boldsymbol{\ell}}=\phi _{1}(\mathbf{r})+\phi _{2}(\mathbf{r})}\phi_t(\mathbf{r})=-\int_\infty^\mathbf{r}\mathbf{E}_t\cdot\mathrmhyxjpf1ac\boldsymbol{\ell}=-\int_\infty^\mathbf{r}(\mathbf{E}_1+\mathbf{E}_2)\cdot\mathrmhyxjpf1ac\boldsymbol{\ell}=\phi_1(\mathbf{r})+\phi_2(\mathbf{r})。


      所以,電勢也遵守疊加原理。當(dāng)計算一組電荷分布所產(chǎn)生的電勢時,只需要知道在電荷分布的每個源位置的單獨電荷所產(chǎn)生在檢驗位置的電勢,就可以應(yīng)用積分運算,得到整個電荷分布所產(chǎn)生在檢驗位置的電勢。


    電勢的微分方程


      應(yīng)用積分符號內(nèi)取微分方法,電勢的梯度為


      {\displaystyle\mathbf{\nabla}\phi(\mathbf{r})=-\mathbf{\nabla}\int _{\infty}^{\mathbf{r}}\mathbf{E}(\mathbf{r}')\cdot\mathrmhyxjpf1ac{\boldsymbol{\ell}}^{\,\prime}=-\mathbf{E}(\mathbf{r})}\mathbf{\nabla}\phi(\mathbf{r})=-\mathbf{\nabla}\int_\infty^\mathbf{r}\mathbf{E}(\mathbf{r}')\cdot\mathrmhyxjpf1ac\boldsymbol{\ell}^{\,\prime}=-\mathbf{E}(\mathbf{r})。


      所以,電場與電勢之間的關(guān)系為


      {\displaystyle\mathbf{E}(\mathbf{r})=-\mathbf{\nabla}\phi(\mathbf{r})}\mathbf{E}(\mathbf{r})=-\mathbf{\nabla}\phi(\mathbf{r})。


      根據(jù)高斯定律的方程,


      {\displaystyle\mathbf{\nabla}\cdot\mathbf{E}=\rho/\epsilon _{0}}\mathbf{\nabla}\cdot\mathbf{E}=\rho/\epsilon_0;


      其中,{\displaystyle\rho}\rho是電荷密度,{\displaystyle\epsilon _{0}}\epsilon _{0}是電常數(shù)。


      所以,電勢滿足泊松方程:


      {\displaystyle\nabla^{2}\phi=-\rho/\epsilon _{0}}\nabla^2\phi=-\rho/\epsilon_0。


      假設(shè)電荷密度為零,則這方程變?yōu)槔绽狗匠蹋?/span>


      {\displaystyle\nabla^{2}\phi=0}\nabla^2\phi=0。


      請注意,假若{\displaystyle\mathbf{\nabla}\times\mathbf{E}\neq 0}\mathbf{\nabla}\times\mathbf{E}\ne 0,也就是說,電場不具保守性(由于隨時間變化的磁場造成的效應(yīng);參閱麥克斯韋方程組),則不能使用這些方程。


      由于電勢乃是標(biāo)量,而電場是具有三個分量的向量,所以,很多時候,使用電勢來解析問題會省去很多運算工作,帶來很大的便利。


      拉普拉斯方程的解答


      在某空間區(qū)域內(nèi),假設(shè)電荷密度為零,則電勢必須滿足拉普拉斯方程,并且符合所有相關(guān)邊界條件。


    邊界條件


      在靜電學(xué)里,有三種邊界條件:


      狄利克雷邊界條件:在所有邊界,電勢都已良態(tài)給定。具有這種邊界條件的問題稱為狄利克雷問題。


      紐曼邊界條件:在所有邊界,電勢的法向?qū)?shù)都已良態(tài)給定。具有這種邊界條件的問題稱為紐曼問題。


      混合邊界條件:一部分邊界的電勢都已良態(tài)給定,其它邊界的電勢的法向?qū)?shù)也已良態(tài)給定。


      根據(jù)拉普拉斯方程的唯一性定理,對于這些種類的邊界條件,拉普拉斯方程的解答都具有唯一性。所以,只要找到一個符合邊界條件的解答,則這解答必定為正確解答。


    分離變數(shù)法


      應(yīng)用分離變數(shù)法來解析拉普拉斯方程,可以將問題的偏微分方程改變?yōu)橐唤M較容易解析的常微分方程。對于一般問題,通常會采用直角坐標(biāo)系、圓柱坐標(biāo)系或球坐標(biāo)系來分離拉普拉斯方程。但是,對于其它比較特別的問題,另外還有八種坐標(biāo)系可以用來分離拉普拉斯方程。[3]分離之后,找到每一個常微分方程的通解(通常為一組本征方程的疊加),電勢可以表達(dá)為這些通解的乘積。將這表達(dá)式與邊界條件相匹配,就可以設(shè)定一般解的系數(shù),從而找到問題的特解。根據(jù)拉普拉斯方程的唯一性定理,這特解也是唯一的正確解答。


    兩個半平面導(dǎo)體案例

    500px-Two_half_planes_with_different_voltages.svg.png

      被位于{\displaystyle y=0}y=0的絕緣線條分隔為處于y+、y--半平面的兩個導(dǎo)體的電勢分別設(shè)定為{\displaystyle+V}+V、{\displaystyle-V}-V。


      假設(shè)在xy-平面的無限平面導(dǎo)體被一條位于{\displaystyle y=0}y=0的絕緣線條分為兩半,兩個處于y+、y--半平面的導(dǎo)體的電勢分別設(shè)定為{\displaystyle+V}+V、{\displaystyle-V}-V,則計算z+-半空間任意位置的電勢這問題,由于邊界條件的幾何形狀適合用直角坐標(biāo)來描述,可以以直角坐標(biāo){\displaystyle(x,y,z)}(x,y,z)將拉普拉斯方程表示為:


      {\displaystyle\nabla^{2}\phi={\frac{\partial^{2}\phi}{\partial x^{2}}}+{\frac{\partial^{2}\phi}{\partial y^{2}}}+{\frac{\partial^{2}\phi}{\partial z^{2}}}=0}\nabla^2\phi=\frac{\partial^2\phi}{\partial x^2}+\frac{\partial^2\phi}{\partial y^2}+\frac{\partial^2\phi}{\partial z^2}=0。


      因為這案例與x-坐標(biāo)無關(guān),方程可以簡化為


      {\displaystyle\nabla^{2}\phi(y,z)={\frac{\partial^{2}\phi}{\partial y^{2}}}+{\frac{\partial^{2}\phi}{\partial z^{2}}}=0}\nabla^2\phi(y,z)=\frac{\partial^2\phi}{\partial y^2}+\frac{\partial^2\phi}{\partial z^2}=0。


      應(yīng)用分離變數(shù)法,猜想解答的形式為


      {\displaystyle\phi(y,z)=Y(y)Z(z)}\phi(y,z)=Y(y)Z(z)。


      將這公式代入拉普拉斯方程,則可得到


      {\displaystyle{\frac{1}{Y(y)}}\{\frac{\mathrmhyxjpf1ac^{2}Y(y)}{\mathrmhyxjpf1acy^{2}}}+{\frac{1}{Z(z)}}\{\frac{\mathrmhyxjpf1ac^{2}Z(z)}{\mathrmhyxjpf1acz^{2}}}=0}\frac{1}{Y(y)}\\frac{\mathrmhyxjpf1ac^2 Y(y)}{\mathrmhyxjpf1acy^2}+\frac{1}{Z(z)}\\frac{\mathrmhyxjpf1ac^2 Z(z)}{\mathrmhyxjpf1acz^2}=0。


      注意到這方程的每一個項目都只含有一個變量,并且跟其它變量無關(guān)。所以,每一個項目都等于常數(shù):


      {\displaystyle{\frac{1}{Y(y)}}\{\frac{\mathrmhyxjpf1ac^{2}Y(y)}{\mathrmhyxjpf1acy^{2}}}=C}\frac{1}{Y(y)}\\frac{\mathrmhyxjpf1ac^2 Y(y)}{\mathrmhyxjpf1acy^2}=C、


      {\displaystyle{\frac{1}{Z(z)}}\{\frac{\mathrmhyxjpf1ac^{2}Z(z)}{\mathrmhyxjpf1acz^{2}}}=-C}\frac{1}{Z(z)}\\frac{\mathrmhyxjpf1ac^2 Z(z)}{\mathrmhyxjpf1acz^2}=-C。


      這樣,一個二次偏微分方程被改變?yōu)閮蓚€簡單的二次常微分方程。解答分別為


      {\displaystyle Y(y)=A_{1}e^{iky}+A_{2}e^{-iky}}Y(y)=A_1 e^{iky}+A_2 e^{-iky}、


      {\displaystyle Z(z)=B_{1}e^{kz}+B_{2}e^{-kz}}Z(z)=B_1 e^{kz}+B_2 e^{-kz};


      其中,{\displaystyle A_{1}(k)}A_1(k)、{\displaystyle A_{2}(k)}A_2(k)、{\displaystyle B_{1}(k)}B_1(k)、{\displaystyle B_{2}(k)}B_2(k)都是系數(shù)函數(shù)。


      當(dāng){\displaystyle z}z趨向于無窮大時,{\displaystyle Z(z)}Z(z)趨向于零,所以,{\displaystyle B_{1}=0}B_1=0。綜合起來,電勢為


      {\displaystyle\phi(y,z)=\int _{0}^{\infty}(A_{1}e^{iky}+A_{2}e^{-iky})e^{-kz}\mathrmhyxjpf1ack}\phi(y,z)=\int_0^{\infty}(A_1 e^{iky}+A_2 e^{-iky})e^{-kz}\mathrmhyxjpf1ack。


      由于在{\displaystyle z=0}z=0,y+、y--半平面的電勢分別為{\displaystyle+V}+V、{\displaystyle-V}-V,所以,


      當(dāng){\displaystyle y&gt;0}y&gt;0時,{\displaystyle\int _{0}^{\infty}(A_{1}e^{iky}+A_{2}e^{-iky})\mathrmhyxjpf1ack=+V}\int_0^{\infty}(A_1 e^{iky}+A_2 e^{-iky})\mathrmhyxjpf1ack=+V、


      當(dāng){\displaystyle y&lt;0}y&lt;0時,{\displaystyle\int _{0}^{\infty}(A_{1}e^{iky}+A_{2}e^{-iky})\mathrmhyxjpf1ack=-V}\int_0^{\infty}(A_1 e^{iky}+A_2 e^{-iky})\mathrmhyxjpf1ack=-V。


      應(yīng)用傅里葉變換,可以得到


      {\displaystyle A_{1}(k)={\frac{V}{2\pi}}\left(\int _{0}^{\infty}e^{-iky'}\mathrmhyxjpf1acy'-\int _{-\infty}^{0}e^{-iky'}\mathrmhyxjpf1acy'\right)}A_1(k)=\frac{V}{2\pi}\left(\int_0^{\infty}e^{-iky'}\mathrmhyxjpf1acy'-\int_{-\infty}^0 e^{-iky'}\mathrmhyxjpf1acy'\right)、


      {\displaystyle A_{2}(k)={\frac{V}{2\pi}}\left(\int _{0}^{\infty}e^{iky'}\mathrmhyxjpf1acy'-\int _{-\infty}^{0}e^{iky'}\mathrmhyxjpf1acy'\right)}A_2(k)=\frac{V}{2\pi}\left(\int_0^{\infty}e^{iky'}\mathrmhyxjpf1acy'-\int_{-\infty}^0 e^{iky'}\mathrmhyxjpf1acy'\right)。


      所以,由{\displaystyle A_{1}(k)}A_1(k)項目貢獻(xiàn)出的電勢為


      {\displaystyle{\begin{aligned}\phi _{1}&={\frac{V}{2\pi}}\int _{0}^{\infty}\mathrmhyxjpf1ack\left\{\int _{0}^{\infty}e^{ik(y-y')-kz}\mathrmhyxjpf1acy'-\int _{-\infty}^{0}e^{ik(y-y')-kz}\mathrmhyxjpf1acy'\right\}\\&=-\{\frac{V}{2\pi}}\int _{0}^{\infty}{\frac{\mathrmhyxjpf1acy'}{i(y-y')-z}}+\{\frac{V}{2\pi}}\int _{-\infty}^{0}{\frac{\mathrmhyxjpf1acy'}{i(y-y')-z}}\\\end{aligned}}}\begin{align}\phi_1&=\frac{V}{2\pi}\int_0^{\infty}\mathrmhyxjpf1ack\left\{\int_0^{\infty}e^{ik(y-y')-kz}\mathrmhyxjpf1acy'-\int_{-\infty}^0e^{ik(y-y')-kz}\mathrmhyxjpf1acy'\right\}\\


      &=-\\frac{V}{2\pi}\int_0^{\infty}\frac{\mathrmhyxjpf1acy'}{i(y-y')-z}+\\frac{V}{2\pi}\int_{-\infty}^0\frac{\mathrmhyxjpf1acy'}{i(y-y')-z}\\


      \end{align}。


      類似地,由{\displaystyle A_{2}(k)}A_2(k)項目貢獻(xiàn)出的電勢為


      {\displaystyle{\begin{aligned}\phi _{2}&={\frac{V}{2\pi}}\int _{0}^{\infty}\mathrmhyxjpf1ack\left\{\int _{0}^{\infty}e^{-ik(y-y')-kz}\mathrmhyxjpf1acy'-\int _{-\infty}^{0}e^{-ik(y-y')-kz}\mathrmhyxjpf1acy'\right\}\\&=-\{\frac{V}{2\pi}}\int _{0}^{\infty}{\frac{\mathrmhyxjpf1acy'}{-i(y-y')-z}}+\{\frac{V}{2\pi}}\int _{-\infty}^{0}{\frac{\mathrmhyxjpf1acy'}{-i(y-y')-z}}\\\end{aligned}}}\begin{align}\phi_2&=\frac{V}{2\pi}\int_0^{\infty}\mathrmhyxjpf1ack\left\{\int_0^{\infty}e^{-ik(y-y')-kz}\mathrmhyxjpf1acy'-\int_{-\infty}^0e^{-ik(y-y')-kz}\mathrmhyxjpf1acy'\right\}\\


      &=-\\frac{V}{2\pi}\int_0^{\infty}\frac{\mathrmhyxjpf1acy'}{-i(y-y')-z}+\\frac{V}{2\pi}\int_{-\infty}^0\frac{\mathrmhyxjpf1acy'}{-i(y-y')-z}\\


      \end{align}。


      總電勢為[4]


      {\displaystyle{\begin{aligned}\phi&={\frac{Vz}{\pi}}\int _{0}^{\infty}{\frac{\mathrmhyxjpf1acy'}{(y-y')^{2}+z^{2}}}-\{\frac{Vz}{\pi}}\int _{-\infty}^{0}{\frac{\mathrmhyxjpf1acy'}{(y-y')^{2}+z^{2}}}\\&={\frac{2V}{\pi}}\\arctan{\left({\frac{y}{z}}\right)}\\\end{aligned}}}\begin{align}\phi&=\frac{Vz}{\pi}\int_0^{\infty}\frac{\mathrmhyxjpf1acy'}{(y-y')^2+z^2}-\\frac{Vz}{\pi}\int_{-\infty}^0\frac{\mathrmhyxjpf1acy'}{(y-y')^2+z^2}\\


      &=\frac{2V}{\pi}\\arctan{\left(\frac{y}{z}\right)}\\


      \end{align}。

    泊松方程的解答


    電荷分布所產(chǎn)生的電勢


      根據(jù)庫侖定律,一個源位置為{\displaystyle\mathbf{r}'}\mathbf{r}'的點電荷{\displaystyle q}q,所產(chǎn)生在任意位置{\displaystyle\mathbf{r}}\mathbf{r}的電場為


      {\displaystyle\mathbf{E}(\mathbf{r})={\frac{q}{4\pi\epsilon _{0}}}\{\frac{(\mathbf{r}-\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|^{3}}}}\mathbf{E}(\mathbf{r})=\frac{q}{4\pi\epsilon_0}\\frac{(\mathbf{r}-\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|^3}。


      對于一群點電荷,應(yīng)用疊加原理,總電場等于每一個點電荷所產(chǎn)生的電場的疊加。體積區(qū)域{\displaystyle\mathbb{V}'}\mathbb{V}'內(nèi)部電荷密度為{\displaystyle\rho(\mathbf{r}')}\rho(\mathbf{r}')的電荷分布,在檢驗位置{\displaystyle\mathbf{r}}\mathbf{r}所產(chǎn)生的電場為


      {\displaystyle\mathbf{E}(\mathbf{r})={\frac{1}{4\pi\epsilon _{0}}}\int _{\mathbb{V}'}\rho(\mathbf{r}'){\frac{(\mathbf{r}-\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|^{3}}}\\mathrmhyxjpf1ac^{3}r'}\mathbf{E}(\mathbf{r})=\frac{1}{4\pi\epsilon_0}\int_{\mathbb{V}'}\rho(\mathbf{r}')\frac{(\mathbf{r}-\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|^3}\\mathrmhyxjpf1ac^3 r';


      其中,{\displaystyle\mathrmhyxjpf1ac^{3}r'}\mathrmhyxjpf1ac^3 r'是微小體積元素。


      應(yīng)用一條向量恒等式,


      {\displaystyle\nabla{\frac{1}{|\mathbf{r}-\mathbf{r}'|}}=-\{\frac{(\mathbf{r}-\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|^{3}}}}\nabla\frac{1}{|\mathbf{r}-\mathbf{r}'|}=-\\frac{(\mathbf{r}-\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|^3},


      可以得到


      {\displaystyle\mathbf{E}(\mathbf{r})=-\{\frac{1}{4\pi\epsilon _{0}}}\nabla\int _{\mathbb{V}'}{\frac{\rho(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}}\\mathrmhyxjpf1ac^{3}r'}\mathbf{E}(\mathbf{r})=-\\frac{1}{4\pi\epsilon_0}\nabla\int_{\mathbb{V}'}\frac{\rho(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}\\mathrmhyxjpf1ac^3 r'。


      設(shè)定在無窮遠(yuǎn)的電勢為參考值0,則在任意位置的電勢為


      {\displaystyle\phi(\mathbf{r})={\frac{1}{4\pi\epsilon _{0}}}\int _{\mathbb{V}'}{\frac{\rho(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}}\\mathrmhyxjpf1ac^{3}r'}\phi(\mathbf{r})=\frac{1}{4\pi\epsilon_0}\int_{\mathbb{V}'}\frac{\rho(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}\\mathrmhyxjpf1ac^3 r';(1)


      應(yīng)用一則關(guān)于狄拉克δ函數(shù)的向量恒等式


      {\displaystyle\nabla^{2}\left({\frac{1}{|\mathbf{r}-\mathbf{r}'|}}\right)=-4\pi\delta(\mathbf{r}-\mathbf{r}')}\nabla^2\left(\frac{1}{|\mathbf{r}-\mathbf{r}'|}\right)


      =-4\pi\delta(\mathbf{r}-\mathbf{r}'),


      假設(shè)檢驗位置{\displaystyle\mathbf{r}}\mathbf{r}在積分體積{\displaystyle\mathbb{V}'}\mathbb{V}'內(nèi),則可得到泊松方程:


      {\displaystyle\nabla^{2}\phi(\mathbf{r})={\frac{1}{4\pi\epsilon _{0}}}\int _{\mathbb{V}'}\nabla^{2}\left({\frac{\rho(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}}\right)\\mathrmhyxjpf1ac^{3}r'=-\{\frac{1}{\epsilon _{0}}}\int _{\mathbb{V}'}\rho(\mathbf{r}')\delta(\mathbf{r}-\mathbf{r}')\\mathrmhyxjpf1ac^{3}r'=-\{\frac{\rho(\mathbf{r})}{\epsilon _{0}}}}\nabla^2\phi(\mathbf{r})=\frac{1}{4\pi\epsilon_0}\int_{\mathbb{V}'}\nabla^2\left(\frac{\rho(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}\right)


      \\mathrmhyxjpf1ac^3 r'=-\\frac{1}{\epsilon_0}\int_{\mathbb{V}'}\rho(\mathbf{r}')\delta(\mathbf{r}-\mathbf{r}')\\mathrmhyxjpf1ac^3 r'


      =-\\frac{\rho(\mathbf{r})}{\epsilon_0}。


      所以,電勢的方程(1)為泊松方程的解答。


    邊界條件


      電勢的方程(1)只考慮到一群電荷分布所產(chǎn)生的電勢。假若遭遇邊界條件為電勢的靜電學(xué)問題,就不能使用方程(1),必需使用更具功能的方法。


      根據(jù)格林第二恒等式,對于任意良態(tài)函數(shù){\displaystyle\phi(\mathbf{r})}\phi(\mathbf{r})與{\displaystyle\psi(\mathbf{r})}\psi(\mathbf{r}),[5]


      {\displaystyle\int _{\mathbb{V}}\left(\phi\nabla^{2}\psi-\psi\nabla^{2}\phi\right)\\mathrmhyxjpf1ac^{3}r=\oint _{\mathbb{S}}\left(\phi{\partial\psi\over\partial n}-\psi{\partial\phi\over\partial n}\right)\\mathrmhyxjpf1ac^{2}r}\int_{\mathbb{V}}\left(\phi\nabla^2\psi-\psi\nabla^2\phi\right)\\mathrmhyxjpf1ac^3 r=\oint_{\mathbb{S}}\left(\phi{\partial\psi\over\partial n}-\psi{\partial\phi\over\partial n}\right)\\mathrmhyxjpf1ac^2 r;


      其中,{\displaystyle\mathbb{V}}\mathbb{V}是積分體積,{\displaystyle\mathbb{S}}\mathbb{S}是包住{\displaystyle\mathbb{V}}\mathbb{V}的閉表面,{\displaystyle\mathrmhyxjpf1ac^{2}r}\mathrmhyxjpf1ac^2 r是微小面元素,{\displaystyle\partial\phi\over\partial n}\partial\phi\over\partial n或{\displaystyle\partial\phi\over\partial n}\partial\phi\over\partial n都是取垂直于閉表面{\displaystyle\mathbb{S}}\mathbb{S}的法向?qū)?shù),都是從積分體積{\displaystyle\mathbb{V}}\mathbb{V}朝外指出。


      設(shè)定{\displaystyle\phi(\mathbf{r}')}\phi(\mathbf{r}')為在{\displaystyle\mathbf{r}'}\mathbf{r}'的電勢,{\displaystyle\psi={\frac{1}{|\mathbf{r}-\mathbf{r}'|}}}\psi=\frac{1}{|\mathbf{r}-\mathbf{r}'|}為{\displaystyle\mathbf{r}'}\mathbf{r}'與{\displaystyle\mathbf{r}}\mathbf{r}之間的距離。應(yīng)用泊松方程{\displaystyle\nabla^{2}\phi(\mathbf{r})=-\rho/\epsilon _{0}}\nabla^2\phi(\mathbf{r})=-\rho/\epsilon_0,則可得到


      {\displaystyle\int _{\mathbb{V}'}\left[\phi(\mathbf{r}')\nabla^{2}\left({\frac{1}{|\mathbf{r}-\mathbf{r}'|}}\right)+{\frac{\rho(\mathbf{r}')}{\epsilon _{0}|\mathbf{r}-\mathbf{r}'|}}\right]\mathrmhyxjpf1ac^{3}r'=\oint _{\mathbb{S}'}\left[\phi\{\partial\over\partial n'}\left({\frac{1}{|\mathbf{r}-\mathbf{r}'|}}\right)-\left({\frac{1}{|\mathbf{r}-\mathbf{r}'|}}\right){\partial\phi\over\partial n'}\right]\mathrmhyxjpf1ac^{2}r'}\int_{\mathbb{V}'}\left[\phi(\mathbf{r}')\nabla^2\left(\frac{1}{|\mathbf{r}-\mathbf{r}'|}\right)+\frac{\rho(\mathbf{r}')}{\epsilon_0|\mathbf{r}-\mathbf{r}'|}\right]\mathrmhyxjpf1ac^3 r'=\oint_{\mathbb{S}'}\left[\phi\{\partial\over\partial n'}\left(\frac{1}{|\mathbf{r}-\mathbf{r}'|}\right)-\left(\frac{1}{|\mathbf{r}-\mathbf{r}'|}\right){\partial\phi\over\partial n'}\right]\mathrmhyxjpf1ac^2 r'。


      再應(yīng)用向量恒等式


      {\displaystyle\nabla^{2}\left({\frac{1}{|\mathbf{r}-\mathbf{r}'|}}\right)=-4\pi\delta(\mathbf{r}-\mathbf{r}')}\nabla^2\left(\frac{1}{|\mathbf{r}-\mathbf{r}'|}\right)=-4\pi\delta(\mathbf{r}-\mathbf{r}')。


      假設(shè)檢驗位置{\displaystyle\mathbf{r}}\mathbf{r}在積分體積{\displaystyle\mathbb{V}'}\mathbb{V}'內(nèi),則可得到


      {\displaystyle\phi(\mathbf{r})={\frac{1}{4\pi\epsilon _{0}}}\int _{\mathbb{V}'}{\frac{\rho(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}}\\mathrmhyxjpf1ac^{3}r'+{\frac{1}{4\pi}}\oint _{\mathbb{S}'}\left[\left({\frac{1}{|\mathbf{r}-\mathbf{r}'|}}\right){\partial\phi\over\partial n'}-\phi\{\partial\over\partial n'}\left({\frac{1}{|\mathbf{r}-\mathbf{r}'|}}\right)\right]\mathrmhyxjpf1ac^{2}r'}\phi(\mathbf{r})=\frac{1}{4\pi\epsilon_0}\int_{\mathbb{V}'}\frac{\rho(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}\\mathrmhyxjpf1ac^3 r'+\frac{1}{4\pi}\oint_{\mathbb{S}'}\left[\left(\frac{1}{|\mathbf{r}-\mathbf{r}'|}\right){\partial\phi\over\partial n'}-\phi\{\partial\over\partial n'}\left(\frac{1}{|\mathbf{r}-\mathbf{r}'|}\right)\right]\mathrmhyxjpf1ac^2 r'。


      這方程右手邊的體積分就是電勢的方程(1),而面積分就是因為邊界條件而添加的項目。這是{\displaystyle\mathbb{V}'}\mathbb{V}'體內(nèi)與體外之間的邊界曲面。面積分的第一個項目要求給定在邊界曲面的法向電場,即{\displaystyle E_{n'}=-{\partial\phi\over\partial n'}}E_{n'}=-{\partial\phi\over\partial n'},也就是面感應(yīng)電荷密度{\displaystyle\sigma=\epsilon _{0}E_{n'}}\sigma=\epsilon_0 E_{n'}。面積分的第二個項目要求給定在邊界曲面的電勢{\displaystyle\phi}\phi。假若能夠知道積分體積內(nèi)的電荷密度、在閉曲面的面電荷密度與電勢,就可以計算出在積分體積內(nèi)任意位置的電勢。


      根據(jù)柯西邊界條件,有時候,給定在邊界曲面的法向電場與電勢,可能會因為給定過多邊界條件,而造成無法計算出一致的電勢的狀況。實際而言,只要給定法向電場或電勢,兩者之一,就可以計算出電勢。[5]


      假若積分體積為無窮大空間,當(dāng){\displaystyle r'}r'趨向于無窮大時,則面積分的被積分項目會以{\displaystyle 1/r'^{3}}1/r'^3速率遞減,而積分面積會以{\displaystyle r'^{2}}r'^2速率遞增,所以,面積分項目會趨向于零,這方程約化為先前的電勢方程(1)。


    格林函數(shù)


      包括函數(shù){\displaystyle 1/|\mathbf{r}-\mathbf{r}'|}1/|\mathbf{r}-\mathbf{r}'|在內(nèi),有一類函數(shù){\displaystyle G(\mathbf{r},\mathbf{r}')}G(\mathbf{r},\mathbf{r}'),稱為格林函數(shù),能夠滿足方程


      {\displaystyle\nabla^{2}G(\mathbf{r},\mathbf{r}')=-4\pi\delta(\mathbf{r}-\mathbf{r}')}\nabla^2 G(\mathbf{r},\mathbf{r}')=-4\pi\delta(\mathbf{r}-\mathbf{r}')。


      另外,假設(shè)函數(shù){\displaystyle H(\mathbf{r},\mathbf{r}')}H(\mathbf{r},\mathbf{r}')滿足拉普拉斯方程


      {\displaystyle\nabla^{2}H(\mathbf{r},\mathbf{r}')=0}\nabla^2 H(\mathbf{r},\mathbf{r}')=0,


      則函數(shù){\displaystyle G'(\mathbf{r},\mathbf{r}')=G(\mathbf{r},\mathbf{r}')+H(\mathbf{r},\mathbf{r}')}G'(\mathbf{r},\mathbf{r}')=G(\mathbf{r},\mathbf{r}')+H(\mathbf{r},\mathbf{r}')也是格林函數(shù)。


      應(yīng)用這靈活性質(zhì),可以更嚴(yán)格地規(guī)定格林函數(shù):[5]


      對于狄利克雷問題,當(dāng)源位置{\displaystyle\mathbf{r}'}\mathbf{r}'在邊界表面{\displaystyle{\mathbb{S}'}}{\mathbb{S}'}時,規(guī)定格林函數(shù){\displaystyle G_{D}(\mathbf{r},\mathbf{r}')=0}G_D(\mathbf{r},\mathbf{r}')=0。這樣,從格林第二恒等式,設(shè)定{\displaystyle\phi(\mathbf{r}')}\phi(\mathbf{r}')為在{\displaystyle\mathbf{r}'}\mathbf{r}'的電勢,{\displaystyle\psi(\mathbf{r},\mathbf{r}')=G_{D}(\mathbf{r},\mathbf{r}')}\psi(\mathbf{r},\mathbf{r}')=G_D(\mathbf{r},\mathbf{r}'),則可得到


      {\displaystyle\phi(\mathbf{r})={\frac{1}{4\pi\epsilon _{0}}}\int _{\mathbb{V}'}\rho(\mathbf{r}')G_{D}(\mathbf{r},\mathbf{r}')\\mathrmhyxjpf1ac^{3}r'-\{\frac{1}{4\pi}}\oint _{\mathbb{S}'}\phi(\mathbf{r}')\{\partial G_{D}(\mathbf{r},\mathbf{r}')\over\partial n'}\mathrmhyxjpf1ac^{2}r'}\phi(\mathbf{r})=\frac{1}{4\pi\epsilon_0}\int_{\mathbb{V}'}\rho(\mathbf{r}')G_D(\mathbf{r},\mathbf{r}')\\mathrmhyxjpf1ac^3 r'


      -\\frac{1}{4\pi}\oint_{\mathbb{S}'}\phi(\mathbf{r}')\{\partial G_D(\mathbf{r},\mathbf{r}')\over\partial n'}\mathrmhyxjpf1ac^2 r'。(2)


      對于滿足紐曼問題,當(dāng)源位置{\displaystyle\mathbf{r}'}\mathbf{r}'在邊界表面{\displaystyle{\mathbb{S}'}}{\mathbb{S}'}時,規(guī)定格林函數(shù){\displaystyle\oint _{\mathbb{S}'}{\frac{\partial G_{D}(\mathbf{r},\mathbf{r}')}{\partial n'}}\mathrmhyxjpf1ac^{2}r'=-{\frac{4\pi}{S}}}\oint_{\mathbb{S}'}\frac{\partial G_D(\mathbf{r},\mathbf{r}')}{\partial n'}\mathrmhyxjpf1ac^2 r'=-\frac{4\pi}{S}。


      這兩種規(guī)定都能夠唯一地設(shè)定格林函數(shù)。注意到格林函數(shù)是一個幾何函數(shù),與整個系統(tǒng)的電荷分布無關(guān)。對于任何系統(tǒng),只要計算出適合其幾何形狀的格林函數(shù),則不論系統(tǒng)的電荷分布為何,都可以使用同樣的格林函數(shù)。


    無限平面導(dǎo)體案例

    400px-Method_Of_Images_1_electric_charge.svg.png

      位于xy-平面的是一個接地的無限平面導(dǎo)體。其上方的點電荷{\displaystyle q}q的直角坐標(biāo)是{\displaystyle(0,\,0,\,a)}(0,\,0,\,a)。


      假設(shè)xy-平面是接地的無限平面導(dǎo)體,則對于z+半空間、滿足狄利克雷邊界條件的格林函數(shù)為


      {\displaystyle{\begin{matrix}G_{D}(\mathbf{r},\mathbf{r}')={\cfrac{1}{\sqrt{(x-x')^{2}+(y-y')^{2}+(z-z')^{2}}}}\\\qquad\qquad\qquad-\{\cfrac{1}{\sqrt{(x-x')^{2}+(y-y')^{2}+(z+z')^{2}}}}\\\end{matrix}}}\begin{matrix}G_D(\mathbf{r},\mathbf{r}')=\cfrac{1}{\sqrt{(x-x')^2+(y-y')^2+(z-z')^2}}\\


      \qquad\qquad\qquad-\\cfrac{1}{\sqrt{(x-x')^2+(y-y')^2+(z+z')^2}}\\


      \end{matrix};


      其中,{\displaystyle(x,y,z)}(x,y,z)、{\displaystyle(x',y',z')}(x',y',z')分別是檢驗位置{\displaystyle\mathbf{r}}\mathbf{r}、源位置{\displaystyle\mathbf{r}'}\mathbf{r}'的直角坐標(biāo)。


      由于接地導(dǎo)體的電勢為零,方程(2)的面積分項目等于零,方程(2)變?yōu)?/span>


      {\displaystyle\phi(\mathbf{r})={\frac{1}{4\pi\epsilon _{0}}}\int _{\mathbb{V}'}\rho(\mathbf{r}')G_{D}(\mathbf{r},\mathbf{r}')\\mathrmhyxjpf1ac^{3}r'}\phi(\mathbf{r})=\frac{1}{4\pi\epsilon_0}\int_{\mathbb{V}'}\rho(\mathbf{r}')G_D(\mathbf{r},\mathbf{r}')\\mathrmhyxjpf1ac^3 r'


      。


      假設(shè)在位置{\displaystyle(0,0,a)}(0,0,a)有點電荷{\displaystyle q}q,則在z+半空間任意位置的電勢為


      {\displaystyle{\begin{aligned}\phi(\mathbf{r})&={\frac{1}{4\pi\epsilon _{0}}}\int _{\mathbb{V}'}\rho(\mathbf{r}')\left({\frac{1}{\sqrt{x^{2}+y^{2}+(z-a)^{2}}}}-{\frac{1}{\sqrt{x^{2}+y^{2}+(z+a)^{2}}}}\right)\\mathrmhyxjpf1ac^{3}r'\\&={\frac{1}{4\pi\epsilon _{0}}}\left({\frac{q}{\sqrt{x^{2}+y^{2}+(z-a)^{2}}}}-{\frac{q}{\sqrt{x^{2}+y^{2}+(z+a)^{2}}}}\right)\\\end{aligned}}}\begin{align}


      \phi(\mathbf{r})&=\frac{1}{4\pi\epsilon_0}\int_{\mathbb{V}'}\rho(\mathbf{r}')\left(\frac{1}{\sqrt{x^2+y^2+(z-a)^2}}-\frac{1}{\sqrt{x^2+y^2+(z+a)^2}}\right)\\mathrmhyxjpf1ac^3 r'\\


      &=\frac{1}{4\pi\epsilon_0}\left(\frac{q}{\sqrt{x^2+y^2+(z-a)^2}}-\frac{q}{\sqrt{x^2+y^2+(z+a)^2}}\right)\\


      \end{align}。


      仔細(xì)檢察這方程,右手邊第一個項目,是在沒有平面導(dǎo)體的狀況時,點電荷{\displaystyle q}q所產(chǎn)生的電勢;右手邊第二個項目,是使用鏡像法時,鏡像電荷{\displaystyle-q}-q所產(chǎn)生的電勢。請參閱鏡像法條目的點電荷與無限平面導(dǎo)體段落。


    導(dǎo)引


      已知函數(shù){\displaystyle 1/|\mathbf{r}-\mathbf{r}'|}1/|\mathbf{r}-\mathbf{r}'|為格林函數(shù){\displaystyle G(\mathbf{r},\mathbf{r}')}G(\mathbf{r},\mathbf{r}'),滿足方程


      {\displaystyle\nabla^{2}G(\mathbf{r},\mathbf{r}')=-4\pi\delta(\mathbf{r}-\mathbf{r}')}\nabla^2 G(\mathbf{r},\mathbf{r}')=-4\pi\delta(\mathbf{r}-\mathbf{r}')。


      在三維無限空間里,{\displaystyle 1/|\mathbf{r}-\mathbf{r}'|}1/|\mathbf{r}-\mathbf{r}'|的傅里葉級數(shù)為[6]


      {\displaystyle{\begin{aligned}{\frac{1}{|\mathbf{r}-\mathbf{r}'|}}&\equiv{\frac{1}{2\pi^{2}}}\int _{-\infty}^{\infty}\int _{-\infty}^{\infty}\int _{-\infty}^{\infty}\mathrmhyxjpf1ac^{3}k{\frac{e^{i\mathbf{k}\cdot(\mathbf{r}-\mathbf{r}')}}{k^{2}}}\\&={\frac{1}{2\pi^{2}}}\int _{-\infty}^{\infty}\int _{-\infty}^{\infty}\mathrmhyxjpf1ack_{x}\\mathrmhyxjpf1ack_{y}e^{ik_{x}(x-x')+ik_{y}(y-y')}\int _{-\infty}^{\infty}\\mathrmhyxjpf1ack_{z}{\frac{e^{ik_{z}(z-z')}}{k_{x}^{2}+k_{y}^{2}+k_{z}^{2}}}\\\end{aligned}}}\begin{align}\frac{1}{|\mathbf{r}-\mathbf{r}'|}


      &\equiv\frac{1}{2\pi^2}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}\mathrmhyxjpf1ac^3 k\frac{e^{i\mathbf{k}\cdot(\mathbf{r}-\mathbf{r}')}}{k^2}\\


      &=\frac{1}{2\pi^2}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}\mathrmhyxjpf1ack_x\\mathrmhyxjpf1ack_y e^{ik_x(x-x')+ik_y(y-y')}\int_{-\infty}^{\infty}\\mathrmhyxjpf1ack_z\frac{e^{ik_z(z-z')}}{k_x^2+k_y^2+k_z^2}\\


      \end{align}。


      現(xiàn)在,必需找到格林函數(shù){\displaystyle G_{D}(\mathbf{r},\mathbf{r}')=G(\mathbf{r},\mathbf{r}')+H(\mathbf{r},\mathbf{r}')}G_D(\mathbf{r},\mathbf{r}')=G(\mathbf{r},\mathbf{r}')+H(\mathbf{r},\mathbf{r}'),滿足狄利克雷邊界條件{\displaystyle G_{D}((x,y,0),\mathbf{r}')=0}G_D((x,y,0),\mathbf{r}')=0,同時,函數(shù){\displaystyle H(\mathbf{r},\mathbf{r}')}H(\mathbf{r},\mathbf{r}')滿足拉普拉斯方程


      {\displaystyle\nabla^{2}H(\mathbf{r},\mathbf{r}')=0}\nabla^2 H(\mathbf{r},\mathbf{r}')=0。


      對于z+半空間,{\displaystyle H(\mathbf{r},\mathbf{r}')}H(\mathbf{r},\mathbf{r}')以傅里葉級數(shù)擴(kuò)張為


      {\displaystyle{\begin{aligned}H(\mathbf{r},\mathbf{r}')&={\frac{1}{2\pi^{2}}}\int _{-\infty}^{\infty}\int _{-\infty}^{\infty}\\mathrmhyxjpf1ack_{x}\\mathrmhyxjpf1ack_{y}e^{ik_{x}(x-x')+ik_{y}(y-y')}\int _{-\infty}^{\infty}\\mathrmhyxjpf1ack_{z}\left[B(\mathbf{k},z')e^{ik_{z}z}+C(\mathbf{k},z')e^{-ik_{z}z}\right]\\\end{aligned}}}\begin{align}H(\mathbf{r},\mathbf{r}')


      &=\frac{1}{2\pi^2}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}\\mathrmhyxjpf1ack_x\\mathrmhyxjpf1ack_y e^{ik_x(x-x')+ik_y(y-y')}\int_{-\infty}^{\infty}\\mathrmhyxjpf1ack_z\left[B(\mathbf{k},z')e^{ik_z z}+C(\mathbf{k},z')e^{-ik_z z}\right]\\


      \end{align}。


      對于x-坐標(biāo)與對于y-坐標(biāo)的傅里葉級數(shù)擴(kuò)張,{\displaystyle H}H函數(shù)與{\displaystyle G}G函數(shù)的形式相同。這是因為對于無限空間案例與無限平面導(dǎo)體案例,兩種案例的x-邊界條件與y-邊界條件都相同,只有z-邊界條件稍有改變。將{\displaystyle H}H函數(shù)的方程代如,{\displaystyle G_{D}(\mathbf{r},\mathbf{r}')}G_D(\mathbf{r},\mathbf{r}')變?yōu)?/span>


      {\displaystyle G_{D}(\mathbf{r},\mathbf{r}')={\frac{1}{2\pi^{2}}}\int _{-\infty}^{\infty}\int _{-\infty}^{\infty}\\mathrmhyxjpf1ack_{x}\\mathrmhyxjpf1ack_{y}e^{ik_{x}(x-x')+ik_{y}(y-y')}\int _{-\infty}^{\infty}\\mathrmhyxjpf1ack_{z}\left[{\frac{e^{ik_{z}(z-z')}}{k_{x}^{2}+k_{y}^{2}+k_{z}^{2}}}+B(\mathbf{k},z')e^{ik_{z}z}+C(\mathbf{k},z')e^{-ik_{z}z}\right]}G_D(\mathbf{r},\mathbf{r}')


      =\frac{1}{2\pi^2}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}\\mathrmhyxjpf1ack_x\\mathrmhyxjpf1ack_y e^{ik_x(x-x')+ik_y(y-y')}\int_{-\infty}^{\infty}\\mathrmhyxjpf1ack_z\left[\frac{e^{ik_z(z-z')}}{k_x^2+k_y^2+k_z^2}+B(\mathbf{k},z')e^{ik_z z}+C(\mathbf{k},z')e^{-ik_z z}\right];


      其中,{\displaystyle B(\mathbf{k},z')}B(\mathbf{k},z')與{\displaystyle C(\mathbf{k},z')}C(\mathbf{k},z')都是系數(shù)函數(shù)。


      由于{\displaystyle G_{D}((x,y,0),\mathbf{r}')=0}G_D((x,y,0),\mathbf{r}')=0,對于任意{\displaystyle\mathbf{k}}\mathbf{k}與{\displaystyle z'}z',{\displaystyle B(\mathbf{k},z')}B(\mathbf{k},z')與{\displaystyle C(\mathbf{k},z')}C(\mathbf{k},z')之間的關(guān)系為


      {\displaystyle{\frac{e^{-ik_{z}z'}}{k_{x}^{2}+k_{y}^{2}+k_{z}^{2}}}+B(\mathbf{k},z')+C\mathbf{k},z')=0}\frac{e^{-ik_z z'}}{k_x^2+k_y^2+k_z^2}+B(\mathbf{k},z')+C\mathbf{k},z')=0、


      {\displaystyle B(\mathbf{k},z')={\frac{B_{0}e^{-ik_{z}z'}}{k_{x}^{2}+k_{y}^{2}+k_{z}^{2}}}}B(\mathbf{k},z')=\frac{B_0 e^{-ik_z z'}}{k_x^2+k_y^2+k_z^2}、


      {\displaystyle C(\mathbf{k},z')={\frac{C_{0}e^{-ik_{z}z'}}{k_{x}^{2}+k_{y}^{2}+k_{z}^{2}}}}C(\mathbf{k},z')=\frac{C_0 e^{-ik_z z'}}{k_x^2+k_y^2+k_z^2};


      其中,{\displaystyle B_{0}}B_0與{\displaystyle C_{0}}C_{0}都是系數(shù)常數(shù),而且,{\displaystyle B_{0}+C_{0}=-1}B_0+C_0=-1


      將這些公式代入{\displaystyle G_{D}}G_D,可以得到


      {\displaystyle G_{D}(\mathbf{r},\mathbf{r}')={\frac{1}{2\pi^{2}}}\int _{-\infty}^{\infty}\int _{-\infty}^{\infty}\\mathrmhyxjpf1ack_{x}\\mathrmhyxjpf1ack_{y}e^{ik_{x}(x-x')+ik_{y}(y-y')}\int _{-\infty}^{\infty}\\mathrmhyxjpf1ack_{z}\left\{{\frac{(1+B_{0})}{k^{2}}}\left[e^{ik_{z}(z-z')}-e^{ik_{z}(z+z')}\right]\right\}}G_D(\mathbf{r},\mathbf{r}')


      =\frac{1}{2\pi^2}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}\\mathrmhyxjpf1ack_x\\mathrmhyxjpf1ack_y e^{ik_x(x-x')+ik_y(y-y')}\int_{-\infty}^{\infty}\\mathrmhyxjpf1ack_z\left\{\frac{(1+B_0)}{k^2}\left[e^{ik_z(z-z')}-e^{ik_z(z+z')}\right]\right\}。


      為了滿足方程{\displaystyle\nabla^{2}G_{D}(\mathbf{r},\mathbf{r}')=-4\pi\delta(\mathbf{r}-\mathbf{r}')}\nabla^2 G_D(\mathbf{r},\mathbf{r}')=-4\pi\delta(\mathbf{r}-\mathbf{r}'),必需設(shè)定{\displaystyle B_{0}=0}B_0=0。所以,


      {\displaystyle{\begin{aligned}G_{D}(\mathbf{r},\mathbf{r}')&={\frac{1}{2\pi^{2}}}\int _{-\infty}^{\infty}\int _{-\infty}^{\infty}\\mathrmhyxjpf1ack_{x}\\mathrmhyxjpf1ack_{y}e^{ik_{x}(x-x')+ik_{y}(y-y')}\int _{-\infty}^{\infty}\\mathrmhyxjpf1ack_{z}\left\{{\frac{1}{k^{2}}}\left[e^{ik_{z}(z-z')}-e^{ik_{z}(z+z')}\right]\right\}\\&={\frac{1}{|\mathbf{r}-\mathbf{r}'|}}-{\frac{1}{|\mathbf{r}-\mathbf{r}''|}}\\&={\cfrac{1}{\sqrt{(x-x')^{2}+(y-y')^{2}+(z-z')^{2}}}}-\{\cfrac{1}{\sqrt{(x-x')^{2}+(y-y')^{2}+(z+z')^{2}}}}\\\end{aligned}}}\begin{align}G_D(\mathbf{r},\mathbf{r}')&=\frac{1}{2\pi^2}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}\\mathrmhyxjpf1ack_x\\mathrmhyxjpf1ack_y e^{ik_x(x-x')+ik_y(y-y')}\int_{-\infty}^{\infty}\\mathrmhyxjpf1ack_z\left\{\frac{1}{k^2}\left[e^{ik_z(z-z')}-e^{ik_z(z+z')}\right]\right\}\\


      &=\frac{1}{|\mathbf{r}-\mathbf{r}'|}-\frac{1}{|\mathbf{r}-\mathbf{r}''|}\\


      &=\cfrac{1}{\sqrt{(x-x')^2+(y-y')^2+(z-z')^2}}-\\cfrac{1}{\sqrt{(x-x')^2+(y-y')^2+(z+z')^2}}\\


      \end{align};


      其中,{\displaystyle\mathbf{r}''=(x',y',-z')}\mathbf{r}''=(x',y',-z')是鏡像電荷的位置。


    兩個半平面導(dǎo)體案例


      假設(shè)在xy-平面的無限平面導(dǎo)體被一條位于{\displaystyle y=0}y=0的絕緣線條分為兩半,兩個處于y+、y--半平面的導(dǎo)體的電勢分別設(shè)定為{\displaystyle+V}+V與{\displaystyle-V}-V,則由于{\displaystyle\rho(\mathbf{r}')=0}\rho(\mathbf{r}')=0,方程(2)變?yōu)?/span>


      {\displaystyle\phi(\mathbf{r})=-\{\frac{1}{4\pi}}\oint _{\mathbb{S}'}\phi(\mathbf{r}')\{\partial G_{D}(\mathbf{r},\mathbf{r}')\over\partial n'}\mathrmhyxjpf1ac^{2}r'}\phi(\mathbf{r})=-\\frac{1}{4\pi}\oint_{\mathbb{S}'}\phi(\mathbf{r}')\{\partial G_D(\mathbf{r},\mathbf{r}')\over\partial n'}\mathrmhyxjpf1ac^2 r'。(3)


      注意到{\displaystyle\mathbb{V}'}\mathbb{V}'是z+-半空間,xy-平面是其邊界閉曲面的一部分,格林函數(shù)在xy-平面的法向?qū)?shù)的方向是朝著負(fù)z方向:


      {\displaystyle{\begin{aligned}{\partial G_{D}\over\partial n'}&=-\{\partial G_{D}\over\partial z'}\\&=-\{\cfrac{z-z'}{[(x-x')^{2}+(y-y')^{2}+(z-z')^{2}]^{3/2}}}\-\{\cfrac{z+z'}{[(x-x')^{2}+(y-y')^{2}+(z+z')^{2}]^{3/2}}}\\&=-\{\cfrac{2z}{[(x-x')^{2}+(y-y')^{2}+z^{2}]^{3/2}}}\\\end{aligned}}}\begin{align}{\partial G_D\over\partial n'}&=-\{\partial G_D\over\partial z'}\\


      &=-\\cfrac{z-z'}{[(x-x')^2+(y-y')^2+(z-z')^2]^{3/2}}\-\\cfrac{z+z'}{[(x-x')^2+(y-y')^2+(z+z')^2]^{3/2}}\\


      &=-\\cfrac{2z}{[(x-x')^2+(y-y')^2+z^2]^{3/2}}\\


      \end{align}。


      {\displaystyle\mathbb{V}'}\mathbb{V}'的邊界閉曲面在無窮遠(yuǎn)位置的電勢為0,所以,只需要計算xy-平面給出的貢獻(xiàn),就可以得到在{\displaystyle\mathbb{V}'}\mathbb{V}'內(nèi)部任意位置的電勢。將上述方程代入方程(3):[4]


      {\displaystyle{\begin{aligned}\phi(\mathbf{r})&={\frac{2z}{4\pi}}\left\{\int _{0+}^{\infty}\int _{-\infty}^{\infty}{\cfrac{V\mathrmhyxjpf1acx'\mathrmhyxjpf1acy'}{[(x-x')^{2}+(y-y')^{2}+z^{2}]^{3/2}}}+\int _{-\infty}^{0-}\int _{-\infty}^{\infty}{\cfrac{-V\mathrmhyxjpf1acx'\mathrmhyxjpf1acy'}{[(x-x')^{2}+(y-y')^{2}+z^{2}]^{3/2}}}\right\}\\&=\{\frac{zV}{\pi}}\left\{\int _{0+}^{\infty}{\frac{\mathrmhyxjpf1acy'}{(y-y')^{2}+z^{2}}}-\int _{-\infty}^{0-}{\frac{\mathrmhyxjpf1acy'}{(y-y')^{2}+z^{2}}}\right\}\\&={\frac{2V}{\pi}}\\arctan{\left({\frac{y}{z}}\right)}\\\end{aligned}}}\begin{align}\phi(\mathbf{r})&=\frac{2z}{4\pi}\left\{\int_{0+}^{\infty}\int_{-\infty}^{\infty}\cfrac{V\mathrmhyxjpf1acx'\mathrmhyxjpf1acy'


      }{[(x-x')^2+(y-y')^2+z^2]^{3/2}}+\int_{-\infty}^{0-}\int_{-\infty}^{\infty}\cfrac{-V\mathrmhyxjpf1acx'\mathrmhyxjpf1acy'


      }{[(x-x')^2+(y-y')^2+z^2]^{3/2}}\right\}\\


      &=\\frac{zV}{\pi}\left\{\int_{0+}^{\infty}\frac{\mathrmhyxjpf1acy'}{(y-y')^2+z^2}


      -\int_{-\infty}^{0-}\frac{\mathrmhyxjpf1acy'}{(y-y')^2+z^2}\right\}\\


      &=\frac{2V}{\pi}\\arctan{\left(\frac{y}{z}\right)}\\


      \end{align}。


    推廣至電動力學(xué)


      假設(shè)磁場含時間(每當(dāng)電場含時間,則此假設(shè)成立。逆過來亦成立),則不能簡單地以標(biāo)勢{\displaystyle\phi}\phi描述電場。因為根據(jù)法拉第電磁感應(yīng)定律,{\displaystyle\mathbf{\nabla}\times\mathbf{E}=-\{\frac{\partial\mathbf{B}}{\partial t}}\neq 0}\mathbf{\nabla}\times\mathbf{E}=-\\frac{\partial\mathbf{B}}{\partial t}\neq 0,電場不再具有保守性,{\displaystyle\int\mathbf{E}\cdot\mathrmhyxjpf1ac{\boldsymbol{\ell}}}\int\mathbf{E}\cdot\mathrmhyxjpf1ac\boldsymbol{\ell}跟路徑有關(guān)。


      替代地,在定義標(biāo)勢時,必須引入磁矢勢{\displaystyle\mathbf{A}}\mathbf{A},定義為


      {\displaystyle\mathbf{B}\{\stackrel{def}{=}}\\mathbf{\nabla}\times\mathbf{A}}\mathbf{B}\\stackrel{def}{=}\\mathbf{\nabla}\times\mathbf{A};


      其中,{\displaystyle\mathbf{B}}\mathbf{B}是磁場。


      根據(jù)亥姆霍茲定理[7](Helmholtz theorem),假設(shè)一個向量函數(shù){\displaystyle\mathbf{F}(\mathbf{r})}\mathbf{F}(\mathbf{r})滿足以下兩條件:


      {\displaystyle\nabla\cdot\mathbf{F}(\mathbf{r})=D(\mathbf{r})}\nabla\cdot\mathbf{F}(\mathbf{r})=D(\mathbf{r})、


      {\displaystyle\nabla\times\mathbf{F}(\mathbf{r})=\mathbf{C}(\mathbf{r})}\nabla\times\mathbf{F}(\mathbf{r})=\mathbf{C}(\mathbf{r});


      其中,{\displaystyle D(\mathbf{r})}D(\mathbf{r})是個標(biāo)量函數(shù),{\displaystyle\mathbf{C}(\mathbf{r})}\mathbf{C}(\mathbf{r})是個向量函數(shù)。


      再假設(shè){\displaystyle D(\mathbf{r})}D(\mathbf{r})和{\displaystyle\mathbf{C}(\mathbf{r})}\mathbf{C}(\mathbf{r}),在無窮遠(yuǎn)處都足夠快速地趨向0,則{\displaystyle\mathbf{F}(\mathbf{r})}\mathbf{F}(\mathbf{r})可以用方程表達(dá)為


      {\displaystyle\mathbf{F}(\mathbf{r})=-\nabla\left({\frac{1}{4\pi}}\int _{\mathbb{V}'}{\frac{D(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}}\mathrmhyxjpf1ac^{3}r'\right)+\nabla\times\left({\frac{1}{4\pi}}\int _{\mathbb{V}'}{\frac{\mathbf{C}(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}}\mathrmhyxjpf1ac^{3}r'\right)}\mathbf{F}(\mathbf{r})=-\nabla\left(\frac{1}{4\pi}\int_{\mathbb{V}'}\frac{D(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}\mathrmhyxjpf1ac^3 r'\right)+\nabla\times\left(\frac{1}{4\pi}\int_{\mathbb{V}'}\frac{\mathbf{C}(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}\mathrmhyxjpf1ac^3 r'\right);


      在這里,{\displaystyle\nabla}\nabla只作用于{\displaystyle\mathbf{r}}\mathbf{r},體積分的體積為{\displaystyle\mathbb{V}'}\mathbb{V}'。


      采用庫侖規(guī)范(Coulomb gauge),則磁矢勢{\displaystyle\mathbf{A}}\mathbf{A}遵守


      {\displaystyle\mathbf{\nabla}\cdot\mathbf{A}=0}\mathbf{\nabla}\cdot\mathbf{A}=0。


      所以,


      {\displaystyle\mathbf{A}(\mathbf{r})=\nabla\times\left({\frac{1}{4\pi}}\int _{\mathbb{V}'}{\frac{\mathbf{B}(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}}\mathrmhyxjpf1ac^{3}r'\right)={\frac{1}{4\pi}}\int _{\mathbb{V}'}\mathbf{B}(\mathbf{r}')\times{\frac{\mathbf{r}-\mathbf{r}'}{|\mathbf{r}-\mathbf{r}'|^{3}}}\mathrmhyxjpf1ac^{3}r'}\mathbf{A}(\mathbf{r})=\nabla\times\left(\frac{1}{4\pi}\int_{\mathbb{V}'}\frac{\mathbf{B}(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}\mathrmhyxjpf1ac^3 r'\right)=\frac{1}{4\pi}\int_{\mathbb{V}'}\mathbf{B}(\mathbf{r}')\times\frac{\mathbf{r}-\mathbf{r}'}{|\mathbf{r}-\mathbf{r}'|^3}\mathrmhyxjpf1ac^3 r'。


      注意到,以上這些推導(dǎo),并沒有涉及時間參數(shù)。加入時間參數(shù){\displaystyle t}t,結(jié)果也成立。所以,永遠(yuǎn)可以找到磁矢勢{\displaystyle\mathbf{A}}\mathbf{A}:


      {\displaystyle\mathbf{A}(\mathbf{r},\,t)={\frac{1}{4\pi}}\int _{\mathbb{V}'}\mathbf{B}(\mathbf{r}',\,t)\times{\frac{\mathbf{r}-\mathbf{r}'}{|\mathbf{r}-\mathbf{r}'|^{3}}}\mathrmhyxjpf1ac^{3}r'}\mathbf{A}(\mathbf{r},\,t)=\frac{1}{4\pi}\int_{\mathbb{V}'}\mathbf{B}(\mathbf{r}',\,t)\times\frac{\mathbf{r}-\mathbf{r}'}{|\mathbf{r}-\mathbf{r}'|^3}\mathrmhyxjpf1ac^3 r'。


      根據(jù)法拉第電磁感應(yīng)定律,向量場{\displaystyle\mathbf{G}=\mathbf{E}+\partial\mathbf{A}/\partial t}\mathbf{G}=\mathbf{E}+\partial\mathbf{A}/\partial t是一個保守場:


      {\displaystyle\nabla\times\mathbf{G}=\nabla\times\mathbf{E}+\nabla\times\partial\mathbf{A}/\partial t=0}\nabla\times\mathbf{G}=\nabla\times\mathbf{E}+\nabla\times\partial\mathbf{A}/\partial t=0。


      所以,必定可以找到標(biāo)勢{\displaystyle\phi}\phi,滿足{\displaystyle\mathbf{G}=-\nabla\phi}\mathbf{G}=-\nabla\phi。因此,下述方程成立:


      {\displaystyle\mathbf{E}=-\mathbf{\nabla}\phi-{\frac{\partial\mathbf{A}}{\partial t}}}\mathbf{E}=-\mathbf{\nabla}\phi-\frac{\partial\mathbf{A}}{\partial t}。


      靜電勢只是這含時定義的一個特別案例,在這案例里,磁矢勢{\displaystyle\mathbf{A}}\mathbf{A}不含時間。從另一方面來說,對于含時向量場,電場的路徑積分與靜電學(xué)的結(jié)果大不相同:


      {\displaystyle\int _{a}^\mathbf{E}\cdot\mathrmhyxjpf1ac{\boldsymbol{\ell}}\neq\phi(b)-\phi(a)}\int_a^b\mathbf{E}\cdot\mathrmhyxjpf1ac\boldsymbol{\ell}\neq\phi(b)-\phi(a)。



    關(guān)于電勢,小編為大家就分享這些。歡迎聯(lián)系我們合運電氣有限公司,以獲取更多相關(guān)知識。

    上一篇:電感

    下一篇:混合信號集成電路

    相關(guān)新聞

    首頁 產(chǎn)品 手機(jī) 頂部
    在線客服
    聯(lián)系方式

    熱線電話

    15588886921

    400熱線

    400-0886921

    上班時間

    周一到周五

    郵箱地址

    [email protected]

    二維碼

      <object id="bg47o"></object>

      
      

    1. 欧美大香蕉免费看 | 大香蕉人人 | 国产乱人伦久久免费 | 成人黄色一级片 | 丁香五月偷拍 | 五月丁香亚洲无码 | 伊人中文字幕 | 大屌肏屄视频在线播放。 | 国产美女被干网站 | 色之综合天天综合色天天棕色 | av婷婷免费在线观看 | 后入极品美女视频在线 | 天天综合网永久入口 | 国产三级视频在线 | 五月丁香综合在线 | 亚洲AV成人无码久久精品麻豆 | 亚洲狼友 | 北条麻妃 视频 | 精品视频天天在线免费 | 亚洲欧美高清国产精品 | 国产免费AV片在线无码免费看 | 国产粉嫩小泬18p | 日本成人片在线看 | 在线视频二区 | 人妻日韩精品中文字幕 | w黄视频欧美精品韩日 | 黄色片网站在线观看 | 免费无遮挡 视频网乱码 | 精品无码少妇一区二区三区 | 日韩草比视频 | 黄笔毛片大全免费观看 | 一区二区三区在线观看 | 91九色啦 | 大香av| 黄视频免费| 啪啪啪网站大全 | 东京热一讴二区三区 | 日本久久人体视频 | 黄片免费观看 | 日日夜夜狠狠 | 欧美操逼毛片 | 天天摸天天干 | 日日干日韩 | 中曰韩一级片 | 天天干在线影院 | 黄色视频在线免费看电影 | 操你吧欧美性爱 | 亚洲中文字幕无码爆乳av | 亚洲欧洲在线视频 | 国产毛片AV一区二区三区牛牛影视 | 靠逼视频网站在线观看 | 手机在线看a | 爆操视频在线观看免费 | 人人操人人爱超碰 | 日本成人在线不卡视频 | 熊猫成人视频 | 国产高清成人a片 | 求一个做爱视频网站免费在线观看 | 成人精品在线观看 | 4438激情网 | 国产一级电影 | 韩国三级av | 巨乳全裸家证父 | 国产精品内射无码 | 国产1234在线电影 | 中文字幕五码 | 久久亚洲精品成人777 | 成人A片无码水蜜桃免费网站软件 | 国产性爱第一页 | 一级a毛片免费看 | 天天干天天日天天干天天日 | 欧美性爱第四页 | 91人妻精 | 亚洲欧洲免费看 | 翔田千里黑人一区在线观看 | 肏屄视频在线免费观看 | 在线黄色片视频 | 中文人妻av | 男人天堂2025手机版 | 大吊操逼视频 | 国产无遮挡裸体色视频免费观看 | 欧美成人精品首页 | 91久久婷婷亚洲精品成人 | 亚洲黄色免费电影 | 久草青娱乐 | 成人久久 挤奶 | 亚洲福利网站 | 日韩在线观看一区 | 国产亚洲色婷婷久久99精品 | 激情婷婷色色 | 俺也去色网络 | 开心色播五月天 | 欧美日韩丁香五月天色网 | 中文字幕一区免费 | 成人一级电影 | 午夜色综合 | 久久九| 黄色片视频日韩 | 欧美一性一乱一交一视频 | 五月天AV电影在线 | 大伊香蕉在线播放 | 国产男女啪啪网站 | 夜夜爽AV | 亚洲精品成人片在线播放波多野吉 | 免费看一区二区三区A片 | www.婷婷综合 | 欧美日本黄色一级视频 | a天堂中文在线 | 青娱乐在线视频网站 | 日韩福利一区二区 | 美女视频黄a视频全免费不卡 | 俺去也俺就去www色情网最新更新内容 | 偷拍日本大香蕉 | 操逼视频。 | 在线亚洲黄色视频 | 日韩欧美第一页 | 日产精品久久久久 | 亚洲成人无码专区在线 | 91久久精品国自产合 | 67194亚洲 | 一级黄色录像带 | 欧美精品一卡二卡 | 色婷婷国产精品视频一区 | 欧美美女后进式插逼视频 | 欧美一级颜情免费视频 | A片黄色电影网站 | 国产精品豆花视频www | 国产久久久久 | 欧美视频手机在线观看 | 国产精品77 | 你操我操综合 | 国产无码激情视频 | 操逼视频播放 | www.xxx操逼网 | 成人网站在线免费 | 国产一级a一级a看免费视频 | 天天添夜夜干 | 撸一撸综合 | 操逼的网 | 成人视频久久 | 大香蕉久久久 | 99精品视频99 | 国产婬| 精品黄色 | 三级片网站在线观看视频 | 日韩无码一二三四区 | 日韩一级一级 | 欧美操逼视频免费观看 | 99婷婷五月天 | 18禁精品 | 麻豆精品秘 国产 | 日本爱爱一区二区视频 | 天天爽夜夜爽夜夜爽 | 操操操电影| 伊人视频网站 | 天天撸一撸在线免费观看 | 丁香五月激情中文 | 熟女国产精品 | 五月婷婷综合网 | 亚洲福利视频网站 | 夜色视频在线播放 | 亚洲V免费看 | 五月丁香月 | 黄色电影在线视频 | 欧美熟女一区 | 日本高清色情视频 | 国产精品久久电影网 | 经典三级先锋影音 | www.99热这里只有精品 www国产夜插内射视频网站 | 九九九九在线视频 | 久久精品观看 | 国产精彩视频免费观看 | 国产黄色免费观看 | 亚洲一级日比毛片 | 高清一区在线观看视频 | 国产精品成人无码a无码 | 天天日天天操天天插 | 黄色电影A片电影 | 玖草福利 | 久久狼人综合香蕉 | 国产1234在线观看 | 三级片大香蕉 | 亚洲人妻蜜桃 | 爱爱视频天天 | 免费观看亚洲视频 | 国产成人电影久久 | 熟女三区 | 91豆花视频 | 一级操逼 | 果冻传媒性爱操逼电影三级片 | 老色批影院 | www.青青操 | 手机免费A V | 日韩无码黄色电影网站 | 依人大香蕉视频网站 | 久久熟女 | 波多野结衣在线网站 | 台湾成人在线观看 | 特黄特色大片勉费看 | 久久精品视频在线观看 | 91久久精品人人槡人妻人 | 国产美女操逼网站 | 免费a√在线 | 香蕉在线观看 | 激情射av | 夜夜干天天撸时时操 | xxxxx无码 | 亚洲视频在线播放免费 | 影音先锋在线成人视频 | 激情A学生妹 | 久艹久艹 | 无码肉厚美臂小早川怜子 | 麻豆精品无码视频 | 国产无码久久久久 | 日本大乳高潮视频在线观看 | 亚洲AV无码久久久久网站飞鱼 | 日韩av一卡电影在线观看 | 久久久久久久久久一级 | 伊人大香蕉在线观看视频 | 在线免费观看小黄片 | 小视频一区 | 波多野吉衣中文字幕在线 | 影音先锋AV成人资源网 | 欧美性爱视频精品 | 91久色| 亚洲AV成人漫画无码精品网站 | 大香蕉操骚逼 | 欧美一级视频在线免费观看 | 毛毛毛毛毛毛毛片123 | 在线毛片网址 | 国产高清无码视频在线播放 | 美女草逼视频 | 黄色成人视频免费 | 久久综合色气网站 | 午夜性爱网站 | 欧美成人性爱在线视频 | 不卡日本视频 | 欧美日产久久 | 国产一級A片免费看 | 国产精品久久影院 | 人妻日逼网| 丁香婷婷在线 | 天堂极品mv | 欧美性爱视频福利 | 青娱乐国产一区二区三区 | 亚洲AV蜜桃永久无码精品XXX | 含羞草自拍原创一区二区三区视频 | 天天久天操天天淫 | 欧美熟妇乱伦 | 91无码国产成人精品 | www第一页 | 久久精品AV无码夜色 | 欧美三级高清在线观看 | 成人才看的在线视频 | 日本一级免费视频 | 樱桃综合网 | 亚洲理论在线观看 | 99这里只有 | 青青激情| 欧洲精品无码一区二区在线 | 91久久久久久久久 | 99热在线只有精品 | 男女动态拍拍片 | 欧美日韩国产三级 | 水蜜桃精品视频 | 性爱无码网站 | 国产精品18久久久 | 亚洲天堂在线看 | 在线观看国产色情 | 成人精品91 | 欧美成人性爱网址 | 草久网 | 青青草人人艹 | 暗呦网一区二区三区 | 欧美在线va | 久久网站又改成什么 | 97人人色 | 骚无码 | 久久婷婷婬片A片AAA | 成人黄页网 | 无码在线中文字幕 | 伊人婷婷狠狠 | 色婷婷在线无码精品 | 欧美俺来也 | 美女扒开尿口让 | 综合在线第一页 | www.免费看黄色 | 黄色电影视频看看 | 国产特黄 | 成人视频网站在线看 | 一区日韩 | 一区二区电影网 | 国内自拍夫妻视频 | 免费草逼网站 | 精品在线免费视频 | 亚洲国产一级黄片 | 中文字幕精品无码网站人口 | 国产传媒精品 | 日韩欧美18禁 | 91成人精品一区在线播放 | 加勒比很很操免费视频 | 五月婷婷欧美激情 | 操屄视频免费观看 | 国产精品一区二区三区四区 | 天天干天天干天天天干 | 超碰在线人人操 | 91精品国产91久久久久 | 91麻豆精品国产 | 亚洲乱码国产乱码精品天美传媒 | 亚洲一区二区三区在线 | 国产激情视频久久久久久 | 亚洲黄色自拍 | 精品人妻一区二区三区含羞草 | 无码 操| 特黄特色特刺激免费播放 | 在线免费看黄 | 成人伊人影视在线 | 亚洲无码高清视频在线 | 欧美性爱亚洲性爱 | 欧美肏屄视频网站 | 日本成人黄色视频 | 日韩无码高清视频看看看看 | 国产婷婷在线视频 | www视频网站在线操 | 永久黄色WWW.99 | 亚洲免费日本精品在线观看 | 老熟女 码A片 | 影音先锋自拍 | 国产人妻成人免费无码 | 一级a黄色电影片 | 欧美日韩国产成人 | 在线a黄网站 | 黄频免费看 | 久久久午夜视频 | 国产精品视频免费在线观看 | 天天天干夜夜夜 | 天天色天天日天天射 | 成人精品视频 | 天堂综合精品无码 | 亚洲视频免费观看视频 | 人高清无码在线播放 | 国产精品乱码69一区二区三区 | 国产精品永久成人免费 | 天天干天天插天天 | 天堂成人在线视频 | 国产夫妻视频 | 一级A黄色片 | 青青草无码在线观看 | 大香蕉综合在线 | 蜜臀久久99精品久久宅男 | 日韩操逼毛片 | 中国大香蕉黄色在线视频 | 天天射日日舔 | 色色色五月天婷婷 | 国产三级三级三级看三级囯产亚洲 | 精品午夜人妻 | 日本一区二区三区在线观看视频 | 看片一区 | 91日韩欧美在线观看 | 国产成人丁香五月 | 插逼视频免费看 | 91精品国产福利 | 亚洲天堂视频网 | 欧美一区二区三区的 | 波多野结衣一区二区 | 台湾无码一区二区 | 中文字幕在线一区观看 | 处一女一级a一片 | 天天爽夜夜爽人人爽 | 北条麻妃一区二区三区 | 日韩无码流出 | 97国产成人无码精品久久久 | 在线欧美视频 | 7799日日夜精品日日夜精品 | 99Re视频官网 | 奶水饱胀一区二区三区 | 亚洲ww| а√天堂资源在线 | 日韩色逼| 国产裸体视频网站女69 | 99日韩无码| 日韩精品一区二区三区四区苍老师 | 竹菊视频一区二区三区 | 欧美黑人大鸡巴操逼视频 | 亚洲性爱网络 | 麻豆操逼毛片视频 | 日批视频免费播放 | se01亚洲 | 亚洲国产精品成人综合久 | 樱桃视频一区二区 | 欧美一级a视频免费放 | 我想看中国台湾特级黄色录像1级特黄特黄的 | 国产人人人操 | 婷婷五月天在线观看 | 51ⅴ精品国产91久久久久久 | 粉嫩在线视频 | 九九九精品在线 | 欧美精品一区二区三区成人片在线 | 久久一区二区三区在线 | 欧美高清在线 | 日韩国产在线一区 | 欧美激情成人 | 欧美日韩在线播放 | 久久久精品影院av | 国产网红女主播操逼 | 草肉肉XXXXHD劉亦菲 | 亚洲中文字幕在线免费观看视频 | 免费看成人747474九号视频在线观看 | 欧美日高清视频免费在线播放 | 亚洲欧美精品suv | 国人免费无码区久久久免费 | 亚洲1区2区3区4区5区 | 大香蕉网欧美 | 天堂在线1 | 久久性| 久热在线精品视频 | 九一精品网站 | 99亚州| 日韩理论在线 | 人妻无码久久精品人妻成人 | 亚洲高清视频在线播放 | 黄片视频链接 | 久久视频网站 | 黄色小电影在线看 | 国产精品久久久久久久久久久久午夜片 | 一级黄色电影免费看 | 波多野结衣视频在线 | 黄色电影视频在线精品 | 国产乱婬A∨片 | 成人在线三级片 | 超碰在线大香蕉 172.86.93.25 | 色婷婷五月天在线 | 国产丝袜人妖CD露出 | 色婷婷激情五月天 | 亚洲综合色网 | 青青青亚州视频在线 | 淫秽视频在线观看 | 亚洲二视频 | 又黄又爽视频 | 日本黄色片免费看 | 狠狠操狠狠爽 | 亚洲一级黄片免费观看视频洲免费观看 | 亚洲第一页综合 | 在线观看成人自拍 | 亚洲免费无码视频 | 精品韩无码台湾成人在线 | 大大鸡吧轻轻操在线视频 | 嫩逼网| www.99热这里只有精品 www国产夜插内射视频网站 | 成年人视频在线播放 | 九一免费观看网站 | 麻豆成人在线 | 欧美日本AA | 国产精品自拍偷拍 | 欧美性爱免费网址 | 人人草人人爽 | 中文字幕一区第一页 | 日韩AV性爱 | 亚洲性爱电影网站 | 爱情岛黄色 | 一级A片一毛片大全 | 色激情五月天 | 操骚逼黄色网址 | 欧美特级AAAAAA特黄高清 | 91av在线观看视频五月天 | 正在播放:欧美系列黑人老哥干大奶少妇 | 欧美性猛交XXXXXX | www.8x8x | 做爱黄色视频网站 | 东京热先锋影音 | 免费一级片电影网站 | 京熱大亂交无碼大亂交在线 | 91福利视频在线观看 | 日韩视频――中文字幕 | 亚洲va高清 | 苍井空在线一区二区 | 老阿姨的丝丝波涛胸涌诱惑 | 四虎8848精品成人免费网站 | 777国产盗摄偷窥精品0000 | 女人操女人一级毛片 | 国内亚洲欧美精品九九九久久免费 | 色综合五月天 | 天堂а√在线中文在线新版 | 99精品无码视频 | 日韩A片在线免费观看 | 第四色大香蕉 | 久久精品国产亚洲7777 | 欧美闷骚影院 | 国家产在线观看一二 | jk无码| 污污污午夜 | 国内自拍在线观看 | AV软件在线免费观看 | 色婷婷粉嫩精品综合在线 | 欧美日韩一区二区黄 | 在线免费无码 | 色婷婷中文在线观看视频 | xxxxx无码 | 操逼手机高清免费看 | 大香蕉美女啪啪啪视频 | 欧美精品人妻视频 | 日韩电影无码一区二区 | 亚洲日韩小说图片视频首页 | 麻豆国产精品无码人妻无码 | 亚洲精品成人在线视频久久 | 大香蕉欧美伊人 | 日韩a毛片 | 婷婷国产夫妻 | 国产熟女视频 豆花视频 | 自拍偷拍成人在线视频 | 免费A∨视频 | 一级片电影网站 | 成人18毛片| 黄色免费国产电影 | 亚韩一区二区三区 | 久久亚洲中文字幕 | 你懂的在线观看视频 | 亚洲理论视频 | 国产露脸91国语对白 | 神尾舞av| 久久夜天天天天 | aa色黄视频 | 日屄视频免费观看 | 婷婷久久福利 | 一级a性色毛片 | 精品久久久久中 | 日韩特级视频 | 99热在线只有精品 | 97人妻精品一区二区三区 | 羞羞色院91蜜桃 | 婷婷狠狠干 | 国际精品一级视频 | 中文字幕在线永久 | 国产精品久久久久永久免费看 | 91丨国产亚洲丨精品白丝 | 影音先锋AV男人资源站 | 久久超碰大香蕉 | 超碰自拍中文字幕 | 99自拍网 | 东京热99 | 香蕉久久a毛片 | 国产精品小电影 | 91精品老司机 | 狠狠干狠狠干 | 国产欧美鸡巴 | 美女激晴一级播放在线观看 | 无码av人| 欧美性成人 | 大香蕉综合一二 | 青娱乐亚洲高清视频 | 亚洲无码家庭乱论小说区 | 亚洲AV无码乱码国产精品牛牛 | 俺去俺来也在线www色 | 一道无码| 欧美精品综合 | 欧美强开小嫩苞 | 久久艹骚逼| 欧洲免费观看成品视频 | 欧美老熟妇乱大交XXXXX动漫 | 天天爽天天摸天天爱 | 艹骚逼| 国内精品三级视屏 | 又大又粗免费视频 | 狼人综合免费视频 | 亚洲黄色电影网站 | 在线无码免费观看 | 六区九区一区在线 | 逼逼爱插| 精品无码一区二区三区爱与 | 囯产精品宾馆在线精品酒店 | 看着撸影音先锋资源 | 免费在线观看亚洲视频 | 三级视频手机在线播放 | 天天日天天弄 | 亚洲福利视频网 | 91北条麻妃二区 | 九色丨蝌蚪丨成人 | 在线播放高清无码 | 日韩成人黄色电影网站 | 少妇人妻无码AV片在线蜜芽 | 青娱乐av免费观看 | 狠狠狠狠狠狠狠狠狠狠狠狠狠狠狠狠 | 午夜成人性爱 | 亚洲娱乐在线 | 视频一区三区 | 三级日本三级网站三级网站在线 | 国产一级乱伦片 | 欧美日韩一区二区三区四区论理片 | 91久久成人视频 | 免费无码婬片AAAA片视频 | 美女扒开尿口让男生桶爽视频 | 免费一级黄色片 | 农村嫩苞一区二区三区… | 国产国内精品灰色夫妻网 | 免费看一区二区三区四区 | 国产伦精品一区二区三区视频痴汉 | 午夜精品少妇 | 狠狠干AV | 加勒比综合 | 色五月 丁香 婷婷 | 亚洲免费观看高清在线观看 | 大香蕉网大香蕉网 | 豆花传剧高清在线看 | 国产一卡二卡三卡免费观看 | 精品无码一区三区三区 | 中文字幕日韩视频 | 黄色小说片五月 | 北条麻妃网址大全 | 国产激情视频免费在线看 | 爱婷婷五月丁香 | 51精品一区二区三区 | 吖v在线视频免费观看免费观看 | 熟女乱伦图片 | 一区二区三区免费看 | 中文字幕一区二区三区四区五区人 | 国产免费aaa高清 | 亚洲一区二区免费视频 | 欧美成人精品一级乱黄 | 国产偷拍自拍在线观看 | 看黄色大片网站操逼 | 操逼图片视频 | 91嫩草欧美久久久九九九 | 天堂中文在线资源的 | 一级黄色片中文字幕 | 国产最新网址 | 韩国精品无码电影 | 韩国成人三级 | 青娱乐精品在线视频 | 三级黄色天天天天 | a√天堂资源中文8 | 大香蕉大香蕉视频网 | 黄色一级大片在线观看 | 久久婷婷丁香五月综合 | 欧美成人精品三级在线网站 | 日韩专区第一页。日韩中文字幕在线亚洲 | 亚洲色综合 | 蜜臀VA | 日韩在线黄| 国产精品又污又黄又爽污污 | 欧美 日韩 国产 中文 | 欧美成人一区二区 | 亚洲AV毛片在线观看 | 亚洲精品少妇 | 日韩精品熟妇 | 开心色情五月综合激情 | 大色吻97综合 | 五月婷丁香 | 色大香蕉色大香蕉色大香蕉色大香蕉色色 | 黄色视频在线观看免费 | 靠逼网站免费观看 | 国产精品高潮呻吟久久AⅤu的 | 曰韩一级A片 | 天天草天天干 | 黄色做爱日本动漫网站 | 久草青青草 | 色婷婷一区二区三区久久午夜成人 | 久久久久久69精品 | 色五月婷婷综合 | 爱操逼网| www.欧美成人 | 国产乱子伦-区二区三区熟睡91 | 天堂久久网 | 天天爽天天日 | 九九九九精品视频 | 国产很很鲁 | 私人女仆扫地偷懒被主人颜色吃现在被喷尿洗脸 | 欧美韩三级片 | 国产精品人妻人伦a 6 2v久软件 特级西西444www无码视频免费看 | 免费观看国产一卡二卡电影 | 天天骚天天骚 | 天堂色色 | 黄片一区 | 国产无码一级视频 | 亚洲天堂91 | 中文字幕无码视频在线 | 干BB视频 | 日韩一级电影在线观看 | 中文字幕一区二区久久人妻网站 | 久久在线精品 | 毛片毛片毛片毛片 | 人妻免费视频 | 亚洲乱伦视频网址 | 91视频久久 | 成人精品免费视频 | 日韩欧美精品 | 久久艹伊人 | 天天爽天天爽成人A片影院 | 国产成人MV | 无码中文字幕在线视频 | AV无码免费一区二区三区不卡 | 操大逼 | www.想要xx | 嗯嗯无码| 大香蕉在线精品视频 | 全部在线A片免费播放 | 成人网站日皮视频 | 一级片日日爱66 | 亚洲乱伦一区二区三区 | 日韩成人片 | 888精品成人片 | 久久精品国产亚洲A | 美女黄18禁 | 特级日本欧美日韩亚洲精品综合免费在线 | 91麻豆成人精品国产免费网站 | 日皮视频在线观看免费 | 日韩精品777视频 | 国内毛片毛片毛片毛片毛片毛片毛片毛片毛片毛片毛片毛片 | 国产日逼逼 | 大香蕉最新视频精品 | 2018天天日天天射 | 日日三级网 | 国产成人综合久久久久久 | 久久精品黄色视频 | 国产麻豆成人免费视频 | 精品少妇在线刺激 | 大屌肏屄视频在线播放。 | 国产一级A片久久久免费看快餐 | 中文字幕乱伦 | 精品三级AV无码 | 婷婷一区二区三区 | 日韩一级电影观看 | 久久久久久国产精品高清 | 成人导航网站 | 亚洲免费观看A∨中文 | 亚洲一区影音先锋 | 日韩欧美在线视频 | 夜夜嗨AⅤ一区二区三区 | 成人视频网 | 1204人妻一区二区三区 | 日批视频免费看 | 9999热这里只有精品 | 九九九九精美视频 | 国产精品久久久豆花视频 | 婷婷五月天丁香 | 美腿丝袜亚洲综合 | av黄色电影一区天堂一区二区三区 | 天天干天天爱av 天天干天天骑AV 天天干天天色免费 | 中国婬乱a—级毛片多女 | 东京热亚洲无码 | 青草色天堂 | 99热在线日韩精品免费 | 奇米一区| 国产九色| aaa天堂网 | 伊人大香蕉www | 欧美午夜成人视频 | 就要操妓女网 | 91视频 久久 | 一级特黄视频 | 日本一级A片在线观看 | 欧美色图9 | 亚洲成人网站免费在线观看 | 操逼高清 | 欧美久久久久久久久久 | 少妇操屄视频 | 精品产在线播放内射 | 亚洲黄色片 | 丁香婷婷色五月激情深爱 | 欧美一区二区三区四区精品 | 逼中逼啪视频 | 大鸡巴在线看 | 免费a黄色电影 | 97自拍视频 | 新人妻中文字幕网站 | 大香蕉国产毛片儿操逼申影 | 多人无码| 久久成人无码电影 | 苍井空肏屄乱伦 | 日韩一区二区免费视频 | 777.av | 一级二级黄色电影 | 黄色毛片一级 | 日本欧美一级色情片免费观看 | 国产传媒天美果冻精品亚洲一区二区三区直播 | 青春草免费在线视频 | 九九九九九精品视频 | 成人午夜在线 | 91精品免费久久 | 久久大鸡八三级片 | 免费一级特黄大学生毛片 | 一级二级三级在线观看 | 欧美大香蕉在线观看免费一区二区三区 | 国产成人 亚洲 日韩 | 高清国产AV | www.撸| 日本成人五月天 | 黄色强奸免费小视频网站 | 一插综合 | 香蕉国产成人毛片 | 豆花视频理论在线播放 | 亚洲无码精品久久 | 国国外苗毛片 | 蜜桃视频成人网站入口 | x88AV~熟女人妻 | 秋霞福利在线 | 91探花精品偷拍在线播放 | 波多野结衣中文字幕一区二区 | 大香蕉日| 99在线精品免费视频 | 一区二区三区无马亚 | a黄色在线视频无 | 日逼大全 | 日本中出中文字幕 | 天天噪天天射天天拍 | 狠狠操五月天 | 自拍偷拍成人视频 | 亚洲综合 p | 欧美videos办公室丝袜长腿 | 五月天无码视频 | www.操操操.com | 开心五月深深爱婷婷和丁香 | 青青草天天干 | 国产特级毛片AAAAAA | 天天曰天天操 | 男插女视频网站 | 无码三级视频在线观看 | 国产黑料视频你懂的 | 99国产在线观看 | 逼特逼视频免费看 | 欧美 日韩 国产在线观看 | 五月天色婷婷综合 | 午夜偷拍精品 | 一级大黄a片 | 色婷婷综合av | 小泽玛利亚av在线 | 性日熟妇 | 高清无码视频在线免费观看 | 色婷婷亚洲色 | 成人大香蕉综合电影 | 欧美一区二区三区四区视频 | 亚洲视频在线免费看 | 一区二区三区四区五区精品无码 | 欧美成人手机免费在线观看 | 久操激情网 | 激情乱伦91网站 | 日日av影院 | 日逼片| 99国产在线视频 | 久热在线观看视频 | 国产一区二区三区四区五区入口 | 欧美在线黄色 | 青青青青青青草草草草草草草视频 | 一级淫乱视频 | 国产波霸巨爆乳无码视频在线 | 久久99精品久久久水蜜桃 | 久久伊人大香蕉 | 亚洲成人免费AV | 国产精品婷婷久久久 | 色色色色五月天 | 精品国产污污污免费入口15 | 亚洲中文字幕不卡 | 影音先锋成人网 | 亚洲精品黄色电影网站 | 成人午夜色情无码 | 十八禁网站在线看 | 成人三级电影在线观看 | 大肠浣肠调教一区二区三区在线 | 免费做爱视频网站 | 日本激情视频网站免费 | 加勒比蜜桃日韩成人 | 欧美日韩免费看 | 国产靠逼视频网站 | 日韩成人激情 | 色播视频在线观看 | 88无码| 黄色视频高清无码 | 天天日天天射天天 | 豆花视频成人版WWW18 | 韩国三级久久 | 久久无码一区 | 亚洲乱伦无码视频 | 国产久久久久久久 | 国产AV黄片 | gg在线精品视频观看 | 天天日狠狠操 | 视频一区二区三区免费 | 成 人 在线偷拍视频 | 天天天天射天天天搞天天要 | 欧美伊人 | 激情图区| 天天玩天天操 | 亚洲日韩中文字幕视频在线 | 99re99 | 99久久精品国产精品有折扣吗 | 亚洲成人社区网站 | 91人妻人人澡人人爽人人精品乱 | 日韩一级在线播放免费观看 | 学生妹无遮挡毛片 | 欧美视频天天干 | 国产精品久久欠久久久久久九秃 | 亚欧精品久久久久久久久久久 | 中文字幕av一区二区三区 | 中文无码熟妇人妻 | 天天干,天天射免费视频 | 北条麻妃无码一区二区三区视频 | 大香蕉久久久久久久久久久 | 男女激情91视频网站 | 亚洲在线视频第一页 | 日批在线观看 | 国产黄色在线看 | 日韩午夜| 成人AV中文字幕 | 亚洲影院之台湾 | 免费黄色视频。 | 操逼精品免费视频 | 成人无码影音先锋 | 色婷婷激情综合网 | 97碰人人操 | 水蜜桃成人A片视频网站 | 日日干夜夜爽 | 婷婷毛片 | 少妇喷水视频 | 影音先锋久久久久久 | 久久婷婷五月综合伊人 | 伊人久热亚洲精品 | 久久中文字幕网 | 五月天婷婷影院 | 大香蕉福利视频 | 天天爽天天爽天天 | 日本一级片免费看 | 四虎影院最新地址 | 第四色五月天 | 欧美交换乱婬A片免费 | 天天做爱视频 | 亚洲色图88 | 亚洲成人欧美 | 大鸡巴操我骚逼 | 激情视频福利一区在线看久 | 日本视频3号级 | 午夜抽插 | 亚洲视频第1页 | 爱搞搞视频 | 亚洲无码1000 | 欧美干在线观看 | 东京热先锋影音 | 成人电影91 | 青青欧美,青青美女视频 | 日韩欧美伊人 | 成人三级无码视频 | 一级A片黄片 | 操出水视频在线观看网站国产 | 欧美色18在线 | 国产精品伦理一区 | 亚洲色大成 | 九九久久网站 | 亚洲精彩视频 | 手机在线观看免费视频人 | 青青操手机在线观看 | 黄色在线免费一级视频 | 91chaopeng在线 | 激情无码青青草 | 9l视频自拍九色9l视频成人 | 以及视频片又粗又猛 | 西西444WWW无码大胆 | 黄色片在线免费 | 99爱精品视频在线观看 | 中文字幕北条麻妃在线 | 暗呦网一区二区三区 | 成全在线观看高清的 | 青青操拍拍视频 | 99青娱乐 | 操逼喷水视频 | 五月婷婷综合视频 | 欧美一区二区三区成人片下载 | 爱AV操 | 欧性猛交XXXX乱大交 | 日韩无码三级片 | 成人性生交大片免费看黄103 | 午夜影院在线 | 亚洲精品乱码久久久久久蜜桃不卡 | 加勒比精品 | 久久久麻豆 | 国产草逼视频 | 亚洲AV成人精品一区二区三区_ | 欧美性猛交XXXX乱大交 | www.狠狠鲁 | xxxxx免费视频 | 免费观看一区二区三区四区五区 | 国产大屁股 | 在线观看的黄色小视频 | 玖玖精品视频在线观看 | 激情色五月天 | 成人做爱网站视频| 风间精品一区二区三区 | 在线内射视频 | 中国操逼毛片 | 91高清无码在线观看 | 免费看黄色一级片 | 青青草草草在线视频资源站 | 天天操中文字幕 | 日本猛少妇色XXXXX猛叫 色秘 乱码一区二区三区唱戏 | 成人三级在线看 | 自拍一区不卡 | 欧美一级免费A片 | 爱爱网站无码 | a级黄色的网站在线观看 | 性爱免费视频毛片 | 黄色片视频日本 | 国产精品成人久久免费 | 日本理论片一道本 |