邊界層流中雙棟鄰近建築物之受風效應研究
73
0
0
全文
(2) 095301070000G2027. 邊界層流中雙棟鄰近建築物之受 風效應研究. 研究主持人：黎益肇. 內政部建築研究所自行研究報告 中華民國 95 年 12 月.
(3) 目次. 目次 目次························································································· I 表次······················································································III 圖次······················································································ IV 摘. 要·················································································VII. ABSTRACT ········································································ IX 第一章 緒. 論 ······································································ 1. 第一節 研究緣起與背景 ·················································1 第二章 文獻回顧 ·································································· 5 第一節 相關研究 ·····························································5 第二節 耐風設計規範 ·····················································7 第三章 研究方法 ································································ 11 第一節 實驗方法 ···························································11 第二節 量測分析項目 ···················································20 第三節 數值方法 ···························································21 第四章 結果討論 ································································ 25 第一節 第二節 第三節 第四節 第五節. 模式驗證 ···························································25 雷諾數效應 ·······················································27 流場結果 ···························································32 風壓係數 ···························································39 風力係數 ···························································44. 第五章 結論與建議 ···························································· 44. I.
(4) 邊界層流中雙棟鄰近建築物之受風效應研究. 第一節 結論 ···································································44 第二節 建議 ···································································44 附錄一 期中簡報會議紀錄與回應 ···································· 44 附錄二 電子式壓力掃瞄器操作說明 ································ 44 參考書目·············································································· 44.
(5) 表次. 表次 表 2-1 地況相關參數 ························· 8 表 2-2 牆之平均外風壓係數 ··················· 9 表 3-1 實驗用流場條件 ······················ 20 表 4-1 實驗與數值結果比較表 ················ 26 表 4-2 平均風力係數表 ······················ 44 表 4-3 擾動風力係數表 ······················ 44. III.
(6) 邊界層流中雙棟鄰近建築物之受風效應研究. 圖次 圖 3-1 邊界層流中雙建築物模型 ·············· 12 圖 3-2 模型孔位配置圖 ······················ 12 圖 3-3 電子式壓力掃瞄器 ···················· 13 圖 3-4 熱線探針 ···························· 14 圖 3-5 熱線測速儀 ·························· 14 圖 3-6 皮托管 ······························ 16 圖 3-7 薄膜式壓力計 ························ 16 圖 3-8 資料擷取系統 ························ 17 圖 3-9 研究個案示意圖 ······················ 18 圖 3-10 來流邊界層特性剖面(地況 B) ········· 19 圖 3-11 來流邊界層特性剖面(地況 A) ········· 19 圖 3-12 透式邊界條件說明圖 ················· 23 圖 3-13 不透式邊界條件說明圖 ··············· 24 圖 4-1 邊界層流通過建築物之流場示意圖 ······ 26 圖 4-2 柱體中心線表面風壓比較圖 ············ 27 圖 4-3 平均風壓係數分布圖(Re=9.1×10 4 ) ······ 28 圖 4-4 平均風壓係數分布圖(Re=1.35×10 5 ) ····· 28 圖 4-5 平均風壓係數分布圖(Re=1.89×10 5 ) ····· 29. IV.
(7) 圖次. 圖 4-6 擾動風壓係數分布圖(Re=9.1×10 4 ) ······ 29 圖 4-7 擾動風壓係數分布圖(Re=1.35×10 5 ) ····· 30 圖 4-8 擾動風壓係數分布圖(Re=1.89×10 5 ) ····· 30 圖 4-9 平均風壓係數分布圖(數值結果) ········ 31 圖 4-10 擾動風壓係數分布圖(數值結果) ······· 32 圖 4-11 速度等值圖(S/D=0.5、y/H=1/3) ······· 33 圖 4-12 速度等值圖(S/D=0.5、y/H=2/3) ······· 33 圖 4-13 速度等值圖(S/D=1、y/H=1/3) ········· 34 圖 4-14 速度等值圖(S/D=1、y/H=2/3) ········· 34 圖 4-15 速度等值圖(S/D=2、y/H=1/3) ········· 35 圖 4-16 速度等值圖(S/D=2、y/H=2/3) ········· 35 圖 4-17 速度等值圖(S/D=3、y/H=1/3) ········· 36 圖 4-18 速度等值圖(S/D=3、y/H=2/3) ········· 36 圖 4-19 速度等值圖(S/D=4、y/H=1/3) ········· 37 圖 4-20 速度等值圖(S/D=4、y/H=2/3) ········· 37 圖 4-21 速度等值圖(S/D=5、y/H=1/3) ········· 38 圖 4-22 速度等值圖(S/D=5、y/H=2/3) ········· 38 圖 4-23 平均風壓係數分布圖(S/D=0.5) ········ 40 圖 4-24 平均風壓係數分布圖(S/D=1) ·········· 40. V.
(8) 邊界層流中雙棟鄰近建築物之受風效應研究. 圖 4-25 平均風壓係數分布圖(S/D=2) ·········· 41 圖 4-26 平均風壓係數分布圖(S/D=3) ·········· 41 圖 4-27 平均風壓係數分布圖(S/D=4) ·········· 42 圖 4-28 平均風壓係數分布圖(S/D=5) ·········· 42 圖 4-29 擾動風壓係數分布圖(S/D=0.5) ········ 43 圖 4-30 擾動風壓係數分布圖(S/D=1) ·········· 44 圖 4-31 擾動風壓係數分布圖(S/D=2) ·········· 44 圖 4-32 擾動風壓係數分布圖(S/D=3) ·········· 44 圖 4-33 擾動風壓係數分布圖(S/D=4) ·········· 44 圖 4-34 擾動風壓係數分布圖(S/D=5) ·········· 44. VI.
(9) 摘. 摘. 要. 要. 關鍵詞：邊界層、風洞實驗、表面風壓、數值模擬 一、研究緣起 單 一 建 築 物 於 風 域 中 之 流 場 變 化，常 為 研 究 建 築 物 風 荷 重 之 基 本 情 況，前 人 研 究 亦 多 著 眼 於 此。然 就 實 際 情 況 而 言，都 會 建 築 皆 為 毗 臨 而 立 而 非 單 一 存 在，甚 至 於 集 合 式 住 宅 之 設 計 亦 常 出 現 此 情 形。因 此，相鄰建築物間於風場中之交互作用，便為本研究之主要課題。 與 風 洞 試 驗 比 較，應 用 數 值 模 擬 的 方 法 解 析 此 類 問 題 則 可 較 經 濟 地 獲 取 更 為 完 備 之 整 場 設 計 資 料。本 研 究 以 數 值 方 法 為 主，並 同 步 進 行風洞實驗，對本數值模式之準確性予以驗證。 二、研究方法 為 深 入 瞭 解 風 域 中 雙 矩 形 建 物 受 鄰 近 風 效 應 之 極 端 情 況，本 研 究 針對邊界層來流中雙棟建築物探討其受風效應。藉由橫列 (Side-by-side)之 方 形 建 築 物 間 距 改 變，量 測 建 築 物 的 表 面 風 壓、阻 昇 力 以 及 流 場 特 性，以 探 討 建 物 表 面 氣 動 力 作 用 之 機 制，並 以 此 試 驗 結 果 作 為 提 昇 數 值 計 算 準 確 性 之 比 對 基 準。主 要 工 作 係 以 風 洞 模 型試 驗 量 測 與 非 恆 定 紊 流 流 場 數 值 模 擬 發 展 二 者 並 行，結 合 兩 種 方 式 獲 致 之結果，交叉分析建築物受風效應與相應紊流風場間相互之變化機 制，據以獲得綜合性的結論。 三、重要發現 (1) 數 值 結 果 與 實 驗 數 據 有 相 當 不 錯 的 吻 合 度 ， 可 作 為 建 築 物 受 風 壓 初步評估之工具。且在檢視流場與結構體受風力之相關性時，為 一有利之工具。 (2) 平 均 阻 力 係 數 在 間 距 比 0.5 時 有 最 大 值 1.015， 明 顯 高 於 單 柱 之. VII.
(10) 邊界層流中雙棟鄰近建築物之受風效應研究. 0.828， 故 在 此 情 況 時 需 慎 重 考 慮 其 風 力 設 計 參 數 。 (3) 當 間 距 比 大 於 4 之 後 ， 建 築 物 各 項 風 力 特 性 已 近 似 單 柱 之 情 形 。 (4) 間 距 比 從 0.5 變 化 至 5 時 ， 迎 風 面 最 大 壓 力 逐 漸 由 由 上 角 轉 移 至 中心線，且流場分布亦逐漸分離成單一柱體之型態。. 四、主要建議事項 (1) 本 所 於 本 年 度 添 購 一 套 SUN V440 FDS 伺 服 器 主 機 ， 對 於 數 值 計 算 能量有實質之幫助。未來將近一步添購相關運算軟體，以利本模 式使用。 (2) 對 於 氣 動 力 實 驗 所 採 用 之 電 子 式 壓 力 掃 描 系 統 ， 需 建 立 一 套 自 行 校正程序，以確保實驗數據之準確性。 (3) 在 規 劃 氣 動 力 模 型 試 驗 時，銳 緣 附 近 壓 力 孔 位 之 解 析 度 應 提 高 (或 壓 力 變 化 較 大 之 位 置 )， 有 助 於 繪 製 壓 力 分 布 圖 或 尋 找 局 部 極 端 值。 (4) 後 續 應 規 劃 更 多 雙 棟 建 築 物 之 排 列 組 合 狀 況 ， 以 確 立 邊 界 層 流 中 雙棟建築物之風力特性。. VIII.
(11) 摘. 要. ABSTRACT Keywords: boundary layer、 wind tunnel experimental、 numerical simulation. Many. investigators. have. discussed. the. wind. field. distinction of single building, which must often be taken into study the wind load of structure. In fact, the architectures not stand alone in the city, even in the design of aggregating residence. Therefore, in this research, the major object is the wind field interaction behavior of neighbor building。 Compared to wind tunnel tests, the application of numerical simulations can yield complete design data in a more economical way. Besides, since the resulting wind data are represented in time series, it becomes more applicable to the building designs. For deeply realizing of the peak situation of wind around the twin rectangular building, this article investigates the wind effect of twin building in the boundary layer incoming flow. According to the conversion of distance between two square buildings in side-by-side arrangement, we measure the surface pressure, drag force, lift force and flow characteristic, investigate. building. surface. aerodynamic. behavior.. The. accuracy of the numerical predictions will be guided and improved based on the measurement results.. IX.
(12)
(13) 第一章 緒. 第一章 緒. 論. 論. 第一節 研究緣起與背景. 壹、研究緣起 矩形建築物在風域中紊流流場之特性分析一向是建築研究領域 中 重 要 的 一 環。為 了 要 獲 取 高 層 建 築 物 承 受 之 時 變 風 荷 重 以 推 算 其動 態 反 應，進 而 決 定 建 築 物 住 民 舒 適 性 之 設 計 指 標，現 有 方 法 往 往 採用 風 洞 試 驗，或 經 由 實 場 監 測 之 整 合 資 料，以 頻 率 域 表 徵 之 風 函 數 為輸 入 依 據 進 行 結 構 之 動 態 反 應 評 估。但 當 建 築 物 之 造 型 特 殊 或 建 物 鄰 近 區 之 地 況 複 雜 時，此 方 法 之 應 用 將 產 生 誤 差，甚 至 導 致 不 保 守 之 設計 結果。 單 一 高 層 建 築 物 於 風 域 中 之 流 場 變 化，常 為 研 究 建 築 物 風 荷 重 之 基 本 情 況，前 人 研 究 亦 多 著 眼 於 此。然 就 實 際 情 況 而 言，都 會 建 築皆 為毗臨而立而非單一存在，甚至於集合式住宅之設計亦常出現此情 形。因 此，相 鄰 建 築 物 間 於 風 場 中 之 交 互 作 用，便 為 本 研 究 之 主 要課 題。 與 風 洞 試 驗 比 較，應 用 數 值 模 擬 的 方 法 解 析 此 類 問 題 則 可 較 經 濟 地 獲 取 更 為 完 備 之 整 場 設 計 資 料。本 研 究 以 數 值 方 法 為 主，並 同 步 進 行風洞實驗，對本數值模式之準確性予以驗證。. 貳、研究目的 本研究之目標有如下四項： (1) 系 統 化 地 探 討 來 流 情 況 與 建 築 物 幾 何 條 件 之 改 變 對 受 風 效 應 與 流 場機制之影響，以提供工程分析之通則。 (2) 細 部 檢 視 建 築 物 幾 何 條 件 與 鄰 近 風 場 間 之 相 互 關 係 ， 藉 系 統 化 模. 1.
(14) 邊界層流中雙棟鄰近建築物之受風效應研究. 擬結果獲取相關設計之重要參考資料。 (3) 數 值 計 算 模 式 預 測 精 度 之 提 昇 ， 將 可 提 供 更 為 精 確 的 風 場 與 結 構 受風效應時序列資料，除了可以成為高層建築風力設計之分析工 具外，更可以配合模型試驗量測，達到相輔相成之效果。 (4) 增 進 本 所 同 仁 相 關 實 驗 設 計 與 安 裝 量 測 流 程 之 熟 析 度 。 本 研 究 針 對 邊 界 層 來 流 的 情 況 下，使 用 建 築 物 橫 向 排 列 方 式，並 藉 由 其 來 流 風 況 與 間 距 比 改 變，分 析 建 築 物 各 項 氣 動 力 參 數 與 鄰 近建 物之相關性。. 參、研究方法 本 研 究 之 工 作 包 括 風 洞 模 型 試 驗 與 數 值 模 擬 計 算 兩 個 部 份，茲 分 述如後： (1) 風 洞 模 型 試 驗 1.來 流 風 況 除 均 勻 來 流 外，本 研 究 擬 於 針 對 代 表 平 坦 市 郊 及 大 型 都 會 中 心 兩 種 地 形 的 大 氣 邊 界 層 形 式 之 來 流 風 場 中 進 行 實 驗 量 測，目 標 邊 界 層平 均 風 速 剖 面 以 指 數 律 表 示 之 指 數 分 別 約 為 0.21(地 況 C)及 0.32(地況 A)， 邊 界 層 厚 度 約 為 1.5 米 。 2. 建 物 模 型 建 築 模 型 以 壓 克 力 板 粘 合 製 作 ， 採 用 高 度 (H)120 公 分 、 厚 度 (D) 為 20 公 分 之 方 柱 ， 相 應 之 阻 塞 比 低 於 4%。 3.間 距 比 模 型 配 置 之 間 距 比 採 用 0.5、1、2、3、4、5(相 應 距 離 為 10、20、 40、 60、 80、 100 公 分 )， 並 與 單 柱 結 果 作 比 對 分 析 。 (2) 數 值 模 擬 研 究 中 流 場 之 計 算 係 採 用 微 可 壓 縮 流 (Weakly-compressible-flow Method [1])方 法 ， 以 模 擬 建 築 物 鄰 近. 2.
(15) 第一章 緒. 論. 區 域 之 三 維 (Three-dimensional)非 恆 定 (Unsteady)風 場 。 數 值 模 式 中 另 應 用 大 渦 模 擬 (Large Eddy Simulation)配 合 以 次 網 格 紊 流 模 型 (Subgrid-scale Turbulence Model)， 以 反 映 出 真 實 風 場 中 之 紊 流 特 性。 (3) 研 究 步 驟 1.風 洞 模 型 試 驗 ○ 1 試驗設施規劃與設計。 ○ 2 模型與試驗機構製作與建構。 ○ 3 模型試驗與量測。 2.數 值 模 擬 ○ 1 數值模式發展。 ○ 2 模式測試。 ○ 3 數值計算結果比對與驗證。 ○ 4 資料整合分析。 ○ 5 結論與成果報告。. 3.
(16) 邊界層流中雙棟鄰近建築物之受風效應研究. 4.
(17) 第二章 文獻回顧. 第二章 文獻回顧 單 一 高 層 建 築 物 於 風 域 中 之 流 場 變 化，常 為 研 究 建 築 物 風 荷 重 之 基 本 情 況，前 人 研 究 亦 多 著 眼 於 此。然 就 實 際 情 況 而 言，都 會 建 築皆 為毗臨而立而非單一存在，甚至於集合式住宅之設計亦常出現此情 形。因 此，相 鄰 建 築 物 間 於 風 場 中 之 交 互 作 用，便 為 本 研 究 之 主 要課 題。. 第一節 相關研究. 壹、二維研究 Okajima(1984)採 用 風 洞 量 測 作 用 於 橫 向 排 列 雙 方 柱 表 面 壓 力 分 佈及渦流脫離頻率，以探討高雷諾數下不同間距比對柱體昇阻力係 數、背壓係數與史特赫數的變化。實驗結果發現當兩柱間距微小時， 其 流 場 行 為 近 似 於 單 一 矩 柱 之 情 況；當 間 距 比 為 1 至 1.5 時，雙 方 柱 後之尾流呈現間歇性不穩定的偏向形成與雙柱對稱及偏向兩不同流 況，而 於 雙 柱 尾 流 對 稱 之 流 場 較 偏 向 有 較 小 之 背 壓 係 數，顯 示 其 受較 大 之 阻 力。而 隨 著 間 距 增 大，不 穩 定 偏 向 發 生 則 大 為 降 低 說 明 了 彼此 互制效應的減弱。 Kobayasi(1976)採 用 風 洞 量 測 均 勻 流 來 流 下 在 雷 諾 數 為 1.2~7.6 ×10 5 於 不 同 間 距 比 作 用 在 橫 列 雙 方 柱 之 昇 力 係 數 及 阻 力 係 數 分 佈。實 驗 結 果 顯 示 於 兩 方 柱 極 為 接 近 （ S/D≦ 0.1） 時 C D 有 最 大 值 ， 而 隨 著 間 距 比 之 增 加，方 柱 相 應 之 阻 力 係 數 呈 現 漸 減 之 趨 勢。此 外，兩 方柱 之昇力係數於所有間距比情況下皆不為零。 Sarode 等（ 1981）於 實 驗 研 究 中 發 現，紊 流 場 中 橫 向 排 列 之 雙 方 柱 其 平 均 阻 力 係 數 （ C D ） 皆 比 單 一 柱 體 時 為 高 。 而 當 間 距 比 (S/D)小. 5.
(18) 邊界層流中雙棟鄰近建築物之受風效應研究. 於 3 時 ， CD 值 隨 間 距 之 增 大 而 呈 現 上 昇 的 趨 勢 。 Sun 等 (1995)於 雷 諾 數 為 1.03×10 5 時，量 測 不 同 寬 深 比 之 雙 矩 柱 (含 雙 方 柱 )於 各 間 距 比 下 所 相 應 之 流 場 特 性。實 驗 中 發 現，當 2 個 方 柱 為 橫 向 排 列 時，在 微 小 間 距 比 (S/D＜ 0.1)其 壓 力 分 布 近 似 於 寬 深比 為 0.5 之 單 一 矩 柱。而 其 間 距 比 小 於 2 時，雙 方 柱 之 中 間 間 隙 流 系 成 穩 定 地 偏 向 其 中 一 柱 。 導 致 上 下 方 雙 柱 相 應 之 平 均 阻 力 係 數 ( C D )不 同 ， 並 隨 著 間 距 之 加 大 兩 柱 之 CD 值 趨 於 一 致 ， 於 間 距 比 大 於 3 後 其 值 約 與 單 一 方 柱 相 同。而 於 雙 柱 縱 向 排 列 時 發 現，在 間 距 比 為 3.5 時， 後 柱 之 CD 值 由 負 值 轉 為 正 值 。. 貳、三維研究 Tutar 等 (2002)在 探 討 雙 棟 並 排 建 築 物 流 場 的 數 值 模 擬 研 究 成 果 指 出，對 於 同 樣 高 度 的 雙 棟 並 排 建 築 物 而 言，隨 著 兩 建 築 物 間 距 的增 大 ， 所 產 生 氣 流 渦 漩 的 強 度 反 而 減 小 。 LES 模 擬 方 法 中 的 RNG 次 網格 尺 寸 模 型 運 用 在 雙 棟 並 排 建 築 物 的 大 氣 流 場 時，大 致 而 言 具 有 較 佳的 模 擬 成 果 且 與 實 驗 結 果 相 符。但 該 模 擬 方 法 在 預 測 該 並 排 建 築 物 間通 道入口端處的流場時，則仍待檢討改進。 卓 勇 志 (2001)探 討 雙 棟 並 排 矩 形 建 築 物 (B/D=2、H/D=7.5)，在 不 同 間 距 比 之 相 鄰 側 風 面 (confronted side face)之 表 面 風 壓 ， 建 築 物 迎 風 面 皆 為 窄 迎 風 面 模 式。研 究 中 發 現，相 鄰 面 靜 態 壓 力 分 布：最大 負 壓 之 CP 均 發 生 在 x/D=0.125、z/D=0.87~0.93 處，分 布 趨 勢 與 單 棟 建 築 側 風 面 相 近。 CP 值 隨 著 間 距 比 縮 小 而 增 加 負 壓。而 史 特 赫 數 隨 間 距 比 增 大 而 減 小，顯 示 兩 棟 建 築 物 越 靠 近，所 產 生 的 渦 漩 逸 散 頻 率越 大 。 此 外 ， 間 隙 風 速 隨 間 距 比 減 小 而 增 大 的 趨 勢 ， 在 S/D < 0.5 時 轉 為相反，間距更小時風速反而減緩。 綜 括 上 述 文 獻 研 究 成 果 指 出，對 於 雙 棟 橫 向 並 排 的 建 築 物 環 境 風 場，無 論 是 二 維 流 況 或 是 三 維 流 況 的 探 討，均 能 發 現 隨 著 建 築 物 之 間. 6.
(19) 第二章 文獻回顧. 的 間 距 加 大，渦 流 溢 散 的 程 度 將 逐 漸 減 弱，甚 至 該 間 距 增 大 至 一 定程 度 時，兩 建 築 物 對 流 場 的 影 響，將 不 再 產 生 互 制 效 應，亦 即 形 成 兩個 獨立的單棟建築物環境風場。. 第二節 耐風設計規範. 參、風速之垂直分佈 風 速 隨 距 地 面 高 度 增 加 而 遞 增，與 地 況 種 類 有 關，依 下 列 指 數 律 公式計算之： Vz ⎛ z ⎞ =⎜ ⎟ V10 ⎝ 10 ⎠. α. ； 0≦ z≦ z g. (2-1). 其中， Vz. ：. 高 度 z 處 之 風 速 (m/sec)。. V 10 ：. 10 公 尺 高 之 風 速 (m/sec)。. α. 相 對 於 10 分 鐘 平 均 風 速 之 垂 直 分 布 法 則 的 指 數，與 地 況 種. ：. 類 有 關 ， 見 表 2-1。 zg. ：. 梯 度 高 度 (m)， 與 地 況 種 類 有 關 ， 見 表 2-1。. 地 況 種 類 依 建 築 物 所 在 位 置 及 其 附 近 地 表 特 性 而 定，分 成 以 下 三 類： (1) 地 況 A： 大 城 市 市 中 心 區 ， 至 少 有 50%之 建 築 物 高 度 大 於 20 公 尺 者。 建 築 物 迎 風 向 之 前 方 至 少 800 公 尺 或 建 築 物 高 度 10 倍 的 範 圍 （ 兩 者 取 大 值 ） 係 屬 此 種 條 件 下 ， 才 可 使 用 地 況 A。 (2) 地 況 B：大 城 市 市 郊、小 市 鎮 或 有 許 多 像 民 舍 高 度（ 10~20 公 尺 ）， 或較民舍為高之障礙物分布其間之地區者。建築物迎風向之前方 至 少 500 公 尺 或 建 築 物 高 度 10 倍 的 範 圍（ 兩 者 取 大 值 ）係 屬 此 種. 7.
(20) 邊界層流中雙棟鄰近建築物之受風效應研究. 條 件 下 ， 方 可 使 用 地 況 B。 (3) 地 況 C： 平 坦 開 闊 之 地 面 或 草 原 或 海 岸 或 湖 岸 地 區 ， 其 零 星 座 落 之 障 礙 物 高 度 小 於 10 公 尺 者 。 若附近地況為介於地況 A 與地況 B 間或地況 B 與地況 C 間之過渡 地 況，原 則 上 應 採 用 會 產 生 較 大 風 力 之 地 況，但 也 可 利 用 可 信 賴 之合 理分析法，決定此一過渡地況之風速垂直分布。. 表2-1 地況相關參數 地況. α. z g (m) b. c. λ(m). ε. z m i n (m). A. 0.32. 500. 0.45. 0.45. 55. 0.5. 18. B. 0.25. 400. 0.62. 0.30. 98. 0.33. 9. C. 0.15. 300. 0.94. 0.20. 152. 0.20. 4.5. 肆、風速壓 各 種 不 同 用 途 係 數 之 建 築 物 在 不 同 地 況 下，離 地 面 z 公 尺 高 之 風 2 速 壓 q ( z )依 下 式 計 算 ， 其 單 位 為 kgf/m 。. q(z) = 0.06 K(z) K zt [I V 10 (C)] 2. (2-2). 式 中 ， K(z) 稱 為 風 速 壓 地 況 係 數 ， 此 值 為 離 地 面 z 公 尺 之 風 速 壓 與 標 準 風 速 壓 （ 地 況 C ， 離 地 面 10 公 尺 處 ） 之 比 值 ， 依 下 式 計 算 ：. 8. 2α. ⎛ z K(z) = 2.774 ⎜ ⎜z ⎝ g. ⎞ ⎟ ； z ＞ 5m ⎟ ⎠. ⎛ 5 = 2.774 ⎜ ⎜z ⎝ g. ⎞ ⎟ ； z ≦ 5m ⎟ ⎠. 2α. (2-3).
(21) 第二章 文獻回顧. 各 種 地 況 種 類 之 α 值 及 梯 度 高 度 z g ， 照 風 力 規 範 2.3 節 規 定 ， 見 表. 2-1 。. 伍、風壓係數與風力係數 計算建築物或地上獨立結構物主要風力抵抗系統之設計風力 時 ， 其 所 使 用 之 風 壓 係 數 Cp （ 封 閉 式 或 部 分 封 閉 式 建 築 物 用 ） 見 表 2-2。. 表2-2 牆之平均外風壓係數 所屬牆面. L/B. Cp. 使用的風速壓. 迎風面. 所有值. 0.8. q(z). 背風面. 0-1 2 ≧4. -0.5 -0.3 -0.2. q(h). 側風面. 所有值. -0.7. q(h). 註 ： L： 平 行 於 風 向 之 結 構 物 水 平 尺 寸 ， m B： 垂 直 於 風 向 之 結 構 物 水 平 尺 寸 ， m G： 陣 風 反 應 因 子 q(h)GC p q(z)GC p. 風向 B. q(h)GC p. q(h)GC p L. 9.
(22) 邊界層流中雙棟鄰近建築物之受風效應研究. 10.
(23) 第三章 研究方法. 第三章 研究方法. 第一節 實驗方法. 壹、試驗設施與量測設備 邊界層流中雙建築物模型試驗於本所之風洞實驗室進行，風洞 本 體 為 一 垂 直 向 的 封 閉 迴 路 系 統 ， 總 長 度 為 77.9m ， 最 大 寬 度 為 9.12m，最 大 高 度 為 15.9m。另 顧 慮 到 未 來 如 進 行 污 染 擴 散 試 驗 或 煙 霧視流試驗可能對風洞本體及工作氣體造成污染，原封閉迴路風洞 切換為開放式風洞。 建築模型以壓克力板粘合製作，以垂直方式固定於第一測試區 截 面 (4m x 2.6m)之 第 二 轉 盤 處 (如 見 圖 3-1)， 而 均 勻 來 流 實 驗 則 固 定 於 第 一 轉 盤 處。模 型 採 用 高 度 (H)120 公 分、厚 度 (D)為 20 公 分 之 方 柱 ， 相 應 之 阻 塞 比 低 於 4%。 模 型 的 四 個 面 皆 均 勻 預 埋 設 壓 力 孔 ， 單 面 壓 力 孔 配 置 64 點，配 置 方 式 如 圖 3-2，編 號 方 式 由 左 下 角 至 右 上角。透過壓力傳輸管線系統與量測儀器連接，再利用電子式壓力 掃 瞄 器 (Scanivalve； ZOC 33)量 測 風 壓 資 料 (如 圖 3-3)。 風 速 剖 面 係 以 熱 線 (Hot Wire) 探 針 配 合 恆 溫 流 速 儀 (Constant Temperature Anemometer; DANTEC 9090N10101) 進 行 量 測 ( 如 圖 3-4、圖 3-5)，而 參 考 風 速 則 採 用 直 式 皮 托 管 (如 圖 3-6)配 合 薄 膜 式 壓 力 計 (VALIDYNE Differential Pressure Transducer； DP103)。 本實驗所使用之量測設備介紹如下：. 11.
(24) 邊界層流中雙棟鄰近建築物之受風效應研究. 圖 3-1 邊 界 層 流 中 雙 建 築 物 模 型. 圖 3-2 模 型 孔 位 配 置 圖. 12.
(25) 第三章 研究方法. (1) 電 子 式 壓 力 掃 瞄 器 表 面 風 壓 量 測 採 用 電 子 式 壓 力 掃 瞄 器 (ZOC 33/64 PX ； 如 圖 3-3) ， 該 系 統 每 個 單 一 模 組 有 64 個 壓 力 輸 入 管 (Pneumatic Inputs)， 對 應 64 個 壓 電 式 壓 力 感 應 器 ， 每 一 個 壓 力 感 應 器 皆 可 單 獨 校 正，輸 入 管 藉 由 內 徑 1mm PVC 管 連 接 至 模 型 量 測 點 以 量 測 壓 力。 而 各 模 組 皆 連 接 至 壓 力 訊 號 處 理 系 統 (RAD BASE 3200)， 此 系 統 可支援類比數位訊號之轉換，最高可支援 8 個模組，其解析度達 16bit， 最 大 採 樣 頻 率 500Hz， 傳 輸 介 面 為 USB， 具 網 路 控 制 及 傳 輸 功能。擷取資料轉換完成之後由此系統傳至個人電腦儲存分析。. 圖 3-3 電 子 式 壓 力 掃 瞄 器 (2) 熱 線 流 速 儀 在 必 要 的 風 場 量 測 上 ， 本 實 驗 使 用 DANTEC 之 熱 線 流 速 儀 (constant temperature anemometer， 如 圖 3-4、 圖 3-5)。 熱 線 探 針的率定配合已完成率定的薄膜式壓力轉換器，在風洞入口處地面 架 上 皮 托 管 (平 行 中 心 線 )， 將 其 動 靜 壓 接 於 薄 膜 式 壓 力 轉 換 器 上 ， 孔 口 位 置 距 離 地 面 40 cm 高 ， 而 熱 線 探 針 利 用 移 動 機 構 架 於 與 皮 托. 13.
(26) 邊界層流中雙棟鄰近建築物之受風效應研究. 管孔口等高、平行、盡量接近的位置，最後利用熱線探針所得的電 壓值與薄膜式壓力轉換器轉換後的風速得出一條四次方多項式的迴 歸曲線。. 圖 3-4 熱 線 探 針. 圖 3-5 熱 線 測 速 儀. 14.
(27) 第三章 研究方法. (3) 皮 托 管 本實驗中採用皮托管進行來流平均風速之量測，由皮托管所量 測 到 的 壓 力 差 值 ， 利 用 伯 努 利 方 程 式 (Bernoulli equation)， 即 依 據後式計算出相應之風速。. U=. 2 Δp. ρ air. (3-1). (4) 壓 力 轉 換 器 本 研 究 採 用 的 壓 力 轉 換 器 為 薄 膜 式 壓 力 轉 換 器 (very-low pressure transducer， VALIDYNE DP103-18， 參 見 圖 3-1)， 具 有 堅 固之金屬外殼，其內部包有一壓電膜片。當受到外部壓力時會導致 金 屬 薄 片 變 形，致 使 產 生 電 壓 變 化，再 經 由 訊 號 放 大 器 讀 出 電 壓 值。 壓 力 轉 換 器 若 與 皮 托 管 (pitot tube)連 接 ， 經 率 定 後 可 用 以 量 測 流 場平均速度。 薄膜式壓力轉換器率定應配合壓力轉換器內部的壓電膜片的受 壓 範 圍，依 照 其 膜 片 可 承 受 範 圍，利 用 壓 力 校 正 器 (DPI 610)連 接 兩 條短油管傳輸壓力給薄膜式轉換器之動壓與靜壓。壓力由小至大， 直 到 可 承 受 之 最 大 壓 力，透 過 資 料 擷 取 系 統 (取 樣 頻 率 為 250Hz，取 樣 時 間 為 70 秒 )將 所 測 之 電 壓 值 轉 換 存 檔 後 ， 其 迴 歸 率 定 曲 線 呈 線 性型態。. 15.
(28) 邊界層流中雙棟鄰近建築物之受風效應研究. 圖 3-6 皮 托 管. 圖 3-7 薄 膜 式 壓 力 計 (5) 資 料 擷 取 系 統. 16.
(29) 第三章 研究方法. 實 驗 所 量 得 之 類 比 訊 號 係 經 由 IOTECH ADC-488/8SA 擷 取 後 作 類 比 數 位 (analog-digital)轉 換 (圖 3-8)。 本 系 統 共 有 8 組 輸 入 端 ， 最 高 採 樣 頻 率 為 100 kHz，具 有 16-bit 之 解 析 度，精 確 度 (accuracy) 高 達 0.02%。 數 位 化 的 訊 號 以 大 於 200 kb/s 的 速 度 經 由 IEEE-488 界面傳至電腦，進行資料儲存與統計運算。. 圖 3-8 資 料 擷 取 系 統. 貳、試驗狀況 試 驗 所 採 用 的 主 要 變 數 有 間 距 比 (S/D)以 及 邊 界 層 特 性 (採 用 均 勻 、 郊 區 、 都 市 )， 建 築 物 高 寬 比 (H/D)為 6， 本 研 究 之 參 數 定 義 與 實 驗 配 置 如 圖 3-9所 示 。 規 劃 測 試 風 速 約 為 8.21m/s、 12.61m/s 以 及 17.06m/s， 雷 諾 數 分 布 約 在 9.1×10 4 至 1.89×10 5 之 間 。 模 型 配 置 之 間 距 比 採 用 0.5、 1、 2、 3、 4、 5(相 應 距 離 為 10、 20、 40、 60、 80、 100 公 分 )， 並 與 單 柱 結 果 作 比 對 分 析 。. 17.
(30) 邊界層流中雙棟鄰近建築物之受風效應研究. 研 究 中 邊 界 層 來 流 欲 採 用 均 勻 (Uniform)、郊 區 與 都 市 等 3 種 特 性 之 地 況 ， 相 應 之 風 速 與 紊 流 強 度 (Turbulence Intensity)剖 面 如 見 圖 3-10與 圖 3-11， 邊 界 層 流 場 特 性 如 表 3-1所 示 。 D. H. U0 S. 圖 3-9 研 究 個 案 示 意 圖. 18.
(31) 第三章 研究方法. Turbulence Intensity (%) 250. 0. 10. 20. 30. 40. 50. Wind Speed. Height (cm). 200. Turbulence Intensity. 150. 100. 50. 0. 0. 5. 10. 15. 20. 25. 30. Wind speed (m/s). 圖 3-10 來 流 邊 界 層 特 性 剖 面 (地 況 B). Turbulence Intensity (%) 250. 0. 10. 20. 30. 40. 50. Wind Speed. Height (cm). 200. Turbulence Intensity. 150. 100. 50. 0. 0. 5. 10. 15. 20. 25. 30. Wind speed (m/s). 圖 3-11 來 流 邊 界 層 特 性 剖 面 (地 況 A). 19.
(32) 邊界層流中雙棟鄰近建築物之受風效應研究. 表 3-1 實 驗 用 流 場 條 件 地況. 地況 B. 地況 A. 邊 界 層 厚 度 δ (cm). 138.82. 139.19. 層 緣 風 速 U δ (m/s). 14.21. 15.35. 冪數率α. 0.21. 0.32. 0.64 6.01 × 10 - 2. 1.10. 參數. *. 剪 力 速 度 u (m/s) 粗 糙 長 度 z o (cm) 地 表 阻 力 係 數 Cd. 2.0 × 10. -3. 0.83 5.02 × 10 - 3. 第二節 量測分析項目. 由單點風壓歷時資料可計算得單點之平均與擾動性風壓係數， 配 合 歷 時 資 料 的 合 成 可 計 算 的 平 均 與 擾 動 性 之 阻 力 /昇 力 係 數。各 係 數定義如下： 平 均 風 壓 係 數 CP =. 擾 動 風 壓 係 數 CP ' =. 平 均 阻 力 係 數 CL =. 擾 動 阻 力 係 數 CD ' =. 20. P ρU 2. (3-2). P '2 2 1 2 ρU. (3-3). 1 2. FX ρU 2 BD. (3-4). FX '2 2 1 2 ρU BD. (3-5). 1 2.
(33) 第三章 研究方法. 平 均 昇 力 係 數 CL =. 擾 動 昇 力 係 數 CL ' =. FY ρU 2 BD. (3-6). FY ' 2 2 1 2 ρU BD. (3-7). 1 2. 第三節 數值方法. 本 研 究 中 流 場 之 模 擬 係 採 用 微 可 壓 縮 流 (Weakly-Compressible-Flow Method[1]，WCF)的 動 力 計 算 方 法。對 無 滯 性 流 (Inviscid Flow)而 言 ， 其 相 應 之 連 續 及 動 量 方 程 式 分 別 為：. ∂ρ + ∇ ⋅ ( ρV ) = 0 ∂t. (3-8). ∂V 1 + V ⋅ ∇ V = − ∇p ∂t ρ. (3-9). 式中 ρ 、 V 、 p 、 t 分別代表流體密度、速度向量、壓力及時間， 而 ν 為 運 動 黏 滯 度 (Kinematic Viscosity)。 在密度僅為壓力的函數的假設條件(即 Barotropic 假定)下，對低馬赫 數流(Low-Mach-Number Flows)而言，空間平均形式紊流(Turbulent Flow) 控制方程式可近似為如下之指標符號形式(" "意指平均(Average)，"’"表攝 動(Fluctuation)部份， ui = ui + ui′ ， p = p + p′ )：. ∂p ∂ + ( ρa 2 u j ) = 0 ∂t ∂x j. (3-10). 21.
(34) 邊界層流中雙棟鄰近建築物之受風效應研究. ∂ ui ∂ ui u j ∂ ( p / ρ ) + = ∂t ∂x j ∂x j ∂ + ∂x j. ⎧⎪ ∂ ui ⎫⎪ ⎨− ui′u ′j − ui′u ′j − ui u ′j − (ui u j − ui u j ) + ν ⎬ ∂x j ⎪⎭ ⎪⎩. (3-11). 由雷諾平均的假設，亦即： ui′u j − ui u′j + (ui u j − ui u j ) = 0. (3-12). (3-11) 式 成 為 ：. ∂ ui ∂ ui u j ∂ p* ∂ ( )+ + =− ∂t ∂x j ∂x j ρ ∂x j. ⎧⎪ ∂ ui ⎫⎪ 1 ⎨− (ui′u ′j − ui′u ′jδ ij ) + ν ⎬ 3 ∂x j ⎪⎭ ⎪⎩. 其 中 δ ij 為 Kronecker delta 函 數 ， p * = p +. ρ 3. (3-13). u i′u ′j 。. 根據渦流粘滯性之假設(Eddy-Viscosity Hypothesis)，次網格紊流剪 應力項可表示為：. 1 − (ui′u′j − ui′u′jδ ij ) = ν t Sij ， 3. Sij = (. ∂ u j ∂ ui + ) ∂xi ∂x j. (3-14). 其 中 ν t 為 紊 流 粘 滯 度 (Turbulence Diffusivity)。而 (3-13) 式 可 重 寫 為：. ∂ ui ∂ ui u j ∂ p* ∂ τ ij + =− ( )+ ( ) ∂t ∂x j ∂x j ρ ∂x j ρ. (3-15). 式 中 τ ij 乃 片 流 粘 性 剪 應 力 以 及 次 網 格 紊 流 剪 應 力 的 和，而 紊 流 擴 散 係 數 則 依 據 Germano 等 [2]建 議 之 動 態 次 網 格 型 式 而 得。(3-10) 式 及 (3-15)式 即 為 本 數 值 模 擬 中 流 場 解 析 之 依 據 。 模 擬 中 針 對 各 計 算 網 格 採 用 MacCormack[3] 之 預 測 步 ⎯ 校 正 步 (Predictor-Corrector)，以顯式(Explicit)的方法處理網格邊界之通量 (Fluxes)進而計算出下一時間之結果。. 22.
(35) 第三章 研究方法. 在數值邊界條件之處理上，對一個固定邊界(Stationary Boundary)而 言，計算中將於鏡射點(Phantom Point)給定適當物理量的方式，以正確地反 應邊界之特性。於速度方面，入流截面之速度為均勻或冪數律(Power Law) 分佈之給定值以定義入流形態；建築物與地表等其他不可穿透邊界(上、下邊 界與出流截面)之鏡射點(如圖 3-12)則分別依據不可滑動(No-Slip)及零梯 度之物理特性給定之(如圖 3-13)。在壓力方面，除將出流截面之平均壓力給 為零值以定義出流場參考壓力外，其餘邊界之壓力給定均基於零梯度之物理 特性。. C D. A. B. UC=UD VC=VD PC=PD. F. Prism. UA=U0 V A= V B PA=PB. E. UE=UF VE=VF ∑PE=0. H G. UG=UH VG=VH PG=PH. 圖 3-12 透 式 邊 界 條 件 說 明 圖. 23.
(36) 邊界層流中雙棟鄰近建築物之受風效應研究. 圖 3-13 不 透 式 邊 界 條 件 說 明 圖. 24.
(37) 第四章 結果討論. 第四章 結果討論. 第一節 模式驗證. 壹、前人實驗驗證 為 確立 本 模式之 正 確性，故 先針對 Murakami[11]的立 方柱 體的實 驗 結果 進 行 比 對 。 本 實 驗 採 用之邊 界 層厚度 (δ /D)為 5.2，冪數 率 為 0.25， 相 應雷諾 數(Reynolds Number)為 7×10 4 。 當邊界 層 流場通 過 方形 柱 體 時 (如 圖 4-1)， 頂 部 銳 緣 處 (Leading edge) 產 生 分 離 區 (Separation) 致 使 往 下 游 處 發 生 再 接 觸 現 象 (Reattachment)，而 在背風 面 處 亦 產 生 一 迴 旋 區 (Recirculation)。在 比對再 接 觸點長 度(Reattachment length)、回 流 區 長 度 (Reverse flow length)以 及迴旋 區 長度(Recirculation Length)後發現 (見 表 4-1)， 數值結 果 與實驗 值 頗為吻 合 ， 其 相 應 誤 差 小 於 11%。 另 比 對 順 流 向 之 柱 體 中 心 線 平 均 表 面 風 壓 發 現 (圖 4-2)， 壓 力 最 大 值 發 生 在 迎 風 面 上 緣 處 (近 B 點 )， 而 在 B 點 附 近 因 強 大 的 逆 壓 力 梯 度(Adverse Pressure Gradient)致 使 柱 體 頂 部 產 生 分 離 區 ， 並 使 壓 力 逐漸回升，在通過背風面之後的迴旋區使得壓力變化趨於穩定。整體 而 言，數 值 結 果 亦 有 不 錯 的 吻 合 度，但 在 柱 頂 銳 緣 (B、C 點)附 近 的 預 測較差，須做進一步修正。. 25.
(38) 邊界層流中雙棟鄰近建築物之受風效應研究. 圖4-1 邊界層流通過建築物之流場示意圖. 表4-1 實驗與數值結果比較表. Murakami[11] 數值結果. 26. 再接觸點長度. 回流區長度. 迴旋區長度. 0.67D 0.75D. 0.56D 0.63D. 1.27D 1.43D.
(39) 第四章 結果討論. 1.5 Experimental, Murakami and Mochida (1988) Numemical. 1. 0.5. Cp. B. C. A. D. Wind. 0. -0.5. -1. -1.5. A. B. C. D. x/D. 圖4-2 柱體中心線表面風壓比較圖. 第二節 雷諾數效應. 本節首先針對均勻來流狀況下，單一柱體柱在不同風速下測試雷 諾 數 效 應 (Reynolds Number Effect)之影 響，採 用 風速 為 8.21m/s、 4 12.61m/s 以及 17.06m/s，相 應 雷 諾 數 (以 Re 表 示 )為 9.1×10 、1.35×. 10 5 、1.89×10 5 ， 其 平 均 風 壓 係 數 ( C p )結 果 如 圖 4-3至 圖 4-5所 示 。 圖 中 a、b、c、d 標 號 分 別 顯 示 迎 風 面、右 側 風 面、左 側 風 面 以 及 背 風 面 。 平均表面風壓係數結果顯示，最大風壓係數發生在迎風面之中心 線，最小值發生在側風面與背風面上緣。迎風面之壓力分布隨著雷諾 數增加，其對稱性越顯著，左右側風面亦呈現同樣之趨勢，且風壓係 數分布變化越趨穩定，可見雷諾數效應已不顯著，故後續之實驗均採 用 雷 諾 數 1.89×10 5 來 進 行 。. 27.
(40) 邊界層流中雙棟鄰近建築物之受風效應研究. -1 -0.8 -0.6 -0.4 -0.2 0. (a). (b). (c). 0.2 0.4 0.6 0.8. 1. (d). 4. 圖4-3 平均風壓係數分布圖(Re=9.1×10 ). -1 -0.8 -0.6 -0.4 -0.2 0. (a). (b). (c). 0.2 0.4 0.6 0.8. 1. (d). 5. 圖4-4 平均風壓係數分布圖(Re=1.35×10 ). 28.
(41) 第四章 結果討論. -1 -0.8 -0.6 -0.4 -0.2 0. (a). (b). (c). 0.2 0.4 0.6 0.8. 1. (d). 5. 圖4-5 平均風壓係數分布圖(Re=1.89×10 ). 0 0.04 0.08 0.12 0.16 0.2 0.24 0.28 0.32 0.36 0.4. (a). (b). (c). (d). 4. 圖4-6 擾動風壓係數分布圖(Re=9.1×10 ). 29.
(42) 邊界層流中雙棟鄰近建築物之受風效應研究. 0 0.04 0.08 0.12 0.16 0.2 0.24 0.28 0.32 0.36 0.4. (a). (b). (c). (d). 5. 圖4-7 擾動風壓係數分布圖(Re=1.35×10 ). (a). 0 0.04 0.08 0.12 0.16 0.2 0.24 0.28 0.32 0.36 0.4. (b). (c). (d). 5. 圖4-8 擾動風壓係數分布圖(Re=1.89×10 ). 30.
(43) 第四章 結果討論. 另 於 數 值 計 算 之 單 柱 平 均 表 面 風 壓 係 數 (圖 4-9)中發 現 ， 最 大 正 壓發生 在 迎風面 中 心線上 緣 處 ， 約 為 1.5。 最 大 負 壓 值 發 生 在 側 風 面 頂 端 前 緣 處(Leading Edge)， 其值約 為-1.2。 在背 風面部 份 ，柱體 下 2/3 處 壓 力 有 回 升 之 狀 況 ， 此 應 是 馬 蹄 型 渦 流 破 壞 順 風 向 之 渦 散 現 象 所造成。 檢視單 柱 擾動表 面 風壓係 數 結 果 (圖 4-10)可 發 現，在 柱 體 中 段 位 置，擾動 風壓值 有 最大值，此 處應 為 渦散(Vortex Shedding)現 象較 完 整的地方，渦流剝離的狀況亦較接近二維特性之位置。而迎風面因不 受 到 渦 散 及 馬 蹄 渦 直 接 影 響 ， 其 擾 動 值 較 小 ， 約 為 0.02。 整體而言，實驗與數值在單棟建築物之平均與擾動風壓係數分布 頗 為 接 近 ， 相 對 位 置 風 壓 係 數 誤 差 值 小 於 10%。. (a). -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1. (b). (c). (d). 圖4-9 平均風壓係數分布圖(數值結果). 31.
(44) 邊界層流中雙棟鄰近建築物之受風效應研究. (a). 0 0.04 0.08 0.12 0.16 0.2 0.24 0.28 0.32 0.36 0.4. (b). (c). (d). 圖4-10 擾動風壓係數分布圖(數值結果). 第三節 流場結果. 檢 視 來 流 邊 界 層 為 地 況 C 時，不 同 間 距 比 (S/D)的 情 況 下，在 高 度 1 分別 為 1/3H(H 為建 築高度 )、 2/3H 的流 場速度 等 值圖 (圖 4-11至 圖. 4-22)。 整 體 而 言，在 建 築 物 後 方 各 自 產 生 獨 立 尾 流 區，僅 在 間 距 比 (S/D) 為 0.5 時 會 有 合 併 之 現 象 。 此外， 當 S/D=0.5 時，雙 柱 中間間 隙 流速 最 高 ， 並 隨 著 間 距 增 大 而 變 小 。 而 當 S/D 大 於 4 後， 間 隙 流 速 反 有 增 大之現象，且流場之整體對稱性亦越佳。 在各個間距比狀況下檢視不同高度流場結果得知，雙棟建築物之 流場特性基本上是對稱的，因此本研究僅針對面向來流左側建築進行 探討。. 1. 取 順 、 橫 風 向 速 度 合 之 純 量 值 ； uv = u 2 + v 2 。. 32.
(45) 第四章 結果討論. 0.2. 0.4. 0.6. 0.8. 1. 1.2. 1.4. 1.6. Y. 0. X. 圖4-11 速度等值圖(S/D=0.5、y/H=1/3). 0.2. 0.4. 0.6. 0.8. 1. 1.2. 1.4. 1.6. Y. 0. X. 圖4-12 速度等值圖(S/D=0.5、y/H=2/3). 33.
(46) 邊界層流中雙棟鄰近建築物之受風效應研究. 0.2. 0.4. 0.6. 0.8. 1. 1.2. 1.4. 1.6. Y. 0. X. 圖4-13 速度等值圖(S/D=1、y/H=1/3). 0.2. 0.4. 0.6. 0.8. 1. 1.2. 1.4. 1.6. Y. 0. X. 圖4-14 速度等值圖(S/D=1、y/H=2/3). 34.
(47) 第四章 結果討論. 0.2. 0.4. 0.6. 0.8. 1. 1.2. 1.4. 1.6. Y. 0. X. 圖4-15 速度等值圖(S/D=2、y/H=1/3). 0.2. 0.4. 0.6. 0.8. 1. 1.2. 1.4. 1.6. Y. 0. X. 圖4-16 速度等值圖(S/D=2、y/H=2/3). 35.
(48) 邊界層流中雙棟鄰近建築物之受風效應研究. 0.2. 0.4. 0.6. 0.8. 1. 1.2. 1.4. 1.6. Y. 0. X. 圖4-17 速度等值圖(S/D=3、y/H=1/3). 0.2. 0.4. 0.6. 0.8. 1. 1.2. 1.4. 1.6. Y. 0. X. 圖4-18 速度等值圖(S/D=3、y/H=2/3). 36.
(49) 第四章 結果討論. 0.2. 0.4. 0.6. 0.8. 1. 1.2. 1.4. 1.6. Y. 0. X. 圖4-19 速度等值圖(S/D=4、y/H=1/3). 0.2. 0.4. 0.6. 0.8. 1. 1.2. 1.4. 1.6. Y. 0. X. 圖4-20 速度等值圖(S/D=4、y/H=2/3). 37.
(50) 邊界層流中雙棟鄰近建築物之受風效應研究. 0.2. 0.4. 0.6. 0.8. 1. 1.2. 1.4. 1.6. Y. 0. X. 圖4-21 速度等值圖(S/D=5、y/H=1/3). 0.2. 0.4. 0.6. 0.8. 1. 1.2. 1.4. 1.6. Y. 0. X. 圖4-22 速度等值圖(S/D=5、y/H=2/3). 38.
(51) 第四章 結果討論. 第四節 風壓係數. 為探討間距比對雙棟建築物之受風效應，本節將針對平均及擾動 風 壓 係 數 在 柱 體 之 迎 風 面、背 風 面、外 側 風 面、內 側 風 面 (靠 鄰 近 柱 面 ) 進行討論。 (1) 平 均 風 壓 係 數 檢 視間 距 比為 0.5 之表面 風 壓結果 (如 圖 4-23)，圖由 左 至右(標 號 a、 b、 c、 d)分 別 代 表 迎 風 面 、 背 風 面 、 外 側 風 面 、 內 側 風 面 。 結 果 發 現 ， 迎 風 面 壓 力 分 布 明 顯 向 右 偏 斜 ， 最 大 壓 力 值 ( C P =0.67)發 生 在 y/D=0.7、z/D=5.0 處。在 背 風 面 壓 力 部 份，其 壓 力 分 布 並 無 明 顯 變 化。而外側風面與內側風面之壓力分布有顯著的不同，在內側風面的 表面壓力，由於收縮效應的影響，其壓力會小於外側風面的壓力值。 且其壓力變化是從上游往下游、頂部往根部逐漸遞增。 當間距 比 增加時 (如 圖 4-24至 圖 4-28)，迎 風 面 風 壓 分 布 越 接 近 左 右對稱形式，而外側與內側風面之壓力分布形式亦越接近。當間距比 大於 4 之後，內側風壓近似外側風壓，表示此時之表面壓力分布受到 鄰近柱體的影響越小。. 39.
(52) 邊界層流中雙棟鄰近建築物之受風效應研究. -1. (a). -0.8 -0.6 -0.4 -0.2. (b). 0. (c). 0.2. 0.4. (d). 0.6. 0.8. 1. 圖4-23 平均風壓係數分布圖(S/D=0.5). -1. (a). -0.8 -0.6 -0.4 -0.2. (b). 0. (c). 0.2. 0.4. 0.6. 0.8. 1. (d). 圖4-24 平均風壓係數分布圖(S/D=1). 40.
(53) 第四章 結果討論. -1. (a). -0.8 -0.6 -0.4 -0.2. (b). 0. (c). 0.2. 0.4. (d). 0.6. 0.8. 1. 圖4-25 平均風壓係數分布圖(S/D=2). -1. (a). -0.8 -0.6 -0.4 -0.2. (b). 0. (c). 0.2. 0.4. (d). 0.6. 0.8. 1. 圖4-26 平均風壓係數分布圖(S/D=3). 41.
(54) 邊界層流中雙棟鄰近建築物之受風效應研究. (a). -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 (b) (c) (d). 1. 圖4-27 平均風壓係數分布圖(S/D=4). -1. (a). -0.8 -0.6 -0.4 -0.2. (b). 0. (c). 0.2. 0.4. 0.6. 0.8. 1. (d). 圖4-28 平均風壓係數分布圖(S/D=5). (2) 擾 動 風 壓 係 數. 42.
(55) 第四章 結果討論. 檢 視 各 間 距 比 之 擾 動 風 壓 係 數 分 布 (圖 4-29至 圖 4-34)後 發 現 ， S/D 為 0.5 時，迎 風面之 C P' 值 約在 0.08 至 0.25 間，但 在 S/D 增加 至 1 時 ， C P' 值 上 升 至 0.28。而 後 隨 著 間 距 比 增 大 ， 表 面 之 C P' 值 逐 漸 下 降 至 穩 定 ， 約 為 0.02。 在 內 側 風 面 部 份，其 不 穩 定 度 均 高 於 外 側 風 面，僅在 S/D 為 5 時， 恢 復 對 稱 之 狀 況。而 在 S/D≦2 時， C P' 最 大 值 發 生 在 柱 體 前 緣 上 方 處； S/D≧3 後 ， C P' 最 大 值 則 落 在 柱 體 中 下 段 。. 0 0.04 0.08 0.12 0.16 0.2 0.24 0.28 0.32 0.36 0.4. (a). (b). (c). (d). 圖4-29 擾動風壓係數分布圖(S/D=0.5). 43.
(56) 邊界層流中雙棟鄰近建築物之受風效應研究. (a). 0 0.04 0.08 0.12 0.16 0.2 0.24 0.28 0.32 0.36 0.4 (b) (c) (d). 圖4-30 擾動風壓係數分布圖(S/D=1). 0 0.04 0.08 0.12 0.16 0.2 0.24 0.28 0.32 0.36 0.4. (a). (b). (c). (d). 圖4-31 擾動風壓係數分布圖(S/D=2). 44.
(57) 第四章 結果討論. 0 0.04 0.08 0.12 0.16 0.2 0.24 0.28 0.32 0.36 0.4. (a). (b). (c). (d). 圖4-32 擾動風壓係數分布圖(S/D=3). 0 0.04 0.08 0.12 0.16 0.2 0.24 0.28 0.32 0.36 0.4. (a). (b). (c). (d). 圖4-33 擾動風壓係數分布圖(S/D=4). 45.
(58) 邊界層流中雙棟鄰近建築物之受風效應研究. (a). 0 0.04 0.08 0.12 0.16 0.2 0.24 0.28 0.32 0.36 0.4 (b) (c) (d). 圖4-34 擾動風壓係數分布圖(S/D=5). 第五節 風力係數 各間距 比 與單柱 之 平均風 力 係數結 果 如表 4-2所示。 當 S/D=0.5 時，平均 阻 力 係 數 值 ( C D ) 為 最 大(1.015)，明 顯高 於單柱 之 實驗結 果。 而 隨 著 間 距 比 增 大 ， C D 值 逐 漸 減 小 並 回 穩 至 單 柱 結 果 (0.828)。 由於鄰近柱體的交互影響，流場整體向內側偏移，因此平均升力 係 數 值 ( C L )在 S/D=0.5 時為 最 小，並 為 負值(-0.226)。且隨 著間距 比 增 大 ， C L 值 逐 漸 趨 於 零 值 (單 柱 結 果 )。. 46.
(59) 第四章 結果討論. 表4-2 平均風力係數表 係數 CD. S/D. CL. 0.5. 1.015. -0.226. 1. 0.865. -0.079. 2. 0.785. 0.008. 3. 0.764. 0.021. 4. 0.757. 0.042. 5. 0.838. 0.026. 單柱. 0.851. -. 在 擾 動 風 力 係 數 部 份 (表 4-3)， 擾 動 阻 力 係 數 值 ( C D' )僅 在 S/D 略 大 (0.092)外，其 餘 間 距 比 之 擾動阻 力 係數極 為 穩定，可 見其受 間 距比 影響並不顯著。 因 渦 流 剝 離 現 象 受 鄰 近 建 築 影 響，故 擾 動 升 力 係 數 ( C L' )在 小間距 比 時 (S/D≦ 4)， 其 C L' 遠 小 於 單 柱 之 結 果 (0.236)； S/D=4 時 ， 其 C L' 已 接 近 單 柱 之結果 (0.23)。. 47.
(60) 邊界層流中雙棟鄰近建築物之受風效應研究. 表4-3 擾動風力係數表 係數. C D'. C L'. 0.5. 0.042. 0.044. 1. 0.092. 0.096. 2. 0.034. 0.040. 3. 0.015. 0.068. 4. 0.029. 0.054. 5. 0.047. 0.105. 單柱. 0.053. 0.161. S/D. 48.
(61) 第五章 結論與建議. 第五章 結論與建議. 第一節 結論. (1) 數 值 結 果 與 實 驗 數 據 有 相 當 不 錯 的 吻 合 度 ， 可 作 為 建 築 物 受 風 壓 初步評估之工具。且在檢視流場與結構體受風力之相關性時，為 一有利之工具。 (2) 平 均 阻 力 係 數 在 間 距 比 0.5 時 有 最 大 值 1.015， 明 顯 高 於 單 柱 之 0.828， 故 在 此 情 況 時 需 慎 重 考 慮 其 風 力 設 計 參 數 。 (3) 當 間 距 比 大 於 4 之 後 ， 建 築 物 各 項 風 力 特 性 已 近 似 單 柱 之 情 形 。 (4) 間 距 比 從 0.5 變 化 至 5 時 ， 迎 風 面 最 大 壓 力 逐 漸 由 由 上 角 轉 移 至 中心線，且流場分布亦逐漸分離成單一柱體之型態。. 第二節 建議. (1) 本 所 於 本 年 度 添 購 一 套 SUN V440 FDS 伺 服 器 主 機 ， 對 於 數 值 計 算 能量有實質之幫助。未來將近一步添購相關運算軟體，以利本模 式使用。 (2) 對 於 氣 動 力 實 驗 所 採 用 之 電 子 式 壓 力 掃 描 系 統 ， 需 建 立 一 套 自 行 校正程序，以確保實驗數據之準確性。 (3) 在 規 劃 氣 動 力 模 型 試 驗 時，銳 緣 附 近 壓 力 孔 位 之 解 析 度 應 提 高 (或 壓 力 變 化 較 大 之 位 置 )， 有 助 於 繪 製 壓 力 分 布 圖 或 尋 找 局 部 極 端 值。. 49.
(62) 邊界層流中雙棟鄰近建築物之受風效應研究. (4) 後 續 應 規 劃 更 多 雙 棟 建 築 物 之 排 列 組 合 狀 況 ， 以 確 立 邊 界 層 流 中 雙棟建築物之風力特性。. 50.
(63) 附錄一 期中簡報會議紀錄與回應. 附錄一 期中簡報會議紀錄與回應. 內 政 部 建 築 研 究 所 95 年 度 研 究 計 畫 「 單 一 建 築 物 樓 頂 風 場 特 性 之 研 究 」、「 低 層 廠 棚 式 建 築 物 風 載 重 評 估 之 研 究 」、「 邊 界 層 流 中 雙 棟 鄰近建築物之受風效應研究」等研究案期中審查會議紀錄 一、日. 期 ： 95 年 7 月 24 日 （ 星 期 一 ） 上 午 9 時 30 分. 二、地. 點 ： 本 大 樓 15 樓 第 一 會 議 室. 三、主 持 人：李主任秘書玉生 四、出席人員：（詳簽到單） 五、記. 錄：張恭銘、曹源暉、黎益肇. 審查意見回應對照表 審查委員. 審查意見. 回應. 王技師亭復. 1. 在 兩 棟 建 築 物 相 鄰 內 側. 1. 在 期 末 報 告 說 明. 及 各 自 外 側 端 區 (End Zone) 之 ； 圖 比 例 均 相 同 。 或 端 邊 帶 (Edge Strips) 面. 2. 待 案 例 分 析 完 整. 寬 之 10%兩 端 區 之 u/Uo 對. 後，將 與 現 行 之 風 力. z/ Ｄ 與 x/ Ｄ 高 度 變 化 如. 規範做整合分析。. 何 ？ 其 比 例 是 否 與 圖 11 至. 3.遵 照 辦 理 。. 圖 13 相 同 。 2. 建 議 將 研 究 成 果 整 合 納 入 風 力 規 範 中，對 於 業 界 有 實質之幫助。 3. 兩 棟 建 築 物 相 鄰 內 側 面 之 速 度 亦 請 比 照 圖 11 做 比 較。. 51.
(64) 邊界層流中雙棟鄰近建築物之受風效應研究. 方教授富民. 1. 本 研 究 中 數 值 計 算 是 主. 1.已 於 11 月 底 添 購. 要 工 作 之 一，然 以 建 研 所 風. 大型運算主機。. 洞 實 驗 室 而 言，尚 無 計 算 流. 2.遵 照 辦 理 。. 體 動 力 學 (CFD)試 驗 室 ， 以 提 供 所 需 之 計 算 能 量。建 議 建 研 所 積 極 成 立 CFD 試 驗 室，以 達 國 際 風 洞 試 驗 室 之 水準。 2.在 檢 討 雙 棟 間 距 比 (S/D) 之影響時，應與單柱結果 (S/D→ 無 窮 大 )比 較 ， 建 議 增加單柱相應結果之呈現。 朱教授佳仁. 1.報 告 中 所 指 之 「 人 行 風 」 1. 行 人 風 高 度 依 縮 一 般 稱 為 「 行 人 風 」 (Pedestrian. 尺 不 同 而 變 化，於 此. level 列 入 整 個 風 速 剖 面 。. wind)， 報 告 中 未 說 明 其 高. 使 用 LES Model 並 無. 度為何？計算模式是否有. Wall Function 問 題。. 使 用 「 邊 牆 函 數 」 (Wall 2.感 謝 建 議，並 納 入 Function) ？ 計 算 結 果 是 否 會 受 使 用 (或 不 使 用 )邊 牆 函數的影響？ 2.建 議 間 距 比 S/D 可 以 再 加 一. 組. S/D<0.5( 譬. 如. S/D=0.25) 與 S/D=0.5 比 較，其 結 果 可 以 應 用 於 集 合 式 住 宅 大 樓，或 防 火 巷 寬 度 (間 距 )的 決 定 。. 52. 期末報告方向。.
(65) 附錄一 期中簡報會議紀錄與回應. 蕭教授葆羲. 1. 期 中 進 度 皆 已 達 到 預 期. 1.感 謝 肯 定 。. 之工作進度。. 2. 於 期 末 報 告 中 納. 2. 未 來 應 將 擾 動 風 壓 均 方. 入。. 根值之分布納入探討。. 3. 在 雙 棟 建 築 物 對. 3. 關 於 平 均 風 壓 係 數 分. 稱 排 列 狀 況 下，其 平. 布，於 迎 風 面、背 風 面、內. 均流場為對稱之狀. 側、外 側 均 僅 繪 出 一 棟，是. 況，其 說 明 會 在 期 末. 否意味該平均風壓係數之. 報告中加入。. 分佈變化為對稱？請於報 告之結果中註明。 中華民國土木. 1. 本 研 究 案 既 在 探 討 對 建. 1.感 謝 指 正，並 依 審. 技師公會全聯. 築 物 之 受 風 效 應，建 議 應 與. 查委員意見修改。. 會. 目前耐風設計規範內使用. 2.感 謝 建 議，並 納 入. 鍾技師肇滿. 之專有名詞及內容相呼應。 期末報告方向。 第 1 頁 有 部 分 錯 別 字，請 注 意修正。 2. 建 議 本 案 研 究 結 果 可 在 實務上對結構體設計之影 響，提 出 應 予 注 意 之 處，以 利實務上之應用。. 53.
(66) 邊界層流中雙棟鄰近建築物之受風效應研究. 中華民國建築. 1. 本 案 採 取 雙 棟 方 柱 平 行. 1.個 案 頗 多，會 在 下. 師公會全聯會. 横 列 之 配 置，且 僅 改 變 高 度. 年度計劃或延伸研. 溫建築師卓炫. 與 間 距 值 進 行 測 試 (主 要 變. 究時列入考量。. 數 為 高 寬 比 H/D 及 間 距 比. 2.遵 照 辦 理 。. S/D)，針 對 兩 棟 之 B/D(深 寬. 3.遵 照 辦 理 。. 比 )或 不 平 行 橫 列 (Side By 4.感 謝 指 正 。 Side)時 ， 則 未 測 試 。 請 往. 5. 模 擬 個 案 完 整 後. 後能繼續補充上述尺寸規. 進行。. 模變化時之探討。 2.第 4 頁 之 平 面 請 以 剖 面 圖 表 示，另 行 人 風 示 意 圖 部 分 也請加入報告中。 實驗用流場條件僅列「市 郊 」及「 都 市 」地 況，請 依 「 耐 風 設 計 規 範 」三 種 地 況 A、 B、 C 三 類 進 行 實 驗 。 3. 請 考 慮 不 同 風 向 的 影 響，使 用 風 洞 試 驗 測 點 之 設 計 風 力 及 風 力 效 應 時，應 請 遵 照 耐 風 設 計 規 範 5.3 節 之 模 擬 準 則，並 考 慮 不 同 方 向 之 影 響，及 其 合 理 組 合 中 最 不利的風力載重。 4. 請 補 充 符 號 說 明 及 名 詞 定義，並與「耐風設計規 範 」契 合，同 時 編 列 目 錄 及 頁碼。. 54.
(67) 附錄一 期中簡報會議紀錄與回應. 5.試 驗 數 值 請 與「 耐 風 設 計 規 範 」相 關 列 表 中 之 值 進 行 比較與評估。. 55.
(68) 邊界層流中雙棟鄰近建築物之受風效應研究. 附錄二 電子式壓力掃瞄器操作說明 壹、硬體安裝 P×A→ P×A cal→ cal C→ 校 正 (校 正 時 接 壓 力 校 正 器 ， 量 測 時 接 參 考 壓 力 ) R→ 背 景 值 (參 考 壓 力 ， 接 皮 托 管 靜 壓 值 或 側 壁 靜 壓 ). 貳、開機程序 1.Power supply 2.RPM 1000 3.scanivalve 主 機. 參、指令說明 首 先 將 RAD 程 式 (C:\PROGRAM FILE\SCANIVALVE\RAD.EXE)打 開 ， 接 著 輸 入 所 需 之 指 令 ,各 指 令 與 其 所 含 意 義 說 明 如 下 ：. 一、設定所需要之測點數： 1. set chan1 0 (將 channel 1 的 測 點 設 定 先 歸 零 ) 2. set chan1 1-1..1-27. (假 設 第 一 個 module 使 用 1 至 27 點 ). 3. set chan1 1-30..1-64. (假 設 第 一 個 module 使 用 30 至 64 點 ). 4. set chan1 2-1..2-64. (假 設 第 二 個 module 使 用 1 至 64 點 ). 二、設定擷取資料值與輸出型式 1. set fps1 2000. (設 定 擷 取 資 料 筆 數 2000). 2. set period 500 (設 定 掃 描 率 500， 也 就 是 取 樣 頻 率 500HZ). 56.
(69) 附錄一 期中簡報會議紀錄與回應. 3. set format 1. (設 定 位 置 輸 出 型 式 ， 共 有 0、 1、 2 三 種 型 式 ). 4. save. 三、設定偵測數據顯示於螢幕上或自動存檔至資料夾 1. set conout 2. (設 定 偵 測 數 據 顯 示 於 螢 幕 上 但 不 存 檔 ). 2. set conout 3. (設 定 偵 測 數 據 自 動 存 檔 至 資 料 夾 ). 四、再確認設定群組與各項設定值 1. list sg 1 2. list I. (確 認 設 定 群 組 ,主 要 確 認 fps1 與 chan1). (確 認 設 定 值 ,主 要 確 認 format 與 conout). 五、設定為校正模式 1. dout 2 1. (充 氮 氣 ). 2. dout 1 0 3. dout 3 0 4. calz (利 用 氮 氣 清 除 管 內 壓 力 及 殘 留 氣 體 ) 5. scan (執 行 偵 測 ) 6. close. (在 所 設 定 之 偵 測 時 間 過 後 輸 入 close,關 閉 偵 測 指 令 ). 附 註 ： 一 般 氣 體 壓 力 控 制 應 小 於 65 pa. 六、設定為正式操作模式 1. dout 1 1 2. dout 2 0 3. dout 3 0 4. scan (執 行 偵 測 ) 5. close. (在 所 設 定 之 偵 測 時 間 過 後 輸 入 close,關 閉 偵 測 指 令 ). 57.
(70) 邊界層流中雙棟鄰近建築物之受風效應研究. 七、實驗程序 (1)設 定 與 校 正 1. set fps1 取 樣 數 目 2. set period 取 樣 頻 率 3. set chan1 0 4. set chan1 1-1..1-64. (自 行 調 整 ， 未 用 孔 位 請 封 住 ). 5. dout 2 1 6. dout 1 0 7. calz. 二、調整為偵測模式 1. dout 1 1 2. dout 2 0. 三、資料擷取 1. set conout 3(螢 幕 輸 出 時 採 用 set conout 2) 2. scan 3. close. 58. (此 指 令 下 達 後 ， 方 完 成 存 檔 工 作 ).
(71) 參考書目. 參考書目 1. Song, C. and Yuan, M. (1988), "A Weakly Compressible Flow Model and Rapid Convergence Method," ASME Journal of Fluids Engineering, Vol. 110, No. 4, pp.441-455.. 17 2. Germano, U., Piomelli, P. and Cabot, W.H. (1991), “A Dynamic Subgrid-scale Eddy Viscosity Model,” Physics of Fluids, A(3), 1760-1765. 3. MacCormack,. R.. (1969),. "The. Effect. of. Viscosity. in. Hyper-velocity Impact Cratering," AIAA paper No. 69-354, American Institute of Aeronautics and Astronautics, Washington D.C. 4. Sakamoto, H., Haniu, H. and Obata, Y. (1987), “Fluctuating Force Acting on Two Square Prisms in a Tandem Arrangement,” Journal of Wind Engineering and Industrial Aerodynamic, Vol. 26, pp. 85-113. 5. Tutar, M. and Oguz, G. (2002), “Large Eddy Simulation of Wind Flow Around Parallel Buildings with Varying Configurations,” Fluid Dynamics Research, Vol.31, pp. 289-315. 6. Luo, S.C., Li, L.L. and Shah, D.A. (1999), “Aerodynamic Stability of the Downstream of Two Tandem Square-section Cylinders,” Journal of Wind Engineering and Industrial Aerodynamic., Vol. 79, pp. 79-103. 7. Sun, T. F. and Gu, Z. F. (1995), “Interference Between Wind Loading on Group of Structure,” Journal of Wind Engineering and Industrial Aerodynamic., Vol. 54-55, pp. 213-225. 8. Liu, C. H. and Chen, J. M. (2002), “Observations of hysteresis in flow around two square cylinders in a tandem arrangement,”. 59.
數據
+2
Outline
相關文件
Randall,D.,te al,(1997),Water Supply Planning Simulation Model Using Mixed-Integer Linear Programming, Jouenal of Water Resources Planning and
and Shinmoto, Y.,” Effects of Dynamic Stall on Propulsive Efficiency and Thrust of Flapping Airfoil “, AIAA JOURNAL Vol. Liou, “Numerical Simulation of Dynamics Stall Using Upwind
Therefore, the purpose of this study is to perform a numerical analysis on the thermal effect of shape-stabilized PCM plates as inner linings on the indoor air temperature
Y., (1999), “Examining the technology acceptance model using physician acceptance of telemedicine technology,” Journal of Management Information Systems, Vol. and Baroudi,
(1985), “A Conceptual Model of Service Quality and Its Implications for Future Research,” Journal of Marketing, Vol. (1988), “SERVQUAL: A Multiple-Item Scale for Measuring
In Type I, one performed comparison of 2D&3D model, impact of intrinsic stress of devices, package thermal effect, and transistor location effect.. On the other hand,
Shinar (2002), “Effect of Uncertainty, Transmission Type, Driver Age and Gender on Brake Reaction and Movement Time,” Journal of Safety Research,
Harmelink (2001), ”Comparison of Linear Scheduling Model (LSM) And Critical Path Method (CPM).” Journal of Construction Engineering and Management, Vol. (1993), “Project
