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

    基于LBM 的鋁微滴斜柱沉積水平偏移研究1)

    2021-11-09 08:46:26任彥霖劉趙淼
    力學(xué)學(xué)報(bào) 2021年6期
    關(guān)鍵詞:斜柱液滴軸線

    任彥霖 劉趙淼,?,2) 逄 燕,? 王 翔

    *(北京工業(yè)大學(xué)材料與制造學(xué)部,北京 100124)

    ?(北京工業(yè)大學(xué)先進(jìn)制造重點(diǎn)實(shí)驗(yàn)室,北京 100124)

    引言

    在傳統(tǒng)增材制造技術(shù)中,微管結(jié)構(gòu)[1-2]、多孔結(jié)構(gòu)[3-4]因支撐材料難以去除而無法生產(chǎn).金屬液滴沉積制造技術(shù)基于逐點(diǎn)成型原理,通過按需噴射方式[5-6],進(jìn)行無支撐斜柱沉積[7],實(shí)現(xiàn)空間網(wǎng)絡(luò)結(jié)構(gòu)的生產(chǎn)[8-9].鋁液滴是增材制造技術(shù)中常用的工質(zhì),具有較好的界面張力和黏度等材料特性,被廣泛應(yīng)用于液滴沉積制造技術(shù)的研究中[10-11].然而,先沉積液滴非受限表面凝固后的球帽形貌會(huì)影響后沉積液滴運(yùn)動(dòng)行為,使液滴實(shí)際沉積位置與其生成位置存在偏差.這種偏差使液滴沉積行為產(chǎn)生的偏移距離,會(huì)降低斜柱沉積形貌精度.精確分析影響水平偏移程度的控制因素,成為斜柱沉積過程亟待解決的關(guān)鍵問題.

    沉積液滴在球帽形凝固表面上的運(yùn)動(dòng)過程,是一個(gè)復(fù)雜的流動(dòng)、傳熱過程,其實(shí)際沉積位置受到?jīng)_擊[12]、震蕩[13]、共同凝固過程[14]等因素影響.Zhang等[7]在液滴連續(xù)沉積實(shí)驗(yàn)中,通過調(diào)整融合率引起實(shí)際沉積位置偏移,并影響斜柱傾角;進(jìn)一步固定融合率,卻依然存在沉積位置偏移,影響傾角的控制精度[15].Graham 等[16]和Dalili 等[17]實(shí)驗(yàn)觀測(cè)了水滴的相互融合過程,表明沉積位置的偏移由界面張力作用引起.Li 等[18]探究了石蠟滴水平連續(xù)沉積過程,發(fā)現(xiàn)液滴間距會(huì)影響位置偏移程度.這些研究多集中于液滴融合過程產(chǎn)生的偏移運(yùn)動(dòng),而缺少金屬液滴凝固形成的球帽結(jié)構(gòu)對(duì)偏移運(yùn)動(dòng)的影響.Ju 等[19]研究了液滴在球面上的沉積運(yùn)動(dòng),表明球面曲率影響液滴沉積運(yùn)動(dòng),但由于液滴與球面同軸,而未產(chǎn)生滑落.Tian 等[20]實(shí)驗(yàn)研究了液滴在傾斜表面的沉積過程,提出了沉積液滴動(dòng)力學(xué)行為的能量分析法.對(duì)于鋁液滴在球帽形凝固表面上的偏移運(yùn)動(dòng)過程,前、后液滴間軸線距離的影響及其受力機(jī)理尚未明確,對(duì)斜柱沉積過程參數(shù)的優(yōu)化設(shè)計(jì)有待研究.

    為分析偏移運(yùn)動(dòng)和受力機(jī)理,鋁液滴的連續(xù)沉積過程需要采用格子玻爾茲曼方法(LBM),進(jìn)行數(shù)值模擬研究.LBM 是一種介觀尺度的模擬算法,在兩相流動(dòng)方面采用S-C 偽勢(shì)模型統(tǒng)一處理離散相和分散相[21],并引入C-S 狀態(tài)方程計(jì)算分子碰撞有效質(zhì)量使其滿足熱力學(xué)一致性[22],對(duì)平衡速度進(jìn)行修正[23];在傳熱相變方面采用焓的顯式法求解傳熱過程,通過雙分布方程使其與流動(dòng)耦合[24],并對(duì)固液相變前沿進(jìn)行無滑移邊界處理[25-26].LBM 已經(jīng)具備解決實(shí)際工程中傳熱、相變問題的能力[27-31],并應(yīng)用于單液滴沉積過程的研究中[32-34],有助于實(shí)現(xiàn)金屬液滴連續(xù)沉積過程的建模.

    本文將基于LBM 建立鋁液滴斜柱沉積模型,采用自編程序研究液滴水平偏移過程和受力機(jī)理.根據(jù)液滴能量演化趨勢(shì),對(duì)沉積過程進(jìn)行階段劃分.結(jié)合質(zhì)心位移曲線,分析水平偏移產(chǎn)生的主要階段和主要推動(dòng)力.根據(jù)偏移距離的演化規(guī)律,對(duì)掃描步距進(jìn)行優(yōu)化,以實(shí)現(xiàn)斜柱的均勻沉積和傾角的精確控制.

    1 物理模型

    斜柱沉積制造是多個(gè)液滴按相同間距連續(xù)沉積的過程,以其中兩液滴先后下落的物理過程,作為主要研究對(duì)象.如圖1 所示,密度、初始焓、直徑分別為ρH,HA,Dd的兩個(gè)鋁液滴,從同一高度h,水平相距L的位置先后下落.先沉積的液滴接觸溫度為TB的恒溫基板后凝固,后沉積的液滴接觸前者表面后凝固.外界環(huán)境氣體密度為ρ0,初始溫度為T0(T0=TB).鋁液滴材料參數(shù)如表1 所示.

    圖1 多液滴沉積物理過程簡圖Fig.1 Physical process of multi-droplet deposition

    表1 鋁液滴材料參數(shù)Table 1 Material properties of the aluminum droplets

    2 數(shù)值方法

    2.1 兩相流動(dòng)方程

    液滴沉積過程是一個(gè)兩相流動(dòng)過程,需要準(zhǔn)確追蹤氣液界面.本文采用偽勢(shì)模型,通過流體粒子間相互作用力,使其自發(fā)地形成兩相界面,以滿足熱力學(xué)一致性和第一性原理.首先對(duì)流動(dòng)過程建立密度的格子玻爾茲曼方程

    式中,i表示格子方向的索引編號(hào),ei表示i方向上的單位格子速度,τf表示密度松弛時(shí)間[35].浸沒邊界法是LBM 中處理曲面邊界的常用方法[36].本文采用Noble 等[26]提出的方法處理相變界面,以B表示固體邊界判別函數(shù),表示附加碰撞項(xiàng).

    在流動(dòng)方程的基礎(chǔ)上,引入相鄰位點(diǎn)間粒子的長程相互作用力,實(shí)現(xiàn)流體的兩相分離

    式中,G表示相互作用強(qiáng)度,ψ 表示碰撞有效質(zhì)量.為了保證兩相分離狀態(tài)具有熱力學(xué)一致性,碰撞有效質(zhì)量通過C-S 狀態(tài)方程計(jì)算[22]

    在流固界面上,近壁面粒子所受作用力為

    式中,ρw表示壁面粒子假想密度,用于調(diào)節(jié)固液表面浸潤度.采用Shan-Doolen 平衡態(tài)分布函數(shù)的速度校正法,對(duì)平衡速度和物理速度進(jìn)行修正[23]

    2.2 傳熱相變方程

    液滴在沉積過程中傳熱并發(fā)生凝固,需要引入焓的格子玻爾茲曼方程求解

    式中,τg表示焓的松弛時(shí)間,表示焓的平衡分布函數(shù),并采用Huang 等[34]提出的顯式方程表示

    式中,ωi表示i方向上的加權(quán),cs表示格子聲速,Cp表示定壓熱容.根據(jù)流體焓不僅可以計(jì)算出當(dāng)?shù)販囟?還可以追蹤固液相變界面

    式中,HS表示流體固化焓,HL表示流體液化焓,TS表示固相線,TL表示液相線.對(duì)于相變界面上的無滑移邊界條件,采用Noble 等[26]提出的移動(dòng)固體邊界處理方法,并通過固體判別函數(shù)和附加碰撞項(xiàng)修正流動(dòng)方程

    式中,k表示運(yùn)動(dòng)方向與ei相反的分布函數(shù)索引,US表示固體的運(yùn)動(dòng)速度.

    2.3 網(wǎng)格獨(dú)立性驗(yàn)證

    通過模擬單液滴下落過程中的形心高度變化,對(duì)模型進(jìn)行網(wǎng)格獨(dú)立性驗(yàn)證.計(jì)算域上下表面設(shè)為半步反彈及恒溫邊界條件,四周設(shè)置為循環(huán)邊界條件,模擬單液滴下落過程.如圖2 所示,當(dāng)液滴直徑為Dd=36 lu 時(shí),計(jì)算域總網(wǎng)格量超過7.9×105個(gè)單元,液滴形心運(yùn)動(dòng)趨勢(shì)不再有明顯變化.在LBM中,計(jì)算域在無特殊格式處理時(shí),通常采用均勻網(wǎng)格劃分,網(wǎng)格精度由宏觀尺寸在格子上的分辨率體現(xiàn).根據(jù)液滴實(shí)際直徑Dd=0.6 mm,和網(wǎng)格獨(dú)立性驗(yàn)證得到的網(wǎng)格量大小,得到宏觀空間尺度的國際單位分辨率為1 m=6.0×104lu,并在保證計(jì)算收斂性的基礎(chǔ)上,進(jìn)一步確定時(shí)間、溫度、質(zhì)量等基本單位在格子上的分辨率(如表2 所示).

    圖2 網(wǎng)格獨(dú)立性驗(yàn)證Fig.2 Grid independence verification

    表2 國際基本單位與格子單位的換算關(guān)系Table 2 Relationship between the SI units and the lattice units

    3 結(jié)果與討論

    對(duì)水平偏移距離的準(zhǔn)確控制,需要根據(jù)液滴沉積過程中的能量變化,對(duì)液滴運(yùn)動(dòng)過程進(jìn)行劃分,探尋水平偏移的主要產(chǎn)生階段和主要推動(dòng)力.

    3.1 能量演化及運(yùn)動(dòng)階段劃分

    沉積液滴與環(huán)境氣體、凝固液滴形成三相接觸,其表面能直接反應(yīng)液滴形態(tài)變化和運(yùn)動(dòng)過程.Tian等[20]根據(jù)表面能的充放過程,對(duì)液滴運(yùn)動(dòng)過程進(jìn)行階段劃分,并建立了表面能計(jì)算公式

    式中,Sgl表示氣液界面面積,Sls表示固液界面面積,θe表示液滴靜態(tài)接觸角.沉積液滴除了表面能外,還具有重力勢(shì)能GP和動(dòng)能KE,這兩種能量的計(jì)算模型分別表示為

    式中,ρ 表示流體微元密度,g表示重力加速度,Δz表示流體微元運(yùn)動(dòng)的相對(duì)高度,u表示流體微元速度.沉積液滴接觸凝固表面后,由于壁面剪切作用產(chǎn)生黏性耗散VD,使表面能、重力勢(shì)能與動(dòng)能之和降低.因此,液滴沉積過程中產(chǎn)生的黏性耗散可推算為

    圖3 提取了浸潤度為ξ=0.42、兩液滴軸線距離為L*=0.33 下的液滴沉積能量變化.其中Fo是無量綱時(shí)間尺度,表示為

    式中,υ 是鋁的運(yùn)動(dòng)黏度,t表示時(shí)間,Dd為液滴的特征長度,取液滴的初始直徑.根據(jù)SE的充放變化,可以將沉積液滴運(yùn)動(dòng)過程分為下落階段、快速擴(kuò)張階段、慢速擴(kuò)張階段、回彈階段[37](如圖3).在下落階段(Fo<0.17),液滴在界面張力作用下維持球狀,SE不發(fā)生明顯改變.重力作用使勢(shì)能下降|ΔGP|=2.32 mu·lu2/ts2,并等額轉(zhuǎn)化為動(dòng)能KE,形成接觸凝固表面前的沖擊速度.

    圖3 后沉積液滴運(yùn)動(dòng)過程中的能量演化曲線Fig.3 Tendency of energy in the motion of deposit droplet

    液滴接觸凝固表面后,與氣體、凝固表面形成三相接觸,進(jìn)入快速擴(kuò)張階段(0.17 <Fo≤0.22).表面能降低|ΔS E|=9.77 mu·lu2/ts2,重力勢(shì)能減小|ΔGP|=3.23 mu·lu2/ts2,表明毛細(xì)力和重力在快速擴(kuò)張運(yùn)動(dòng)中起推動(dòng)作用.此時(shí)液滴中央截面狀態(tài)如圖3(b) 所示,毛細(xì)力帶動(dòng)接觸線擴(kuò)張,并使液滴質(zhì)心下降.過流斷面的減小使黏性剪切作用增強(qiáng),使黏性耗散增大|ΔVD|=13.01 mu·lu2/ts2,動(dòng)能減小|ΔKE|=0.01 mu·lu2/ts2.

    在慢速擴(kuò)張階段中(0.22 <Fo≤0.28).表面能增大|ΔS E|=3.09 mu·lu2/ts2,重力勢(shì)能降低|ΔGP|=0.94 mu·lu2/ts2.該結(jié)果表明毛細(xì)力在慢速擴(kuò)張過程中起阻礙作用,而重力起推動(dòng)作用.慢速擴(kuò)張階段的初始狀態(tài)如圖3(c) 所示,接觸線達(dá)到毛細(xì)平衡狀態(tài),并開始在重力和流動(dòng)慣性作用下受迫擴(kuò)張.流速因氣液界面阻礙而降低,使壁面剪切強(qiáng)度降低,黏性耗散僅為|ΔVD|=0.91 mu·lu2/ts2,同時(shí)動(dòng)能降低|ΔKE|=3.06 mu·lu2/ts2.

    進(jìn)入回彈階段后 (Fo> 0.28),表面能下降|ΔS E|=4.44 mu·lu2/ts2,同時(shí)重力勢(shì)能提高|ΔGP|=0.13 mu·lu2/ts2,表明毛細(xì)力在接觸線回彈過程中起推動(dòng)作用、重力起阻礙作用,但二者作用效果均低于快速擴(kuò)張階段.此時(shí)液滴中央截面狀態(tài)如圖3(d)所示,毛細(xì)力帶動(dòng)接觸線回縮,并使質(zhì)心位置升高.過流斷面的增大和流速的減小使黏性耗散增大|ΔVD|=4.43 mu·lu2/ts2,并使動(dòng)能降低|ΔKE|=0.12 mu·lu2/ts2.

    盡管毛細(xì)力在擴(kuò)張運(yùn)動(dòng)中,先后起到推動(dòng)作用和抑制作用,但由于流動(dòng)形式相似,所以下文被統(tǒng)稱為擴(kuò)張階段.圖4(a) 表明擴(kuò)張階段中液滴內(nèi)部流場(chǎng)以慣性力主導(dǎo)的下落運(yùn)動(dòng)為主,而凝固表面上的流場(chǎng)表現(xiàn)為擴(kuò)散運(yùn)動(dòng).如圖4(b)所示的回彈階段中,上游接觸點(diǎn)切線傾角較小,重力與毛細(xì)力同向,促使接觸線進(jìn)行回彈運(yùn)動(dòng);下游接觸點(diǎn)切線傾角較大,重力與毛細(xì)力反向,反而使接觸線繼續(xù)向下游運(yùn)動(dòng).由此導(dǎo)致液滴整體向下游運(yùn)動(dòng),形成了滾動(dòng)流動(dòng)狀態(tài).

    圖4 擴(kuò)張和回彈階段中的流場(chǎng)形式Fig.4 The flow field in spreading and rebound stage

    由于前沉積液滴的凝固表面具有軸對(duì)稱的球面形貌,使后沉積液滴在各方向上的運(yùn)動(dòng)相同.為分析水平偏移運(yùn)動(dòng)過程,后文以偏移方向作為水平x軸正方向,并提取三維結(jié)果中液滴質(zhì)心運(yùn)動(dòng)距離和速度進(jìn)行分析.

    3.2 水平偏移主要階段及推動(dòng)力

    圖5 展示了兩液滴軸線距離為L*=0.33 的偏移運(yùn)動(dòng)過程.其中L*和δL*分別為無量綱軸線間距和無量綱偏移距離,表示為

    圖5 沖擊液滴偏移過程Fig.5 Horizontal displacements of the deposit droplets

    是液滴形心在水平方向的無量綱偏移速度

    沉積液滴在擴(kuò)張階段中發(fā)生加速偏移,并在Fo=0.29 時(shí)偏移速度達(dá)到峰值=1.32,而偏移距離僅為δL*=0.05.在隨后的回彈階段中,偏移運(yùn)動(dòng)開始減速,在Fo=0.52 時(shí)偏移速度降低至=0.06,偏移距離增大至δL*=0.21.由此可見,水平偏移運(yùn)動(dòng)主要發(fā)生在回彈階段,而加速過程發(fā)生在擴(kuò)張階段.根據(jù)3.1 節(jié)對(duì)擴(kuò)張階段的受力分析,表明重力、毛細(xì)力是水平偏移運(yùn)動(dòng)的主要推動(dòng)力.

    3.3 掃描步距優(yōu)化

    為實(shí)現(xiàn)對(duì)偏移距離的準(zhǔn)確控制,需要在實(shí)際工況中,選擇與主要推動(dòng)力對(duì)應(yīng)的控制條件,并對(duì)其影響程度進(jìn)行標(biāo)定.

    鋁液滴與氣體、凝固表面形成的三相接觸過程中,固液界面浸潤度是控制毛細(xì)力作用的主要參數(shù).液滴間軸線距離決定了兩液滴的接觸點(diǎn)位置,而在半球形凝固表面上,接觸點(diǎn)的位置與該點(diǎn)切線斜率直接相關(guān),并決定了切線方向上的重力分量.液滴軸線距離和固液界面浸潤度成為本文偏移距離研究的主要控制參數(shù).

    如圖6 所示,水平偏移距離隨軸線距離的演化趨勢(shì)呈現(xiàn)階段化特性:L*<0.10 時(shí),偏移距離隨著軸線距離的增大而增大;而L*>0.10 后,偏移距離隨著軸線距離的增大而減小.在不同固液浸潤度條件下對(duì)比,表明偏移程度隨著浸潤度的增大而減小,并且不影響偏移程度隨軸線距離的演化趨勢(shì).由于固液浸潤度是材料屬性,需要根據(jù)材料浸潤度以及所需的偏移距離δL,選擇對(duì)應(yīng)的軸線距離L.以浸潤度為ξ=0.42 的材料工況為例,對(duì)偏移距離進(jìn)行數(shù)據(jù)擬合得到對(duì)勾函數(shù)關(guān)系

    圖6 液滴間軸線距離對(duì)水平偏移距離的影響Fig.6 Influence of axes distance between droplets on the horizontal deviation

    偏移程度隨軸線間距演化的階段化特征,表明液滴的偏移運(yùn)動(dòng)過程存在競(jìng)爭機(jī)制.圖7 對(duì)比了不同軸線距離下,液滴水平偏移速度的演化趨勢(shì):在接觸凝固表面后開始加速,進(jìn)入回彈階段后減速,并在表面能完全釋放后速度趨于平緩.演化趨勢(shì)差異,體現(xiàn)在加速段時(shí)長和速度極大值的變化.

    如圖7 所示,隨著軸線間距由L*=0.03 增大到L*=0.14,速度極大值從=0.66 增大到=1.44,隨后小幅減小.液滴接觸凝固表面后(Fo>Fo0),速度達(dá)到極大值前的加速段時(shí)長Foa=Fo-Fo0,隨傾角的增大而縮短.隨著軸線間距由L*=0.03 增大到L*=0.33,加速段時(shí)長從Foa=0.25 減小到Foa=0.17.對(duì)速度極大值和加速段時(shí)長的模擬結(jié)果進(jìn)行數(shù)據(jù)擬合,得到與軸線間距的關(guān)系式

    圖7 軸線間距對(duì)偏移運(yùn)動(dòng)的影響Fig.7 Influence of axis distance on deviate movement

    該結(jié)果表明加速段時(shí)長與軸線間距呈單調(diào)減小關(guān)系,而速度極大值隨軸線間距先增大后小幅減小.二者對(duì)偏移程度起著相反的作用:加速段時(shí)長的縮短會(huì)提高偏移程度,而速度極大值的增大則與之相反.這種競(jìng)爭關(guān)系是導(dǎo)致偏移程度隨軸線間距表現(xiàn)出階段化特征的根本原因.

    作為液滴沉積的控制參數(shù)之一,沉積高度的增大會(huì)提高沖擊速度,也會(huì)使液滴溫度提早降低,所以該工況的研究范圍被控制在h*=1.67~2.00 內(nèi).圖8展示了h*=1.67,1.72,1.83 三組工況的對(duì)比結(jié)果,表明偏移程度隨軸線間距具有相似的演化趨勢(shì),最大偏移距離與沉積高度負(fù)相關(guān),分別為δL*=0.41,0.38,0.35.這是由于液滴接觸凝固表面時(shí),內(nèi)部流體以凝固表面上的擴(kuò)散運(yùn)動(dòng)為主,沖擊速度的增大使流體克服重力作用產(chǎn)生的流量分配差異,使液滴在擴(kuò)張階段中的偏移距離和速度減小.該結(jié)果表明,沉積高度的增大對(duì)偏移程度起抑制作用.

    圖8 軸線間距對(duì)水平偏移程度的影響Fig.8 Influence of axial distance on horizontal deviation

    為實(shí)現(xiàn)斜柱的均勻沉積,后續(xù)液滴沉積時(shí)的掃描步距既要包括液滴間的軸線距離,還要考慮前沉積液滴的偏移距離

    在保證每個(gè)液滴掃描步距一致的基礎(chǔ)上,根據(jù)Zhang 等[7]提出的斜柱傾角理論公式,可以得到

    式中ε 表示液滴間的融合率,是沉積后的液滴間形心距離與初始直徑之比.本文采用固液浸潤度為ξ=0.42 的鋁液滴,從高度為h*=1.5 下落.通過測(cè)量液滴形心間距,得到融合率為ε=0.63,并將其代入式(29)、與式(25)和式(28)聯(lián)立,得到液滴軸線距離與傾角關(guān)系式為

    如圖9 所示,根據(jù)式(28) 計(jì)算得到掃描步距分別為w*=0.21,0.29,0.34,0.37 的工況下,多液滴斜柱沉積模擬結(jié)果具有均勻的形貌結(jié)構(gòu),并且傾角與理論結(jié)果一致[7].

    圖9 斜柱傾角與掃描步距關(guān)系及其沉積形貌Fig.9 Relationship between the tilt angle and the scanning step

    4 結(jié)論

    本文通過建立LBM 液滴連續(xù)沉積模型,探究了鋁液滴沉積過程中水平偏移的主要發(fā)生階段及其影響因素.液滴沉積運(yùn)動(dòng)從能量角度被分為下落、快速擴(kuò)張、慢速擴(kuò)張和回彈4 個(gè)階段.擴(kuò)張階段是偏移加速運(yùn)動(dòng)的主要階段,而回彈階段產(chǎn)生了主要偏移距離.偏移運(yùn)動(dòng)主要推動(dòng)力源于毛細(xì)力、重力和黏性剪切力在擴(kuò)張階段中的共同作用.此外,偏移距離隨軸線距離演化趨勢(shì)呈階段化特性,擬合結(jié)果呈對(duì)勾函數(shù)關(guān)系.這種階段化特性是偏移運(yùn)動(dòng)中加速段時(shí)長和速度極大值間的競(jìng)爭關(guān)系導(dǎo)致.沉積高度和浸潤度的減小均不影響該演化趨勢(shì),但會(huì)使偏移程度增大.結(jié)合斜柱幾何特性,掃描步距與斜柱傾角的關(guān)系式得到建立.在該模型下選擇掃描步距,進(jìn)行的多液滴連續(xù)沉積模擬,具有均勻的斜柱形貌,并且傾角與理論結(jié)果吻合程度較好.本文所得結(jié)果有助于偏移距離的精確控制,提高斜柱沉積制造的形貌質(zhì)量.

    猜你喜歡
    斜柱液滴軸線
    某鋼框架斜柱結(jié)構(gòu)設(shè)計(jì)分析
    超高混凝土獨(dú)立斜柱關(guān)鍵施工技術(shù)的研究
    曲軸線工件劃傷問題改進(jìn)研究
    液滴間相互碰撞融合與破碎的實(shí)驗(yàn)研究
    噴淋液滴在空氣環(huán)境下的運(yùn)動(dòng)特性
    混凝土斜柱施工質(zhì)量控制
    綠色科技(2019年4期)2019-03-20 06:57:14
    基于回歸分析的水電機(jī)組軸線曲折預(yù)判斷分析
    行書章法淺析(十五)書寫應(yīng)把握行軸線
    建筑工程中帶斜柱轉(zhuǎn)換結(jié)構(gòu)設(shè)計(jì)分析
    凸輪軸孔軸線與止推面垂直度超差問題研究
    河南科技(2014年16期)2014-02-27 14:13:21
    欧美激情国产日韩精品一区| 国产精品国产三级专区第一集| 美女大奶头视频| 亚洲18禁久久av| 一级毛片我不卡| 免费看日本二区| 亚洲精华国产精华液的使用体验| 身体一侧抽搐| 天堂√8在线中文| 日韩av在线大香蕉| 免费看日本二区| 色5月婷婷丁香| ponron亚洲| 亚洲av.av天堂| 在线观看免费高清a一片| 国产69精品久久久久777片| 日韩欧美一区视频在线观看 | 亚洲最大成人av| 久久久精品欧美日韩精品| 淫秽高清视频在线观看| 亚洲欧美一区二区三区黑人 | 色5月婷婷丁香| 精品欧美国产一区二区三| 日本午夜av视频| 春色校园在线视频观看| 免费观看在线日韩| 熟妇人妻不卡中文字幕| 国产av国产精品国产| 亚洲精品久久久久久婷婷小说| 男女那种视频在线观看| 老司机影院成人| 日韩三级伦理在线观看| 黄色日韩在线| 国产日韩欧美在线精品| 成人av在线播放网站| 天堂√8在线中文| 日韩一本色道免费dvd| 夫妻性生交免费视频一级片| 人人妻人人看人人澡| 欧美极品一区二区三区四区| 在线 av 中文字幕| 国产av国产精品国产| 国产精品人妻久久久影院| 永久免费av网站大全| 国产探花极品一区二区| 日韩 亚洲 欧美在线| 国内精品宾馆在线| 久久6这里有精品| 亚洲成人精品中文字幕电影| 在线天堂最新版资源| 国产激情偷乱视频一区二区| 99久久人妻综合| 国产老妇伦熟女老妇高清| 综合色丁香网| 99热网站在线观看| 中文字幕av成人在线电影| 搞女人的毛片| 亚洲av在线观看美女高潮| 成人鲁丝片一二三区免费| 亚洲精品影视一区二区三区av| 国产亚洲91精品色在线| 国产免费福利视频在线观看| 一级毛片aaaaaa免费看小| 麻豆精品久久久久久蜜桃| 亚洲婷婷狠狠爱综合网| 国产高清有码在线观看视频| 天堂√8在线中文| 欧美精品国产亚洲| 少妇丰满av| 成年女人看的毛片在线观看| 国产精品.久久久| 18+在线观看网站| 亚洲av电影不卡..在线观看| 性色avwww在线观看| 激情 狠狠 欧美| 91午夜精品亚洲一区二区三区| 亚洲丝袜综合中文字幕| 在线免费观看不下载黄p国产| 国产黄色小视频在线观看| 亚洲国产精品专区欧美| 日本-黄色视频高清免费观看| 伦理电影大哥的女人| 人妻制服诱惑在线中文字幕| 大香蕉久久网| 亚洲成人av在线免费| 欧美最新免费一区二区三区| 免费观看在线日韩| 三级经典国产精品| 国产伦一二天堂av在线观看| 国产又色又爽无遮挡免| 亚洲在线观看片| 午夜爱爱视频在线播放| 亚洲四区av| 久久午夜福利片| 久久97久久精品| 我要看日韩黄色一级片| 久久久久久久久久黄片| 亚洲欧洲国产日韩| 亚洲真实伦在线观看| 欧美日韩综合久久久久久| 成人高潮视频无遮挡免费网站| 小蜜桃在线观看免费完整版高清| 丝瓜视频免费看黄片| 天堂俺去俺来也www色官网 | 国产精品久久视频播放| 69av精品久久久久久| 日本欧美国产在线视频| 亚洲精品成人av观看孕妇| 久久久久九九精品影院| 日韩亚洲欧美综合| 一区二区三区四区激情视频| 免费看美女性在线毛片视频| 大香蕉久久网| 久久久久久久大尺度免费视频| 久久精品夜色国产| 69av精品久久久久久| 久久久精品94久久精品| 免费观看性生交大片5| 成人特级av手机在线观看| 国产亚洲最大av| 精品久久国产蜜桃| 天堂√8在线中文| 九色成人免费人妻av| 国产精品国产三级国产av玫瑰| 只有这里有精品99| 亚洲激情五月婷婷啪啪| 成人毛片a级毛片在线播放| 免费看光身美女| 成人一区二区视频在线观看| 日韩电影二区| videos熟女内射| 又爽又黄a免费视频| 特大巨黑吊av在线直播| 狂野欧美白嫩少妇大欣赏| 午夜免费激情av| 欧美日韩视频高清一区二区三区二| 国产av不卡久久| 亚洲国产精品sss在线观看| 日韩av在线免费看完整版不卡| 日本黄大片高清| 国产高清不卡午夜福利| 91狼人影院| 大又大粗又爽又黄少妇毛片口| 99热这里只有精品一区| 三级国产精品片| 国产69精品久久久久777片| 三级毛片av免费| 成人性生交大片免费视频hd| 国产老妇伦熟女老妇高清| 免费观看性生交大片5| 丝袜美腿在线中文| 天美传媒精品一区二区| 18禁在线播放成人免费| 国产高清国产精品国产三级 | 国产乱来视频区| 99热全是精品| 成年免费大片在线观看| 秋霞在线观看毛片| 亚洲成人一二三区av| 欧美丝袜亚洲另类| 少妇熟女欧美另类| 欧美3d第一页| 少妇熟女欧美另类| 麻豆乱淫一区二区| 又爽又黄无遮挡网站| 精品久久久久久成人av| 亚洲精品乱码久久久v下载方式| 亚洲美女搞黄在线观看| 直男gayav资源| 亚洲成人一二三区av| 99久国产av精品| 大片免费播放器 马上看| 日韩欧美 国产精品| 午夜福利网站1000一区二区三区| 永久免费av网站大全| 日本午夜av视频| 国产亚洲午夜精品一区二区久久 | 男人舔奶头视频| 欧美最新免费一区二区三区| 国产午夜精品论理片| 亚洲国产日韩欧美精品在线观看| 波野结衣二区三区在线| 1000部很黄的大片| 国产三级在线视频| 精品久久久久久电影网| 国产一区二区亚洲精品在线观看| 一级毛片久久久久久久久女| 日本爱情动作片www.在线观看| 观看美女的网站| 欧美xxxx黑人xx丫x性爽| 精品久久久久久久末码| 午夜精品在线福利| 国产麻豆成人av免费视频| 久久97久久精品| 卡戴珊不雅视频在线播放| 国产老妇伦熟女老妇高清| 狠狠精品人妻久久久久久综合| 亚洲av成人精品一区久久| 特级一级黄色大片| 人体艺术视频欧美日本| 亚洲av中文av极速乱| 黄片无遮挡物在线观看| 少妇被粗大猛烈的视频| 高清日韩中文字幕在线| 日本猛色少妇xxxxx猛交久久| 久久久精品欧美日韩精品| 能在线免费看毛片的网站| 三级毛片av免费| 国产乱人偷精品视频| 青春草视频在线免费观看| 中文资源天堂在线| 日韩,欧美,国产一区二区三区| 国产精品一区二区三区四区免费观看| 国产精品久久久久久久电影| 中文字幕av成人在线电影| eeuss影院久久| 国产一区亚洲一区在线观看| 日韩av免费高清视频| 嫩草影院精品99| av在线播放精品| 毛片女人毛片| 人妻系列 视频| 日本三级黄在线观看| 国产高清三级在线| 亚洲成人一二三区av| 日本三级黄在线观看| 欧美+日韩+精品| 午夜福利在线观看免费完整高清在| 亚洲乱码一区二区免费版| 国产精品国产三级国产av玫瑰| 成人漫画全彩无遮挡| 97热精品久久久久久| 青春草国产在线视频| 国产大屁股一区二区在线视频| 亚洲av中文字字幕乱码综合| 91精品一卡2卡3卡4卡| 黄色配什么色好看| 免费观看av网站的网址| 男的添女的下面高潮视频| 又大又黄又爽视频免费| 亚洲国产精品成人综合色| 免费看日本二区| 免费看不卡的av| av天堂中文字幕网| 人人妻人人澡欧美一区二区| 高清日韩中文字幕在线| 日本免费a在线| 欧美丝袜亚洲另类| 免费看不卡的av| 大话2 男鬼变身卡| 嫩草影院精品99| 国产亚洲精品久久久com| 午夜精品在线福利| 亚洲美女搞黄在线观看| 中文字幕av在线有码专区| 免费看光身美女| 国产老妇女一区| 91aial.com中文字幕在线观看| 99久久精品热视频| 久久久精品欧美日韩精品| 丝袜美腿在线中文| 最近2019中文字幕mv第一页| 亚洲av电影不卡..在线观看| 国产淫语在线视频| 国产精品一区二区三区四区免费观看| 身体一侧抽搐| 国产黄频视频在线观看| 亚洲精品乱久久久久久| 一级二级三级毛片免费看| 国内精品宾馆在线| 狠狠精品人妻久久久久久综合| 国产激情偷乱视频一区二区| 99视频精品全部免费 在线| 午夜日本视频在线| 亚洲av.av天堂| 一级二级三级毛片免费看| 亚洲国产色片| 熟妇人妻不卡中文字幕| 国产一级毛片七仙女欲春2| 久久精品综合一区二区三区| 国产69精品久久久久777片| 国产激情偷乱视频一区二区| 中文天堂在线官网| 在线免费观看的www视频| 久久午夜福利片| 18禁动态无遮挡网站| 老司机影院毛片| 亚洲欧洲日产国产| 久久久久性生活片| 老师上课跳d突然被开到最大视频| 欧美高清成人免费视频www| 大陆偷拍与自拍| 亚洲乱码一区二区免费版| 亚洲国产欧美在线一区| 国产黄片视频在线免费观看| 国产精品一区www在线观看| av线在线观看网站| 欧美日本视频| 伦理电影大哥的女人| 日本熟妇午夜| 国产精品福利在线免费观看| 国产成人福利小说| 亚洲精品国产av成人精品| 建设人人有责人人尽责人人享有的 | 色综合色国产| 国产有黄有色有爽视频| 亚洲精品影视一区二区三区av| 免费观看的影片在线观看| 蜜桃亚洲精品一区二区三区| 人人妻人人看人人澡| 三级国产精品欧美在线观看| 亚洲成人中文字幕在线播放| 久久99热这里只频精品6学生| 亚洲精品第二区| 精品一区二区三卡| 成年版毛片免费区| 亚洲欧美日韩无卡精品| 久久久久久国产a免费观看| 一本久久精品| 亚洲美女搞黄在线观看| 搞女人的毛片| 街头女战士在线观看网站| 午夜激情久久久久久久| 午夜免费激情av| 精品人妻一区二区三区麻豆| 黄色一级大片看看| 亚洲激情五月婷婷啪啪| 亚洲18禁久久av| 一级毛片我不卡| 亚洲精品一区蜜桃| 国产精品1区2区在线观看.| 欧美日韩在线观看h| 一级二级三级毛片免费看| 午夜激情福利司机影院| 久久这里有精品视频免费| 身体一侧抽搐| av又黄又爽大尺度在线免费看| 一级毛片我不卡| 国产乱人偷精品视频| 亚洲色图av天堂| 亚洲精品成人久久久久久| 99九九线精品视频在线观看视频| 麻豆精品久久久久久蜜桃| 国产在线男女| 午夜福利在线观看免费完整高清在| 欧美bdsm另类| 欧美日韩视频高清一区二区三区二| 麻豆乱淫一区二区| 亚洲精品一二三| 日日啪夜夜撸| 搡老乐熟女国产| 久久精品国产亚洲网站| 日韩欧美一区视频在线观看 | 高清视频免费观看一区二区 | 男人爽女人下面视频在线观看| 99热这里只有是精品50| 国产老妇女一区| 国产精品一及| 春色校园在线视频观看| 国产大屁股一区二区在线视频| 欧美日韩综合久久久久久| 最近最新中文字幕免费大全7| 美女内射精品一级片tv| 国产在线一区二区三区精| 大话2 男鬼变身卡| 色吧在线观看| 丰满少妇做爰视频| 乱码一卡2卡4卡精品| 熟女电影av网| 亚州av有码| 久久久色成人| 夜夜爽夜夜爽视频| 日本黄大片高清| 亚洲综合色惰| 免费观看的影片在线观看| 亚洲av成人精品一区久久| 亚洲欧美一区二区三区黑人 | 丝袜喷水一区| 九九在线视频观看精品| 热99在线观看视频| 国产老妇伦熟女老妇高清| 熟妇人妻久久中文字幕3abv| 国产 亚洲一区二区三区 | 两个人视频免费观看高清| 精品国内亚洲2022精品成人| 午夜福利网站1000一区二区三区| 嫩草影院精品99| 熟妇人妻不卡中文字幕| 中文天堂在线官网| 91在线精品国自产拍蜜月| 老女人水多毛片| 日韩精品青青久久久久久| 免费观看性生交大片5| 成人二区视频| .国产精品久久| 国产黄a三级三级三级人| 国产精品国产三级国产专区5o| 国产成人freesex在线| 国产成年人精品一区二区| 少妇人妻精品综合一区二区| 国产精品99久久久久久久久| 亚洲乱码一区二区免费版| 综合色丁香网| 成人鲁丝片一二三区免费| 大香蕉久久网| 黑人高潮一二区| 少妇人妻精品综合一区二区| 黄色一级大片看看| 亚洲伊人久久精品综合| 国产淫片久久久久久久久| 日韩成人伦理影院| 哪个播放器可以免费观看大片| 一区二区三区乱码不卡18| 成人午夜精彩视频在线观看| 日韩av在线免费看完整版不卡| 乱人视频在线观看| 日日啪夜夜爽| 国产乱人视频| 女人十人毛片免费观看3o分钟| 成年av动漫网址| 欧美xxⅹ黑人| 男女国产视频网站| 亚洲不卡免费看| 国产有黄有色有爽视频| 亚洲av电影在线观看一区二区三区 | 亚洲综合色惰| 春色校园在线视频观看| 午夜精品一区二区三区免费看| 久久久午夜欧美精品| 一二三四中文在线观看免费高清| 免费黄频网站在线观看国产| 久久99热这里只频精品6学生| 免费看日本二区| 国产有黄有色有爽视频| 久久久久久国产a免费观看| 中文字幕亚洲精品专区| 国产熟女欧美一区二区| 国产精品日韩av在线免费观看| 麻豆乱淫一区二区| 波野结衣二区三区在线| 精品久久国产蜜桃| 婷婷色av中文字幕| 日本-黄色视频高清免费观看| 91久久精品国产一区二区三区| 精品99又大又爽又粗少妇毛片| 国产免费福利视频在线观看| 69人妻影院| 国产精品麻豆人妻色哟哟久久 | 男女视频在线观看网站免费| 国产一区有黄有色的免费视频 | 亚洲自偷自拍三级| 99热6这里只有精品| 亚洲在久久综合| 午夜激情久久久久久久| 99热这里只有精品一区| 国产永久视频网站| 在线观看免费高清a一片| 中文字幕久久专区| 黄片无遮挡物在线观看| 国产在视频线精品| 日韩欧美一区视频在线观看 | 亚洲乱码一区二区免费版| 久久99热这里只频精品6学生| 中文字幕久久专区| 精品国产三级普通话版| 国产美女午夜福利| 熟女电影av网| 国国产精品蜜臀av免费| 国产精品蜜桃在线观看| 国产免费福利视频在线观看| av在线播放精品| 赤兔流量卡办理| 亚洲精品日韩av片在线观看| 在线免费十八禁| .国产精品久久| 女人十人毛片免费观看3o分钟| 看黄色毛片网站| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 18禁在线播放成人免费| 亚洲av.av天堂| 日韩av在线大香蕉| 久久久欧美国产精品| 日韩大片免费观看网站| 亚洲精品一区蜜桃| 少妇丰满av| 蜜桃亚洲精品一区二区三区| 人人妻人人澡欧美一区二区| 成人性生交大片免费视频hd| 国产白丝娇喘喷水9色精品| 久久久久久久午夜电影| 午夜久久久久精精品| 黄片无遮挡物在线观看| 午夜精品在线福利| 少妇高潮的动态图| 国产精品日韩av在线免费观看| 亚洲图色成人| 亚洲国产精品sss在线观看| 99热6这里只有精品| 欧美日韩一区二区视频在线观看视频在线 | 午夜福利成人在线免费观看| 日韩人妻高清精品专区| 91久久精品国产一区二区成人| 中国美白少妇内射xxxbb| 亚洲av.av天堂| 欧美成人午夜免费资源| 黄片无遮挡物在线观看| 三级国产精品片| 国产亚洲5aaaaa淫片| 久久久久精品性色| 国产女主播在线喷水免费视频网站 | 久久精品夜色国产| 午夜福利高清视频| 午夜福利视频1000在线观看| 亚洲色图av天堂| 卡戴珊不雅视频在线播放| 最近中文字幕2019免费版| 亚洲婷婷狠狠爱综合网| 免费无遮挡裸体视频| 成人午夜高清在线视频| 日韩,欧美,国产一区二区三区| 在线观看美女被高潮喷水网站| 不卡视频在线观看欧美| 亚洲精品久久午夜乱码| 波多野结衣巨乳人妻| 亚洲av福利一区| 亚洲久久久久久中文字幕| 国产男女超爽视频在线观看| 一级毛片电影观看| 伊人久久国产一区二区| 亚洲第一区二区三区不卡| 国产黄片视频在线免费观看| 尾随美女入室| 久久久精品欧美日韩精品| 免费av观看视频| 嫩草影院精品99| 久久久欧美国产精品| 久久人人爽人人爽人人片va| 亚洲图色成人| 欧美日韩精品成人综合77777| 99热全是精品| 日韩人妻高清精品专区| 国产精品一区www在线观看| 精品国产一区二区三区久久久樱花 | 五月玫瑰六月丁香| 亚洲久久久久久中文字幕| 亚洲熟女精品中文字幕| 99热这里只有是精品在线观看| 在线观看av片永久免费下载| 大香蕉97超碰在线| 欧美区成人在线视频| 欧美成人a在线观看| 黄色配什么色好看| 搡老妇女老女人老熟妇| 97超碰精品成人国产| 国产成人福利小说| 精品一区二区三卡| 97超视频在线观看视频| 亚洲av在线观看美女高潮| 精华霜和精华液先用哪个| 99热这里只有是精品在线观看| 天天躁夜夜躁狠狠久久av| 一本一本综合久久| 插逼视频在线观看| 午夜福利在线观看吧| 网址你懂的国产日韩在线| 丝袜美腿在线中文| 亚洲经典国产精华液单| 全区人妻精品视频| 一级毛片久久久久久久久女| 小蜜桃在线观看免费完整版高清| 2018国产大陆天天弄谢| 国产日韩欧美在线精品| 亚洲一级一片aⅴ在线观看| 久热久热在线精品观看| 91aial.com中文字幕在线观看| 久久99热6这里只有精品| 神马国产精品三级电影在线观看| 99久久精品国产国产毛片| 在线天堂最新版资源| 日韩欧美精品免费久久| 免费观看av网站的网址| 一本久久精品| 插逼视频在线观看| 亚洲综合精品二区| 国产熟女欧美一区二区| 丰满少妇做爰视频| 高清午夜精品一区二区三区| 成人无遮挡网站| 最近中文字幕2019免费版| 亚洲国产最新在线播放| 国产精品一区二区在线观看99 | 亚洲精品第二区| 神马国产精品三级电影在线观看| 国产黄色视频一区二区在线观看| 久久久亚洲精品成人影院| 国产成人aa在线观看| 91午夜精品亚洲一区二区三区| 最近2019中文字幕mv第一页| 男女啪啪激烈高潮av片| 在线免费观看不下载黄p国产| 亚洲精品视频女| 久久久久性生活片| 久久99蜜桃精品久久| 亚洲精品中文字幕在线视频 | 亚洲国产日韩欧美精品在线观看| 日韩强制内射视频| 国产伦理片在线播放av一区| 综合色丁香网| 搞女人的毛片| 一二三四中文在线观看免费高清| 国产白丝娇喘喷水9色精品| 国产精品久久视频播放| 亚洲成人av在线免费| 春色校园在线视频观看| 麻豆久久精品国产亚洲av|