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

    壩基面受力狀態(tài)的有限元計(jì)算誤差影響分析

    2017-02-27 00:30:55程帥李守義司政杜占科李萌郝曉飛
    南水北調(diào)與水利科技 2017年1期
    關(guān)鍵詞:誤差分析有限元法

    程帥+李守義+司政+杜占科+李萌+郝曉飛

    摘要:有限元數(shù)值計(jì)算在水利工程中運(yùn)用廣泛,其計(jì)算精度和計(jì)算誤差是需要考慮和控制的重要問(wèn)題,對(duì)其進(jìn)行研究可為有限元法在水利設(shè)計(jì)和復(fù)核中提供技術(shù)支撐。以重力壩為研究載體,根據(jù)壩體荷載合力值與有限元計(jì)算的壩基面內(nèi)力值(壩基面應(yīng)力積分值)的平衡條件,分析有限元應(yīng)力計(jì)算誤差的影響因素,提出減小誤差的措施。研究結(jié)果表明:在同一受力狀態(tài)下,提取應(yīng)力的對(duì)象不同計(jì)算誤差不同,需針對(duì)所研究的問(wèn)題及誤差控制要求合理選擇相應(yīng)提取應(yīng)力的對(duì)象;將緊接壩體的薄層地基保持與壩體相同的模型尺寸和約束條件,可充分減小壩基面有限元應(yīng)力計(jì)算誤差;合理選擇模型網(wǎng)格尺寸、地基模擬范圍、應(yīng)力提取路徑的條數(shù)、提取點(diǎn)間距,可在不消耗過(guò)量的計(jì)算資源和時(shí)間下滿足精度要求。研究成果可為控制類似結(jié)構(gòu)的有限元計(jì)算誤差提供參考。

    關(guān)鍵詞:水工結(jié)構(gòu)工程;有限元法;誤差分析;壩基面;受力狀態(tài)

    中圖分類號(hào):TV314文獻(xiàn)標(biāo)志碼:A文章編號(hào):1672-1683(2017)01-0179-07

    Abstract:Finite element numerical calculation is widely used in hydraulic engineering.Its accuracy and error is one of the important problems that need to be considered and controlled.Studies on this issue can provide technical support for the FEM in the design and check of hydraulic engineering.Gravity dam was used as the object of research.The factors influencing the error in finite element stress calculation were analyzed by comparing the resultant load value of the gravity dam with the internal force value of the dam base surface from finite element calculation (integral value of stress of dam base surface).The results showed that extracting the stress from different objects results in different finite element calculation error,even under the same stress state.It is necessary to select the proper object to extract the stress from according to the research goal and error control requirements.The error can be greatly reduced by making the thin-layer foundation adjacent to the dam the same size and in the same constraint condition as the dam body.Proper determination on mesh size,simulation range of foundation,number of paths to extract the stress,and spacing of extracting points can meet the required precision without overly consuming resources or time.The results can provide reference for controlling the finite element calculation error in similar projects.

    Key words:hydraulic structure engineering;finite element method;error analysis;dam base surface;forced state

    隨著計(jì)算機(jī)技術(shù)及數(shù)值分析方法的迅速發(fā)展,有限元數(shù)值計(jì)算在各個(gè)領(lǐng)域得到了越來(lái)越廣泛的應(yīng)用[1-2]。近幾年,有限元法在水利工程中正發(fā)揮著重要作用,特別是大中型水利工程中有限元的應(yīng)用已經(jīng)相當(dāng)廣泛[3-5]。目前,對(duì)于復(fù)雜水利工程中的有限元計(jì)算結(jié)果,在沒(méi)有試驗(yàn)數(shù)據(jù)和標(biāo)準(zhǔn)答案的情況下,如何判斷其正確性,如何用所掌握的信息對(duì)設(shè)計(jì)進(jìn)行完善和修正,如何控制其計(jì)算精度和計(jì)算誤差,這些問(wèn)題沒(méi)有得到很好的解決。這也是為何至今利用有限元軟件進(jìn)行有限元分析的水工建筑物或水力模型僅作為水利工程設(shè)計(jì)參考和復(fù)核,并未給出設(shè)計(jì)或檢驗(yàn)標(biāo)準(zhǔn)的原因之一。

    影響有限元計(jì)算精度的因素很多,Jin[4],Mirzabozorg[5],Kennedy[6],Natarajan[7]等都針對(duì)提高有限元計(jì)算精度開(kāi)展了很多研究。國(guó)內(nèi)已有學(xué)者對(duì)其中的部分因素進(jìn)行了探討。虞皓影[8]、張國(guó)新[9]等在利用有限元分析結(jié)構(gòu)穩(wěn)定性中指出邊界計(jì)算范圍、應(yīng)力積分方法等因素對(duì)計(jì)算結(jié)果影響較為敏感,并各自提出了減小有限元計(jì)算誤差的方法。鄒超英等[10]通過(guò)對(duì)拱壩進(jìn)行應(yīng)力分析,研究了有限元單元類型、自由度數(shù)量、單元階次等因素與能量誤差的關(guān)系。網(wǎng)格密度也是一個(gè)影響計(jì)算精度的因素,李濤等[11]、章春亮等[12]等指出,有限元求解模型中所有關(guān)鍵性區(qū)域的網(wǎng)格劃分水平?jīng)Q定整個(gè)有限元模型的計(jì)算結(jié)果精度,理論上單元數(shù)量的增多,結(jié)果的精確度就會(huì)提高,并收斂于真實(shí)解。單元類型及網(wǎng)格劃分方法對(duì)有限元計(jì)算精度也有一定影響,姬賀炯等[13]、Cavin等[14]、張麗媛等[15]分別對(duì)這些因素進(jìn)行了研究,為相應(yīng)結(jié)構(gòu)的有限元計(jì)算誤差提供了理論依據(jù)。然而,對(duì)于大體積水工結(jié)構(gòu)的模擬,網(wǎng)格劃分過(guò)密,會(huì)消耗大量的資源和時(shí)間,需確定相應(yīng)網(wǎng)格尺寸來(lái)兼顧計(jì)算效率和精度。此外,運(yùn)用有限元分析結(jié)構(gòu)受力狀態(tài)、穩(wěn)定性等問(wèn)題時(shí),通常需要選擇整體結(jié)構(gòu)中的不同對(duì)象來(lái)提取計(jì)算截面的應(yīng)力,選擇不同的提取應(yīng)力對(duì)象,同一受力位置的有限元計(jì)算誤差不同,而在既定的應(yīng)力提取對(duì)象下應(yīng)力路徑的范圍及疏密程度也會(huì)影響其計(jì)算精度。目前,在大體積水工結(jié)構(gòu)的有限元計(jì)算中,針對(duì)上述幾點(diǎn)影響因素及誤差影響分析方面的研究甚少。

    壩基面作為壩體的薄弱控制面,該部位的受力狀態(tài)對(duì)分析壩體穩(wěn)定性等問(wèn)題有很大影響。本文以重力壩為研究載體,根據(jù)力學(xué)平衡條件得出壩基面有限元計(jì)算結(jié)果的誤差,對(duì)同一位置不同提取應(yīng)力對(duì)象下有限元計(jì)算值的偏差及其原因進(jìn)行了分析,提出了減小誤差的措施,考察了網(wǎng)格尺寸大小、地基模擬范圍、應(yīng)力路徑條數(shù)、應(yīng)力點(diǎn)間距等因素對(duì)壩基面有限元計(jì)算誤差的影響。研究結(jié)果可為分析導(dǎo)墻、閘墩、廠房、擋土墻等類似工程結(jié)構(gòu)的有限元計(jì)算誤差提供依據(jù)。

    1 研究資料

    1.1 研究對(duì)象概況

    本文研究主要依托某水利樞紐工程中一重力壩壩段開(kāi)展。該工程屬Ⅱ等大(2)型工程,攔河壩為混凝土重力壩,最大壩高75.97 m,壩頂寬度7.0 m,上游壩面鉛直,下游壩坡1∶0.8,最大壩底寬度60.74 m。正常蓄水位450 m,死水位440 m,下游正常尾水位405 m,壩前淤沙高程389.3 m。

    1.2 荷載施加

    本文在進(jìn)行有限元計(jì)算時(shí),考慮的荷載有:自重、上下游靜水壓力、揚(yáng)壓力、浪壓力、泥沙壓力。

    1.3 材料參數(shù)

    不同材料的計(jì)算參數(shù)取值見(jiàn)表1。

    2.1 有限元模型

    本文采用ANSYS有限元軟件對(duì)某一重力壩壩段進(jìn)行分析,壩段寬18 m,整體有限元計(jì)算模型包括壩體與基礎(chǔ)巖體?;鶐r的模擬范圍取矩形區(qū)域,基巖邊界值:深度方向均取152 m,約為2倍壩高,左、右岸方向與壩段同寬,上下游方向取2倍壩高(在研究地基范圍影響時(shí)分類考慮)。模型整體坐標(biāo)系的原點(diǎn)設(shè)在壩踵右側(cè)底部,沿水流方向指向下游為X軸正方向;沿高度方向鉛直向上為Y軸正方向;垂直水流方向指向右岸為Z軸正方向。整體模型網(wǎng)格剖分基本采用8結(jié)點(diǎn)六面體實(shí)體單元,壩基緊接壩體10 m范圍內(nèi)網(wǎng)格尺寸與壩體相同,壩基其它區(qū)域網(wǎng)格尺寸較大,中間10 m厚度的槽型區(qū)域通過(guò)四面體單元進(jìn)行過(guò)渡。整體有限元計(jì)算模型見(jiàn)圖1。

    2.2 計(jì)算假定

    根據(jù)重力壩、廠房等類似結(jié)構(gòu)的特點(diǎn)及其相關(guān)研究成果[16-18],本文在三維有限元靜力計(jì)算中,做如下假定。

    (1)壩段間設(shè)置有橫縫,因此各方案下壩段獨(dú)立承受荷載,壩段間無(wú)相互作用。

    (2)充分考慮壩體的材料分區(qū),在靜力計(jì)算中假定壩體混凝土是線彈性的,地基為均勻無(wú)質(zhì)量彈性體,二者皆滿足彈性力學(xué)基本假定。

    (3)由于橫縫作用,在計(jì)算模型中壩體兩側(cè)無(wú)約束,地基巖體上下游面及左右側(cè)面考慮法向約束,基巖底面考慮三向約束(在研究邊界約束影響時(shí)分類考慮)。

    2.3 研究方法

    分析截面受力狀態(tài),需對(duì)應(yīng)力進(jìn)行曲面積分以得到建基面截面內(nèi)力[19-21]。由以上計(jì)算模型和資料,施加荷載,得到有限元應(yīng)力計(jì)算結(jié)果。選取上部壩體結(jié)構(gòu)或下部地基結(jié)構(gòu)作為應(yīng)力提取對(duì)象,沿著壩軸線方向取N條路徑,并選定每條路徑的提取點(diǎn)間距,提取建基面上的法向應(yīng)力和切向應(yīng)力,按式(1)沿建基面進(jìn)行積分后得到法向和切向力的合力。

    式中:σy,τxy為提取的各點(diǎn)法向和切向應(yīng)力;Fy,F(xiàn)s為法向和切向力的合力值。

    根據(jù)力學(xué)平衡原理計(jì)算出外荷載合力值與有限元計(jì)算得出的內(nèi)力值作對(duì)比,按式(2)分析其相對(duì)誤差。

    式中:F0為外荷載合力值;Fi為有限元計(jì)算的內(nèi)力值。

    采取一系列方案得到有限元計(jì)算結(jié)果,根據(jù)各方案相對(duì)誤差分析壩基面受力狀態(tài)的有限元計(jì)算誤差影響因素,研究減小誤差的措施。

    3 研究?jī)?nèi)容與結(jié)果分析

    3.1 提取應(yīng)力對(duì)象的影響

    分別選取上部壩體結(jié)構(gòu)和下部地基結(jié)構(gòu)為應(yīng)力提取對(duì)象,提取壩基面上的法向和切向應(yīng)力,積分得合力值與外荷載合力值對(duì)比,計(jì)算其相對(duì)誤差見(jiàn)表2。

    由表2可知,有限元計(jì)算結(jié)果中,提取應(yīng)力的對(duì)象不同,法向力和切向力都有所變化,且誤差不同,法向力誤差對(duì)于應(yīng)力提取對(duì)象的敏感度小于切向力誤差。選擇下部地基結(jié)構(gòu)作為應(yīng)力提取對(duì)象計(jì)算誤差偏大,法向力誤差超過(guò)5%,切向力誤差甚至達(dá)到12.3%。在實(shí)際分析中,需合理選擇提取應(yīng)力的對(duì)象,才能保證計(jì)算結(jié)果的精確性。由于壩基面部位對(duì)分析壩體穩(wěn)定性等問(wèn)題極為關(guān)鍵,且某些問(wèn)題中必須以下部地基結(jié)構(gòu)作為提取應(yīng)力的對(duì)象,以下分析選擇下部地基作為提力對(duì)象時(shí)誤差過(guò)大的原因。

    由地基與壩體交接面的特點(diǎn)可知,這種上下偏差較大的現(xiàn)象可能是由于材料屬性變化、模型尺寸突變、邊界約束條件的不同等因素造成的?,F(xiàn)針對(duì)以上幾種可能的影響因素做以下六種方案進(jìn)行研究:方案①:原始計(jì)算模型,材料、尺寸、邊界約束都不變;方案②:將地基的材料屬性變?yōu)榕c壩體相一致,其它條件不變;方案③:將地基范圍的X方向長(zhǎng)度取為與壩體相等,其它條件不變;方案④:地基范圍不變,在緊接壩體部位一個(gè)單元厚度處地基與壩體X方向長(zhǎng)度相等,且約束不變;方案⑤:在地基的原約束基礎(chǔ)上,將建基面向下一個(gè)單元厚度的地基約束去掉,使其與壩體的約束條件相一致,地基其他邊界約束不變;方案⑥:地基范圍不變,在緊接壩體部位一個(gè)單元厚度處地基與壩體X方向長(zhǎng)度相等,且無(wú)約束,其它條件不變。各方案的模型、材料及邊界約束示意圖見(jiàn)圖2。

    根據(jù)圖2中各種方案,建立相應(yīng)的有限元計(jì)算模型,施加荷載,分別以上部壩體結(jié)構(gòu)和下部地基結(jié)構(gòu)作為提取應(yīng)力的對(duì)象,壩基面上的法向和切向力的相對(duì)誤差計(jì)算結(jié)果見(jiàn)圖3。

    由圖3可知,選擇下部地基結(jié)構(gòu)誤差較大,其中法向力誤差在2.8%~6.5%之間,切向力誤差在0.3%~27.9%之間;選取上部壩體結(jié)構(gòu)計(jì)算相對(duì)較精確,其中法向力誤差在1.7%~3.5%之間,切向力誤差在0.7%~1.5%之間。切向力誤差變幅較法向力更大,其中下部切向力誤差對(duì)方案的選取最為敏感,變幅最大。

    由方案①和②可知,雖然將壩體與地基的材料屬性取為一致,但其上下切向應(yīng)力值相差依然很大,下部切向力誤差都在10%以上,故造成這種誤差的原因并非材料屬性的變化所致。由方案③和④可知,兩種不同地基范圍下,在壩基面處模型尺寸均無(wú)突變,但上下切向應(yīng)力相差依然很大,下部切向力誤差達(dá)到28%,可知造成這種誤差的原因并非模型尺寸突變所致。方案⑤和⑥驗(yàn)證是否為邊界約束條件的影響,方案⑤中解除了建基面向下一個(gè)單元厚度的地基約束,誤差值并未減小,下部切向應(yīng)力依然有12%,各項(xiàng)誤差與方案①相近。方案⑥中緊接壩體部位一個(gè)單元厚度處地基與壩體X方向長(zhǎng)度相等,且無(wú)約束,使此處一小部分地基的約束條件和模型尺寸皆與壩體相一致,這種情況下,地基面上、下部的法向力和切向力誤差都達(dá)到最小,上部和下部提取的力幾乎相等,法向誤差減小至1.7%,切向誤差減小至0.3%。

    上述六種方案充分說(shuō)明:造成壩基上部與下部的所提取力的較大偏差的原因,不是材料屬性變化所致,而是模型尺寸突變和邊界約束條件變化的雙重原因所致。故為充分減小有限元計(jì)算誤差,建議將緊接壩體很小厚度范圍(一個(gè)單元尺寸)內(nèi)的薄層地基保持與壩體同寬且不做任何約束,即保持與壩體具有相同的模型尺寸和約束條件。

    3.2 網(wǎng)格尺寸的影響

    為研究壩體網(wǎng)格尺寸與地基網(wǎng)格尺寸對(duì)計(jì)算精度的影響,現(xiàn)設(shè)置A、B兩組方案見(jiàn)表3。其中A組方案中地基網(wǎng)格尺寸固定為10 m,研究壩體網(wǎng)格尺寸由10 m減小至1 m時(shí)的計(jì)算誤差;B組方案中壩體網(wǎng)格尺寸固定為1 m,研究地基網(wǎng)格尺寸由10 m減小至1 m時(shí)的計(jì)算誤差。計(jì)算結(jié)果見(jiàn)圖4。

    由圖4可知,網(wǎng)格尺寸的大小會(huì)影響計(jì)算精度。壩體網(wǎng)格尺寸的變化對(duì)計(jì)算誤差的影響更為敏感,結(jié)果與文獻(xiàn)[11]中關(guān)鍵區(qū)域的網(wǎng)格劃分水平?jīng)Q定整個(gè)有限元計(jì)算結(jié)果精度的觀點(diǎn)相吻合。當(dāng)壩體與地基的網(wǎng)格尺寸同為10 m時(shí),法向力誤差達(dá)到7.45%,切向力誤差達(dá)到2.61%,隨著壩體網(wǎng)格尺寸的減小,計(jì)算誤差也隨之減小,計(jì)算精度越精確。當(dāng)網(wǎng)格尺寸減小到一定程度以后,對(duì)計(jì)算結(jié)果準(zhǔn)確性的提高就很小,且會(huì)導(dǎo)致計(jì)算耗時(shí)大幅度增加。當(dāng)壩體網(wǎng)格尺寸為2 m時(shí),法向誤差降至2.41%,切向誤差降至0.63%。地基的網(wǎng)格尺寸對(duì)計(jì)算結(jié)果影響甚微,當(dāng)壩體網(wǎng)格尺寸固定為1 m時(shí),地基網(wǎng)格尺寸由10 m減小至1 m,計(jì)算誤差都維持在一個(gè)相對(duì)穩(wěn)定的水平。綜合考慮計(jì)算機(jī)耗時(shí)和計(jì)算精度的影響,在水利工程中大體積結(jié)構(gòu)模擬中,地基網(wǎng)格尺寸取為10 m,結(jié)構(gòu)網(wǎng)格尺寸取為2 m,即可滿足分析要求。

    3.3 地基范圍的影響

    為分析地基模擬范圍對(duì)計(jì)算精度的影響,分C、D兩組方案進(jìn)行研究,各方案取值如表4所示(表中H為壩高)。其中C組方案中地基深度方向(Y向)保持2倍壩高不變,研究上下游方向(X向)地基范圍由壩踵和壩趾分別向上、下游延伸不同長(zhǎng)度時(shí)壩基面合力的計(jì)算誤差;D組方案中X向地基范圍分別向上、下游延伸2倍壩高不變,研究Y向地基不同深度范圍下的壩基面合力值的計(jì)算誤差。計(jì)算結(jié)果見(jiàn)圖5。

    由圖5可知:地基范圍的不同會(huì)導(dǎo)致計(jì)算精度的不同。X方向地基范圍的變化對(duì)計(jì)算誤差的影響更為敏感。Y方向地基深度保持2倍壩高不變,隨著X方向地基范圍的增大,計(jì)算誤差整體上呈減小趨勢(shì),法向誤差在X方向地基范圍增大時(shí),也會(huì)出現(xiàn)略微的增大現(xiàn)象,當(dāng)X方向地基范圍由壩踵和壩趾分別向上游和下游延伸2倍壩高及以上時(shí),法向誤差會(huì)控制在2.8%以下,切向誤差會(huì)控制在0.85%以下。地基Y方向的范圍對(duì)計(jì)算結(jié)果的誤差影響較小,當(dāng)X方向地基范圍由壩踵和壩趾分別向上游和下游延伸2倍壩高不變時(shí),地基深度范圍由0.2倍壩高增加至1倍壩高時(shí),隨深度的增加誤差略有減小,當(dāng)深度大于1倍壩高時(shí),計(jì)算誤差雖有所偏差,但浮動(dòng)較小,法向誤差在2.8%左右,切向誤差在0.85%左右。

    考慮到地基范圍過(guò)大會(huì)造成模型單元與節(jié)點(diǎn)數(shù)目的成倍增加,需要更多的計(jì)算資源和時(shí)間。故建議在類似工程的模擬中,X方向地基范圍由壩踵和壩趾分別向上游和下游延伸2倍壩高,Y方向地基深度取為2倍壩高。

    3.4 應(yīng)力路徑條數(shù)及應(yīng)力點(diǎn)間距的影響

    以上是針對(duì)有限元應(yīng)力計(jì)算的誤差影響分析。在確定的應(yīng)力計(jì)算結(jié)果下,壩基面的合力計(jì)算值也會(huì)有不同的精度,這是由于計(jì)算合力時(shí)應(yīng)力積分路徑條數(shù)與應(yīng)力點(diǎn)間距的選取不同而造成的,現(xiàn)通過(guò)E、F兩組方案分析這兩個(gè)因素對(duì)壩基面合力值計(jì)算的誤差影響,各方案設(shè)置見(jiàn)表5(表中B為壩段寬度,L為路徑總長(zhǎng)度)。其中,E組方案在相同的應(yīng)力點(diǎn)間距下(0.01L),沿著壩軸線方向取不同的應(yīng)力路徑條數(shù),分析應(yīng)力路徑條數(shù)對(duì)計(jì)算精度的影響;F組方案固定路徑條數(shù)(B/2),分析不同應(yīng)力點(diǎn)間距下計(jì)算誤差的影響。各方案計(jì)算結(jié)果見(jiàn)圖6。

    由圖6可知,應(yīng)力提取路徑的條數(shù),以及每條路徑提取點(diǎn)的間距,對(duì)計(jì)算精度有一定的影響。同一提取點(diǎn)間距,隨著路徑條數(shù)的增多,計(jì)算誤差會(huì)減小,當(dāng)條數(shù)大于B/2時(shí),法向力誤差可控制在2.6%以內(nèi),切向力誤差可控制在1.0%以內(nèi),且繼續(xù)增多路徑精度的提高效果不明顯。在同一提取路徑條數(shù)(B/2)下,當(dāng)提取應(yīng)力點(diǎn)間距過(guò)大時(shí),為負(fù)偏差(標(biāo)準(zhǔn)值小于計(jì)算值),間距過(guò)小時(shí),為正偏差(標(biāo)準(zhǔn)值大于計(jì)算值)。隨著應(yīng)力提取點(diǎn)間距的減小,計(jì)算誤差會(huì)由負(fù)到正逐漸增大,最后會(huì)處于一個(gè)穩(wěn)定狀態(tài),當(dāng)間距小于0.01 L時(shí),法向力誤差可控制在2.0%以內(nèi),切向力誤差可控制在0.9%以內(nèi)。

    由以上結(jié)果可知,應(yīng)力路徑條數(shù)取B/2即可滿足分析要求。綜合考慮計(jì)算精度和計(jì)算工作量,提取應(yīng)力點(diǎn)間距取為0.006 25L~0.02L為宜。

    4 結(jié)論

    (1) 在不同提取應(yīng)力的對(duì)象下,壩基面受力的有限元計(jì)算值不同,且計(jì)算誤差不同。選擇上部壩體結(jié)構(gòu)作為研究對(duì)象計(jì)算誤差較小,地基面上的法向應(yīng)力誤差對(duì)于提力對(duì)象的敏感度小于切向應(yīng)力誤差。在類似工程中,需針對(duì)所研究的結(jié)構(gòu)及誤差控制要求合理選擇相應(yīng)提取應(yīng)力的對(duì)象。

    (2) 在有限元計(jì)算模型中,將緊接壩體很小厚度范圍(一個(gè)單元尺寸)內(nèi)的薄層地基保持與壩體相同的模型尺寸和約束條件,可充分減小有限元計(jì)算誤差。其中法向誤差可控制在2.0%以內(nèi),切向誤差可控制在0.5%以內(nèi)。

    (3)對(duì)類似大體積水工結(jié)構(gòu)進(jìn)行模擬并分析其受力狀態(tài)時(shí),綜合考慮計(jì)算精度、計(jì)算資源和時(shí)間、工作效率等因素,結(jié)構(gòu)網(wǎng)格尺寸取為2 m,地基網(wǎng)格尺寸取為10 m,地基范圍在上下游方向向兩側(cè)延伸2倍建筑物高度,深度方向取為2倍建筑物高,即可在不消耗過(guò)量的計(jì)算資源和時(shí)間下滿足精度要求;應(yīng)力路徑條數(shù)取B/2即可滿足分析要求,應(yīng)力點(diǎn)間距取為0.006 25L~0.02L為宜。

    考慮到影響有限元計(jì)算精度的因素很多,如模型本身的結(jié)構(gòu)誤差、階次誤差,荷載分布的簡(jiǎn)化,材料本構(gòu)模型選取等,各個(gè)因素是相互關(guān)聯(lián),共同作用的。此外,亦可通過(guò)其他方法評(píng)價(jià)有限元計(jì)算結(jié)果精度,如對(duì)比試驗(yàn)數(shù)據(jù)、觀察應(yīng)力等值線圖、核算關(guān)鍵節(jié)點(diǎn)位移等。這些規(guī)律還有待于進(jìn)一步研究。本文以經(jīng)典彈性力學(xué)為基礎(chǔ),視混凝土和地基材料為均質(zhì)連續(xù)體,通過(guò)對(duì)比有限元應(yīng)力計(jì)算值與外荷載合力值,論述了部分影響壩基面受力狀態(tài)計(jì)算誤差的因素,所得結(jié)論適用于水利工程中的大體積混凝土結(jié)構(gòu),可滿足實(shí)際工程的分析需求,對(duì)類似工程具有一定的指導(dǎo)意義。

    參考文獻(xiàn)(References):

    [1] 曾攀.工程有限元方法[M].北京:科學(xué)出版社,2010.(ZENG Pan.The finite element method in engineering[M].Beijing:Science Press,2010.(in Chinese))

    [2] 何本國(guó).ANSYS土木工程應(yīng)用實(shí)例[M].北京:中國(guó)水利水電出版社,2011.(HE Ben-guo.Application examples of ANSYS civil engineering[M].Beijing:China Water Power Press,2011.(in Chinese))

    [3] 王樹(shù)仿,李向明.三維有限元靜力分析在小型引水式電站廠房中的應(yīng)用[J].浙江水利科技,2014,42(4):49-53.(WANG Shu-fang,LI Xiang-ming.A small diversion-type power station plant three-dimensional finite element static analysis[J].Zhejiang Hydrotechnics,2014,42(4):49-53.(in Chinese))DOI:10.13641/j.cnki.33-1162/tv.2014.04.020

    [4] Jin B,Lazarov R,Pasciak J,et al.Error analysis of a finite element method for the space-fractional parabolic equation[J].SIAM Journal on Numerical Analysis,2014,52(5):2272-2294.DOI:10.1137/13093933X

    [5] Mirzabozorg H,Hariri-Ardebili M A,Heshmati M,et al.Structural safety evaluation of Karun III Dam and calibration of its finite element model using instrumentation and site observation[J].Case Studies in Structural Engineering,2014,1:6-12.DOI:10.1016/j.csse.2014.02.001

    [6] Kennedy F E,Colin F,F(xiàn)loquet A,et al.Improved techniques for finite element analysis of sliding surface temperatures[J].Developments in Numerical and Experimental Methods Applied to Tribology,2014:138-150.DOI:10.1016/B978-0-408-22164-1.50025-9

    [7] Natarajan S,Ooi E T,Chiong I,et al.Convergence and accuracy of displacement based finite element formulations over arbitrary polygons:Laplace interpolants,strain smoothing and scaled boundary polygon formulation[J].Finite Elements in Analysis and Design,2014,85:101-122.DOI:10.1016/j.finel.2014.03.006

    [8] 虞皓影,何建平.影響有限元法分析滑坡穩(wěn)定性的因素研究[J].河南工程學(xué)院學(xué)報(bào):自然科學(xué)版,2009,21(1):33-35.(YU Hao-ying,HE Jian-ping.Research on factors of element method for slope stability analysis[J].Journal of Henan Institute of Engineering,2009,21(1):33-35.(in Chinese))

    [9] 張國(guó)新,劉毅.壩基穩(wěn)定分析的有限元直接反力法[J].水力發(fā)電,2006,32(12):30-32.(GHANG Guo-xin,LIU Yi.Finite element direct force method in stability analysis of dam foundation [J].Water Power,2006,32(12):30-32.(in Chinese))

    [10] 鄒超英,蘇志敏,曾海軍,等.拱壩有限元分析網(wǎng)格剖分方案研究[J].南水北調(diào)與水利科技,2011,9(1):47-49.(ZOU Chao-ying,SU Zhi-min,ZENG Hai-jun,et al.Mesh generation of finite element method for arch dam[J].South-to-North Water Diversion and Water Science&Technology,2011,9(1):47-49.(in Chinese))DOI:10.3724/SP.J.1201.2001.01047

    [11] 李濤,左正興,廖日東.結(jié)構(gòu)仿真高精度有限元網(wǎng)格劃分方法[J].機(jī)械工程學(xué)報(bào),2009,45(6):304-308.(LI Tao,ZUO Zheng-xing,LIAO Ri-dong.Meshing method of high precision FEM in structural simulations[J].Chinese Journal of Mechanical Engineering,2009,45(6):304-308.(in Chinese))DOI:10.3901/JME.2009.06.304

    [12] 章春亮.有限元分析中的單元?jiǎng)澐仲|(zhì)量與計(jì)算精度[J].輕工機(jī)械,2002 (1):27-30.(ZHANG Chun-liang.Meshing quality and calculation precision of finite element analysis[J].Light Industry Machinery,2002 (1):27-30.(in Chinese))

    [13] 姬賀炯,白長(zhǎng)青,韓省亮.輸流管道動(dòng)力有限元建模及實(shí)驗(yàn)研究[J].應(yīng)用力學(xué)學(xué)報(bào),2013,30(3):422-428.(JI He-jiong,BAI Chang-qing,HAN Sheng-liang.Dynamic finite element modeling and experimental research of the fluid-filled pipeline[J].Chinese Journal of Applied Mechanics,2013,30(3):422-428.(in Chinese))DOI:10.11776/cjam.30.03.C059

    [14] Cavin P,Gravouil A,Lubrecht A A,et al.Efficient FEM calculation with predefined precision through automatic grid refinement[J].Finite Elements in Analysis and Design,2005,41(11):1043-1055.DOI:10.1016/j.finel.2004.09.004

    [15] 張麗媛,胡昱,李慶斌,等.基于平面布爾運(yùn)算的重力壩3D有限元建模[J].水力發(fā)電學(xué)報(bào),2012,31(5):209-215.(ZHANG Li-yuan,HU Yu,LI Qing-bin,et al.Hex-meshing method for gravity dam based on planar Boolean operations[J].Journal of Hydroelectric Engineering,2012,31(5):209-215.(in Chinese))

    [16] 袁自立,馬福恒,焦延濤.石漫灘碾壓混凝土重力壩溫度效應(yīng)分析[J].南水北調(diào)與水利科技,2013(5):61-64.(YUAN Zi-li,MA Fu-heng,JIAO Yan-tao.Analysis of temperature effects of Shimantan RCC gravity dam[J].South-to-North Water Transfers and Water Science & Technology,2013(5):61-64.(in Chinese))DOI:10.3724/SP.J.1201.2013.05061

    [17] 童偉,張連明,范書(shū)立.官地碾壓混凝土重力壩的抗震分析[C]中國(guó)碾壓混凝土筑壩技術(shù).2010.(TONG Wei,ZHANG Lian-ming,F(xiàn)AN Shu-li.Antiseismic analysis for Guandi RCC gravity dam[C].RCC Damming Technology of China,2010.(in Chinese))

    [18] 李守義,肖靜.河床式水電站廠房橫縫止水布置形式研究[J].水力發(fā)電學(xué)報(bào),2014,33(5):165-168.(LI Shou-yi,XIAO Jing.Study on waterstop arrangements of transverse joints of hydropower house in river channel[J].Journal of Hydroelectric Engineering,2014,33(5):165-168.(in Chinese))

    [19] 朱伯芳,高季章,陳祖煜,等.拱壩設(shè)計(jì)與研究[M].北京:中國(guó)水利水電出版社,2002.(ZHU Bo-fang,GAO Ji-zhang,CHEN Zu-yu,et al.Design and research of arch dam[M].Beijing:China Water & Power Press,2002.(in Chinese))

    [20] 楊強(qiáng),劉福深,周維垣.基于直接內(nèi)力法的拱壩建基面等效應(yīng)力分析[J].水力發(fā)電學(xué)報(bào),2006,25(1):19-23.(YANG Qiang,LIU Fu-shen,ZHOU Wei-yuan.Equivalent stress analysis for the base surfaces of arch dams by direct internal force method[J].Journal of Hydroelectric Engineering,2006,25(1):19-23.(in Chinese))

    [21] 覃海藝,馬寧.薄壁彎曲結(jié)構(gòu)有限元計(jì)算精度研究[J].應(yīng)用力學(xué)學(xué)報(bào),2015,32(1):139-145.(QIN Hai-yi,MA Ning.Finite element calculating precision of thin-wall bending structure[J].Chinese Journal of Applied Mechanics,2015,32(1):139-145.(in Chinese))DOI:10.11776/cjam.32.01.B109

    猜你喜歡
    誤差分析有限元法
    正交各向異性材料裂紋疲勞擴(kuò)展的擴(kuò)展有限元法研究
    氣壓高度計(jì)的測(cè)量誤差分析及修正方法
    基于“比值比較法”確定最佳實(shí)驗(yàn)方案的研究
    云錫老廠錫礦貫通測(cè)量分析研究
    電學(xué)計(jì)量的誤差分析及不確定度理論探究分析
    車載雷達(dá)陣面舉升裝置的位置誤差
    分光計(jì)調(diào)節(jié)方法改進(jìn)及實(shí)驗(yàn)誤差分析
    三維有限元法在口腔正畸生物力學(xué)研究中發(fā)揮的作用
    集成對(duì)稱模糊數(shù)及有限元法的切削力預(yù)測(cè)
    有限元法在機(jī)械設(shè)計(jì)方向中的教學(xué)實(shí)踐
    国产在线观看jvid| 操美女的视频在线观看| 国产精品,欧美在线| 亚洲成人久久性| 99久久国产精品久久久| 亚洲性夜色夜夜综合| 啦啦啦免费观看视频1| 一二三四在线观看免费中文在| 精品第一国产精品| 欧美黄色片欧美黄色片| 国产野战对白在线观看| 精品不卡国产一区二区三区| 亚洲午夜精品一区,二区,三区| 久热爱精品视频在线9| 亚洲成人免费电影在线观看| 欧美+亚洲+日韩+国产| 成人国产一区最新在线观看| 少妇被粗大的猛进出69影院| 亚洲成av片中文字幕在线观看| 99国产综合亚洲精品| 亚洲精华国产精华精| 99精品久久久久人妻精品| 女警被强在线播放| 亚洲精品国产一区二区精华液| 国产日韩一区二区三区精品不卡| 精品久久久久久久人妻蜜臀av | 狂野欧美激情性xxxx| 成人18禁在线播放| 亚洲国产日韩欧美精品在线观看 | 一边摸一边抽搐一进一出视频| 中文字幕人妻丝袜一区二区| 国产av一区二区精品久久| 精品一区二区三区视频在线观看免费| 午夜福利影视在线免费观看| 两人在一起打扑克的视频| 久久九九热精品免费| 老鸭窝网址在线观看| 极品教师在线免费播放| 啦啦啦观看免费观看视频高清 | 制服诱惑二区| 久久性视频一级片| 久久久久国内视频| 国产成人精品久久二区二区91| 琪琪午夜伦伦电影理论片6080| 在线永久观看黄色视频| 亚洲激情在线av| 国产真人三级小视频在线观看| 亚洲自拍偷在线| 国产精品综合久久久久久久免费 | 亚洲精华国产精华精| 欧美黄色片欧美黄色片| 麻豆久久精品国产亚洲av| 亚洲国产精品合色在线| 亚洲情色 制服丝袜| 十分钟在线观看高清视频www| 一边摸一边抽搐一进一小说| 免费人成视频x8x8入口观看| 国产成人系列免费观看| 乱人伦中国视频| 悠悠久久av| 国产精品亚洲美女久久久| 天堂影院成人在线观看| 亚洲精品国产色婷婷电影| 亚洲成a人片在线一区二区| 成人手机av| 午夜福利,免费看| 美女午夜性视频免费| 午夜精品国产一区二区电影| 熟妇人妻久久中文字幕3abv| 欧美国产日韩亚洲一区| 成年人黄色毛片网站| 精品国产一区二区三区四区第35| 老司机午夜十八禁免费视频| or卡值多少钱| 国产精品免费视频内射| 9191精品国产免费久久| 成人18禁高潮啪啪吃奶动态图| 99国产精品一区二区蜜桃av| 91精品三级在线观看| 乱人伦中国视频| 激情在线观看视频在线高清| 午夜久久久在线观看| 国产精品美女特级片免费视频播放器 | 一进一出好大好爽视频| av视频免费观看在线观看| 999久久久精品免费观看国产| 色老头精品视频在线观看| 成人国产综合亚洲| 啦啦啦观看免费观看视频高清 | 一级a爱视频在线免费观看| 18禁裸乳无遮挡免费网站照片 | 性少妇av在线| 一区二区日韩欧美中文字幕| 日本免费一区二区三区高清不卡 | 久久人人97超碰香蕉20202| 欧美激情久久久久久爽电影 | 国产成人av激情在线播放| 亚洲欧洲精品一区二区精品久久久| 国产区一区二久久| av福利片在线| 十八禁人妻一区二区| 99久久国产精品久久久| 久久 成人 亚洲| 国产精品乱码一区二三区的特点 | 亚洲少妇的诱惑av| 色综合亚洲欧美另类图片| 99在线人妻在线中文字幕| 两个人视频免费观看高清| 国产成人av教育| av天堂久久9| 欧美大码av| 成人永久免费在线观看视频| 又黄又爽又免费观看的视频| 亚洲精品中文字幕在线视频| 国产1区2区3区精品| 欧美绝顶高潮抽搐喷水| 国产精品综合久久久久久久免费 | 欧美日本视频| 侵犯人妻中文字幕一二三四区| 此物有八面人人有两片| 十八禁网站免费在线| 欧美在线一区亚洲| 国产av又大| 亚洲色图 男人天堂 中文字幕| 满18在线观看网站| 男女下面插进去视频免费观看| 亚洲av第一区精品v没综合| 18美女黄网站色大片免费观看| 高清毛片免费观看视频网站| 高清毛片免费观看视频网站| 一区二区三区精品91| 免费在线观看黄色视频的| 亚洲五月婷婷丁香| bbb黄色大片| 久久久国产成人精品二区| 亚洲免费av在线视频| 国内精品久久久久久久电影| 国产欧美日韩综合在线一区二区| 免费看美女性在线毛片视频| 国产成人精品在线电影| 日韩欧美三级三区| www.999成人在线观看| 国产xxxxx性猛交| 波多野结衣av一区二区av| 国产精品久久久久久亚洲av鲁大| 亚洲五月色婷婷综合| 在线观看免费视频日本深夜| 一边摸一边做爽爽视频免费| 可以在线观看毛片的网站| 男男h啪啪无遮挡| 老司机在亚洲福利影院| 久久久国产欧美日韩av| 久久中文字幕一级| 亚洲成a人片在线一区二区| av天堂久久9| 欧美不卡视频在线免费观看 | 老司机在亚洲福利影院| 国产成+人综合+亚洲专区| 99久久99久久久精品蜜桃| 中文字幕av电影在线播放| 岛国在线观看网站| 欧美一区二区精品小视频在线| 神马国产精品三级电影在线观看 | 亚洲精品美女久久久久99蜜臀| 国产一区在线观看成人免费| 久久天堂一区二区三区四区| 中文字幕人妻丝袜一区二区| 欧美日韩亚洲综合一区二区三区_| 99久久精品国产亚洲精品| 女人精品久久久久毛片| www.熟女人妻精品国产| 精品无人区乱码1区二区| 黄色片一级片一级黄色片| 午夜免费激情av| 国产成人一区二区三区免费视频网站| 欧美丝袜亚洲另类 | 中文亚洲av片在线观看爽| 校园春色视频在线观看| 波多野结衣巨乳人妻| 91字幕亚洲| 日韩精品免费视频一区二区三区| 国产一区二区三区在线臀色熟女| 真人一进一出gif抽搐免费| 久久久国产成人精品二区| 亚洲avbb在线观看| 日韩欧美国产在线观看| 欧美乱码精品一区二区三区| 激情视频va一区二区三区| 久久久久久免费高清国产稀缺| 亚洲最大成人中文| 午夜两性在线视频| 波多野结衣高清无吗| 无遮挡黄片免费观看| 啦啦啦免费观看视频1| 一二三四社区在线视频社区8| 日韩视频一区二区在线观看| 深夜精品福利| 波多野结衣高清无吗| 免费看a级黄色片| 亚洲国产欧美日韩在线播放| 男女午夜视频在线观看| 身体一侧抽搐| 精品国产美女av久久久久小说| 国产蜜桃级精品一区二区三区| 亚洲av日韩精品久久久久久密| 日韩 欧美 亚洲 中文字幕| 久久久久久久午夜电影| 成人免费观看视频高清| 午夜福利视频1000在线观看 | 亚洲全国av大片| 久久亚洲精品不卡| 性色av乱码一区二区三区2| 久久精品国产亚洲av高清一级| 自拍欧美九色日韩亚洲蝌蚪91| 女同久久另类99精品国产91| 精品乱码久久久久久99久播| 亚洲黑人精品在线| 两人在一起打扑克的视频| 亚洲中文日韩欧美视频| 变态另类成人亚洲欧美熟女 | 最近最新中文字幕大全免费视频| 大香蕉久久成人网| 亚洲欧美日韩高清在线视频| 国产精品1区2区在线观看.| 啦啦啦免费观看视频1| 老鸭窝网址在线观看| 亚洲欧美激情综合另类| 中文字幕高清在线视频| 国产精品综合久久久久久久免费 | 国产精品电影一区二区三区| 精品第一国产精品| 好看av亚洲va欧美ⅴa在| 午夜激情av网站| 黄色 视频免费看| 黑人巨大精品欧美一区二区蜜桃| 亚洲精品一区av在线观看| 久久人人97超碰香蕉20202| 精品少妇一区二区三区视频日本电影| 老熟妇乱子伦视频在线观看| 国产一区二区三区视频了| 久久人人精品亚洲av| 男女做爰动态图高潮gif福利片 | 一区二区三区精品91| 国产精品自产拍在线观看55亚洲| 69av精品久久久久久| 久久九九热精品免费| bbb黄色大片| 嫁个100分男人电影在线观看| 亚洲人成电影免费在线| 可以在线观看的亚洲视频| 性欧美人与动物交配| 免费在线观看亚洲国产| 桃色一区二区三区在线观看| svipshipincom国产片| 国产精品一区二区在线不卡| 国产成人影院久久av| av在线播放免费不卡| 午夜免费观看网址| 夜夜爽天天搞| 日韩大码丰满熟妇| 1024香蕉在线观看| 免费看美女性在线毛片视频| 91麻豆精品激情在线观看国产| 久久久久久久午夜电影| 国产成人精品无人区| 后天国语完整版免费观看| 18禁黄网站禁片午夜丰满| 久久久久精品国产欧美久久久| 国产精品一区二区免费欧美| 欧美日韩瑟瑟在线播放| 久久人人97超碰香蕉20202| 欧美国产精品va在线观看不卡| 午夜激情av网站| 午夜福利18| 露出奶头的视频| 老司机午夜十八禁免费视频| 在线播放国产精品三级| 99国产精品99久久久久| 亚洲成人免费电影在线观看| 97碰自拍视频| 国产精品国产高清国产av| 伊人久久大香线蕉亚洲五| 两个人免费观看高清视频| 一二三四社区在线视频社区8| 久久九九热精品免费| 无人区码免费观看不卡| 国产精品一区二区在线不卡| 亚洲 欧美一区二区三区| 91精品国产国语对白视频| 咕卡用的链子| 欧美亚洲日本最大视频资源| 国产蜜桃级精品一区二区三区| 69av精品久久久久久| 欧美成人性av电影在线观看| 精品欧美国产一区二区三| 非洲黑人性xxxx精品又粗又长| 法律面前人人平等表现在哪些方面| 一级毛片高清免费大全| 九色亚洲精品在线播放| 涩涩av久久男人的天堂| 日韩欧美三级三区| 在线观看免费日韩欧美大片| 精品午夜福利视频在线观看一区| 正在播放国产对白刺激| 色婷婷久久久亚洲欧美| 激情视频va一区二区三区| 国产精品影院久久| 一二三四社区在线视频社区8| 69av精品久久久久久| 亚洲成国产人片在线观看| 精品国产乱子伦一区二区三区| 成人国产综合亚洲| 99国产极品粉嫩在线观看| 免费高清视频大片| 久久久久久久精品吃奶| 操美女的视频在线观看| 日韩一卡2卡3卡4卡2021年| 欧美一级毛片孕妇| 怎么达到女性高潮| 亚洲精品国产区一区二| 少妇粗大呻吟视频| 一区二区日韩欧美中文字幕| 午夜福利免费观看在线| 桃红色精品国产亚洲av| 亚洲av电影不卡..在线观看| 精品第一国产精品| 一本大道久久a久久精品| 天天躁夜夜躁狠狠躁躁| 国产精品久久久久久精品电影 | 在线观看免费视频日本深夜| 亚洲电影在线观看av| 最好的美女福利视频网| 在线观看一区二区三区| 午夜免费观看网址| 久久精品aⅴ一区二区三区四区| 黄频高清免费视频| av片东京热男人的天堂| 国内精品久久久久精免费| 91九色精品人成在线观看| 91字幕亚洲| 国产精品一区二区三区四区久久 | 不卡一级毛片| 欧美老熟妇乱子伦牲交| 亚洲欧美一区二区三区黑人| tocl精华| 可以免费在线观看a视频的电影网站| 最新美女视频免费是黄的| 成人永久免费在线观看视频| 午夜福利高清视频| 免费女性裸体啪啪无遮挡网站| 日韩欧美在线二视频| 精品日产1卡2卡| 国产一级毛片七仙女欲春2 | 国产野战对白在线观看| 中文字幕另类日韩欧美亚洲嫩草| 久久人妻av系列| 97碰自拍视频| 午夜福利在线观看吧| 日日摸夜夜添夜夜添小说| 亚洲激情在线av| 在线av久久热| 日本免费一区二区三区高清不卡 | 亚洲国产高清在线一区二区三 | 18禁国产床啪视频网站| 日本vs欧美在线观看视频| 精品少妇一区二区三区视频日本电影| 91av网站免费观看| 久久人妻av系列| 国产成人系列免费观看| 国产亚洲av嫩草精品影院| 香蕉久久夜色| 日韩欧美一区二区三区在线观看| 日韩一卡2卡3卡4卡2021年| 欧美日韩精品网址| 91麻豆精品激情在线观看国产| 一级作爱视频免费观看| 在线十欧美十亚洲十日本专区| 日韩欧美三级三区| 热99re8久久精品国产| 国产熟女xx| 人人妻,人人澡人人爽秒播| 大香蕉久久成人网| 国产成+人综合+亚洲专区| 一进一出好大好爽视频| 久久久精品国产亚洲av高清涩受| 亚洲一区中文字幕在线| 日韩欧美免费精品| 在线播放国产精品三级| 麻豆久久精品国产亚洲av| 亚洲中文av在线| 美女 人体艺术 gogo| 国产免费男女视频| 91麻豆精品激情在线观看国产| 国产成人一区二区三区免费视频网站| 亚洲专区字幕在线| 免费少妇av软件| 久久久精品国产亚洲av高清涩受| 日韩国内少妇激情av| 国产亚洲精品久久久久久毛片| 日韩欧美一区视频在线观看| 亚洲人成伊人成综合网2020| 午夜激情av网站| 亚洲av日韩精品久久久久久密| 国产人伦9x9x在线观看| 日本五十路高清| 久久天躁狠狠躁夜夜2o2o| av电影中文网址| 国产av精品麻豆| 亚洲情色 制服丝袜| 亚洲精品在线观看二区| 日本一区二区免费在线视频| 侵犯人妻中文字幕一二三四区| 精品免费久久久久久久清纯| 丁香六月欧美| 18禁国产床啪视频网站| 欧美丝袜亚洲另类 | 欧美不卡视频在线免费观看 | 久久精品亚洲熟妇少妇任你| 欧美日本中文国产一区发布| 亚洲第一欧美日韩一区二区三区| 亚洲电影在线观看av| 免费女性裸体啪啪无遮挡网站| 69av精品久久久久久| 国产精品美女特级片免费视频播放器 | 香蕉国产在线看| 欧美绝顶高潮抽搐喷水| 国产精品一区二区三区四区久久 | 又大又爽又粗| 少妇熟女aⅴ在线视频| 99在线人妻在线中文字幕| 精品日产1卡2卡| 亚洲全国av大片| 免费看十八禁软件| 啦啦啦韩国在线观看视频| 在线观看免费视频网站a站| 69精品国产乱码久久久| 又紧又爽又黄一区二区| 黄频高清免费视频| 亚洲国产欧美网| 99国产精品免费福利视频| 久久伊人香网站| 久久国产精品男人的天堂亚洲| 99在线视频只有这里精品首页| 两性午夜刺激爽爽歪歪视频在线观看 | 黄色a级毛片大全视频| 亚洲av成人av| 精品国产国语对白av| 亚洲九九香蕉| 午夜激情av网站| 视频在线观看一区二区三区| 欧美成人午夜精品| 亚洲人成电影免费在线| 亚洲 国产 在线| 免费在线观看亚洲国产| 精品国产一区二区三区四区第35| 亚洲天堂国产精品一区在线| 亚洲熟妇中文字幕五十中出| 成人精品一区二区免费| 亚洲一码二码三码区别大吗| 一级黄色大片毛片| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲欧美日韩无卡精品| 国产成人啪精品午夜网站| 亚洲,欧美精品.| 窝窝影院91人妻| 国产av精品麻豆| 女性被躁到高潮视频| 18美女黄网站色大片免费观看| 一区二区三区激情视频| 欧美日韩福利视频一区二区| 日本五十路高清| www.精华液| 91麻豆精品激情在线观看国产| 亚洲色图综合在线观看| 日本免费a在线| 亚洲午夜精品一区,二区,三区| 久久性视频一级片| 亚洲伊人色综图| 成人三级做爰电影| 亚洲激情在线av| 亚洲色图综合在线观看| 国产高清有码在线观看视频 | 黄片播放在线免费| 在线天堂中文资源库| 制服人妻中文乱码| 久久久国产精品麻豆| 日韩国内少妇激情av| 免费看十八禁软件| 波多野结衣高清无吗| 欧美成人免费av一区二区三区| 男女做爰动态图高潮gif福利片 | 黄片小视频在线播放| 国产1区2区3区精品| 搡老妇女老女人老熟妇| 在线观看www视频免费| 国产精品美女特级片免费视频播放器 | 真人一进一出gif抽搐免费| 神马国产精品三级电影在线观看 | 久久久国产成人精品二区| 精品人妻在线不人妻| 免费在线观看日本一区| 日本免费a在线| 黑人巨大精品欧美一区二区mp4| 一级作爱视频免费观看| 中文字幕人成人乱码亚洲影| 久久久久久久久免费视频了| 黄色视频,在线免费观看| 一二三四在线观看免费中文在| 一区在线观看完整版| 欧美av亚洲av综合av国产av| 欧美大码av| 天天添夜夜摸| 在线观看66精品国产| 亚洲人成电影观看| 久久久久亚洲av毛片大全| 色综合欧美亚洲国产小说| 岛国视频午夜一区免费看| 精品久久久精品久久久| 久久久久国内视频| av中文乱码字幕在线| 亚洲av电影不卡..在线观看| 黄色视频不卡| 欧美乱色亚洲激情| 午夜免费激情av| av中文乱码字幕在线| 校园春色视频在线观看| 99在线人妻在线中文字幕| 搞女人的毛片| 欧美乱妇无乱码| 国产亚洲精品第一综合不卡| 精品久久久久久成人av| 黄色视频不卡| 777久久人妻少妇嫩草av网站| 午夜福利视频1000在线观看 | 国产亚洲欧美精品永久| 国产欧美日韩精品亚洲av| 熟妇人妻久久中文字幕3abv| tocl精华| 久久国产乱子伦精品免费另类| 国产亚洲欧美精品永久| 精品人妻在线不人妻| 成人欧美大片| 中亚洲国语对白在线视频| 老汉色∧v一级毛片| www国产在线视频色| 男女下面进入的视频免费午夜 | 国产精品亚洲美女久久久| 90打野战视频偷拍视频| 一二三四在线观看免费中文在| 国产精品秋霞免费鲁丝片| 国产蜜桃级精品一区二区三区| 看片在线看免费视频| 最近最新免费中文字幕在线| 亚洲全国av大片| 国产精品亚洲av一区麻豆| 亚洲第一青青草原| 中文字幕精品免费在线观看视频| 波多野结衣一区麻豆| 亚洲第一av免费看| 人人澡人人妻人| 欧美日韩乱码在线| 999久久久精品免费观看国产| 欧美一级毛片孕妇| АⅤ资源中文在线天堂| 国产乱人伦免费视频| 精品无人区乱码1区二区| 国产一卡二卡三卡精品| 搞女人的毛片| 中文字幕人妻熟女乱码| 婷婷六月久久综合丁香| 久久久久久久久中文| 十八禁人妻一区二区| 亚洲欧美日韩无卡精品| 99精品在免费线老司机午夜| 亚洲精华国产精华精| 免费看a级黄色片| 精品一区二区三区av网在线观看| 亚洲av熟女| 妹子高潮喷水视频| 亚洲精品国产一区二区精华液| videosex国产| 成人特级黄色片久久久久久久| 久久久久国产一级毛片高清牌| 亚洲av成人av| 国产精品98久久久久久宅男小说| 男女下面插进去视频免费观看| www.自偷自拍.com| 国产精品久久电影中文字幕| 国产一区二区三区在线臀色熟女| 激情视频va一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 国产亚洲欧美在线一区二区| 亚洲色图av天堂| 桃色一区二区三区在线观看| 亚洲精品美女久久久久99蜜臀| 免费无遮挡裸体视频| 欧美色欧美亚洲另类二区 | 国产亚洲精品av在线| 如日韩欧美国产精品一区二区三区| 一区福利在线观看| 女人爽到高潮嗷嗷叫在线视频| 69精品国产乱码久久久| 免费av毛片视频| 天天一区二区日本电影三级 | 一区二区日韩欧美中文字幕| 女人爽到高潮嗷嗷叫在线视频| 久9热在线精品视频| 一区福利在线观看| 好看av亚洲va欧美ⅴa在| 久久久国产成人免费| 欧美在线黄色| 亚洲av片天天在线观看| 如日韩欧美国产精品一区二区三区| 国产高清激情床上av| 国产成+人综合+亚洲专区|