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

    考慮折射影響的攝影測量法土工三軸模型精度分析

    2023-01-12 07:05:50唐怡懷黃少染牟春梅
    中國測試 2022年12期
    關鍵詞:土樣徑向修正

    唐怡懷, 黃少染, 牟春梅,2

    (1. 桂林理工大學土木與建筑工程學院,廣西 桂林 541004; 2. 廣西巖土力學與工程重點實驗室,廣西 桂林 541004)

    0 引 言

    常規(guī)三軸試驗由于試驗裝置及方法的局限,只能獲取土樣整體的平均變形[1],無法詳盡地反映土樣的變形特征,且難以用于非飽和土體的變形測量[2]。針對以上問題,國內(nèi)外學者將數(shù)字圖像測量技術用于土工三軸試驗中[3],通過圖像測量手段實現(xiàn)對三軸試驗過程中土樣的變形測量?;趫D像測量的三軸試驗方法不受土樣必須完全飽和的條件限制,還可以獲取土樣的局部化應變。目前常見的技術手段有數(shù)字圖像分析技術(digital image analysis technology,DIA)[4-5]、數(shù)字圖像相關技術 (digital image correlation technology,DIC)[6-7]、邵龍?zhí)堆芯繄F隊的三軸試樣全表面變形數(shù)字圖像測量系統(tǒng)[8-10]以及Li 等[11-15]的攝影測量技術?;贒IA技術的三軸試驗是通過識別圖像中土樣的邊緣輪廓,獲取土樣軸向變形與徑向變形。該方法假設條件在實際情況中難以滿足,并且必須準確知曉壓力室及土樣與相機的相對方位[15],實際的測量精度難以達到預期水準?;?D-DIC技術的三軸試驗通過關聯(lián)處理土樣變形前后的圖像,得到不同圖像中同名點的變化實現(xiàn)土樣表面局部化變形測量。但此方法僅局限于土樣局部化平面測量,無法實現(xiàn)對土樣的全表面立體觀測。隨著圖像技術的發(fā)展,在2D-DIC技術的基礎上,結(jié)合雙目立體視覺原理形成了3DDIC技術[16],運用雙相機從不同空間位置同時觀測目標物表面一點,根據(jù)其在兩幅同時刻的圖像中像素點以及雙相機的內(nèi)外參數(shù)[17-18],確定其在空間中的三維坐標。然而,由于3D-DIC技術不能有效解決光線在壓力室和水的折射問題,難以應用于土工三軸試驗中。全表面變形數(shù)字圖像測量系統(tǒng)的三軸試驗利用亞像素點檢測原理,追蹤土樣表面各角點的位置變化,實現(xiàn)土樣表面任意時刻的應變場繪制。然而此方法需將圓筒狀壓力室改制成平面板狀從而不適用于常規(guī)三軸試驗,并且改制后的平面板狀壓力室仍無法消除壓力室與水產(chǎn)生的折射影響。攝影測量法的三軸試驗一定程度克服了上述圖像測量方法的局限性,無需對常規(guī)三軸試驗儀器進行任何改造,可在常規(guī)三軸試驗基礎上進行。該方法基于近景攝影測量原理,手持一部普通數(shù)碼相機(單鏡頭)采集土樣全表面二維圖像,隨后利用多視角下的三軸土樣圖像重構土樣的三維空間形態(tài),實現(xiàn)對三軸土樣的軸向應變、徑向應變、局部應變及體變的全過程測量。此外,攝影測量法建立的折射修正模型,考慮了壓力室受圍壓作用下出現(xiàn)鼓脹變形導致折射界面的改變對折射修正結(jié)果產(chǎn)生的影響,可實現(xiàn)高精度的折射修正。

    本文將近景攝影測量法用于土工三軸試驗,選用表面附有編碼點的硬質(zhì)圓柱作為三軸模型。通過基于攝影測量的三軸模擬試驗,建立折射修正模型對折射誤差進行校正,獲取空氣中、壓力室并注水及折射修正后三種情況下模型軸向值與徑向值。將以上測量值與空氣中三軸模型游標卡尺測量值進行比較,評估在三軸試驗中攝影測量軸向與徑向的測量精度。同時比較三軸模型折射修正前后的形態(tài)異同,探究了折射對三軸模型成像的影響及軸向與徑向成像產(chǎn)生差異的原因。最后,為進一步提高攝影測量法三軸土樣變形測量的精度提供一定借鑒。

    1 攝影測量法的基本原理

    攝影測量法是基于近景攝影測量的原理,利用相機從至少兩個不同方位拍攝被測目標點的圖像,通過軟件識別不同圖像中的同名點,獲取不同圖像中同名點的圖像像素坐標,根據(jù)不同圖像的視覺交匯,確定該目標點的唯一空間點。具體的方法[11-12]如下:以常規(guī)三軸試驗為基礎,在加載框架、壓力室外表面及橡皮膜上張貼標記點RAD(ringed automatically detected)編碼點(詳見圖 1),每個編碼點都賦有唯一的序列編號,可被PhotoModelerScanner圖像處理軟件自動識別;利用校正表和軟件檢校模塊校正測試鏡頭,手持校正后的單鏡頭相機對壓力室內(nèi)部被攝物體進行環(huán)繞式拍攝,將獲取的圖像分組輸入PhotoModelerScanner圖像處理軟件中進行識別與匹配;運用攝影測量分析技術確定拍攝相機的方位以及壓力室的形態(tài)與位置,隨后利用光線追蹤技術處理多重介質(zhì)(水、有機玻璃壓力室、空氣)引起的折射問題,并運用最小二乘法確定被攝物體表面編碼點的最終三維坐標,重建加載框架、壓力室和被攝物體的3D坐標模型,進而根據(jù)編碼點位置的變化計算被攝物體的體積變化,實現(xiàn)被攝物體的局部與整體變形測量。

    在三軸試驗攝影測量過程中,成像光線在多重介質(zhì)的交界處折射產(chǎn)生偏移。為確定光線折射后的方向,在折射界面處運用折射定理即可根據(jù)入射光的入射角計算出光線折射角。但是在三軸試驗過程中,有機玻璃壓力室受圍壓作用下發(fā)生蠕變呈現(xiàn)微鼓脹形態(tài),此時的折射界面與標準圓柱狀壓力室的折射界面不符。因此,如何擬合與壓力室實際形態(tài)相吻合的函數(shù)界面對恢復光線傳播路徑至關重要。針對上述問題,攝影測量法通過建立基于光線追蹤的折射修正模型校正折射誤差,獲取趨于被攝物體原型的三維模型。折射修正模型是基于壓力室玻璃材質(zhì)均勻的假設前提下進行,其過程主要分成兩部分,一是有機玻璃壓力室形態(tài)和位置的確定,二是編碼點坐標的修正。圓柱形壓力室的形態(tài)和位置由含有9個未知參數(shù)的數(shù)學模型表征,模型參數(shù)根據(jù)壓力室表面編碼點的三維點云坐標使用最小二乘法求得最優(yōu)的參數(shù)組合,擬合與實際相吻合的折射界面。編碼點坐標的修正即利用多條光線追蹤技術和最小二乘法消除折射對圖像測量結(jié)果的影響。光線追蹤利用光線的可逆性,通過像點逆向追蹤成像光線的傳播路徑確定物點坐標[19-20],但由于外部因素的不可控性,多條逆向光線難以按照理想狀態(tài)收斂成一點。為此,運用最小二乘法確定逆向光線最終交匯點,實現(xiàn)對編碼點坐標的準確定位。本文采用Li[11-12]所提出的基于光線追蹤的折射修正模型,并在此基礎上,將上述折射處理過程集合成程序,通過試驗分析攝影測量法建立的折射修正模型處理折射問題的準確性。

    2 試驗方案及試驗數(shù)據(jù)獲取

    2.1 試驗方案

    本次試驗選擇一個形態(tài)規(guī)則的硬質(zhì)圓柱作為三軸土樣模型(以下稱為三軸模型),模擬三軸土樣在空氣中以及充水后壓力室內(nèi)部的狀態(tài)。分別采用精度為0.02 mm的游標卡尺和攝影測量法對空氣中和壓力室內(nèi)部的三軸模型進行測量,將兩種方法所得軸向高度H和徑向直徑D進行對比分析。

    2.2 試驗準備

    2.2.1 模型標記點的設置

    從PhotoModelerScanner圖像處理軟件中導出RAD編碼點并打印在A4紙上,隨后將若干RAD編碼點粘貼在三軸模型表面(呈9行×15列分布,共135個編碼點)。另外對每行、列進行編號,以便于后續(xù)試驗的圖像采集及結(jié)果分析,如圖1所示。

    圖1 三軸模型標記點的設置

    2.2.2 試驗系統(tǒng)的建立

    在攝影空間內(nèi)定義統(tǒng)一的世界坐標系與模型比例對重構三維模型具有十分重要意義,可實現(xiàn)不同重構模型之間的變形對比。為此,在三軸儀兩側(cè)加載框架上分別粘貼一列編碼點,世界坐標系通過加載框架上編碼點定義,構建以水平方向為X軸、豎直方向為Y軸的世界坐標系O-XYZ,模型比例通過游標卡尺測得加載架上點間距確定。隨后在壓力室外表面粘貼四行和四列編碼點,并將附有編碼點的三軸模型置于壓力室中心,試驗系統(tǒng)設置如圖2所示。

    圖2 試驗系統(tǒng)的設置

    2.2.3 相機校正

    本次試驗采用手機(f/2.2光圈,800萬像素5鏡式鏡頭)進行攝影測量。同時選用手機具備:1)高分辨率、防抖、輕巧等特點,適用于室內(nèi)開展圖像采集;2)其后置鏡頭為單攝像鏡頭,適用于攝影測量領域。然而,普通相機由于內(nèi)方位元素未知以及鏡頭存在光學畸變,無法直接用于攝影測量。為了獲取相機內(nèi)方位元素以及多種畸變參數(shù),首先運用測試鏡頭環(huán)繞拍攝等距網(wǎng)格校正表,共12幅圖像,然后通過PhotoModelerScanner圖像處理軟件的檢校模塊完成上述參數(shù)的求解[21]。等距校正表及手機拍攝方位如圖3所示,相機校正結(jié)果詳見表1。

    圖3 等距校正表及手機拍攝方位圖

    表1 手機校正前后參數(shù)

    從表1可見,除相機的橫向像素M、縱向像素N以及傳感器尺寸Fy外,其余參數(shù)與相機的標定參數(shù)皆存在明顯的差異。K1、K2表示徑向畸變參數(shù),P1和P2表示切向畸變參數(shù),在校正前相機的畸變參數(shù)均不為0,說明原圖像存在一定程度的圖像畸變,通過上述畸變參數(shù)可將呈現(xiàn)桶狀畸變或枕狀畸變的圖像恢復正常。相機的校正質(zhì)量由像點的殘差值決定,通過軟件檢校模塊處理12幅圖像獲取校正表的空間點云圖,見圖3(b)。當最大像素殘差值在1.0 pixel及以下時,表明相機校正成功。

    2.3 試驗數(shù)據(jù)的獲取

    利用游標卡尺對三軸模型附有編碼點區(qū)域內(nèi)的徑向直徑、軸向高度進行測量,按照編號統(tǒng)計相應測量結(jié)果。隨后進行系統(tǒng)設置并在三軸儀底座上安裝三軸模型,手持校正后的相機拍攝獲取空氣中、安裝壓力室并注水兩種情況下的三軸模型圖像。具體操作步驟如下:1)首先正對附有編碼點的荷載架從不同方位拍攝5~6幅圖像,隨后環(huán)繞三軸模型拍攝25~35幅圖像。2)拍攝完成后,安裝壓力室并注滿水,以相同的方式環(huán)繞壓力室拍攝,獲取包含加載框架編碼點、壓力室外壁編碼點及模型表面編碼點的圖像,一次拍攝50幅圖像左右。此外,相鄰圖像間的重疊率應大于70%,三軸模型表面每個編碼點至少在5幅圖像上出現(xiàn)。

    運用PhotoModelerScanner圖像處理軟件對三軸模型圖像進行分組處理,輸出所有相機光心的三維坐標、壓力室外壁編碼點的三維坐標、空氣中三軸模型的三維坐標以及壓力室內(nèi)部三軸模型的三維坐標等信息。根據(jù)相機光心、壓力室表面編碼點的三維坐標與壓力室內(nèi)部三軸模型的三維坐標建立基于光線追蹤的折射修正模型,消除折射引起的放大效應,確定三軸模型表面編碼點的真實三維坐標。為便于執(zhí)行折射修正過程,在Matlab軟件編寫了相應折射修正的程序,實現(xiàn)折射誤差的快速修正。

    將三軸模型在空氣中以及壓力室內(nèi)修正前后的編碼點三維坐標按行列編號進行分類統(tǒng)計。軸向上根據(jù)三軸模型表面首行與末行的兩RAD編碼點三維坐標z值之差確定,而徑向上提取模型RAD編碼點的三維坐標中x、y值,通過擬合圓函數(shù)獲取徑向直徑,并計算平均值。

    3 試驗結(jié)果分析

    試驗結(jié)果可分為三部分,一是評估空氣中進行攝影測量的準確性,二是評估模擬試驗中攝影測量的準確性,三是折射修正前后試驗結(jié)果對比分析。

    3.1 空氣中的準確性驗證

    利用卡尺和攝影測量法對編碼點間距進行對比測量,初步評估攝影測量法的準確性。如圖4所示,將分辨率為0.02 mm的游標卡尺放在印有編碼點的A4紙上,測量編碼點(編號1~10)兩兩之間的距離;使用經(jīng)校正后的手機從不同方位拍攝,共6幅圖像,然后運用PhotoModelerScanner圖像處理軟件對圖像進行處理,獲取編號1~10相鄰兩編碼點的間距值,兩種方法測量結(jié)果見表2。

    圖4 空氣中的精度驗證試驗

    表2 游標卡尺和攝影測量法分析結(jié)果對比

    從表2可見,以卡尺測量結(jié)果作為點間距的“基準值”,攝影測量結(jié)果與卡尺測量結(jié)果基本一致,點間距測量值的最大絕對誤差和最大相對誤差分別為0.025 mm和0.381%。在9個測量值中,7個測量值的絕對誤差值小于0.01 mm。由此表明,攝影測量法用于空氣中具有較高的測量精度,將其用于變形測量具有一定可行性。

    3.2 三軸試驗的模擬試驗

    在三軸試驗攝影測量過程中,由于有機玻璃壓力室和水所產(chǎn)生的折射會引起攝影測量結(jié)果的失準,故分析時需要考慮折射對測量結(jié)果產(chǎn)生的影響,其測量精度需要進一步分析。通過基于攝影測量的三軸模擬試驗,分別獲取空氣中三軸模型的尺寸測量值(H空氣、D空氣)和注水壓力室內(nèi)三軸模型的尺寸測量值(H修正、D修正),將以上結(jié)果與空氣中三軸模型的游標卡尺測量值(H卡尺、D卡尺)進行比較,相關結(jié)果見表3、表4。

    表3 三軸模型軸向高度測量結(jié)果分析

    表4 三軸模型平均徑向直徑測量結(jié)果

    根據(jù)表3三軸模型軸向高度的測量結(jié)果分析,以卡尺測量H卡尺值作為軸向“基準值”,在空氣中的攝影測量軸向測量值與卡尺軸向測量值相近,其中最大誤差為-0.075 mm,最小誤差為-0.001 mm,平均誤差為-0.038 mm;而修正后的攝影測量軸向測量值H修正與軸向“基準值”H卡尺存在明顯差異,最大誤差為 0.193 mm,最小誤差為 0.085 mm,平均誤差為0.133 mm,相較之下H修正值的誤差比H空氣值大。同時結(jié)合表2可知,同是在空氣中進行攝影測量,三軸模擬試驗中攝影測量的準確性卻不如在空氣中精度驗證試驗中的表現(xiàn),空氣中三軸模型高度的測量誤差大于點間距的測量誤差。這是因為在處理空氣中三軸模型圖像時,各點坐標通過識別不同圖像中同名點確定,編碼點的殘差值并不完全小于1.0,故獲取的編碼點坐標存在一定的偏差;而在空氣中精度驗證試驗中,相鄰圖像間的重疊率達到100%,編碼點識別的準確性相對更高,故得到點的位置坐標更準確。

    根據(jù)表4三軸模型徑向平均直徑的測量結(jié)果分析,以卡尺測量D卡尺值作為徑向“基準值”,D空氣值的誤差為-0.005 mm,而D修正值的誤差為 0.201 mm,在徑向上折射修正后的測量誤差明顯大于攝影測量空氣中的測量誤差。造成徑向修正結(jié)果誤差較大的原因可能有兩個:一是人為因素造成行向編碼點不完全處于一條直線上,從而影響了此次圓函數(shù)的擬合;二是攝影測量在折射修正過程中目標點的定位存在偏差所致。三軸模型上每個編碼點至少在5幅圖像上出現(xiàn),理論上通過不同圖像產(chǎn)生的逆向光線形成的視覺交匯即為該編碼點的坐標位置,然而由于外部因素影響,光線無法按照理想狀態(tài)匯聚成一點。因此,運用最小二乘法解決光線無法交匯確定物點的問題,其基本原理是:假設物點為光源,那么所有光線都由物點發(fā)散,此時物點與所有光線存在最短距離即物點至逆向光線的距離di的平方和最小。在處理過程中,當di值大于臨界值0.15 mm時,即認為光線偏離距離較大,確定的物點誤差較大,應剔除對應的光線,隨后將剩余光線進行重新計算。

    圖5展示了其中5個修正點所使用的圖像數(shù)量及具體的di值。從圖上可見,折射修正過程使用的圖像數(shù)量最少有8幅,最多有11幅,獲得的編碼點坐標將更加精確。同時,此次抽選5個點的di值絕大部分處在0.04 mm以內(nèi),表明通過最小二乘法確定編碼點的最終位置相對準確。

    圖5 模型上5個修正點所用的圖像數(shù)量以及具體的di值

    3.3 折射修正前后對比分析

    為了進一步分析折射影響產(chǎn)生的放大現(xiàn)象,將折射修正前后三軸模型的形態(tài)繪制在同一圖形上,并根據(jù)點云坐標獲取折射修正前后的徑向測量平均值與軸向測量平均值,如圖6所示。同時,將修正前后的三軸模型軸向高度(H未修正、H修正)與卡尺測量值(H卡尺)進行比較,如圖7所示。

    圖6 三軸模型修正前后形態(tài)對比

    圖7 三軸模型軸向高度修正前后與卡尺測量值對比

    從圖6可見,由于多重介質(zhì)的折射影響,充水后的壓力室內(nèi)三軸模型徑向與軸向均呈現(xiàn)出折射放大現(xiàn)象。但折射對三軸模型軸向與徑向的影響不一致,相對于軸向,徑向受折射影響的作用更為顯著,與圖2實際觀察的結(jié)果相吻合。分析其原因,在徑向上壓力室玻璃罩為一弧形曲面,使得充水后的壓力室如同一面透鏡,因此三軸模型的徑向測量不僅受多重介質(zhì)的折射影響,還受到壓力室弧形曲面影響,導致成像光線在弧形曲面處發(fā)生明顯的曲折[22],成像時三軸模型的徑向直徑呈現(xiàn)明顯的放大變形現(xiàn)象,徑向測量平均值從 31.961 mm(D修正)放大至40.311 mm(D未修正),差值為 8.35 mm,其平均誤差增加26.29%;在軸向上壓力室玻璃罩如同一面平面玻璃,由于平面狀玻璃受到折射的影響較小[8],并且水、有機玻璃壓力室及空氣三者的折射率相近[15],成像光線受不同介質(zhì)的折射的曲折程度較小,使得三軸模型軸向測量時呈現(xiàn)的折射放大現(xiàn)象相對不明顯,軸向測量平均值從57.321 mm()放大至59.091 mm(),差值僅為 1.77 mm,其平均誤差增加3.09%。同時,根據(jù)圖7模型軸向高度修正前后與卡尺測量值對比的結(jié)果分析,從圖上可見,盡管軸向受折射的影響作用相對要小,但其測量值與三軸模型實際測量值仍存在明顯差異。在進行折射修正前H未修正值均大于H卡尺值,此時每列軸向高度測量值波動幅度較大,從58.9 mm附近變化至59.3 mm上下波動,而折射修正后H修正測量值曲線變化幅度較小,與H卡尺測量值曲線近似平行,二者差值明顯減小。

    由此表明,折射修正模型可以較為準確地修正折射所產(chǎn)生的放大變形,獲取趨近三軸模型原狀的三維模型。此外,成像時三軸模型軸向高度受折射影響并非呈一定比例放大,而是在某一范圍內(nèi)波動,產(chǎn)生上述現(xiàn)象的原因主要是試驗受到相機拍攝的方位與角度、光照環(huán)境穩(wěn)定程度以及玻璃材質(zhì)的均勻程度等影響,三軸模型各方向受到的折射影響不一致,折射對三軸模型成像的影響存在波動。結(jié)合表3、表4數(shù)據(jù)分析,以空氣中卡尺測量的平均值作為“基準值”,折射修正后的軸向測量平均誤差為0.23%,徑向測量平均誤差為0.63%,這與折射修正前的軸向測量平均誤差減少3.09%,徑向測量平均誤差減少26.29%。上述試驗結(jié)果表明,通過基于光線追蹤的折射修正模型可以將折射所致誤差控制在試驗允許的范圍,將考慮折射因素的攝影測量法用于土工三軸試驗,基本滿足三軸土樣變形測量的精度要求。

    4 結(jié)束語

    本文通過空氣中的精度驗證試驗初步評估攝影測量法在空氣中拍攝測量的準確性。以硬質(zhì)圓柱體作為三軸土樣模型,通過攝影測量的三軸模擬試驗對三軸試驗中攝影測量的精度進行驗證,分析壓力室與水引起的折射對三軸模型的成像影響以及折射修正模型處理折射問題的效果。同時探討將考慮折射影響的攝影測量法用于三軸試驗的適用性,得出以下結(jié)論:

    1)在空氣中進行精度驗證試驗,攝影測量結(jié)果與游標卡尺測量結(jié)果基本吻合,點間距測量值絕對誤差在0.025 mm以內(nèi),攝影測量法在空氣中具有較高的測量精度。

    2)不同介質(zhì)的折射對三軸模型的成像有重要影響,基于光線追蹤的折射修正模型可以有效減少折射引起的測量誤差,經(jīng)修正后的徑向平均誤差為0.63%,軸向平均誤差為0.23%,這與折射修正前的徑向測量平均誤差減少26.29%,而軸向測量平均誤差減少3.09%。

    3)在折射修正前后,攝影測量徑向直徑平均值從31.961 mm放大至40.311 mm,而攝影測量軸向高度平均值從57.321 mm放大至59.091 mm,從數(shù)值變化可知,徑向受折射的影響顯著大于軸向。此外,三軸模型各方向受到的折射影響各不相同,成像時三軸模型各向軸向高度在某一范圍內(nèi)波動,并非按照一定比例放大變形。

    目前攝影測量法應用研究尚且不足,測量系統(tǒng)中還存在一系列誤差,如成像系統(tǒng)誤差、折射修正參數(shù)誤差、外部環(huán)境誤差等,有待進一步完善。從以下方面提出減少相應誤差的建議與措施:1)升級相機成像采集系統(tǒng)。引起相機成像系統(tǒng)誤差的主要因素為相機參數(shù)與圖像分辨率。今后可配備具有高精度參數(shù)的相機與高圖像分辨率的相機成像采集系統(tǒng),從而提高攝影測量的精度以及圖像處理的效率。2)完善折射修正模型。折射修正模型參數(shù)是影響攝影測量法在三軸試驗測量精度的主要因素之一,該誤差是模型參數(shù)與實際值存在偏差引起的。故而今后對折射修正模型進行完善,以獲取準確性更高的模型參數(shù)組合,從而構建更精確的折射修正模型。3)研發(fā)自動化補光模塊。拍攝過程中光照強度過強或過弱會導致拍攝的圖像曝光過度或曝光不足,影響圖像的成像質(zhì)量。因此,有必要研發(fā)配備自動化補光模塊減小外部環(huán)境所致誤差。

    猜你喜歡
    土樣徑向修正
    灌區(qū)渠道基土工程水敏性試驗研究
    檸檬酸對改良紫色土中老化銅的淋洗研究
    Some new thoughts of definitions of terms of sedimentary facies: Based on Miall's paper(1985)
    修正這一天
    快樂語文(2021年35期)2022-01-18 06:05:30
    淺探徑向連接體的圓周運動
    RN上一類Kirchhoff型方程徑向?qū)ΨQ正解的存在性
    基于PID+前饋的3MN徑向鍛造機控制系統(tǒng)的研究
    重型機械(2020年3期)2020-08-24 08:31:40
    合同解釋、合同補充與合同修正
    法律方法(2019年4期)2019-11-16 01:07:28
    一類無窮下級整函數(shù)的Julia集的徑向分布
    膨脹土干濕交替作用下殘余強度試驗方案分析
    治淮(2018年6期)2018-01-30 11:42:44
    午夜免费观看网址| 欧美黄色片欧美黄色片| 99热这里只有是精品50| 日本a在线网址| 一区福利在线观看| 小说图片视频综合网站| 制服人妻中文乱码| 午夜精品久久久久久毛片777| 亚洲美女视频黄频| 亚洲av五月六月丁香网| 午夜福利成人在线免费观看| 91在线精品国自产拍蜜月 | av在线天堂中文字幕| 午夜影院日韩av| 久久久久国产精品人妻aⅴ院| 国产在视频线在精品| 男插女下体视频免费在线播放| 琪琪午夜伦伦电影理论片6080| 国产午夜精品久久久久久一区二区三区 | 看片在线看免费视频| 老司机午夜十八禁免费视频| 叶爱在线成人免费视频播放| 热99在线观看视频| 亚洲电影在线观看av| 看免费av毛片| 在线免费观看不下载黄p国产 | 99国产综合亚洲精品| av天堂在线播放| 757午夜福利合集在线观看| 国产av一区在线观看免费| 国产中年淑女户外野战色| 国产精品一及| 免费av不卡在线播放| 757午夜福利合集在线观看| 丁香六月欧美| 精品人妻一区二区三区麻豆 | 国产蜜桃级精品一区二区三区| 国产亚洲精品一区二区www| 久久精品国产综合久久久| 午夜视频国产福利| 国产三级中文精品| 一个人免费在线观看的高清视频| 久久精品国产99精品国产亚洲性色| 久久久久久久久中文| 大型黄色视频在线免费观看| 99国产综合亚洲精品| xxxwww97欧美| 两性午夜刺激爽爽歪歪视频在线观看| 啦啦啦免费观看视频1| 两个人视频免费观看高清| 欧美中文日本在线观看视频| 又黄又粗又硬又大视频| 色噜噜av男人的天堂激情| 国产精品自产拍在线观看55亚洲| 免费看十八禁软件| 免费看美女性在线毛片视频| 少妇丰满av| 国产高潮美女av| 91九色精品人成在线观看| 欧美成人一区二区免费高清观看| 麻豆国产97在线/欧美| 国产男靠女视频免费网站| 日韩人妻高清精品专区| 国产高清视频在线观看网站| 国产亚洲av嫩草精品影院| 亚洲国产精品久久男人天堂| netflix在线观看网站| 99视频精品全部免费 在线| 亚洲精品成人久久久久久| 成人特级av手机在线观看| 国产精品久久久久久人妻精品电影| 亚洲aⅴ乱码一区二区在线播放| 国产精品美女特级片免费视频播放器| 成人18禁在线播放| 亚洲午夜理论影院| 波野结衣二区三区在线 | av国产免费在线观看| 国产高清有码在线观看视频| 精品久久久久久成人av| 国内精品美女久久久久久| 十八禁网站免费在线| 中亚洲国语对白在线视频| 成人无遮挡网站| 免费看美女性在线毛片视频| 久久精品影院6| 国产精品日韩av在线免费观看| 九九久久精品国产亚洲av麻豆| 午夜影院日韩av| 国产真实伦视频高清在线观看 | 成人午夜高清在线视频| 久久精品影院6| 男人舔女人下体高潮全视频| 国产av在哪里看| 久久香蕉精品热| 又黄又爽又免费观看的视频| 国产伦在线观看视频一区| 男女那种视频在线观看| e午夜精品久久久久久久| 在线看三级毛片| 国产美女午夜福利| 一级黄片播放器| av女优亚洲男人天堂| 久久九九热精品免费| 久久亚洲真实| 亚洲av熟女| 欧美成人性av电影在线观看| 男女床上黄色一级片免费看| 深夜精品福利| 无遮挡黄片免费观看| 精品一区二区三区视频在线观看免费| 一区二区三区国产精品乱码| 一本精品99久久精品77| 日本撒尿小便嘘嘘汇集6| 一级黄片播放器| 国产69精品久久久久777片| 国产美女午夜福利| 亚洲国产精品成人综合色| 亚洲内射少妇av| 国产真人三级小视频在线观看| 国产av在哪里看| 国产三级中文精品| 亚洲专区中文字幕在线| 国产精品久久久久久久电影 | 亚洲,欧美精品.| 国产真人三级小视频在线观看| 久久国产乱子伦精品免费另类| 欧美+亚洲+日韩+国产| 欧美三级亚洲精品| 精品福利观看| 亚洲一区二区三区不卡视频| 国产精品野战在线观看| 午夜精品一区二区三区免费看| 国内久久婷婷六月综合欲色啪| 国产97色在线日韩免费| aaaaa片日本免费| 在线观看免费午夜福利视频| 操出白浆在线播放| 色精品久久人妻99蜜桃| 欧美日韩乱码在线| 免费看美女性在线毛片视频| or卡值多少钱| 亚洲狠狠婷婷综合久久图片| 19禁男女啪啪无遮挡网站| 哪里可以看免费的av片| 中文字幕人成人乱码亚洲影| 亚洲自拍偷在线| 国产一区在线观看成人免费| 亚洲国产精品999在线| 女人高潮潮喷娇喘18禁视频| 亚洲真实伦在线观看| 日本免费一区二区三区高清不卡| 一区二区三区激情视频| 婷婷精品国产亚洲av| 国产毛片a区久久久久| 成人精品一区二区免费| 高潮久久久久久久久久久不卡| 蜜桃久久精品国产亚洲av| 夜夜夜夜夜久久久久| 亚洲国产精品久久男人天堂| 美女被艹到高潮喷水动态| 脱女人内裤的视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 日韩中文字幕欧美一区二区| 亚洲专区国产一区二区| 欧美区成人在线视频| 国产极品精品免费视频能看的| 国产精品1区2区在线观看.| 全区人妻精品视频| 两个人的视频大全免费| 欧美一区二区精品小视频在线| 国产精品影院久久| 丰满的人妻完整版| 欧美黄色片欧美黄色片| 淫秽高清视频在线观看| 亚洲av日韩精品久久久久久密| 天天添夜夜摸| 国产精品一区二区三区四区久久| АⅤ资源中文在线天堂| 天堂√8在线中文| or卡值多少钱| 国产伦人伦偷精品视频| 久久亚洲真实| 久久性视频一级片| 国产真实乱freesex| 最新在线观看一区二区三区| 亚洲一区高清亚洲精品| 熟女少妇亚洲综合色aaa.| 久久人妻av系列| 亚洲av中文字字幕乱码综合| 波多野结衣巨乳人妻| 国产精品国产高清国产av| 欧美日韩瑟瑟在线播放| 我的老师免费观看完整版| 少妇人妻一区二区三区视频| 国产三级在线视频| 九九久久精品国产亚洲av麻豆| 国产成年人精品一区二区| 精品久久久久久久毛片微露脸| 国产黄色小视频在线观看| 最后的刺客免费高清国语| 真人一进一出gif抽搐免费| 操出白浆在线播放| 桃色一区二区三区在线观看| 亚洲久久久久久中文字幕| 国产成人a区在线观看| 日韩欧美在线乱码| 国产老妇女一区| 午夜影院日韩av| 国产精品一区二区三区四区久久| av视频在线观看入口| x7x7x7水蜜桃| 国产色爽女视频免费观看| 欧美极品一区二区三区四区| 99久久精品一区二区三区| 岛国在线观看网站| 午夜福利高清视频| 99精品欧美一区二区三区四区| 桃红色精品国产亚洲av| 亚洲av五月六月丁香网| 制服丝袜大香蕉在线| 国产精品一区二区免费欧美| avwww免费| 狂野欧美白嫩少妇大欣赏| h日本视频在线播放| 国产亚洲欧美在线一区二区| 国产在线精品亚洲第一网站| 日本 av在线| 51国产日韩欧美| 一级黄色大片毛片| 亚洲欧美日韩卡通动漫| 在线播放无遮挡| 丰满人妻一区二区三区视频av | 手机成人av网站| 久久精品综合一区二区三区| 色尼玛亚洲综合影院| 国产高潮美女av| 久久久久久大精品| 国产精品久久久久久人妻精品电影| 十八禁人妻一区二区| 好男人在线观看高清免费视频| 国产日本99.免费观看| 国内精品久久久久精免费| 国语自产精品视频在线第100页| 国产一级毛片七仙女欲春2| 国产真实伦视频高清在线观看 | 日韩有码中文字幕| 亚洲精品456在线播放app | 国产久久久一区二区三区| 97碰自拍视频| 1024手机看黄色片| 99视频精品全部免费 在线| 观看美女的网站| 人妻丰满熟妇av一区二区三区| avwww免费| 老司机午夜十八禁免费视频| 午夜福利高清视频| 亚洲欧美日韩高清专用| 欧美日韩综合久久久久久 | 亚洲精品成人久久久久久| 亚洲人成电影免费在线| 国产老妇女一区| 18禁美女被吸乳视频| 美女被艹到高潮喷水动态| 色综合亚洲欧美另类图片| 女人被狂操c到高潮| 国产精品综合久久久久久久免费| 亚洲av免费高清在线观看| 白带黄色成豆腐渣| 观看美女的网站| 天堂av国产一区二区熟女人妻| 久久久成人免费电影| 免费大片18禁| 91九色精品人成在线观看| 国产国拍精品亚洲av在线观看 | 亚洲精品成人久久久久久| 一进一出抽搐动态| 国产精品99久久99久久久不卡| 99国产综合亚洲精品| 窝窝影院91人妻| 国产91精品成人一区二区三区| 人妻久久中文字幕网| 日本五十路高清| 别揉我奶头~嗯~啊~动态视频| 变态另类成人亚洲欧美熟女| 天天躁日日操中文字幕| 成年版毛片免费区| 欧美在线一区亚洲| 18+在线观看网站| 十八禁网站免费在线| 国产在视频线在精品| 日韩 欧美 亚洲 中文字幕| 精品无人区乱码1区二区| 身体一侧抽搐| 国产精华一区二区三区| 亚洲人成网站在线播放欧美日韩| 色精品久久人妻99蜜桃| av福利片在线观看| 69人妻影院| 午夜视频国产福利| 国产精品日韩av在线免费观看| 99热这里只有是精品50| 岛国视频午夜一区免费看| 精品99又大又爽又粗少妇毛片 | 大型黄色视频在线免费观看| 午夜精品久久久久久毛片777| 国产精品久久久久久久电影 | 久久久久久九九精品二区国产| 最后的刺客免费高清国语| 激情在线观看视频在线高清| 丝袜美腿在线中文| 欧美不卡视频在线免费观看| 欧美激情在线99| 久久久久免费精品人妻一区二区| 露出奶头的视频| 窝窝影院91人妻| 国产高清视频在线观看网站| 高潮久久久久久久久久久不卡| 国产欧美日韩一区二区三| 色综合欧美亚洲国产小说| 日韩精品青青久久久久久| 久9热在线精品视频| 亚洲人成网站在线播放欧美日韩| 高清毛片免费观看视频网站| 亚洲国产高清在线一区二区三| 久久精品91无色码中文字幕| 日本成人三级电影网站| 欧美成人一区二区免费高清观看| 麻豆成人av在线观看| 免费人成视频x8x8入口观看| 免费无遮挡裸体视频| 搡老妇女老女人老熟妇| 国内精品一区二区在线观看| 国产精品久久久久久人妻精品电影| 两性午夜刺激爽爽歪歪视频在线观看| 51午夜福利影视在线观看| 久久精品国产综合久久久| 90打野战视频偷拍视频| 在线国产一区二区在线| 欧美性猛交╳xxx乱大交人| 中文字幕av在线有码专区| 有码 亚洲区| 中文字幕久久专区| 国产精品三级大全| 看黄色毛片网站| 国产一区二区激情短视频| 国内精品久久久久精免费| 美女高潮喷水抽搐中文字幕| 一级毛片女人18水好多| 久久久久久国产a免费观看| 成人一区二区视频在线观看| 国产亚洲欧美在线一区二区| 欧美高清成人免费视频www| 亚洲精品粉嫩美女一区| 看黄色毛片网站| 久久久国产精品麻豆| 成人一区二区视频在线观看| 亚洲av日韩精品久久久久久密| 看片在线看免费视频| 国产一区二区亚洲精品在线观看| 亚洲国产欧美网| 亚洲av日韩精品久久久久久密| 真人做人爱边吃奶动态| 嫩草影院精品99| 国产亚洲精品av在线| or卡值多少钱| 搞女人的毛片| 日韩av在线大香蕉| av福利片在线观看| 免费观看的影片在线观看| 成年女人看的毛片在线观看| 别揉我奶头~嗯~啊~动态视频| 两个人的视频大全免费| 国产成人av教育| 亚洲av成人不卡在线观看播放网| www.www免费av| 亚洲无线观看免费| 色在线成人网| 亚洲精华国产精华精| 内地一区二区视频在线| 成人国产综合亚洲| 国产免费av片在线观看野外av| 亚洲欧美日韩高清专用| 动漫黄色视频在线观看| 真人一进一出gif抽搐免费| 麻豆久久精品国产亚洲av| 久9热在线精品视频| 国产高清视频在线播放一区| 亚洲人成伊人成综合网2020| 成人三级黄色视频| 亚洲不卡免费看| 婷婷六月久久综合丁香| 日本一本二区三区精品| 亚洲精品国产精品久久久不卡| 窝窝影院91人妻| 欧美性猛交黑人性爽| 一区二区三区国产精品乱码| 国产精品女同一区二区软件 | 丁香六月欧美| 深爱激情五月婷婷| 国产精品国产高清国产av| 一个人观看的视频www高清免费观看| 一个人看的www免费观看视频| 一区二区三区免费毛片| 亚洲 欧美 日韩 在线 免费| 久9热在线精品视频| 亚洲专区中文字幕在线| 成人鲁丝片一二三区免费| 精华霜和精华液先用哪个| 成年免费大片在线观看| 国产精品香港三级国产av潘金莲| 国产aⅴ精品一区二区三区波| 美女高潮喷水抽搐中文字幕| 真人做人爱边吃奶动态| 欧美日韩乱码在线| 亚洲精品一卡2卡三卡4卡5卡| 一本久久中文字幕| 国内少妇人妻偷人精品xxx网站| 99热6这里只有精品| 亚洲自拍偷在线| 中出人妻视频一区二区| 99久久精品热视频| 天堂网av新在线| 国产一区二区亚洲精品在线观看| 99久久99久久久精品蜜桃| 一个人免费在线观看电影| 亚洲午夜理论影院| 女生性感内裤真人,穿戴方法视频| 亚洲欧美日韩高清在线视频| 精品久久久久久久毛片微露脸| 人妻丰满熟妇av一区二区三区| 久久婷婷人人爽人人干人人爱| 在线播放国产精品三级| 国产精品一区二区三区四区久久| 日韩欧美在线乱码| 一个人免费在线观看电影| 亚洲色图av天堂| 国模一区二区三区四区视频| 每晚都被弄得嗷嗷叫到高潮| 91麻豆精品激情在线观看国产| 欧美大码av| 国产成人福利小说| 一本综合久久免费| 色精品久久人妻99蜜桃| 12—13女人毛片做爰片一| 免费在线观看亚洲国产| 在线观看一区二区三区| 99热这里只有是精品50| 亚洲av电影在线进入| 在线十欧美十亚洲十日本专区| 亚洲精华国产精华精| 午夜福利在线观看免费完整高清在 | 又黄又爽又免费观看的视频| 欧美成人一区二区免费高清观看| 成人av一区二区三区在线看| 女人被狂操c到高潮| 久久精品国产亚洲av香蕉五月| 亚洲在线观看片| 亚洲无线在线观看| 欧美日韩国产亚洲二区| 欧美一区二区亚洲| 韩国av一区二区三区四区| 国产亚洲欧美在线一区二区| 亚洲人成网站高清观看| 色综合亚洲欧美另类图片| 91麻豆av在线| 国产野战对白在线观看| 一本一本综合久久| 两性午夜刺激爽爽歪歪视频在线观看| 国产欧美日韩精品亚洲av| 91麻豆av在线| 搡老岳熟女国产| 久久久久性生活片| 亚洲片人在线观看| 中文字幕久久专区| 观看美女的网站| 在线观看66精品国产| 日本一二三区视频观看| 在线天堂最新版资源| 欧美xxxx黑人xx丫x性爽| 一个人看视频在线观看www免费 | 少妇的丰满在线观看| 日韩免费av在线播放| 久久精品91无色码中文字幕| 在线观看舔阴道视频| 88av欧美| 欧美日韩亚洲国产一区二区在线观看| 一区二区三区高清视频在线| 99久久无色码亚洲精品果冻| 午夜视频国产福利| 久久久久久久久大av| 老司机福利观看| 精品一区二区三区视频在线观看免费| 观看美女的网站| 可以在线观看毛片的网站| 男人舔女人下体高潮全视频| 亚洲人与动物交配视频| 久久国产精品影院| 日韩欧美免费精品| 怎么达到女性高潮| 午夜免费观看网址| 亚洲avbb在线观看| 精华霜和精华液先用哪个| 亚洲一区高清亚洲精品| 国产真实伦视频高清在线观看 | av天堂中文字幕网| 欧美日韩乱码在线| 精品日产1卡2卡| 又紧又爽又黄一区二区| av在线天堂中文字幕| 精品日产1卡2卡| 欧美三级亚洲精品| 啪啪无遮挡十八禁网站| 国产午夜精品论理片| 老汉色av国产亚洲站长工具| av黄色大香蕉| 午夜福利欧美成人| 一级黄片播放器| 成人av一区二区三区在线看| 亚洲精华国产精华精| 两性午夜刺激爽爽歪歪视频在线观看| 深爱激情五月婷婷| 特级一级黄色大片| 欧美zozozo另类| 亚洲精品一卡2卡三卡4卡5卡| 18禁裸乳无遮挡免费网站照片| 国产精品99久久99久久久不卡| 丝袜美腿在线中文| 一个人看视频在线观看www免费 | 亚洲天堂国产精品一区在线| 国产高清激情床上av| 久久精品国产99精品国产亚洲性色| 国产精品久久久久久亚洲av鲁大| 男人和女人高潮做爰伦理| 在线观看66精品国产| 亚洲七黄色美女视频| netflix在线观看网站| 国产精品日韩av在线免费观看| 99久久精品一区二区三区| 欧美黑人欧美精品刺激| 在线播放国产精品三级| av欧美777| 91麻豆av在线| 亚洲自拍偷在线| 在线看三级毛片| 91久久精品电影网| 午夜福利视频1000在线观看| 国产乱人伦免费视频| 一夜夜www| 香蕉久久夜色| 淫秽高清视频在线观看| av专区在线播放| 亚洲18禁久久av| 级片在线观看| 久久精品夜夜夜夜夜久久蜜豆| 日本a在线网址| 亚洲国产色片| 91久久精品电影网| 精品一区二区三区视频在线 | 日韩欧美免费精品| 国产亚洲精品综合一区在线观看| 老汉色av国产亚洲站长工具| 丰满人妻熟妇乱又伦精品不卡| 久久精品国产自在天天线| 亚洲av电影不卡..在线观看| 亚洲国产精品999在线| 狂野欧美激情性xxxx| 欧美日本亚洲视频在线播放| 非洲黑人性xxxx精品又粗又长| 国产精品久久久久久亚洲av鲁大| 国产日本99.免费观看| ponron亚洲| 国产欧美日韩一区二区精品| 国产高潮美女av| 一进一出抽搐动态| 好男人电影高清在线观看| 国产精品美女特级片免费视频播放器| 国产亚洲精品久久久久久毛片| 男人和女人高潮做爰伦理| 一进一出抽搐gif免费好疼| 黄色片一级片一级黄色片| 女人高潮潮喷娇喘18禁视频| 成年女人毛片免费观看观看9| 最新中文字幕久久久久| 国产精品久久电影中文字幕| 在线观看一区二区三区| 熟女少妇亚洲综合色aaa.| 少妇丰满av| 精品99又大又爽又粗少妇毛片 | 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 日韩精品中文字幕看吧| 欧美zozozo另类| 手机成人av网站| 天美传媒精品一区二区| 亚洲无线观看免费| av在线天堂中文字幕| 国产免费男女视频| 欧美黑人巨大hd| 亚洲av二区三区四区| 欧美成人a在线观看| 小说图片视频综合网站| 久久久久国内视频| 欧美日韩精品网址| 男女视频在线观看网站免费| 久久国产精品影院| 十八禁人妻一区二区| 亚洲成人免费电影在线观看| 日韩欧美精品免费久久 | 亚洲人成网站高清观看| 色吧在线观看| 精品一区二区三区av网在线观看| 美女免费视频网站| 成人18禁在线播放| 免费高清视频大片| 韩国av一区二区三区四区| 亚洲av第一区精品v没综合|