本文主旨在結合電腦輔助設計軟體(Pro Engineer)及有限元素分析軟體(LS-DYNA)分析不同壁厚及外部尺寸的高爾夫球桿頭之碰撞行為，並根據其運動模式設計出較佳的桿頭尺寸。
利用Pro Engineer軟體的參數式特性，迅速地設計不同尺寸外形的桿頭，並提供LS-DYNA作有限元素法的碰撞分析。本研究將設計不同厚度及外部尺寸並在不違反高爾夫球桿頭設計規範內作碰撞分析，並藉由這些尺寸變數所產生與高爾夫球碰撞後運動狀態（速度、角速度等）根據其碰撞結果之物理性質及機械性質，比較不同外形設計對高爾夫球碰撞物理性質影響，並模擬高爾夫球在偏心撞擊所顯現的物理性質，收集統計所有結果資訊，歸納出影響擊球的外部尺寸，進而設計出高擊球飛行初速及寬廣甜蜜點（Sweet Spot）之桿頭。

The purpose of this study is to integrate the efficiency of the computer-aided design software Pro Engineer and the finite element analysis software LS-DYNA on the structural design of the golf club head. By means of the parametric design characteristics of Pro Engineer, various thickness and different shapes of the golf club head have been studied. The impact model will extract the ball velocities and spins for a given hit location on the club face which is analyzed by the software LS-DYNA. After that, the factors for influencing the dynamic analysis, such as the thickness and the shape of the hitting face, the off-center hitting, the ball velocity and spin, the stress distribution, and stress propagation in the head are further identified and measured. Finally, the effects of the factors described above have been analyzed quantitatively to investigate the maximum velocity of the ball and the sweet spot of the golf club head.

中 文 摘 要…………………………………………………i
英 文 摘 要………………………………………………ii
內 容 目 錄………………………………………………iii
圖 目 錄…………………………………………………vi
表 目 錄…………………………………………………ix
第一章、緒論………………………………………………1
1-1 前言…………………………………………………1
1-2 文獻回顧……………………………………………2
1-3 研究動機與目的………………………………………4
1-4 使用之軟體介紹………………………………………7
1-4-1使用之軟體介紹………………………………7
1-4-2前處理器部分（Pre-Processor）………………7
Pro Engineer (release 20.0) …………………7
HyperMesh （Version 3.0）……………………8
FEMB （Version 27）……………………………8
1-4-3求解器部分（Processor）………………………8
LS-DYNA (Version 950) …………………………8
1-4-4後處理器部分（Post-Processor）………………9
LS-TAURUS………………………………………9
PostGL及PostGL-Graph…………………………9
1-5高爾夫桿頭設計要點及規範………………………10
1-6研究方法……………………………………………12
1-7進行步驟………………………………………………14
1-8工作項目及具體成果…………………………………15
第二章、理論分析………………………………………16
2-1前言………………………………………………16
2-2運動方程式的時間積分公式………………………17
2-3 LS-DYNA軟體於有限元素法之運用…………………19
2-4中央差分法(The Central Difference Method) 22
2-4-1中央差分法的穩定性條件與時間增量控制……23
2-4-2摩擦係數……………………………………24
第三章、碰撞模型之建構…………………………………25
3-1前言………………………………………………25
3-2基本假設…………………………………………25
3-3材料性質…………………………………………26
3-4有限元素模型建立過程……………………………27
3-5有限元素模型性質及元素型式說明…………………30
3-6邊界及初始條件……………………………………31
3-7各項設定值…………………………………………34
第四章、案例分析………………………………………34
4-1前言………………………………………………34
4-2案例一平板傾斜角0度與靜摩擦係數及動摩擦係數之關係…34
4-2-1基本假設……………………………………34
4-2-2結果統計與分析……………………………35
4-2-3案例分析與討論……………………………39
4-3案例二.平板傾斜角11度與靜摩擦係數及動摩擦係數之關係40
4-3-1結果統計與分析……………………………40
4-3-2案例分析與討論……………………………45
4-4案例三不同的平板傾斜角之碰撞效應………………47
4-4-1結果統計與分析……………………………47
4-4-2案例分析與討論……………………………52
4-5案例四不同的撞擊物體形狀之碰撞效應……………53
4-5-1基本假設……………………………………53
4-5-2結果統計與分析……………………………54
4-5-3案例分析與討論……………………………58
4-6案例五不同的高爾夫桿頭打擊面厚度之碰撞效應…60
4-6-1結果統計與分析……………………………60
4-6-2案例分析與討論……………………………66
4-7案例六打擊面厚度設計於田口式直交表法之運用…68
4-7-1直交表法設計參數規劃………………………69
4-7-2結果統計與分析……………………………69
4-7-3案例分析與討論……………………………78
4-8案例七甜蜜點之斷定及最佳桿頭擊球區域之探討…80
4-8-1結果統計與分析……………………………80
4-8-2案例分析與討論……………………………89
水平偏心撞擊………………………………89
垂直偏心撞擊………………………………90
第五章、結論……………………………………………93
參考文獻…………………………………………………96

1.Takuzo IWATSUBO, Nobuki KONISHI, and Tetsuo YAMAGUCHI, "Research on Optimum Design of a 高爾夫球桿頭 Club", 日本機械學會論文集(c), Vol. 56, pp. 68-73,1990.
2.Takuzo IWATSUBO, Noritoshi NAKAGAWA, Masakatsu AKAO, Ichiro TOMINNAGA, and Tetsuo YAMAGUCHI, "Optimum Design of a 高爾夫球桿頭 Club Head", 日本機械學會論文集(c), Vol. 56, pp.203-209, 1990.
3.Tetsuo YAMAGUCHI, Ichiro TOMINAGA, Takuzo IWATSUBO, Nortitoshi NAKAGAWA, and Masakatsu AKAO, "Transfer Characteristics in Collision of Two Viso-Elastomers under Mechanical Impedance", Theoretical and Appl. Mech., Vol. 34, pp. 153-166,1985.
4.Charles E. Smith, Pao-Pao Liu, "Coefficient of Restitution", J. Appl. Mech., Vol. 59, pp. 963-969,1992.
5.A.C. Eringen, "Transverse Impact on Beam and Plates", J. Appl. Mech.,1953.
6.Y. Lee, J. F. Hamilton, and J. W. Sullivan, "The Lumped Parameter Method for Elastic Impact Problems", J.Appl. Mech., Vol. 50, pp.823-827,1983.
7.C. Zener, and H. Feshbach, "A method of Calculating Energy Losses During Impact", J. Appl. Mech., June, pp. A67-A70, 1939.
8.J. B. Keller, "Impact With Friction", J. Appl. Mech., Vol. 53, pp.1-4, 1986.
9.Brylawski, A. M(1994) "An investigation of three dimension deformation of a golf club during downswing". Science and golf II: Proceedings of the 1944 World Scientific Congress of golf. pp.265-170.
10.Miline, R.D.(1990)"What is the role of the shaft in the golf swing?" Science and Golf: Proceedings of the 1944 World Scientific Congress of golf. pp.252-257.
11.Winfield,D.C.and Tan,T.E.(1994)"Optimization of clubhead loft and swing elevation angles for maximum distance of a golf drive". Computers & Structures Vol.53,nl,pp19-25.
12.Kane, T K. and Levision, D.A.,(1985)Dynamices: Theory and Applications. McGraw-Hill, NewYork.
13.Takuzo IWATSUBO(1998)"Analysis of Golf Impact Phenomenon and Ball Trajection."JSME Int. Journal., Series C, Vol.41,No.4,pp.823-828.