Please use this identifier to cite or link to this item: http://hdl.handle.net/11455/2536
DC FieldValueLanguage
dc.contributor游明輝zh_TW
dc.contributor蕭飛賓zh_TW
dc.contributor.advisor郭正雄zh_TW
dc.contributor.advisorCheng-Hsiung Kuoen_US
dc.contributor.author張修豪zh_TW
dc.contributor.authorChang, Hsiu-Haoen_US
dc.contributor.other中興大學zh_TW
dc.date2011zh_TW
dc.date.accessioned2014-06-05T11:43:30Z-
dc.date.available2014-06-05T11:43:30Z-
dc.identifierU0005-2508201011152200zh_TW
dc.identifier.citation1. Alam, M. M. and Zhou, Y., Strouhal numbers,forces and flow structures around two tandem cylinders of different diameters. Journal of Fluids and Structures, 2008. 24, pp. 505-526. 2. Auteri, F., Belan, M., Gibertini, G., and Grassi, D., Normal flat plates in tandem: An experimental investigation. Journal of Wind Engineering and Industrial Aerodynamics 2008. 96, pp. 872-879. 3. Bloor, M. S., The transition to turbulence in the wake of a circular cylinder. Journal of Fluid Mechanics, 1964. 19, pp. 290-304. 4. Igarashi, T., Characteristics of a flow around two circular cylinders of dofferent diameters arranged in tandem. Bulletin of the JSME, 1982. 25, pp. 349-357. 5. Kiya, M., Vortex shedding from two circular cylinders in staggered arrangement. Journal of Fluids Engineering, Transactions of the ASME, 1980. 102, pp. 166-173. 6. Ko, N. W. M., Wong, P. T. Y., and Leung, R. C. K., interaction of flow structures within bistable flow behind two circular cylinders of different diameters. Experimental Thermal and Fluid Science, 1996. 12, pp. 33-44. 7. Spivack, H. M., Vortex frequency and flow pattern in the wake of two parallel cylinder at naried spacing normal to an air stream. Aeronautical Sciences, 1946. 13, pp. 289-297. 8. Sumner, D., Price, S. J., PA, Iuml, and DOUSSIS, M. P., Flow-pattern identification for two staggered circular cylinders in cross-flow. Journal of Fluid Mechanics, 2000. 411, pp. 263-303. 9. Sumner, D., Richards, M. D., and Akosile, O. O., Strouhal number data for two staggered circular cylinders. Journal of Wind Engineering and Industrial Aerodynamics, 2008. 96, pp. 856-871. 10. Sumner, D., Wong, S. S. T., and Price, S. J., Fluid behaviour of side-by-side circular cylinders in steady cross-flow. Journal of Fluids and Structures, 1999. 13, pp. 309-338. 11. Wang, Z. J. and Zhou, Y., Vortex interactions in a two side-by-side cylinder near-wake. International Journal of Heat and Fluid Flow, 2005. 26, pp. 362-377. 12. Xu, S. J., Zhou, Y., and So, R. M. C., Reynolds number effects on the flow structure behind two side-by-side cylinder. Physics of Fluids, 2003. 15, pp. 1214-1219. 13. Zdravkovich, M. M., Review of flow interference between two circular cylinder in various arrangements.ASME Trans. Journal of Fluid Engineering, 1977. 26, pp. 618-633. 14. Zhou, Y. and Antonia, R. A., Effect of initial condition on structures in a turbulent near-wake. AIAA, 1994, pp. 1207-1213. 15. Zhou, Y. and Antonia, R. A., The effect of Reynolds number on a turbulent far-wake. Experiments in Fluids, 1998. 25, pp. 118-125. 16. ProVISION-XS User Manual, 2007. 17. 李宜儒, 運用小波方法分析並列雙圓柱尾流的長時間特性,中興大學 機械工程系 碩士論文,2006. 18. 廖景晨, 縱列圓柱-垂直平板間自激式振盪流場及下游水平板影響之實驗研究,中興大學 機械工程系 碩士論文,2009.zh_TW
dc.identifier.urihttp://hdl.handle.net/11455/2536-
dc.description.abstract本論文係研究具有兩倍直徑比之錯列雙圓柱下游流場之特性。雷諾數為1000時,錯列圓柱之中心間距在L/D=-2.0~1.0,T/D=1.0~2.0間變化。本實驗在低速循環式水槽中進行,以流場可視化作定性的觀察,再以雷射都普勒測速儀、質點影像測速系統作定量量測。並且運用小波分析以及環流量之變化,進一步了解具有兩倍直徑比之錯列雙圓柱下游流場之特性。實驗結果顯示:當錯列雙圓柱之中心間距較大時,大、小圓柱下游流場會形成兩對平行流向的渦街結構,其特徵頻率會趨近單一大、小圓柱者。當錯列雙圓柱之中心間距逐漸縮小,因大、小圓柱相對位置之不同,導致間隙流偏斜,形成模式1(Mode 1)與模式2(Mode 2)。當錯列雙圓柱之中心間距持續縮小時,寬域兩側形成剪力層特性,往下游演化過程中,可能會產生渦漩配對,其特徵頻率會逐漸變為近域特徵頻率的 倍,在更下游處則會形成一交替脫離的渦列;而窄域下游尾流則形成交替脫離的渦漩,其特徵頻率為高頻。當錯列雙圓柱之中心間距很小時,圓柱後方為寬域者,其外側之剪力層特性更明顯,僅在下游處形成無渦漩配對低頻的渦漩結構;圓柱後方為窄域者,在下游則先形成尾流特性。因三個渦漩形成的位置相當靠近,且渦漩形成長度很短,尾流特性之範圍很小,往下游演化過程中,陸續結合成一大尺度渦漩,其特徵頻率接近寬域外側剪力層者。在更下游處,與寬域外側之渦漩逐漸形成一個具有低頻的交替脫離渦漩列。zh_TW
dc.description.abstractThis study investigates the flow characteristics behind two staggered circular cylinders of diameter ratio two by experiment. The Reynolds number was kept 1000; and the streamwise and the transverse distances between the cylinder centers are located within the range of L/D=-2.0~1.0, T/D=1.0~2.0, respectively. All the experiments are performed in a low-speed recirculation water channel. The qualitative flow structures were observed by the dye flow visualization technique, and the quantitative velocity measurements were performed by the FLDV system. Some of the data are analyzed by the cross Wavelet transformation. Further, using the PIV system acquire the instantaneous vorticity and circulation of the related flow structures between these two staggered cylinders. Some important results are summarized as follows: When the distances (either streamwise or transverse) between the centers of staggered circular cylinders are large, the gap flow becomes nearly parallel; and the flow structure forms two parallel vortex streets shed alternatively in the downstream regions. In these cases, the characteristic frequencies are very close to those behind a single large and small circular cylinders. While the center spacing decreases, two modes of the flow structures are defined, mode 1 and 2, depending upon the phases of the gap vortices. At smaller gap, the flow of the narrow wake shows the characteristic with much higher frequency. However, the two outer shear layers of the wide wake sometimes show the pairing of the vortices, forming the 1/2 subharmonic component of the frequency of the narrow wake during the evolution process. In the further downstream region, the wake character may be formed with a low characteristic frequency. Further decrease of the distances (L/D and T/D) between two cylinders much shorter vortex formation length is found behind the narrow wake. On the other hand, the outer shear layers of the wide wake can only form the vortices of much low frequency due to instability. The widely spread shear layer interacts strongly with the vortex street of the narrow wake and will form a much large scale vortex in the farther downstream region. The characteristic frequencies of the two shear layers are close, and may form a vortex street shed alternatively with much lower frequency.en_US
dc.description.tableofcontents目錄..............................................................................................................Ⅰ 圖表目錄………………………………………………………….…....….Ⅲ 符號說明…………………………………………………………….….…Ⅸ 第一章 緒論 1-1 研究動機…………………………………………………...…...1 1-2 文獻回顧…………………………………………………...…...1 第二章 實驗設備及模型 2-1 低速循環式水槽……………………………….………..…...…8 2-2 光纖雷射都普勒雷射測速儀………………………….…...…..8 2-3 PIV量測系統……………………………….………….…..…....9 2-4 精密三維移行機構……………………………….…….….……9 2-5 實驗模型………………………………….…………….….……9 第三章 實驗方法與數據分析 3-1 流場可視化………………………………………….…….....…10 3-2 LDV速度量測步驟與原理……………………...………......…10 3-3 PIV質點影像測速儀量測原理……………………….….....….11 3-4 小波分析………………………………………………….....….11 3-5 渦度與環流量……………………………………………....…..12 3-6 訊號分析與測試………………………………………………..13 第四章 實驗結果與討論 14 4-1 不同中心間距之錯列圓柱流場結構特性 14 4-1-1 模式1(Mode 1)的流場結構可視化…………….…. 15 4-1-2 過度流場結構可視化………………………...……......16 4-1-3 模式2(Mode 2)的流場結構可視化...............................16 4-1-4 流場可視化之小結…………………...………………..17 4-2 尾流流場結構之特徵頻率分析 18 4-2-1 模型1(Mode 1)之流場結構…………….……..…........18 4-2-2 模型2(Mode 2)之流場結構…………….….……....….19 4-3 全域流場觀察 4-3-1 模型1(Mode 1)之渦度分佈及演化………………..…20 4-3-2 模型2(Mode 2)之渦度分佈及演化………………...…21 4-4 環流量之變化 4-4-1 單一圓柱之渦漩結構之演化…………...………...….....22 4-4-2 模型1(Mode 1)之渦漩結構演化…………………….…22 4-4-3 模型2(Mode 2)之渦漩結構演化………………….……24 4-5 錯列圓柱下游流場長時間平均及擾動能量之分析 4-5-1 交替脫離之渦列尾流特性……………...……….….......26 4-5-2 剪力層流場特性……………...........................26 4-5-3 錯列圓柱流場……………………………26 4-6 直徑比為2錯列圓柱下游之流場特徵與分類 4-6-1 模型1之特徵………………………………………34 4-6-2 模型2之特徵………………………………………35 第五章 結論…………………………………………………37 文獻回顧……………………………………………………………39 圖1-1:史卓赫數隨著淨間距變化的分布[7]。………………...............40 圖1-2:(a) 渦漩A、B、C在下游處會結合在一起。 (b) 渦漩B、C會結合在一起;而渦漩A 會在下游處潰散[11]。…………………………………………40 圖1-3:固定雷諾數,改變中心間距其無因次頻率的分布[12]。……41 圖1-4:壓力波動係數圖可觀察聯接、跳動、雙穩態之區域[4]。......41 圖1-5:不同雷諾數對下游圓柱阻力係數之影響[4]。……………......42 圖1-6:不同圓柱間距所產生不同的流場結構[4]。………………......42 圖1-7:錯列圓柱排列之中心間距以及夾角示意圖[5]。…………......43 圖1-8:在雷諾數為850~1900,以兩圓柱間的間距變化為縱軸, 以兩圓夾角變化為橫軸,又可細分出九種不同流場[8]。…….43 圖1-9:觀察直徑比為2:1的雙圓柱下游的流場結構,可分成 兩種基本模式:間隙流偏向大圓柱,且小圓柱後方的 回流區比大圓柱後方的回流區大為模式1;相反則為 模式2[6]。……………………………………………….…..…..44 圖2-1:低速循環式水槽示意圖:(a)上視圖;(b)側視圖。…….....45 圖2-2:FLDV速度量測系統示意圖。……………………………........46 圖2-3:PIV量測系統示意圖。……………………………………...….47 圖2-4:錯列圓柱擺設模型示意配置圖。………………………….......48 圖3-1:FLDV量測系統之干涉條紋示意圖。…………………………49 圖3-2:PIV測速系統分析向量示意圖。………………………….…..50 圖3-3:小波母函數經調整尺度因子及平移因子後的示意圖: (a)尺度因子a為1,平移因子b=0;(b)尺度因子 a為2,平移因子b=0;(c)尺度因子a為0.5,平移因 子b=0;(d)平移因子b=1。………………………………..…51 圖3-4:(a)模擬速度訊號U(t)=5sin(2π×2×t), 0<t≦10; 5sin(2π×1×t), 0<t≦20(b)能量密度頻譜圖 (c)小波轉換後的等高線圖。………….…………………....52 圖4-1:直徑比為2之錯列圓柱下游流場結構之演化特徵 與分類。…………………………………………………………53 圖4-2:錯列雙圓柱之中心間距為L/D=-2.5,T/D=3.0之流 場可視化照片。……………………………………..………….54 圖4-3:錯列雙圓柱之中心間距為L/D=-1.0,T/D=1.5之流 場可視化照片。………………………………………..…….….55 圖4-4:錯列雙圓柱之中心間距為L/D=-2.0,T/D=1.0之流 場可視化照片。………….……………….…………………….56 圖4-5:錯列雙圓柱之中心間距為L/D=-0.5,T/D=1.0之流 場可視化照片。.…………………………………………….….58 圖4-6:錯列雙圓柱之中心間距為L/D=-0.5,T/D=2.0時, 間隙流產生上、下不穩定擺動的流場現象。………………….59 圖4-7:錯列雙圓柱之中心間距為L/D=0.0,T/D=1.0之流 場可視化照片。……………………………………….……….60 圖4-8:錯列雙圓柱之中心間距為L/D=1.0,T/D=2.0之流 場可視化照片。…………………………………………..……61 圖4-9:錯列雙圓柱之中心間距為L/D=1.0,T/D=4.0之流 場可視化照片。………………………………………….…….62 圖4-10:當雷諾數為1000時,單一圓柱其尾流擾動速度之 能量密度頻譜圖;(a)單一圓柱直徑為D,(b)單一 圓柱直徑為d。……………………………………………....63 圖4-11:當雷諾數為1000,錯列雙圓柱之中心間距為 L/D=-2.0,T/D=2.0時,其尾流擾動速度之 能量密度頻譜圖。……………………………..……………..64 圖4-12 當雷諾數為1000時,錯列雙圓柱之中心間距為 L/D=-2.0,T/D=1.0,其尾流擾動速度之能量密 度頻譜圖。………………………………………………..…….65 圖4-13:當雷諾數為1000,錯列雙圓柱之中心間距為 L/D=-1.0,T/D=1.5時,其尾流擾動速度之能量 密度頻譜圖。……………………………………….………......67 圖4-14:當雷諾數為1000,錯列雙圓柱之中心間距為 L/D=-0.5,T/D=1.0時,其尾流擾動速度之能量 密度頻譜圖。……………………………………………..…….69 圖4-15:當雷諾數為1000,錯列雙圓柱之中心間距為 L/D=0.0,T/D=1.0時,其流擾動速度之能量密 度頻譜圖。…………………………………………………..…71 圖4-16:當雷諾數為1000,錯列雙圓柱之中心間距為 L/D=1.0,T/D=2.0時,其尾流擾動速度之能量 密度頻譜圖。…………..……………………………………....73 圖4-17:當雷諾數為1000時,錯列雙圓柱之中心間距 為L/D=-2.0,T/D=2.0下游渦漩A與B交替渦漩 脫離,以及渦漩C與渦漩D交替渦漩脫離,形成 兩對渦街的尾流流場渦度圖。……………………………..…74 圖4-18:當雷諾數為1000時,錯列雙圓柱之中心間距 為L/D=-2.0,T/D=1.0下游渦漩A+B+C之結合, 與獨自脫離D渦漩流場渦度圖。.……………………………75 圖4-19:當雷諾數為1000時,錯列雙圓柱之中心間距為 L/D=-1.0,T/D=1.5下游渦漩A與B渦漩脫離,以 及渦漩Cn配對與渦漩Dn配對渦漩脫離的尾流流場 渦度圖。…………………………………………………….….76 圖4-20:當雷諾數1000時,錯列雙圓柱之中心間距為 L/D=-0.5,T/D=1.0,下游渦漩A+B+C之結合 與獨自脫離的Dn配對流場渦度圖。…………………….……77 圖4-21:當雷諾數為1000時,錯列雙圓柱中心間距為 L/D=0.0,T/D=1.0,下游渦漩B+C+D之結合與 獨自脫離的Dn配對流場渦度圖。………………………….…78 圖4-22:當雷諾數為1000時,錯列雙圓柱之中心間距 為L/D=1.0,T/D=2.0,下游渦漩A與B交替渦 漩脫離以及渦漩C與渦漩D交替渦漩脫離形成 兩對渦街的尾流流場渦度圖。…………………………………79 圖4-23:當雷諾數為1000時,單一大(左)、小(右)圓柱 後方尾流之(a)環流量的週期性變化圖,(b)環流量 的能量密度頻譜圖,(c)利用LDV量測之能量密度 頻譜圖。………………………..……………………………….80 圖4-24:當雷諾數為1000時,間距為L/D=-2.0,T/D=2.0, 錯列圓柱後方上半部尾流之環流量圖(左);以及環流量 的能量密度頻譜圖(右)。下半部尾流之環流量圖;以及 環流量的能量密度頻譜圖。…………………………….……..81 圖4-25:當雷諾數為1000時,間距為L/D=-2.0,T/D=1.0, (a-c)錯列圓柱後方上半部尾流之環流量圖(左); 以及環流量的能量密度頻譜圖(右)。(d-f) 錯列圓柱 後方下半部尾流之環流量圖;以及環流量的能量密 度頻譜圖。……………………………………………....…..83 圖4-26:當雷諾數為1000時,間距為L/D=-1.0,T/D=1.5, 錯列圓柱後方上半部尾流之環流量圖(左);以及 環流量的能量密度頻譜圖(右)。下半部尾流之環 流量圖;以及環流量的能量密度頻譜圖。………………..86 圖4-27:當雷諾數為1000時,間距為L/D=-0.5,T/D=1.0 錯列圓柱後方上半部尾流之環流量圖(左) ;以及 環流量的能量密度頻譜圖(右)。下半部尾流之 環流量圖;以及環流量的能量密度頻譜圖。……………..88 圖4-28:當雷諾數為1000時,間距為L/D=0.0,T/D=1.0, 錯列圓柱後方上半部尾流之環流量圖(左);以及 環流量的能量密度頻譜圖(右)。下半部尾流之 環流量圖;以及環流量的能量密度頻譜圖。……………...90 圖4-29:當雷諾數為1000時,間距為L/D=1.0,T/D=2.0 ,錯列圓柱後方上半部尾流之環流量圖(左); 以及環流量的能量密度頻譜圖(右)。下半部 尾流之環流量圖;以及環流量的能量密度頻譜圖。.……..92 圖4-30;將環流量進行FFT計算得到近、遠域之特徵頻率 除以單一大圓柱之尾流特徵頻率(a)T/D為1.0時 大小圓柱近、遠域之無因次化特徵頻率(b) T/D為 1.5時大小圓柱近、遠域之無因次化特徵頻率(c) T/D為2.0時大小圓柱近、遠域之無因次化特徵頻率。..…94 圖4-31:雷諾數為1000時,(a-b)單一圓柱下游的尾流擾動速 度的Urms與Vrms等高圖,(c)沿著Y軸在X/D=3.0取出 Urms與Vrms之分布圖,(d-e) 沿著X軸在Y/D=0.5、 -0.5取出Urms與Vrms之分布圖。………………….………..95 圖4-32:雷諾數為1000時,(a-b)剪力層流場下游的尾流 擾動速度的Urms與Vrms等高圖,(c-d)沿著Y軸 在X/D=2.0、3.0取出Urms與Vrms之分布圖,(e) 沿著X軸在Y/D=1.4取出Urms與Vrms之分布圖。.……….96 圖4-33:雷諾數為1000時,(a-b)L/D=-2.0,T/D=2.0錯列 圓柱下游的尾流擾動速度的Urms與Vrms等高圖, (c-d)沿著Y軸在X/D=0.8、1.8取出Urms與Vrms之 分布圖,(e-f) 沿著X軸在Y/D=-0.2、2.5取出 Urms與Vrms之分布圖。……………………………..………97 圖4-34:(a-b)沿著Y軸在X/D=1.0、3.0取出特徵頻率為 0.65Hz、0.9Hz能量強度之分佈圖,(c-d) 沿著 X軸在Y/D=-0.5、2.8取出特徵頻率為0.65Hz、 0.9Hz能量強度之分佈圖。………………………….…..…..98 圖4-35:雷諾數為1000時,(a-b)L/D=-2.0,T/D=1.0錯列 圓柱下游的尾流擾動速度的Urms與Vrms等高圖, (c-d)沿著Y軸在X/D=0.4、2.6取出Urms與Vrms 之分布圖,(e-f) 沿著X軸在Y/D=-0.2、1.8取出 Urms與Vrms之分布圖。…………………………..……..…..99 圖4-36:雷諾數為1000時,在(a-b)L/D=-1.0,T/D=1.5 錯列圓柱下游的尾流擾動速度的Urms與Vrms 等高圖,(c)沿著Y軸在X/D=1.8取出Urms與 Vrms之分布圖,(d-e) 沿著X軸在Y/D=-0.4、 2.0取出Urms與Vrms之分布圖。………………………...….100 圖4-37:(a-b)沿著Y軸在X/D=1.0、4.0取出特徵頻率為 0.45Hz、0.9Hz能量強度之分佈圖,(c-d) 沿著 X軸在Y/D=-0.5、2.4取出特徵頻率為0.45Hz、 0.9Hz能量強度之分佈圖。………………………………..101 圖4-38:雷諾數為1000時,(a-b)L/D=-0.5,T/D=1.0錯列圓柱下游的尾流擾動速度的Urms與Vrms等高圖,(c-d)沿著Y軸在X/D=1.2、2.2取出Urms與Vrms之分布圖,(e-f) 沿著X軸在Y/D=-0.4、1.8取出Urms與Vrms之分布圖,(g-h) 下游遠域的尾流擾動速度的Urms與Vrms等高圖,(i) 沿著Y軸在X/D=7.4取出Urms與Vrms之分布圖,(j-k) 在遠域沿著X軸在Y/D=-0.4、1.4取出Urms與Vrms之分布圖。……………………………………………………….…..102 圖4-39:(a-b)沿著Y軸在X/D=1.5、4.5取出特徵頻率為0.3Hz、0.6Hz、1.2Hz能量強度之分佈圖,(c-d)沿著X軸在Y/D=1.6、-0.5取出特徵頻率為0.3Hz、0.6Hz、1.2Hz能量強度之分佈圖。.……………..….104 圖4-40:雷諾數為1000時,(a-b)L/D=0.0,T/D=1.0錯列圓柱下游的尾流擾動速度的Urms與Vrms等高圖,(c-d)沿著Y軸在X/D=1.8、5.0取出Urms與Vrms之分布圖,(e-f) 沿著X軸在Y/D=-0.5、1.8取出Urms與Vrms之分布圖,(g-h) 下游遠域的尾流擾動速度的Urms與Vrms等高圖,(i) 沿著Y軸在X/D=7.4取出Urms與Vrms之分布圖,(j-k) 在遠域沿著X軸在Y/D=-0.4、1.8取出Urms與Vrms之分布圖。………..………105 圖4-41:(a-b)沿著Y軸在X/D=2.0、5.0取出特徵頻率為0.35Hz、0.7Hz、1.4Hz能量強度之分佈圖,(c-d) 沿著X軸在Y/D=1.8、-0.5取出特徵頻率為0.35Hz、0.7Hz、1.4Hz能量強度之分佈圖。………………………………………………………….107 圖4-42:雷諾數為1000時,(a-b)L/D=1.0,T/D=2.0錯列圓柱下游的尾流擾動速度的Urms與Vrms等高圖,(c-d)沿著Y軸在X/D=1.8、4.4取出Urms與Vrms之分布圖,(e-f) 沿著X軸在Y/D=-0.2、2.5取出Urms與Vrms之分布圖。………………….…………….……108 圖4-43:(a-b)沿著Y軸在X/D=2.0、5.0取出特徵頻率為0.5Hz、1.3Hz能量強度之分佈圖,(c-d) 沿著X軸在Y/D=2.4、-0.5取出特徵頻率為0.5Hz、1.3Hz能量強度之分佈圖。…………………………..….….109 圖4-44:錯列圓柱下游直徑比為2之流場結構特徵之分類圖。…….110zh_TW
dc.language.isoen_USzh_TW
dc.publisher機械工程學系所zh_TW
dc.relation.urihttp://www.airitilibrary.com/Publication/alDetailedMesh1?DocID=U0005-2508201011152200en_US
dc.subjectstaggereden_US
dc.subject錯列zh_TW
dc.subjectbluff bodyen_US
dc.subjectgap flowen_US
dc.subjecttwo circular cylindersen_US
dc.subject鈍體zh_TW
dc.subject間隙流zh_TW
dc.subject雙圓柱zh_TW
dc.title具有兩倍直徑比之錯列雙圓柱下游流場特性之研究zh_TW
dc.titleStudy on Flow Characteristics behind Two Staggered Circular Cylinder of Diameter Ratio Twoen_US
dc.typeThesis and Dissertationzh_TW
Appears in Collections:機械工程學系所
文件中的檔案:

取得全文請前往華藝線上圖書館



Items in DSpace are protected by copyright, with all rights reserved, unless otherwise indicated.