坡地建物耐震設計對於位處地震頻繁的我國而言，相當重要，但有關的文獻資料卻付之闕如，故本研究針對此一課題進行探討。
    本研究使用PLAXIS有限元素軟體，並參考台北地層鑽探資料，建置邊坡模型，輸入921地震歷時，並以建築耐震設計規範所訂的30年、475年、2500年再現週期所對應的最大地表加速度0.12g、0.29g、0.51g進行邊坡受震行為分析。其模擬分為三部分：(1)自然邊坡受震行為模擬 (2)建物與邊坡互制關係 (3)人工擋土邊坡受震行為觀察。望藉由模擬，了解邊坡受震所產生的反應及趨勢，進而提供邊坡設計之參考使用。
    研究結果顯示：(1)邊坡受震而尚達破壞前，最大位移發生在坡面，坡頂次之。邊坡破壞後，坡趾位移量最大。(2)邊坡受震時，若坡高不變，改變坡角，將直接影響坡頂受震反應，使坡頂位移增量大於坡面位移增量。(3)若坡角固定，改變坡高時，將放大邊坡受震反應，隨著坡高增加，邊坡位移量及位移增量隨之增加。(4)在邊坡頂緣及趾緣建置建物進行測試可知，在坡趾地面上加置建物，受震時對邊坡及建物的影響不大；但坡頂地面上之建物退縮距離小於10m時，明顯影響邊坡及建物受震行為。(5)建物的加重、加寬或降低樓層高度，均有助於減緩建物於受震時的反應，降低建物位移量的產生。(6)將邊坡進行挖除設置人工擋土設施時，將導致坡頂及擋土設施頂端位移量明顯增加；坡面位移則反而較原邊坡為低。(7)在人工邊坡頂緣後方加設建物，受震時，擋土設施牆頂為最大位移發生處，在人工邊坡頂緣加置建物，將增加邊坡與擋土設施位移及危險性。(8)擋土設施高度增加時，將放大其受震反應，可知擋土設施高度的增加，首先促使擋土牆頂產生受震反應放大現象。(9)當固定挖方，改變挖方位置時，擋土設施及開挖位置距坡頂越遠，將有助於受震時擋土設施及坡頂的穩定。

Slope stability problems have been studied considerably in the past years. Relevant design codes and specifications were established based on the learning. In engineering practice, more people use FDM or FEM analysis while the physical modeling becomes more attractive for particular research studies. In this paper, 2D FEM analysis was conducted using computer program PLAXIS. Numerical examples of a homogeneous earth slope (10m height and 30o inclination) were studied extensively for the design earthquake acceleration (0.29g) under a return period of 475 years in Taipei. Considering the geographic and geological conditions of the local seismic stations and with the past significant earthquakes into the measurement, the acceleration time histories recorded at a particular seismic station located at Chungho District in New Taipei City were calibrated to target PGA for the inputs. Varying the PGA, the soil parameters and the slope geometry, the maximum horizontal displacements of the slope were computed and compared to make the observations and discussions. Slopes with and without the building on the top or at the toe were studied to see their interactions. In general, it is found that local Design Specification on minimum distance of the building from a slope is appropriate for the seismic concerns. Potential defects will be discussed with more elaborations.
    The observation are summarized as follows：(1) Before the slope collapsed due the earthquake, the maximum displacement will occur on the slope surface where the crest displacement is also found critical. (2) As the slope angle increases, the maximum displacement of the slope will move to the crest. (3) By changing the height of the slope with the same inclination angle, the displacement will increase. (4) A building around the toe will not affect the response very much, however a building on the crest will affect significantly the response especially within a certain distance, say 10m. (5) By increasing the weight and width of the building or reducing the height will reduce the building responses. (6) An excavation on the slope will make some amplification for the slope on top of the excavation, however the rest parts of the slope are de-amplified. (7) If a building is built behind the man-made (say retaining wall exist) slope, the construction will enlarge the slope displacements. (8) The height of excavation or the retaining wall is important, its increase will make the slope responses worse especially for slope on top of the excavation. (9) With the same amount of soils in the excavations, the lower the excavation is, the smaller the response will be; which indicates that when making the land developments on the sloping ground, the higher the location is, the more attentions must be paid.

目錄
第一章 緒論	1
1-1 研究動機與目的	1
1-2 研究方法與內容	2
第二章 文獻回顧	4
2-1 前言	4
2-2 靜態分析方法	9
2-2-1擬靜力分析 (Pseudo static Analysis)	9
2-2-2 Newmark分析理論	14
2-3數值分析方法	18
2-4坡地受震試驗	21
2-4-1離心機試驗	24
2-4-2振動台試驗	26
2-5國內外坡地相關規範	32
2-5-1國內坡地規範	32
2-5-2國外坡地規範	37
2-6台灣地區地震危害度	41
第三章  PLAXIS介紹及初步模擬結果	45
3-1 PLAXIS 簡介	45
3.2 PLAXIS程序組成及指令說明	49
3-2-1 PLAXIS基本操作	49
3-2-2 PLAXIS材料組構模式	57
3-2-3莫爾庫倫土壤模型	59
3-2-4 PLAXIS元素類型	65
3-2-5 PLAXIS邊界條件	67
3-2-6 PLAXIS動態模擬	70
3-2-7 PLAXIS與EQWEAP比較結果	72
3-3有限元素模型分析探討	73
3-3-1地震力選取	73
3-3-2基線修正法	80
3-3-4小結	93
第四章  自然邊坡與建物	95
4-1自然邊坡模型設置	95
4-1-1土壤凝聚力影響	100
4-1-2土壤摩擦角影響	105
4-1-3土壤剪力波速影響	107
4-1-4坡角影響	109
4-1-5坡高影響	112
4-1-6地表加速度對邊坡之影響	113
4-2 建物與邊坡	115
4-2-1坡地建物模型	115
4-2-2上坡建物對邊坡之影響	117
4-2-3下坡建物對邊坡之影響	119
4-2-4邊坡建物與自然邊坡比較	122
4-2-5建物高對邊坡之影響	125
4-2-6 建物重對邊坡之影響	127
4-2-7 建物寬對邊坡之影響	128
第五章 人工邊坡	130
5-1 前言	130
5-2人工邊坡	130
5-2-1 擋土牆設置模型	130
5-2-2 人工邊坡與自然邊坡比較	132
5-2-3 擋土牆高度對邊坡之影響	133
5-2-4 挖方位置對人工建物的影響	135
5-2-5 地表加速度對人工邊坡的影響	136
5-2-6 上坡建物對人工邊坡的影響	138
5-2-7  加設地錨對人工邊坡的影響	140
第六章 結論與建議	142
6-1 結論	143
6-2 建議	146
6-3 未來展望	147
參考文獻	148

表目錄
表2-1  靜態分析方法模擬邊坡受震之相關文獻一覽	36
表2-2  數值分析方法模擬邊坡受震之相關文獻一覽	6
表2-3  坡地受震試驗模擬邊坡受震之相關文獻一覽	7
表2-4  坡地退縮距離相關規範及文獻一覽	8
表2-5  自然河岸高度超過五公尺範圍者(摘自建築技術規則，2002)	36
表2-6  各種土壤內摩擦角(陳世芳，1979)	40
表2-7  台灣主要都市再現週期對應之 IM 地表加速度值	43
表3-1  PLAXIS (2011)六節點形狀函數	52
表3-2  PLAXIS (2011)十五節點形狀函數	53
表3-3  PLAXIS所提供之土壤材料模式	58
表3-4  台北地區921地震各測站記錄(整理自中央氣象局網站)	76
表3-5  台北地區331地震各測站記錄(整理自中央氣象局網站)	77
表3-6  PLAXIS元素密度	92
表3-7  土壤模擬參數	94
表4-1  改變凝聚力於各點位移增加量表	101
表4-2  改變摩擦角於各點位移增加量表	106
表4-3  改變剪力波速於各點位移增加量表	108
表4-4  改變坡角於各點位移增加量表	110
表4-5  改變坡角於各點位移增加量表	111
表4-6  改變坡高於各點位移增加量表	113
表4-7  改變PGA於各點位移增加量表	114
表4-8  建築物剛度性質表	116
表4-9  改變建物於頂緣位置對各點位移影響增加量表	118
表4-10 改變建物於址緣位置對各點位移影響增加量表	120
表4-11 改變建物高度對各點位移影響增加量	126
表4-12 改變建物重量對各點位移影響增加量	128
表4-13 改變建物寬對各點位移影響增加量	129
表5-1  人工擋土設施之材料性質	131
表5-2  改變擋土高度對人工邊坡之位移增加量表	134
表5-3  改變挖方位置對人工邊坡之位移增加量表	136
表5-4  改變地表加速度對人工邊坡之位移增加量表	137

圖目錄
圖1-1    研究流程圖	3
圖2-1    擬靜力分析滑動情況下內力作用圖(摘自羅曉萍，2006)	9
圖2-2    潛在滑動坡體之擬靜力分析力素(摘自Seed，1979)	12
圖2-3    Newmark位移法二次積分之示意圖	16
圖2-4    平面型地滑示意圖(重繪自張石角，1997)	33
圖2-5    擋土建物與安全距離示意圖(重繪自建築技術規則，2002)	34
圖2-6    日本陡坡地崩坍危險區域範圍基準示意圖(重繪自日本急傾斜地崩壞防止工事技術指針，1998)	38
圖2-7    邊坡崩坍最大波及距離示意圖(重繪自林澄清，2004)	39
圖2-8    不同地震危害度曲線比較	44
圖2-9    台灣主要都是地震危害度曲線	44
圖3-1    平面應變與軸對稱示意圖(摘自PLAXIS參考手冊)	49
圖3-2    土壤元素節點與應力點位置(重繪自 PLAXIS 2D Reference，2011)	51
圖 3 3   PLAXIS 6節點三角形元素示意圖	52
圖 3 4   PLAXIS 15節點三角形元素示意圖	52
圖3-5    Mohr-Coulomb 完全塑性之應力－應變關係(摘自PLAXIS 2D Reference，2011)	62
圖3-6    有效主應力空間中之不規則六角錐降伏函數	65
圖3-7    PLAXIS阻尼邊界示意圖	68
圖3-8    EQWEAP與PLAXIS模擬水平地表受震之位移反應圖	73
圖3-9    積穗國小位置圖	74
圖3-10   積穗國小921地震歷時	75
圖3-11   積穗國小331地震歷時	75
圖3-12   基線修正前之位移歷時圖	82
圖3-13   基線修正後之位移歷時圖	82
圖3-14   水平測試標準模型	83
圖3-15   PLAXIS使用動態分析示意圖	87
圖3-16   PLAXIS設定附加計算步示意圖	88
圖3-17   PLAXIS設定動態步數示意圖	88
圖3-18   PLAXIS輸入地震資料示意圖	89
圖3-19   地層厚度對水平地表位移影響之位移圖	90
圖3-20   寬幅對水平地表位移影響之位移圖	91
圖3-21   網格疏密度對水平地表位移影響之位移圖	93
圖4-1    標準邊坡模型	96
圖4-2    自然邊坡各記錄點示意圖	96
圖4-3    坡頂30m處速度歷時	97
圖4-4    坡頂20m處速度歷時	97
圖4-5    坡頂10m處速度歷時	97
圖4-6    坡頂速度歷時	97
圖4-7    坡面速度歷時	97
圖4-8    坡趾速度歷時	97
圖4-9    坡趾10m處速度歷時	97
圖4-10   坡趾20m處速度歷時	97
圖4-11   坡趾30m處速度歷時	97
圖4-12   坡頂30m處速度歷時	97
圖4-13   坡頂20m處速度歷時	97
圖4-14   坡頂30m修正結果(最大位移量23.5cm)	98
圖4-15   坡頂修正結果(最大位移量29.9cm)	….98
圖4-16   坡面修正結果(最大位移量33.5cm)	….98
圖4-17   坡趾修正結果(最大位移量25.3cm)	….98
圖4-18   坡趾10m修正結果(最大位移量24.8cm)	.98
圖4-19   坡趾20m修正結果(最大位移量25cm)	.98
圖4-20   坡趾30m修正結果(最大位移量25cm)	.98
圖4-21   邊坡受震時921及331地表位移比較圖	99
圖4-22   邊坡受震時凝聚力對地表位移之影響圖	100
圖4-23   C=0.01 kN/m2時土壤強度折減分析圖	101
圖4-24   C=1 kN/m2時土壤強度折減分析圖	102
圖4-25   C=2 kN/m2時土壤強度折減分析圖	102
圖4-26   C=4 kN/m2時土壤強度折減分析圖	103
圖4-27   C=8時 kN/m2土壤強度折減分析圖	103
圖4-28   C=10時 kN/m2時土壤強度折減分析圖	104
圖4-29   C=20時 kN/m2時土壤強度折減分析圖	104
圖4-30   邊坡受震時摩擦角對地表位移之影響圖	105
圖4-31   邊坡受震時剪力波速對地表位移之影響圖	107
圖4-32   邊坡受921地震時坡角對地表位移之影響圖	109
圖4-33   邊坡受331地震時坡角對地表位移之影響圖	110
圖4-34   邊坡受震時坡高對地表位移之影響圖	112
圖4-35   邊坡受震時坡高對地表位移之影響圖	114
圖4-36   邊坡頂緣設置建物標準模型示意圖	116
圖4-37   邊坡受震時頂緣建物後退距離對地表位移之影響圖	118
圖4-38   邊坡受震時址緣建物後退距離對地表位移之影響圖	120
圖4-39   邊坡受震時建物後退對地表位移之影響圖	121
圖4-40   自然邊坡與坡頂加設建物比較圖	122
圖4-41   自然邊坡土壤強度折減分析圖	123
圖4-42   建物距頂緣5m之土壤強度折減分析圖	123
圖4-43   建物距頂緣10m之土壤強度折減分析圖	124
圖4-44   建物距頂緣20m之土壤強度折減分析圖	124
圖4-45   邊坡受震時建物高度對地表位移之影響圖	126
圖4-46   邊坡受震時建物重量對地表位移之影響圖	127
圖4-47   邊坡受震時建物寬對地表位移之影響圖	129
圖5-1    人工擋土設施標準模型示意圖	131
圖5-2    邊坡受震時人工邊坡與自然邊坡比較圖	132
圖5-3    邊坡受震時改變擋土高度對人工邊坡之影響圖	134
圖5-4    邊坡受震時改變剪力波速對人工邊坡之影響圖	135
圖5-5    邊坡受震時改變地表加速度對人工邊坡之影響圖	137
圖5-6    人工邊坡受震時加設上坡建物前後比較圖	138
圖5-7    人工擋土設施加設地錨之標準模型示意圖	141
圖5-8    加設地錨後人工邊坡受震時位移圖	142
圖5-9    人工邊坡受震時加設地錨前後比較圖	142

1.Anteneh Biru Tsegaye (2010), PLAXIS Liquefaction model.
2.Arias, A. (1970), “A Measure of Earthquake Intensity”, Seismic Design for Nuclear Power Plants, Massachusetts Institute of Techology Press, Cambridge, MA, pp.438-483
3.Aurelian C. Trandafir (2008) “ Earthquake-induced displacements of gravity retaining walls and anchor-reinforced slopes” soil dynamics and earthequake eng
4.Chang , M. H., Mitchell, J.K. and Seed, R. B.(1999) “Model Test on the 1988 Kettleman Hill Landfill Slope Failure.”Geotechnial Testing Journal, 22,61-66
5.Chen, H. T., Hung, W. Y., Chang, C. C., Chen, Y. J. and Lee, C. J. (2007) “Centrifuge Modeling Test of a Geotextile-Reinforced Wall with a Very Wet Clayey Backfill.”Geotextiles and Geomembranes, 25, 346-259.
6.Clough R W. Pirtz D. (1956)” Earthquake Resisitance of Rockfill Dams” jounal for the Soil Mechanics and Foundations Division”, ASCE. Vol.89, No.1
7.George D. Bouckovalas (2005) ” Numerical evaluation of slope topography effects on seismic ground motion”, Soil Dynamics and Earthquake Engineering 25,pp. 547–558
8.Hong, Y. S., Chen, R. H., Wu, C. S. and Chen, J. R. (2005) “Shaking Table Tests and Stability Analysis of Steep Nailed Slopes.”Canadian Geotechnical Journal, 42,1264-1279
9.Hynes-Griffin, M.E. and Franklin, A.G. (1984) “Rationalizing the seismic coefficient method” ,U.S. Army corps of Engineering Waterways experiment Station, Misscellaneous Paper GL-84-13
10.Iai, S. (1989),“Similitude for Shaking Table Tests on Soil-Structure
-Fluid Model in 1-g Gravitational Field.” Soils and Foundations, 29, 105-118.
11.Jibson, R.W. (1993 ) “Predicting Earthquake-Induced Landslide Displacements Using Newmark’s Sliding Block Analysis” , Transportation Research Record, 1411,pp.9-17
12.Jibson, R.W.,Harp, E.L. and Michael, J.A. (2000), “Method for Producing Digital Probabilistic Seismic Landslide Hazard Maps”, Enginneering Geology, v25, n3-4, pp.271-289
13.Kagawa, T. (1978) “On the Similitude in Model Vibration Tests of Earth Structures.”Proceedings of Japan Society of Civil Engineering, 275, 69-77.
14.Keefer, D.K. and R.C. Wilson, (1989), “Predicting Earthquake-Induced Landslides with Emphasis on Arid and Semi-arid Environment”,In Landslides in a Semi-Arid Environment, Inland Geological Survey Society, Riverside, California , Vol. 2, pp. 118-149
15.Keefer, D.K. (2000)“Statistical Analysis of an Earthquak-Induced Landslides Distribution 1989 Loma Prieta ,California Event”,Geol Soc.  Am. Bull., pp.231-249 
16.Khazai, B., and Sitar, N., (2000). “Assessment of Seismic Slope Stability Using GIS Modeling”, Geographic Information Sciences, Vol. 6, No.2, 121-128. 
17.Kutter, B.L.(1982), Centrifugal Modeling of the Response of Clay Embankments to Earthquakes, Ph.D Thesis, Cambridge University.
18.Arai K. (2009) “Model Test and Numerical Analysis of Slope Stability during Earthquake” Proceedings of the International Conference on Performance-Based Design in Earthquake Geotechnical Engineering ,pp.873-880.
19.Lin, M. L. and Wang, K. L. (2006)“Seismic Slope Behavior in a Large Shaking Table Model Test.”Engineering Geology, 86, 118-133.
20.Lin, M.L. and Wang, K.L. (2005), “Seismic Slope Behavior In A Large-Scale Shaking Table Model Test.”, Engineering Geology , Volume 86, Issues 2–3, pp. 118–133. 
21.Low, B.K. and Tang, W.H. (2004) “Reliability Assessment of Slope Reliability for Effective Landslide Hazard Management," Geotechnique, February,Vol. LV, No. 1, pp85-94.
22.Newmark, N.W. (1965) “Effects of Earthquake on Dams and Embankments”,Geotechnique, Vol.15, No.2 pp. 139-160.
23.Okamoto , S. (1980) Introduction to Earthquake Engineering, University of Tokyo press, Tokyo 629p
24.PLAXIS (2011), PLAXIS 2D Material Models Manual 2010
25.Proceedings of the 1988 Conference on Geotechnaical Earthquake Engineering and Soil Dynamics III, pp. 779-789.
26.Roark, Raymond J.(1954) “Formulas for Stress and Strain, 3rd Edition”, McGraw-Hill, New York
27.Roscoe, K. (1968)“Soils and Model Tests.”Journal of Strain Analysis,3,57-64
28.Rocha, M. (1957) “The possibility of Solving Soil Mechanics Problems by the Use of Models.”4 th International Conference on Soil Mechanics and Foundation Engineering, London, U.K.,1,183-188.
29.Seed H B, Clough, R W. (1963)”Earthquake resisitance of sloping coredam. “ ASCE Vol.89.No. SM1,
30.Seed, H. B. (1979)“Considerations in the earthquake-resistant design of earthquake and rockfill dams.”Geotechnique, Vol. 29, No.3, pp. 215-263.
31.Taylor, R.N. (1995)”Geotechnical centrifuge technology” Blackie Academic & Professional, London 296p.
32.Terzaghi, K., 1950. Mechanism of landslides , Application of Geology to Engineering Practice (Berkey Volume), Geological Society of America, New York, pp. 83-123.
33.Wartman, J., Riemer, M.F.,k Bray, J. D., and Seed, R. B. (1998), “Newmark Analyses of a Shaking Table Slope stability Experiment. ” Geotechnical Earthquake Engineering and Soil Dynamic. pp. 778-789
34.Wartman, J., Bray, J. D., and Seed, R. B. (2001), “Physical Model　Studies of Seismically Induced Deformations In Slopes”, GeoEngineering Volume: 61-03, Section: B, page: 1541
35.Wartman, J., Bray, J. D., and Seed, R. B. (2005) Shaking Table Modeling of Seismically Induced Deformation in Slopes. Journal of Geotechnical and Geoenvironmental Engineering, Vol. 131, No. 5, pp. 610-622
36.Wakai, A ., Ugai, K,.Sato, M. and Tazo, T (2001), “Numerical Analysis for Seismic Behavior of a Slope Based on a Simple Cyclic Loading model”, Proceedings Fourth International Conference on Recent Advances in Geotechnaical Earthquake Engineering and Soil Dynamics and Symposium in Honor of Professor W.D. Liam Finn, Paper No. 5.06
37.Wieczorek, G.F., Wilson, R.C. and Harp, E.L. (1985), “Map Showing Slope Sstability During Earthquakes in San Mateo County, California”, Miscellaneous Investigations Map I-1257-E., U.S. Geological Survey
38.Yiannis TSOMPANAKIS, Prodromos N. PSARROPOULOS, and Stefanos TSIMPOURAKIS (2007)“Dynamic interaction between retaining wallsand retained structures, 4th International Conference on Earthquake Geotechnical Engineering, 25-28, 
39.Yuzhen Yu .,Lijun Deng.,Xun Sun and He Lu (2008)” Centrifuge modeling of a dry sandy slope response to earthquake loading” Eng Bull Earthquake Eng ,pp:447–461
40.內政部營建署 (2006)，“建築物耐震設計規範及解說”，營建雜誌社。
41.內政部營建署(2002)，”建築技術規則”，營建雜誌社。
42.日本建設省河川局砂防部 (1998). 斜面崩壞防止工事設計の實例—本篇(急傾斜地崩壞防止工事技術指針），第109～131頁。
43.王世榮(2002)，「以位移法分析自然邊坡在地震作用下的平面式破壞」，碩士論文，成功大學資源工程學研究所，台灣，台南。
44.王彥誌(2012) 「以波動方程和有限元素分析樁基礎受震行為之比較」，碩士論文，淡江大學土木工程研究所，台灣，淡水。
45.石塚忠範、高梨和行、草野慎一(2004)，「日本土砂災害防止法介紹」，中華水土保持學報，第35期，295-308頁。
46.台灣世曦工程顧問股份有限公司 (2011)，“特二號道路高架橋梁樁基礎場址地盤鑽探報告”
47.吳宗達(2007)，「不同退縮距離對坡地建物之應力傳遞與變形探討」，碩士論文，台灣科技大學，台灣，台北。
48.朱聖心(2001)，「應用地理資訊系統製作地震及降雨所引致之山崩危險圖」，碩士論文，國立台灣大學土木工程研究所，台灣，台北。
49.余鍾琦(2009) ，「臨海深開挖災變實例數值分析之探討」，碩士論文，高雄第一科技大學營建工程系，台灣，高雄。
50.李姿穎(2009)，「沉箱式碼頭受震動力反應分析之研究」，碩士論文，台灣大學土木研究所，台灣，台北。
51.林澄清(2004)，「 坡地建築物退縮邊坡安全合理距離之研究」中華民國建築學會建築學報，第46期，135~152頁。
52.林美聆(2002)，「 動態作用下之邊坡滑移與破壞分析」，國家科學委員會專題研究計畫成果報告。
53.林美聆(2004)，「區域性邊坡受震反應分析及振動台模型試驗」，國家科學委員會專題研究計畫成果報告。
54.林美聆(2007)，「大型粒狀土壤邊坡振動台模型實驗」，國家科學委員會專題研究計畫成果報告。
55.林美聆(2009)「利用小型振動台模擬邊坡受震情形之研究」，國家科學委員會專題研究計畫成果報告。
56.香港土木工程署土力工程處《Highway Slope Manual》(2000)
57.翁作新(1985)，「 動態土壓之計算與應用」地工技術雜誌，第9期，55~60頁。
58.國家地震工程研究中心(2004)，“耐震與性能設計規範研究（一）”。
59.陳世芳(1979)，”理論土壤力學與實用基礎工程(上)”，文笙書局。
60.陳重光(2011)，「以RIDO及PLAXIS程式模擬支撐擋土開挖行為初探」，碩士論文，國立宜蘭大學，台灣，宜蘭。
61.陳時祖、溫郁菁、彭文飛、蘇容瑩(2002)，「以位移法分析自然邊坡破壞行為之研究」，國科會工程處地震尖端研究計畫群—地震衍生之邊坡破壞行為研究及應用—九十一學年度期中研究成果研討會論文集。
62.陳嬑璇(2002)，「溪頭地區潛感圖製作研究」，碩士論文，國立台灣大學土木工程研究所，台灣，台北。
63.陳榮河，紀柏全(2010)，「模型邊坡試驗之因次分析」地工技術雜誌，第125期，7-14頁。
64.彭文飛(2001)，「以位移法分析自然邊坡在地震時之破壞行為的初步探討」， 碩士論文，國立成功大學資源工程學研究所，台灣，台南。
65.黃鈺珊(2008)，「載重位置對崩積土邊坡位移影響之研究」 ，碩士論文，華梵大學環境與防災設計研究所，台灣，台北。
66.曾惠蒂(2004)，「變形分析法應用於邊坡穩定分析之研究」，碩士論文，國立成功大學資源工程學系研究所，台灣，台南。
67.楊子逸(2009)，“基樁耐震性能分析之初探”，碩士論文，淡江大學土木工程研究所，台灣，淡水。
68.廖南華(2003)，「土壤經驗參數於數值分析之應用」，碩士論文，成功大學土木工程研究所，台灣，台南。
69.鄭世豪(2009)，「聲坡監測應用於地錨結構破壞之模擬和訊號分析」，博士論文，淡江大學土木工程博士班，台灣，淡水。
70.鄭柏威(2008)，「應用PLAXIS程式分析SMW擋土壁特性與開挖變形之關係」 ，碩士論文，台灣科技大學研究所，台灣，台北。
71.鄭錦桐(2002)，「台灣地區地震危害度的不確定性分析與參數拆解」，博士論文，國立中央大學地球物理研究所，台灣，桃園。
72.毆章煜(2009)，“深開挖工程-分析設計理論與實務”，科技圖書。
73.劉哲瑋(2008)，「應用不連續變形分析法模擬受荷重下邊坡之穩定性分析」，碩士論文，高雄第一科技大學，台灣，高雄。
74.羅曉萍(2006)，「九份地區邊坡受震之位移分析」，碩士論文，國立臺北科技大學土木與防災研究所，台灣，台北。