|
 |
 |
廣州市萬(wàn)通通風(fēng)設(shè)備有限公司
廣州市天河區(qū)長(zhǎng)興街岑村松崗小區(qū)五巷5號(hào)101鋪
電 話:13570323867 劉生
網(wǎng) 址:www.st0r.cn |
|
 |
 |
玻璃鋼軸流風(fēng)機(jī),
玻璃鋼離心風(fēng)機(jī),
塑料離心風(fēng)機(jī),
酸霧廢氣凈化處理塔,
玻璃鋼負(fù)壓風(fēng)機(jī),
玻璃鋼柜式風(fēng)機(jī),
玻璃鋼屋頂風(fēng)機(jī),
風(fēng)機(jī)配套產(chǎn)品,
煤氣加壓風(fēng)機(jī),
大型離心風(fēng)機(jī),
玻璃鋼斜流風(fēng)機(jī),
|
|
|
風(fēng)機(jī)常識(shí) | Common |
|
|
完全可逆軸流風(fēng)機(jī)的設(shè)計(jì)與優(yōu)化(離心風(fēng)機(jī),環(huán)保風(fēng)機(jī),空氣凈化器,酸霧塔,油水分離器) |
發(fā)布時(shí)間:2017/4/18 10:42:12
|
摘要:針對(duì)可逆軸流風(fēng)機(jī)S型葉片正反向運(yùn)行損失較大,效率難以大幅度提高的問(wèn)題,本文分析了組合葉柵相對(duì)于S型葉片的優(yōu)勢(shì),采用組合葉柵的構(gòu)想設(shè)計(jì)組合葉片可逆風(fēng)機(jī);采用可控渦扭向規(guī)律的設(shè)計(jì)方法設(shè)計(jì)出高效軸流風(fēng)機(jī)單轉(zhuǎn)子;為了探索可逆風(fēng)機(jī)組合葉片設(shè)計(jì)中前后排葉片的相對(duì)位置對(duì)風(fēng)機(jī)總性能的影響,對(duì)組合葉片設(shè)計(jì)進(jìn)行了基于不同軸向和周向位置的布局探索,共建立了16套不同的轉(zhuǎn)子模型,并進(jìn)行詳細(xì)的數(shù)值計(jì)算和分析對(duì)比,得出了可逆風(fēng)機(jī)組合葉片設(shè)計(jì)前后排葉片布局的最佳方案,最終得到的完全可逆風(fēng)機(jī)的性能滿足要求,正反風(fēng)效率均在85%以上。同時(shí),由于前后排葉片的完全對(duì)稱布置,保證了風(fēng)機(jī)正反風(fēng)性能的完全可逆性。
關(guān)鍵詞:軸流式通風(fēng)機(jī);可逆式;組合葉片;設(shè)計(jì)與優(yōu)化;數(shù)值模擬
中圖分類號(hào):TH432.1 文獻(xiàn)標(biāo)志碼:A
Design and Optimization of Fully Reversible Axial-flow Fan
Abstract: Aiming at the problem of greater losses in the forward and reverse operating of reversible axial-flow fan S-blade and problem of efficiency difficult to significantly improved, this paper analyzed the advantage of combined cascade relative to the S-shaped blade, so the concept of adopting combined cascade to design and combine reversible fan blades was pointed out, and the method of adopting controllable vortex torsional law was used to design high-efficient axial-flow fan with single rotor. Furthermore, in order to explore the influence of the relative positions of front and rear blades in the design of reversible fan combined blades on the overall performance of fan, the layout exploration was carried out on the combined blade’s design based on different axial and circumferential positions. Meanwhile, 16 sets of different rotor models were established, which were calculated numerically and analyzed in detail. So the optimal solution of designing the layout of front and rear blades for reversible fan was obtained. Finanlly, the obtained performance of fully reversible fan met the requirements, the efficiency in positive and negative position is above 85%. Meanwhile, the completely symmetrical layout of front and rear blades had ensured the fan performance was completely reversible.
Key words: axial-flow fan; reversible; combined blade; design and optimization; numerical simulation
引言
許多場(chǎng)合要求軸流風(fēng)機(jī)能完全反轉(zhuǎn)反風(fēng)[1-3] ,如地鐵、隧道、礦井和地下工程等?赡骘L(fēng)機(jī)承擔(dān)著地下空間與外部自然空間的氣體交換工作,特別在事故工況下,例如列車阻塞和突發(fā)火災(zāi)時(shí),需要風(fēng)機(jī)進(jìn)行應(yīng)急反向送風(fēng),而且要求風(fēng)量、風(fēng)壓與正向工作時(shí)相當(dāng),因此,地下工程通風(fēng)系統(tǒng)配置的基本都是可逆轉(zhuǎn)式軸流風(fēng)機(jī)。隨著各大城市地鐵建設(shè)的不斷發(fā)展,可逆風(fēng)機(jī)的需求越來(lái)越廣泛,研發(fā)一種正反向通風(fēng)效率都較高的可逆風(fēng)機(jī)成為很多學(xué)者致力實(shí)現(xiàn)的目標(biāo)。李俊超[4-7]等在國(guó)內(nèi)率先提出了一種采用普通翼型反向搭接構(gòu)成可逆風(fēng)機(jī)新翼型和一種S形反向?qū)ΨQ翼型,并采用優(yōu)化方法開發(fā)出了可逆風(fēng)機(jī)。隨著經(jīng)濟(jì)發(fā)展,我國(guó)的地鐵和公共隧道里程大大增加,對(duì)各種可逆軸流風(fēng)機(jī)的需求量也相應(yīng)增大。原有的可逆地鐵風(fēng)機(jī)設(shè)計(jì)技術(shù)不能完全滿足要求,相應(yīng)的改進(jìn)也不斷出現(xiàn)。
S翼型在可逆風(fēng)機(jī)的設(shè)計(jì)中由于葉型的后半段無(wú)論是正向還是反向流動(dòng)都不符合流動(dòng)規(guī)律,因此其效率難以得到大幅度的提高,而組合葉柵相對(duì)于S翼型在升力系數(shù)和失速攻角范圍都有了很大程度的改善。當(dāng)前后排葉片的軸向重合度為30%~50%時(shí),組合葉柵的氣動(dòng)性能明顯得到改善,這是因?yàn)榍傲腥~片背弧面的分離點(diǎn)明顯后置。由于前列葉柵的環(huán)量誘導(dǎo)作用,氣流產(chǎn)生轉(zhuǎn)折,減小了進(jìn)入后列葉片的攻角,使得后列葉片能夠與前列葉片一起,各自在最大升力系數(shù)附近的攻角下工作,因而,以這種方式將前后葉柵進(jìn)行組合,能夠使前后列葉片各自在最適宜的工作條件下工作,從而使得組合葉柵整體性能的優(yōu)勢(shì)得到充分體現(xiàn)。但是,目前組合葉柵依然存在葉型簡(jiǎn)單、負(fù)荷低、運(yùn)行效率低等問(wèn)題,然而毋庸置疑,在可逆風(fēng)機(jī)的葉型設(shè)計(jì)中它應(yīng)該是一種應(yīng)用潛力非常大的葉型。關(guān)于組合葉柵的優(yōu)化設(shè)計(jì)值得進(jìn)一步花大力氣進(jìn)行研究,本著這樣的目標(biāo),本文開展基于組合葉片的可逆軸流風(fēng)機(jī)的設(shè)計(jì)和優(yōu)化,以期設(shè)計(jì)出一種具有更好可逆性能的完全可逆風(fēng)機(jī)。
1 單轉(zhuǎn)子設(shè)計(jì)與性能模擬
1.1 單轉(zhuǎn)子氣動(dòng)設(shè)計(jì)方法介紹
1.1.1 設(shè)計(jì)方法
為了抑制葉根區(qū)強(qiáng)烈的旋渦流動(dòng)和減少分離損失,使風(fēng)機(jī)效率、風(fēng)量和風(fēng)壓達(dá)到設(shè)計(jì)要求,采用可控渦扭向規(guī)律的設(shè)計(jì)方法,把加功量集中在葉中區(qū)域,而葉中區(qū)域可以視為無(wú)粘區(qū),其二次流損失非常小。按照這樣的思路,可以根據(jù)以下步驟來(lái)設(shè)計(jì)扭向規(guī)律:
1) 首先根據(jù)設(shè)計(jì)流量值,由等環(huán)量規(guī)律計(jì)算出風(fēng)輪扭速Δcu沿葉高的分布。
根據(jù)中、尖部加功量加大,根部加功量小,等環(huán)量值進(jìn)行適當(dāng)?shù)匦拚T瓌t是,葉尖部分增加的少一點(diǎn),因?yàn)檫@里有較大的圓周速度,稍微增加一點(diǎn),產(chǎn)生的效果就相當(dāng)明顯;葉根部分可適當(dāng)多減少一些,因?yàn)檫@里的圓周速度較低;葉中部分增加到合適的程度。
2) 由Δcu(r)的確定方程、基元級(jí)壓升方程以及簡(jiǎn)化徑向平衡方程三方程聯(lián)立,求出相應(yīng)的軸向速度沿徑向的分布Ca(r),來(lái)保證流線的穩(wěn)定性。另外,用流量Q確定這種扭向規(guī)律下Ca(r)方程中的積分常數(shù)。最后用加功量和流量Q來(lái)檢驗(yàn)是否達(dá)到總的壓升和流量的設(shè)計(jì)指標(biāo)。
對(duì)于這些參數(shù)進(jìn)行無(wú)量綱化:
式中:腳標(biāo)t表示葉片尖部的參數(shù);u為圓周速度;p為壓力升。得無(wú)量綱方程組如下:
解上述方程組,其中積分方程的解中有常數(shù)項(xiàng),通過(guò)流量連續(xù)Q確定這種扭向規(guī)律下Ca(r)方程中的積分常數(shù)。
圖1 Δcu沿葉高的分布圖
1.1.2 計(jì)算機(jī)輔助設(shè)計(jì)系統(tǒng)
計(jì)算機(jī)輔助設(shè)計(jì)(CAD-Computer Aided Design)的應(yīng)用,可以大大提高產(chǎn)品的設(shè)計(jì)工作效率和質(zhì)量。因此,為了加快風(fēng)機(jī)的設(shè)計(jì)進(jìn)程和設(shè)計(jì)水平,本設(shè)計(jì)應(yīng)用了風(fēng)機(jī)CAD系統(tǒng),包括氣動(dòng)設(shè)計(jì)、幾何造型和CFD特性網(wǎng)圖計(jì)算。整套的CAD系統(tǒng)流程見圖2。
圖2 CAD設(shè)計(jì)流程圖
1.1.3 氣動(dòng)設(shè)計(jì)和幾何造型
在本氣動(dòng)設(shè)計(jì)的初始階段,采用簡(jiǎn)化徑向平衡方程基礎(chǔ)上的變環(huán)量設(shè)計(jì)方法,即在豐富的經(jīng)驗(yàn)數(shù)據(jù)基礎(chǔ)上,利用對(duì)沿葉高環(huán)量的控制,增加加功效益,提高葉片的氣動(dòng)效率。
1) 根據(jù)流量Q、設(shè)計(jì)全壓p及轉(zhuǎn)速n,由經(jīng)驗(yàn)公式可確定合理的外徑和輪轂比;
2) 在葉片氣動(dòng)設(shè)計(jì)上,考慮風(fēng)機(jī)內(nèi)部流場(chǎng)近似于圓柱面流動(dòng),因此采用簡(jiǎn)單徑向平衡方程設(shè)計(jì)方法,扭向規(guī)律按照可控渦設(shè)計(jì),控制環(huán)量沿葉高分布的初值,則按經(jīng)驗(yàn)?zāi)P徒o出;
3) 根據(jù)Δcu(r),Cu(r),Δca(r)的算出,可由速度三角形關(guān)系式計(jì)算出各截面的氣動(dòng)參數(shù)。由葉柵稠度的選定(查低速葉柵特性圖)以及經(jīng)驗(yàn)落后角公式,并經(jīng)過(guò)反復(fù)迭代計(jì)算,最終得到葉柵氣動(dòng)結(jié)構(gòu)參數(shù)和葉型幾何參數(shù)。
1.2 葉輪基本設(shè)計(jì)參數(shù)選取
以某一可逆風(fēng)機(jī)設(shè)計(jì)指標(biāo)為例,設(shè)計(jì)指標(biāo):全壓320~820Pa,流量6.8~17.5m3/h;直徑為0.9m、正反向的風(fēng)量全壓和功率基本一致,正反向效率大于76%;考慮到組合葉片無(wú)論正反向工作時(shí)后排葉片不可避免的損失,單轉(zhuǎn)子設(shè)計(jì)指標(biāo):壓升900Pa,流量10m3/h,效率90%;轉(zhuǎn)速1 450r/min,外徑900mm。
具體設(shè)計(jì)參數(shù)選擇如下。
1) 最優(yōu)輪轂比選取:d=0.65
2) Δcu的選取:
等環(huán)量設(shè)計(jì)——12.2(尖);14.78(中);18.76(根)
可控渦設(shè)計(jì)——12.2(尖);15(中);17.8(根)
3) 稠度: 0.463(尖); 0.663(中);1.02(根)
4) 積疊方式:按葉型前緣線進(jìn)行彎掠積疊
5) 葉片數(shù):11
三維造型見圖3。
圖3 單轉(zhuǎn)子葉片圖
1.3 單轉(zhuǎn)子氣動(dòng)校核
數(shù)值模擬程序采用商用軟件ANSYS/CFX求解三維流場(chǎng),求解器是CFX-Solver Manager,可以完成二維/三維的歐拉方程/N-S方程粘性求解;趯(duì)低速壓縮機(jī)內(nèi)部流動(dòng)的充分認(rèn)識(shí),數(shù)值模擬中采用了軸對(duì)稱假設(shè),對(duì)通風(fēng)機(jī)轉(zhuǎn)子進(jìn)行單通道三維粘性定常計(jì)算。求解方程為雷諾平均的三維粘性N-S方程,差分格式采用了二階精度的迎風(fēng)格式,湍流模型選用k-epsilon模型,前后級(jí)之間的連接面上采用周向平均的數(shù)據(jù)傳遞方式。
1) 計(jì)算網(wǎng)格與邊界條件
計(jì)算網(wǎng)格由NUMECA軟件中的AutoGrid5模塊生成,圖4和圖5為計(jì)算網(wǎng)格的子午面流道及葉輪轉(zhuǎn)子壁面網(wǎng)格,葉輪帶有1mm的葉尖間隙,網(wǎng)格數(shù)共約22萬(wàn),網(wǎng)格正交性>15°,延展比<3,長(zhǎng)寬比<3 000,均滿足計(jì)算要求。
邊界條件:進(jìn)口給定標(biāo)準(zhǔn)大氣壓;出口給定平均靜壓;計(jì)算域左右兩邊界面給定周期性邊界條件;葉輪機(jī)匣為絕對(duì)坐標(biāo)系下靜止的固體壁面;葉片和葉片輪轂為相對(duì)坐標(biāo)系下靜止的固體壁面。
圖4 子午面流道圖
圖5 單通道網(wǎng)格圖
2) 單轉(zhuǎn)子氣動(dòng)性能曲線見圖6~圖10。
圖6 單轉(zhuǎn)子全壓性能曲線圖
圖7 單轉(zhuǎn)子效率性能曲線圖
圖8 單轉(zhuǎn)子功率壓升性能曲線圖
圖9 轉(zhuǎn)子出口軸向速度分布圖
圖10 轉(zhuǎn)子出口總壓展向分布曲線圖
從性能曲線上我們可以看出,所設(shè)計(jì)的單轉(zhuǎn)子設(shè)計(jì)點(diǎn)流量10.16m3/s,效率90.62%,總壓升770.56Pa。
2 基于組合葉片設(shè)計(jì)的可逆風(fēng)機(jī)布局優(yōu)化
2.1 設(shè)計(jì)方案
為了描述前后排葉片的相對(duì)位置,定義周向弧度系數(shù)和軸向重合度如下:周向弧度系數(shù),即前列某一葉片與距離其最近的后列葉片的弧度夾角與同排兩相鄰葉片弧度夾角的比值;軸向相對(duì)位置,即用前后列葉片在70%葉高處沿弦長(zhǎng)方向重合距離與該處弦長(zhǎng)之比來(lái)描述。轉(zhuǎn)子模型的建立采用軸向相對(duì)位置為70%;葉高弦長(zhǎng)處重合度分別為10%、30%、50%、70%;周向相對(duì)位置采用周向弧度系數(shù)分別為25%、40%、50%、75%,建立了不同布局方案的轉(zhuǎn)子模型。上下游分別延長(zhǎng)了3個(gè)柵距,葉尖間隙1mm,對(duì)不同重合度和相對(duì)柵距系數(shù)下布局方案進(jìn)行模擬計(jì)算。采用UG軟件建立轉(zhuǎn)子模型(見圖11)。
圖11 轉(zhuǎn)子模型示意圖
2.2 不同布局方案的數(shù)值模擬
在本方案數(shù)值計(jì)算時(shí),選用商用CFD軟件NUMECA中的AutoGrid5進(jìn)行網(wǎng)格生成,采用ANSYS旗下的CFX-Solver軟件進(jìn)行計(jì)算,后處理采用的是CFX-Post軟件。網(wǎng)格數(shù)45萬(wàn),上下游分別延長(zhǎng)了三個(gè)柵距,組合葉片的子午道流面圖如圖12所示。
圖12 子午道流面圖
湍流模型采用k-elpson,邊界條件:進(jìn)口給定標(biāo)準(zhǔn)大氣壓;出口給定平均靜壓;計(jì)算域左右兩邊界面給定周期性邊界條件;葉輪機(jī)匣為絕對(duì)坐標(biāo)系下靜止的固體壁面;葉片和葉片輪轂為相對(duì)坐標(biāo)系下靜止的固體壁面。
對(duì)不同的軸向和周向布局方案數(shù)值模擬所得轉(zhuǎn)子效率、功率、流量特性變化見圖13~圖15。圖13為不同布局組合葉片效率變化曲線,此處定義不同布局方案下組合葉片轉(zhuǎn)子效率與設(shè)計(jì)點(diǎn)效率之比;圖14為不同布局組合葉片壓升變化曲線,此處定義不同布局方案下組合葉片轉(zhuǎn)子壓升與設(shè)計(jì)點(diǎn)壓升之比;圖15為不同布局組合葉片流量變化曲線,此處定義不同布局方案下組合葉片轉(zhuǎn)子流量與設(shè)計(jì)點(diǎn)流量之比。
圖13 不同布局組合葉片效率性能曲線圖
圖14 不同布局組合的葉片壓升性能曲線圖
圖15 不同布局組合葉片流量曲線圖
從不同布局形式的效率曲線來(lái)看,效率最高點(diǎn)為軸向重合度10%,周向弧度系數(shù)為75%布局形式,最高點(diǎn)效率系數(shù)為1.008,效率為85.56%,正反向效率完全相同。通過(guò)分析可以看出,周向弧度系數(shù)的改變對(duì)于固定周向重合度的轉(zhuǎn)子葉片性能影響不大,與壓氣機(jī)和風(fēng)扇設(shè)計(jì)的串列葉片不同,串列葉片的優(yōu)勢(shì)體現(xiàn)在后排葉片能減少前排葉片吸力面的分離,從而改善整體性能。而分離發(fā)生嚴(yán)重的工況大多是非設(shè)計(jì)工況,對(duì)于風(fēng)機(jī)的設(shè)計(jì)情況則不同,風(fēng)機(jī)幾乎所有的工作時(shí)間都是設(shè)計(jì)工況,前排葉片的吸力面并沒有發(fā)生明顯的分離現(xiàn)象。所以對(duì)于組合葉片來(lái)說(shuō),后排葉片的插入并沒有明顯的改善轉(zhuǎn)子的整體性能,隨著軸向重合度的不斷增加,組合葉輪轉(zhuǎn)子葉根前后排葉片有大范圍的重疊,造成流動(dòng)堵塞,使得正反向轉(zhuǎn)子性能同時(shí)下降,周向重合度對(duì)于組合葉片的正反向性能具有較大影響。
2.3 最佳布局方案轉(zhuǎn)子性能
圖16~圖18分別給出了最高效率點(diǎn)組合葉輪轉(zhuǎn)子的效率、功率和壓升特性曲線,組合葉輪轉(zhuǎn)子正反向工作效率均達(dá)85%以上,全壓升800Pa,流量9m3/h,均在可逆風(fēng)機(jī)設(shè)計(jì)指標(biāo)范圍內(nèi),且具有較高的效率。這對(duì)于常年日工作時(shí)間在12小時(shí)以上的可逆風(fēng)機(jī)來(lái)說(shuō),節(jié)能效益非常明顯。
圖16 最佳布局組合葉片的壓升特性圖
圖17 最佳布局組合葉片的效率特性圖
圖18 最佳布局組合葉片的功率特性圖
3 結(jié)論
1) 基于組合葉片設(shè)計(jì)的可逆風(fēng)機(jī)性能優(yōu)劣主要取決于構(gòu)成組合葉片的單轉(zhuǎn)子性能,布局優(yōu)化只能在一定范圍內(nèi)減少后列葉片。由于反向布置而產(chǎn)生的損失,采用可控渦設(shè)計(jì)方法對(duì)于軸流風(fēng)機(jī)單轉(zhuǎn)子的設(shè)計(jì)是行之有效的一種設(shè)計(jì)方法;
2) 組合葉片的前后排葉片相對(duì)位置直接影響著組合葉片的整體性能,本文通過(guò)對(duì)不同布局形式的組合葉片數(shù)值模擬,分析對(duì)比可知,前后排葉片的周向相對(duì)位置對(duì)風(fēng)機(jī)總性能的影響不大,最佳周向位置為周向弧度系數(shù)為10%。由于風(fēng)機(jī)主要在設(shè)計(jì)工況工作,前列葉片在設(shè)計(jì)工況下并沒有明顯的流動(dòng)分離,所以后排葉片的插入未能明顯改善前排葉片的流動(dòng),而軸向相對(duì)位置對(duì)組合葉片性能影響較為明顯。結(jié)果表明,軸向重合度在30%位置處,組合葉輪具有較好的流動(dòng)性能;
3) 該設(shè)計(jì)方法所得可逆風(fēng)機(jī)轉(zhuǎn)子的壓升和流量均達(dá)到性能要求,且具有較高的效率,同時(shí)由于組合葉片的前后排對(duì)稱布局設(shè)計(jì),結(jié)構(gòu)正反向?qū)ΨQ,能夠保證風(fēng)機(jī)性能的完全可逆而且都在較高的效率點(diǎn)工作;
4) 由于時(shí)間關(guān)系,本文并未對(duì)后列葉片的逆向損失機(jī)理進(jìn)行研究,下一步工作將基于逆流損失機(jī)理進(jìn)行探索,如若能更好的把握逆流損失并進(jìn)行有效控制,該設(shè)計(jì)方案將得到進(jìn)一步的優(yōu)化。
參 考 文 獻(xiàn)
[1] 李國(guó)慶. 城市軌道交通通風(fēng)空調(diào)系統(tǒng)技術(shù)發(fā)展新趨勢(shì)[J].都市快軌交通,2004,17(6):5- 7.
[2] Gabay D. Fire safety: a short history in the Paris subway. 2002ITA Open Session: Fire and Life Safety [J]. Tunnelling andUnderground space Technology 2002, 17:139- 143.
[3] Carvel R O, Beard A N and Jowitt P W. The influence of longitudinal ventilation systems on fires in tunnels [J]. Tunnelling and underground space technology 2001, 16:3- 21.
[4] 李超俊,趙德文,魏百鎖. 雙向軸流風(fēng)機(jī)的優(yōu)化設(shè)計(jì)與實(shí)驗(yàn)研究[J]. 機(jī)械工程學(xué)報(bào),1991,27(5):54- 60.
[5] 李超俊,趙德文,魏百鎖.直接反轉(zhuǎn)反風(fēng)軸流風(fēng)機(jī)的理論與實(shí)驗(yàn)研究[J].西安交通大學(xué)學(xué)報(bào),1991,25(3):9- 13.
[6] 李超俊,魏百鎖.一種用于可逆軸流風(fēng)機(jī)的新翼型:雙圓弧S形翼型的研究[J].機(jī)械工程報(bào),1992,28(4):37- 41.
[7] 田琳,蘇莫明.雙向軸流通風(fēng)機(jī)設(shè)計(jì)方法的研究[J].風(fēng)機(jī)技術(shù),2010(5):32-34,39.
更多環(huán)保風(fēng)機(jī)可用我司產(chǎn)品類型:玻璃鋼風(fēng)機(jī),玻璃鋼軸流風(fēng)機(jī),玻璃鋼離心風(fēng)機(jī),負(fù)壓風(fēng)機(jī),玻璃鋼斜流風(fēng)機(jī),玻璃鋼柜式風(fēng)機(jī)。
|
[ 關(guān)閉本頁(yè) ] |
|
 |
|
|