日韩av黄I国产麻豆传媒I国产91av视频在线观看I日韩一区二区三区在线看I美女国产在线I麻豆视频国产在线观看I成人黄色短片

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

怎么作非线性拟合_值得收藏ansys非线性收敛总结

發布時間:2024/9/30 编程问答 50 豆豆
生活随笔 收集整理的這篇文章主要介紹了 怎么作非线性拟合_值得收藏ansys非线性收敛总结 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

ansys計算非線性時會繪出收斂圖,其中橫坐標是cumulative iteration number 縱坐標是absolute convergence norm。他們分別是累積迭代次數和絕對收斂范數,用來判斷非線性分析是否收斂。

ansys在每荷載步的迭代中計算非線性的收斂判別準則和計算殘差。其中計算殘差是所有單元內力的范數,只有當殘差小于準則時,非線性疊代才算收斂。

ansys的收斂是基于力的收斂的,以力為基礎的收斂提供了收斂量的絕對值,而以位移為基礎的收斂僅提供表現收斂的相對量度。一般不單獨使用位移收斂準則,否則會產生一定偏差,有些情況會造成假收斂.(ansys非線性分析指南--基本過程Page.6)?。因此ansys官方建議用戶盡量以力為基礎(或力矩)的收斂誤差,如果需要也可以增加以位移為基礎的收斂檢查。

ANSYS缺省是用L2范數控制收斂。其它還有L1范數和L0范數,可用CNVTOL命令設置。在計算中L2值不斷變化,若L2

由于ANSYS缺省的criterion計算是全部變量的平方和開平方(SRSS)*valuse(你設置的值),所以crition也有小小變化。如有需要,也可自己指定crition為某一常數, CNVTOL,F,10000,0.0001,0就指定力的收斂控制值為10000*0.0001=1。

另外,非線性計算中用到的一個開關是SOLCONTROL,如關閉SOLCONTROL 選項,那么軟件默認收斂準則:力或彎矩的收斂容差是0.001,而不考慮位移的收斂容差;如果打開SOLCONTROL 選項,同樣的默認收斂準則:力或彎矩的收斂容差是0.005,而位移收斂容差是0.05。

非線性收斂非常麻煩,與網格精度、邊界條件、荷載步等一系列因素有關,單元的特點對收斂的影響很大,單元的性態不好收斂則困難些;合理的步長可以使求解在真解周圍不至于振蕩,步長過小,計算量太大,步長過大,會由于過大的荷載步造成不收斂。網格密度適當有助于收斂,網格太密計算量太大,當然太稀計算結果會有較大的誤差。究竟多少往往要針對問題進行多次試算。

1、如果不收斂,可以考慮一下方法改進

1.放松非線性收斂準則。

(CNVTOL #Sets convergence values for nonlinear analyses).

2.增加荷載步數。

(NSUBST #Specifies the number of substeps to be taken this load step)

3.增加每次計算的迭代次數(默認的25次)

(NEQIT #Maximum number of equilibrium iterations allowed each substep)

4 重新劃分單元試試,后續會得到不同的答案。

就我的經驗,前兩種方法效果好一些,后一種方法不一定奏效。但應注意到,放大收斂準則其實是在降低計算精度的條件下得到近似解,并且放大的收斂準則是否與實際相符或有實際意義應仔細考慮。

2、問非線性計算的收斂和速度

Q:我在計算一個大型結構,地震荷載,BEAM188 計算時間太長一個小時可能計算了1秒總共40秒而且越來越慢,不小心早上還停了電如何能使計算加快?或者怎么才能即使突然結束以后還能繼續算?

A:調整優化非線性計算的收斂和速度可以說幾乎是一種藝術, 即沒有固定的可循規則。

我的經驗是, 你的結構的"非線性"越小, 非線性的變化越規則, 就越容易收斂. 想象一下如果你是手算這個非線性問題, 對你來講較容易的, 對ANSYS的相應算法也會容易些.

可以把你的地震時程分析拿出幾點, 做一下靜態的非線性分析, 同時調整模型看看分析出來的結果是否合理. 如果這一步還沒有做, 那花大量時間做出的時程分析是廢品的可能性十分之大.

一定要記住有限元分析是一個"簡化"問題的過程. 建立一個模型一定要由淺到深. 線性的模型沒有搞透不要貿然進攻非線性. 靜態沒有搞透不要碰時程分析.

A:影響非線性收斂穩定性及其速度的因素很多,我們可以看看這幾點:

a、影響非線性收斂穩定性及其速度的因素:

模型——主要是結構剛度的大小。

對于某些結構,從概念的角度看,我們可以認為它是幾何不變的穩定體系。但如果結構相近的幾個主要構件剛度相差懸殊,或者懸索結構的索預應力過小(即它的剛度不夠大),在數值計算中就可能導致數值計算的較大誤差,嚴重的可能會導致結構的幾何可變性——忽略小剛度構件的剛度貢獻。如果還不能理解,我們可以進一步說:我們有一種通用的方法判斷結構的幾何可變性,即det(K)=0。在數值計算中,要得到det(K)恒等于零是不可能的,我們也就只能讓它較小時即認為結構是幾何可變的。對于上述的結構,他們的K值是很小的,故而也可判斷為幾何可變體系。事實上這類結構在實際工程中 也的確是非常危險的。

為此,我們要看看模型有沒有問題。如出現上述的結構,要分析它,就得降低剛度很大的構件單元的剛度,可以加細網格劃分,或著改用高階單元(BEAM->SHELL,SHELL->SOLID)。構件的連接形式(剛接或鉸接)等也可能影響到結構的剛度。

線性算法(求解器)

ANSYS中的非線性算法主要有:稀疏矩陣法(SPARSE DIRECT SOLVER)、預共軛梯度法(PCG SOLVER)和波前法(FRONT DIRECT SLOVER)。

稀疏矩陣法是性能很強大的算法,一般默認即為稀疏矩陣法(除了子結構計算默認波前法外)。預共軛梯度法對于3-D實體結構而言是最優的算法,但當結構剛度呈現病態時,迭代不易收斂。為此推薦以下算法:

1)BEAM單元結構,SHELL單元結構,或以此為主的含3-D SOLID的結構,用稀疏矩陣法;

2)3-D SOLID的結構,用預共軛梯度法;

3)當你的結構可能出現病態時,用稀疏矩陣法;

4)當你不知道用什么時,可用稀疏矩陣法。

3)非線性逼近技術

在ANSYS里還是牛頓-拉普森法和弧長法。牛頓-拉普森法是我們常用的方法,收斂速度較快,但也和結構特點和步長有關。弧長法常被某些人推崇備至,它能算出力加載和位移加載下的響應峰值和下降響應曲線。但也發現:在峰值點,弧長法仍可能失效,甚至在非線性計算的線性階段,它也可能會無法收斂。

為此,我們盡量不要從開始即激活弧長法,還是讓程序自己激活為好(否則出現莫名其妙的問題)。子步(時間步)的步長還是應適當,自動時間步長也是很有必要的。

b.如何加快計算速度

在大規模結構計算中,計算速度是一個非常重要的問題。下面就如何提高計算速度作一些建議:

充分利用ANSYS MAP分網和SWEEP分網技術,盡可能獲得六面體網格,這一方面減小解題規模,另一方面提高計算精度。在生成四面體網格時,用四面體單元而不要用退化的四面體單元。比如95號單元有20節點,可以退化為10節點四面體單元,而92號單元為10節點單元,在此情況下用92號單元將優于95號單元。

選擇正確的求解器。對大規模問題,建議采用PCG法。此法比波前法計算速度要快10倍以上(前提是您的計算機內存較大)。對于工程問題,可將ANSYS缺省的求解精度從1E-8改為1E-4或1E-5即可。

????設置直接影響到收斂。應該注意以下幾點:

1)足夠大的荷載步(將MAXMIUM SUBSTEP=1000000),可以更容易收斂,避免發散的出現(nsub,nsbstp,nsbmx,nsbmn);

2)足夠大的平衡迭代步數,默認為25,可以放大到很大(100)(eqit,eqit);

3)將收斂準則調整,以位移控制時調整為0.05,以力控制為0.01(CNVTOL,lab,value,toler,norm,minref)。

4)對于線性單元和無中間節點的單元(SOLID65和SOLID45),關閉EXTRA DISPLACEMENTS OPTIONS(在OPTIONS中)。

5)對于CONCRETE材料,可以關閉壓碎功能,將CONCRETE中的單軸抗壓強度設置為-1(tadata,mat,shrcf-op,shrcf-cl,UntensSt,UnCompSt(-1))。

c.非線性計算無法收斂原因

非線性計算無法收斂是非常頭疼的問題,往往不同的分析過程會有不同的提示,總結不收斂的原因,主要可從以下幾大類著手:

1)檢查模型是否存在不完全約束或剛體位移的情況:這主要是通過施加合理的約束方法來誰解決。如果剛體位移不可避免,可以人為施加合理的約束使其收斂,workbench就是通過施加弱彈簧來消除剛體位移的。

2)檢查網格:尤其是錯誤信息提示有“單元出現嚴重扭曲”的語句時,這種情況發生在原模型計算前存在嚴重網格畸變的情形,可以通過 check,esle,warn 選擇存在警告的單元位于何處,找到位置后再在該處細化。

3)檢查載荷子步和時間步設置:一般是增加子步數或者減少時間步長,這點不用多說了,很實用很有效的方法。

4)檢查材料參數設置:這一錯誤通常不容易發現,然而非常重要。材料模型不正確意味著不合理的應力應變關系,在施加載荷后往往出現不合理的結構響應,導致自由度位移過大而不收斂。具體說就是材料的EX,PRXY,金屬材料的雙線性和多線性參數等是否正確。

5)檢查結構是否有不穩定:如果我們分析的結構是一個局部或全局不穩定結果,這一點必須考慮。比如在結構變形過程中出現了屈曲、剛度突變的情況,這個錯誤是非線性不收斂的一個重要錯誤。

6)檢查接觸的設置:接觸是一個狀態非線性問題,ansys中可選的參數很多,對于不同的接觸其中的參數需要個別設置,其中主要的參數有接觸行為方式,法向罰剛度因子,pingball域,初始穿透等等。

7)檢查積分參數的設置:在瞬態分析采用完全求解方法時,默認的方法是Newmark方法,四個積分參數由γ控制,在一些特殊情況下,考慮數值衰減效果更加理想的HHT方法,或者手動設置四個積分參數(但是這必須在你非常清楚積分參數的影響和保證無條件收斂的情況下才可以)。

8)檢查非線性求解器的選擇:ansys默認的求解方法是iterative,即迭代法,即Newton-Raphson迭代,很多情況下該方法是可行的,但有時候直接迭代法direct能更好的收斂。

9)變形分析中約束方程的設置,一旦約束方程的節點的位移不滿足方程存在的條件的話就會產生變形不合理、不收斂的情況。

????,最好的方法就是從out信息中尋找突破點。

d.ANSYS的非線性收斂準則描述:

計算收斂過程圖

CNVTOL, Lab, VALUE, TOLER, NORM, MINREF. U5 b

ANSYS中,非線性收斂準則主要有力的收斂,位移的收斂,彎矩的收斂和轉角的收斂。一般用力的控制加載時,可以使用殘余力的2-范數控制收斂;而位移控制加載時,最好用位移的范數控制收斂。

When SOLCONTROL,ON, TOLER Defaults to 0.005 (0.5%) for force and moment, and 0.05 (5%) for displacement when rotational DOFs are not present.

When SOLCONTROL,OFF, defaults to 0.001 (0.1%) for force and moment.

收斂精度一般可放寬至 5%,以提高收斂速度。

加快收斂的方法有一下幾種

1)可以增大荷載子步數,nsubst,nsbstp,nsbmn,carry

2)修改收斂準則,cnvtol,lab,value,toler,norm,minref

3)打開優化的非線性默認求解設置和某些強化的內部求解算法,solcontrol, key1,key2,key3,vtol(一般情況下,默認是打開的)

4)重新劃分網格,網格的單元不宜太大或太小, 一般在5~10厘米左右

5)檢查模型的正確性

如何看非線性計算時繪出的收斂圖

1.圖片含義:

非線性計算是一個迭代計算的過程,曲線表示兩次迭代之間的誤差,圖中分別表示力和位移在迭代過程中的每次迭代之間的誤差

CRIT表示收斂準則,不同的收斂準則對應不同的顏色的曲線

L2表示收斂殘差,相應的有F L2(力收斂殘差)、M L2(彎矩收斂殘差

2.收斂的判斷依據:

  • L2曲線均在CRIT曲線以下;

  • Time=1(time=1,只是一個完成的標志,這里默認的認為總時間是1,說明到達最后一步了!這個也可以不稱作判據)

  • 滿足以上條件,則非線性分析收斂。

    ansys計算非線性時會繪出收斂圖,該圖是對你計算過程的一個記錄,是對計算過程的直觀認識。其中橫坐標是cumulative iteration number ,是累積迭代次數;縱坐標是absolute convergence norm,絕對收斂范數,由他們來判斷非線性分析是否收斂。

    ansys在每荷載步的迭代中計算非線性的收斂判別準則和計算殘差。其中計算殘差是所有單元內力的范數,只有當殘差小于準則時,非線性疊代才算收斂。ansys默認的收斂準則是將不平衡力的SRSS與VALUEYOLER的值進行比較,對力或力矩進行收斂檢查。也可以添加位移收斂準則。ansys的位移收斂是基于力的收斂的,以力為基礎的收斂提供了收斂量的絕對值,而以位移為基礎的收斂僅提供表現收斂的相對量度。一般不單獨使用位移收斂準則,否則會產生一定偏差,有些情況會造成假收斂。因此ansys官方建議用戶盡量以力(或力矩)為基礎的收斂準則,如果需要也可以增加以位移為基礎的收斂檢查。

    對于多自由度體系的收斂檢查,ANSYS提供了三種不同的矢量范數用于收斂檢查:無限范數、L1范數、L2范數。ANSYS缺省是用L2范數控制收斂。在圖形右邊表示的就是你這次計算使用的收斂準則——crit、L2分別是按照兩種收斂準則計算出來的誤差量。F CRIT-收斂標準;F L2-誤差范數,按照L2收斂準則計算出來的力的誤差量(迭代計算中的概念)。M表示力矩。TIME是與定義的子步與子步數的過程反映。由TIME 和NSUBST命令控制。

    3.繪制收斂圖

    結束后會顯示收斂圖,如果你把它弄消失了,就無法再查看了。只有再算一次。所有小心操作!不過收斂圖不怎么重要,只是模型計算過程的記錄。

    紫的是殘差力,即({F} - {Fnr}),籃的是收斂準則,當殘差在準則以下時,求解收斂

    Ansys 牛頓-拉普森法用如下方程迭代到一個收斂解:

    [KT]{deta u} = {F} - {Fnr}

    以增量形式逐漸施荷加載。

    在每一載荷增量中完成平衡迭代來使得增量求解達到平衡。

    求解平衡方程[KT]{Deta u} = {F} - {Fnr}

    [KT]?=切線剛度矩陣

    {deta u}?=位侈增量

    {F}?=外部載荷向量

    {Fnr}?=內部力向量

    迭代進行,直到{F} - {Fnr}在允許誤差范圍內。程序反復求解此方程(二分載荷增量), 直到殘差(失衡力), {F} - {Fnr},小到可以接受的程度.最大可接受值稱為力收斂準則

    當L2處于crit曲線下時收斂,但是很多時候,L2曲線是在CRIT曲線上下波動,并不都處于CRIT曲線之下

    關于ansys中收斂準則(cnvtol)理解

    ansys中依據缺省的收斂準則,程序將對不平衡力SRSS與VALUE*TOLER的值進行比較;而VALUE的缺省值是在SRSS和MINREF中取較大值。現假如TOLER的缺省值是0.1的話,這個準則是不是可以理解成后一次的SRSS是前一次的SRSS的01倍就收斂啦?

    請指點

    我是這樣理解的例如下面的命令流:

    convergence value 是收斂值,convergence norm是收斂準則。ansys可以用cnvtol命令,如:

    cnvtol,f,5000,0.0005,0

    cnvtol,u,10,0.001,2

    其中f是指采用力結果,5000是收斂絕對值,0.0005是收斂系數,0是收斂范數。如果不平衡力(獨立的檢查每一個自由度)小于等于5000*0.0005(也就是2.5),并且如果位移的變化小于等于10*0.001時,認為子步是收斂的。

    ANSYS中收斂準則,程序默認力與位移共同控制,并且收斂的控制系數好像是0.001。這樣的收斂精度一般很難使塑性分析收斂,對于一般的塑性分析收斂問題,前幾個荷載步(彈性階段)用力與位移共同控制,進入塑性后用力控制或位移控制,也可以先用力后用位移控制(位移控制比較容易收斂),至于控制系數取多少,自己根據需要逐步放大直至收斂!也有人建議最后用能量來控制收斂

    收斂準則應該是指選取那種結果進行收斂判定,通常有三種選擇,分別是力(f),位移(u)、和能量。當然這三種形式可以單獨使用也可以聯合使用。收斂準則的另一層意思應該是選取什么范數形式(1、2、3范數)。一般結構通常都選取2范數格式。

    而收斂值只是收斂準則中的一部分,如cnvtol命令中的收斂絕對值與收斂系數的乘積就應該是你所指的收斂值(convergence value)

    ansys 使用收斂準則有L1,L2,L~~(無窮大)三個收斂準則。

    在工程中,一般使用收斂容差(0.05)就可以。

    建議使用位移收斂準則( cnvtol,u,0.05,,, )與力收斂準則( cnvtol,f,0.05,,, )。因為僅僅只使用一個收斂準則,會存在較大的誤差。

    假如你只能是使用一個收斂準則,建議你提高收斂容差(0.01以下)。

    ANSYS缺省是用L2范數控制收斂。其它還有L1范數和L0范數,可用CNVTOL命令設置。在計算中L2值不斷變化,若L20 E; s, U; n??T3 I' r'

    由于ANSYS缺省的criterion計算是你全部變量的平方和開平方(SRSS)*valuse(你設置的值),所以crition也有小小變化。如有需要,也可自己指定crition為某一常數, CNVTOL,F,10000,0.0001,0?就指定力的收斂控制值為10000*0.0001=1。

    另外,非線性計算中用到的一個開關是SOLCONTROL如關閉SOLCONTROL 選項,那么軟件默認收斂準則:力或彎矩的收斂容差是0.001,而不考慮位移的收斂容差;如果打開SOLCONTROL 選項,同樣的默認收斂準則:力或彎矩的收斂容差是0.005,而位移收斂容差是0.05。

    注:文章來源于網絡

    總結

    以上是生活随笔為你收集整理的怎么作非线性拟合_值得收藏ansys非线性收敛总结的全部內容,希望文章能夠幫你解決所遇到的問題。

    如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。

    国产五月 | 久久tv| 99久久99久国产黄毛片 | av导航福利 | 97成人在线观看视频 | 91成人午夜| 黄色片亚洲| 天天干天天干天天干 | 亚洲激情在线视频 | 午夜精品电影一区二区在线 | 国产精品成人av在线 | 久久99热精品 | 亚洲国产97在线精品一区 | 超碰公开在线观看 | 国产精品女同一区二区三区久久夜 | 在线观看免费福利 | 久草香蕉在线视频 | 99久久99久久精品免费 | 五月天免费网站 | 伊人久久精品久久亚洲一区 | 91精品国自产在线偷拍蜜桃 | 婷色在线 | 久久视频免费在线观看 | 日本精品va在线观看 | 亚洲精品免费视频 | 国产九九热视频 | 日韩综合一区二区 | 亚洲男男gaygay无套同网址 | 亚洲天堂精品 | 欧美日韩调教 | 一区二区三区韩国免费中文网站 | www.天天射.com | 黄色网址av| 久久视频精品在线 | 国语自产偷拍精品视频偷 | 国产又粗又长的视频 | 黄色大全视频 | 九九视频免费观看视频精品 | 国产成人亚洲在线观看 | 激情五月视频 | 久久免费视频在线观看30 | 国产成人综合在线观看 | 国产美女黄网站免费 | 男女啪啪网站 | 激情开心 | 成人av网页| v片在线看 | 韩国av永久免费 | 国产一区福利在线 | 久久999久久 | 欧美日韩高清在线一区 | 中文字幕 国产 一区 | 91.精品高清在线观看 | 91九色视频国产 | 92av视频| 亚洲精品在线一区二区 | 久久久国产在线视频 | 日韩在线视频免费播放 | 九九九毛片 | 久久综合婷婷综合 | 欧美日韩国产mv | 精品色999 | 久久免费看毛片 | 欧美激情视频一区 | 在线播放精品一区二区三区 | 国产永久免费观看 | 天天综合网入口 | 久久久久久97三级 | 中文字幕一区二区三区精华液 | 久草影视在线观看 | 久草免费在线观看视频 | 最新中文字幕在线观看视频 | 四虎成人在线 | 亚洲精品视频免费在线 | 黄色小说视频网站 | 四虎在线免费观看 | 亚洲夜夜爽 | 97超碰免费在线观看 | 久久精品视频在线看 | 久久99影院 | 91精品一区二区在线观看 | 国产一区二区在线看 | 国产精品 日本 | 日韩在线观看你懂的 | 国产精品成人一区二区三区 | 国产日韩欧美在线一区 | 综合色亚洲 | 亚洲一区二区观看 | 中文在线8新资源库 | 亚洲欧洲国产视频 | 91精品久久久久久久久久久久久 | 日韩试看 | 欧美日韩二区三区 | 五月天婷婷综合 | 日韩专区av| 国产精品久久久av久久久 | 日韩二区在线观看 | 999久久久精品视频 日韩高清www | 五月天av在线 | 国产字幕在线观看 | 深夜男人影院 | 黄色亚洲免费 | 久久久国产视频 | 99精品一区 | 天天躁日日 | 91成人精品观看 | 国产资源精品 | 国产91精品一区二区绿帽 | 久久久久久久久综合 | 欧美日本不卡视频 | 成人av一区二区在线观看 | 亚洲国产精彩中文乱码av | 丁香婷婷在线观看 | 亚洲精品视频在线免费播放 | 国产精品视频999 | 免费国产在线精品 | 久久久久久久久久久久国产精品 | www.操.com| 成人午夜精品 | 探花视频在线版播放免费观看 | 亚洲在线视频播放 | 黄色成年片 | 久久久久久久综合色一本 | 免费观看9x视频网站在线观看 | 色亚洲网 | 免费日韩一区二区 | 久久国色夜色精品国产 | 日本乱码在线 | 国产精品精品国产 | 色欧美视频 | 区一区二区三在线观看 | www视频免费在线观看 | 最近2019好看的中文字幕免费 | 久久精彩视频 | 六月天色婷婷 | 奇米影视8888在线观看大全免费 | www.亚洲| 国产精品色婷婷视频 | 91人人澡 | 亚州成人av在线 | 精品国产成人av在线免 | 国产又粗又硬又爽视频 | 亚洲va欧美va人人爽 | 婷婷网在线 | 欧美日韩国产一区二区三区在线观看 | 国产精品原创在线 | 亚州av网站 | 婷婷色六月天 | 日韩欧美在线高清 | 青青色影院 | 操操操天天操 | av女优中文字幕在线观看 | 久久精品男人的天堂 | 日日夜夜精品免费观看 | 深夜免费福利在线 | 又大又硬又黄又爽视频在线观看 | www操操| 亚洲成人av片在线观看 | 国产综合小视频 | 精品99免费视频 | 国产亚州精品视频 | 91精彩视频 | 久久综合狠狠 | 激情综合啪 | 日本中文字幕高清 | 91热这里只有精品 | 国产精品亚州 | 中文字幕在线看视频国产 | 国产精品久久久久永久免费 | 久久久精品在线观看 | 午夜av网站| 看av免费网站 | 成人黄在线观看 | 免费观看完整版无人区 | 精品一区 在线 | 成年人视频免费在线 | 97在线播放视频 | 六月丁香婷 | 在线观看韩日电影免费 | 黄色小说免费在线观看 | 一区二区三区免费在线观看 | 97激情影院 | 国产99久久久国产精品 | 一区二区视频网站 | www.午夜视频| 天堂在线成人 | 亚洲午夜小视频 | 日韩色爱| 热久久国产精品 | av成人资源 | 婷婷深爱| 中文字幕在线观看免费 | 一区二区欧美日韩 | 国产特级毛片aaaaaa | 日韩毛片在线播放 | va视频在线观看 | 亚洲伊人婷婷 | 久久乐九色婷婷综合色狠狠182 | 99精品视频在线看 | 欧美日韩国产综合一区二区 | 国产精品免费视频网站 | 欧美日韩国产一区 | 色播激情五月 | 国产.精品.日韩.另类.中文.在线.播放 | 欧美一级xxxx | 久久久久亚洲天堂 | 青春草免费在线视频 | 欧美 日韩 视频 | 九九九国产 | 日韩素人在线观看 | 国产美女视频 | 91日韩在线| 婷婷在线网站 | 香蕉影视 | 国产午夜精品一区二区三区四区 | 欧美极度另类 | 成片免费观看视频大全 | 色伊人网 | 中文字幕精品一区二区三区电影 | 久久91网| 国产亚洲精品久 | 香蕉色综合 | 久久久久久久久亚洲精品 | 免费欧美高清视频 | 黄色一区二区在线观看 | 激情婷婷网 | 国产精品理论片在线观看 | 亚洲激情网站免费观看 | 狠狠色伊人亚洲综合网站野外 | 国产成人免费网站 | 日韩免费一级a毛片在线播放一级 | 精品免费久久久久 | 五月婷婷丁香 | 国产黄色看片 | 精品国产观看 | 天天操天天弄 | 久久精品亚洲综合专区 | 久久免费成人精品视频 | 久久这里有 | 国产精品日韩精品 | 国产精品在线看 | 国产欧美日韩精品一区二区免费 | 91丨porny丨九色 | 天天爽综合网 | 免费在线黄色av | 亚洲黄色高清 | 免费看的黄网站 | 91看毛片 | 欧美精品久久久久久久久老牛影院 | 香蕉免费 | 天堂v中文 | 国产精品高潮呻吟久久av无 | 天堂网一区二区 | 伊人午夜| 国偷自产视频一区二区久 | 亚洲少妇久久 | 久久精品综合视频 | 日本中文字幕影院 | 久草在线精品观看 | 精品国产电影一区二区 | 久久亚洲区 | 一区电影 | 91在线亚洲 | 亚洲午夜av | 黄色免费观看网址 | 色夜视频 | 九草视频在线观看 | 97视频在线 | 国产精品大片免费观看 | 国产在线高清 | 久久毛片视频 | 亚洲人毛片 | 91久久精品一区二区二区 | www中文在线 | 中文字幕在线观看亚洲 | av片子在线观看 | 人人草网站 | 亚洲激情电影在线 | 国产网站在线免费观看 | 婷婷性综合| www天天干com| 久久久久久久久久免费视频 | 久久精品这里都是精品 | 狠狠色丁香婷婷综合久小说久 | 久久99国产精品久久99 | 亚洲成人资源 | 91中文字幕永久在线 | 97成人免费 | 久久久久久毛片 | 九九热免费在线观看 | 精品一区二区三区香蕉蜜桃 | 蜜臀av麻豆| 国产高清在线观看av | 999国内精品永久免费视频 | 日韩精品久久一区二区三区 | 亚洲国产精品va在线看黑人动漫 | av成人在线看| 中文字幕在线观看的网站 | 久久男女视频 | 国产一区二区日本 | 久久高清毛片 | 最近中文字幕mv免费高清在线 | 婷婷六月综合亚洲 | 在线视频一二三 | 日韩影视在线 | 欧美午夜性生活 | 一级片视频在线 | 日韩欧美精品在线观看视频 | 国产视频一区二区三区在线 | 日本久久成人 | 久99久精品视频免费观看 | 国产欧美在线一区二区三区 | 色国产精品一区在线观看 | 在线日本v二区不卡 | 欧美成人在线网站 | 99久久99视频只有精品 | 91av电影 | 国内精品久久久久影院男同志 | 91精品国产综合久久久久久久 | 欧美国产三区 | 五月情婷婷 | 99视频免费观看 | 成人免费视频网站 | 精品国产诱惑 | 狠狠干我 | 成人在线免费看 | 欧美日韩xx| 国产精品久久久久av | 欧美精品久久久久久久免费 | 久久深爱网 | 青青河边草免费 | 精品视频资源站 | 97成人资源 | 久久免费av电影 | 国产看片 色 | 99精品在这里| 亚洲国产成人久久 | 一级黄色片在线免费观看 | 国产精品自产拍在线观看中文 | 在线高清 | 久久久久久国产精品免费 | 欧美日韩1区 | 婷婷丁香色| 日韩精品视频免费专区在线播放 | 国产高清 不卡 | 欧美视频99 | 成人av在线直播 | 免费在线观看日韩视频 | 91精品国产麻豆 | 日本视频高清 | 69欧美视频| 一区二区三区四区五区六区 | 国产精品久久久久久久久蜜臀 | 婷婷社区五月天 | 成人国产精品av | 久久久久亚洲精品成人网小说 | av天天色 | 中文字幕色播 | 久精品一区 | 91福利视频免费 | 中文字幕亚洲在线观看 | 亚洲国产精品久久久 | 麻豆视频国产在线观看 | 91精品啪 | 精品国产乱码久久久久久久 | 丁香午夜婷婷 | 蜜臀av性久久久久av蜜臀妖精 | 91香蕉国产在线观看软件 | 九九爱免费视频在线观看 | 91精品国自产在线 | 黄色大片入口 | 久久久免费视频播放 | 一二区电影 | 国产色网 | 日韩免费观看一区二区三区 | 91视频免费国产 | 久久久久北条麻妃免费看 | 丁香视频| 国产精品毛片一区视频 | 狠狠狠色 | 五月天婷婷在线观看视频 | 91亚洲在线 | 免费福利在线观看 | 久久字幕精品一区 | 岛国av在线免费 | 久久久久久伊人 | 日本黄色大片免费看 | 日日射天天射 | 91一区啪爱嗯打偷拍欧美 | 91久久久国产精品 | 欧美一级特黄高清视频 | 中文字幕av在线 | 在线观看色网 | 欧美一区二区在线免费看 | 黄色小说视频在线 | 欧美日韩国产网站 | 草在线| 综合色亚洲 | 成人精品电影 | 最新av免费 | 精品视频专区 | 美女视频永久黄网站免费观看国产 | 中文字幕av在线播放 | av免费在线免费观看 | 视频在线一区 | 亚洲一二三区精品 | 日韩欧美在线综合网 | 蜜臀久久99静品久久久久久 | 日韩免费视频在线观看 | 激情综合网五月 | 天天操天天射天天舔 | 99久热在线精品视频观看 | 久久久久国产一区二区三区四区 | 国色天香第二季 | 久久亚洲影视 | 日本精品久久久久中文字幕5 | 国产日产在线观看 | 午夜在线免费观看视频 | 成人av一区二区兰花在线播放 | 日韩电影在线观看一区 | 亚洲国产精品va在线 | 激情五月色播五月 | 99成人精品 | 天天射射天天 | 亚洲综合在线发布 | 在线 高清 中文字幕 | 一级黄色大片在线观看 | 久久综合久久伊人 | 亚洲成a人片在线观看中文 中文字幕在线视频第一页 狠狠色丁香婷婷综合 | 国产欧美三级 | 国产黄网在线 | 天天草天天摸 | 黄色片网站av | 久久人人97超碰精品888 | 探花视频在线观看免费 | 亚洲精品欧美精品 | 国产麻豆精品在线观看 | 中文字幕高清免费日韩视频在线 | 色老板在线 | 一区二区三区在线观看中文字幕 | 天天伊人网 | 免费观看9x视频网站在线观看 | 激情视频二区 | 日韩动态视频 | 在线观看免费av网 | www.国产在线 | 黄色av网站在线免费观看 | 欧美福利片在线观看 | 亚洲自拍偷拍色图 | 久久影院一区 | 久久精品之 | 综合久久五月天 | 亚洲影音先锋 | 视频福利在线观看 | 亚洲视频aaa| 在线视频 91| 免费亚洲电影 | 午夜精品久久久久久久久久久久 | 欧美性猛片, | 99色资源| 99久久久国产精品美女 | 色窝资源 | 999男人的天堂 | 国偷自产视频一区二区久 | 亚洲欧洲av| 国产高清第一页 | 精品成人国产 | 日韩性色| 91天天操 | 91尤物国产尤物福利在线播放 | 狠狠插狠狠操 | 国产精品初高中精品久久 | 国产在线不卡视频 | 天堂网中文在线 | 欧美激情精品久久久久 | 国产精品视频专区 | 久久久麻豆精品一区二区 | 成人影视免费 | a视频在线看 | 91视频传媒| 国产日韩精品在线观看 | 中文字幕在线观看视频一区 | 久在线观看 | 日韩丝袜视频 | 97超碰在线免费观看 | 日本久久精 | 欧美一区二区三区四区夜夜大片 | 久久免费黄色网址 | 日韩电影在线一区二区 | 国产91在线观 | 国产中文在线观看 | 精品视频在线免费观看 | 久久精品国产精品亚洲 | 狠狠综合久久av | 久久久久久久久久久电影 | av大全在线看 | 丁香花中文在线免费观看 | 免费在线精品视频 | 蜜臀av性久久久久蜜臀aⅴ四虎 | 91精选 | 国产精品高清一区二区三区 | 91精品国产91久久久久久三级 | 最近中文字幕大全 | 亚洲精品乱码久久久久久蜜桃动漫 | 视频国产精品 | 91视频免费| 97超碰人人| 人人澡人| 一区二区三区在线观看免费 | www.久久免费 | 在线观看v片 | 天天爱天天操天天射 | 久久久久久久久免费视频 | 欧洲视频一区 | 亚洲va欧洲va国产va不卡 | 日韩影视精品 | 亚洲闷骚少妇在线观看网站 | 日韩免费一二三区 | www久久九 | 日韩伦理片hd | 亚洲国产精品成人综合 | 国产亚洲va综合人人澡精品 | 成人av高清在线 | 玖玖爱在线观看 | 欧美成人猛片 | 国产成人a亚洲精品 | 久久综合精品国产一区二区三区 | 在线播放91| 免费看十八岁美女 | 国产传媒中文字幕 | 2024国产精品视频 | 色婷丁香| 亚洲人xxx| 国产精品嫩草影院99网站 | 日日爽天天操 | 手机av电影在线 | 天天干天天天 | 久久国产亚洲视频 | 免费看成人av | 中文字幕精品在线 | 久久伊人精品天天 | 欧美日韩国产二区三区 | 国产成人在线一区 | a黄色大片| 国产91粉嫩白浆在线观看 | 久久视频中文字幕 | 狠狠操导航 | 国产在线理论片 | 免费看特级毛片 | 夜夜视频资源 | 中文有码在线 | 国产精品美乳一区二区免费 | 99精品视频在线播放观看 | 果冻av在线 | 久久久www成人免费毛片 | 久久免费99精品久久久久久 | 欧美黄污视频 | 在线视频观看成人 | 成人av电影在线播放 | 国产精品18久久久久久首页狼 | 久久久精品国产一区二区三区 | 欧美日韩亚洲一 | 日本激情中文字幕 | 免费a v在线| 国产不卡网站 | 在线观看完整版免费 | 久草电影免费在线观看 | 国产成人一区二区三区免费看 | 一本一本久久a久久精品综合小说 | 欧美va天堂在线电影 | 国产精品99久久久久久宅男 | 国产色资源 | 色视频国产直接看 | 99在线视频免费观看 | 国产手机视频精品 | 国产精品va视频 | 国产精品伦一区二区三区视频 | 狠狠色丁香九九婷婷综合五月 | 香蕉在线视频观看 | 在线观看视频国产一区 | 成人永久在线 | 91av视屏 | 99视频免费看 | 2019精品手机国产品在线 | 国产日韩欧美精品在线观看 | 国产亚洲免费的视频看 | 绯色av一区 | 午夜资源站| 欧美大香线蕉线伊人久久 | 婷婷色5月| 日本黄色免费在线 | 又湿又紧又大又爽a视频国产 | 日日天天av| 国产精品成人国产乱一区 | 激情av在线播放 | 国产中文自拍 | 久久精品一区八戒影视 | 成人免费视频a | 国产第一页福利影院 | 四虎国产永久在线精品 | 欧美va电影| 国产专区在线 | 精品一区免费 | 久久精品99久久 | 久草视频手机在线 | 国产精品免费久久久久久久久久中文 | 国产视频亚洲 | 91视频久久久久 | 日日干日日操 | 日韩一区二区三区高清在线观看 | 国产午夜精品一区二区三区欧美 | a在线观看免费视频 | 精品高清视频 | 久久久久国产精品厨房 | 久久久久久久久久久影院 | 黄色特一级片 | 亚洲精品免费看 | 天堂av免费观看 | 精品中文字幕在线播放 | 中文乱码视频在线观看 | 又粗又长又大又爽又黄少妇毛片 | 亚洲精品成人av在线 | 日韩乱色精品一区二区 | 国产99免费视频 | 亚欧日韩av | 亚洲日本在线一区 | 久久艹人人 | 国产精品久久久久久av | 免费电影播放 | 国产黄网站在线观看 | 久久天堂影院 | 日韩一区二区三区视频在线 | 极品久久久 | 国产精品96久久久久久吹潮 | 91视频链接 | 超碰人人超 | 欧美精品国产综合久久 | 一区二区三区四区久久 | 狠狠干婷婷 | 91九色精品 | 欧美日韩一区三区 | 最近中文字幕第一页 | 一级a毛片高清视频 | 午夜视频在线观看一区二区 | 国产精品99蜜臀久久不卡二区 | 日韩av一卡二卡三卡 | 久久人人爽 | 国产精品久久久久久模特 | 亚洲开心激情 | 精品国产久 | 开心激情久久 | 中文字幕在线观看视频网站 | 一区二区三区四区久久 | 91九色蝌蚪在线 | 国产高清视频免费观看 | 国产视频精选 | a黄色片在线观看 | 特黄一级毛片 | 精品字幕| 九九九电影免费看 | 在线看污网站 | 99热亚洲精品 | 国产成人一区二区啪在线观看 | 日韩在线 一区二区 | 国产精品成人aaaaa网站 | 国产久草在线观看 | 久久欧美在线电影 | 五月天中文在线 | 国产剧情一区二区 | 日本一区二区三区免费看 | 99免费观看视频 | 波多野结衣视频一区二区三区 | 久久免费的精品国产v∧ | 国产91精品在线播放 | 午夜精品久久久久久中宇69 | 国产九九九九九 | 免费av小说 | 午夜影院先 | 人人舔人人射 | 在线久草视频 | 亚洲一区日韩 | 国产精品美女网站 | 在线观看91久久久久久 | 在线视频 国产 日韩 | 日韩一区二区三区免费视频 | 四季av综合网站 | 国产小视频在线看 | av一区二区在线观看中文字幕 | 色综合久久88色综合天天人守婷 | 国产精品第2页 | 亚洲欧美在线综合 | 国产在线久草 | 久久久久久麻豆 | 粉嫩av一区二区三区四区在线观看 | 国产精品一区二区你懂的 | 天天插综合 | 久久精品久久精品 | 天天拍天天操 | av一区二区三区在线 | 亚洲男男gaygayxxxgv | 色欧美88888久久久久久影院 | av夜夜操| 日韩三级中文字幕 | 综合色狠狠 | 美女免费视频网站 | 黄污视频网站 | 国产高清免费在线观看 | 在线视频观看成人 | 亚洲高清视频在线 | 亚洲男男gaygay无套 | 在线观看中文字幕亚洲 | 一区二区三区免费在线播放 | 亚洲一级理论片 | 福利一区二区 | 绯色av一区 | 中文字幕91视频 | 91九色在线观看 | 亚洲国产日韩一区 | 91精品蜜桃| 在线播放一区二区三区 | 日日夜夜天天综合 | 一区二区三区在线观看免费视频 | 亚洲一级电影在线观看 | 久久综合九色综合久久久精品综合 | 亚洲精品网站在线 | 国产美女视频免费观看的网站 | 91激情小视频 | 久久精品中文字幕一区二区三区 | 久久久久久电影 | 亚洲视频在线观看免费 | 中文字幕在线观看一区 | 91高清免费| 欧美一级久久久 | 在线成人中文字幕 | 西西大胆免费视频 | 国产美女视频 | 亚洲永久精品一区 | 美女视频黄色免费 | 日韩在线免费视频 | 欧美一区二区三区在线看 | 日韩大片在线观看 | 亚洲国产精品第一区二区 | 在线91视频 | 国产精品久久久久久久久毛片 | 在线电影a | 国内外成人免费在线视频 | 色com网| 伊人五月天 | 日韩在线观看第一页 | 2019中文最近的2019中文在线 | 国产精品国产亚洲精品看不卡 | 日韩午夜视频在线观看 | 午夜精品久久久久99热app | 色天天中文 | 五月天天色 | av丁香花| 久草亚洲视频 | 福利视频一区二区 | 国产精品无av码在线观看 | 五月婷婷色丁香 | 欧美日韩在线视频观看 | 97人人艹| www.com操| 深夜男人影院 | 日本在线观看一区 | 狠狠色狠狠色综合日日小说 | 久久免费播放视频 | 亚洲成人av片在线观看 | 五月婷婷一级片 | 欧美福利片在线观看 | 狠狠干网址| 九九热免费在线视频 | 人人看人人做人人澡 | 成年人免费观看在线视频 | 日韩精品中文字幕有码 | 99视频免费看 | 亚洲国产高清在线 | 国产一卡二卡在线 | 99热99 | 国产精品一区在线观看 | 日本中文一区二区 | 国产亚洲激情视频在线 | 91成品人影院 | 久久99在线 | 亚洲激情p| 亚洲综合色视频在线观看 | 婷婷丁香九月 | 911免费视频| 精品一区二区在线免费观看 | 国产精品成人av久久 | 黄色aaa毛片 | 成人免费在线观看av | 日韩免费三级 | 亚洲亚洲精品在线观看 | 亚洲电影久久 | a视频在线观看 | 国产婷婷色 | 在线观看成年人 | 天天操狠狠操夜夜操 | 国产又粗又猛又色又黄网站 | 国产一二三区av | 午夜精品一区二区国产 | 色狠狠狠| 欧美性网站 | 欧美色精品天天在线观看视频 | 国产精品入口传媒 | 高潮久久久久久久久 | 国产特级毛片 | 激情久久久久久久久久久久久久久久 | 婷婷在线观看视频 | 亚洲电影久久 | 亚洲 成人 欧美 | 天天插狠狠插 | 亚洲精品国产精品国 | 久草免费在线视频观看 | 日日干av | 黄色小网站免费看 | 国产精品18久久久久久久 | 日韩欧美国产精品 | 国产96在线 | 国内精品久久久久影院日本资源 | 国产色视频一区 | 69夜色精品国产69乱 | 91丨九色丨国产丨porny精品 | 久久国产精品精品国产色婷婷 | 免费一级片久久 | 伊人狠狠干| 久久在线视频精品 | 国产综合精品久久 | 亚洲天堂精品视频 | 91久久精品日日躁夜夜躁国产 | 欧美极品少妇xxxxⅹ欧美极品少妇xxxx亚洲精品 | 免费高清看电视网站 | 99爱在线观看 | 久久国产一区二区三区 | 亚洲精品美女在线观看播放 | av黄色国产 | 91精品国产99久久久久久红楼 | 最近免费观看的电影完整版 | 狠狠色伊人亚洲综合网站野外 | 91亚洲成人| 免费成人短视频 | 国产小视频福利在线 | 成人黄色中文字幕 | 国产这里只有精品 | 久草免费色站 | 成人免费中文字幕 | 特级西西www44高清大胆图片 | 婷婷在线色 | 夜夜摸夜夜爽 | 日韩在线观看av | 欧美在线观看视频一区二区三区 | 高清av免费一区中文字幕 | 免费看片日韩 | 欧美精品一区二区性色 | 人人爱爱人人 | 国产一级一片免费播放放 | 99精品欧美一区二区三区 | 丁香婷婷在线 | 亚洲日本精品视频 | 黄网在线免费观看 | 在线观看免费成人av | 亚洲成av人电影 | 亚洲精品国产综合99久久夜夜嗨 | 免费看特级毛片 | 午夜久久美女 | 成人久久免费 | 久久精品成人热国产成 | 狠狠狠狠干 | 欧美另类性| 欧美一区二区三区免费观看 | 超碰97在线资源 | 91在线91 | 国产亚洲视频在线 | 精品久久久久久久久久岛国gif | 97av影院| 久久久久久国产一区二区三区 | 国产综合福利在线 | 国产伦理久久精品久久久久_ | 99视频免费看 | 啪啪小视频网站 | 国产精品一区二区av日韩在线 | 免费看的毛片 | 欧美性黑人 | 国产成人精品综合 | 国产涩涩在线观看 | 91人人揉日日捏人人看 | 五月天色中色 | 国偷自产视频一区二区久 | 狠狠色噜噜狠狠狠狠2022 | 国产精品一区二区中文字幕 | 免费麻豆视频 | 亚洲永久国产精品 | 婷婷色网| 毛片网站免费 | 一区二区高清在线 | 国产中文字幕网 | 亚洲成人一区 | 在线免费观看视频你懂的 | 国产一二区视频 | 国产在线97 | 国产一区成人在线 | 国产一区二区在线播放 | 国产午夜精品一区二区三区欧美 | 久一网站 | 久久成人国产精品一区二区 | 精品久久久久久久久久久久久 | 中文字幕视频网站 | 国产美女视频网站 | 天天操夜操视频 | 中文字幕在线播放一区 | 久久久午夜电影 | 亚洲高清视频在线观看 | 亚洲免费观看在线视频 | 黄色成人免费电影 | 亚洲一区久久 | 亚洲91精品 | 欧美一二三视频 | 99视频精品免费视频 | 超碰97公开 | 欧美性色网站 | 欧美日韩免费一区二区 | 91免费观看视频在线 | 天天操天天干天天 | 久久99久久99免费视频 | 97超级碰碰碰碰久久久久 | 97成人精品视频在线播放 | 久久免费视频3 | 国产在线观看国语版免费 | 久草在线91 | 精品视频在线播放 | 91精品国产自产老师啪 | 国产精品久久久久久久久久白浆 | 久久综合五月婷婷 | 日韩xxxbbb | 成人免费电影 | 成人毛片网 | 国内丰满少妇猛烈精品播放 | 国产精品片 | 亚洲国产日韩在线 | 国产理伦在线 | www色,com| 超碰97在线资源 | 久久久久国产精品厨房 | 黄色日批网站 | 亚洲欧美在线观看视频 | 国产黄色片免费看 | 色综合天天射 | 久久久一本精品99久久精品66 | 四虎永久免费网站 | 97超碰中文字幕 | 日本黄色免费网站 | 在线a人片免费观看视频 | 国产美女免费视频 | 亚洲一二区视频 | 91人人澡人人爽人人精品 | 黄色毛片视频免费观看中文 | 开心激情五月网 | 国产高清专区 | 国产精品免费人成网站 | 国产 日韩 在线 亚洲 字幕 中文 | 成人在线播放网站 | 四虎影视成人 | 五月婷婷毛片 | a级片久久久 | 波多野结衣在线观看一区二区三区 | 国产一区二区在线观看免费 | 欧美a免费| 激情深爱 | 色视频在线免费 | 爱情影院aqdy鲁丝片二区 | 久久不卡免费视频 | 一区二区三高清 | 青青射| 成人黄色小说在线观看 | 国产视频高清 | 五月开心婷婷 | www天天干| 国产亚洲精品久久网站 | 日本精品中文字幕在线观看 | 在线国产一区二区三区 | av观看在线观看 | 五月激情婷婷丁香 | 久久久久色 | 国产精品午夜免费福利视频 | 国产一区二区三精品久久久无广告 | 欧美日韩精品在线观看 | 日本三级久久 | 美女一级毛片视频 | 丁香花在线视频观看免费 | 成人免费亚洲 | 日韩av在线免费看 | 欧美黑人性猛交 | 国产一级免费在线观看 | 亚洲国产成人在线播放 | 日本性生活免费看 | 五月婷婷天堂 | 亚洲午夜av电影 |