《電子技術(shù)應(yīng)用》
您所在的位置:首頁(yè) > 可編程邏輯 > 設(shè)計(jì)應(yīng)用 > 基于FPGA的像素探測(cè)器數(shù)據(jù)緩存設(shè)計(jì)
基于FPGA的像素探測(cè)器數(shù)據(jù)緩存設(shè)計(jì)
來(lái)源:電子技術(shù)應(yīng)用2011年第9期
楊 萌, 王祖強(qiáng)
(山東大學(xué) 信息科學(xué)與工程學(xué)院, 山東 濟(jì)南104221)
摘要: 針對(duì)粒子對(duì)撞機(jī)像素探測(cè)器的數(shù)據(jù)緩存,提出了基于FPGA的解決方法。通過(guò)時(shí)序轉(zhuǎn)換匹配模塊,可以使用定制FIFO核實(shí)現(xiàn)對(duì)探測(cè)器模型的數(shù)據(jù)進(jìn)行緩存。經(jīng)過(guò)FPGA驗(yàn)證,在功能上達(dá)到了像素探測(cè)器模型的要求,對(duì)于粒子對(duì)撞機(jī)像素探測(cè)器的低成本工程驗(yàn)證具有顯著的參考價(jià)值。
關(guān)鍵詞: FPGA FIFO 像素探測(cè)器 CycloneII
中圖分類號(hào): TP333
文獻(xiàn)標(biāo)識(shí)碼: B
文章編號(hào): 0258-7998(2011)09-092-05
Design and realization of pixel detector data buffer based on FPGA
Yang Meng, Wang Zuqiang
Information Science and Technology Institute, Shangdong University, Jinan 104221, China
Abstract: In order to solving data buffer problem of partical collider pixel detector model, this paper introduces a method based on FPGA device. By using time sequence switch module, it could use custom-built FIFO IP core for buffering data from pixel detector output. After testing on FPGA development board, this method could satisfy the requirement of pixel detector model, and it has remarkable reference value to experimentation of partical collider pixel detector project.
Key words : FIFO; FPGA; pixel detector


 粒子對(duì)撞機(jī)[1]是建立在高能同步加速器上的一種粒子對(duì)撞裝置。對(duì)撞機(jī)運(yùn)行時(shí),能量強(qiáng)大的粒子在其內(nèi)部經(jīng)過(guò)加速后對(duì)撞,以達(dá)到一定的相互反應(yīng)速率。對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行分析研究,可以幫助實(shí)驗(yàn)者了解相應(yīng)粒子的形態(tài)等特性,從而推動(dòng)各種理論或技術(shù)的發(fā)展。
 像素探測(cè)器[2]是粒子對(duì)撞機(jī)內(nèi)部用于探測(cè)質(zhì)子撞擊后所形成的離子軌跡的裝置,是一種專用傳感器,由像素探測(cè)器探測(cè)得出的數(shù)據(jù)結(jié)果即可繪制出撞擊后的離子軌跡。目前,像素探測(cè)器主要采用以下幾種設(shè)計(jì)方法:ASIC設(shè)計(jì)方法[3],其優(yōu)點(diǎn)是集成度高且功耗低,但是設(shè)計(jì)成本高、周期長(zhǎng)且風(fēng)險(xiǎn)較大;分立元件采用微處理器的設(shè)計(jì)方法,其優(yōu)點(diǎn)是成本低、設(shè)計(jì)周期短且風(fēng)險(xiǎn)小,但是速度慢且體積大。針對(duì)這幾種設(shè)計(jì)方法,本文結(jié)合FPGA[4]提出了一種基于粒子對(duì)撞機(jī)像素探測(cè)器的數(shù)據(jù)緩存方法設(shè)計(jì)。相對(duì)ASIC設(shè)計(jì)方法,具有設(shè)計(jì)投入低、參數(shù)設(shè)置靈活、開(kāi)發(fā)風(fēng)險(xiǎn)低而且開(kāi)發(fā)周期短的優(yōu)勢(shì);相對(duì)于分立元件采用微處理器的設(shè)計(jì)方法,具有體積小、功耗低且速度快的優(yōu)勢(shì)。
 由于像素探測(cè)器的數(shù)據(jù)輸出具有保持時(shí)間短(最少保持時(shí)間為25 ns)、數(shù)據(jù)間隔時(shí)間不確定等特征,所以對(duì)于數(shù)據(jù)的存儲(chǔ)必須經(jīng)過(guò)時(shí)序轉(zhuǎn)換電路將數(shù)據(jù)事先存入緩存,再寫(xiě)入存儲(chǔ)器。本文的數(shù)據(jù)緩存設(shè)計(jì)由FIFO輸入時(shí)序模塊、FIFO模塊和FIFO輸出時(shí)序模塊三部分組成,通過(guò)驗(yàn)證可以實(shí)現(xiàn)由探測(cè)器數(shù)據(jù)輸出端到存儲(chǔ)器的數(shù)據(jù)緩存功能。
1 總體設(shè)計(jì)
    FPGA是一種新型的可編程邏輯器件,可以通過(guò)硬件描述語(yǔ)言生成所需的硬件,相對(duì)于傳統(tǒng)的ASIC設(shè)計(jì),具有設(shè)計(jì)周期短、靈活性高、開(kāi)發(fā)費(fèi)用少、設(shè)計(jì)風(fēng)險(xiǎn)低等優(yōu)點(diǎn),是ASIC的一種高效低成本解決方案。本設(shè)計(jì)使用的FPGA是Altera公司的CYCLONE2芯片,最高時(shí)鐘頻率100 MHz。總體架構(gòu)如圖1所示。

 由于探測(cè)器數(shù)據(jù)輸出需要符合FIFO數(shù)據(jù)輸入的時(shí)序,F(xiàn)IFO的數(shù)據(jù)輸出需要符合存儲(chǔ)器輸入的時(shí)序,因此設(shè)計(jì)了FIFO輸入時(shí)序模塊和FIFO輸出時(shí)序模塊,使數(shù)據(jù)能夠正確緩存進(jìn)而存儲(chǔ)。下面對(duì)具體的時(shí)序規(guī)范及各個(gè)模塊進(jìn)行說(shuō)明。

 


2 模塊設(shè)計(jì)及說(shuō)明
2.1 FIFO模塊

 FIFO即先進(jìn)先出緩存器。FIFO與其他存儲(chǔ)器的最大區(qū)別就是沒(méi)有外部讀寫(xiě)的地址線,因此應(yīng)用起來(lái)十分方便簡(jiǎn)單,但是相應(yīng)的缺點(diǎn)就是只能順序?qū)懭牖蛘唔樞蜃x取,其尋址方式由內(nèi)部指針自加自減完成,不能對(duì)特定的地址進(jìn)行讀寫(xiě)。FIFO通常使用在不同速度的接口數(shù)據(jù)交換中,通過(guò)FIFO使得時(shí)序分別符合兩個(gè)接口的特定速度,以進(jìn)行數(shù)據(jù)傳輸。
 Altera的FIFO按照驅(qū)動(dòng)時(shí)鐘分類有兩種:?jiǎn)螘r(shí)鐘FIFO(SCFIFO),即讀和寫(xiě)FIFO用的同一個(gè)時(shí)鐘信號(hào);雙時(shí)鐘FIFO(DCFIFO),即讀和寫(xiě)FIFO分別使用不同時(shí)鐘信號(hào)。其中雙時(shí)鐘FIFO還根據(jù)不同的端口數(shù)據(jù)寬度分為DCFIFO和DCFIFO_MIXED_WIDTH,所謂DCFIFO是輸入數(shù)據(jù)和輸出數(shù)據(jù)具有相同的數(shù)據(jù)寬度;而DCFIFO_MIXED_WIDTH是指輸入輸出端口可以使用不同的數(shù)據(jù)寬度。
 使用Altera的FIFO MegaWizard inteface launched對(duì)FIFO進(jìn)行構(gòu)建,像素探測(cè)器模型中有16路傳感器輸出,每路輸出需要記錄100個(gè)數(shù)據(jù)。根據(jù)模型需要,構(gòu)建FIFO的主要參數(shù)如下:
    dcfifo_component.intended_device_family
    = "Cyclone II"dcfifo_component.lpm_hint"MAXIMIZE_
        SPEED=7  RAM_BLOCK_TYPE=M4K"
    dcfifo_component.lpm_numwords = 128
    dcfifo_component.lpm_showahead = "OFF"
    dcfifo_component.lpm_type = "dcfifo"
    dcfifo_component.lpm_width = 16
    dcfifo_component.lpm_widthu = 7
    dcfifo_component.overflow_checking="ON" dcfifo_compo-
        nent.rdsync_delaypipe = 5
    dcfifo_component.underflow_checking="ON"
    dcfifo_component.use_eab = "ON"
    dcfifo_component.write_aclr_synch = "OFF"
    dcfifo_component.wrsync_delaypipe = 5;
 構(gòu)建完FIFO后,對(duì)其進(jìn)行時(shí)序仿真,仿真結(jié)果如圖2所示。

 參照時(shí)序仿真圖對(duì)各個(gè)端口及時(shí)序進(jìn)行說(shuō)明:
 aclr: 異步清零端,1 bit。清零所有輸出狀態(tài)端口,對(duì)于DCFIFO,3個(gè)wrclk時(shí)鐘上升沿后清零wrfull端口,清零rdfull端口。如果輸出端口定義reg類型的則會(huì)被清零;否則會(huì)保持輸出值。
 data: 數(shù)據(jù)輸入端口, 16 bit。當(dāng)寫(xiě)請(qǐng)求wrreq有效時(shí),保持?jǐn)?shù)據(jù)直到數(shù)據(jù)被寫(xiě)入FIFO。當(dāng)使用手動(dòng)定義FIFO時(shí),其數(shù)據(jù)寬度用參數(shù)LPM_WIDTH定義。
 q: 數(shù)據(jù)輸出端口, 16 bit。當(dāng)有數(shù)據(jù)請(qǐng)求時(shí)(rdreq有效時(shí)),輸出數(shù)據(jù)。對(duì)于DCFIFO,輸出數(shù)據(jù)的寬度可以與輸入數(shù)據(jù)data端口寬度不同,具體用參數(shù)LPM_WIDTH_R定義。
 rdcld: 上升沿出發(fā)時(shí)鐘,1bit。用來(lái)同步以下信號(hào):q、dreq、dfull、dempty、rdusedw。
 rdempty: 輸出數(shù)據(jù)為零時(shí)輸出高電平,1 bit。不管目標(biāo)設(shè)備是什么,在讀請(qǐng)求發(fā)送前必須查詢r(jià)dempty信號(hào)是否為高電平,以避免錯(cuò)誤指令發(fā)出。
 rdreq: 讀請(qǐng)求信號(hào)端口,1 bit。當(dāng)需要從FIFO中讀數(shù)據(jù)時(shí),向rdreq端發(fā)送讀請(qǐng)求,讀取數(shù)據(jù)個(gè)數(shù)與rdreq持續(xù)的時(shí)鐘數(shù)相同。需要注意的是,當(dāng)rdempty有效時(shí)不能發(fā)送rdreq信號(hào)。對(duì)于這種情況,可以開(kāi)啟空保護(hù)功能,通過(guò)設(shè)置參數(shù)UNDERFLOW_CHECKING高電平實(shí)現(xiàn),當(dāng)rdempty為高電平時(shí)rdreq信號(hào)自動(dòng)被置為無(wú)效。
 rdusedw: 輸出數(shù)據(jù)顯示FIFO中可讀數(shù)據(jù)的數(shù)據(jù)量,7 bit。在DCFIFO中,其端口寬度要與手動(dòng)設(shè)置參數(shù)LPM_WIDTHU相等。需要注意的是,對(duì)于Cyclone系列的FPGA,當(dāng)顯示數(shù)據(jù)滿時(shí)實(shí)際FIFO的存儲(chǔ)數(shù)據(jù)量有可能并沒(méi)有達(dá)到存儲(chǔ)的最大值,因此應(yīng)該參考full或者wrfull端口來(lái)執(zhí)行正確的寫(xiě)操作,參考empty或者rdempty端口來(lái)執(zhí)行正確的讀操作。
 wrclk: 寫(xiě)數(shù)據(jù)時(shí)鐘,上升沿觸發(fā)有效,1 bit。用于同步以下端口:data、wrreq、wrusedw、 wrfull和wrempty。
 wrfull: 寫(xiě)數(shù)據(jù)滿信號(hào),1 bit。當(dāng)此端口電平有效時(shí),F(xiàn)IFO已經(jīng)被寫(xiě)滿。其注意事項(xiàng)同rdfull。總體而言,rdfull信號(hào)要比wrfull信號(hào)有所延遲,因此,應(yīng)該通過(guò)wrfull信號(hào)的電平來(lái)判斷是否可以發(fā)送寫(xiě)請(qǐng)求信號(hào)wrreq。
 wrreq: 寫(xiě)請(qǐng)求信號(hào)。當(dāng)需要向FIFO寫(xiě)入數(shù)據(jù)時(shí),向wrreq端發(fā)送讀請(qǐng)求,讀取數(shù)據(jù)個(gè)數(shù)與wrreq持續(xù)的時(shí)鐘數(shù)相同。需要注意的是,當(dāng)wrfull有效時(shí)不能發(fā)送wrreq信號(hào)。對(duì)于這種情況,可以開(kāi)啟溢出保護(hù)功能,通過(guò)設(shè)置參數(shù)OVERFLOW_CHECKING高電平實(shí)現(xiàn),當(dāng)wrfull為高電平時(shí)wrreq信號(hào)自動(dòng)被置為無(wú)效。同時(shí)在取消aclr信號(hào)時(shí)不應(yīng)輸入wrreq信號(hào),否則aclr信號(hào)的下降沿和wrreq信號(hào)置為高電平后寫(xiě)數(shù)據(jù)的上升沿會(huì)產(chǎn)生競(jìng)爭(zhēng)冒險(xiǎn)現(xiàn)象。對(duì)于CYCLONE系列的FPGA中DCFIFO器件,可以選擇添加同步電路同步aclr信號(hào)和wrclk信號(hào),在手動(dòng)設(shè)置中也可以通過(guò)設(shè)置參數(shù)WRITE_ACLR_SYNCH有效實(shí)現(xiàn)同步。
 wrusedw: 輸出數(shù)據(jù)顯示FIFO中寫(xiě)入數(shù)據(jù)的數(shù)據(jù)量,7 bit。在DCFIFO中,其端口寬度要與手動(dòng)設(shè)置參數(shù)LPM_WIDTHU相等。需要注意的是,對(duì)于Cyclone系列的FPGA,當(dāng)顯示數(shù)據(jù)滿時(shí)實(shí)際FIFO的存儲(chǔ)數(shù)據(jù)量有可能并沒(méi)有達(dá)到存儲(chǔ)的最大值,因此應(yīng)該參考full或者wrfull端口來(lái)執(zhí)行正確的寫(xiě)操作,參考empty或者rdempty端口來(lái)執(zhí)行正確的讀操作。
  FIFO工作時(shí)的狀態(tài)轉(zhuǎn)換圖如圖3所示。


2.2 FIFO輸入時(shí)序模塊
    FIFO輸入時(shí)序模塊用于使TIME COUNTER_WITHID模塊(用于對(duì)信號(hào)進(jìn)行計(jì)時(shí))的輸出數(shù)據(jù)時(shí)序符合FIFO輸入數(shù)據(jù)時(shí)序的接口要求。FIFO的寫(xiě)請(qǐng)求信號(hào)wrreq的長(zhǎng)度所包含的周期數(shù)為寫(xiě)入FIFO數(shù)據(jù)的個(gè)數(shù),并且寫(xiě)請(qǐng)求信號(hào)有效時(shí)輸入數(shù)據(jù)端口的數(shù)據(jù)即被寫(xiě)入FIFO,幾乎沒(méi)有延遲,因此輸入數(shù)據(jù)長(zhǎng)度應(yīng)符合FIFO的時(shí)鐘信號(hào)周期長(zhǎng)度,并且寫(xiě)請(qǐng)求信號(hào)長(zhǎng)度也需要符合FIFO時(shí)鐘信號(hào)周期長(zhǎng)度,兩者在時(shí)間上需要達(dá)到同步。
 為了測(cè)試輸出數(shù)據(jù)長(zhǎng)度符合FIFO接口的時(shí)序要求,分別進(jìn)行了圖4(a)、圖4(b)、圖4(c)所示的波形仿真,輸入信號(hào)的持續(xù)長(zhǎng)度分別為10 ns、20 ns、30 ns。如圖所示三種輸入的輸出皆為10 ns保持時(shí)間, 同時(shí)輸出10 ns write_en使能信號(hào),符合FIFO接口要求。

 參照時(shí)序仿真圖對(duì)各個(gè)端口及時(shí)序進(jìn)行說(shuō)明:
 clk: FIFO輸入時(shí)序模塊時(shí)鐘輸入,1 bit。用于同步datain16、dataout16、write_en、complete、enable、usedw信號(hào),與FIFO中wrclk信號(hào)相頻相同。
    complete:數(shù)據(jù)接收完成信號(hào),1 bit。當(dāng)從timecounter_withid模塊接收數(shù)據(jù)完成時(shí)發(fā)送給timecounter_withid模塊,持續(xù)時(shí)長(zhǎng)1周期,timecounter_withid接收到信號(hào)時(shí)對(duì)內(nèi)部數(shù)據(jù)進(jìn)行初始化,做好準(zhǔn)備對(duì)新的信號(hào)輸入進(jìn)行計(jì)時(shí)。
 datain16: 16位數(shù)據(jù)輸入端口,16 bit。用于輸入16位數(shù)據(jù)。
 dataout16: 16位數(shù)據(jù)輸出端口,16 bit。用于輸出16位數(shù)據(jù)。
 enable: 數(shù)據(jù)寫(xiě)入輸入時(shí)序模塊使能,1 bit。當(dāng)有數(shù)據(jù)要寫(xiě)入模塊時(shí),首先在此端口輸入高電平,然后輸入數(shù)據(jù)方可被正確接收。
    usedw: FIFO中可用數(shù)據(jù)端口,7 bit。用于查詢FIFO中可用數(shù)據(jù)量,以此判斷是否有剩余空間,進(jìn)而確定是否向FIFO中寫(xiě)數(shù)據(jù)。
 write_en: 寫(xiě)使能信號(hào),1bit。當(dāng)要向FIFO中寫(xiě)入數(shù)據(jù)時(shí)為高電平,dataout16進(jìn)行數(shù)據(jù)輸出。
 模塊流程圖如圖5所示。整個(gè)模塊以posedge clk為同步時(shí)鐘進(jìn)行循環(huán)判斷。當(dāng)寫(xiě)入模塊使能enable有效且寫(xiě)入完畢信號(hào)complete無(wú)效時(shí),首先對(duì)FIFO的存儲(chǔ)情況進(jìn)行判斷,當(dāng)FIFO中字節(jié)小于7 FH時(shí),說(shuō)明FIFO未被寫(xiě)滿,可以向其輸入數(shù)據(jù)。數(shù)據(jù)輸入FIFO的同時(shí),寫(xiě)使能write_en置為有效,對(duì)于計(jì)數(shù)模塊輸出的接收完成信號(hào)complete信號(hào)置為有效。對(duì)一個(gè)周期的計(jì)數(shù)寄存器counter進(jìn)行判斷,如果counter為1即計(jì)時(shí)一周期到,則寄存器清零。counter的數(shù)值是在每次posedge clk信號(hào)到來(lái)時(shí)對(duì)write_en判斷后進(jìn)行加1,即當(dāng)write_en為有效電平時(shí),counter才被加1用來(lái)計(jì)時(shí)。

2.3 FIFO輸出時(shí)序模塊
 FIFO輸出時(shí)序模塊用于使FIFO的輸出信號(hào)符合Flash讀寫(xiě)時(shí)序規(guī)范。由FIFO的時(shí)序仿真圖可以看出,當(dāng)讀請(qǐng)求信號(hào)rdreq發(fā)出后,數(shù)據(jù)要延遲15 ns左右才會(huì)輸出。如果Flash控制器發(fā)出讀請(qǐng)求信號(hào)后立即讀FIFO的數(shù)據(jù),則會(huì)造成差錯(cuò)。因此本模塊可以銜接FIFO和Flash控制器的端口時(shí)序。
 時(shí)序仿真圖如圖6所示。由時(shí)序仿真圖可以看出,sendready信號(hào)相對(duì)于empty信號(hào)延遲了13 ns左右,dataout信號(hào)相對(duì)于sendready信號(hào)延遲了20 ns左右。可以滿足FIFO時(shí)序要求。

 參照時(shí)序仿真圖對(duì)各個(gè)端口及時(shí)序進(jìn)行說(shuō)明:
 clk: FIFO輸入時(shí)序模塊時(shí)鐘輸入,1 bit。用于同步datain16、dataout16、empty、complete、sendready信號(hào),與FIFO中wrclk信號(hào)相頻相同。
    complete: 數(shù)據(jù)接收完成信號(hào),1 bit。當(dāng)存儲(chǔ)器從模塊接收數(shù)據(jù)完成時(shí)發(fā)送給FIFO輸出時(shí)序模塊,持續(xù)時(shí)長(zhǎng)1周期,F(xiàn)IFO輸出時(shí)序模塊接收到此信號(hào)時(shí)對(duì)內(nèi)部數(shù)據(jù)進(jìn)行初始化,開(kāi)始對(duì)新的信號(hào)輸入進(jìn)行計(jì)時(shí)。
 datain16: 16位數(shù)據(jù)輸入端口,16 bit。用于輸入16位數(shù)據(jù)。
 dataout16: 16位數(shù)據(jù)輸出端口,16 bit。用于輸出16位數(shù)據(jù)。    
 sendready: 數(shù)據(jù)準(zhǔn)備就緒信號(hào),1 bit。當(dāng)FIFO收到讀數(shù)據(jù)請(qǐng)求信號(hào)并且有數(shù)據(jù)輸出時(shí)sendready發(fā)送有效信號(hào),存儲(chǔ)器讀取輸出數(shù)據(jù)。
 write_en: 寫(xiě)使能信號(hào),1 bit。當(dāng)要向FIFO中寫(xiě)入數(shù)據(jù)時(shí)為高電平,dataout16進(jìn)行數(shù)據(jù)輸出。
 模塊流程圖如圖7所示。

3 設(shè)計(jì)驗(yàn)證
3.1  驗(yàn)證

    將程序下載到CycloneII FPGA芯片中,并且用按鍵作為輸入信號(hào)進(jìn)行了測(cè)試,在控制臺(tái)上打印 FIFO的輸出結(jié)果,結(jié)果如圖8所示。

3.2  驗(yàn)證結(jié)果分析
 當(dāng)FIFO中的數(shù)據(jù)經(jīng)過(guò)時(shí)序轉(zhuǎn)換模塊可以向存儲(chǔ)器寫(xiě)入時(shí),ready信號(hào)為1,同時(shí)數(shù)據(jù)寫(xiě)入寄存器并在控制臺(tái)打印,然后再將寄存器中的數(shù)據(jù)寫(xiě)入存儲(chǔ)器,并向時(shí)序轉(zhuǎn)換模塊返回完成信號(hào)complete高電平1,等待從FIFO中讀取新數(shù)據(jù)。經(jīng)過(guò)驗(yàn)證,本設(shè)計(jì)可以將輸入數(shù)據(jù)進(jìn)行緩存并且輸出給NiosII CPU,符合像素探測(cè)器數(shù)據(jù)緩存要求。
    本文結(jié)合FPGA對(duì)粒子對(duì)撞機(jī)像素探測(cè)器的數(shù)據(jù)緩存提出了解決方法并進(jìn)行了設(shè)計(jì)和驗(yàn)證。相對(duì)于傳統(tǒng)的ASIC構(gòu)建方法,具有高效率、低投入的優(yōu)勢(shì),提高了探測(cè)器升級(jí)和參數(shù)設(shè)置等方面的靈活性;相對(duì)于分立元件的設(shè)計(jì)方法,具有高速、低功耗和集成度高的優(yōu)勢(shì)。通過(guò)時(shí)序仿真和FPGA驗(yàn)證,能夠滿足像素探測(cè)器的設(shè)計(jì)需要。
參考文獻(xiàn)
[1] 王直華. 粒子對(duì)撞機(jī)與宇宙大爆炸[J]. 科技潮,2008(10).
[2] 孟祥承.新型半導(dǎo)體探測(cè)器發(fā)展和應(yīng)用[J]. 核電子學(xué)與探測(cè)技術(shù),2004,24(1):87-96.
[3] GAGLIARD G. The ATLAS pixel detector electronics.Nuclear Instruments and Methods in Physics Research,2001:275-281.
[4] 楊海鋼,孫嘉斌,王慰.FPGA器件設(shè)計(jì)技術(shù)發(fā)展綜述[J].電子與信息學(xué)報(bào),2010,32(3):714-727.

此內(nèi)容為AET網(wǎng)站原創(chuàng),未經(jīng)授權(quán)禁止轉(zhuǎn)載。
亚洲一区二区欧美_亚洲丝袜一区_99re亚洲国产精品_日韩亚洲一区二区
国产精品美女久久久久久免费 | 久久综合伊人77777蜜臀| 亚洲一二三区在线观看| 亚洲人成网站777色婷婷| 久久精品欧美| 久久精品视频在线播放| 欧美一区二区高清在线观看| 亚洲欧美国产精品桃花| 亚洲专区免费| 亚洲一二三区在线观看| 亚洲香蕉网站| 亚洲伊人色欲综合网| 一区二区免费在线视频| 91久久精品国产91久久性色tv | 亚洲人成欧美中文字幕| 亚洲三级免费| 日韩亚洲视频| 夜夜嗨av色一区二区不卡| 一区二区三区高清在线观看| 一区二区三区鲁丝不卡| 宅男噜噜噜66国产日韩在线观看| 夜夜嗨av一区二区三区中文字幕| 亚洲国产高清视频| 亚洲福利免费| 亚洲人体偷拍| 一区二区三区回区在观看免费视频| 亚洲日本一区二区| 99国产精品私拍| 亚洲欧美国产77777| 欧美影院精品一区| 亚洲人成人99网站| 中文国产亚洲喷潮| 亚洲欧美中文字幕| 久久久一本精品99久久精品66| 久久久五月天| 欧美国产另类| 欧美午夜激情小视频| 国产精品国产自产拍高清av王其| 国产精品综合av一区二区国产馆| 国产一区视频在线观看免费| 亚洲大片精品永久免费| 99re热精品| 香蕉精品999视频一区二区 | 91久久国产精品91久久性色| 一区二区精品在线| 欧美一区影院| 农夫在线精品视频免费观看| 欧美日韩一区二区免费在线观看 | 国内精品久久久久影院薰衣草| 有坂深雪在线一区| 亚洲精品乱码| 亚洲欧美日韩系列| 亚洲人成在线观看一区二区| 亚洲一区二区视频在线| 久久久久久一区| 欧美精品一区二区三区蜜桃| 国产精品久久久久久久久借妻| 国产日韩欧美高清免费| 亚洲高清免费视频| 亚洲一区精品视频| 亚洲国产91精品在线观看| 在线视频精品一| 久久精品在线视频| 欧美成年人视频网站| 国产精品黄视频| 在线精品观看| 亚洲直播在线一区| 99精品福利视频| 久久久999精品免费| 欧美日韩精品二区| 国产视频欧美视频| 夜夜嗨av一区二区三区| 亚洲第一网站免费视频| 亚洲一区二区在线看| 久久中文字幕一区| 国产精品日韩欧美一区二区| 亚洲经典自拍| 久久高清福利视频| 亚洲欧美日韩综合aⅴ视频| 欧美国产欧美亚洲国产日韩mv天天看完整 | 亚洲午夜精品一区二区三区他趣| 亚洲高清资源| 久久aⅴ国产欧美74aaa| 欧美美女bbbb| 1769国产精品| 欧美一区视频在线| 午夜国产精品视频| 欧美男人的天堂| 亚洲高清三级视频| 欧美在线观看一区二区三区| 亚洲一区二区三区影院| 欧美黄色aaaa| 一区二区在线视频| 欧美一区精品| 欧美亚洲日本网站| 国产精品va| 99精品国产在热久久婷婷| 亚洲精品日韩一| 久久久噜噜噜久久中文字免| 国产精品三级视频| 亚洲深夜福利| 亚洲一区二区黄| 欧美日本在线播放| 亚洲激情网址| 亚洲精品日韩激情在线电影 | 在线看视频不卡| 久久国产日韩欧美| 久久国产精品99精品国产| 国产精品毛片va一区二区三区 | 久久一区欧美| 国内精品**久久毛片app| 亚洲欧美日韩在线高清直播| 亚洲免费视频在线观看| 欧美日韩一区二区三区四区五区| 亚洲国产精品一区二区三区| 亚洲国产精品一区制服丝袜| 久久婷婷蜜乳一本欲蜜臀| 国内精品视频一区| 欧美在线播放高清精品| 久久精品国产99精品国产亚洲性色 | 欧美亚洲综合网| 久久国产夜色精品鲁鲁99| 国产欧美日本| 欧美一区二区福利在线| 久久不射中文字幕| 国产日产高清欧美一区二区三区| 亚洲一区二区高清视频| 亚洲欧美在线另类| 国产日韩成人精品| 性高湖久久久久久久久| 久久久久久夜| 亚洲高清在线| 一本色道**综合亚洲精品蜜桃冫| 欧美日韩国产免费| 一区二区三区欧美成人| 欧美亚洲日本国产| 国产又爽又黄的激情精品视频| 欧美一区二区在线看| 久久婷婷国产综合尤物精品| 1769国产精品| 夜夜精品视频| 国产精品区一区二区三| 午夜视频一区| 男男成人高潮片免费网站| 亚洲精品系列| 香蕉久久夜色| 国语自产精品视频在线看抢先版结局| 亚洲第一久久影院| 欧美日本精品在线| 一级日韩一区在线观看| 久久福利毛片| 亚洲国产精品久久久| 亚洲一区二区在线免费观看视频| 国产精品三区www17con| 亚洲高清视频一区| 欧美日本一区二区三区| 亚洲欧美激情诱惑| 欧美xx69| 亚洲一级一区| 美国成人毛片| 亚洲天堂成人在线观看| 久久久欧美精品sm网站| 亚洲精品国产无天堂网2021| 午夜日韩激情| 在线精品亚洲| 亚洲欧美日韩在线高清直播| 黑丝一区二区三区| 中文国产亚洲喷潮| 国内精品免费午夜毛片| 中国成人黄色视屏| 国产午夜精品久久| 一本到高清视频免费精品| 国产日韩在线播放| 一本久道久久综合婷婷鲸鱼| 国产亚洲毛片在线| 中文在线不卡| 韩国av一区二区三区在线观看| 一本色道久久综合狠狠躁篇的优点| 国产美女一区| 99在线精品观看| 国内揄拍国内精品久久| 亚洲一区二区在线免费观看| 在线观看一区| 欧美一区日韩一区| 亚洲人体偷拍| 久久久久久久尹人综合网亚洲| 亚洲精品一二三| 六月婷婷一区| 午夜在线电影亚洲一区| 欧美日本一区二区三区| 亚洲高清自拍| 国产视频精品免费播放| 一区二区三欧美| 伊人久久噜噜噜躁狠狠躁 | 国产亚洲a∨片在线观看| 一区二区三区日韩欧美精品| 尹人成人综合网| 欧美在线亚洲综合一区| 99这里有精品| 欧美激情综合网|