• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    一種三對角矩陣的相似對角化及其應(yīng)用

    2021-01-13 05:04張興剛曹磊
    關(guān)鍵詞:特征向量特征值

    張興剛 曹磊

    基金項(xiàng)目:國家自然科學(xué)基金資助項(xiàng)目(11965007);貴州省科技合作計劃資助項(xiàng)目(20157641)

    摘 要:針對二維筒倉中隨機(jī)堆積的顆粒體系,建立了一個由吸收系數(shù)p和側(cè)向傳遞系數(shù)q決定的格點(diǎn)系統(tǒng)模型。導(dǎo)出了筒倉底部的壓力分布與頂部壓力分布的關(guān)系式,該關(guān)系式的解析計算涉及到三對角形式的傳遞系數(shù)矩陣A(p,q)的相似對角化。通過相似變換將A(p,q)對稱化,并且提出基于二階差分方程的方法討論這種對稱三對角矩陣的特征值和特征向量,只有在特殊情況下才能解析地求出A(p,q)的所有特征值并將其相似對角化。最后將相關(guān)的數(shù)學(xué)結(jié)果用于計算有效質(zhì)量隨顆??傎|(zhì)量的變化關(guān)系,結(jié)果表明格點(diǎn)系統(tǒng)模型可以得到與Janssen模型類似的結(jié)果。

    關(guān)鍵詞:糧倉效應(yīng);格點(diǎn)模型;三對角矩陣;特征值;特征向量

    中圖分類號:O151,O312? 文獻(xiàn)標(biāo)志碼:A

    顆粒介質(zhì)的靜力學(xué)問題是顆粒物質(zhì)研究中基礎(chǔ)且重要的方面,它與實(shí)際應(yīng)用密切相關(guān)(例如,工程中地基的承載能力、砂石骨料的受力結(jié)構(gòu),材料科學(xué)中顆粒材料的力學(xué)性能,等等),對于其它復(fù)雜體系(例如,懸浮液、泡沫、膠體等)的理解也有重要的意義[1-2]。人們在很早的時候就發(fā)現(xiàn),筒倉中的糧食與水在靜力學(xué)方面有很不相同的表現(xiàn),糧倉底部所受壓力會隨糧食堆積高度的增加而趨于飽和,這就是顆粒物質(zhì)靜力學(xué)中一個著名的效應(yīng)——糧倉效應(yīng)[1,3]。早在1895年,德國工程師Janssen就提出了基于一定假設(shè)的連續(xù)介質(zhì)模型,對糧倉效應(yīng)進(jìn)行解釋。然而,顆粒物質(zhì)是一種典型的復(fù)雜體系,目前還沒有關(guān)于顆粒介質(zhì)靜力學(xué)的系統(tǒng)、成熟的理論,許多基本問題仍需深入探究。因此,自上世紀(jì)九十年代的顆粒物質(zhì)研究熱潮以來,仍有許多研究者對糧倉效應(yīng)及相關(guān)的問題進(jìn)行研究[4-13]。例如,最近有關(guān)反糧倉效應(yīng)[11]和磁性顆粒體系糧倉效應(yīng)[12]的工作。已有的研究工作以實(shí)驗(yàn)和計算機(jī)模擬為主,理論方面也主要是對Janssen模型、OSL模型等宏觀唯象模型[4]的討論,對于糧倉效應(yīng)產(chǎn)生的微觀機(jī)理以及唯象參數(shù)與微觀結(jié)構(gòu)的聯(lián)系還不是很清楚。為了深入認(rèn)識顆粒介質(zhì)所服從的力學(xué)規(guī)律,很有必要加強(qiáng)顆粒體系的微觀結(jié)構(gòu)、力網(wǎng)、力傳遞等問題的研究,基于微觀的幾何和力學(xué)分析建立體系的微觀特征與宏觀性質(zhì)之間的聯(lián)系[13-15]。

    對于圓盤的晶格堆積,文獻(xiàn)[6]中通過結(jié)構(gòu)和力學(xué)平衡的分析得到了相鄰兩層顆粒間力傳遞的規(guī)則,并且利用力傳遞規(guī)則得到了隨著堆積層數(shù)增加容器底部總壓力會趨向飽和的結(jié)果。對于非晶格的堆積結(jié)構(gòu),不借助計算機(jī)模擬很難從理論上直接進(jìn)行結(jié)構(gòu)和力學(xué)的分析。標(biāo)量q模型是研究顆粒堆積中力的幾率分布的一個著名的模型[15]。為了便于從理論上研究筒倉底部壓力分布與頂部壓力分布的關(guān)系,可以采用q模型的基本思想,將筒倉中隨機(jī)堆積的顆粒體系抽象簡化為格點(diǎn)系統(tǒng)。本文以二維筒倉中隨機(jī)堆積的顆粒體系為物理背景,建立了一個由吸收系數(shù)p和側(cè)向傳遞系數(shù)q決定的格點(diǎn)系統(tǒng)模型,該模型的解析求解涉及到三對角形式的傳遞系數(shù)矩陣A(p,q)的相似對角化。文中通過相似變換將A(p,q)對稱化,并且提出基于二階差分方程的方法對矩陣的特征值和特征向量進(jìn)行分析和計算;最后,將相關(guān)的數(shù)學(xué)結(jié)果用于討論糧倉效應(yīng)的有關(guān)問題。

    1 力傳遞的數(shù)學(xué)模型

    為討論問題方便,這里主要研究二維的情形,將糧倉效應(yīng)中的物理體系抽象為圖1所示的帶有邊界的格點(diǎn)系統(tǒng)。在格點(diǎn)系統(tǒng)中,一個格點(diǎn)等效于一小團(tuán)顆粒,第m行n列的格點(diǎn)用符號P(m,n)表示,格點(diǎn)的總行數(shù)M≥2,總列數(shù)N≥5,并且每個格點(diǎn)都受到重力m的作用。類似于標(biāo)量q模型,這里主要關(guān)注每個作用力在豎直方向上的分量,因此在下面的討論中,所說的作用力一般是指力在豎直方向上分力的大小。頂部各個格點(diǎn)所受的壓力可用行向量f=(f, f,…,f)描述,底部各個格點(diǎn)對容器底的壓力用f=(f,f,…,f)描述。為了給出底部壓力分布f與頂部壓力分布f的關(guān)系,需要結(jié)合實(shí)際物理體系的特征定義恰當(dāng)?shù)南噜弮蓪痈顸c(diǎn)間力傳遞的規(guī)則。這里以隨機(jī)堆積為物理背景,定義了圖1所示的一種典型的力傳遞規(guī)則。規(guī)定同一行的格點(diǎn)間沒有豎直方向的分力,不相鄰的兩行格點(diǎn)間也沒有相互作用。對于任意格點(diǎn)P(m,n),稱其所受的重力、頂部對其產(chǎn)生的壓力、以及第m-1行格點(diǎn)對其產(chǎn)生的作用力為輸入力;將P(m,n)對第m+1行格點(diǎn)、及其對容器的作用力稱為輸出力;對于每一個格點(diǎn),只有上述五種豎直方向分力的存在,為了保證格點(diǎn)在豎直方向的力平衡,總輸入力必須等于總輸出力。圖1中用實(shí)線箭頭比較完整地畫出了第一行格點(diǎn)的各個輸入力及輸出力的情況(重力沒有畫),圖中的虛線箭頭上也標(biāo)出了每種輸出力對應(yīng)的系數(shù)。

    是一個三對角矩陣[16]。傳遞系數(shù)矩陣只由p和q這兩個參量確定,因此為了表述方便,稱這個模型為pq模型。模型中的p稱為吸收系數(shù),反映了由于容器側(cè)壁與顆粒的摩擦所導(dǎo)致的格點(diǎn)總輸入力被器壁吸收的程度,它滿足0≤p<1;若容器側(cè)壁與顆粒無摩擦,則p=0,有摩擦?xí)r0? ? ?由圖1中定義的力傳遞規(guī)則,有fM=wM;因此由(1)、(2)兩式可以導(dǎo)出fM與f0的關(guān)系。假設(shè)存在可逆矩陣Q使A=QΛQ-1,并且假設(shè)A-E可逆(E是N階單位矩陣),不難導(dǎo)出

    利用Janssen模型[1,4]可以得到Me=Ms(1-e-Mt/Ms),其中Ms叫作飽和質(zhì)量;顯然(33)式的結(jié)果與Janssen模型得到的結(jié)果很接近,這說明本文中的pq模型也可以解釋糧倉效應(yīng),它在一定程度上正確反映了筒倉系統(tǒng)中力傳遞的特征。利用(29)式,還可以比較方便地討論容器底部的應(yīng)力分布、點(diǎn)載荷的響應(yīng)等方面的力學(xué)問題。當(dāng)然,如果要深入了解p和q對這些力學(xué)問題的影響,就需要求解方程(25);根據(jù)前面的數(shù)學(xué)討論,在一般情況下應(yīng)該結(jié)合數(shù)值計算才能解決這個問題。

    4 結(jié)論

    為了從微觀上理解糧倉效應(yīng),我們針對二維隨機(jī)堆積建立了一個格點(diǎn)系統(tǒng)模型——pq模型。該模型的解析求解涉及到三對角形式的傳遞系數(shù)矩陣A(p,q)的相似對角化,這也是本文的重點(diǎn)。我們通過相似變換將A(p,q)對稱化,并且提出基于二階差分方程的方法比較有效地解決了該對稱矩陣的特征值和特征向量的計算。對于p=1/2、q=1這種典型的情況,我們解析地導(dǎo)出了糧倉底部的有效質(zhì)量隨顆粒總質(zhì)量的變化關(guān)系;其它p,q取值的情況也可轉(zhuǎn)化為數(shù)值求解高次代數(shù)方程(25)的問題。本文討論的pq模型只是格點(diǎn)系統(tǒng)模型中一種可能的力傳遞形式,對于不同的堆積結(jié)構(gòu)和環(huán)境條件,可能具有不同的力傳遞規(guī)則,在數(shù)學(xué)上一般表現(xiàn)為不同的傳遞系數(shù)矩陣。本文提出的理論方法對于其它傳遞系數(shù)矩陣的情形有借鑒的價值,研究的結(jié)果也能加深對微觀力傳遞規(guī)則與宏觀壓力分布之間關(guān)系的認(rèn)識。參考文獻(xiàn):

    [1]陸坤權(quán), 劉寄星. 顆粒物質(zhì)(上)[J]. 物理, 2004, 33(9): 629-635.

    [2] 孫其誠, 厚美瑛, 金峰, 等. 顆粒物質(zhì)物理與力學(xué)[M]. 北京: 科學(xué)技術(shù)出版社, 2011.

    [3] DE-GENNES P G. Granular matter: a tentative view[J]. Rev Mod Phys, 1999,71(2):374-382.

    [4] VANEL L, CLAUDIN P, BOUCHAUD J P, et al. Stresses in silos: comparison between theoretical models and new experiments[J]. Phys Rev Lett, 2000, 84(7): 1439-1442.

    [5] BERTHO Y, FRDRIQUE G D, HULIN J P. Dynamical Janssen effect on granular packing with moving walls[J]. Phys Rev Lett, 2003, 90(14): 144301.

    [6] 張興剛, 胡林, 隆正文. 糧倉效應(yīng)的圓盤堆積模型[J]. 計算物理, 2006, 23(6): 642-646.

    [7] WAMBAUGH J F, HARTLEY R R, BEHRINGER R P. Force networks and elasticity in granular silos[J]. Eur Phys J E, 2010, 32(2): 135.

    [8] 李智峰, 彭政, 蔣亦民. 糧倉內(nèi)顆粒壓力的測量:Janssen行為及其偏差[J]. 物理學(xué)報, 2014, 63(10): 292-299.

    [9] 宗謹(jǐn), 周志剛, 王文廣, 等. 顆粒固體應(yīng)力轉(zhuǎn)向比的光彈法探測[J]. 物理學(xué)報, 2017, 66(10): 177-184.

    [10]BLANCO-RODRGUEZ R, PREZ-NGEL G. Stress distribution in two-dimensional silos[J]. Phys Rev E, 2018, 97(1): 12903.

    [11]MAHAJAN S, TENNENBAUM M, PATHAK S N, et al. Reverse janssen effect in narrow granular columns[J]. Phys Rev Lett, 2020, 124(12): 128002.

    [12]THORENSL, MLY K J, BOURGOIN M, et al. Magnetic Janssen effect[J]. Nat Commun, 2021, 12(1):2486.

    [13]ZHANG X G, DAI D. Exact solutions for stress distribution in rhombic disk packings[J]. Granular Matter, 2021, 23: 51.

    [14]GENDELMAN O, POLLACK Y G, PROCACCIA I. What determines force chains in granular media?[J]. Phys Rev Lett ,2016, 116(2): 078001.

    [15]COPPERSMITH S N, LIU C, MAJUMDAR S, et al. Model for force fluctuations in bead packs[J]. Phys Rev E, 1996, 53(5): 4673-4685.

    [16]楊勝良. 三對角矩陣的特征值及其應(yīng)用[J]. 數(shù)學(xué)的實(shí)踐與認(rèn)識, 2010, 40(3): 155-160.

    (責(zé)任編輯:于慧梅)

    作者簡介:張興剛(1980—),男,副教授,博士,研究方向:統(tǒng)計物理和凝聚態(tài)理論的研究,E-mail:xgzhang@gzu.edu.cn.

    通訊作者:張興剛,E-mail:xgzhang@gzu.edu.cn.

    Similarity Diagonalization of a Tridiagonal Matrix and its Applications

    ZHANG Xinggang CAO Lei

    (1.School of Physics,Guizhou University,Guiyang 550025,China;

    2.Powerchina Guiyang Engineering Co. Ltd., Guiyang 550081,China)

    Abstract: A lattice model which is determined by the absorption coefficient p and the lateral transfer coefficient q is proposed for the two dimensional disordered granular columns. The pressure distribution relation between the top and the bottom of the container is deduced, and the analytic calculation of this mathematical relation involves the similarity diagonalization of tridiagonal transfer matrix A(p,q). A similarity transformation is proposed to symmetrize A(p,q), and the method based on the second order difference equations is proposed to obtain the eigenvalues and eigenvectors of the symmetric matrix. Finally, these mathematical results are used to discuss Janssen effect. Calculations show that results similar to Janssen model can be produced by our lattice model.

    Key words: janssen effect; lattice model; tridiagonal matrix; eigenvalue; eigenvector

    猜你喜歡
    特征向量特征值
    高中數(shù)學(xué)特征值和特征向量解題策略
    三個高階微分方程的解法研究
    求矩陣特征值的一個簡單方法
    利用“降階法”求解歐拉方程
    一道全國大學(xué)生數(shù)學(xué)競賽試題的推廣
    球殼區(qū)域上二階橢圓特征值問題的一種高精度數(shù)值逼近
    氨基酸序列特征向量提取方法的探討
    基于鼠標(biāo)行為的電子商務(wù)中用戶異常行為檢測
    一種改進(jìn)的峰均功率比判源方法
    一類非線性矩陣方程組性質(zhì)的研究
    日本黄大片高清| 听说在线观看完整版免费高清| 久久热精品热| 少妇熟女aⅴ在线视频| 久久99热这里只频精品6学生| av免费在线看不卡| av网站免费在线观看视频 | 亚洲国产精品专区欧美| 免费看不卡的av| 秋霞在线观看毛片| 亚洲精品日本国产第一区| 少妇裸体淫交视频免费看高清| av免费观看日本| 亚洲真实伦在线观看| 精品熟女少妇av免费看| 伦理电影大哥的女人| 亚洲无线观看免费| 97超碰精品成人国产| 日韩视频在线欧美| 高清在线视频一区二区三区| 非洲黑人性xxxx精品又粗又长| 春色校园在线视频观看| 亚洲精品乱码久久久久久按摩| 亚洲国产精品国产精品| 免费观看无遮挡的男女| 九九久久精品国产亚洲av麻豆| 在线天堂最新版资源| 国产单亲对白刺激| 建设人人有责人人尽责人人享有的 | 久久午夜福利片| 在线观看美女被高潮喷水网站| 亚洲欧洲日产国产| 日本猛色少妇xxxxx猛交久久| av又黄又爽大尺度在线免费看| 男女啪啪激烈高潮av片| 亚洲av日韩在线播放| 国产伦在线观看视频一区| 欧美xxxx性猛交bbbb| 午夜老司机福利剧场| 亚洲欧美精品自产自拍| 日韩欧美国产在线观看| 一本一本综合久久| 亚洲真实伦在线观看| 特级一级黄色大片| 51国产日韩欧美| 精品久久久精品久久久| 国产精品一区二区三区四区久久| 在线观看一区二区三区| 91av网一区二区| 五月伊人婷婷丁香| 精品人妻视频免费看| 人人妻人人澡欧美一区二区| 亚洲精品自拍成人| 国产黄色免费在线视频| 欧美日韩在线观看h| 亚洲第一区二区三区不卡| 看免费成人av毛片| 久久久久免费精品人妻一区二区| 国产亚洲精品av在线| 中文字幕久久专区| 国产精品伦人一区二区| 国产69精品久久久久777片| 久久精品国产亚洲网站| 免费黄色在线免费观看| 成年人午夜在线观看视频 | 国产高清有码在线观看视频| 少妇裸体淫交视频免费看高清| 综合色丁香网| 啦啦啦中文免费视频观看日本| 看非洲黑人一级黄片| 亚洲国产高清在线一区二区三| 精品久久久久久久人妻蜜臀av| 欧美一级a爱片免费观看看| 日本黄色片子视频| 男人和女人高潮做爰伦理| 美女内射精品一级片tv| 久久午夜福利片| 99久久精品国产国产毛片| 久久久色成人| 亚洲人与动物交配视频| 日韩av在线大香蕉| av在线天堂中文字幕| 国产欧美另类精品又又久久亚洲欧美| 日韩人妻高清精品专区| 91久久精品电影网| 欧美变态另类bdsm刘玥| 秋霞伦理黄片| 91狼人影院| 日韩欧美一区视频在线观看 | 亚洲av二区三区四区| 国产淫语在线视频| 亚洲人成网站高清观看| 国产乱人偷精品视频| 欧美成人午夜免费资源| 老师上课跳d突然被开到最大视频| 欧美高清性xxxxhd video| 久久久久久久久久久丰满| 免费黄频网站在线观看国产| 国产日韩欧美在线精品| 免费人成在线观看视频色| 久久久欧美国产精品| 亚洲成人久久爱视频| 嫩草影院入口| 国产一区二区亚洲精品在线观看| 久久热精品热| 搞女人的毛片| 久久久精品免费免费高清| 亚洲人与动物交配视频| 如何舔出高潮| 女人被狂操c到高潮| 99热这里只有是精品在线观看| 国产69精品久久久久777片| 国产精品久久久久久久久免| 日本三级黄在线观看| 99热这里只有是精品50| 久久这里有精品视频免费| 亚洲av成人精品一二三区| 国产爱豆传媒在线观看| 精品国产一区二区三区久久久樱花 | 成年人午夜在线观看视频 | 午夜福利视频1000在线观看| 听说在线观看完整版免费高清| 久久久国产一区二区| 亚洲av成人精品一二三区| 97人妻精品一区二区三区麻豆| 一区二区三区高清视频在线| 男的添女的下面高潮视频| 国产高清不卡午夜福利| 黄色一级大片看看| 成年版毛片免费区| 十八禁网站网址无遮挡 | 美女高潮的动态| 简卡轻食公司| 国产v大片淫在线免费观看| 伊人久久国产一区二区| 免费人成在线观看视频色| 亚洲,欧美,日韩| 看十八女毛片水多多多| 国内少妇人妻偷人精品xxx网站| 亚洲电影在线观看av| 欧美高清成人免费视频www| 国产色爽女视频免费观看| 夜夜爽夜夜爽视频| 国产免费又黄又爽又色| 尾随美女入室| 久久久久久久久中文| 色5月婷婷丁香| 欧美最新免费一区二区三区| 两个人的视频大全免费| 免费看a级黄色片| 少妇被粗大猛烈的视频| 亚州av有码| 亚洲精品国产成人久久av| 嫩草影院新地址| 男插女下体视频免费在线播放| 国产精品美女特级片免费视频播放器| 晚上一个人看的免费电影| 99久久精品国产国产毛片| 久久国内精品自在自线图片| 欧美另类一区| 久久精品夜夜夜夜夜久久蜜豆| freevideosex欧美| 在线天堂最新版资源| 国产亚洲精品av在线| 国产伦理片在线播放av一区| 中文资源天堂在线| 女人久久www免费人成看片| 日本免费在线观看一区| 免费看a级黄色片| 亚洲国产av新网站| 久久精品综合一区二区三区| 日日摸夜夜添夜夜爱| 国产片特级美女逼逼视频| 婷婷色av中文字幕| 国产 一区 欧美 日韩| 亚洲精华国产精华液的使用体验| av专区在线播放| 又黄又爽又刺激的免费视频.| 国产亚洲精品av在线| 亚洲色图av天堂| 亚洲精品成人av观看孕妇| 亚洲国产成人一精品久久久| 久久国内精品自在自线图片| 美女黄网站色视频| 丰满人妻一区二区三区视频av| 国产黄a三级三级三级人| 一二三四中文在线观看免费高清| 久久综合国产亚洲精品| 噜噜噜噜噜久久久久久91| 激情 狠狠 欧美| 97在线视频观看| 国产精品精品国产色婷婷| 日本免费a在线| 国产免费视频播放在线视频 | av黄色大香蕉| 又爽又黄a免费视频| 久久久a久久爽久久v久久| 国产成人精品一,二区| 秋霞伦理黄片| 免费看不卡的av| 精品少妇黑人巨大在线播放| 国产精品一区二区在线观看99 | 亚洲精品亚洲一区二区| 色吧在线观看| 国精品久久久久久国模美| 国产成人aa在线观看| 亚洲不卡免费看| 亚洲国产日韩欧美精品在线观看| 99视频精品全部免费 在线| 亚洲激情五月婷婷啪啪| 少妇被粗大猛烈的视频| 久久久成人免费电影| 国产三级在线视频| 国产精品av视频在线免费观看| 亚洲国产欧美人成| 日韩欧美精品免费久久| 99热全是精品| 国产v大片淫在线免费观看| 免费少妇av软件| 99久久九九国产精品国产免费| 中文字幕av成人在线电影| 国产精品av视频在线免费观看| 亚洲最大成人av| 97热精品久久久久久| 一级av片app| 日韩欧美一区视频在线观看 | 国产成人一区二区在线| 国产 一区 欧美 日韩| a级一级毛片免费在线观看| 午夜久久久久精精品| 精品久久久噜噜| 国产 一区 欧美 日韩| 亚洲精品一区蜜桃| 天堂av国产一区二区熟女人妻| 国产 一区 欧美 日韩| 欧美变态另类bdsm刘玥| 黄色一级大片看看| 亚洲欧美一区二区三区黑人 | 免费黄网站久久成人精品| 欧美高清成人免费视频www| 永久网站在线| 日韩制服骚丝袜av| 国产精品女同一区二区软件| 男插女下体视频免费在线播放| 亚洲欧美一区二区三区国产| 特大巨黑吊av在线直播| 精品一区二区三卡| 欧美最新免费一区二区三区| 日本-黄色视频高清免费观看| 日日啪夜夜撸| 日韩伦理黄色片| 又粗又硬又长又爽又黄的视频| 一级二级三级毛片免费看| 国产免费视频播放在线视频 | 成年版毛片免费区| 1000部很黄的大片| 久久国产乱子免费精品| 久久久久久久久久黄片| 2021少妇久久久久久久久久久| 国产日韩欧美在线精品| 日本黄大片高清| 97在线视频观看| 五月天丁香电影| 亚洲av电影不卡..在线观看| 精品久久久久久久久久久久久| 午夜精品一区二区三区免费看| www.av在线官网国产| 国产高潮美女av| 亚洲一级一片aⅴ在线观看| 国产高清国产精品国产三级 | 中文字幕av在线有码专区| 天堂中文最新版在线下载 | 全区人妻精品视频| 男人舔女人下体高潮全视频| 一级二级三级毛片免费看| 99热全是精品| 亚洲人成网站在线观看播放| 免费观看av网站的网址| 波多野结衣巨乳人妻| 你懂的网址亚洲精品在线观看| 国产欧美另类精品又又久久亚洲欧美| 熟妇人妻久久中文字幕3abv| 三级经典国产精品| 成人国产麻豆网| 亚洲av电影在线观看一区二区三区 | 一个人看视频在线观看www免费| 成年女人看的毛片在线观看| 哪个播放器可以免费观看大片| 综合色av麻豆| 婷婷色综合www| 真实男女啪啪啪动态图| 在线a可以看的网站| 97人妻精品一区二区三区麻豆| 精品久久久久久久久av| 一二三四中文在线观看免费高清| 免费看日本二区| 人妻一区二区av| 亚洲av国产av综合av卡| 亚洲欧洲日产国产| 国产片特级美女逼逼视频| 免费黄频网站在线观看国产| 国产黄片美女视频| 久久精品夜夜夜夜夜久久蜜豆| 超碰av人人做人人爽久久| 亚洲精品第二区| 99久久精品一区二区三区| 久久久久久九九精品二区国产| 亚洲欧美成人综合另类久久久| av免费在线看不卡| 亚洲最大成人av| 免费观看性生交大片5| 欧美xxxx黑人xx丫x性爽| 99热全是精品| 亚洲精品国产av蜜桃| 自拍偷自拍亚洲精品老妇| 亚洲国产精品sss在线观看| av免费在线看不卡| 毛片一级片免费看久久久久| 在线免费观看不下载黄p国产| 欧美一级a爱片免费观看看| 大话2 男鬼变身卡| 床上黄色一级片| 国产成人aa在线观看| 亚洲av中文字字幕乱码综合| 成年女人看的毛片在线观看| av免费在线看不卡| 久久国内精品自在自线图片| 国产色爽女视频免费观看| av网站免费在线观看视频 | 亚洲国产欧美人成| 99热这里只有是精品在线观看| 狂野欧美白嫩少妇大欣赏| 国产在线一区二区三区精| 免费看光身美女| 欧美日韩精品成人综合77777| 九九久久精品国产亚洲av麻豆| a级毛片免费高清观看在线播放| av卡一久久| 欧美日韩国产mv在线观看视频 | 22中文网久久字幕| 亚洲自偷自拍三级| 国模一区二区三区四区视频| 特大巨黑吊av在线直播| 亚洲国产日韩欧美精品在线观看| 亚洲精品乱码久久久久久按摩| 亚洲内射少妇av| 禁无遮挡网站| 中文字幕免费在线视频6| 特级一级黄色大片| 日日啪夜夜撸| 少妇的逼好多水| 能在线免费看毛片的网站| 久久久久久久国产电影| 国产伦在线观看视频一区| 午夜福利网站1000一区二区三区| 七月丁香在线播放| 中文乱码字字幕精品一区二区三区 | 尾随美女入室| 热99在线观看视频| 一级毛片我不卡| 国产成人精品一,二区| 亚洲精品自拍成人| 亚洲aⅴ乱码一区二区在线播放| 久久久久久久大尺度免费视频| 黄色配什么色好看| 在线观看一区二区三区| 校园人妻丝袜中文字幕| 熟妇人妻不卡中文字幕| 色网站视频免费| 一级毛片 在线播放| 国产精品嫩草影院av在线观看| 欧美一级a爱片免费观看看| 老司机影院毛片| 九九久久精品国产亚洲av麻豆| 插阴视频在线观看视频| 国产色婷婷99| 久久精品人妻少妇| 国产黄色小视频在线观看| 日本三级黄在线观看| 亚洲成人av在线免费| 色综合亚洲欧美另类图片| 高清欧美精品videossex| 免费观看精品视频网站| 国产乱人视频| 国产综合精华液| 精品99又大又爽又粗少妇毛片| 日本免费在线观看一区| av在线观看视频网站免费| 五月玫瑰六月丁香| 毛片一级片免费看久久久久| 干丝袜人妻中文字幕| 秋霞伦理黄片| 日本-黄色视频高清免费观看| 男女下面进入的视频免费午夜| 一级二级三级毛片免费看| 日韩欧美国产在线观看| 国产乱人视频| 亚洲综合色惰| 日韩制服骚丝袜av| 99久久精品国产国产毛片| 伦精品一区二区三区| 久久久久久久久久久丰满| 国产精品不卡视频一区二区| www.色视频.com| 日韩 亚洲 欧美在线| 赤兔流量卡办理| 人妻夜夜爽99麻豆av| 床上黄色一级片| 亚洲一区高清亚洲精品| 亚洲在久久综合| 少妇被粗大猛烈的视频| 性色avwww在线观看| 精品一区二区三卡| 亚洲乱码一区二区免费版| 欧美激情在线99| av在线老鸭窝| 久久99精品国语久久久| 一区二区三区四区激情视频| 亚洲图色成人| 天天躁夜夜躁狠狠久久av| 天美传媒精品一区二区| 观看免费一级毛片| 国产乱人视频| 亚洲精品影视一区二区三区av| 午夜免费观看性视频| 欧美日韩精品成人综合77777| 一区二区三区乱码不卡18| 亚洲丝袜综合中文字幕| 亚洲乱码一区二区免费版| videossex国产| 青春草视频在线免费观看| 日日啪夜夜撸| 两个人视频免费观看高清| 嘟嘟电影网在线观看| 精品国产露脸久久av麻豆 | 亚洲熟妇中文字幕五十中出| 99热6这里只有精品| av在线亚洲专区| 亚洲欧美日韩无卡精品| 丰满少妇做爰视频| 中文字幕亚洲精品专区| 亚洲精品日韩在线中文字幕| 在线观看美女被高潮喷水网站| 色尼玛亚洲综合影院| freevideosex欧美| 亚洲四区av| 日日摸夜夜添夜夜添av毛片| 久久99热这里只有精品18| 亚洲人与动物交配视频| 九草在线视频观看| 九色成人免费人妻av| 国产在线男女| 日日啪夜夜撸| 午夜福利在线观看吧| 一级爰片在线观看| 免费观看a级毛片全部| 美女cb高潮喷水在线观看| 日韩成人伦理影院| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 午夜福利在线观看免费完整高清在| 国产精品人妻久久久影院| 韩国高清视频一区二区三区| 老司机影院毛片| 乱人视频在线观看| 老师上课跳d突然被开到最大视频| 国产精品福利在线免费观看| 国模一区二区三区四区视频| 免费人成在线观看视频色| 久久99精品国语久久久| 亚洲高清免费不卡视频| 久久亚洲国产成人精品v| 久久这里只有精品中国| 国产乱来视频区| 国产高潮美女av| 久久韩国三级中文字幕| 免费看不卡的av| 国产 亚洲一区二区三区 | 男插女下体视频免费在线播放| 国产伦在线观看视频一区| 五月玫瑰六月丁香| 欧美极品一区二区三区四区| 在线观看av片永久免费下载| 欧美日韩国产mv在线观看视频 | 激情五月婷婷亚洲| 日日啪夜夜爽| 亚洲av国产av综合av卡| 日本爱情动作片www.在线观看| 天堂网av新在线| 亚洲国产精品sss在线观看| 在线免费十八禁| 看十八女毛片水多多多| 中国国产av一级| 老师上课跳d突然被开到最大视频| 少妇人妻一区二区三区视频| 一本一本综合久久| 99久久九九国产精品国产免费| 男女国产视频网站| 真实男女啪啪啪动态图| 欧美bdsm另类| 成年免费大片在线观看| 纵有疾风起免费观看全集完整版 | 久久人人爽人人片av| 大陆偷拍与自拍| 综合色丁香网| 成人一区二区视频在线观看| 美女国产视频在线观看| 嘟嘟电影网在线观看| freevideosex欧美| 国产精品嫩草影院av在线观看| 看十八女毛片水多多多| 听说在线观看完整版免费高清| 欧美97在线视频| 久久久久久久久久黄片| 久久精品国产亚洲网站| 亚洲成色77777| 听说在线观看完整版免费高清| 99热这里只有精品一区| 国产爱豆传媒在线观看| 亚洲欧美日韩卡通动漫| 精品人妻一区二区三区麻豆| 人妻夜夜爽99麻豆av| 一个人免费在线观看电影| 真实男女啪啪啪动态图| 只有这里有精品99| 国产在视频线精品| 色5月婷婷丁香| 国产精品三级大全| 精品99又大又爽又粗少妇毛片| 蜜臀久久99精品久久宅男| 男人和女人高潮做爰伦理| av天堂中文字幕网| 国产免费视频播放在线视频 | 中国国产av一级| 久久亚洲国产成人精品v| 热99在线观看视频| av一本久久久久| 精品久久久噜噜| 中文精品一卡2卡3卡4更新| 精品久久久噜噜| 亚洲自偷自拍三级| 日韩欧美精品免费久久| 国产在视频线在精品| 最新中文字幕久久久久| 中文字幕制服av| 五月伊人婷婷丁香| 久久99热6这里只有精品| 午夜日本视频在线| 人妻少妇偷人精品九色| 永久免费av网站大全| 午夜福利视频1000在线观看| 波野结衣二区三区在线| 久久6这里有精品| 国内精品一区二区在线观看| 国产精品女同一区二区软件| 看黄色毛片网站| 日韩强制内射视频| 秋霞伦理黄片| 又粗又硬又长又爽又黄的视频| 高清在线视频一区二区三区| 美女主播在线视频| 亚洲精品中文字幕在线视频 | 女的被弄到高潮叫床怎么办| 亚洲av日韩在线播放| 国产在视频线精品| 成人亚洲欧美一区二区av| 91在线精品国自产拍蜜月| 在线播放无遮挡| 国国产精品蜜臀av免费| 黄色欧美视频在线观看| 日韩三级伦理在线观看| 精品人妻一区二区三区麻豆| 日韩精品有码人妻一区| 成人二区视频| av播播在线观看一区| 国产精品国产三级国产专区5o| 男的添女的下面高潮视频| 国语对白做爰xxxⅹ性视频网站| 亚洲av二区三区四区| 丝袜美腿在线中文| 黄色配什么色好看| 最近中文字幕高清免费大全6| 伦精品一区二区三区| 美女xxoo啪啪120秒动态图| 亚洲人成网站在线观看播放| 国产精品日韩av在线免费观看| 秋霞伦理黄片| 久久精品综合一区二区三区| 99久久精品一区二区三区| 国产国拍精品亚洲av在线观看| 国产片特级美女逼逼视频| 国产免费福利视频在线观看| 午夜精品在线福利| 亚洲精品aⅴ在线观看| 亚洲国产色片| 男女啪啪激烈高潮av片| 97精品久久久久久久久久精品| 男女边摸边吃奶| 久久精品国产亚洲网站| 日产精品乱码卡一卡2卡三| 日韩亚洲欧美综合| 成年女人看的毛片在线观看| 性色avwww在线观看| 插逼视频在线观看| 久久国内精品自在自线图片| 国产免费福利视频在线观看| 麻豆成人午夜福利视频| 久久热精品热| 人妻系列 视频| 国产 一区精品| 久久久久久久大尺度免费视频| 日韩中字成人| 国产久久久一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 直男gayav资源| 国产成年人精品一区二区|