本研究針對隅角開放型青光眼的生成與發病機制建立出完整的數值模型以及進行模擬力學分析探討。隅角開放型青光眼是因為小樑網的病變導致眼睛內房水排出受阻，使眼球前房前端眼壓上升影響到在後房進而導致壓力堆積，後房壓力經由傳遞壓迫到視神經造成視野上的缺陷與模糊。
本研究使用了兩套軟體,3D繪圖軟體PROE與有限元素法分析軟體ANSYS,來進行模型的建立與隅角開放型青光眼相關模擬,此模擬分析可以了解到青光眼內部的壓力分佈情況以及不同杯盤比時應力的聚集之處。 
本研究針對了隅角開放型青光眼不同的杯盤比做應力分佈的探討,針對杯盤比程度0.3-0.7做每個不同模型的研究,在相同條件下進行分析模擬,發現凹陷從0.3-0.7為一種集中靠近視神經盤趨勢,並研究探討每個不同凹陷應力集中之處也不一樣,也針對不同集中之處的不同數值大小做分析探討。
本研究內容對於隅角開放型青光眼的模擬分析結果,可以令青光眼之患者更了解自己相關發病情況,因為杯盤比大小是造成青光眼發病的主要原因之一,能提供未來醫學上一個資訊與具有參考價值的模擬數據。

This study explores open-angle glaucoma and its Pathogenesis to create the whole Numerical model and Mechanical analysis. Open-angle glaucoma is caused by pathological changes of trabecular meshwork. Intraocular aqueous fluid will be stuck and intraocular pressure will increase in both anterior and posterior chamber. Then, intraocular pressure of posterior chamber diffuses and oppress optic disc, causes the blurring and disability of sight. 
This study apply two software, 3D Design software PROE and Finite element method simulation software ANSYS, to create mode and related simulation of open-angle glaucoma. This simulation analysis helps us to understand pressure distribution of glaucoma and where the stress gathers in different Cup-disc ratio.
    This study analyzes stress distribution in different Cup-disc ratio of open-angle glaucoma. I do different researches based on Cup-disc ratio degree from 0.3 to 0.7. Under the same condition, I find from degree 0.3 to 0.7, surface depression shows that stress tends to get close to optic disc. Besides, where the stress gathers differs in different surface depression, so I also conduct numerical analysis based on different stress-gathering section. 
The result of analyzing open-angle glaucoma can help the patient understand more about the cause of open-angle glaucoma. In addition, because Cup-disc ratio degree is the main cause of open-angle glaucoma, this study can also offer simulated data of reference value in medical use in the future.

第一章	序論	1
1-1	前言	1
1-2	研究背景	3
1-3	文獻回顧	8
第二章	研究方法	15
2-1	有限元素法概論	15
2-2	應用軟體	16
2-3	研究過程	18
第三章	分析理論方程式	20
3-1	連續性方程式(CONTINUITY EQUATION)	20
3-2	動量方程式(MOMENTUM EQUATION)	20
3-3	可壓縮能量方程式(COMPRESSIBLE ENERGY EQUATION)	21
3-4	壓力(PRESSURE)	22
3-5	收斂檢測	23
第四章	分析元素、材料參數與邊界條件設定	25
4-1	定義分析元素	25
4-2	設定材料參數	27
4-3	眼球實體模型建立	30
4-4	初始條件設定	35
4-5	負載條件設定	36
第五章	數值模型分析結果	38
5-1	結果分析	38
5-1-1.	流體分析結果-眼球前後房流體模擬分析	38
5-1-2.	結構體分析結果	42
5-2	分析結果驗證	46
5-3	結果討論	57
5-3-1白種人、亞洲人、非洲人之視神經盤結構差異比較探討	57
5-3-2 白種人、亞洲人、非洲人之不同杯盤比探討	62
5-3-3 白種人、亞洲人、非洲人之不同材料參數下比較探討	68
第六章	結論與未來展望	77
6-1	結論	77
6-2	未來展望	79
參考文獻	80

圖目錄
      圖 1-1青光眼視野損毀過程	3
圖 1-2房水流動示意圖[7]	5
圖 1-3 原發性隅角開放型青光眼示意圖[9]	6
圖2-1 ANSYS分析作業流程圖	17
圖2-2研究過程與流程圖	19
圖4-1 FLUID142元素圖[45]	25
圖4-2 SOLID45元素圖[45]	26
圖4-3 眼球前後房模型圖	31
圖4-4 眼球前後房網格圖	31
圖4-5 本研究之眼球模型剖面圖	33
圖4-6 完整眼球之3D模型圖	33
圖4-7完整眼球模型之3D網格圖	34
圖4-8房水流場設定圖	35
圖4-9 眼球結構示意圖[54]	37
圖4-10 邊界條件設定示意圖	37
圖5-1 眼球前後房壓力分佈圖	39
圖5-2 眼球前後房流場分佈圖	39
圖5-3 眼球前後房流線分佈圖	40
圖5-4 隅角開放型青光眼前後房壓力分佈圖	41
圖5-5 隅角開放型青光眼杯盤比0.3之應力分佈	43
圖5-6隅角開放型青光眼杯盤比0.35之應力分佈	43
圖5-7隅角開放型青光眼杯盤比0.4之應力分佈	43
圖5-8隅角開放型青光眼杯盤比0.45之應力分佈	44
圖5-9隅角開放型青光眼杯盤比0.5之應力分佈	44
圖5-10隅角開放型青光眼杯盤比0.55之應力分佈	44
圖5-11隅角開放型青光眼杯盤比0.6之應力分佈	45
圖5-12隅角開放型青光眼杯盤比0.65之應力分佈	45
圖5-13隅角開放型青光眼杯盤比0.7之應力分佈	45
圖5-14隅角開放型青光眼之前後房壓力分佈圖	48
圖5-15隅角開放型青光眼與文獻[21]比較之結構體應力分佈圖	48
圖5-16 文獻[49]MIN 應變峰值對照圖  	52
圖5-17 文獻[49]MIN 應力峰值對照圖	52
圖5-18 文獻[49]MIN 應變均值對照圖	53
圖5-19 文獻[49]MIN 應力均值對照圖	53
圖5-20文獻[49]MAX 應變峰值對照圖	54
圖5-21 文獻[49]MAX 應力峰值對照圖	54
圖5-22文獻[49]MAX 應變均值對照圖	55
圖5-23文獻[49]MAX 應力均值對照圖	55
圖5-24 文獻[51]分析應力圖	56
圖5-25 為白種人、亞洲人、非洲人杯盤比0.3之眼球結構體應力分佈圖	59
圖5-26為白種人、亞洲人、非洲人之眼球結構體應變分佈圖	60
圖5-27 為白種人杯盤比0.3-0.7應力分佈圖	64
圖5-28 為亞洲人杯盤比0.3-0.7應力分佈圖	65
圖5-29 為非洲人杯盤比0.3-0.7應力分	66
圖5 30白種人之少年、中年、老年杯盤比0.3時應力分佈圖	69
圖5 31亞洲人之少年、中年、老年杯盤比0.3時應力分佈圖	70
圖5 32非洲人之少年、中年、老年杯盤比0.3時應力分佈圖	71
圖5 33白種人之少年、中年、老年應變分佈圖	72
圖5 34亞洲人之少年、中年、老年應變分佈圖	73
圖5 35非洲人之少年、中年、老年應變分佈圖	74

1.	http://hospital.kingnet.com.tw/essay/essay.html?pid=18796
2.	K. D.Kallaur, M. Fingeret, Medicla management of the glaucomas. Clinical Ocular Pharmacology, 4th ed., Edited by J. D. Bartlett, S. D. Jaanus, Butterworth Heinemann, Woburn, 2001, pp.831-883.
3.	T. L. Lewis, J. T. Wing, Anatomy and physiology of the optic nerve. Primary care of the glaucomas, Edited by M. Fingeret, T. L. Lewis, 2ed ed., McGraw-Hill, Palatino, Italy, 2001, pp.56-58.
4.	R. L. Stamper, M. F. Lieberman, M. V. Drake, Diagnosis and therapy of the glaucomas, 7th ed., Mosby, Missouri, 1999, pp.177-187.
5.	http://cht.a-hospital.com/w/%E9%9D%92%E5%85%89%E7%9C%BC
6.	http://www.uho.com.tw/sick.asp?aid=950
7.	http://www.dennislameyecenter.com/glaucoma/46.html
8.	Hyun Kyung Cho, Changwon Kee, “Population-based glaucoma prevalence studies in Asians”, SurvOphthalmol. 2014 Jul-Aug;59(4):434-47
9.	http://www.masseyeandear.org/research/ophthalmology/vision-resear
10.	https://www.ntuh.gov.tw/OPH/DocLib10/%E7%9C%BC%E5%A3%93%EF%BC%8D%E5%BD%B1%E9%9F%BF%E7%9A%84%E5%9B%A0%E7%B4%A0%E5%8F%8A%E6%B8%AC%E9%87%8F%E7%9A%84%E6%96%B9%E6%B3%95.aspx
11.	R.C. Tripathi, B.J. Tripathi, “Functional anatomy of the anterior chamber angle”, Biomedical Foundations of Ophthalmology Philadelphia, pp. 197, 1982
12.	MarkJohnson and Kristine Erickson “Mechanisms and Routes of Aqueous Humor Drainage”Principles and Practice of Ophthalmo- logy, Vol. 1,  2000
13.	Heys, Jeffrey J., Victor H. Barocas, and Michael J. Taravella. "Modeling passive mechanical interaction between aqueous humor and iris." TRANSACTIONS-AMERICAN SOCIETY OF MECHANICAL ENGINEERS JOURNAL OF BIOMECHANICAL ENGINEERING 123.6 (2001): 540-547.
14.	JostBruno Jonas, et al. “Optic Disc, Cup and Neuroretinal Rim Size, Configuration and Correlations in Normal Eyes” Investigative Ophthalmology & Visual Science, Vol. 29, No. 7, July 1988
15.	Timothy W. Olsen, et al. “Human Sclera: Thickness and Surface Area” American Journal of Ophthalmology 1998; 125: 237-241
16.	Power, Erik D., et al. "A nonlinear finite element model of the human eye for large deformation loading." Proceedings of the 25th Annual Meeting of Biomechanics, San Diego. 2001.
17.	AdanVillamarin, et al. “3D simulation of the aqueous flow in the human eye” Medical Engineering & Physics 34 (2012) 1462-1470
18.	D. B. Yan, et al. “Deformation of the lamina cribrosa by elevated intraocular pressure” British Fournal of Ophthalmology 1994; 78:643-648
19.	Claude F. Burgoyne, et al. “Measurement of Optic Disc Compliance by Digitized Image Analysis in the Normal Monkey Eye” Ophthalmology 1995; 102: 1790-1799
20.	EiichiUchio,ShigeakiOhno, JojuKudoh, Koki Aoki, Lech Tomasz Kisielewicz “Simulation model of an eyeball based on finite element analysis on a supercomputer”Br J Ophthalmol 1999;83:1106–1111
21.	Antbony J. Bellezza, et al. “The Optic Nerve Head as a Biomechanical Structure. “ Invest Ophthalmol Vis Sci. 2000;41:2991–3000
22.	Esther M. Hoffmann, MD, Linda M. Zangwill, PhD, Jonathan G. Crowston, MD, PhD, and Robert N. Weinreb, MD, "Optic Disk Size and Glaucoma", SurvOphthalmol. 2007 ; 52(1): 32–49.
23.	CamilaZangalli, Shelly R Gupta, George L Spaeth,"The disc as the basis of treatment for glaucoma", Wills Eye Institute, 840 Walnut Street, Suite 1110, Philadelphia, PA 19107, USA.
24.	MichaëlJ.A.Girard, et al. “Peripapillary andPosterior scleral mechanics-Part II: Experimental and inverse finite element characterization” Journal of Biomechanical Engineering: MAY 2009, vol. 131
25.	Richard E. Norman, et al. “Finite element modeling of the human sclera: influence on optic nerve head biomechanics and connections with glaucoma” Experimental Eye Research 93 (2011) 4-12
26.	Younes Agoumi, et al. “Laminar and prelaminar tissue displacement during intraocular pressure elevation in glaucoma patients and healthy controls” Ophthalmology 2011;118:52-59
27.	Sigal, Ian A., et al. "Finite element modeling of optic nerve head biomechanics."Investigative ophthalmology & visual science 45.12 (2004): 4378-4387.
28.	黃士豪，「眼球前後房房水壓力數值模擬分析研究」，淡江大學機械與機電工程研究所，2008。
29.	張又仁，「隅角開放型青光眼之眼球流場與壓力分析」，淡江大學機械與機電工程研究所，2011。
30.	Copt, René-Pierre, Ravi Thomas, and André Mermoud. "Corneal thickness in ocular hypertension, primary open-angle glaucoma, and normal tension glaucoma." Archives of ophthalmology 117.1 (1999): 14.
31.	徐瑋廷，「隅角開放型青光眼之眼球流場與壓力分析」，淡江大學機械與機電工程研究所，2013。
32.	Sander Schutte, et al. “A finite-element analysis model of orbital biomechanics” Vision Research 46 (2006) 1724-1731
33.	李易融，「隅角開放型青光眼之視神經壓力分佈探討」，淡江大學機械與機電工程學研究所，2014。
34.	Jing Wu, M. A. Nasserib, M. Ederc, M. Azqueta Gavaldonc, C. P. Lohmannd and Alois Knollc,"The 3D Eyeball FEA Model with Needle Rotation", APCBEE Procedia 7 ( 2013 ) 4 – 10
35.	Michail E. Kavousanakis, Nikolaos G. Kalogeropoulos, Dimitrios T. Hatziavramidis,"Computational modeling of drug delivery to the posterior eye", Chemical Engineering Science108(2014)203–212
36.	Stitzel, Joel D., et al. "A nonlinear finite element model of the eye with experimental validation for the prediction of globe rupture." SAE CONFERENCE PROCEEDINGS P. SAE; 1999, 2002.
37.	Sigal, Ian A. "Interactions between geometry and mechanical properties on the optic nerve head." Investigative ophthalmology & visual science 50.6 (2009): 2785-2795.
38.	Lai, J. S., and Rita A. Gangwani. "Medication-induced acute angle closure attack." Hong Kong Med J 18.2 (2012): 139-45.
39.	Heys, Jeffrey J. ,Barocas, Victor H. “A boussinesq model of natural convection in the human eye and the formation of Krukenbergs spindle” Annals of Biomedical Engineering. 30(3): 392-401
40.	趙蘭英, 「生醫材料的研究與發展」，工業技術研究院經資中心生醫組, 2002。
41.	Uchio, Eiichi, et al. "Simulation model of an eyeball based on finite element analysis on a supercomputer." British Journal of Ophthalmology 83.10 (1999): 1106-1111.
42.	Von Noorden G K : Summary of the gross anatomy of the extraocular muscles. In Binocular Vision and Ocular Motility, 4thed, Klein EA and Falk KH, editors. St. Louis, C.V. Mosby,1990, pp. 49-50.
43.	Yaacov Shauly, Benjamin Miller “ Tenon’s Capsule: Ultrastructure of Collagen Fibrils in Normals and Infantile Esotropia” Investigative Ophthalmology & Visual Science, Vol. 33, No. 3, March 1992.
44.	Amini R1, Whitcomb JE, Al-Qaisi MK, Akkin T, Jouzdani S, Dorairaj S, Prata T, Illitchev E, Liebmann JM, Ritch R, Barocas VH, "The posterior location of the dilator muscle induces anterior iris bowing during dilation, even in the absence of pupillary block.", Invest Ophthalmol Vis Sci. 2012 Mar 9;53(3):1188-94
45.	http://www.ansys.com/
46.	Ian A. Sigal, et al. “IOP-induced lamina cribrosa deformation and scleral canal expansion: independent or related?” Invest Ophthalmol Vis Sci. 2011;52:9023-9032
47.	Joel D. Stitzel, et al. “A nonlinear finite element model of the eye with experimental validation for the prediction of globe rupture”Stapp Car Crash Journal, Vol. 46 (November 2002)
48.	李世豪，「不同人種及年齡對隅角開放型青光眼之視神經盤壓力影響探討」，淡江大學機械與機電工程學研究所，2015
49.	E. Corona et al “Digital Stereo Image Analyzer for Generating Automated 3-D Measures of Optic Disc Deformation in Glaucoma” IEEE, VOL. 21 No.10 (2002)
50.	I A. Sigal et al “Factors Influencing Optic Nerve Head Biomechanics” Investigative Ophthalmology & Visual Science, Vol. 46, No. 11(2005)
51.	I. A. Sigal et al “A Method to Estimate Biomechanics and Mechanical Properties of Optic Nerve Head Tissues From Parameters Measurable Using Optical Coherence Tomography” IEEE VOL. 33, NO. 6, (2014)