高爾夫球頭的重心變化，影響球體碰撞後，球體的初速度與球的旋轉量，本文以高爾夫球頭容積變化，固定球頭重量為200g的情況下，改變球頭的打擊面水平弧度半徑、砝碼配重，藉此探討球頭水平弧度半徑、砝碼配重，對於球體發射後的影響，歸納出不同情況下，對於球體三維飛行軌跡與偏移距離。
對於球頭容積為400 cc，其球頭的最佳水平半徑弧度為11 in，430 cc宜選用水平半徑弧度為12 in，而460 cc的球頭可以選用水平半徑弧度為13 in 的球頭。
砝碼配重可以有效改善打者習慣，讓打者的擊出的球得到適當的調整。
本文以有限元素軟體分析球受撞擊後，得到球體一系列的撞擊後球體的物理行為，將初始值代入三維飛行運動方程式，以數值分析得到飛行軌跡，針對不同球頭，可用本文的方法作一般球體三維飛行軌跡預測，作為球頭設計時之依據。

Variation of the center of gravity of a golf club head will influence the initial velocity and rotation of speed of a ball after the golf ball is struck by golf club head. After fixing the weight of 200g of a golf head, the researcher changes the volume of golf head and the horizontal curvature of radius. He also distribute counterpoise to investigates the effect of launching of a golf ball. This thesis summarizes the ball of three-dimensional flight trajectory and offset distance.
For the volume of the golf head is 400 cc, the best level of the radius of horizontal curvature is 11 in, in the 430 cc should use a radius of horizontal curvature of 12 in, and the 460 cc head club can chose a radius of horizontal curvature of 13 in. The distribution of counterpoise can effectively improve the play’s habits, so the trajectory of a golf ball can be appropriately adjusted.
By finite element method, the physical behavior of a series of the lunching ball can be predicted. The trajectory of golf ball can be measured by substituting the inertial value of ball into the three-dimension equations of motion. According to the trajectory of golf ball flight by this study, this study provides the characteristics for designing a golf club head.

目錄
論文審定書 i
誌謝 ii
中文摘要 iii
Abstract iv
目錄 v
表次 viii
圖次 ix
第一章 緒論 1
1.1 前言 1
1.2 文獻回顧 1
1.3 研究動機 3
1.4 研究方法 3
第二章 理論基礎 4
2.1 顯式分析與隱式分析 4
2.1.1 顯式 4
2.2 中心差分法 4
2.3 Time收斂準則 5
2.4 高爾夫球飛行經驗理論 6
2.5 RK4解微分方程式 10
2.6 高爾夫球頭尺寸定義與名稱 12
2.7 齒輪效應 16
2.8 開口閉口 19
第三章 有限元素建模與設定 20
3.1 前言 20
3.2 基礎假設 20
3.3 模型之材料參數 21
3.4 高爾夫球頭模型設計 23
3.5 網格收斂分析 26
3.6 模擬座標表示 28
3.7 模型參數設計 29
3.8 有限元素與模型 30
第四章 結果與討論 33
4.1 球頭參數打擊面水平弧度半徑與側旋量 關係 33
4.2 球頭容積改變對甜蜜區的影響 36
4.3 球體撞擊後之飛行軌跡 42
4.4 砝碼配重(重心左右、重心向上)1g、5.5g、10g 48
4.5 最佳設計 52
第五章 結論 55
參考文獻 56

參考文獻
[1]A.Francavilla and O.C.Zienkiewicz , “A Note On Numerical Computation Of Elastic Contact Problem”,Int. J. Num. Mech Engng ,Vol.9,pp.913-924,1975.
[2](1981)T. D. Sachdeva ,C. V. Ramakrishnan and R. Natarnjan ,“ A Finite Element Method For The Elastic Contact Problem” , Trans. ASME J. Engng. Industry Vol. 103,pp. 456-461,1981.
[3]G.R. Johnson ,“High Velocity Impact Calculation In Three Dimension”,Tran.ASME J.Appl. Mech. Vol. 44, pp. 95-106,1977.
[4]N.Okamoto and M. Nakazawa , “Finite Element Incremental Contact Analysis With Various Frictional Conditions” ,Int. J. Num. Meth. Engng ,Vol. 14,pp. 337-357,1979.
[5]G.R. Johnson , “Dynamic Responsion Of Axisymmetric Solids Subjected To Impact And Spin”,AIAA J.,Vol. 17,No. 9,pp. 975-979,1979.
[6]年Y.Lee,J. F. Hamiltion, and J.W.Sullivan,“The Lumped Parameter Method for Elastic Impact Problems”, J. Appl. Mech., Vol.53, p.1-4, 1986.
[7]Ｋ.J. Bathe and A. Chaudhary,“A Solution Method For Planar And Axisymmetric Contact Problems”, Int. J. Num. Meth. Engng.,Vol.21,pp. 65-88,1985.
[8]H.Parisch , “A Consistent Tangent Stiffness Matrix For Three-Dimensional Non-Linear Contact Problems” ,Int. J. Num. Meth. Engng., Vol.28,pp. 1803-1812,1989.
[9]J.E. Mottershead ,S.K.Pascoe and R.G.English ,“A General Finite Element Approach For Contact Stress Analysis” ,Int. J. Num. Meth. Engng.,Vol. 33,pp.765-779,1992.
[10]N.Chandrasekaran ,W.E.Haisler and R.E.Goforth ,“A Finite Element Solution Method For Contact Problems With Friction” , Int. J. Num. Meth. Engng.,Vol. 24,pp.477-495.1987.
[11]Y.Kanto and G. Yagawa ,“ A Dynamic Contact Buckling Analysis By The Penalty Finite Element Method” , Int. J.Num. Meth. Engng., Vol.29,pp. 755-774,1990.
[12]Ghosh ,“Arbitrary Lagrangian-Eulerian Finite Element Analysis Of Large Deformation In Contacting Bodies” ,Int. J. Num. Meth. Engng., Vol.33,pp.1891-1925,1992.
[13]J.C. Simo and T.A Lanrsen ,“An Augrmented Lagrangian Treatment Of Contact Problems Involving Friction”,Computers & Structures , Vol.42,No.1,pp.97-116,1992.
[14]F.Paris ,A Foces and J. A. Garrido ,“Application Of Boundary Element Method To Solve Three-Dimensional Elastic Contact Problems Without Friction” ,Computers & Structures ,Vol.45,No.1,pp.19-33,1992.
[15]A. Gaknaya , D. Lambert and A. Cardou ,“A boundary Element And Mathematical Programming Approach For Frictional Contact Problems”, Computers & Structures ,Vol. 42,No.3,pp. 341-353,1992.
[16]R.J. Neal and L. M. Hubinger“Effect Of Groove Shape On Spin Of A Golf Ball”,Biomechanics J.,Vol.22,No.10,pp.1062,1989.
[17]T.Iwatsubo, N. Nakagawa, M. Akao, I. Tominaga and T.Yamaguchi,“Optimum Design Of A Golf Club Head ”, C Hen/Trans. JSME Vol. 56,No. 524 pp.1053-1059,1990.
[18]Y.Nakao and E.Urata,’’Shape Measurement Of A Three Dimension Body By A Cylindricak Probe Contacting With It Side Surface’’, C Hen/Trans. JSME,Vol. 57,No.544, pp. 4008-4015, 1991.
[19]蕭翰林，高爾夫碰撞現象之分析與杆頭設計，國立中山大學論文，2000。
[20]黃信凱，碳纖維之高爾夫球杆之動態分析，國立中山大學論文，2010。
[21]K. Tanaka1, F. Sato2, H. Oodaira1, Y. Teranishi3, F. Sato3, and S. Ujihashi1, “Construction of the finite-element models of golf balls and simulations of their collisions”, Journal of Materials, Design And Application, Vol. 220, Part I, 2005.
[22]Takeshi NARUO and Taketo MIZOTA, Fukuika Institute Of Technology。
[23]Maltby, Ralph, Golf club design, fitting, alteration & repair: the principles & procedures , second edition