高爾夫球桿頭的重心位置，影響撞擊後球的發射角度、發射初速與球體的旋轉量。本文以重心深度、重心高度與重心距離描述重心位置，改變桿頭薄殼厚度調整重心落點，變更打擊面傾角做桿頭外型變化，最後以三種桿頭初速進行撞擊模擬。藉此探討重心對球體擊發後行為的影響，並且歸納出不同桿頭速度、打擊面傾角、桿頭重心配置與飛行距離的關係。
對於揮桿速度較小的打者，宜選用大深度、低重心、打擊面傾角較大的桿頭，以增加起飛角度、彈道高度與滯空時間，利於加長飛行距離；揮桿速度快的打者，則需降低球的起飛角度，不宜使用打擊面傾角過大的桿頭，且桿頭的重心深度稍淺、重心高度稍高。
本文以有限元素法模擬球頭撞擊球體的行為，得出一系列撞擊後球體之物理特性結果，並以這些數值作為初值代入飛行運動方程式，以數值方法得到飛行軌跡。針對不同打擊習性的打者，可用本文方法做一飛行射程預測，或作為球頭選用之依據。

The center of gravity (CG) in the golf club head is crucial to the launch angle, launch velocity and spin of golf balls after the golf impact. The CG locations in this study refer to the depth of CG, height of CG and distance of CG, which were determined by the change of shell thickness of a golf club head. By means of the finite element method (FEM), the researcher analyzed the impact process for golf balls and golf head clubs. Then he discussed the relationship between club heads and ball behaviors after impact. He also simulated the trajectory by the numerical method on the basis of the impact results of the golf ball.
To sum up, a longer carry was generalized by deeper CG, lower CG and higher loft angle for a golfer with a slow club head speed. A shallower CG, higher CG and lower loft on a club head for a golfer with a fast club head speed were also verified in this study. Finally, a series of impacts produced by FEM and numerical method were generalized. Specifically, the method can be used to predict the flying distance of golf ball. It also can help predict the proper CG locations for golf club head designing.

謝詞 i
中文摘要 ii
Abstract iii
目錄 iv
表次 vi
圖次 vii
第一章 緒論 1
1-1 前言 1
1-2 文獻回顧 2
1-3 研究動機與目的 3
1-4 研究方法 3
第二章 基礎理論 4
2-1 LS-DYNA顯式分析理論 4
2-2 高爾夫球飛行射程推算 10
2-3 以Runge Kutta Method解微分方程式 13
2-4 高爾夫球頭設計尺寸名稱與定義 15
2-5 高爾夫球桿頭重心位置定義 18
2-6 高爾夫球桿頭板殼名稱 20
第三章 有限元素模型與設計參數 21
3-1 碰撞模擬基本假設與有限元素模型 21
3-2 高爾夫球頭模型設計 23
3-3 模型之材料參數 25
3-4 桿頭模型重心控制 26
3-5 有限元素模型 31
第四章 碰撞結果與討論 33
4-1 球體發射速度(launch velocity)與重心關係 33
4-2 球體發射角度(launch angle)與重心關係 42
4-3 後旋量(ωx)與重心關係 48
4-4 側旋量(ωy)與重心關係 54
4-5 射程與球體擊發後行為之關係 60
4-6 射程與重心之關係 66
第五章 結論 79
參考文獻 81
附錄一 83
附錄二 91

[1]Tetsuo Yamaguchi, Ichiro Tominaga, Takuzo Iwatsubo, Nortitoshi Nakagawa, and Masakatsu Akao, “Transfer Characteristic in Collision of Two Viso-Elastomers under Mechanical Impedance,” Theoretical and Appl. Mech., Vol. 34, pp. 153-166, 1985.
[2]Y. Lee, J. F. Hamiltion, and J. W. Sullivan, “The Lumped Parameter Method for Elastic Impact Problems,” J. Appl. Mech., Vol. 53, pp. 1-4, 1986.
[3]Theodore Jorgensen, Jr., “On the Dynamics of a Golf Club,” American journal of physics, Vol. 38, pp. 644-651, 1970.
[4]M. A. Lampsa, “Maximizing Distance of the Golf Drive: An Optimal Control Study,” Transactions of the ASME, pp. 362-367, 1975.
[5]S.K. Cheong, K. W. Kang, S. K. Jeong, “Evaluation of the mechanical performance of golf shafts,” Engineering Failure Analysis, Vol. 13, pp. 464-473, 2006.
[6]Takuzo Iwatsubo, Nobuki Konishi, and Tetsuo Yamaguchi, “Research on Optimum Design of a Golf Club,” Transactions of the Japan Society of Mechanical Engineers, Vol. 56, pp.2386-2392, 1990.
[7]Takuzo Iwatsubo, Shozo Kawamura, Kazuyoshi Miyamoto and Tetsuo Yamaguchi, “Analysis of Golf Impact Phenomenon and Ball Trajectory,” JSME international journal. Series C, Vol. 41, pp. 822-828, 1998.
[8]邱佑宗，高爾夫球頭甜區計算，中華民國第二十二屆全國力學會議，頁 395-400，1998。
[9]Frank Thomas, “The benefits of bulge and roll,” Frank Talk, pp. 68, 2005.
[10]Takuzo Iwatsubo, Shozo Kawamura, Kazuyoshi Miyamoto, and Tetsuo Yamaguchi, “Numerical Analysis of Golf Club Head and Ball at Various Impact Points,” Sports Engineering, pp. 195-204, 2000.
[11]K. Arakawa, T. Mada, H. Komatsu, T. Shimizu, M. Satou, K. Takehara and G. Etoh, “Dynamic Contact Behavior of a Golf Ball during an Oblique Impact,” Experimental Mechanics, Vol. 46, pp. 691–697,2007
[12]John C. Adams Jr., “Flight Mechanics of a Spinning Dimpled Spheroid,” ASME seminar, 2004.
[13]K.R. Williams and B.L. Sih, “Three Dimensional Analysis of Club Movement and Clubhead Orientation at Impact During the Golf Swing,” Journal of Biomechanics, Vol. 40, pp. 317, 2007.
[14]Takeshi Naruo, Taketo Mizota, “Aerodynamic Force Measurement of a Golf Ball and 3D Trajectory Analysis,” Journal of Japan Society of Fluid Mechanics, Vol. 23(3), pp. 203-211, 2003.
[15]張庭彰，高爾夫球桿頭形狀變化之撞擊結果分析，國立中山大學碩士論文，2004。
[16]謝昇宏，高爾夫球桿頭打擊面曲率變化之擊球效果分析，國立中山大學碩士論文，2005。
[17]J. L. Meriam and L. G.. Kraige, Egineering Mechanics Dynamics, 5th ed., Wiley, USA, 2001.
[18]John O. Hallquist, LS-DYNA Theoretical Manual, Livermore Software Technology Corporation, Livermore, 1998.
[19]Daryl L. Logan, A First Course in the Finite Element Method, 4th ed., Thomson Learning, Canada, 2006.
[20]Saeed Moaveni, Finite Element Analysis: Theory and Application with ANSYS, 2nd ed., Upper Saddle River, New Jersey, 2003.
[21]趙海歐，LS-DYNA動力分析指南，兵器工業出版社，2003。