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

    球面凸類圖形Delaunay三角剖分再分算法及其收斂性分析

    2018-01-08 07:33:57李映華
    計算機應(yīng)用 2017年12期
    關(guān)鍵詞:小圓圓心角剖分

    夏 俊,李映華

    (昆明理工大學 理學院,昆明 650500)

    球面凸類圖形Delaunay三角剖分再分算法及其收斂性分析

    夏 俊*,李映華

    (昆明理工大學 理學院,昆明 650500)

    在計算曲面Ricci Flow時,會因為三角網(wǎng)格中存在過小的角而出現(xiàn)不收斂的情況。針對這種不收斂的問題,提出一種提高最小角角度的球面凸類圖形Delaunay三角剖分再分算法。首先,給出球面凸類圖形Delaunay三角剖分再分算法。它的核心操作有兩個:1)如果某條Delaunay劣弧被“侵占”, 通過添加Delaunay劣弧中點分割Delaunay劣弧;2)如果存在“瘦”球面三角形, 通過添加球面三角形外接球面小圓圓心分解球面三角形。然后,利用局部特征尺度探索出所提算法的收斂條件并給出輸出頂點的一個上界公式。根據(jù)實驗輸出的網(wǎng)格驗證,所提算法網(wǎng)格生成的球面三角形沒有狹小的角,適合用來計算Ricci Flow。

    球面;Delaunay三角剖分;劣??;局部特征尺度;收斂

    0 引言

    Delaunay三角剖分在圖形學、地理學以及有限元分析等領(lǐng)域是一項極為重要的預處理技術(shù)。隨著計算機技術(shù)的發(fā)展以及分析對象的復雜性,Delaunay三角剖分在這些領(lǐng)域扮演著重要角色。例如,在有限元分析中,需要把研究的區(qū)域分解成三角形單位元。如果剖分的三角形角度過小,那么有限元算法的收斂性得不到保證。因為Delaunay三角剖分可以最大化最小角,所以Delaunay三角剖分是一種有效的剖分方法。

    目前,平面上的Delaunay三角剖分的研究已經(jīng)相當成熟。平面上的Delaunay三角剖分分為點集Delaunay三角剖分和約束Delaunay三角剖分。 對于平面點集Delaunay 三角剖分來說,Lee等[1]給出了點集Delaunay三角剖分的幾何性質(zhì)和算法。Dwyer[2]利用分治合并算法構(gòu)建了點集Delaunay三角剖分。近年來,對于點集Delaunay三角剖分的研究都是不斷降低時間復雜度,插入點混合定位算法[3]有效降低算法的時間復雜度。對于平面約束Delaunay三角剖分來說,Chew[4]提出了分治合并算法,Lee等[5]對平面直線圖形作約束Delaunay三角剖分。近幾年,關(guān)于平面約束Delaunay三角剖分的研究也都是不斷降低復雜度,如帶斷層約束的Delaunay三角剖分混合算法[6],它在生成網(wǎng)格質(zhì)量和時間效率上都優(yōu)于傳統(tǒng)算法。隨著計算機技術(shù)的發(fā)展,約束Delaunay三角剖分已經(jīng)不能滿足生產(chǎn)需求。在有限元分析中,網(wǎng)格的質(zhì)量對有限元方程的收斂性有著直接的影響。Ruppert[7]通過添加“Steiner”點提出Delaunay Refinement算法以及Shewchuk[8]提出的Delaunay Refinement算法增大了輸出三角形的最小角角度,這些網(wǎng)格都適用于有限元方程。但是Ruppert’s Refinement算法復雜度的證明由Barbic等[9]完成。Rand[10]提出的“Off-Centers”不僅提高輸出三角形的質(zhì)量,還添加更少的輸出點。

    近年來,出現(xiàn)了非歐空間的Delaunay三角剖分,如Caroli等[11]給出了球面點集Delaunay三角剖分和Zeng等[12]提出的雙曲曲面Delaunay Refinement三角剖分。在計算Ricci Flow方程時,一旦出現(xiàn)某個退化的球面三角形,就有可能出現(xiàn)不收斂的情況,文獻[13-14]研究就需要高質(zhì)量的球面三角形。因此,必須提高球面三角形的質(zhì)量。本文通過添加Delaunay劣弧中點以及球面三角形外接測地圓圓心,有效地提高整體球面三角形的質(zhì)量,避免狹小角的出現(xiàn),有效解決Ricci Flow方程不收斂的問題。

    1 球面凸類圖形Delaunay三角剖分再分算法

    首先,給出幾個相關(guān)概念。

    定義1 球面上兩點之間劣弧的長度,稱之為球面上兩點之間的距離, 簡稱距離。

    定義2 球面n邊形A1A2…An總在它的每一邊所在大圓的半個球面內(nèi), 那么稱這個球面多邊形為球面凸n邊形。

    定義3 球面凸n邊形以及包含在球面凸n邊形內(nèi)的點和劣弧稱為球面凸類圖形。

    因為不同半徑的球面之間可以建立一一映射,所以不妨在單位球面上對輸入的球面凸類圖形進行Delaunay三角剖分再分,如果沒有特殊說明,下文中的球面指的是單位球面。

    為了方便,假設(shè)球面凸類圖形中球面凸n邊形兩條相交劣弧夾角角度至少是90°,這個限制根據(jù)文獻[7]中的“Corner-Lopping” 可以解除,具體不再贅述;同時,記球面凸類圖形為X。另外本文中提到的角指的是球面角。

    球面凸類圖形Delaunay三角剖分再分算法是在球面點集Delaunay三角剖分算法基礎(chǔ)上生成的,球面點集Delaunay三角剖分(Sphere Delaunay Triangulation(V))簡記為SDT(V)。球面點集Delaunay三角剖分指的是給定球面上一系列點,由球面點集Delaunay三角剖分算法產(chǎn)生的所有球面三角形最短邊對應(yīng)的外接球面小圓圓心角達到最大,也就是最大化最小圓心角。如果球面上不存在四點共圓,那么這個SDT(V)是唯一的。如果存在四點共圓,兩條對角線任取一條不影響最小圓心角的大小,不妨取對角線最短的作為邊。具體的算法實現(xiàn)見參考文獻[11]。 為了方便下文的分析,本文只需要對球面凸類圖形內(nèi)部進行Delaunay三角剖分再分,故球面凸類圖形外的球面三角形全部刪除。

    把球面凸類圖形中的劣弧以及剖分產(chǎn)生的劣弧叫作Delaunay劣弧,用集合G表示。球面上任意一條劣弧AB記為AB。為了方便,也可用AB表示球面上A和B兩點之間的距離, 那么球面上AB的取值范圍為(0,π]。同時,把球面凸類圖形中的點以及剖分添加的點叫作頂點,用集合V表示。球面上任意一個點叫作點。對于這個算法,初始時,V就是所有輸入的頂點,G就是所有的輸入劣弧。

    生成球面點集Delaunay三角剖分之后,利用本文算法繼續(xù)剖分,添加點的原因有兩個:保證G中的Delaunay劣弧是點集Delaunay三角剖分產(chǎn)生的測地線;整體提高球面三角形最短邊外接球面小圓圓心角。

    對Delaunay劣弧g來說,以g為測地直徑的球面小圓稱為g的測地直徑圓。如果有一個頂點A在g的測地直徑圓內(nèi),稱頂點A“侵占”Delaunay劣弧g。如圖1所示, 輸入的球面凸類圖形是圖中的實Delaunay劣弧,球面點集Delaunay三角剖分是圖中除Delaunay劣弧AB外,其他實Delaunay劣弧和虛劣弧組成的圖形。對于球面凸類圖形來說,這不是它的Delaunay三角剖分,因為AB不是點集Delaunay三角剖分中的測地線且頂點A“侵占”Delaunay劣弧CD。

    本文算法的核心操作有兩個:第一,如果某條Delaunay劣弧被“侵占”,通過添加Delaunay劣弧中點分割Delaunay劣弧;第二,如果存在“瘦”球面三角形,通過添加球面三角形外接球面小圓圓心分解球面三角形,“瘦”球面三角形是球面三角形外接球面小圓中,球面三角形最短邊對應(yīng)的圓心角的角度小于φ(φ是預先給定的值)的球面三角形。

    圖1 非法球面Delaunay三角剖分Fig. 1 Delaunay triangulation of illegal sphere

    下面給出球面凸類圖形Delaunay三角剖分再分算法:

    輸入 球面凸類圖形X和參數(shù)φ。

    BEGIN

    初始化 計算初始球面點集Delaunay三角剖分SDT(V)

    while Delaunay劣弧g被“侵占”:

    SpliltGeo(Delaunay劣弧g)

    if “瘦”球面三角形t的外接球面小圓圓心P“侵占”Delaunay劣弧g1,g2,…,gk

    fori=1 tok:

    SplitGeo(Delaunay劣弧gi)

    else SplitStri(球面三角形t)

    end if

    until 沒有Delaunay劣弧被“侵占”和沒有球面三角形最短邊所對應(yīng)的圓心角<φ

    END

    子算法SplitStri(球面三角形t)。

    輸入 球面三角形t。

    BEGIN

    添加球面三角形t的外接球面小圓圓心到V中,更新SDT(V)

    END

    子算法SplitGeo(Delaunay劣弧g)。

    輸入 Delaunay劣弧g。

    BEGIN

    添加Delaunay劣弧g的中點到V中,更新SDT(V),將g從G中 移除,將g分成的兩個子段g1和g2添加到G中

    END

    2 算法收斂性分析

    在本章證明本文算法在執(zhí)行有限的步驟之后是可以終止的,這與輸入球面凸類圖形的局部特征尺度有關(guān),給出局部特征尺度的概念。

    定義4 給定一個球面凸類圖形X,以球面上任意一點P為圓心的與X中兩個和P不相關(guān)的頂點或Delaunay劣弧相交的最小球面小圓的測地半徑叫作點P的局部特征尺度,記為lfs(P)。

    在球面上,根據(jù)定義1知,lfs(P)的取值范圍為(0,π)。圖2是對這個定義的直觀解釋,lfs(A)、lfs(B)和lfs(C)分別對應(yīng)圖2中三個球面小圓的測地半徑,其中l(wèi)fs(C)有一部分頂點或Delaunay劣弧相交。對于給定的輸入球面凸類圖形X,lfs(X)是連續(xù)函數(shù)且滿足下列引理1。

    引理1 局部特征尺度是Lipschitz函數(shù),即任意給定一個球面凸類圖形,在球面上任取兩個點P和Q,有:

    lfs(Q)≤lfs(P)+PQ

    其中PQ是球面上兩點P和Q的距離。

    證明 根據(jù)lfs()的定義,lfs(P)是以點P為圓心的與球面凸類圖形中兩個和P不相關(guān)的頂點或Delaunay劣弧相交的最小球面小圓的測地半徑,記這個球面小圓為D1。

    1)若lfs(P)+PQ≥π,因為lfs()<π,顯然成立。

    2)若lfs(P)+PQ<π,現(xiàn)在以點Q為圓心,r=lfs(P)+PQ為測地半徑的球面小圓記為D2,那么D1包含在D2中,因此:

    lfs(Q)≤r

    故得到:

    lfs(Q)≤r=lfs(P)+PQ

    每當一個點添加之后,它不改變lfs()的性質(zhì),lfs()只和輸入的球面凸類圖形有關(guān),下面的引理2說明頂點的密度完全是由lfs()決定。

    引理2 1)初始時,對于每個輸入的頂點P,P到它最近的頂點的距離至少是lfs(P)。

    2)當點P是球面“瘦”球面三角形的外接球面小圓圓心時,點P到它最近頂點的距離至少是lfs(P)/CT(點P屬于這個三角剖分中的頂點或者是“侵占”某條Delaunay劣弧)。

    3)當頂點P是某條Delaunay劣弧的中點時,點P到它最近的頂點的距離至少為lfs(P)/CG。

    證明 當頂點P是輸入的點時,根據(jù)lfs()的定義,頂點P到它最近的頂點的距離至少是lfs(P)。對于后來加進去的點,假設(shè)這個引理對之前所有的頂點都成立。

    1)當點P是球面“瘦”三角形t的外接球面小圓圓心時,那么點P到球面“瘦”三角形t的頂點最近。設(shè)點P到球面“瘦”三角形t的頂點A、頂點B和頂點C的距離等于r1,那么r1的取值范圍為(0,π/2)。如圖3所示,不妨假設(shè)弧AB是球面△ABC的最短邊,那么弧AB所對應(yīng)的圓心角最小,設(shè)為ψ。

    不失一般性,首先考慮點A是在點B之后加進去的或者兩者都是輸入的頂點。當點A是輸入的頂點時,則lfs(A)≤AB;當點A是球面三角形外接球面小圓圓心時,設(shè)這個球面小圓的測地半徑為rA,根據(jù)這個引理有l(wèi)fs(A)≤CTrA≤CTAB;當點A是某條Delaunay劣弧的中點時,在A點利用這個引理有l(wèi)fs(A)≤CGAB。

    假設(shè)給定條件CG≥CT≥1,則上面幾種情況綜合起來為:

    lfs(A)≤CGAB

    如圖4所示,先把球面△PAB提出來,再由點P作AB的垂線,根據(jù)球面正弦定理有:

    由此得到:

    AB=2arcsin(sin(ψ/2)·sinr1)

    令ρ(ψ,r1)=(AB)/(2r1),因為:

    lfs(P)≤lfs(A)+r1

    則ρ(ψ,r1)在r1∈(0,π/2)上的最大值的上確界為sin(ψ/2),得:

    這里只需要取CT≥2sin(ψ/2)CG+1。

    圖3 球面三角形及其圓心角Fig. 3 Spherical triangle and its central angle

    圖4 球面三角形及其垂線Fig. 4 Spherical triangle and its vertical

    2)如圖5所示,當頂點P是Delaunay劣弧g的中點時,如果存在頂點A或某個球面“瘦”三角形外接球面小圓圓心A在以g為測地直徑的球面小圓內(nèi)時,Delaunay劣弧g需要添加中點并假設(shè)這個球面小圓的測地半徑為r2,易知,r2的取值范圍為(0,π/2)。

    圖5 點A“侵占”劣弧BCFig. 5 Vertex A “encroaches upon” minor arc BC

    a)點A位于Delaunay劣弧h上并且和Delaunay劣弧g不相交,之前假設(shè)所有的輸入球面凸類圖形的角度至少是90°,所以必然存在兩個不相交的Delaunay劣弧,一個包含點P,另一個包含點A,使得它們的距離比r2小。因此,得到:

    lfs(P)≤r2

    根據(jù)CG≥1知,這種情況是成立的。

    b)點A是球面“瘦”三角形外接球面小圓圓心,但是這個點在以Delaunay劣弧g為測地直徑的球面小圓內(nèi),所以這個點是不存在于三角剖分內(nèi)。那么存在一個以點A為圓心,測地半徑為rA的球面小圓CA使得圓內(nèi)沒有其他點,在A點應(yīng)用這個引理得到:

    rA≥lfs(A)/CT

    同時,點B和點C在CA之外,故當rA最大時有:

    cosrA=cosr2×cosr2

    令ρ(r2)=rA/r2,根據(jù)引理1有:

    lfs(P)≤lfs(A)+r2

    根據(jù)上面的分析,得到如下不等式:

    通過解這個方程組,得到:

    所以ψ的一個上界為φ=41.4°。

    引理2表明當點P加進去之后,沒有其他頂點在以P為圓心、lfs(P)/CG為測地半徑的球面小圓內(nèi)。下面的這個定理表明當點加進去之后,頂點P到其他頂點的距離存在下界。

    定理1 給定輸出的球面凸類圖形的Delaunay三角剖分再分,對于任意一個頂點P,頂點P到它最近的頂點的距離至少是lfs(P)/(CG+1)。

    證明 引理2解決了除P在Q之后加入的情況,利用這個引理2對Q進行處理,得到:

    PQ≥lfs(Q)/CG

    結(jié)合引理1得到:

    PQ≥lfs(P)/(CG+1)

    有了上面的定理之后,就能證明輸出頂點存在上界。

    定理2 Delaunay三角剖分再分輸出球面三角網(wǎng)格中,頂點的數(shù)量至多是:

    其中:B是球面凸類圖形的內(nèi)部;C1是具體的常數(shù)。

    根據(jù)引理1,在DP中,lfs()的最大值是lfs(P)+rP,設(shè)DP的面積為SDP,故:

    下面對rP分情況討論:

    1)當rP≤π/2時,因為SDP=4π sin2(rP/2),所以:

    綜合起來得:

    2)當π/2

    因此,只需取C1≥(2CG+3)2π即可。

    根據(jù)以上兩種情況知,只需取C1≥(2CG+3)2π。

    由上面的分析知,當輸入的參數(shù)φ≤41°時,本文算法在執(zhí)行過程中添加的點是有限的,也就是本文算法是收斂的。

    3 算法實現(xiàn)及分析

    球面凸類圖形的Delaunay三角剖分再分算法比歐氏空間的Delaunay Refinement Algorithm復雜。Delaunay Refinement Algorithm收斂條件由剖分的三角形最小角限定,本文算法收斂條件是由剖分三角形最短邊對應(yīng)的球面圓心角限定。這是因為球面三角形內(nèi)角和不是固定常數(shù),導致無法通過球面三角形最小角限定收斂條件。在具體算法實現(xiàn)過程中,球面凸類圖形Delaunay三角剖分再分算法比Delaunay Refinement Algorithm更復雜。Delaunay Refinement Algorithm可以由Matlab直接實現(xiàn)并且可視化,但是球面凸類圖形Delaunay三角剖分再分算法只能通過程序語言輸出邊和點的數(shù)據(jù),再由其他繪圖軟件可視化。

    本文算法的核心操作和平面的情況是相同的,具體的證明過程見參考文獻[9],所以本文算法的時間復雜度和平面情況相同,都是O(n2)。本文算法的程序?qū)崿F(xiàn)如圖6所示。

    圖6 程序?qū)崿F(xiàn)圖(最小圓心角為41.4°)Fig. 6 Program implementation diagram (minimum central angle is 41.4°)

    本文算法還可以適當提高最小圓心角的角度,可以在本文算法中插入一個限制球面三角形大小的語句,不妨以球面三角形外接測地圓半徑為限制。這種方法可以在球面上作近似均勻剖分,生成盡量同樣大小的三角形,并且可以加密剖分,滿足許多應(yīng)用的要求。

    4 結(jié)語

    近年來,有關(guān)Delaunay三角剖分研究重點仍然是在歐氏空間進行。對于非歐空間的Delaunay三角剖分研究得還不夠深入。在求解Ricci Flow方程時,Ricci Flow的收斂性是和剖分的三角形相關(guān)的,一旦出現(xiàn)某個角度過小的三角形,Ricci Flow就有可能出現(xiàn)不收斂情況。本文中的三角網(wǎng)格通過提高最短邊對應(yīng)的圓心角提高三角形最小角角度,有效解決不收斂的問題。本文算法中添加的是球面三角形外接測地圓的圓心,下一步可以研究添加“Off-center Steiner Vertices”,或許這種方法輸出的點更少,且可以突破本文收斂條件41.4°的限制,提高到更大的角度。

    References)

    [1] LEE D T, SCHACHTER B J. Two algorithms for constructing Delaunay triangulations [J]. International Journal of Computer and Information Sciences, 1980, 9(3): 219-242.

    [2] DWYER R A. A faster divide-and-conquer algorithm for constructing Delaunay triangulations [J]. Algorithmica, 1987, 2(1/2/3/4): 137-151.

    [3] 高莉.改進的Deluanay三角剖分算法研究[D].蘭州:蘭州交通大學,2015:30-33.(GAO L. An improved Deluanay triangulation algorithm [D]. Lanzhou: Lanzhou Jiaotong University, 2015: 30-33.)

    [4] CHEW L P. Constrained Delaunay triangulation [J]. Algorithmica, 1989, 4(1/2/3/4): 97-108.

    [5] LEE D T, LIN A K. Generalized Delaunay triangulation for planar graphs [J]. Discrete & Computational Geometry, 1986, 1(3): 201-217.

    [6] 張群會,解子毅.帶斷層約束的Delaunay三角剖分混合算法[J].西安科技大學學報,2014,34(1):52-56.(ZHANG Q H, XIE Z Y. Mixed algorithm of Delaunay triangular subdivision with fault constraint [J]. Journal of Xi’an University of Science and Technology, 2014, 34(1): 52-56.)

    [7] RUPPERT J. A Delaunay refinement algorithm for quality 2-dimensional mesh generation [J]. Journal of Algorithms, 1995, 18(3): 548-585.

    [8] SHEWCHUK J R. Delaunay refinement algorithms for triangular mesh generation [J]. Computational Geometry, 2002, 22(1/2/3): 21-74.

    [9] BARBIC J, MILLER G. A quadratic running time example for ruppert’s refinement algorithm [EB/OL]. [2017- 04- 16]. http://www-bcf.usc.edu/~jbarbic/BarbicMiller-RupperQuadraticExampleTechReport.pdf.

    [10] RAND A. Where and how Chew’s second Delaunay refinement algorithm works [C/OL]// Proceedings of the 2013 23rd Annual Canadian Conference on Computational Geometry, [2017- 04- 16]. http://www.w.cccg.ca/proceedings/2011/papers/paper91.pdf.

    [11] CAROLI M, DE CASTRO P M M, LORIOT S, et al. Robust and efficient Delaunay triangulations of points on or close to a sphere [C]// Proceedings of the 2010 9th International Conference on Experimental Algorithms. Berlin: Springer, 2010: 462-473.

    [12] ZENG W, SHI R, GU X F. Global surface remeshing using symmetric Delaunay triangulation in uniformization spaces [C]// Proceedings of the 2011 Eighth International Symposium on Voronoi Diagrams in Science and Engineering. Piscataway, NJ: IEEE, 2011: 160-169.

    [13] GU X F, GUO R, LUO F, et al. A discrete uniformization theorem for polyhedral surfaces II [EB/OL]. [2017- 04- 16]. http://math.oregonstate.edu/~guoren/docs/u10-10.pdf.

    [14] GU X F, GUO R, LUO F, et al. A discrete uniformization theorem for polyhedral surfaces I [EB/OL]. [2017- 04- 16]. https://arxiv.org/pdf/1401.4594.pdf.

    [15] PREPARATA F P, SHAMOS M I. Computational Geometry — an Introduction [M]. New York: Springer, 1985: 241-276.

    This work is partially supported by the Natural Science Foundation of Kunming University of Science and Technology (KKSY201507066).

    XIAJun, born in 1991, M. S. candidate. His research interests include computational geometry.

    LIYinghua, born in 1978, Ph. D., lecturer. His research interests include computational geometry.

    Delaunaytriangulationsubdivisionalgorithmofsphericalconvexgraphanditsconvergenceanalysis

    XIA Jun*, LI Yinghua

    (SchoolofScience,KunmingUniversityofScienceandTechnology,KunmingYunnan650500,China)

    When calculating curved Ricci Flow, non-convergence emerges due to the existence of undersized angles in triangular meshes. Concerning the problem of non-convergence, a Delaunay triangulation subdivision algorithm of spherical convex graph of enhancing the minimum angle was proposed. First of all, the Delaunay triangulation subdivision algorithm of spherical convex graph was given. The proposed algorithm had two key operations: 1) if a Delaunay minor arc was “encroached upon”, a midpoint of the Delaunay minor arc was added to segment the Delaunay minor arc; 2) if there was a “skinny” spherical triangle, it was disassembled by adding the center of minor circle of its circumscribed sphere. Then, the convergence criteria of the proposed algorithm was explored on local feature scale and an upper-bound formula of the output vertex was given. The grids based on the output of experiment show that the spherical triangle generated by the grids of the proposed algorithm has no narrow angle, so it is suitable for calculating Ricci Flow.

    sphere; Delaunay triangulation; minor arc; local feature scale; convergence

    2017- 06- 29;

    2017- 09- 04。

    昆明理工大學自然科學基金資助項目(KKSY201507066)。

    夏俊(1991—),男,安徽合肥人,碩士研究生,主要研究方向:計算幾何; 李映華(1978—),男,廣東梅州人,講師,博士,主要研究方向:計算幾何。

    1001- 9081(2017)12- 3558- 05

    10.11772/j.issn.1001- 9081.2017.12.3558

    (*通信作者電子郵箱yinghuali@kmust.edu.cn)

    TP301.6

    A

    猜你喜歡
    小圓圓心角剖分
    誰的周長長
    我高估了別人對我的關(guān)注度
    利用問題鏈導向深度學習的設(shè)計與實施
    各種各樣的扇形
    基于重心剖分的間斷有限體積元方法
    二元樣條函數(shù)空間的維數(shù)研究進展
    “圓心角”度數(shù):弧長計算的關(guān)鍵點
    一種實時的三角剖分算法
    復雜地電模型的非結(jié)構(gòu)多重網(wǎng)格剖分算法
    不同圓心角對小半徑曲線橋受力性能影響的研究
    国语自产精品视频在线第100页| 亚洲精品粉嫩美女一区| 亚洲av免费高清在线观看| 成人亚洲精品av一区二区| 国产伦在线观看视频一区| 中文字幕av成人在线电影| 亚洲国产欧洲综合997久久,| 国内精品宾馆在线| 亚洲欧美日韩高清在线视频| 国产精品一区二区性色av| 国产私拍福利视频在线观看| 亚洲最大成人手机在线| 国产aⅴ精品一区二区三区波| 国产毛片a区久久久久| 99国产极品粉嫩在线观看| 欧美精品国产亚洲| 亚洲狠狠婷婷综合久久图片| 少妇的逼好多水| 韩国av一区二区三区四区| 午夜免费男女啪啪视频观看 | 亚洲成a人片在线一区二区| 亚洲成a人片在线一区二区| 干丝袜人妻中文字幕| 色av中文字幕| 91精品国产九色| 成人二区视频| 久久久久精品国产欧美久久久| 十八禁国产超污无遮挡网站| 成人一区二区视频在线观看| 欧美高清成人免费视频www| 日日夜夜操网爽| av在线亚洲专区| av天堂中文字幕网| 成年女人看的毛片在线观看| 97碰自拍视频| 久久精品国产亚洲av香蕉五月| 国产国拍精品亚洲av在线观看| 啪啪无遮挡十八禁网站| 别揉我奶头 嗯啊视频| 亚洲国产精品合色在线| 一级a爱片免费观看的视频| 欧美一区二区亚洲| 欧美中文日本在线观看视频| h日本视频在线播放| avwww免费| 91在线精品国自产拍蜜月| 免费看日本二区| 日韩欧美精品免费久久| 性欧美人与动物交配| av在线亚洲专区| 久久精品人妻少妇| 久久久精品欧美日韩精品| x7x7x7水蜜桃| 成人高潮视频无遮挡免费网站| 国产午夜精品久久久久久一区二区三区 | 日本黄色视频三级网站网址| 久久久久久久久大av| 国产精品久久视频播放| 亚洲欧美日韩高清在线视频| 亚洲av中文av极速乱 | 国产亚洲精品av在线| 国产精品久久久久久亚洲av鲁大| 观看美女的网站| 啦啦啦啦在线视频资源| 白带黄色成豆腐渣| 丰满的人妻完整版| 1024手机看黄色片| 波多野结衣巨乳人妻| 97碰自拍视频| 亚洲av熟女| 99热这里只有精品一区| 国产一区二区三区视频了| 亚洲国产欧洲综合997久久,| 午夜精品久久久久久毛片777| 国产精品电影一区二区三区| 婷婷精品国产亚洲av在线| 欧美性猛交╳xxx乱大交人| 99热6这里只有精品| 联通29元200g的流量卡| 网址你懂的国产日韩在线| 免费人成在线观看视频色| 级片在线观看| 中国美白少妇内射xxxbb| av专区在线播放| 国产精品亚洲一级av第二区| av国产免费在线观看| 色综合站精品国产| 一本一本综合久久| 女人被狂操c到高潮| 人人妻人人看人人澡| 老熟妇乱子伦视频在线观看| 国产av在哪里看| 欧美精品国产亚洲| 亚洲性夜色夜夜综合| 在线免费十八禁| 国产 一区 欧美 日韩| 白带黄色成豆腐渣| 一区福利在线观看| 欧美日韩乱码在线| 在线看三级毛片| 久久久久久久亚洲中文字幕| 人妻制服诱惑在线中文字幕| 亚洲欧美日韩无卡精品| 欧美日韩中文字幕国产精品一区二区三区| 给我免费播放毛片高清在线观看| 直男gayav资源| 少妇熟女aⅴ在线视频| 亚洲精品亚洲一区二区| 亚洲一级一片aⅴ在线观看| 久久6这里有精品| 亚洲精华国产精华精| 99在线视频只有这里精品首页| 国产日本99.免费观看| 精品人妻一区二区三区麻豆 | 能在线免费观看的黄片| 亚洲久久久久久中文字幕| 色综合婷婷激情| 99视频精品全部免费 在线| 国产真实伦视频高清在线观看 | 国产色爽女视频免费观看| 搡老熟女国产l中国老女人| 亚洲一区二区三区色噜噜| 3wmmmm亚洲av在线观看| 琪琪午夜伦伦电影理论片6080| 中国美女看黄片| 亚洲四区av| 日日撸夜夜添| 琪琪午夜伦伦电影理论片6080| 听说在线观看完整版免费高清| 狠狠狠狠99中文字幕| 国产精华一区二区三区| 一级黄片播放器| 亚洲第一电影网av| 欧美高清成人免费视频www| 全区人妻精品视频| 精品人妻视频免费看| 我的女老师完整版在线观看| 免费观看人在逋| 亚洲欧美日韩高清专用| 日韩欧美国产在线观看| 欧美激情国产日韩精品一区| 国产精品综合久久久久久久免费| 九九久久精品国产亚洲av麻豆| 免费看美女性在线毛片视频| x7x7x7水蜜桃| 午夜免费男女啪啪视频观看 | 国产 一区 欧美 日韩| 啦啦啦啦在线视频资源| 中文资源天堂在线| 村上凉子中文字幕在线| 淫秽高清视频在线观看| 欧美成人一区二区免费高清观看| 在线播放国产精品三级| 亚洲专区国产一区二区| 午夜福利欧美成人| 亚洲av第一区精品v没综合| 亚洲精品456在线播放app | 久久精品综合一区二区三区| 国产精华一区二区三区| 国产大屁股一区二区在线视频| 看十八女毛片水多多多| 国产乱人视频| 亚洲 国产 在线| 色哟哟·www| 别揉我奶头~嗯~啊~动态视频| 日本三级黄在线观看| 日日摸夜夜添夜夜添av毛片 | 国内精品宾馆在线| 亚洲天堂国产精品一区在线| 韩国av一区二区三区四区| 国产精品久久久久久久电影| 精品久久久久久久久亚洲 | 看黄色毛片网站| 国产视频一区二区在线看| 一区福利在线观看| 免费大片18禁| 欧美xxxx性猛交bbbb| 99热精品在线国产| 91久久精品电影网| 丰满的人妻完整版| а√天堂www在线а√下载| 久久99热6这里只有精品| 99精品在免费线老司机午夜| 身体一侧抽搐| 九九热线精品视视频播放| 老司机福利观看| 国产精品久久久久久久电影| 一区二区三区四区激情视频 | 久久久久久久久久久丰满 | 日本一二三区视频观看| 天天一区二区日本电影三级| 国产精品日韩av在线免费观看| 国产精品伦人一区二区| 成人毛片a级毛片在线播放| 深爱激情五月婷婷| 欧美极品一区二区三区四区| 22中文网久久字幕| 91久久精品国产一区二区成人| 国产高清三级在线| 一区福利在线观看| 中文字幕精品亚洲无线码一区| 少妇熟女aⅴ在线视频| 国产精品免费一区二区三区在线| 啦啦啦啦在线视频资源| 国产精品一区二区性色av| 噜噜噜噜噜久久久久久91| 别揉我奶头 嗯啊视频| 在线观看av片永久免费下载| 高清毛片免费观看视频网站| 特级一级黄色大片| 九色成人免费人妻av| 老司机福利观看| 亚洲中文字幕一区二区三区有码在线看| 听说在线观看完整版免费高清| 精品久久久久久久久久久久久| 精品人妻偷拍中文字幕| av视频在线观看入口| 十八禁国产超污无遮挡网站| 一区福利在线观看| 18禁黄网站禁片午夜丰满| 欧美日韩瑟瑟在线播放| 欧美日韩中文字幕国产精品一区二区三区| 久久精品国产鲁丝片午夜精品 | 中出人妻视频一区二区| av.在线天堂| 最近视频中文字幕2019在线8| 色尼玛亚洲综合影院| 午夜福利高清视频| videossex国产| 日韩在线高清观看一区二区三区 | 国产日本99.免费观看| av福利片在线观看| 日韩欧美免费精品| 久久久精品大字幕| 国产 一区精品| 久久国内精品自在自线图片| 亚洲国产精品合色在线| 变态另类丝袜制服| 精品久久久久久久久久久久久| 午夜免费成人在线视频| 午夜精品在线福利| 亚洲专区国产一区二区| 在线观看美女被高潮喷水网站| 国产在线精品亚洲第一网站| 久久这里只有精品中国| 国产69精品久久久久777片| 欧美三级亚洲精品| 日本 欧美在线| 免费av观看视频| www.色视频.com| 很黄的视频免费| 中文亚洲av片在线观看爽| 超碰av人人做人人爽久久| 精品久久久久久久久久免费视频| 少妇被粗大猛烈的视频| 午夜激情欧美在线| 免费在线观看影片大全网站| 欧美zozozo另类| 久久久久久久久中文| 色吧在线观看| 成人鲁丝片一二三区免费| 精品一区二区三区视频在线观看免费| 免费搜索国产男女视频| 老司机深夜福利视频在线观看| 88av欧美| 亚洲天堂国产精品一区在线| 国产主播在线观看一区二区| 国产精品久久久久久久久免| 色在线成人网| 好男人在线观看高清免费视频| 亚洲色图av天堂| 欧美性猛交╳xxx乱大交人| 观看免费一级毛片| 久久久成人免费电影| 久久久久性生活片| 成人毛片a级毛片在线播放| 中亚洲国语对白在线视频| 尾随美女入室| 国产男人的电影天堂91| 免费黄网站久久成人精品| 久久欧美精品欧美久久欧美| 99在线视频只有这里精品首页| 国产精品98久久久久久宅男小说| 99九九线精品视频在线观看视频| 天堂网av新在线| 国产高清激情床上av| 国产精品美女特级片免费视频播放器| 熟女电影av网| 99久久无色码亚洲精品果冻| 日韩,欧美,国产一区二区三区 | 人妻制服诱惑在线中文字幕| 极品教师在线免费播放| 日韩一区二区视频免费看| 国产高清激情床上av| x7x7x7水蜜桃| 一区二区三区高清视频在线| 真人一进一出gif抽搐免费| 国内精品一区二区在线观看| 99九九线精品视频在线观看视频| 赤兔流量卡办理| 亚州av有码| 校园春色视频在线观看| 99国产极品粉嫩在线观看| 久久久色成人| 亚洲一区二区三区色噜噜| 狠狠狠狠99中文字幕| 尾随美女入室| 最近视频中文字幕2019在线8| 美女被艹到高潮喷水动态| 精品乱码久久久久久99久播| 国产aⅴ精品一区二区三区波| 男女下面进入的视频免费午夜| 久久99热6这里只有精品| 老熟妇仑乱视频hdxx| 久久久精品大字幕| 88av欧美| 亚洲在线观看片| 国产伦精品一区二区三区视频9| 亚洲精品456在线播放app | 中出人妻视频一区二区| 国产精品久久久久久亚洲av鲁大| bbb黄色大片| 不卡一级毛片| 久久久色成人| 亚洲中文字幕一区二区三区有码在线看| 一区二区三区免费毛片| 国内精品一区二区在线观看| 欧美zozozo另类| 国产伦在线观看视频一区| 国产精品爽爽va在线观看网站| 免费观看精品视频网站| 国产成年人精品一区二区| 天堂影院成人在线观看| 男人舔女人下体高潮全视频| 欧美日韩亚洲国产一区二区在线观看| av国产免费在线观看| 美女高潮的动态| 波多野结衣高清无吗| 1024手机看黄色片| 老熟妇仑乱视频hdxx| 午夜精品在线福利| 麻豆久久精品国产亚洲av| 午夜免费激情av| 亚洲经典国产精华液单| 亚洲aⅴ乱码一区二区在线播放| 色5月婷婷丁香| 久久久久久大精品| 欧美日本视频| 国产精品精品国产色婷婷| 少妇的逼水好多| 国产国拍精品亚洲av在线观看| 久久久精品大字幕| 99热只有精品国产| 国产色爽女视频免费观看| 日韩欧美一区二区三区在线观看| 亚洲图色成人| 99热网站在线观看| 国产 一区 欧美 日韩| 99久久中文字幕三级久久日本| 亚洲不卡免费看| avwww免费| 亚洲欧美精品综合久久99| 深夜a级毛片| 国产亚洲欧美98| 人人妻,人人澡人人爽秒播| 午夜激情欧美在线| 久久人人精品亚洲av| 欧美色欧美亚洲另类二区| 日韩欧美国产一区二区入口| 亚洲人成伊人成综合网2020| 人人妻,人人澡人人爽秒播| 51国产日韩欧美| 精品人妻偷拍中文字幕| 可以在线观看的亚洲视频| 日本免费一区二区三区高清不卡| 好男人在线观看高清免费视频| 永久网站在线| 能在线免费观看的黄片| 成人国产一区最新在线观看| 国产麻豆成人av免费视频| 麻豆久久精品国产亚洲av| 久久久色成人| 嫩草影视91久久| 少妇人妻精品综合一区二区 | 欧美潮喷喷水| 久久精品国产清高在天天线| 97超级碰碰碰精品色视频在线观看| 校园春色视频在线观看| 欧美最黄视频在线播放免费| 深夜a级毛片| 男人舔奶头视频| 俄罗斯特黄特色一大片| 真人做人爱边吃奶动态| 亚洲18禁久久av| 日韩欧美在线二视频| 色在线成人网| 亚洲欧美清纯卡通| 又黄又爽又免费观看的视频| 我要看日韩黄色一级片| 天堂网av新在线| 亚洲成人中文字幕在线播放| 日本欧美国产在线视频| 黄色一级大片看看| 长腿黑丝高跟| 久久久久精品国产欧美久久久| 色在线成人网| 春色校园在线视频观看| 免费看光身美女| 国产黄a三级三级三级人| 精品午夜福利视频在线观看一区| 国产免费一级a男人的天堂| 69av精品久久久久久| 小蜜桃在线观看免费完整版高清| 久久久久久伊人网av| 久久99热6这里只有精品| 一区二区三区激情视频| 观看美女的网站| 国产在线男女| 日韩欧美三级三区| 午夜福利成人在线免费观看| 亚洲精品亚洲一区二区| 18禁裸乳无遮挡免费网站照片| 亚洲va日本ⅴa欧美va伊人久久| 国产精品日韩av在线免费观看| 女的被弄到高潮叫床怎么办 | 99国产精品一区二区蜜桃av| 国产亚洲精品av在线| 亚洲av成人av| 九色国产91popny在线| 国产精品人妻久久久影院| 欧美潮喷喷水| 免费观看人在逋| 国内精品一区二区在线观看| 18禁裸乳无遮挡免费网站照片| 国内精品美女久久久久久| 一进一出抽搐动态| 国产三级中文精品| 亚洲精品456在线播放app | 真实男女啪啪啪动态图| 色哟哟·www| 亚洲国产日韩欧美精品在线观看| 嫩草影院新地址| 婷婷精品国产亚洲av| 性色avwww在线观看| 1000部很黄的大片| 成年免费大片在线观看| 亚洲精品成人久久久久久| 久久精品国产亚洲av天美| 国产私拍福利视频在线观看| 日日摸夜夜添夜夜添小说| 99九九线精品视频在线观看视频| 可以在线观看的亚洲视频| 一区二区三区高清视频在线| 国产69精品久久久久777片| 婷婷色综合大香蕉| 男女那种视频在线观看| 有码 亚洲区| 国产午夜精品论理片| 岛国在线免费视频观看| 99视频精品全部免费 在线| 日韩亚洲欧美综合| 2021天堂中文幕一二区在线观| 欧美极品一区二区三区四区| 88av欧美| 午夜免费男女啪啪视频观看 | 在线免费观看不下载黄p国产 | 成熟少妇高潮喷水视频| 国产91精品成人一区二区三区| 联通29元200g的流量卡| 欧美性猛交╳xxx乱大交人| 国产免费男女视频| 欧美成人一区二区免费高清观看| 身体一侧抽搐| 国产69精品久久久久777片| 日韩精品有码人妻一区| 国内精品久久久久精免费| 亚洲狠狠婷婷综合久久图片| 网址你懂的国产日韩在线| 性色avwww在线观看| 成人毛片a级毛片在线播放| 日本黄色片子视频| www.色视频.com| av专区在线播放| 在线观看美女被高潮喷水网站| 婷婷精品国产亚洲av| 在现免费观看毛片| 国产午夜精品久久久久久一区二区三区 | 亚洲四区av| 午夜福利在线在线| 蜜桃久久精品国产亚洲av| 亚洲av成人av| 亚洲性久久影院| 一进一出抽搐gif免费好疼| 亚洲一区高清亚洲精品| 国产精品久久久久久精品电影| 成人一区二区视频在线观看| 亚洲av一区综合| 亚洲专区国产一区二区| 一进一出抽搐动态| 日本 av在线| 欧美最黄视频在线播放免费| 91在线精品国自产拍蜜月| 亚洲国产高清在线一区二区三| 国内精品久久久久久久电影| 久久亚洲精品不卡| 国产激情偷乱视频一区二区| 国产黄色小视频在线观看| 国产精品久久久久久久电影| 99热6这里只有精品| 精品国产三级普通话版| 国产欧美日韩精品一区二区| av视频在线观看入口| 午夜福利在线观看免费完整高清在 | 国产一区二区激情短视频| 久久精品国产亚洲av涩爱 | 国产一区二区激情短视频| 日韩国内少妇激情av| 给我免费播放毛片高清在线观看| 亚洲美女黄片视频| 十八禁国产超污无遮挡网站| 国产精品女同一区二区软件 | 美女 人体艺术 gogo| 亚洲精品日韩av片在线观看| 国产高清三级在线| 日韩欧美免费精品| 久久精品91蜜桃| 97超级碰碰碰精品色视频在线观看| 男女做爰动态图高潮gif福利片| 日日摸夜夜添夜夜添小说| 日韩av在线大香蕉| 国产乱人视频| 日韩欧美在线二视频| 亚洲国产欧美人成| 国产乱人伦免费视频| 男插女下体视频免费在线播放| 一区二区三区激情视频| 美女cb高潮喷水在线观看| 蜜桃久久精品国产亚洲av| 大型黄色视频在线免费观看| 在线看三级毛片| 波多野结衣高清无吗| 久久99热6这里只有精品| 2021天堂中文幕一二区在线观| 两人在一起打扑克的视频| 国产精品福利在线免费观看| 干丝袜人妻中文字幕| 免费看a级黄色片| 国产单亲对白刺激| 色噜噜av男人的天堂激情| 免费黄网站久久成人精品| 丰满的人妻完整版| 国产色爽女视频免费观看| 十八禁网站免费在线| a级一级毛片免费在线观看| 淫秽高清视频在线观看| 尾随美女入室| 校园人妻丝袜中文字幕| 成人高潮视频无遮挡免费网站| 欧美成人a在线观看| 亚洲成人免费电影在线观看| 欧美色视频一区免费| 日本免费一区二区三区高清不卡| 非洲黑人性xxxx精品又粗又长| av天堂在线播放| 又爽又黄无遮挡网站| 男人舔女人下体高潮全视频| 色尼玛亚洲综合影院| 丰满人妻一区二区三区视频av| av在线老鸭窝| 变态另类丝袜制服| 久久久久久久久中文| 亚洲色图av天堂| 91久久精品国产一区二区成人| av天堂中文字幕网| 丝袜美腿在线中文| 国语自产精品视频在线第100页| 在线观看午夜福利视频| 身体一侧抽搐| 成人av在线播放网站| av天堂在线播放| 中出人妻视频一区二区| 国产精品av视频在线免费观看| 日韩欧美国产在线观看| АⅤ资源中文在线天堂| 中国美白少妇内射xxxbb| 一区二区三区四区激情视频 | 91久久精品国产一区二区成人| 欧美日韩综合久久久久久 | 中文字幕久久专区| 久久久久久大精品| 两性午夜刺激爽爽歪歪视频在线观看| 免费搜索国产男女视频| 亚洲中文字幕一区二区三区有码在线看| 精品人妻熟女av久视频| 亚洲美女黄片视频| 高清在线国产一区| 亚洲内射少妇av| 亚洲国产日韩欧美精品在线观看| 欧美激情在线99| 全区人妻精品视频| 欧美3d第一页| 女生性感内裤真人,穿戴方法视频| 色综合色国产| 久久99热6这里只有精品| 国产高清视频在线观看网站| 国产黄片美女视频| 两人在一起打扑克的视频| 久久精品影院6| 国产精品98久久久久久宅男小说| 久久99热6这里只有精品| 精品午夜福利在线看| 国产黄片美女视频| 色视频www国产| 九色成人免费人妻av| 国产老妇女一区|