<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. 精品三级在线观看 | 一二三精品区 | 国产无码艹逼 | 最新一区二区在线 | 日韩一级中文字幕 | 欧美成人手机 | 欧美黄色片视频 | 欧美日韩在线观看视频 | 爽 好紧 别夹 喷水网站 | 日本成人三级视频 | 青青久操视频 | 无码中文字幕视频在线观看 | 天天日天天操天天摸 | 在线a视频免费观看 | 一级av黄片 | 久久精品噜噜噜成人88Aⅴ | 大黑屌视频 | 学生妹一区二区 | 亚洲AV无码成人精品区在线播放 | 日本一级持黄 | 蜜桃在线无码视频 | 久久精品无码视频 | 欧美激情一区二区 | 午夜免费爱爱视频 | 影音先锋乱伦电影 | 欧美日韩黄色 | 8x8x成人免费视频入口 | 黄色免费高清视频 | 香蕉视频在线观看网站资源 | 在线看片日韩 | 四虎影库亚洲无码 | 日日操夜夜叫狠狠插 | 狠狠操夜夜爽 | 免费的看污片丝瓜视频 | 亚洲热情 | 日本视频在线观看中文字幕亚洲 | 欧美内射网 | 特级西西人体444.444人体聚色 | AV天堂资源 | 欧美性手机在线 | 人妖操女人 | 国产黄永久 | 中国航空少妇一级毛片 | 日韩欧美中文在线视频 | 不卡的a12在线观看 | 天堂综合网久久 | 黑人大屌轮奸视频播放免费成人 | 天天色天天干天天日 | 日韩一级a免费在线视频 | 欧美在线无码精品秘 蜜桃 | 日批网站在线看 | 日韩性爱网址 | 国产激情在线 | 成人黄色在线看 | 综合插插插| 欧美日韩中文字幕在线观看 | 男人天堂V| 91视频久久久久久久久久久久 | 激情视频乱伦 | 天天色人人 | 免费看又黄又无码的网站 | 国产一级黄色电影 | 中文精品久久久久久 | 欧美日韩中文 | 无码内射国产 | 91av无码 | 无码人妻一区二区三区精品不付款 | 亲子伦视频一区二区 | 欧美一级A片高清免费播放 | 果冻传媒性爱操逼电影三级片 | 亚洲精品粉嫩小泬18p图片 | 青青青视频在线 | av在线日3 | 成人网站观看免费 | 国产精品人妻在线 | 欧美日韩手机看片免费 | 色欲天香天天免费视频 | 欧美18禁成人 | 五月丁香婷婷综合在线 | 狠狠操狠狠色 | 豆花视频在线观看一区二区 | 青青草免费在线公开视频播放 | 大香蕉久久久久久成人草 | 美女高潮视频在线观看免费视频 | 日本黄色视 | 无码网址| 亚洲欧美国产毛片在线 | 在线观看日韩 | 婷婷在线无码 | 乱伦毛片视频 | 老鸭窝毛片 | 国产欧美久久精品 | 亚洲无码电影大全 | 无码视频高清在线不卡 | 久久久久久国产精品三级玉女聊斋 | aa,xx黄色 | www.欧美成人 | 国产精彩 | 99精品成人免费毛片无码 | 国内免费毛片一区二区 | 国产成人无码区免费视频 | 亚洲AV无码成人精品涩涩麻豆 | sm视频网站| 免费尻屄视频 | 久久福利视频精品 | 亚洲精品国产精品乱码桃花 | 成人毛片女人毛片免费96 | sese亚洲 | 外国黑人的一级视频 | 日本在线观看a | 成人视频在线播放 | 欧美黄色免费第一看 | 色九九九九九九 | 色五月综合| 肏逼网站在线观看 | 黄片精品午夜福利在线免费观看豆花视频 | 日本亲子乱淫一级 | 三区视频免费看 | 岛国免费AV | 国产午夜免费视频 | 91爱爱 | 看黄色操逼视频 | 天天舔天天射 | 麻豆激情四射在线播放 | 国产高清porn | 欧美精品 - 91爱爱 | 亚洲黄色电影在线看 | 豆花成人www、C0m | 91av精选 | 91男女操逼的视频 | 我要毛片毛片毛片毛片毛 | 日韩高清无码网站 | AAA片网站 | 国产网站91 | 自拍做爱视频 | 国产精品日韩 | 精品av国产 | 欧美啪啪视频 | 国产色亭亭视频 | 台湾成人久久网站视频 | 久操视频中文在线 | 欧美福利在线视频 | 超碰99免费 | 97AV电影 | 国产黄色片在线看 | 一区性感在线观看 | 国产乱婬AAAA片视频软件 | 亚洲新人天堂中文 | 成人福利午夜A片公司 | 影音先锋男人网 | 影音先锋久久久久AV综合网成人 | 河南首富 越来越富 | 国产午夜精品视频 | 成年人视频中文字幕在线播放 | 精品久久久麻豆 | 中国a一片一级一片 | 91黑人大屌啪啪 | 日韩破解无码片 | 亚洲免费观看高清 | 高清无码免费视频在线观看 | 日韩在线三级 | 天堂网成人在线 | 欧美自拍视频在线 | 人人爽,人人妻,人人操 | 色五月av综合 | 日韩欧美一级特黄大片 | 午夜亚洲AⅤ无码高潮片苍井空 | 91精品少妇一区二区三区蜜桃臀 | 艹逼视频免费看 | 日韩产的人妻AV在线网 | 中文字幕一区二区久久人妻网站 | 日本爱爱高清视频 | 久久久9久91精品夫妻电影 | 国产传媒三级 | 日本天天艹人人艹视频在线观看 | 成人av三级在线 成人AV在线导航 成人国产经典视频 | 欧美日逼逼 | 日韩无码免费高清 | 变态另类TS人妖一区二区 | 欧美成人网站在线 | 亚洲激情一区二区三区 | 中文字幕直接干 | 日本一二三区豆花视频 | 日韩A片网站 | 操逼王| 欧美老女人操逼 | 国产黄色小电影 | 伊人91在线 | 亚欧洲精品在线视频 | 逼综合网| 亚洲无码影院 | 国产一级乱伦视频 | 韩国一区二区三区在线观看 | 久久精品视频在线 | 看看亚洲的黄色网 | 青娱乐cao | 曰本国电影黄色免看费 | 男人先锋资源网 | 日韩人妻久久 | 三级网站看 | 日本成人在线午夜影院 | 毛片操逼 | 亚洲乱伦黄色小说视频网址 | 一区视频 | 婷婷久久综合激情综合 | 无码囯无精品毛片大码 | 欧美一区二区三区系列电影 | 亚洲电影av | 水蜜桃精品视频 | 蜜桃传媒一区二区亚洲AV | 先锋影音av网址导航 | 亚洲日韩一区二区 | 欧美亚洲三级 | 三级片免费观看网址 | 青娱乐在线播放在线观看视频 | 91精品老司机 | 中文字幕国产原创 | 国产精品v欧美精品v日韩精品 | 操动漫美女一区 | 欧美黄片网站 | 欧美性爱无码视频 | 青青草手机免费在线看片 | 一区二区网站 | 女人18毛片水多毛片久久1 | 蜜桃视频ht19.vp | 啪啪视频在线观看入口 | 国产又粗又大又黄视频 | 欧美久久久久久久 | 精品秘 无码一区二区三区老师 | 色婷婷丁香五月亚洲 | 日本成人免费电影一区二区三区 | 国产精品无人区 | 九九成人手机 | 大屌操视频在线 | 精品啪啪 | 色亭亭无码 | 91av视频在线 | av天天av无码av天天爽 Chinese国产人妖TS | 亚洲无圣光豆花 | 九九这里只有精品 | 五月激情婷婷中文字幕 | 高清无码一线逼美女系列 | 日韩黄色电影免费看 | 免费的一级A片 | 欧美日日 | 北条麻妃的无码视频 | 婷婷亚洲综合 | 99热在线观看精品 | 日本电影一区二区 | 波多野结衣黄色视频 | 91插穴| 久久加勒比 | 在线免费看不卡黄色视频 | 男女啪啪啪啪网站 | 成人电影在线观看A | 高清国产无码在线观看 | 俺也去了大香蕉在线 | 翔田千里中文字幕无码 | 精品无码一区二区三区 | 日韩三级视频网站 | caoporm超碰 | 亚洲无码高清在线播放 | 国产精品久久久久久苍井空 | 国产乱婬片视频 | 精品视频免费 | 91麻豆亚洲国产成人久久精品 | 波多野结衣乱码无码视频 | 蜜桃99视频一区二区三区 | 97激情| 特黄特色A| 日日干视频 | 亚洲无码国产乱伦 | 国产久久视频在线观看 | 亚洲一极黄色片 | 欧美在线视频网 | 欧美黑人XXXXX性受苍井空 | 亚洲中文成人娱乐网 | 俺来也俺就去www色情网 | 国产精品久久7777 | 中日韩免费视频 | 日本级婬乱片A片AAA毛片A | 亚洲乱码一二三四区国产 | 亚洲成人黄色视频 | MFYD-013 肉食人妻女上司が部下を誘惑し | 在线免费观看黄a | 天天舔天天插 | 中文日韩欧美 | www.老太色综合 | 成人艳史成人版 | 国产女人18毛片水真多果冻 | 影音先锋每日最新av | 日本激情视频网站免费 | 最新无码在线观看 | 牛牛视频在线观看 | 99曰高清热视频 | 国产蜜臀秘 入口 | 影音先锋AV网址 | 亚洲最新 | 超碰在线观看97 | 91视频观看 | 码人妻免费视频 | av大香 | BB在线视频网站 | 国产伦精品一区二区三区最新章节 | 麻豆国产91 在线播放猎赤 | 日本熟妇在线 | 日韩一级免费电影 | 人妻天天干 | 99 在线免费观看 | 欧美精品99 | 天天高清无码 | 亚洲无码三级视频 | 午夜123区 | 美女插屄视频 | 搡老女的BB视频 | 精品三级视频 | 人人免费操 | 国产一区二区三区四区五区入口 | 久久夜色AV | 亚洲无码精品九九九 | 操逼网络| 日韩黄色电影网址 | 成人在线免费网站 | 婷婷五月丁香激情免费视频 | 欧洲免费观看成品视频 | 国产毛片在线视频 | 色色色天堂 | 亚洲黄色视频在线播放 | 天天日日天天 | 国产操比小视频 | 日屁视频网站 | 精品福利在线观看 | 国产一a毛一a免费观看 | 九九九精品影视 | 色情小电影免费网站观看网址在线播 | 蜜桃成人综合 | 伊人网视频在线观看 | аⅴ天堂中文在线网 | 少妇久久精品 | 成人无码区免费AV毛片 | 欧美成人在线免费视频 | 婷婷五月色网 | 尻屄网站 | 免费网站看sm调教视频 | 字幕一区二区三区四区五区在线看 | 精品一区二区久久久久久久网站 | 欧美级一人在线视频播放 | 三级视频网站在线观看 | 欧美成人精品一区二区三区在线看 | 成年人视频在线播放 | 日本女人一区二区三区 | 欧美69av | 亚洲无码黄| 超碰啪啪啪 | 五月天婷婷久久 | www999女优館xxxcoh | 国产夫妻自拍性爱视频 | AV片网| 大鸡吧操小逼 | 国产黄色影院 | 婷婷综合网 | 天天日,天天干,天天射 | 俺也去五月婷婷 | 欧美aaa日韩aaa国产 | 亚洲男人的天堂网 | 逼逼爱插插网站 | 亚洲精品无码视频在线观看 | 麻豆视频免费在线观看 | 国产婷婷在线视频 | 欧美日韩亚洲视频 | 午夜影院中文字幕 | 国产精品稀缺盗摄盗拍 | 很很撸AV| 精品一区二区三区四区五区六区七区八区九区 | 亚洲无码成人在线 | 婷婷五月天激情小说 | 黄黄片在线看 | 久久噜噜噜久久熟女-久久久噜噜噜 | 美女操逼视频网 | 北条麻妃中文字幕久 | 在线久操| 伊人天堂网在线 | 国产一区黄片 | 欧美色图第五页 | 草草视频在线观看 | 肏屄视频看 | sm调教网站 | 日本无码高潮 | 手机无码视频在线观看 | 国产黄在线观看 | 强开小嫩苞A片黑人 | 亚洲成人在线影视 | 人人爱人人爽 | 欧美日屁| 午夜爽 | 国产一级a毛一级a做免费高清视频 | 黄色国产在线免费看 | 五月激情六月丁香 | 三级aV成人 | 日韩网站黄| 天天操天天干天天 | 91久久久久久 | 国产精品无码专区AV在线播放 | 高清无码视频在线 | 国产真实露脸乱子伦对白高清视频 | 免费看黄色片网站 | 五月天丁香激情 | 黄色成人在线观看 | 翔田千里无码在线 | 国产在线内射 | 天天躁日日躁狠狠躁免费麻豆 | 欧美天天搞 | 中文字字幕一区二区三区四区五区 | 三区无码视频 | 青青草原免费在线视频 | 久久亚洲欧美 | 国产一级a毛一级a看… | 久操视频在线免费观看 | 香蕉视频99 | 亚洲成人人人爽av电影 | 久久国产色 | 麻豆成人免费视频在线观看 | 亚洲资源站 | 在线黄| 日韩视频在线观看 | 国产极品人妖ts91热爆 | 亚洲无码性爱 | 亚洲第一色婷婷 | 久久久精品三级片 | 中文字幕日韩人妻久热 | 中文字字幕在线中文乱码更新时间 | 中国一级特黄毛片大片 | 99精品在线观看视频 | 国产黄色视频免费观看 | 超碰网站在线观看 | 影音先锋在线播放99av | 国产精品福利在线观看 | 一级黄色片免费看在干嘛呢 | 欧美大鸡巴操逼 | 一区二区三区四区视频在线 | 夜色av最新网址 一本无码免费视频 | 人妇一区三区高清乱码视频 | 韩国不卡视频 | 欧美黄色毛片 | 91豆花| 成人网站在线精品国产免费 | 芲井空αv无码一区二区三区 | 一级黄色毛片播放 | 国产三级电影网站 | 国内毛片毛片毛片 | www.av乱伦 | 青青草一区二区三区 | 韩国精品亚洲精品 | 91大庆哥 在线观看 | 极度兴感美女黄色操逼视频 | 免费69成人无码无遮又大 | 亚洲AV日韩AV高潮喷人人爽 | 五月婷久 | 韩国女主播操逼 | 中文字幕在线免费观看 | 日本精品在线播放 | 观看操逼视频 | 性爱视频在线无码播放 | 日韩欧美视频青青 | 日批视频免费播放 | 色综合五月天 | 黄色片免费视频网站 | 日韩高清在线 | 亚洲成人在线导航 | 国产综合久久7777777 | 欧美 自拍 视频 | 高清无码一线逼美女系列 | 久久99精品久久久久久草莓 | 青娱乐91超碰 | 凹凸精品一区二区三区四区在线 | 日韩无码影视 | 淫荡自拍 | 天天舔天天插天天爱 | 中国AAA一级操逼毛片 | 操逼图首页 | 人人爽夜夜爽 | 一级一级一级一级一级一级 | 69精品人人人 | 日日色综合 | 囯产精品久久久久久久久久久久 | 日韩综合天堂视频 | 色婷婷香蕉 | 日本不卡在线视频 | 8050午夜| 人人草超碰 | 国产人伦子伦一级A片下载 | 中文字幕无码人妻在线二区 | 国产麻豆系列 | 黄色欧美在线 | 国产操b电影网站在线观看 | 中文字幕无码在线视频 | 天天操夜夜操狠狠操 | 欧美成人激情 | 日韩无码毛片 | 日韩欧美一级A片久久久一区S图 | 久久国产精品久久久久久电车 | 欧美激情网 | 成人性生活片无码免费 | 五月天六月色婷婷在线 | 一区无码精品 | 青娱乐偷窥成人 | 大香蕉视频成人视频 | 爆操黑丝| 日本成人电影在线观看 | 日韩黄色电影院 | 国产一级二级 | 大香蕉婷婷| 国产视频一区在线 | 丝袜脚丫丁香五月 | 开心五月婷 | av综合人妻 | 亚洲中文成人娱乐网 | 亚洲无码在线精品 | 亚洲成色A片77777在线小说 | 青青草视频免费在线 | 亚洲日韩中文字幕视频在线 | 啊啊啊啊啊啊啊网站 | 91玉足脚交白嫩脚丫 | 成人91av| 神马影音先锋无码视频 | 免费看无码一级A片放24小时 | 91久草视频在线 | 玖玖视频在线观看 | 欧美黄色操逼视频 | 五月天色婷婷97在线视频播放 | 午夜av福利 | 欧美老熟妇xxx | 少妇的呻吟声 | 艹逼图| 黄色天天干 | 特黄av毛片 | 亚洲动漫一区 | 天天射夜夜爽 | 最黄国产操屄视频 | 亚洲色图欧美色图另类 | 国产精品久久久久久久鸭 | 媚药一区二区三区四区 | 免费无码又爽又高潮视频 | 美女被日视频网站 | 亚洲成人精品在线观看 | 成人毛片18毛片女人 | 奇米影视久久 | 久久你懂的| 人人妻人人爽人人操 | 色综合天天综合网天天狠天天 | 69久国产精品无码 | 亚洲综合成人在线视频 | 国产主播在线播放 | 大香蕉免费中文 | 久久麻豆av | 狠狠色综合网站 | 精品无码国产污污污免费网站 | 天天躁日日 | 亚洲黄色一级电影 | 日批在线视频 | A黄色视频网站 | 日本操逼片 | 噜噜吧噜噜久久综合 | 97超碰网 | 在线亚洲免费 | 久久久久无码国产精品一区 | 色情 视频网站入口第一集 | 中文字幕在线观看一区二区 | 国产AV综合播放 | 囯产精品久久久久久久久久98 | 成人国产精品视频 | 思思热免费 | 吖v在线观看 | 最新亚洲在线视频 | 福利国产在线 | 女人三级视屏 | 污污的啪啪网站 | 免费一级A片 | 91老夫老妻干老B | 亚洲天堂手机在线 | 亚洲黄色网扯 | 在线免费观看黄色电影 | 一级黄色片在线观看 | 欧美一级免费视频 | 超清无码在线网上 | 欧美午夜爱爱视频 | 成人欧美日韩 | 在线免费看无霉三级黄片 | 青青草日逼网 | 日韩高清中文字幕 | 欧美经典成人无码 | 能免费看AV的网站 | 天天添天天操 | 国产精品人成A片一区二区 | 狠狠人妻久久久久久综合蜜桃 | 青草青草久久 | 中国熟妇乱伦 | 久草热在线 | 伊人天天色狠狠 | 国产中文人人国际 | 超薄丝袜足交在线视频 | 欧美黄色高清视频免费 | 性爱午夜视频 | 亚洲青草视频 | 亚洲无码中文字幕在线观看视频 | 77777成人 | 国产精品偷窥熟女精品视 | 精品国产AV网站 | 国产无码激情后入 | 免费观看亚洲视频 | 国产AV资源| 欧美三级午夜理伦 | 日韩无码xxxx | 亚洲黄色在线网站 | 99色精品视频 | 在线毛片网站 | 国产精品麻豆三级三级视频 | 操婷婷逼 | 久久综合久色欧美综合狠狠 | 激情深爱五月 | 伊人五月婷婷丁香大香蕉 | 国产精品无码AV无码 | 欧美日韩一级免费电影 | 毛片A级成人片 | 国产在线8| 国产精品成人AV激情在线 | 爱骚逼黄色视频 | 国产亚洲一区二区三区 | 日日干夜夜爽 | 青青草在线免费 | 多人操穴视频在线播放 | 肏屄视频在线免费观看 | 国产v视频 | 一本色道久久88亚洲综合加勒比 | 久久亚洲视频 | 成人黄色性生活视频 | 奇米狠狠色777久久久欧美老妇 | 激情乱伦无码 | 精品在线视频网 | 中文娱乐在线视频 | 欧美操逼在线观看 | 日本黄色三级视频 | 最近中文字幕中文翻译歌词 | 欧美手机在线观看 | 贵族学校极品粉穴清纯学生妹 | 亚洲AV无码精品色午夜红一片 | 国产亚洲欧美日韩高清 | 日韩1页| 骚逼随便操 | 78m成人视频中文音声 | 四虎黄色按摩 | 色婷婷五月天在线观看 | 操逼操逼操逼操 | 欧美午夜成人视频 | 色se| 人人摸人人干人人射 | 日韩AV在线电影 | 精品国产a∨一区天美传媒 | 国产麻豆成人A | 国产插穴 | 艹逼胖子网站免费视频 | 高潮视频在线 | 午夜久久网 | 豆花视频在线看一区二区 | 爆乳一区新 | 三级黄,色毛片 | 大香蕉大香蕉视频影院 | 精品国产久久 | 色老板最新地址 | 啪啪啪视频免费观看 | 亚洲一级二级三级片 | 卡一卡二卡三在线 | 黄片免费在线播放 | 国内一区 | 午夜黄色视频 | 五月丁香六月久久 | 日韩无码一二三四区 | 国内精品久久久久久久久鸭下载 | 亚洲日韩AV电影 | 国产精品中文字幕在线观看 | 奇米久久778 | 欧美精品一卡 | 西西特级444www高清视频 | 小电影毛片 | 黄色录像带免费看 | 又粗又大又黄又爽无遮挡 | 欧美一级a一级a爱片免费 | 亚洲婷婷在线视频 | 精品视频一区二区 | 爱99区区 | 亚洲永久免费黄色电影 | 黄色免费视频网站 | 成人aV无码精品国产一区二区 | 夜色五月丁香久久 | 91豆花网站在线视频 | 丁香五月天婷婷激情 | 乱伦7777 | 成人日皮精品视频 | 污污网站在线免费观看 | 吃瓜偷拍AV一区二区三区 | 中文毛片 | www.蜜桃 | 色婷婷官网 | 91色在线视频 | 高清操逼视频 | 色播视频在线观看 | 91蜜臀在线视频免费 | 五月婷婷丁香网 | 久久久久久国产精品频道 | 免费靠逼视频 | 婷婷色在线 | 黄片在线免费网站 | 在线超碰成人 | 无码中文视频 | 9l视频自拍九色9l视频成人 | 人妻黄色一级 | 大香蕉伊人17c | 成人a一级毛片免费看 | 竹菊国产精品成人竹菊影视 | 操操操逼网 | 国产黄色视频在线观看免费 | 欧美色图片在线观看 | 久久人人妻人人做人人玩精品 | 欧美视频一区二区三区四区五区 | 日日干视频免费 | 黑人久久久 | 全日本爽视频在线观看 | 中文四区| 在线播放亚洲高清 | 欧美成人一级视频 | 青网草久草视频 | 国产成人无码精品色欲天香 | 午影操逼视频 | 国产欧美日韩在线视频 | 热99视频66 | 亚洲天堂影院 | 操必毛片| 无码一道本 | 亚洲欧洲一区二区 | 超碰爆乳| 一级A级黄色片 | 性欧美大战久久久久久久免费观看 | 欧美一级黄色丝袜大片免费 | 欧美三级又粗又硬 | 在线免费黄色视频 | 日韩成人午夜视频 | 天天射天天日天天 | 国产一级片免费看 | 天堂a区 天堂二区 | 国产青娱视频在线观看 | 天天干在线影院 | 免费无码成人片在线观看在线 | 猫咪www成人免费网站无码 | 高潮视频免费看 | 性爱免费视频网站 | www亚洲| 另类TS人妖一区二区三区牧瑶 | 欧美精品一区三区 | 少妇厨房愉情理伦BD在线观 | 免费无码在线 | 91三级成人片 | Av逼逼 | 美女黄片网站 | www.av黄 | 开心激情丁香网 | 久久成人理论电影手机 | 久日韩精品小视频 | 成人三级视频网站 | 麻豆精品三级电影 | 久久乐精品 | 男女污污污 | 十八女人毛片 | 人妻精品无码 | 五月婷婷亚洲综合 | 在线免费亚洲视频 | 99久久精品人妻无码一区二区蜜桃 | 污污网站在线观看 | 日韩久久久性爱 | 艹逼网| 亚洲香焦巨胸女人网视频 | 亚洲视频网站免费观看 | 天天摸天天摸天天 | www日本黄色 | 欧美屄视频| 最新一区二区在线 | 一区一二三视频 | 乱伦中文字幕网 | 一级特黄视频 | 欧美手机精品在线 | 在线aⅴ | 黄色天天影视 | 国产成人视频免费看 | 一本色道久久综合无码人妻 | 中文字幕欧美视频 | 日韩特黄一级片 | 日韩一类A片免费 | 夜夜爽妓女8888视频免费观看 | 毛片一区二区 | 嫩操影院 | 色老板精品视频在线观看 | 日韩无码第四页 | 国产宴妇精品久久久久久 | 色a视频 色操网站 | 欧美精品免费在线观看 | 国产一级一片免费播放 | 国产精品在线看 | 日本在线色视频 | 五月丁香乱伦 | 夜夜操屄| www.国产三级 | 影音先锋AV成人资源站在线播放 | 热久久综合免费频99热 | 欧美成人精品首页 | 在线免费a v | 国产三级精品在线 | 日本黄色一级代 | 精品自拍AV | 国产视频99在线观看 | 免费一区两区三区 | 色影音先锋色资源网站 | 日日撸色丁香五月 | 大操逼逼| 国产久久综合 | 小骚逼操死你视频 | 精品国产天线2024 | 日韩无码专区电影 | 夜夜高潮少妇 | 我要操西西| 欧美黄色电影在线 | 久久成人麻豆午夜电影 | 图片区偷拍区小说区 | 天天草天天日天天干 | 欧美日屄| 久久久午夜福利视频 | 围产精品久久久久久久久久久久 | 97男人的天堂 | 色国产综合免费视频在线播放 | 操逼 88AV | 操熟女91 | 啊啊啊啊操我在线动漫 | 极品美女操逼性 | 69视频网在线观看 | 欧美性猛交XXXX免费看久久久 | 最新一区二区在线 | 成人视频在线免费播放 | av天堂一区 | 色中色综合网 | 翔田千里无码XXXXXX | 影音先锋在线爱爱 | 小骚逼网| 成人综合在线视频 | а√天堂资源官网在线资源 | 日本中文字幕手机在线 | 日韩无码黄色电影 | 免费看欧美三级片 | 俺也来俺也去视频在线看 | 人人人人人人人人操 | av中文无码 | 日本激情视频网站免费 | 欧美成人在线性爱 | 欧美黄色小说 | 日韩AV中文字幕在线免费观看 | 天天艹天天干 | 欧美日韩性爱视频一区二区 | 8050午夜影院 | 人人爱人人摸人人舔 | 欧美性爱一区二区 | 婷婷丁香无码 | 国产女18毛片多18精品 | 亚洲一区二区自拍 | 无码真人操逼 | 亚洲AV无码变态另类在线播放 | 91久久婷婷亚洲精品成人 | 日日国产爽 | 亚洲无码免费观看高清 | 俺也不去网大香蕉 | 射一射在线视频 | 嫩草久久99www亚洲红桃 | 老妇人一区二区三区 | 亚洲成人AV| 午夜精品久久久久久不卡8050 | 北条麻妃一区二区三区 | 欧美操逼网| 台湾午夜成人节目在线播放 | 88AV中文字幕教师 | 夜夜躁狠狠躁日日躁视频 | 青久久| 亚洲人成色777777精品音频 | 日屄视频免费观看 | 婷婷丁香伊人 | 国产精品综合 | 国产女人18水真多18精品 | 丁香激情五月 | 一区二区免费在线 | 人人澡在线视频 | 国产乱婬片A片AAAAPp地址 | 蜜逃视频观看 | 在线色福利 | 草逼网站免费观看 | 黑人精品无码 | 国产成人视频免费看 | 日本黄色一级片日 | 围内精品久久久久久久久久‘变脸 91久久婷婷国产麻豆精品电影 | 亚洲第一中文字幕 | 最近中文字幕在线中文字幕7 | 欧美色图亚洲色图另类 | 性一交一乱一A片久久99蜜桃 | 性高潮视频网站 | 成年人在线视频 | 四虎91| 视频亚洲无码 | 国产成人视频在线播放 | 亚洲成人网站无码在线观看 | 乱伦AV中文字幕 | 国产精品一哟哟 | 日韩中文字幕无码中字字幕 | 亚洲www啪成人一区二区麻豆 | av中文字幕网 | 亚洲粉嫩18p | 黄色一级一级直播间 | doi在线观看 | 男女趴趴网站 | 欧美乱妇日本无乱码特黄大片 | 国产一级乱伦片 | 天天日天天摸天天操 | 夜间福利无码视频无码 | 69精品无码成人久久久久久 | 操丝袜美女的骚逼 | 亚洲天堂婷婷 | 亚洲无码 一道本 | 操肥逼视频| 色先锋av天天操 少妇 后入 视频 少妇精品一区二区 | 日韩精品视频网 | 免费三级成人爱做片 | 精品成人久久久久久久 | 日韩成人免费视频 | 小早川怜子爆乿护士中文 | 欧美成人超碰在线 | 一级黄色电影免费 | 销魂少妇一二三区 | 亚洲AV无码成人专区 | 高清无码在线观看免费 | 淫秽视频在线免费观看 | 豆花无码视频 | 学生妹一级J人片内射视频 | 香蕉毛片 | 很黄很色免费视频 | 亚洲黄色在线视频 | 大香蕉伊人在线观看视频 | 在线豆花Av | 三级网站永久大全 | 五月天Av成人在线播放 | 久久人人爽人人爽人人片aV东京热 | 婷婷五月丁香激情免费视频 | 黑人狂躁女人高潮视频 | 欧美操B视频 | 日韩精品无码一区二区 | 黄色在线播放视频免费 | TS人妖一区二区三区视频 | 国产上婬乱18一级毛 | 大香蕉怡春院日韩精品视频 | www.大香蕉综合网 | 免费在线看的黄片 | 中文字幕一区二区三区四区 | 女人毛片 | 國產精品77777777777 | AV国产在线观看 | 免费黄色在线看 | 蜜桃视频成人 免费 | 性爱777 | 国产免费靠逼视频 | 国产高清无码视频在线播放 | jiZZJIZZ成熟丰满少妇 | 夜色视频网| 亚洲无码五区 | 五月丁香花 | 最新亚洲无码在线观看 | 欧美日屄 | 色色伊人网 | 五月天激情婷婷婷久久 | 操日本美女 | 青青自拍视频免费观看 | 亚洲系列第一页 | 国产免费自拍视频 | 亚洲熟妇性ⅩXXX交潮喷 | 黄色操逼的免费网站 男女操逼的视频免费网站 | 亚洲自拍偷拍视频 | 小蔡头喵喵喵-新年小奶牛 | 美女张开腿让男人捅视频 | 伊人综合操逼网 | 老熟女亂伦88AV | 大香蕉在线啪啪 | 国产一级a毛一级a做免费图 | av婷婷免费 | 天天免费摸天天免费日 | 精品国模 | 操丝袜美女的骚逼 | 国产一级视频在线观看 | 第一页欧美 | ww视频在线观看 |