彭 思, 王漢封,2, 李石清

(1.中南大學(xué) 土木工程學(xué)院,長(zhǎng)沙 410075; 2.高速鐵路建造技術(shù)國(guó)家工程實(shí)驗(yàn)室(中南大學(xué)), 長(zhǎng)沙410075)

頂部定常吸氣對(duì)高層建筑模型氣動(dòng)力的影響

彭 思1, 王漢封1,2, 李石清1

(1.中南大學(xué) 土木工程學(xué)院,長(zhǎng)沙 410075; 2.高速鐵路建造技術(shù)國(guó)家工程實(shí)驗(yàn)室(中南大學(xué)), 長(zhǎng)沙410075)

為減少高層建筑受到的風(fēng)荷載，提高高層建筑抗風(fēng)性能，提出一種頂部吸氣的主動(dòng)控制方法。在一高寬比為H/d=5的正方形截面柱體的高層建筑模型上，通過(guò)在其頂部前邊緣開(kāi)設(shè)狹縫進(jìn)行吸氣來(lái)實(shí)現(xiàn)流動(dòng)控制。通過(guò)風(fēng)洞試驗(yàn)研究了頂部吸氣對(duì)氣動(dòng)力與頂部分離流的影響. 利用流動(dòng)可視化與流場(chǎng)測(cè)試結(jié)果揭示了頂部狹縫吸氣的控制機(jī)理. 實(shí)驗(yàn)對(duì)比了不同吸氣系數(shù)Q(=U/U∞,U為狹縫吸氣速度) 對(duì)氣動(dòng)力控制效果的影響, 并對(duì)Q=0, 1和3的三種工況下風(fēng)壓分布與頂部剪切流進(jìn)行了詳細(xì)對(duì)比. 研究發(fā)現(xiàn)狹縫吸氣改變了頂部分離流特性, 并對(duì)模型所有高度上的氣動(dòng)力均有顯著影響.Q=1時(shí)控制效果最佳, 脈動(dòng)阻力與脈動(dòng)升力分別減小17.8%和45.5%. 此時(shí)頂面分離流被削弱并伴隨有再附現(xiàn)象且湍流強(qiáng)度較大, 最有利于頂部剪切流與尾流間的動(dòng)量交換, 從而削弱柱體展向渦脫落與脈動(dòng)氣動(dòng)力.

空氣動(dòng)力學(xué)；有限長(zhǎng)方柱體；定常吸氣；流動(dòng)分離；流動(dòng)控制

針對(duì)超高層建筑的氣動(dòng)優(yōu)化措施大體可以分為被動(dòng)控制和主動(dòng)控制[1]. 被動(dòng)控制最常見(jiàn)的是結(jié)構(gòu)物的外形優(yōu)化, 如：采用倒角、圓角、凹角或角部擾流片, 柱體貫穿開(kāi)洞, 沿高度改變建筑橫截面面積, 增加棱柱體頂面棱邊和側(cè)面螺旋線的數(shù)量[2-7]. 相比于被動(dòng)控制, 主動(dòng)控制方法以其可控性更能適用于復(fù)雜的風(fēng)環(huán)境. 其中吹、吸氣方法是常見(jiàn)的主動(dòng)控制措施之一. 文獻(xiàn)[8]在高寬比H/d為6的高層建筑模型不同高度上開(kāi)設(shè)吸氣孔, 研究了不同風(fēng)向角下模型側(cè)面和背面定常吸氣的控制效果. 實(shí)驗(yàn)發(fā)現(xiàn)側(cè)面吸氣比背面吸氣效果更好, 且開(kāi)孔距流動(dòng)分離點(diǎn)越近, 控制效果越好. 文獻(xiàn)[9-11]通過(guò)數(shù)值模擬研究了高層建筑側(cè)面全高和分段的吸、吹氣控制, 發(fā)現(xiàn)氣流角度、開(kāi)孔位置和開(kāi)孔寬度對(duì)控制效果都有影響, 但主要受無(wú)量綱的吸氣量系數(shù)控制. 從能耗和折減效率的角度分析, 發(fā)現(xiàn)吸、吹氣有不同的適用范圍, 且分段控制效果不如全高控制. 文獻(xiàn)[12-13]進(jìn)一步探討了被動(dòng)控制和主動(dòng)控制相結(jié)合的方法, 研究了方形、Y形截面形式、錐臺(tái)形式和帶切角或凹角形式的高層建筑模型在0°至45°風(fēng)向角下通過(guò)兩側(cè)吸氣的控制效果. 實(shí)驗(yàn)發(fā)現(xiàn), 除錐臺(tái)模型外, 其他各吸氣模型對(duì)橫風(fēng)向脈動(dòng)力都有很好的抑制效果.

有研究發(fā)現(xiàn), 有限長(zhǎng)柱體頂部形成下掃流對(duì)其氣動(dòng)力有很大影響[14-16], 但是能否通過(guò)控制下掃流實(shí)現(xiàn)對(duì)高層建筑氣動(dòng)力的控制仍有待進(jìn)一步驗(yàn)證. 本文通過(guò)風(fēng)洞試驗(yàn), 在正方形截面的高層建筑簡(jiǎn)化模型頂部施加定常吸氣, 系統(tǒng)地研究了吸氣對(duì)氣動(dòng)力的影響規(guī)律, 并通過(guò)流動(dòng)可視化、流場(chǎng)測(cè)量等方法對(duì)頂部吸氣控制的機(jī)理進(jìn)行了探討.

1 實(shí)驗(yàn)介紹

所有試驗(yàn)在一直流式風(fēng)洞內(nèi)進(jìn)行, 其實(shí)驗(yàn)段截面尺寸為450 mm×450 mm, 風(fēng)速范圍0～40 m/s, 湍流度小于0.5%. 被測(cè)高層建筑簡(jiǎn)化為一正方形截面棱柱, 寬d=40 mm, 高寬比H/d=5. 模型由韌性類(lèi)ABS樹(shù)脂經(jīng)3D打印制成, 剛度及表面光滑度都能滿(mǎn)足實(shí)驗(yàn)要求. 安裝模型前, 對(duì)風(fēng)洞壁面邊界層進(jìn)行標(biāo)定, 其厚度約為0.4d, 即模型0.4d以上高度都處于均勻來(lái)流中. 圖1(a)給出了模型示意圖與坐標(biāo)系的定義. 坐標(biāo)系原點(diǎn)為方柱底面中心. 模型頂面迎風(fēng)側(cè)距邊緣1 mm處開(kāi)設(shè)有一寬1 mm長(zhǎng)36 mm的吸氣細(xì)縫. 模型底部的出氣管與放置在風(fēng)洞外的抽氣泵相連, 并可通過(guò)流量計(jì)監(jiān)測(cè)吸氣流量, 見(jiàn)圖1(b).在z*= 1, 2, 3, 4和4.5等(z*=h/d,h為測(cè)點(diǎn)高度,d為方柱寬度) 不同高度上, 布置五層測(cè)壓孔. 每層各個(gè)面的測(cè)壓點(diǎn)為5個(gè), 沿柱體表面均勻布置.測(cè)壓試驗(yàn)使用美國(guó)Scannivalve公司ZOC33/64PxX2高頻壓力掃描閥進(jìn)行測(cè)量, 測(cè)壓試驗(yàn)采樣頻率為625 Hz, 采樣時(shí)間32 s, 測(cè)壓系統(tǒng)見(jiàn)圖1(b). 實(shí)驗(yàn)中自由來(lái)流風(fēng)速U∞= 9.8 m/s, 基于d與U∞的Re=2.74×104. 通過(guò)流量計(jì)監(jiān)測(cè)吸氣的體積流量, 并依據(jù)吸氣狹縫面積, 估算狹縫吸氣速度U, 定義吸氣系數(shù)Q=U/U∞.試驗(yàn)研究了Q=0～4范圍內(nèi)頂部狹縫吸氣對(duì)模型氣動(dòng)力與頂部剪切流的影響規(guī)律.

(1)

(2)

式中:pi表示各測(cè)點(diǎn)風(fēng)壓,p為風(fēng)洞靜壓,ρ為空氣密度,為測(cè)壓點(diǎn)脈動(dòng)風(fēng)壓的均方根值(RMS值). 根據(jù)測(cè)點(diǎn)代表的面積, 可計(jì)算得到模型總體所受阻力系數(shù)Cd,whole, 脈動(dòng)阻力系數(shù)和脈動(dòng)升力系數(shù).

圖1 實(shí)驗(yàn)裝置

2 結(jié)果分析與討論

2.1 氣動(dòng)力

圖4給出了模型各高度上局部的升力系數(shù)自功率譜. 為方便對(duì)比, 所有工況下各層的功率譜均以無(wú)控制時(shí)z*=1的譜峰值強(qiáng)度作為參考值,按式(3)無(wú)量綱化為dB. 頻率則以來(lái)流速度U∞和方柱寬度d按式(4)無(wú)量綱.

(3)

(4)

式中:pQ=i,z*=j為Q=i,z*=j時(shí)功率譜峰值強(qiáng)度,pQ=0,z*=1為Q=0,z*=1時(shí)功率譜峰值強(qiáng)度,f為頻率.

圖3 Q=0、1和3時(shí)方柱氣動(dòng)力沿展向的分布

Fig.3 Distribution of the aerodynamic forces at different spanwise positions withQ=0, 1 and 3

圖4 方柱在Q=0、1和3時(shí)的各層升力系數(shù)自功率譜

Fig.4 Power spectral density of the lift at different spanwise positions withQ=0,1 and 3

無(wú)論吸氣與否, 模型所有高度上的無(wú)量綱渦脫頻率都為0.116, 這表明頂部吸氣雖然對(duì)脈動(dòng)升力等有明顯削弱作用, 但其不會(huì)改變模型渦脫落頻率.Q=0時(shí), 模型所有高度上的升力系數(shù)譜都具有明顯峰值, 說(shuō)明在不吸氣時(shí), 所有高度上展向渦都

具有強(qiáng)烈的周期性. 此外, 由于模型頂部所形成的下掃流具有削弱展向渦的特性, 因此, 越接近模型頂部功率譜峰值越弱. 當(dāng)Q=1時(shí), 所有高度上功率譜峰值相對(duì)于Q=0時(shí)都明顯減弱了, 且越接近頂部削弱得越顯著. 而當(dāng)Q=3時(shí), 功率譜峰值又有所增強(qiáng).

為進(jìn)一步揭示頂部吸氣對(duì)模型氣動(dòng)力的影響機(jī)理, 圖5給出了模型總體Cd,whole和Cl,whole的時(shí)程曲線以及Cl,whole對(duì)應(yīng)的時(shí)頻分析結(jié)果. 由圖 5可得到：首先,Q=0時(shí),Cl,whole大幅波動(dòng)的情形在時(shí)間上具有明顯的間歇性, 時(shí)頻譜中峰值是不連續(xù)的. 這與文獻(xiàn)[18-20]所觀察到的現(xiàn)象一致, 即有限長(zhǎng)方柱尾流中會(huì)隨機(jī)出現(xiàn)兩種典型狀態(tài): 一種是方柱兩側(cè)展向渦有規(guī)律地交替脫落, 此時(shí)Cl,whole大幅振動(dòng)；另一種是展向渦準(zhǔn)對(duì)稱(chēng)地脫落, 對(duì)應(yīng)Cl,whole振幅很小且無(wú)明顯周期性. 顯然, 前者所對(duì)應(yīng)的升力脈動(dòng)值要顯著大于后者. 其次,Q=1時(shí),Cd,whole和Cl,whole時(shí)程曲線表現(xiàn)得較為平穩(wěn), 這說(shuō)明類(lèi)似卡門(mén)渦街的交替渦脫模態(tài)在很大程度上被抑制了. 而當(dāng)Q=3時(shí),Cl,whole的周期性波動(dòng)再次變強(qiáng).

圖5 模型總體氣動(dòng)力特性

2.2 頂部剪切流

為進(jìn)一步探討頂部剪切流的特性, 進(jìn)行了流動(dòng)顯示實(shí)驗(yàn). 煙線布置見(jiàn)圖1(a), 發(fā)煙材料為石蠟. 實(shí)驗(yàn)風(fēng)速為1.12 m/s, 并以此調(diào)節(jié)吸氣速度, 研究Q=0、1和3工況下方柱頂部流動(dòng)特性. 由于方柱沒(méi)有雷諾數(shù)效應(yīng), 故該風(fēng)速下的實(shí)驗(yàn)結(jié)果仍能為測(cè)壓實(shí)驗(yàn)現(xiàn)象的解釋提供參考. 圖6給出了Q=0、1和3時(shí), 模型頂部剪切流的可視化結(jié)果. 由圖 6可以看出,Q=0時(shí), 剪切流在模型頂部邊緣發(fā)生分離, 分離流直接掠過(guò)模型頂部, 在其下游形成下掃流, 并與模型兩側(cè)形成的展向渦發(fā)生作用. 這一下掃流具有將展向渦向兩側(cè)分開(kāi)的趨勢(shì)并能削弱展向渦的強(qiáng)度, 降低其脫落頻率(相對(duì)于二維柱體)[21-22]. 當(dāng)頂部前邊緣吸氣系數(shù)Q=1時(shí), 頂部分離流被明顯抑制并向模型頂部靠近. 在狹縫吸氣的作用下, 分離流在頂部的下游部分發(fā)生了再附著. 隨著向下游的發(fā)展, 流動(dòng)再次在模型頂部后邊緣發(fā)生分離, 并伴隨有較為顯著的不穩(wěn)定性. 從流動(dòng)可視化圖中可以看出, 該剪切流與模型展向渦在模型后部發(fā)生了較為強(qiáng)烈的相互作用, 這可能是此時(shí)脈動(dòng)升力達(dá)到最小的原因. 當(dāng)Q增大到3時(shí), 模型前邊緣剪切流出現(xiàn)了另一種完全不同的形態(tài). 此時(shí), 頂部前邊緣的流動(dòng)分離被完全抑制, 流動(dòng)緊貼頂部表面向下游運(yùn)動(dòng), 在頂面尾部發(fā)生分離, 且剪切流沒(méi)出現(xiàn)如Q=1時(shí)的不穩(wěn)定擺動(dòng).

圖6 Q=0、1和3時(shí)方柱頂部流動(dòng)顯示

圖7 方柱頂部剪切流測(cè)量結(jié)果

圖 7(b)給出了上述測(cè)點(diǎn)位置上的流向脈動(dòng)速度Iuu.Q=0時(shí), 頂部邊緣形成的分離流的剪切層內(nèi)Iuu最強(qiáng), 但其距離頂面較遠(yuǎn). 前緣吸氣在Q=1時(shí)不僅顯著縮小了分離區(qū), 而且仍保持了頂面附近較強(qiáng)的Iuu. 而當(dāng)Q增大到3時(shí), 分離流被完全抑制后, 來(lái)自柱體前部的自由流緊貼著頂面流過(guò), 此時(shí)頂面附近Iuu與自由流相同, 其強(qiáng)度幾乎為零. 由頂面剪切流的流動(dòng)可視化與測(cè)量結(jié)果可知,Q=0時(shí), 頂面前緣分離流直接掠過(guò)模型頂部, 并會(huì)在頂面上形成分離泡. 而當(dāng)Q=1時(shí), 頂面上的分離流會(huì)發(fā)生再附著并明顯抑制分離泡, 且此時(shí)頂面流動(dòng)具有較高的湍流度, 這有利于模型尾流與外部高速流體進(jìn)行動(dòng)量交換, 并削弱展向漩渦脫落. 而當(dāng)Q增大到3時(shí), 盡管頂部分離流被完全抑制, 但頂面流動(dòng)湍流度幾乎為零, 這可能反而削弱了尾流與外部高速流體的動(dòng)量交換, 所以控制效果反而較Q=1時(shí)有所減弱.

3 結(jié) 論

采用風(fēng)洞試驗(yàn), 研究了H/d=5的正方形截面高層建筑簡(jiǎn)化模型頂部狹縫吸氣對(duì)其氣動(dòng)力與頂部剪切流的控制效果, 得到如下結(jié)論:

1)頂部狹縫吸氣能夠略微減小模型時(shí)均阻力, 并顯著抑制模型脈動(dòng)阻力與脈動(dòng)升力. 就本文吸氣設(shè)置情況而言, 當(dāng)吸氣系數(shù)Q=1時(shí), 氣動(dòng)力控制效果最為顯著, 相對(duì)未控制工況, 時(shí)均阻力減小3.60%, 脈動(dòng)阻力與脈動(dòng)升力分別減小17.85%和45.57%.

2)頂部狹縫吸氣不僅對(duì)頂部附近氣動(dòng)力有影響, 而且對(duì)模型所有高度上的氣動(dòng)力均有顯著控制效果. 頂部吸氣可削弱模型展向旋渦脫落, 但并未改變其渦脫落頻率.

3)狹縫吸氣對(duì)模型頂部分離流有顯著的加速作用, 并顯著改變了分離流特性.Q=1時(shí), 模型頂部流動(dòng)分離被明顯削弱, 并在頂面后部發(fā)生再附著；而當(dāng)Q=3時(shí), 流動(dòng)分離被完全抑制, 且頂面上的邊界層也基本消失.

[1] TAMURA Y, KAREEM A. Advancedstructural wind engineering[M].[S.l.]: Springer Japan, 2013: 377-383. DOI:10.1007/978-4-431-54337-4.

[2] KIM Y C, BANDI E K, YOSHIDA A, et al. Response characteristics of super-tall buildings-Effects of number of sides and helical angle[J]. Journal of Wind Engineering and Industrial Aerodynamics， 2015, 145: 252-262. DOI:10.1016/j.jweia.2015.07.001.

[3] TANAKA H, TAMURA Y, OHTAKE K, et al. Experimental investigation of aerodynamic forces and wind pressures acting on tall buildings with various unconventional configurations[J]. Journal of Wind Engineering and Industrial Aerodynamics， 2012, 107/108: 179-191. DOI:10.1016/j.jweia.2012.04.014.

[4] KAREEM A, KIJEWSKI T, TAMURA Y. Mitigation of motions of tall buildings with specific examples of recent applications[J]. Wind & Structures And International Journal， 1999, 2(3): 201-251.

[5] KWOK K C S, BAILEY P A. Aerodynamicdevices for tall buildings and structures[J]. Journal of Engineering Mechanics， 1987, 113(3): 349-365.

[6] DUTTON R, ISYUMOV N. Reduction of tall building motion by aerodynamic treatments[J]. Journal of Wind Engineering and Industrial Aerodynamics， 1990, 36(1-3): 739-744.

[7] GU Ming, QUAN Yong. Across-wind loads of typical tall buildings[J]. Journal of Wind Engineering and Industrial Aerodynamics， 2004, 92(13)：1147-1165. DOI:10.1016/j.jweia.2004.06.004.

[8] ZHANG Hongfu, XIN Dabo, OU Jinping. Steady suction for controlling across-wind loading of high-rise buildings[J]. The Structural Design of Tall and Special Buildings， 2016, 25(15): 785-800. DOI:10.1002/tal.1283.

[9] ZHENG Chaorong, ZHANG Yaochun. Numerical investigation on the drag reduction properties of a suction controlled high-rise building[J]. Journal of Zhejiang University-Science A， 2010, 11(7): 477-487. DOI:10.1631/jzus.A0900593.

[10]鄭朝榮. 高層建筑風(fēng)荷載吸/吹氣控制的數(shù)值模擬研究[D]. 哈爾濱:哈爾濱工業(yè)大學(xué), 2010.

ZHENG Chaorong. Numerical investigation of wind loads on high-rise buildings controlled by suction/blowing[D]. Harbin: Harbin Institute of Technology, 2010.

[11]鄭朝榮，張耀春. 分段吸氣高層建筑減阻性能的數(shù)值研究[J]. 力學(xué)學(xué)報(bào),2011, 43(2): 372-380.

ZHENG Chaorong, ZHANG Yaohchun. Numerical investigation of drag-reducing property for a high-rise building under subsection suction[J]. Chinese Journal of Theoretical and Applied Mechanics, 2011, 43(2)：372-380.

[12]張繼同. 主動(dòng)吸氣與被動(dòng)氣動(dòng)控制下超高層建筑風(fēng)荷載特性研究[D]. 哈爾濱:哈爾濱工業(yè)大學(xué), 2015.

ZHANG Jitong. Characteristics of wind loads on super high-rise buildings under active suction control and passive aerodynamic control[D].Harbin: Harbin Institute of Technology, 2015.

[13]鄭朝榮，張繼同，張智棟. 凹角與吸氣控制下高層建筑平均風(fēng)荷載特性試驗(yàn)研究[J]. 建筑結(jié)構(gòu)學(xué)報(bào)，2016(10): 125-131.

ZHENG Chaorong, ZHANG Jitong, ZHANG Zhidong. Experimental investigation on characteristics of mean wind loads of high-rise buildings controlled by corner recession and air-suction[J]. Journal of Building Structures, 2016(10): 125-131.

[14]WANG Hanfeng, ZHOU Yu. The finite-length square cylinder near wake[J]. Journal of Fluid Mechanics， 2009, 638: 453-490. DOI:10.1017/S0022112009990693.

[15]WANG Hanfeng, ZHOU Yu, CHAN C K, et al. Effect of initial conditions on interaction between a boundary layer and a wall-mounted finite-length-cylinder wake[J]. Physics of Fluids, 2006, 18(6): 065106. DOI:10.1063/1.2212329.

[16]UFFINGER T, ALI I, BECKER S. Experimental and numerical investigations of the flow around three different wall-mounted cylinder geometries of finite length[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2013, 119: 13-27. DOI: 10.1016/j.jweia.2013.05.006.

[17]CHEN Wenli, LI Hui, HU Hui. An experimental study on a suction flow control method to reduce the unsteadiness of the wind loads acting on a circular cylinder[J]. Experiments in Fluids, 2014, 55：1707. DOI:10.1007/s00348-014-1707-7.

[18]SATTARI P, BOURGEOIS J A, MARTINUZZI R J. On the vortex dynamics in the wake of a finite surface-mounted square cylinder[J]. Experiments in Fluids, 2012, 52(5): 1149-1167. DOI:10.1007/s00348-011-1244-6.

[19]王漢封，徐萌霞. 有限長(zhǎng)正方形棱柱繞流的雙穩(wěn)態(tài)現(xiàn)象[J]. 力學(xué)與實(shí)踐, 2013, 35(3): 46-50.

WANG Hanfeng, XU Mengxia. Bistable phenomenon of the flow around a finite-length square prism[J]. Mechanics in Engineering, 2013, 35(3): 46-50.

[20]王漢封，趙東偉，鄒超，等. 高寬比為5的懸臂正方形棱柱氣動(dòng)力特性[J]. 實(shí)驗(yàn)流體力學(xué),2015, 29(6): 8-15.

WANG Hanfeng, ZHAO Dongwei,ZOU Chao, et al. Aerodynamic forces of a cantilevered square prism with aspect ratio of 5[J]. Journal of Experiments in Fluid Mechanics, 2015, 29(6): 8-15.

[21]SAKAMOTO H, ARIE M. Vortex shedding from a rectangular prism and a circular cylinder placed vertically in a turbulent boundary layer[J]. Journal of Fluid Mechanics, 1983, 126(126): 147-165.

[22]PARK C, LEE S. Free end effects on the near wake flow structure behind a finite circular cylinder[J].Journal of Wind Engineering and Industrial Aerodynamics, 2000, 88(2): 231-246. DOI:10.1016/S0167-6105(00)00051-9.

(編輯趙麗瑩)

Thefree-endsteadysuctioneffectsontheaerodynamicforcesofhigh-risebuildingmodel

PENG Si1, WANG Hanfeng1,2, LI Shiqing1

(1.School of Civil Engineering, Central South University, Changsha 410075, China; 2.National Engineering Laboratory for High-speed Railway Construction (Central South University), Changsha 410075, China)

To reduce wind loads on high-rise buildings, a new active aerodynamic control method named steady suction on the top of the building was introduced. The effects of free-end leading edge steady suction on the aerodynamic forces and free-end separation of a finite-length square cylinder were experimentally investigated in a wind tunnel. Flow visualization and flow field measurement were utilized to reveal the control mechanism. The dependence of aerodynamic forces on the suction ratioQ(=U/U∞) was studied. The pressure distribution and the shear flow above the free end were compared in detail forQ=0, 1 and 3. The results show that the slot suction changes the flow separation on the free end significantly and also the aerodynamic forces on the whole cylinder span. The best control result appears atQ=1, with the fluctuation drag and lift reduced by 17.8% and 45.5%, respectively. At this suction ratio, the flow separation at the leading edge is weakened with reattachment and strong turbulence fluctuation presence on the free end, which results in strong momentum transport between the free-end shear flow and the wake and suppresses the vortex shedding and aerodynamic forces efficiently.

aerodynamic forces; finite-length square cylinder; steady suction; flow separation; flow control

10.11918/j.issn.0367-6234.201703123

O355

A

0367-6234(2017)12-0159-06

2017-03-24

國(guó)家自然科學(xué)基金(11472312)

彭 思(1990—)，男，碩士研究生；

王漢封(1976—)，男，教授,博士生導(dǎo)師

王漢封，wanghf@csu.edu.cn