• <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ì)的研究
    露出奶头的视频| 亚洲色图av天堂| 久久久久九九精品影院| 国产精品免费一区二区三区在线| 欧美亚洲日本最大视频资源| 女同久久另类99精品国产91| 中文字幕av电影在线播放| 国产精品 欧美亚洲| 亚洲成人久久爱视频| 久久亚洲真实| 久久久久久大精品| 少妇被粗大的猛进出69影院| 亚洲成人免费电影在线观看| 久久 成人 亚洲| 校园春色视频在线观看| 在线观看日韩欧美| 大型av网站在线播放| 变态另类丝袜制服| avwww免费| 亚洲精品久久成人aⅴ小说| 草草在线视频免费看| 亚洲精品久久成人aⅴ小说| 久久久国产欧美日韩av| 女人爽到高潮嗷嗷叫在线视频| 亚洲欧美激情综合另类| 欧美日韩中文字幕国产精品一区二区三区| 伊人久久大香线蕉亚洲五| 国产精品电影一区二区三区| 国产日本99.免费观看| а√天堂www在线а√下载| 黄色女人牲交| 亚洲 欧美 日韩 在线 免费| 亚洲人成网站高清观看| 久久精品成人免费网站| 亚洲成人国产一区在线观看| 国产高清视频在线播放一区| 人人妻人人看人人澡| 久久久久国产一级毛片高清牌| 亚洲精华国产精华精| 国产av一区二区精品久久| 国产99白浆流出| 啦啦啦观看免费观看视频高清| 精品熟女少妇八av免费久了| 99国产综合亚洲精品| 国产高清激情床上av| 国产精品亚洲一级av第二区| 国产一级毛片七仙女欲春2 | 白带黄色成豆腐渣| 中文亚洲av片在线观看爽| 国产精品久久久久久亚洲av鲁大| 国产真人三级小视频在线观看| 国产高清激情床上av| 国产精品爽爽va在线观看网站 | 欧美中文日本在线观看视频| 成人18禁高潮啪啪吃奶动态图| 精品高清国产在线一区| 欧美黄色淫秽网站| 18美女黄网站色大片免费观看| 久久久国产成人免费| 桃色一区二区三区在线观看| 亚洲欧洲精品一区二区精品久久久| 欧美一级a爱片免费观看看 | 亚洲精品在线观看二区| 热re99久久国产66热| 天天一区二区日本电影三级| 侵犯人妻中文字幕一二三四区| 色哟哟哟哟哟哟| 91成人精品电影| 亚洲激情在线av| 亚洲av日韩精品久久久久久密| cao死你这个sao货| 国产精品美女特级片免费视频播放器 | 一级a爱视频在线免费观看| 无遮挡黄片免费观看| 最好的美女福利视频网| 90打野战视频偷拍视频| 亚洲自拍偷在线| 身体一侧抽搐| 99久久99久久久精品蜜桃| 亚洲一码二码三码区别大吗| 国产免费av片在线观看野外av| 天天躁狠狠躁夜夜躁狠狠躁| 韩国av一区二区三区四区| x7x7x7水蜜桃| 亚洲aⅴ乱码一区二区在线播放 | 久久久久亚洲av毛片大全| 激情在线观看视频在线高清| 精品久久久久久久久久免费视频| 国产亚洲欧美精品永久| bbb黄色大片| 精品久久久久久久毛片微露脸| 一a级毛片在线观看| 一本综合久久免费| 美女扒开内裤让男人捅视频| 香蕉国产在线看| 久久人妻福利社区极品人妻图片| 国产欧美日韩一区二区三| 亚洲美女黄片视频| 国产精品一区二区精品视频观看| 高清在线国产一区| av片东京热男人的天堂| 男女视频在线观看网站免费 | 亚洲avbb在线观看| 嫁个100分男人电影在线观看| 亚洲av日韩精品久久久久久密| 91麻豆精品激情在线观看国产| 国产高清视频在线播放一区| 日本黄色视频三级网站网址| 国产黄片美女视频| 中文字幕人妻熟女乱码| 日本免费a在线| 亚洲片人在线观看| 久久欧美精品欧美久久欧美| 国产视频一区二区在线看| 欧美乱码精品一区二区三区| 欧美色欧美亚洲另类二区| 99热6这里只有精品| 久久这里只有精品19| 免费搜索国产男女视频| 免费一级毛片在线播放高清视频| 日韩高清综合在线| 国产一区二区在线av高清观看| 天天躁狠狠躁夜夜躁狠狠躁| 熟妇人妻久久中文字幕3abv| 成人18禁高潮啪啪吃奶动态图| 丁香六月欧美| 男男h啪啪无遮挡| 一进一出抽搐动态| 十八禁网站免费在线| 亚洲自拍偷在线| 色av中文字幕| 精品久久久久久久末码| 在线观看一区二区三区| 久热爱精品视频在线9| 国产精品永久免费网站| 免费在线观看成人毛片| 国产日本99.免费观看| 亚洲一区高清亚洲精品| 国产精品 国内视频| a在线观看视频网站| 女性生殖器流出的白浆| 欧美午夜高清在线| 国产又色又爽无遮挡免费看| 欧美av亚洲av综合av国产av| 神马国产精品三级电影在线观看 | 女人被狂操c到高潮| 可以在线观看的亚洲视频| 黄片大片在线免费观看| 母亲3免费完整高清在线观看| а√天堂www在线а√下载| 久久久国产精品麻豆| 久久香蕉国产精品| 国产黄a三级三级三级人| 免费av毛片视频| 男女下面进入的视频免费午夜 | 日韩欧美三级三区| 18禁国产床啪视频网站| 国产蜜桃级精品一区二区三区| 精华霜和精华液先用哪个| 人人妻人人澡欧美一区二区| 伊人久久大香线蕉亚洲五| 老司机靠b影院| 精品国产超薄肉色丝袜足j| 成人特级黄色片久久久久久久| 久热爱精品视频在线9| 国产视频一区二区在线看| 一级黄色大片毛片| 黄网站色视频无遮挡免费观看| 婷婷精品国产亚洲av在线| 国产成人精品久久二区二区免费| 国产成人av激情在线播放| 法律面前人人平等表现在哪些方面| 亚洲第一欧美日韩一区二区三区| 免费在线观看亚洲国产| 精品无人区乱码1区二区| 国产成人一区二区三区免费视频网站| 国产色视频综合| 这个男人来自地球电影免费观看| 老司机午夜十八禁免费视频| 午夜精品在线福利| 国产熟女午夜一区二区三区| 黄色片一级片一级黄色片| 夜夜夜夜夜久久久久| 人人妻,人人澡人人爽秒播| 夜夜看夜夜爽夜夜摸| 亚洲欧美精品综合一区二区三区| 日韩精品中文字幕看吧| 欧美国产精品va在线观看不卡| 欧美日韩瑟瑟在线播放| www日本在线高清视频| 亚洲va日本ⅴa欧美va伊人久久| 亚洲一区二区三区不卡视频| 国产亚洲精品一区二区www| 超碰成人久久| 少妇熟女aⅴ在线视频| 国产精品国产高清国产av| 美女国产高潮福利片在线看| 亚洲色图av天堂| 国内少妇人妻偷人精品xxx网站 | 亚洲,欧美精品.| 免费av毛片视频| 久久香蕉国产精品| 一级片免费观看大全| 亚洲av熟女| 制服诱惑二区| 精品无人区乱码1区二区| 欧美一级毛片孕妇| 成人18禁高潮啪啪吃奶动态图| 美国免费a级毛片| 侵犯人妻中文字幕一二三四区| 中文亚洲av片在线观看爽| 成人av一区二区三区在线看| 国产亚洲精品第一综合不卡| 国产av一区二区精品久久| 亚洲第一青青草原| 变态另类成人亚洲欧美熟女| 免费在线观看视频国产中文字幕亚洲| av免费在线观看网站| 岛国在线观看网站| 精品国产亚洲在线| 搡老岳熟女国产| 88av欧美| 日本免费a在线| 日韩大码丰满熟妇| 大香蕉久久成人网| 好男人在线观看高清免费视频 | 99国产精品99久久久久| www日本在线高清视频| 999精品在线视频| www日本黄色视频网| 变态另类丝袜制服| 久久亚洲真实| 俺也久久电影网| 免费在线观看黄色视频的| 97超级碰碰碰精品色视频在线观看| 成人三级做爰电影| 黄片播放在线免费| 波多野结衣巨乳人妻| 女人爽到高潮嗷嗷叫在线视频| 久久久久免费精品人妻一区二区 | 悠悠久久av| 精品电影一区二区在线| 日本免费一区二区三区高清不卡| 女人高潮潮喷娇喘18禁视频| 精品国产一区二区三区四区第35| 免费av毛片视频| 成人亚洲精品一区在线观看| 最近最新中文字幕大全免费视频| 亚洲精品久久国产高清桃花| 国产欧美日韩一区二区精品| 久久久国产成人精品二区| 中文在线观看免费www的网站 | 制服人妻中文乱码| 法律面前人人平等表现在哪些方面| 最近最新中文字幕大全电影3 | 欧美 亚洲 国产 日韩一| 动漫黄色视频在线观看| 国产精品香港三级国产av潘金莲| 他把我摸到了高潮在线观看| 久久香蕉国产精品| 麻豆av在线久日| 男人操女人黄网站| 国产精品 欧美亚洲| 亚洲国产欧美网| 精品国产国语对白av| 亚洲一区中文字幕在线| 国产视频一区二区在线看| 97碰自拍视频| 一区福利在线观看| 久久草成人影院| 真人一进一出gif抽搐免费| 国产欧美日韩一区二区精品| 美女 人体艺术 gogo| 色尼玛亚洲综合影院| 午夜日韩欧美国产| 亚洲色图av天堂| 伊人久久大香线蕉亚洲五| 男人舔女人下体高潮全视频| 中文字幕另类日韩欧美亚洲嫩草| 亚洲人成电影免费在线| 精品第一国产精品| 亚洲精品一卡2卡三卡4卡5卡| 久久久久久九九精品二区国产 | 色老头精品视频在线观看| 日本黄色视频三级网站网址| 搡老岳熟女国产| 国产一区二区在线av高清观看| 久久天躁狠狠躁夜夜2o2o| 亚洲欧美激情综合另类| 亚洲欧美日韩高清在线视频| 欧美成人一区二区免费高清观看 | 一卡2卡三卡四卡精品乱码亚洲| 欧洲精品卡2卡3卡4卡5卡区| 热re99久久国产66热| 亚洲人成网站在线播放欧美日韩| 亚洲美女黄片视频| 长腿黑丝高跟| 国产成人精品无人区| 国产国语露脸激情在线看| 亚洲人成网站高清观看| 美国免费a级毛片| 国产精品影院久久| 91成年电影在线观看| 亚洲人成伊人成综合网2020| 久久伊人香网站| 搞女人的毛片| 国产精品久久久久久亚洲av鲁大| 日本精品一区二区三区蜜桃| 最近最新免费中文字幕在线| 首页视频小说图片口味搜索| 一级毛片高清免费大全| 久久精品91蜜桃| 日韩精品免费视频一区二区三区| 叶爱在线成人免费视频播放| 日韩有码中文字幕| 黄片播放在线免费| 国产亚洲欧美98| 国产精品野战在线观看| 日韩精品中文字幕看吧| 中文字幕人成人乱码亚洲影| 美女午夜性视频免费| 亚洲成av片中文字幕在线观看| 99久久无色码亚洲精品果冻| 一本久久中文字幕| 在线播放国产精品三级| 激情在线观看视频在线高清| 精华霜和精华液先用哪个| 亚洲精华国产精华精| 亚洲成国产人片在线观看| 免费在线观看视频国产中文字幕亚洲| 国产精品亚洲一级av第二区| 婷婷精品国产亚洲av在线| 99在线人妻在线中文字幕| 欧美人与性动交α欧美精品济南到| 人妻丰满熟妇av一区二区三区| 国产亚洲av高清不卡| 久久精品夜夜夜夜夜久久蜜豆 | 国产成人精品久久二区二区免费| 久久精品91蜜桃| 久久99热这里只有精品18| 男女之事视频高清在线观看| 非洲黑人性xxxx精品又粗又长| 最近最新中文字幕大全电影3 | 国产三级在线视频| 亚洲精品粉嫩美女一区| 日韩高清综合在线| 欧美中文综合在线视频| 两性午夜刺激爽爽歪歪视频在线观看 | 国产又爽黄色视频| 国产亚洲精品久久久久5区| 免费在线观看亚洲国产| 在线国产一区二区在线| 97碰自拍视频| 免费高清在线观看日韩| 首页视频小说图片口味搜索| 色婷婷久久久亚洲欧美| 人人妻,人人澡人人爽秒播| 成人免费观看视频高清| 99在线视频只有这里精品首页| 国内精品久久久久精免费| 变态另类成人亚洲欧美熟女| 久久天堂一区二区三区四区| 国产成人精品久久二区二区91| 两人在一起打扑克的视频| 色播在线永久视频| 久久精品人妻少妇| 久久香蕉国产精品| 久久久国产成人免费| 日韩大尺度精品在线看网址| 精品国产超薄肉色丝袜足j| 18禁国产床啪视频网站| 免费看十八禁软件| 欧洲精品卡2卡3卡4卡5卡区| 成年版毛片免费区| 亚洲第一av免费看| 久久精品夜夜夜夜夜久久蜜豆 | 欧美性猛交╳xxx乱大交人| 午夜激情福利司机影院| 免费在线观看日本一区| 亚洲精品一区av在线观看| 国产精品日韩av在线免费观看| 久99久视频精品免费| 亚洲国产日韩欧美精品在线观看 | 一级黄色大片毛片| 国产一卡二卡三卡精品| www.www免费av| 国产成人av激情在线播放| 12—13女人毛片做爰片一| 久久久久久国产a免费观看| 精品欧美国产一区二区三| 精品国产乱子伦一区二区三区| 不卡一级毛片| a级毛片a级免费在线| 久久草成人影院| 欧美 亚洲 国产 日韩一| 久久精品aⅴ一区二区三区四区| 亚洲七黄色美女视频| 久久久久国产一级毛片高清牌| 中文字幕精品免费在线观看视频| 最好的美女福利视频网| 91av网站免费观看| 久久精品国产99精品国产亚洲性色| 搡老熟女国产l中国老女人| 美女午夜性视频免费| a在线观看视频网站| 99精品在免费线老司机午夜| av天堂在线播放| 69av精品久久久久久| 久久伊人香网站| 成人国产综合亚洲| 国产在线观看jvid| 中国美女看黄片| 久久青草综合色| 亚洲自拍偷在线| 亚洲七黄色美女视频| 亚洲国产中文字幕在线视频| 国产亚洲av高清不卡| АⅤ资源中文在线天堂| 国产日本99.免费观看| 国产精品乱码一区二三区的特点| 一级a爱视频在线免费观看| 中文在线观看免费www的网站 | 黑人操中国人逼视频| 91成人精品电影| 国产伦人伦偷精品视频| 欧美日韩亚洲国产一区二区在线观看| 中文字幕最新亚洲高清| 侵犯人妻中文字幕一二三四区| 女性生殖器流出的白浆| 国产高清videossex| 精品人妻1区二区| av超薄肉色丝袜交足视频| 国产在线精品亚洲第一网站| 91字幕亚洲| 国产精品自产拍在线观看55亚洲| 国产97色在线日韩免费| 女性被躁到高潮视频| 免费电影在线观看免费观看| 桃红色精品国产亚洲av| 亚洲av中文字字幕乱码综合 | 午夜免费激情av| 精品第一国产精品| 久久久久久九九精品二区国产 | 亚洲男人天堂网一区| 热99re8久久精品国产| 一本大道久久a久久精品| 国产私拍福利视频在线观看| 久久香蕉激情| 国产伦在线观看视频一区| 亚洲国产看品久久| av福利片在线| 久久精品亚洲精品国产色婷小说| 1024香蕉在线观看| 亚洲免费av在线视频| 一级a爱视频在线免费观看| 国产欧美日韩精品亚洲av| 国产视频内射| 国产精品乱码一区二三区的特点| 亚洲成人国产一区在线观看| 1024手机看黄色片| 国产精品乱码一区二三区的特点| 国产欧美日韩一区二区三| 国产极品粉嫩免费观看在线| 亚洲国产日韩欧美精品在线观看 | 日本免费一区二区三区高清不卡| 搡老妇女老女人老熟妇| 日韩欧美国产在线观看| 伊人久久大香线蕉亚洲五| 久久婷婷人人爽人人干人人爱| 视频在线观看一区二区三区| 久久 成人 亚洲| 美女大奶头视频| 国产熟女午夜一区二区三区| 欧美精品啪啪一区二区三区| 黄色女人牲交| 性欧美人与动物交配| 老司机靠b影院| 99国产精品一区二区蜜桃av| 日日摸夜夜添夜夜添小说| 草草在线视频免费看| 久99久视频精品免费| 欧美性长视频在线观看| 亚洲国产欧洲综合997久久, | 久久婷婷成人综合色麻豆| 俄罗斯特黄特色一大片| 最近最新中文字幕大全电影3 | 天堂动漫精品| 国产精品香港三级国产av潘金莲| 久久久国产成人免费| 岛国在线观看网站| 99在线人妻在线中文字幕| 成年女人毛片免费观看观看9| 少妇熟女aⅴ在线视频| 在线永久观看黄色视频| 黄片大片在线免费观看| 欧美乱色亚洲激情| 波多野结衣高清作品| 欧美日韩亚洲国产一区二区在线观看| 18禁观看日本| www.www免费av| 欧美国产精品va在线观看不卡| 国产一区二区在线av高清观看| 国产久久久一区二区三区| 人人妻人人澡欧美一区二区| 香蕉国产在线看| 精品福利观看| 18禁美女被吸乳视频| 成人一区二区视频在线观看| 国产欧美日韩一区二区三| 成熟少妇高潮喷水视频| 神马国产精品三级电影在线观看 | 国产精品自产拍在线观看55亚洲| 久久午夜综合久久蜜桃| 午夜影院日韩av| 久久香蕉激情| 国产精品1区2区在线观看.| av在线天堂中文字幕| 怎么达到女性高潮| 草草在线视频免费看| 桃色一区二区三区在线观看| 国产精品亚洲一级av第二区| 国产一区二区在线av高清观看| 丰满人妻熟妇乱又伦精品不卡| 人人妻人人澡欧美一区二区| 欧美色视频一区免费| 国产精品99久久99久久久不卡| 天天添夜夜摸| 免费观看精品视频网站| 国产蜜桃级精品一区二区三区| 日本撒尿小便嘘嘘汇集6| 老汉色∧v一级毛片| 日韩高清综合在线| 白带黄色成豆腐渣| 久久人妻福利社区极品人妻图片| 人妻丰满熟妇av一区二区三区| 日本 av在线| 久久久久国产一级毛片高清牌| 久久欧美精品欧美久久欧美| 精品欧美国产一区二区三| av免费在线观看网站| 狠狠狠狠99中文字幕| 亚洲国产日韩欧美精品在线观看 | 男女做爰动态图高潮gif福利片| 免费看日本二区| 最好的美女福利视频网| 亚洲人成电影免费在线| 亚洲精品国产区一区二| 视频在线观看一区二区三区| 婷婷亚洲欧美| 欧美成人性av电影在线观看| 国产成+人综合+亚洲专区| 91成年电影在线观看| 久久 成人 亚洲| 国产一区二区在线av高清观看| 精品久久久久久久久久免费视频| 一区二区日韩欧美中文字幕| 老司机午夜福利在线观看视频| 男人舔奶头视频| 黄色片一级片一级黄色片| 亚洲成人久久性| 两性夫妻黄色片| 精品国产超薄肉色丝袜足j| 91国产中文字幕| 国产视频一区二区在线看| av超薄肉色丝袜交足视频| 亚洲国产精品合色在线| 欧美日韩亚洲综合一区二区三区_| 精品高清国产在线一区| 女生性感内裤真人,穿戴方法视频| 日韩欧美国产在线观看| 国产蜜桃级精品一区二区三区| √禁漫天堂资源中文www| 亚洲片人在线观看| 两人在一起打扑克的视频| 熟女电影av网| 欧美成人一区二区免费高清观看 | 99久久无色码亚洲精品果冻| 亚洲国产精品合色在线| 日本免费a在线| 欧美绝顶高潮抽搐喷水| 国产精品av久久久久免费| 韩国av一区二区三区四区| 黑丝袜美女国产一区| 久久久久免费精品人妻一区二区 | 成人手机av| 亚洲电影在线观看av| 在线观看免费日韩欧美大片| 亚洲成人免费电影在线观看| 丝袜人妻中文字幕| 在线观看66精品国产| 久久精品成人免费网站| 午夜激情av网站| 欧美黑人巨大hd| 一二三四在线观看免费中文在| 久久亚洲真实| 国产精品98久久久久久宅男小说| 69av精品久久久久久| 欧美不卡视频在线免费观看 | 91成年电影在线观看| 午夜日韩欧美国产| 嫩草影院精品99| 欧美乱码精品一区二区三区| av在线天堂中文字幕| 成人永久免费在线观看视频| 日本免费a在线| 757午夜福利合集在线观看| 黄频高清免费视频| www日本在线高清视频| 亚洲黑人精品在线| 变态另类成人亚洲欧美熟女| 中国美女看黄片| 黄网站色视频无遮挡免费观看| 禁无遮挡网站| 悠悠久久av| 国产成人一区二区三区免费视频网站|