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

    基于響應(yīng)矩陣的RFQ新調(diào)諧算法的驗證與實現(xiàn)

    2023-12-27 02:59:14陸羿行方文程張俊強郭玉森趙振堂
    核技術(shù) 2023年11期
    關(guān)鍵詞:調(diào)諧器四極腔體

    陸羿行 方文程 張俊強 郭玉森,5 趙振堂

    1(中國科學(xué)院上海應(yīng)用物理研究所 上海 201800)

    2(中國科學(xué)院大學(xué) 北京 101408)

    3(中國科學(xué)院上海高等研究院 上海 201204)

    4(上海市低溫超導(dǎo)高頻腔技術(shù)重點實驗室 上海 201800)

    5(上??萍即髮W(xué) 上海 201210)

    射頻四極場加速器(Radio Frequency Qaudrupole,RFQ)[1-4]常用于低能質(zhì)子的加速,可以保證其束流品質(zhì),如今廣泛應(yīng)用于質(zhì)子加速器及相關(guān)裝置中?,F(xiàn)有裝置上的RFQ工作頻率通常在200~400 MHz,而歐洲核子中心(European Organization for Nuclear Research,CERN)研發(fā)的用于質(zhì)子治療裝置的RFQ工作頻率為750 MHz[5],這種高頻率的RFQ具有更高的加速梯度,可以推動質(zhì)子加速器的小型化。近十多年中,上海光源(Shanghai Synchrotron Radiation Facility,SSRF)也在致力于質(zhì)子治療裝置的研發(fā),為瑞金醫(yī)院腫瘤質(zhì)子中心設(shè)計、建造的首臺國產(chǎn)質(zhì)子治療裝置已經(jīng)通過驗收,進入臨床試驗階段[6-8];同時,提出了基于直線加速器的質(zhì)子閃療裝置[9-11],緊湊型RFQ作為裝置中的一部分,也已經(jīng)過了多年的研究。調(diào)諧是緊湊型RFQ的關(guān)鍵技術(shù)之一,傳統(tǒng)的調(diào)諧算法基于腔體理論計算得到的物理特性,用于緊湊型RFQ效果欠佳。CERN在高頻RFQ上利用實際場分布對調(diào)諧棒的響應(yīng)量,即響應(yīng)矩陣,來進行調(diào)諧[12-13],在求解的過程中應(yīng)用了奇異值分解(Singular Value Decomposition,SVD)方法。然而,這種方法仍然存在不足,在求解過程中常常出現(xiàn)調(diào)諧量超出調(diào)諧棒可調(diào)節(jié)范圍的情況,即使對SVD方法進行了改造解決了這個問題,又會造成精度下降,導(dǎo)致調(diào)諧步驟多、二極場分量難以降低等問題。

    針對以上問題,本文提出了一種在響應(yīng)矩陣的基礎(chǔ)上利用最小二乘法(Least Squares,LS)進行求解的調(diào)諧算法。該算法可以限定求解的范圍,得到其中的最優(yōu)解;同時,對腔體的二極場分量賦予權(quán)重,可以有效地對二極場進行調(diào)諧,解決二極場誤差難以降低的問題?;谶@種新型的RFQ調(diào)諧算法,在Matlab下進行了算法的開發(fā),完成了模擬條件下的調(diào)諧實驗,搭建了基于微擾法的測量平臺,編寫了控制和通信的Labview程序,并在RFQ樣機上完成了算法的最終驗證。

    經(jīng)過實驗驗證,新調(diào)諧算法可以在限定范圍內(nèi)得到最優(yōu)解,并有效地解決二極場誤差。該算法具有普適性,可以應(yīng)用到其他頻率的RFQ中,解決更多RFQ的調(diào)諧問題;將來可以促進質(zhì)子加速器的小型化,推動醫(yī)用質(zhì)子裝置的發(fā)展,普及相關(guān)診療手段(質(zhì)子治療、質(zhì)子CT(Computed Tomography)[14-15]、硼中子俘獲療法(Boron Neutron Capture Therapy,BNCT)[16]、醫(yī)用同位素生產(chǎn)[17]等)。

    1 RFQ及其測量方法

    四翼型RFQ的腔體是在諧振腔的腔壁上添加4個極頭形成的,其工作模式是類TE210模。在橫向上,4個極頭間可以激發(fā)出電四極場,腔體中的電場也主要集中在極頭附近;在遠(yuǎn)離極頭的位置,主要為磁場,方向沿軸向(圖1(a))。當(dāng)束流在軸線上通過時,會受到四極電場的聚焦作用,并且該電場的極性會隨時間和粒子的位置變化,因此,每過半個射頻周期粒子就會受到相反的作用。在縱向上,極頭與軸線間的距離最小為a,最大為ma,m為調(diào)制參數(shù)。經(jīng)過調(diào)制后,極頭呈類三角函數(shù)的形狀,形成縱向的電場,實現(xiàn)對束流的聚束和加速作用。極頭形狀變化的周期為βλ,相對極頭的“峰”和“谷”是對應(yīng)的,而相鄰極頭的“峰”和“谷”則是錯開的(圖1(b))。

    圖1 RFQ結(jié)構(gòu)及測量框架圖(a) RFQ橫向形狀及電磁場分布,(b) RFQ極頭縱向形狀,(c) RFQ測量方法框架圖Fig.1 Structure of an RFQ and the framework of the measurement(a) Transverse section of an RFQ and the electromagnetic field on it, (b) Longitudinal section of an RFQ, (c) Framework of measuring an RFQ

    測量腔體中的場分布是RFQ調(diào)諧的第一步,可以利用微擾法實現(xiàn)。RFQ的最小束流孔徑一般都非常小,在1~3 mm,難以測量其軸線上的電場分布,因此,RFQ是通過測量4個象限中的磁場來得到場分布的。用步進電機帶動魚線,魚線又帶動金屬小球在腔體中運動,金屬小球會對腔體中的電磁場分布產(chǎn)生影響,探針測量得到的S12信號也會發(fā)生改變,由此可以推算出小球所在位置的磁場強度,存在關(guān)系[18]:

    式中:B表示磁場強度;ΔP表示小球的擾動引起的相位差。電機控制器和網(wǎng)絡(luò)分析儀可以與計算機進行通信來傳遞數(shù)據(jù)(圖1(c))。

    2 基于響應(yīng)矩陣和最小二乘法的調(diào)諧算法

    加速器調(diào)諧的基本邏輯是利用調(diào)諧器對腔體的擾動,彌補腔體內(nèi)場分布的誤差。通過對腔體的測量,得到場分布及其與理想狀態(tài)的差別,再利用腔體的物理特性,計算得到調(diào)諧器所需的調(diào)諧量。RFQ的主要調(diào)諧方法可以分為兩大類:一種基于理論計算得到的物理特性,例如文獻[19],可以稱為“理論法”;另外一種基于實際測量得到的物理特性,例如文獻[12-13],可以稱為“解析法”。理論法的缺點在于把理想狀態(tài)下計算得到的數(shù)據(jù)應(yīng)用于實際,忽略了實際腔體和理想模型之間的差別;對于敏感度更高的緊湊型RFQ而言,實際情況與理想狀態(tài)相差較大,因此理論法難以適用。解析法則可以消除理論和實際之間的誤差對調(diào)諧造成的影響,可以解決緊湊型RFQ的調(diào)諧問題;解析法的缺點是需要額外測量響應(yīng)矩陣。下面將詳細(xì)介紹基于響應(yīng)矩陣和最小二乘法的新調(diào)諧算法。

    分別用B1、B2、B3、B4來表示各象限中的磁場(圖1(a)),則可以計算得到腔體中的二、四極場:

    式中:Q為四極場;Ds、Dt為兩個方向上的二極場。在實際測量和計算中,以上場量均使用歸一化值,其歸一化誤差定義為:

    將二、四極場組合起來形成場矩陣:

    式中:n為沿軸線取場的點數(shù),并滿足場矩陣中元素的數(shù)量大于等于調(diào)諧器數(shù)量,即3n≥ntuner,因為這樣可以提高調(diào)諧的精度。RFQ初始狀態(tài)的場矩陣記為Fi,接著分別移動每一個調(diào)諧器并重復(fù)測量對應(yīng)的場分布,記為F1,F(xiàn)2,…,F(xiàn)n-tuner,可以得到不同場分量對調(diào)諧器的響應(yīng),即響應(yīng)矩陣:

    式中:d為調(diào)諧器的移動量。對于一個理想的RFQ,需要四極場的值都相等,且二極場大小為0,即目標(biāo)矩陣:

    因此,調(diào)諧器對場分布的補償量應(yīng)為:

    至此,調(diào)諧問題轉(zhuǎn)化為求解方程:

    式中:T即為需要求解的調(diào)諧量矩陣。

    從上文可以發(fā)現(xiàn),響應(yīng)矩陣M通常為一個非正方形矩陣,因此無法直接求解。CERN使用SVD方法來求解,但該方法無法限定調(diào)諧量的范圍,直接求解得到的調(diào)諧量非常大,超出調(diào)諧棒可調(diào)諧的范圍,因此,在求解過程中將特征值向量中的部分元素代換為0。這種改造后的SVD方法雖然解決了結(jié)果的范圍問題,但是舍去部分特征值會影響到結(jié)果的精度,并可能放大二極場誤差難以減小的問題。

    對SVD方法做嘗試時發(fā)現(xiàn)了上述問題。針對結(jié)果的范圍問題,可以利用最小二乘法對調(diào)諧量矩陣T中的每一個元素進行單獨求解,得到這一元素在限定范圍內(nèi)使誤差最小的值,可以用式(9)表示:

    式中:tn為矩陣T中的第n個元素。依次求解矩陣T中的元素,并對這個過程進行迭代,直到所有誤差達(dá)到預(yù)期的目標(biāo)為止。

    在這種方法中,還可以對二極場分量賦予權(quán)重,以解決二極場誤差難以減小的問題,即將原方程變?yōu)榍蠼猓?/p>

    式中:k為給二極場賦予的權(quán)重。

    整個調(diào)諧算法如圖2所示。以上為調(diào)諧實驗過程中針對一些問題給出的新解法,下面從模擬和實驗兩方面闡述對新調(diào)諧算法的驗證。

    圖2 基于響應(yīng)矩陣和最小二乘法的調(diào)諧算法流程圖Fig.2 Flowchart of the tuning algorithm based on response matrix and least squares method

    3 基于模擬軟件的調(diào)諧實驗與算法驗證

    為驗證新調(diào)諧算法,首先利用三維電磁場模擬軟件進行了大量模擬調(diào)諧實驗。RFQ的三維模型基于上海光源設(shè)計的714 MHz RFQ,在正式設(shè)計中,該RFQ的長度為2 m,可以將質(zhì)子從50 keV加速至4 MeV,而用于模擬實驗的模型長為1 m,橫截面形狀與正式設(shè)計相似。RFQ的實驗?zāi)P腿鐖D3(a)所示,為金屬模型,在該腔體上共有16根調(diào)諧棒,如圖中所指部分,其余為耦合器、真空口及探測器預(yù)留位置。為實驗?zāi)P椭匦略O(shè)計了底切,在所有調(diào)諧棒均插入5 mm深的狀態(tài)下,將其頻率調(diào)至714 MHz,且將場分布調(diào)至理想狀態(tài)。讀取腔體4個象限中的磁場分布,在不同象限中讀取的位置是對稱的,其中一個位置如圖3(b)所示。對場分布處理后,得到腔體的二、四極場分布如圖3(c)所示,即為理想狀態(tài)。圖中四極場上的“x”標(biāo)記為后續(xù)調(diào)諧實驗過程中,選取場值用于調(diào)諧算法的位置,即n=6。

    圖3 模擬條件下RFQ的模型及場分布(a) RFQ模型,(b) RFQ截面圖及選取場分布位置示意,(c) RFQ二、四極場分布(Qo為四極場,Ds、Dt分別為水平和豎直方向的二極場)Fig.3 3-D model of the RFQ and its field distribution in simulation(a) Model of the RFQ, (b) Transverse section of the RFQ and the place where the field was picked, (c) Dipole and quadrupole components of the RFQ

    根據(jù)算法在Matlab中編寫了相應(yīng)的程序。在理想狀態(tài)下,將每根調(diào)諧棒的插入深度分別減少和增加3 mm,重復(fù)獲取模型中的場分布,得到每根調(diào)諧棒在6 mm范圍內(nèi)對場的影響,計算出的響應(yīng)矩陣為18×16的矩陣。

    3.1 單調(diào)諧棒擾動模擬調(diào)諧

    首先進行單調(diào)諧棒擾動模擬。在理想狀態(tài)下,將第一根調(diào)諧棒t1的深度增加2 mm,對場產(chǎn)生擾動,破壞理想狀態(tài),通過基于響應(yīng)矩陣和LS方法來對場進行模擬調(diào)諧實驗。

    在單根調(diào)諧棒擾動的情況下,雖然在算法中加入了對調(diào)諧量范圍的限制,但是在實際計算的過程中并未出現(xiàn)解出的調(diào)諧量大于調(diào)諧棒可調(diào)諧范圍的情況。在給二極場誤差加上不同的權(quán)重后,對比得到的結(jié)果可以發(fā)現(xiàn)(圖4),隨著權(quán)重的增加,調(diào)諧棒t1的解越接近-2 mm;而不論權(quán)重如何變化,其他調(diào)諧棒解出的調(diào)諧量絕對值均非常小,都在0.1 mm以內(nèi)。由于在數(shù)值上該結(jié)果就已經(jīng)十分接近擾動量,因此沒有再將該結(jié)果應(yīng)用到模型中進行模擬。基本可以推論,僅通過一次調(diào)諧,這些誤差就已經(jīng)大幅降低。

    圖4 單調(diào)諧棒擾動下不同權(quán)重得到的調(diào)諧量Fig.4 Solution with the different weight for dipole components under single tuners' perturbation

    通過單調(diào)諧棒擾動下的模擬調(diào)諧,初步驗證了調(diào)諧算法的可行性,并且實現(xiàn)了較為精準(zhǔn)的調(diào)諧效果。

    3.2 多調(diào)諧棒擾動模擬調(diào)諧

    單調(diào)諧棒擾動實驗初步驗證了算法的準(zhǔn)確性,但這種誤差情況與實際不符,因此進一步進行了多調(diào)諧棒擾動模擬。通過在[-3 mm,3 mm]范圍內(nèi)生成隨機數(shù),為每一根調(diào)諧棒賦予了“擾動”,作為調(diào)諧的初始狀態(tài)。調(diào)諧過程中的歸一化誤差用下式定義:

    初始的二、四極場分布如圖5(b、c)所示,初始的二極場誤差較大,接近30%。

    圖5 模擬調(diào)諧實驗結(jié)果及對比(a) 不同權(quán)重下得到的二、四極場誤差比較,(b) 調(diào)諧過程中四極場分布比較,(c) 調(diào)諧過程中二極場分布比較Fig.5 Result comparisons of the tuning experiment in simulation(a) Quadrupole and dipole component errors with different weights, (b) Quadrupole field variation during the tuning,(c) Dipole field variation during the tuning

    同樣利用LS調(diào)諧算法來進行調(diào)諧。首先比較了二極場誤差不同權(quán)重下的結(jié)果(圖5(a)),圖5中q為四極場誤差,ds、dt分別為水平和豎直方向的二極場誤差。與單調(diào)諧棒擾動的情況不同,當(dāng)二極場誤差權(quán)重達(dá)到5之后,四極場誤差明顯增長,迅速增長到2%以上,當(dāng)權(quán)重增加到10以后,誤差達(dá)到了6%。為了詳細(xì)了解權(quán)重對結(jié)果的影響,補充了權(quán)重為2、3的結(jié)果,可以發(fā)現(xiàn)這兩種情況下,四極場誤差會有所增長,同時二極場誤差會降低,并且在這個范圍內(nèi),四極場誤差均不超過2%,二極場誤差也均在2%以內(nèi)。幾種權(quán)重相比而言,2倍的權(quán)重可以在保證較小的二極場誤差的同時,讓四極場誤差在1%左右,是比較好的結(jié)果。以上模擬實驗表明,多調(diào)諧棒擾動的情況和單調(diào)諧棒是有明顯區(qū)別的,多調(diào)諧棒擾動下進行模擬調(diào)諧會更加接近實際情況,同時也說明為二極場誤差賦予權(quán)重的方法仍然是有效的。

    進一步地,在權(quán)重為2的調(diào)諧結(jié)果的基礎(chǔ)上進行了第二次調(diào)諧,在第二次調(diào)諧中,賦予二極場誤差的權(quán)重依舊設(shè)為2。兩次調(diào)諧結(jié)果,以及初始場、目標(biāo)場的比較見圖5(b、c)。圖5中,q為四極場誤差,ds、dt分別為水平和豎直方向的二極場誤差,前綴o表示目標(biāo)場,i表示初始場,l1表示第一次調(diào)諧后結(jié)果,l2表示第二次調(diào)諧后結(jié)果。從圖5中可以發(fā)現(xiàn),經(jīng)過一次調(diào)諧后,四極場分布已經(jīng)和目標(biāo)十分接近,同時二極場誤差大大減小。經(jīng)過兩次調(diào)諧后,二、四極場誤差均又在縮小,最終四極場誤差為0.8%,二極場誤差為0.4%,是非常理想的誤差水平。

    經(jīng)過多調(diào)諧棒擾動模擬調(diào)諧實驗,完全驗證了調(diào)諧算法的可行性,驗證了二極場權(quán)重對消除二極場的作用,預(yù)示著調(diào)諧方法在實際調(diào)諧中的優(yōu)秀效果。

    3.3 LS解法與SVD解法比較

    對基于SVD解法的調(diào)諧算法也進行了簡單驗證。不可避免地,通過SVD解法得到的調(diào)諧范圍很大,遠(yuǎn)遠(yuǎn)超出了調(diào)諧棒的可移動范圍,必須利用改造后的SVD方法。根據(jù)算法,一共可以得到16組結(jié)果,其中前幾組結(jié)果存在同樣的問題,而后幾組結(jié)果又存在丟失精度過多的問題,因此挑選了其中第9組來進行調(diào)諧。經(jīng)過調(diào)諧后場分布誤差有明顯改善,四極場誤差在1%以內(nèi),最大二極場誤差約2%;第二次調(diào)諧依舊選用第9組結(jié)果,經(jīng)過第二次調(diào)諧后,二極場誤差減小(0.8%),而四極場誤差增大(2%)。通過SVD方法兩次調(diào)諧后四極場誤差偏大,必須經(jīng)過更多的調(diào)諧次數(shù)來減小誤差;而應(yīng)用二極場誤差權(quán)重加LS方法,如上一節(jié)所述,兩次調(diào)諧后實現(xiàn)了較為理想的場分布。

    最終SVD方法和LS方法得到的二、四極場誤差比較如圖6所示,通過比較結(jié)果,驗證了LS方法可限定求解范圍、高效處理二極場誤差、縮短調(diào)諧過程等優(yōu)點。

    圖6 不同調(diào)諧算法下二、四極場誤差對比Fig.6 Comparison of the quadrupole and dipole component errors with different tuning methods

    4 基于RFQ樣機的調(diào)諧算法驗證

    以上完成了調(diào)諧算法的初步驗證,根據(jù)模擬調(diào)諧實驗的RFQ模型,加工了一段鋁模樣段(圖7(a)),利用RFQ樣機進一步驗證調(diào)諧方法。

    圖7 RFQ樣段調(diào)諧平臺及實驗結(jié)果(a) 714 MHz RFQ樣機,(b) RFQ測量平臺,(c) 電機控制和微波測量程序面板,(d) RFQ樣段初始場,(e) 調(diào)諧過程中二、四極場誤差變化Fig.7 Testbed of the RFQ prototype and the experimental results(a) Prototype of the 714-MHz RFQ, (b) Testbed of the prototype, (c) User interface of the motor control and RF measurement,(d) Initial field of the RFQ prototype, (e) Variation of quadrupole and dipole component errors during the tuning

    RFQ的測量方法如§1所述,具體平臺搭建見圖7(b)。在RFQ樣機的端板上開有4個測量孔,分別在RFQ的4個象限,魚線帶著金屬小球在4個象限中來回穿過,并繞過固定在端板上的滑輪,利用步進電機帶動小球在腔體中勻速運動。通過腔體上預(yù)留的耦合孔、真空抽口以及探測器孔將探針伸入腔體內(nèi),探針的另一頭連接到矢量網(wǎng)絡(luò)分析儀,獲取腔體中的微波參數(shù)。步進電機控制及網(wǎng)絡(luò)分析儀的數(shù)據(jù)傳輸均通過在LabView環(huán)境下編寫的程序來實現(xiàn)(圖7(c))。測量得到的數(shù)據(jù)先進行預(yù)處理,再利用調(diào)諧算法完成調(diào)諧量的計算,數(shù)據(jù)處理及計算均通過Matlab程序?qū)崿F(xiàn)。RFQ樣機上有16個調(diào)諧器,調(diào)諧器中螺桿帶動調(diào)諧頭移動,可以實現(xiàn)較為精準(zhǔn)的調(diào)諧量。

    對RFQ樣段的初始場進行測量。先將每個調(diào)諧器的插入深度都調(diào)節(jié)為5 mm,使用兩根探針測量腔體內(nèi)的微波參數(shù),處理后得到RFQ樣段的初始場分布,其中,二極場誤差在30%以上(圖7(d))。主要原因在于初始狀態(tài)調(diào)諧器的深度不一致,誤差主要來源于兩方面:一是調(diào)諧器與腔體間沒有焊接,導(dǎo)致調(diào)諧器與腔體的間隙不一樣,存在誤差;二是調(diào)諧器的精度未達(dá)到預(yù)想的程度,難以判斷各調(diào)諧器的深度是否一致。以上誤差只會影響初始誤差,不會影響調(diào)諧過程,因此在此樣段上進行調(diào)諧實驗來驗證算法及調(diào)諧平臺的可靠性仍是可行的。

    利用上述調(diào)諧算法對RFQ模型腔進行調(diào)諧。初始狀態(tài)下二、四極場誤差分別為24.09%和1.57%,經(jīng)過5次調(diào)諧后兩者誤差分別降為2.33%和1.39%。此時誤差已經(jīng)達(dá)到理想的水平,對腔體進行了一次頻率調(diào)諧,即將所有調(diào)諧器的深度增加0.6 mm,二、四極場誤差產(chǎn)生輕微變化,為2.52%和1.34%(圖7(e)劃線處)。之后又進行了2次額外的調(diào)諧以驗證調(diào)諧算法的收斂性,結(jié)果顯示誤差在小范圍振蕩,可能是誤差較小的結(jié)果已接近調(diào)諧算法的極限精度。調(diào)諧各步的誤差結(jié)果見圖7(e)。

    5 結(jié)語

    針對高頻率、緊湊型RFQ在調(diào)諧過程中遇到的問題,提出了一種基于響應(yīng)矩陣和最小二乘法的新調(diào)諧算法,可以在調(diào)諧棒可調(diào)節(jié)范圍內(nèi)求出最優(yōu)解,且為二極場分量增加權(quán)重可以高效地減小二極場誤差。在模擬環(huán)境下進行了調(diào)諧實驗,初步驗證了新調(diào)諧算法的可行性,與SVD解法的比較中體現(xiàn)出了新調(diào)諧算法的優(yōu)勢。根據(jù)調(diào)諧算法編寫了Matlab程序,基于Labview開發(fā)了RFQ的測量和調(diào)諧平臺,通過RFQ樣機的調(diào)諧實驗進一步驗證了調(diào)諧算法。這種調(diào)諧算法是針對高頻率、緊湊型RFQ開發(fā)的,同時也具備很強的普適性,適用于其他頻率的RFQ,將來可以促使緊湊型RFQ應(yīng)用在更多質(zhì)子加速器裝置上,推動醫(yī)用質(zhì)子裝置的普及。

    作者貢獻聲明陸羿行負(fù)責(zé)研究的提出及設(shè)計、數(shù)據(jù)的收集和整理、文章的起草和最終版的修訂;方文程負(fù)責(zé)項目的監(jiān)督和管理、文章的修訂;張俊強負(fù)責(zé)實驗數(shù)據(jù)的收集;郭玉森負(fù)責(zé)實驗數(shù)據(jù)的收集;趙振堂負(fù)責(zé)項目的監(jiān)督和管理、文章的修訂。

    猜你喜歡
    調(diào)諧器四極腔體
    黃河文化龍頭城市的“四極”支撐論析
    高鐵復(fù)雜腔體鑄造數(shù)值仿真及控制技術(shù)研究
    高鐵制動系統(tǒng)復(fù)雜腔體鑄造成形數(shù)值模擬
    橡膠擠出裝置
    一種波蕩器段間四極鐵遠(yuǎn)程精密調(diào)節(jié)機構(gòu)的研制
    核技術(shù)(2016年4期)2016-08-22 09:05:20
    開孔金屬腔體場強增強效應(yīng)分析
    SPE凈化-氣相色譜-三重四極桿質(zhì)譜對油脂中多環(huán)芳烴的測定
    芯科新一代硅電視調(diào)諧器Si21x7面世
    同位素稀釋-氣相色譜-三重四極桿串聯(lián)質(zhì)譜法測定植物油中的多環(huán)芳烴
    Silicon Labs針對中國廣闊的電視市場推出新一代電視調(diào)諧器IC
    国产亚洲欧美98| 久久久精品国产亚洲av高清涩受| av网站在线播放免费| 热re99久久国产66热| 亚洲欧美激情在线| 宅男免费午夜| 高清毛片免费观看视频网站 | 亚洲在线自拍视频| 一级黄色大片毛片| 欧美午夜高清在线| 亚洲五月婷婷丁香| 久久久久久久久免费视频了| 午夜两性在线视频| 亚洲一卡2卡3卡4卡5卡精品中文| 看免费av毛片| 久久婷婷成人综合色麻豆| 免费av毛片视频| 交换朋友夫妻互换小说| 大码成人一级视频| 99国产精品一区二区三区| 国产男靠女视频免费网站| 亚洲中文av在线| 麻豆国产av国片精品| 欧美大码av| 精品福利观看| 国产一卡二卡三卡精品| 成人手机av| 搡老乐熟女国产| 欧美精品亚洲一区二区| 亚洲,欧美精品.| 日韩av在线大香蕉| 成人18禁在线播放| 国产视频一区二区在线看| 久热这里只有精品99| 高清黄色对白视频在线免费看| 亚洲成国产人片在线观看| 91九色精品人成在线观看| 91字幕亚洲| 女性生殖器流出的白浆| 日韩有码中文字幕| 久久精品国产综合久久久| 老汉色∧v一级毛片| 99精品欧美一区二区三区四区| 国产成人影院久久av| 国产片内射在线| 一级a爱视频在线免费观看| 一区福利在线观看| 18禁裸乳无遮挡免费网站照片 | 18美女黄网站色大片免费观看| 亚洲精品国产一区二区精华液| 欧美日韩瑟瑟在线播放| 午夜福利一区二区在线看| 电影成人av| 亚洲avbb在线观看| 国产精品一区二区精品视频观看| av视频免费观看在线观看| 老鸭窝网址在线观看| 熟女少妇亚洲综合色aaa.| 欧美精品啪啪一区二区三区| 他把我摸到了高潮在线观看| 亚洲色图 男人天堂 中文字幕| 9色porny在线观看| 老熟妇仑乱视频hdxx| 久久久水蜜桃国产精品网| 欧美日韩一级在线毛片| 91在线观看av| av欧美777| 久久久久久久久久久久大奶| 精品第一国产精品| 久久人人97超碰香蕉20202| 99国产精品免费福利视频| 91国产中文字幕| 免费在线观看黄色视频的| 人人澡人人妻人| 搡老熟女国产l中国老女人| 国产99白浆流出| 男女床上黄色一级片免费看| 成人黄色视频免费在线看| 怎么达到女性高潮| av中文乱码字幕在线| 成在线人永久免费视频| 中出人妻视频一区二区| 成年版毛片免费区| 一进一出好大好爽视频| 欧美日韩福利视频一区二区| 中文字幕高清在线视频| 亚洲国产精品sss在线观看 | 欧美国产精品va在线观看不卡| 日本欧美视频一区| 日韩大码丰满熟妇| 日韩大尺度精品在线看网址 | 国产欧美日韩精品亚洲av| 午夜精品久久久久久毛片777| 日日爽夜夜爽网站| 国产在线观看jvid| 18美女黄网站色大片免费观看| 久久影院123| 精品一区二区三区视频在线观看免费 | 免费在线观看完整版高清| 成年人黄色毛片网站| 人妻丰满熟妇av一区二区三区| 老汉色∧v一级毛片| 视频区欧美日本亚洲| 色婷婷av一区二区三区视频| 国产精品国产高清国产av| 精品一品国产午夜福利视频| 久久精品91无色码中文字幕| 日日夜夜操网爽| 韩国精品一区二区三区| 琪琪午夜伦伦电影理论片6080| 波多野结衣一区麻豆| 91大片在线观看| 麻豆成人av在线观看| 久久天堂一区二区三区四区| 欧美日韩视频精品一区| 99精品在免费线老司机午夜| 国产av在哪里看| 精品久久久久久久毛片微露脸| 很黄的视频免费| 亚洲av片天天在线观看| 男女下面进入的视频免费午夜 | 午夜日韩欧美国产| 午夜免费成人在线视频| 亚洲狠狠婷婷综合久久图片| 免费在线观看亚洲国产| 欧美乱码精品一区二区三区| 涩涩av久久男人的天堂| 午夜亚洲福利在线播放| 高清在线国产一区| 亚洲三区欧美一区| 成年人免费黄色播放视频| 免费在线观看黄色视频的| 国产成人免费无遮挡视频| 好男人电影高清在线观看| 亚洲精品国产精品久久久不卡| 成人手机av| 精品国产乱子伦一区二区三区| 欧美日韩黄片免| 老汉色av国产亚洲站长工具| 国产区一区二久久| 黄色成人免费大全| 精品乱码久久久久久99久播| 国产国语露脸激情在线看| 看黄色毛片网站| 午夜视频精品福利| videosex国产| 国产97色在线日韩免费| 色老头精品视频在线观看| 国产精品国产av在线观看| av超薄肉色丝袜交足视频| 精品日产1卡2卡| 午夜激情av网站| 国产蜜桃级精品一区二区三区| 国产精品久久视频播放| 精品国产乱子伦一区二区三区| 国产精品综合久久久久久久免费 | 欧美在线黄色| 啪啪无遮挡十八禁网站| 50天的宝宝边吃奶边哭怎么回事| 欧美日本亚洲视频在线播放| 满18在线观看网站| 高潮久久久久久久久久久不卡| 少妇粗大呻吟视频| 日本免费a在线| 99久久综合精品五月天人人| 又紧又爽又黄一区二区| 国产91精品成人一区二区三区| 大香蕉久久成人网| 热re99久久国产66热| 精品人妻1区二区| 18禁美女被吸乳视频| 高清av免费在线| 国产亚洲av高清不卡| 在线免费观看的www视频| 久久久久九九精品影院| 99久久99久久久精品蜜桃| 日韩成人在线观看一区二区三区| cao死你这个sao货| 亚洲欧美日韩高清在线视频| 精品人妻1区二区| 国产有黄有色有爽视频| 天堂影院成人在线观看| 乱人伦中国视频| 91av网站免费观看| 久久人人爽av亚洲精品天堂| 国产亚洲精品第一综合不卡| 国产精品久久久久久人妻精品电影| 国产av一区在线观看免费| 妹子高潮喷水视频| 黄频高清免费视频| 国产在线精品亚洲第一网站| 久久人人精品亚洲av| 国产精品久久视频播放| 成人精品一区二区免费| 亚洲av美国av| av网站在线播放免费| 黄色 视频免费看| 中文字幕人妻丝袜制服| 久久久久久久久免费视频了| 亚洲精品久久成人aⅴ小说| 别揉我奶头~嗯~啊~动态视频| 超碰成人久久| 91成年电影在线观看| √禁漫天堂资源中文www| 少妇粗大呻吟视频| 18禁裸乳无遮挡免费网站照片 | 可以免费在线观看a视频的电影网站| 欧美成人午夜精品| 午夜福利在线免费观看网站| 国产成人啪精品午夜网站| 悠悠久久av| 天堂中文最新版在线下载| av片东京热男人的天堂| 九色亚洲精品在线播放| 男女之事视频高清在线观看| 看黄色毛片网站| 国产免费av片在线观看野外av| 国产成人一区二区三区免费视频网站| 中文亚洲av片在线观看爽| av在线天堂中文字幕 | 巨乳人妻的诱惑在线观看| 国产av又大| 亚洲自偷自拍图片 自拍| 精品乱码久久久久久99久播| 性欧美人与动物交配| 热re99久久精品国产66热6| 欧美日韩亚洲综合一区二区三区_| 久久国产精品人妻蜜桃| 在线观看免费日韩欧美大片| 又黄又粗又硬又大视频| 久久久水蜜桃国产精品网| 成人三级做爰电影| 嫩草影视91久久| 丝袜人妻中文字幕| a级毛片黄视频| 色综合站精品国产| 夜夜躁狠狠躁天天躁| 免费少妇av软件| 欧美色视频一区免费| 色婷婷久久久亚洲欧美| 国产av精品麻豆| 9色porny在线观看| avwww免费| 三级毛片av免费| 亚洲熟女毛片儿| 午夜亚洲福利在线播放| 色综合婷婷激情| 一进一出好大好爽视频| 大码成人一级视频| 午夜影院日韩av| 91国产中文字幕| 亚洲午夜精品一区,二区,三区| 欧美日韩中文字幕国产精品一区二区三区 | 无人区码免费观看不卡| 老司机靠b影院| 亚洲五月天丁香| 看片在线看免费视频| 午夜福利一区二区在线看| 两人在一起打扑克的视频| 成人特级黄色片久久久久久久| 一夜夜www| 美女福利国产在线| 国产黄a三级三级三级人| www.熟女人妻精品国产| 久久狼人影院| 老司机午夜福利在线观看视频| 久久久久久人人人人人| 欧美日韩视频精品一区| 亚洲 欧美 日韩 在线 免费| 久久久久久久久久久久大奶| 久久久久国产精品人妻aⅴ院| 女性生殖器流出的白浆| 日日夜夜操网爽| 久久天堂一区二区三区四区| 免费av中文字幕在线| 亚洲成人精品中文字幕电影 | 国产亚洲精品久久久久久毛片| 18禁美女被吸乳视频| 99国产精品99久久久久| 亚洲精品国产精品久久久不卡| 国产av一区在线观看免费| 国产成人系列免费观看| 亚洲成a人片在线一区二区| 超碰97精品在线观看| a级毛片在线看网站| 高清欧美精品videossex| 国产黄a三级三级三级人| 久久午夜综合久久蜜桃| 欧美午夜高清在线| 黑人猛操日本美女一级片| 91成人精品电影| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美亚洲日本最大视频资源| 国产色视频综合| 欧美日韩av久久| 中文字幕另类日韩欧美亚洲嫩草| 在线国产一区二区在线| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲欧美精品综合久久99| 亚洲全国av大片| 成年人黄色毛片网站| 侵犯人妻中文字幕一二三四区| 少妇裸体淫交视频免费看高清 | 国产午夜精品久久久久久| 高清在线国产一区| 在线永久观看黄色视频| 一级片'在线观看视频| 国产三级在线视频| 精品无人区乱码1区二区| 久久婷婷成人综合色麻豆| 国产成人影院久久av| 97人妻天天添夜夜摸| 性色av乱码一区二区三区2| 色播在线永久视频| 久久久久久大精品| 久久久久久久久中文| 91精品三级在线观看| 国产高清国产精品国产三级| 婷婷六月久久综合丁香| 色尼玛亚洲综合影院| 精品一区二区三卡| 大码成人一级视频| 精品一区二区三卡| 美女扒开内裤让男人捅视频| 国产xxxxx性猛交| 在线av久久热| 亚洲成人免费av在线播放| bbb黄色大片| 丝袜在线中文字幕| 69精品国产乱码久久久| 超碰97精品在线观看| 十八禁人妻一区二区| 97人妻天天添夜夜摸| 满18在线观看网站| 欧美日韩av久久| 高清黄色对白视频在线免费看| 色综合婷婷激情| 国产精品日韩av在线免费观看 | 亚洲 欧美 日韩 在线 免费| 久久性视频一级片| 色哟哟哟哟哟哟| 黄片大片在线免费观看| 丰满饥渴人妻一区二区三| 亚洲人成电影观看| 精品国产亚洲在线| tocl精华| 丝袜在线中文字幕| 亚洲色图 男人天堂 中文字幕| 激情视频va一区二区三区| a在线观看视频网站| 国产精品免费视频内射| 夜夜爽天天搞| 黑人巨大精品欧美一区二区mp4| 人成视频在线观看免费观看| 免费日韩欧美在线观看| 日日干狠狠操夜夜爽| 亚洲自偷自拍图片 自拍| 中文字幕另类日韩欧美亚洲嫩草| 99国产精品一区二区三区| 黄频高清免费视频| 美女国产高潮福利片在线看| av超薄肉色丝袜交足视频| 欧美日韩精品网址| 18美女黄网站色大片免费观看| 国产精品一区二区精品视频观看| 又大又爽又粗| 丰满饥渴人妻一区二区三| 亚洲五月婷婷丁香| 老司机靠b影院| 久久香蕉激情| 国产深夜福利视频在线观看| 桃色一区二区三区在线观看| 人人妻,人人澡人人爽秒播| 日本wwww免费看| 男人舔女人的私密视频| bbb黄色大片| 国产深夜福利视频在线观看| 欧美丝袜亚洲另类 | 99国产综合亚洲精品| 搡老熟女国产l中国老女人| 极品教师在线免费播放| 免费日韩欧美在线观看| 夜夜爽天天搞| 黄网站色视频无遮挡免费观看| 亚洲人成电影观看| 激情视频va一区二区三区| 国产色视频综合| 欧美亚洲日本最大视频资源| 99香蕉大伊视频| 老熟妇仑乱视频hdxx| 女警被强在线播放| 国产三级在线视频| 久久久国产一区二区| 国产高清激情床上av| www.精华液| 精品人妻在线不人妻| 多毛熟女@视频| 免费av毛片视频| 精品第一国产精品| 亚洲精品国产精品久久久不卡| 真人做人爱边吃奶动态| 国产亚洲欧美在线一区二区| 女人爽到高潮嗷嗷叫在线视频| 免费观看人在逋| 精品国产超薄肉色丝袜足j| 男人的好看免费观看在线视频 | 亚洲免费av在线视频| 亚洲欧美日韩无卡精品| 国产亚洲欧美精品永久| 少妇被粗大的猛进出69影院| 亚洲情色 制服丝袜| 欧美日本中文国产一区发布| 91在线观看av| 99国产精品一区二区蜜桃av| 少妇 在线观看| 日韩高清综合在线| 老司机靠b影院| 久久香蕉激情| 亚洲精品一区av在线观看| 久久久国产成人精品二区 | 精品国产乱子伦一区二区三区| 黄色视频不卡| 精品一区二区三区av网在线观看| 男女下面进入的视频免费午夜 | 波多野结衣高清无吗| 电影成人av| 不卡av一区二区三区| 亚洲狠狠婷婷综合久久图片| 亚洲 欧美一区二区三区| 久热这里只有精品99| www国产在线视频色| 亚洲自拍偷在线| 国产成人一区二区三区免费视频网站| 久久久久久免费高清国产稀缺| 这个男人来自地球电影免费观看| 激情在线观看视频在线高清| 叶爱在线成人免费视频播放| 中文字幕高清在线视频| 天堂√8在线中文| 久久人人爽av亚洲精品天堂| 久久精品国产清高在天天线| 色婷婷av一区二区三区视频| 男人舔女人的私密视频| 后天国语完整版免费观看| 亚洲精品一卡2卡三卡4卡5卡| 天堂中文最新版在线下载| 亚洲中文日韩欧美视频| 久久久水蜜桃国产精品网| 9色porny在线观看| 欧美黑人精品巨大| 女人精品久久久久毛片| 久久香蕉精品热| 国产伦人伦偷精品视频| 精品卡一卡二卡四卡免费| 亚洲男人天堂网一区| 国内久久婷婷六月综合欲色啪| 成人三级黄色视频| 精品一区二区三区四区五区乱码| √禁漫天堂资源中文www| 久久亚洲精品不卡| 久久热在线av| 国产免费现黄频在线看| 又黄又爽又免费观看的视频| 国产精品成人在线| www.自偷自拍.com| 亚洲欧美精品综合久久99| 精品国产美女av久久久久小说| 免费观看精品视频网站| 久久中文字幕一级| 侵犯人妻中文字幕一二三四区| 久9热在线精品视频| 国产男靠女视频免费网站| 亚洲男人的天堂狠狠| 午夜福利,免费看| 成人精品一区二区免费| cao死你这个sao货| 老汉色∧v一级毛片| 操出白浆在线播放| 少妇粗大呻吟视频| 久9热在线精品视频| 超碰97精品在线观看| 大香蕉久久成人网| 夜夜看夜夜爽夜夜摸 | 久久久国产一区二区| 无遮挡黄片免费观看| 色婷婷av一区二区三区视频| 久久人人97超碰香蕉20202| 在线观看日韩欧美| 妹子高潮喷水视频| 男人的好看免费观看在线视频 | 18禁国产床啪视频网站| 男女高潮啪啪啪动态图| 久久亚洲精品不卡| 免费久久久久久久精品成人欧美视频| 美女午夜性视频免费| 可以免费在线观看a视频的电影网站| 午夜影院日韩av| 夜夜夜夜夜久久久久| 婷婷精品国产亚洲av在线| 很黄的视频免费| 99久久人妻综合| x7x7x7水蜜桃| 亚洲精品一区av在线观看| √禁漫天堂资源中文www| 成年人免费黄色播放视频| 不卡av一区二区三区| 国产精品偷伦视频观看了| 午夜免费成人在线视频| 国产高清videossex| 琪琪午夜伦伦电影理论片6080| 岛国在线观看网站| 国产熟女午夜一区二区三区| 久久人人精品亚洲av| av天堂久久9| 嫩草影院精品99| 久久国产精品影院| 精品无人区乱码1区二区| 久久精品国产亚洲av香蕉五月| 一级毛片女人18水好多| 黑丝袜美女国产一区| 久久国产精品影院| 极品人妻少妇av视频| 国产高清激情床上av| 久久国产亚洲av麻豆专区| 人成视频在线观看免费观看| 日韩免费av在线播放| 校园春色视频在线观看| 久久香蕉国产精品| 水蜜桃什么品种好| 人人妻人人澡人人看| 久久精品aⅴ一区二区三区四区| 人妻丰满熟妇av一区二区三区| 国产伦人伦偷精品视频| 欧美亚洲日本最大视频资源| 婷婷丁香在线五月| 国产欧美日韩精品亚洲av| 在线十欧美十亚洲十日本专区| 男人操女人黄网站| av在线播放免费不卡| 日日爽夜夜爽网站| 亚洲五月天丁香| 久久欧美精品欧美久久欧美| 天堂影院成人在线观看| 久久人妻福利社区极品人妻图片| 精品一区二区三区视频在线观看免费 | 男男h啪啪无遮挡| a级毛片在线看网站| 97碰自拍视频| 成人影院久久| 中文字幕另类日韩欧美亚洲嫩草| 欧美日韩黄片免| 自拍欧美九色日韩亚洲蝌蚪91| 欧美黄色片欧美黄色片| 国产精品久久久久成人av| 成人黄色视频免费在线看| 国产在线观看jvid| 国内久久婷婷六月综合欲色啪| 亚洲,欧美精品.| 色综合婷婷激情| 18禁裸乳无遮挡免费网站照片 | 精品欧美一区二区三区在线| 成年人免费黄色播放视频| 成年版毛片免费区| 很黄的视频免费| 在线免费观看的www视频| 免费在线观看视频国产中文字幕亚洲| 香蕉国产在线看| 国产又色又爽无遮挡免费看| 久久 成人 亚洲| 中国美女看黄片| 亚洲美女黄片视频| 亚洲精品久久午夜乱码| av网站免费在线观看视频| 国产深夜福利视频在线观看| 成人18禁高潮啪啪吃奶动态图| 亚洲九九香蕉| 悠悠久久av| svipshipincom国产片| av在线天堂中文字幕 | 一级,二级,三级黄色视频| 女人高潮潮喷娇喘18禁视频| 国产在线精品亚洲第一网站| 最近最新中文字幕大全免费视频| 大型黄色视频在线免费观看| 天堂影院成人在线观看| 久久人妻熟女aⅴ| 欧美日韩亚洲高清精品| 黄色丝袜av网址大全| 中文字幕另类日韩欧美亚洲嫩草| 亚洲国产欧美日韩在线播放| 18禁黄网站禁片午夜丰满| 国产精品国产av在线观看| 又黄又粗又硬又大视频| 在线观看日韩欧美| 一进一出好大好爽视频| 免费在线观看视频国产中文字幕亚洲| 中文字幕人妻熟女乱码| 在线永久观看黄色视频| 国产免费男女视频| 国产熟女午夜一区二区三区| 少妇被粗大的猛进出69影院| 欧美人与性动交α欧美软件| 国产极品粉嫩免费观看在线| 一区福利在线观看| 好看av亚洲va欧美ⅴa在| 久久性视频一级片| 精品一品国产午夜福利视频| 咕卡用的链子| 后天国语完整版免费观看| 80岁老熟妇乱子伦牲交| 中文字幕精品免费在线观看视频| 亚洲人成电影免费在线| a级片在线免费高清观看视频| 色婷婷久久久亚洲欧美| 精品无人区乱码1区二区|