直昇機之減振設計已是國內外的研究課題之一，在諸多減振法之中，最讓人感興趣的為高階諧模控制法。這些控制輸出端賴於針對某一葉片做轉距角的調整，因此有獨立葉片控制法產生。然而這種利用葉片根部的連桿來改變轉距角的傳統方法已漸漸不被採用了，取而代之的為機翼後緣的伺服襟翼。近代的直昇機學界朝著複合材料之後緣機翼及智慧型結構的應用來發展，而這其中又以分析懸停滯空時之穩定性最具實用價值 (共面運動的緣故)，然而相關後緣襟翼的研究仍停留在實驗的階段，在理論方面，卻是極為缺乏且幾乎為不完整的耦合模式。
    因此，本研究將建立一套整合直昇機主旋翼之非線性複合材料結構葉片、線性ONERA動態失速的空氣動力模式、尾流造成之誘導流及允許旋翼葉片後緣襟翼做任意角度調整之空氣動力流場的耦合模式。並進而分析這一套耦合系統在直昇機懸停時，隨著後緣襟翼角度的改變對於整個主旋翼氣體彈性系統之影響。
    相較於以往之研究，本論文使用完全非線性複材葉片除了能模擬更大的形變量外，亦能從中得到更完整的振動行為，從中也將可以觀察到，在線性結構方程不易被觀察到的現象。
    吾人之耦合系統所使用的比較函數是線性複材葉片方程式之理論解，並假設其邊界條件與非線性複材葉片的邊界條件一致。再利用Galerkin法來求解此耦合系統。本研究將著眼於系統之特徵值及動態響應分析，尤其是升力及扭力矩，亦即求取在不同的後緣襟翼角度之下，整個耦合系統的響應，其結果將可提供做為日後直昇機後緣襟翼設計的資訊及高階諧模控制減振技巧參考之數據。

The dynamic behavior of a long slender beam is of interest in connection with helicopter rotor blades, wind-turbine blades, and other systems that perform large and/or complex motions. When the vibration amplitude is large, various nonlinear effects come into play. Owing to the participation of these nonlinear characteristics, many nonlinear phenomena deserve detailed exam.
    For the past years, the concept of trailing edge flap (TEF) has been applied on the helicopter rotor blade system. Although experimental analyses have made a successfully progress, the theoretical models are still on the way. It is proved that the TEF has great advantages on the vibration reduction and control of the helicopter rotor system, especially for the higher harmonic control (HHC) and individual blade control (IBC). 
    In this research, the fully nonlinear composite rotor blade of a helicopter rotor system is considered. This structural model is further coupled with a linearized ONERA dynamic stall aerodynamic model and the wake induced flow effects are also included. Peters’ generalized deformable blade unsteady aerodynamic theory is used for the formulation of the TEF model on the composite rotor blade. 
    The comparison functions of this aeroelastic system are obtained analytically by a linearized composite blade equation, assuming the boundary conditions is consisting with nonlinear blade. The Galerkin’s method is employed to solve the coupled system. The nonlinear coefficients of these aeroelastic equations are determined by Powell hybrid algorithm and a finite-difference approximation to the Jacobian via IMSL solver. The system eigenvalues and dynamic response due to various TEF locations on the blade and different TEF index angles (pitch angles of TEF) will be studied in this research. Hopefully, several nonlinear dynamics on the composite blade associating with the TEF effects will be discovered in this research.

目錄
中文摘要..............................................Ⅰ
英文摘要..............................................Ⅲ
目錄..................................................Ⅶ
圖目錄................................................Ⅷ
符號說明..........................................ⅩⅩⅦ
第一章  緒論...........................................1
一、1  研究動機........................................1
一、2  文獻回顧........................................4
第二章  複材彈性葉片與尾流耦合模式之建立..............11
二、1  複材葉片結構動力方程之推導.....................11
二、2  葉片上空氣動力表示式之推導.....................16
二、2.1  葉片上的空氣動力表示式之推導.................16
二、2.2  後緣襟翼之空氣動力流場.......................21
二、3  尾流函數.......................................29
二、4  耦合模式的建立.................................31
第三章  複材葉片耦合方程之求解........................33
三、1  非線性耦合方程式之展開.........................33
三、2  複材葉片與尾流耦合方程式.......................34
三、3  非線性耦合方程式之求解.........................35
三、4  動態耦合方程式之求解...........................35
第四章  複材葉片振動模態之建立........................48
四、1  振動模態之建立方法.............................48
四、2  線性等向材料葉片之振動模態.....................55
四、2.1  擺振向之振動模態.............................56
四、2.2  拍撲向之振動模態.............................58
四、2.3  扭曲向之振動模態.............................61
四、3  複合材料葉片之振動模態.........................62
四、3.1  Pai之振動模態................................63
四、3.2  前人之振動模態...............................65
第五章  結果與討論....................................69
五、1  驗證比較之結果.................................70
五、2  外力響應之分析.................................71
五、3  等向均質葉片之非線性與線性之比較...............75
五、4  複材葉片之比較.................................77
五、5  後緣襟翼翼剖面弦向位置之影響...................81
五、6  後緣襟翼展弦鉸接位置之影響.....................84
五、7  後緣襟翼展弦大小之影響.........................88
第六章  結論..........................................93
參考文獻..............................................96
附錄一  NACA 0012葉片翼剖面動態失速之函數參數........102
附錄二  後緣襟翼之相關參數定義.......................103

表目錄
表一  旋翼葉片的參考數據 [8].........................104
表二  環氧化碳複材桿件之材料特性.....................105

圖目錄
圖一  空氣動力的座標.................................106
圖二 (a)  後緣襟翼幾何示意圖-後緣襟翼截面 [14].......107
圖二 (b)  後緣襟翼幾何示意圖-翼展....................107
圖三  T.E.F.在不同弦向位置時的扭力矩與Garrick之比較..108
圖四 (a)  不同襟翼角時，後緣襟翼的扭力矩與實驗值及類線性值之比較圖..........................109
圖四 (b)  不同轉速時，後緣襟翼的扭力矩與實驗值及類線性值之比較圖............................109
圖四 (c)  不同襟翼角時，後緣襟翼的升力與類線性值之比較圖....................................110
圖四 (d)  不同轉速時，後緣襟翼的升力與類線性值之比較圖......................................110
圖五 (a)  等向均質葉片擺振向第一振動模態之外力響應圖........................................111
圖五 (b)  等向均質葉片擺振向第二振動模態之外力響應圖........................................111
圖五 (c)  等向均質葉片擺振向第三振動模態之外力響應圖........................................112
圖五 (d)  等向均質葉片拍撲向第一振動模態之外力響應圖........................................112
圖五 (e)  等向均質葉片拍撲向第二振動模態之外力響應圖........................................113
圖五 (f)  等向均質葉片拍撲向第三振動模態之外力響應圖........................................113
圖五 (g)  等向均質葉片扭曲向第一振動模態之外力響應圖........................................114
圖五 (h)  等向均質葉片扭曲向第二振動模態之外力響應圖........................................114
圖五 (i)  等向均質葉片扭曲向第三振動模態之外力響應圖........................................115
圖六 (a)  Pai複材葉片擺振向第一振動模態之外力響應圖.........................................115
圖六 (b)  Pai複材葉片擺振向第二振動模態之外力響應圖.........................................116
圖六 (c)  Pai複材葉片擺振向第三振動模態之外力響應圖.........................................116
圖六 (d)  Pai複材葉片拍撲向第一振動模態之外力響應圖.........................................117
圖六 (e)  Pai複材葉片拍撲向第二振動模態之外力響應圖.........................................117
圖六 (f)  Pai複材葉片拍撲向第三振動模態之外力響應圖.........................................118
圖六 (g)  Pai複材葉片扭曲向第一振動模態之外力響應圖.........................................118
圖六 (h)  Pai複材葉片扭曲向第二振動模態之外力響應圖.........................................119
圖六 (i)  Pai複材葉片扭曲向第三振動模態之外力響應圖.........................................119
圖七 (a)  D11=0.25,D12=0.15之複材葉片擺振向第一振動模態之外力響應圖.........................120
圖七 (b)  D11=0.25,D12=0.15之複材葉片擺振向第二振動模態之外力響應圖.........................120
圖七 (c)  D11=0.25,D12=0.15之複材葉片擺振向第三振動模態之外力響應圖.........................121
圖七 (d)  D11=0.25,D12=0.15之複材葉片拍撲向第一振動模態之外力響應圖.........................121
圖七 (e)  D11=0.25,D12=0.15之複材葉片拍撲向第二振動模態之外力響應圖.........................122
圖七 (f)  D11=0.25,D12=0.15之複材葉片拍撲向第三振動模態之外力響應圖.........................122
圖七 (g)  D11=0.25,D12=0.15之複材葉片扭曲向第一振動模態之外力響應圖.........................123
圖七 (h)  D11=0.25,D12=0.15之複材葉片扭曲向第二振動模態之外力響應圖.........................123
圖七 (i)  D11=0.25,D12=0.15之複材葉片扭曲向第三振動模態之外力響應圖.........................124
圖八 (a)  D11=0.30,D12=0.15之複材葉片擺振向第一振動模態之外力響應圖.........................124
圖八 (b)  D11=0.30,D12=0.15之複材葉片擺振向第二振動模態之外力響應圖.........................125
圖八 (c)  D11=0.30,D12=0.15之複材葉片擺振向第三振動模態之外力響應圖.........................125
圖八 (d)  D11=0.30,D12=0.15之複材葉片拍撲向第一振動模態之外力響應圖.........................126
圖八 (e)  D11=0.30,D12=0.15之複材葉片拍撲向第二振動模態之外力響應圖.........................126
圖八 (f)  D11=0.30,D12=0.15之複材葉片拍撲向第三振動模態之外力響應圖.........................127
圖八 (g)  D11=0.30,D12=0.15之複材葉片扭曲向第一振動模態之外力響應圖.........................127
圖八 (h)  D11=0.30,D12=0.15之複材葉片扭曲向第二振動模態之外力響應圖.........................128
圖八 (i)  D11=0.30,D12=0.15之複材葉片扭曲向第三振動模態之外力響應圖.........................128
圖九 (a)  D11=0.30,D12=0.20之複材葉片擺振向第一振動模態之外力響應圖.........................129
圖九 (b)  D11=0.30,D12=0.20之複材葉片擺振向第二振動模態之外力響應圖.........................129
圖九 (c)  D11=0.30,D12=0.20之複材葉片擺振向第三振動模態之外力響應圖.........................130
圖九 (d)  D11=0.30,D12=0.20之複材葉片拍撲向第一振動模態之外力響應圖.........................130
圖九 (e)  D11=0.30,D12=0.20之複材葉片拍撲向第二振動模態之外力響應圖.........................131
圖九 (f)  D11=0.30,D12=0.20之複材葉片拍撲向第三振動模態之外力響應圖.........................131
圖九 (g)  D11=0.30,D12=0.20之複材葉片扭曲向第一振動模態之外力響應圖.........................132
圖九 (h)  D11=0.30,D12=0.20之複材葉片扭曲向第二振動模態之外力響應圖.........................132
圖九 (i)  D11=0.30,D12=0.20之複材葉片扭曲向第三振動模態之外力響應圖.........................133
圖十 (a)  等向均質葉片擺振向第一振動模態之外力響應3-D圖.....................................134
圖十 (b)  等向均質葉片擺振向第一振動模態之外力響應2-D圖.....................................134
圖十 (c)  等向均質葉片擺振向第二振動模態之外力響應3-D圖.....................................135
圖十 (d)  等向均質葉片擺振向第二振動模態之外力響應2-D圖.....................................135
圖十 (e)  等向均質葉片擺振向第三振動模態之外力響應3-D圖.....................................136
圖十 (f)  等向均質葉片擺振向第三振動模態之外力響應2-D圖.....................................136
圖十 (g)  等向均質葉片拍撲向第一振動模態之外力響應3-D圖.....................................137
圖十 (h)  等向均質葉片拍撲向第一振動模態之外力響應2-D圖.....................................137
圖十 (i)  等向均質葉片拍撲向第二振動模態之外力響應3-D圖.....................................138
圖十 (j)  等向均質葉片拍撲向第二振動模態之外力響應2-D圖.....................................138
圖十 (k)  等向均質葉片拍撲向第三振動模態之外力響應3-D圖.....................................139
圖十 (l)  等向均質葉片拍撲向第三振動模態之外力響應2-D圖.....................................139
圖十 (m)  等向均質葉片扭曲向第一振動模態之外力響應3-D圖.....................................140
圖十 (n)  等向均質葉片扭曲向第一振動模態之外力響應2-D圖.....................................140
圖十 (o)  等向均質葉片扭曲向第二振動模態之外力響應3-D圖.....................................141
圖十 (p)  等向均質葉片扭曲向第二振動模態之外力響應2-D圖.....................................141
圖十 (q)  等向均質葉片扭曲向第三振動模態之外力響應3-D圖.....................................142
圖十 (r)  等向均質葉片扭曲向第三振動模態之外力響應2-D圖.....................................142
圖十一 (a)  Pai複材葉片擺振向第一振動模態之外力響應3-D圖....................................143
圖十一 (b)  Pai複材葉片擺振向第一振動模態之外力響應2-D圖....................................143
圖十一 (c)  Pai複材葉片擺振向第二振動模態之外力響應3-D圖....................................144
圖十一 (d)  Pai複材葉片擺振向第二振動模態之外力響應2-D圖....................................144
圖十一 (e)  Pai複材葉片擺振向第三振動模態之外力響應3-D圖....................................145
圖十一 (f)  Pai複材葉片擺振向第三振動模態之外力響應2-D圖....................................145
圖十一 (g)  Pai複材葉片拍撲向第一振動模態之外力響應3-D圖....................................146
圖十一 (h)  Pai複材葉片拍撲向第一振動模態之外力響應2-D圖....................................146
圖十一 (i)  Pai複材葉片拍撲向第二振動模態之外力響應3-D圖....................................147
圖十一 (j)  Pai複材葉片拍撲向第二振動模態之外力響應2-D圖....................................147
圖十一 (k)  Pai複材葉片拍撲向第三振動模態之外力響應3-D圖....................................148
圖十一 (l)  Pai複材葉片拍撲向第三振動模態之外力響應2-D圖....................................148
圖十一 (m)  Pai複材葉片扭曲向第一振動模態之外力響應3-D圖....................................149
圖十一 (n)  Pai複材葉片扭曲向第一振動模態之外力響應2-D圖....................................149
圖十一 (o)  Pai複材葉片扭曲向第二振動模態之外力響應3-D圖....................................150
圖十一 (p)  Pai複材葉片扭曲向第二振動模態之外力響應2-D圖....................................150
圖十一 (q)  Pai複材葉片扭曲向第三振動模態之外力響應3-D圖....................................151
圖十一 (r)  Pai複材葉片扭曲向第三振動模態之外力響應2-D圖....................................151
圖十二 (a)  D11=0.25,D12=0.15之複材葉片擺振向第一振動模態之外力響應3-D圖....................152
圖十二 (b)  D11=0.25,D12=0.15之複材葉片擺振向第一振動模態之外力響應2-D圖....................152
圖十二 (c)  D11=0.25,D12=0.15之複材葉片擺振向第二振動模態之外力響應3-D圖....................153
圖十二 (d)  D11=0.25,D12=0.15之複材葉片擺振向第二振動模態之外力響應2-D圖....................153
圖十二 (e)  D11=0.25,D12=0.15之複材葉片擺振向第三振動模態之外力響應3-D圖....................154
圖十二 (f)  D11=0.25,D12=0.15之複材葉片擺振向第三振動模態之外力響應2-D圖....................154
圖十二 (g)  D11=0.25,D12=0.15之複材葉片拍撲向第一振動模態之外力響應3-D圖....................155
圖十二 (h)  D11=0.25,D12=0.15之複材葉片拍撲向第一振動模態之外力響應2-D圖....................155
圖十二 (i)  D11=0.25,D12=0.15之複材葉片拍撲向第二振動模態之外力響應3-D圖....................156
圖十二 (j)  D11=0.25,D12=0.15之複材葉片拍撲向第二振動模態之外力響應2-D圖....................156
圖十二 (k)  D11=0.25,D12=0.15之複材葉片拍撲向第三振動模態之外力響應3-D圖....................157
圖十二 (l)  D11=0.25,D12=0.15之複材葉片拍撲向第三振動模態之外力響應2-D圖....................157
圖十二 (m)  D11=0.25,D12=0.15之複材葉片扭曲向第一振動模態之外力響應3-D圖....................158
圖十二 (n)  D11=0.25,D12=0.15之複材葉片扭曲向第一振動模態之外力響應2-D圖....................158
圖十二 (o)  D11=0.25,D12=0.15之複材葉片扭曲向第二振動模態之外力響應3-D圖....................159
圖十一 (p)  D11=0.25,D12=0.15之複材葉片扭曲向第二振動模態之外力響應2-D圖....................159
圖十二 (q)  D11=0.25,D12=0.15之複材葉片扭曲向第三振動模態之外力響應3-D圖....................160
圖十二 (r)  D11=0.25,D12=0.15之複材葉片扭曲向第三振動模態之外力響應2-D圖....................160
圖十三 (a)  D11=0.30,D12=0.15之複材葉片擺振向第一振動模態之外力響應3-D圖....................161
圖十三 (b)  D11=0.30,D12=0.15之複材葉片擺振向第一振動模態之外力響應2-D圖....................161
圖十三 (c)  D11=0.30,D12=0.15之複材葉片擺振向第二振動模態之外力響應3-D圖....................162
圖十三 (d)  D11=0.30,D12=0.15之複材葉片擺振向第二振動模態之外力響應2-D圖....................162
圖十三 (e)  D11=0.30,D12=0.15之複材葉片擺振向第三振動模態之外力響應3-D圖....................163
圖十三 (f)  D11=0.30,D12=0.15之複材葉片擺振向第三振動模態之外力響應2-D圖....................163
圖十三 (g)  D11=0.30,D12=0.15之複材葉片拍撲向第一振動模態之外力響應3-D圖....................164
圖十三 (h)  D11=0.30,D12=0.15之複材葉片拍撲向第一振動模態之外力響應2-D圖....................164
圖十三 (i)  D11=0.30,D12=0.15之複材葉片拍撲向第二振動模態之外力響應3-D圖....................165
圖十三 (j)  D11=0.30,D12=0.15之複材葉片拍撲向第二振動模態之外力響應2-D圖....................165
圖十三 (k)  D11=0.30,D12=0.15之複材葉片拍撲向第三振動模態之外力響應3-D圖....................166
圖十三 (l)  D11=0.30,D12=0.15之複材葉片拍撲向第三振動模態之外力響應2-D圖....................166
圖十三 (m)  D11=0.30,D12=0.15之複材葉片扭曲向第一振動模態之外力響應3-D圖....................167
圖十三 (n)  D11=0.30,D12=0.15之複材葉片扭曲向第一振動模態之外力響應2-D圖....................167
圖十三 (o)  D11=0.30,D12=0.15之複材葉片扭曲向第二振動模態之外力響應3-D圖....................168
圖十三 (p)  D11=0.30,D12=0.15之複材葉片扭曲向第二振動模態之外力響應2-D圖....................168
圖十三 (q)  D11=0.30,D12=0.15之複材葉片扭曲向第三振動模態之外力響應3-D圖....................169
圖十三 (r)  D11=0.30,D12=0.15之複材葉片扭曲向第三振動模態之外力響應2-D圖....................169
圖十四 (a)  D11=0.30,D12=0.20之複材葉片擺振向第一振動模態之外力響應3-D圖....................170
圖十四 (b)  D11=0.30,D12=0.20之複材葉片擺振向第一振動模態之外力響應2-D圖....................170
圖十四 (c)  D11=0.30,D12=0.20之複材葉片擺振向第二振動模態之外力響應3-D圖....................171
圖十四 (d)  D11=0.30,D12=0.20之複材葉片擺振向第二振動模態之外力響應2-D圖....................171
圖十四 (e)  D11=0.30,D12=0.20之複材葉片擺振向第三振動模態之外力響應3-D圖....................172
圖十四 (f)  D11=0.30,D12=0.20之複材葉片擺振向第三振動模態之外力響應2-D圖....................172
圖十四 (g)  D11=0.30,D12=0.20之複材葉片拍撲向第一振動模態之外力響應3-D圖....................173
圖十四 (h)  D11=0.30,D12=0.20之複材葉片拍撲向第一振動模態之外力響應2-D圖....................173
圖十四 (i)  D11=0.30,D12=0.20之複材葉片拍撲向第二振動模態之外力響應3-D圖....................174
圖十四 (j)  D11=0.30,D12=0.20之複材葉片拍撲向第二振動模態之外力響應2-D圖....................174
圖十四 (k)  D11=0.30,D12=0.20之複材葉片拍撲向第三振動模態之外力響應3-D圖....................175
圖十四 (l)  D11=0.30,D12=0.20之複材葉片拍撲向第三振動模態之外力響應2-D圖....................175
圖十四 (m)  D11=0.30,D12=0.20之複材葉片扭曲向第一振動模態之外力響應3-D圖....................176
圖十四 (n)  D11=0.30,D12=0.20之複材葉片扭曲向第一振動模態之外力響應2-D圖....................176
圖十四 (o)  D11=0.30,D12=0.20之複材葉片扭曲向第二振動模態之外力響應3-D圖....................177
圖十四 (p)  D11=0.30,D12=0.20之複材葉片扭曲向第二振動模態之外力響應2-D圖....................177
圖十四 (q)  D11=0.30,D12=0.20之複材葉片扭曲向第三振動模態之外力響應3-D圖....................178
圖十四 (r)  D11=0.30,D12=0.20之複材葉片扭曲向第三振動模態之外力響應2-D圖....................178
圖十五 (a)  α=0.0°、不同襟翼角時，非線性模式與線性模式整個葉片的升力之比較.................179
圖十五 (b)  α=0.0°、不同襟翼角時，非線性模式與線性模式整個葉片的扭力矩之比較...............179
圖十六 (a)  α=1.2°、不同襟翼角時，非線性模式與線性模式整個葉片的升力之比較.................180
圖十六 (b)  α=1.2°、不同襟翼角時，非線性模式與線性模式整個葉片的扭力矩之比較...............180
圖十七 (a)  α=1.6°、不同襟翼角時，非線性模式與線性模式整個葉片的升力之比較.................181
圖十七 (b)  α=1.6°、不同襟翼角時，非線性模式與線性模式整個葉片的扭力矩之比較...............181
圖十八 (a)  α=3.5°、不同襟翼角時，非線性模式與線性模式整個葉片的升力之比較.................182
圖十八 (b)  α=3.5°、不同襟翼角時，非線性模式與線性模式整個葉片的扭力矩之比較...............182
圖十九  線性穩態空氣動力模式之升力分佈......................................................183
圖二十 (a)  α=0.0°、不同襟翼角時，複材葉片整個葉片的升力之比較 (固定D12，改變D11)..........184
圖二十 (b)  α=0.0°、不同襟翼角時，複材葉片整個葉片的扭力矩之比較 (固定D12，改變D11)........184
圖二十一 (a)  α=3.5、不同襟翼角時，複材葉片整個葉片的升力之比較 (固定D12，改變D11).........185
圖二十一 (b)  α=3.5、不同襟翼角時，複材葉片整個葉片的扭力矩之比較 (固定D12，改變D11).......185
圖二十二 (a)  α=5.0、不同襟翼角時，複材葉片整個葉片的升力之比較 (固定D12，改變D11).........186
圖二十二 (b)  α=5.0、不同襟翼角時，複材葉片整個葉片的扭力矩之比較 (固定D12，改變D11).......186
圖二十三 (a)  α=0.0、不同襟翼角時，複材葉片整個葉片的升力之比較 (固定D11，改變D12).........187
圖二十三 (b)  α=0.0、不同襟翼角時，複材葉片整個葉片的扭力矩之比較 (固定D11，改變D12).......187
圖二十四 (a)  α=3.5、不同襟翼角時，複材葉片整個葉片的升力之比較 (固定D11，改變D12).........188
圖二十四 (b)  α=3.5、不同襟翼角時，複材葉片整個葉片的扭力矩之比較 (固定D11，改變D12).......188
圖二十五 (a)  α=5.0、不同襟翼角時，複材葉片整個葉片的升力之比較 (固定D11，改變D12).........189
圖二十五 (b)  α=3.0、不同襟翼角時，複材葉片整個葉片的扭力矩之比較 (固定D11，改變D12).......189
圖二十六 (a)  等向均質葉片在α=3.5°、不同翼剖面弦向位置時，整個葉片的升力...................190
圖二十六 (b)  等向均質葉片在α=3.5°、不同翼剖面弦向位置時，整個葉片的扭力矩.................190
圖二十六 (c)  等向均質葉片在α=3.5°、不同襟翼角時，整個葉片的升力...........................191
圖二十六 (d)  等向均質葉片在α=3.5°、不同襟翼角時，整個葉片的扭力矩.........................191
圖二十七 (a)  等向均質葉片在α=5.0°、不同翼剖面弦向位置時，整個葉片的升力...................192
圖二十七 (b)  等向均質葉片在α=5.0°、不同翼剖面弦向位置時，整個葉片的扭力矩.................192
圖二十七 (c)  等向均質葉片在α=5.0°、不同襟翼角時，整個葉片的升力...........................193
圖二十七 (d)  等向均質葉片在α=5.0°、不同襟翼角時，整個葉片的扭力矩.........................193
圖二十八 (a)  等向均質葉片在α=3.5°、不同展弦鉸接位置時，整個葉片的升力.....................194
圖二十八 (b)  等向均質葉片在α=3.5°、不同展弦鉸接位置時，整個葉片的扭力矩...................194
圖二十八 (c)  等向均質葉片在α=3.5°、不同襟翼角時，整個葉片的升力...........................195
圖二十八 (d)  等向均質葉片在α=3.5°、不同襟翼角時，整個葉片的扭力矩.........................195
圖二十九 (a)  等向均質葉片在α=5.0°、不同展弦鉸接位置時，整個葉片的升力.....................196
圖二十九 (b)  等向均質葉片在α=5.0°、不同展弦鉸接位置時，整個葉片的扭力矩...................196
圖二十九 (c)  等向均質葉片在α=5.0°、不同襟翼角時，整個葉片的升力...........................197
圖二十九 (d)  等向均質葉片在α=5.0°、不同襟翼角時，整個葉片的扭力矩.........................197
圖三十 (a)  等向均質葉片在α=3.5°、不同展弦大小時，整個葉片的升力...........................198
圖三十 (b)  等向均質葉片在α=3.5°、不同展弦大小時，整個葉片的扭力矩.........................198
圖三十 (c)  等向均質葉片在α=3.5°、不同襟翼角時，整個葉片的升力.............................199
圖三十 (d)  等向均質葉片在α=3.5°、不同襟翼角時，整個葉片的扭力矩...........................199
圖三十一 (a)  等向均質葉片在α=5.0°、不同展弦大小時，整個葉片的升力.........................200
圖三十一 (b)  等向均質葉片在α=5.0°、不同展弦大小時，整個葉片的扭力矩.......................200
圖三十一 (c)  等向均質葉片在α=5.0°、不同襟翼角時，整個葉片的升力...........................201
圖三十一 (d)  等向均質葉片在α=5.0°、不同襟翼角時，整個葉片的扭力矩.........................201
圖三十二 (a)  Pai複材葉片在α=3.5°、不同翼剖面弦向位置時，整個葉片的升力....................202
圖三十二 (b)  Pai複材葉片在α=3.5°、不同翼剖面弦向位置時，整個葉片的扭力矩..................202
圖三十二 (c)  Pai複材葉片在α=3.5°、不同襟翼角時，整個葉片的升力............................203
圖三十二 (d)  Pai複材葉片在α=3.5°、不同襟翼角時，整個葉片的扭力矩..........................203
圖三十三 (a)  Pai複材葉片在α=5.0°、不同翼剖面弦向位置時，整個葉片的升力....................204
圖三十三 (b)  Pai複材葉片在α=5.0°、不同翼剖面弦向位置時，整個葉片的扭力矩..................204
圖三十三 (c)  Pai複材葉片在α=5.0°、不同襟翼角時，整個葉片的升力............................205
圖三十三 (d)  Pai複材葉片在α=5.0°、不同襟翼角時，整個葉片的扭力矩..........................205
圖三十四 (a)  Pai複材葉片在α=3.5°、不同展弦鉸接位置時，整個葉片的升力......................206
圖三十四 (b)  Pai複材葉片在α=3.5°、不同展弦鉸接位置時，整個葉片的扭力矩....................206
圖三十四 (c)  Pai複材葉片在α=3.5°、不同襟翼角時，整個葉片的升力............................207
圖三十四 (d)  Pai複材葉片在α=3.5°、不同襟翼角時，整個葉片的扭力矩..........................207
圖三十五 (a)  Pai複材葉片在α=5.0°、不同展弦鉸接位置時，整個葉片的升力......................208
圖三十五 (b)  Pai複材葉片在α=5.0°、不同展弦鉸接位置時，整個葉片的扭力矩....................208
圖三十五 (c)  Pai複材葉片在α=5.0°、不同襟翼角時，整個葉片的升力............................209
圖三十五 (d)  Pai複材葉片在α=5.0°、不同襟翼角時，整個葉片的扭力矩..........................209
圖三十六 (a)  Pai複材葉片在α=3.5°、不同展弦大小時，整個葉片的升力..........................210
圖三十六 (b)  Pai複材葉片在α=3.5°、不同展弦大小時，整個葉片的扭力矩........................210
圖三十六 (c)  Pai複材葉片在α=3.5°、不同襟翼角時，整個葉片的升力............................211
圖三十六 (d)  Pai複材葉片在α=3.5°、不同襟翼角時，整個葉片的扭力矩..........................211
圖三十七 (a)  Pai複材葉片在α=5.0°、不同展弦大小時，整個葉片的升力..........................212
圖三十七 (b)  Pai複材葉片在α=5.0°、不同展弦大小時，整個葉片的扭力矩........................212
圖三十七 (c)  Pai複材葉片在α=5.0°、不同襟翼角時，整個葉片的升力............................213
圖三十七 (d)  Pai複材葉片在α=5.0°、不同襟翼角時，整個葉片的扭力矩..........................213
圖三十八 (a)  D11=0.25,D12=0.15之複材葉片在α=3.5°、不同翼剖面弦向位置時，整個葉片的升力....214
圖三十八 (b)  D11=0.25,D12=0.15之複材葉片在α=3.5°、不同翼剖面弦向位置時，整個葉片的扭力矩..214
圖三十八 (c)  D11=0.25,D12=0.15之複材葉片在α=3.5°、不同襟翼角時，整個葉片的升力............215
圖三十八 (d)  D11=0.25,D12=0.15之複材葉片在α=3.5°、不同襟翼角時，整個葉片的扭力矩..........215
圖三十九 (a)  D11=0.25,D12=0.15之複材葉片在α=5.0°、不同翼剖面弦向位置時，整個葉片的升力....216
圖三十九 (b)  D11=0.25,D12=0.15之複材葉片在α=5.0°、不同翼剖面弦向位置時，整個葉片的扭力矩..216
圖三十九 (c)  D11=0.25,D12=0.15之複材葉片在α=5.0°、不同襟翼角時，整個葉片的升力............217
圖三十九 (d)  D11=0.25,D12=0.15之複材葉片在α=5.0°、不同襟翼角時，整個葉片的扭力矩..........217
圖四十 (a)  D11=0.25,D12=0.15之複材葉片在α=3.5°、不同展弦鉸接位置時，整個葉片的升力........218
圖四十 (b)  D11=0.25,D12=0.15之複材葉片在α=3.5°、不同展弦鉸接位置時，整個葉片的扭力矩......218
圖四十 (c)  D11=0.25,D12=0.15之複材葉片在α=3.5°、不同襟翼角時，整個葉片的升力..............219
圖四十 (d)  D11=0.25,D12=0.15之複材葉片在α=3.5°、不同襟翼角時，整個葉片的扭力矩............219
圖四十一 (a)  D11=0.25,D12=0.15之複材葉片在α=5.0°、不同展弦鉸接位置時，整個葉片的升力......220
圖四十一 (b)  D11=0.25,D12=0.15之複材葉片在α=5.0°、不同展弦鉸接位置時，整個葉片的扭力矩....220
圖四十一 (c)  D11=0.25,D12=0.15之複材葉片在α=5.0°、不同襟翼角時，整個葉片的升力............221
圖四十一 (d)  D11=0.25,D12=0.15之複材葉片在α=5.0°、不同襟翼角時，整個葉片的扭力矩..........221
圖四十二 (a)  D11=0.25,D12=0.15之複材葉片在α=3.5°、不同展弦大小時，整個葉片的升力..........222
圖四十二 (b)  D11=0.25,D12=0.15之複材葉片在α=3.5°、不同展弦大小時，整個葉片的扭力矩........222
圖四十二 (c)  D11=0.25,D12=0.15之複材葉片在α=3.5°、不同襟翼角時，整個葉片的升力............223
圖四十二 (d)  D11=0.25,D12=0.15之複材葉片在α=3.5°、不同襟翼角時，整個葉片的扭力矩..........223
圖四十三 (a)  D11=0.25,D12=0.15之複材葉片在α=5.0°、不同展弦大小時，整個葉片的升力..........224
圖四十三 (b)  D11=0.25,D12=0.15之複材葉片在α=5.0°、不同展弦大小時，整個葉片的扭力矩........224
圖四十三 (c)  D11=0.25,D12=0.15之複材葉片在α=5.0°、不同襟翼角時，整個葉片的升力............225
圖四十三 (d)  D11=0.25,D12=0.15之複材葉片在α=5.0°、不同襟翼角時，整個葉片的扭力矩..........225
圖四十四 (a)  D11=0.30,D12=0.15之複材葉片在α=3.5°、不同翼剖面弦向位置時，整個葉片的升力....226
圖四十四 (b)  D11=0.30,D12=0.15之複材葉片在α=3.5°、不同翼剖面弦向位置時，整個葉片的扭力矩..226
圖四十四 (c)  D11=0.30,D12=0.15之複材葉片在α=3.5°、不同襟翼角時，整個葉片的升力............227
圖四十四 (d)  D11=0.30,D12=0.15之複材葉片在α=3.5°、不同襟翼角時，整個葉片的扭力矩..........227
圖四十五 (a)  D11=0.30,D12=0.15之複材葉片在α=5.0°、不同翼剖面弦向位置時，整個葉片的升力....228
圖四十五 (b)  D11=0.30,D12=0.15之複材葉片在α=5.0°、不同翼剖面弦向位置時，整個葉片的扭力矩..228
圖四十五 (c)  D11=0.30,D12=0.15之複材葉片在α=5.0°、不同襟翼角時，整個葉片的升力............229
圖四十五 (d)  D11=0.30,D12=0.15之複材葉片在α=5.0°、不同襟翼角時，整個葉片的扭力矩..........229
圖四十六 (a)  D11=0.30,D12=0.15之複材葉片在α=3.5°、不同展弦鉸接位置時，整個葉片的升力......230
圖四十六 (b)  D11=0.30,D12=0.15之複材葉片在α=3.5°、不同展弦鉸接位置時，整個葉片的扭力矩....230
圖四十六 (c)  D11=0.30,D12=0.15之複材葉片在α=3.5°、不同襟翼角時，整個葉片的升力............231
圖四十六 (d)  D11=0.30,D12=0.15之複材葉片在α=3.5°、不同襟翼角時，整個葉片的扭力矩..........231
圖四十七 (a)  D11=0.30,D12=0.15之複材葉片在α=5.0°、不同展弦鉸接位置時，整個葉片的升力......232
圖四十七 (b)  D11=0.30,D12=0.15之複材葉片在α=5.0°、不同展弦鉸接位置時，整個葉片的扭力矩....232
圖四十七 (c)  D11=0.30,D12=0.15之複材葉片在α=5.0°、不同襟翼角時，整個葉片的升力............233
圖四十七 (d)  D11=0.30,D12=0.15之複材葉片在α=5.0°、不同襟翼角時，整個葉片的扭力矩..........233
圖四十八 (a)  D11=0.30,D12=0.15之複材葉片在α=3.5°、不同展弦大小時，整個葉片的升力..........234
圖四十八 (b)  D11=0.30,D12=0.15之複材葉片在α=3.5°、不同展弦大小時，整個葉片的扭力矩........234
圖四十八 (c)  D11=0.30,D12=0.15之複材葉片在α=3.5°、不同襟翼角時，整個葉片的升力............235
圖四十八 (d)  D11=0.30,D12=0.15之複材葉片在α=3.5°、不同襟翼角時，整個葉片的扭力矩..........235
圖四十九 (a)  D11=0.30,D12=0.15之複材葉片在α=5.0°、不同展弦大小時，整個葉片的升力..........236
圖四十九 (b)  D11=0.30,D12=0.15之複材葉片在α=5.0°、不同展弦大小時，整個葉片的扭力矩........236
圖四十九 (c)  D11=0.30,D12=0.15之複材葉片在α=5.0°、不同襟翼角時，整個葉片的升力............237
圖四十九 (d)  D11=0.30,D12=0.15之複材葉片在α=5.0°、不同襟翼角時，整個葉片的扭力矩..........237
圖五十 (a)  D11=0.30,D12=0.20之複材葉片在α=3.5°、不同翼剖面弦向位置時，整個葉片的升力......238
圖五十 (b)  D11=0.30,D12=0.20之複材葉片在α=3.5°、不同翼剖面弦向位置時，整個葉片的扭力矩....238
圖五十 (c)  D11=0.30,D12=0.20之複材葉片在α=3.5°、不同襟翼角時，整個葉片的升力..............239
圖五十 (d)  D11=0.30,D12=0.20之複材葉片在α=3.5°、不同襟翼角時，整個葉片的扭力矩............239
圖五十一 (a)  D11=0.30,D12=0.20之複材葉片在α=5.0°、不同翼剖面弦向位置時，整個葉片的升力....240
圖五十一 (b)  D11=0.30,D12=0.20之複材葉片在α=5.0°、不同翼剖面弦向位置時，整個葉片的扭力矩..240
圖五十一 (c)  D11=0.30,D12=0.20之複材葉片在α=5.0°、不同襟翼角時，整個葉片的升力............241
圖五十一 (d)  D11=0.30,D12=0.20之複材葉片在α=5.0°、不同襟翼角時，整個葉片的扭力矩..........241
圖五十二 (a)  D11=0.30,D12=0.20之複材葉片在α=3.5°、不同展弦鉸接位置時，整個葉片的升力......242
圖五十二 (b)  D11=0.30,D12=0.20之複材葉片在α=3.5°、不同展弦鉸接位置時，整個葉片的扭力矩....242
圖五十二 (c)  D11=0.30,D12=0.20之複材葉片在α=3.5°、不同襟翼角時，整個葉片的升力............243
圖五十二 (d)  D11=0.30,D12=0.20之複材葉片在α=3.5°、不同襟翼角時，整個葉片的扭力矩..........243
圖五十三 (a)  D11=0.30,D12=0.20之複材葉片在α=5.0°、不同展弦鉸接位置時，整個葉片的升力......244
圖五十三 (b)  D11=0.30,D12=0.20之複材葉片在α=5.0°、不同展弦鉸接位置時，整個葉片的扭力矩....244
圖五十三 (c)  D11=0.30,D12=0.20之複材葉片在α=5.0°、不同襟翼角時，整個葉片的升力............245
圖五十三 (d)  D11=0.30,D12=0.20之複材葉片在α=5.0°、不同襟翼角時，整個葉片的扭力矩..........245
圖五十四 (a)  D11=0.30,D12=0.20之複材葉片在α=3.5°、不同展弦大小時，整個葉片的升力..........246
圖五十四 (b)  D11=0.30,D12=0.20之複材葉片在α=3.5°、不同展弦大小時，整個葉片的扭力矩........246
圖五十四 (c)  D11=0.30,D12=0.20之複材葉片在α=3.5°、不同襟翼角時，整個葉片的升力............247
圖五十四 (d)  D11=0.30,D12=0.20之複材葉片在α=3.5°、不同襟翼角時，整個葉片的扭力矩..........247
圖五十五 (a)  D11=0.30,D12=0.20之複材葉片在α=5.0°、不同展弦大小時，整個葉片的升力..........248
圖五十五 (b)  D11=0.30,D12=0.20之複材葉片在α=5.0°、不同展弦大小時，整個葉片的扭力矩........248
圖五十五 (c)  D11=0.30,D12=0.20之複材葉片在α=5.0°、不同襟翼角時，整個葉片的升力............249
圖五十五 (d)  D11=0.30,D12=0.20之複材葉片在α=5.0°、不同襟翼角時，整個葉片的扭力矩..........249

參考文獻
[1] McCloud, J. L. Ⅲ, " An Analytical study of a Multicyclic Controllable Twist Rotor," AHS Forum, 1975.
[2] McCloud, J. L. Ⅲ, and Weisbrich, A. L., " Wind-Tunnel Test Results of a Full-Scale Multicyclic Controllable Twist Rotor," AHS forum, 1978.
[3] McHugh, F. J., and Shaw, J., Jr, " Benefits of Higher-Harmonic Blade Pitch : Vibration Reduction, Blade Load Reduction, and Performance Improvement," American Helicopter Society Mideast Region Symposium on Rotor Technology, Essington, PA, August 1976.
[4] Gates, A. A., Adrezin, R. S., and Wei, F.-S., " Servo-Flap Rotor Power Improvement Though Wind Tunnel Testing," AHS. 55th Annual Forum, Montreal, Quebec, Canada, May 25-27, 1999.
[5] Gates, A. A., Adrezin, R. S., and Wei, F.-S., " Thin Body Shape Optimization Technique Using Wind Tunnel Test Data," AIAA paper, AIAA-99-0836, 1999.
[6] Wei, F.-S., " Advanced Servo-Flap Rotor Using Variable Blade Index Angle Control," 38th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, & Materials, Kissimmee, Florida, April 7-10, 1997.
[7] Wei, F.-S., and Tomashofski, C. A., " HHC and IBC For a Servo-Flap Controlled Main Rotor," the 7th International Workshop on Dynamics and Aeroelastic stability modeling of Rotor craft Systems, St. Louis, Missouri, October 14-16, 1997.
[8] Fulton, M. V., and Ormiston, R. A., " Hover Testing of a Small-Scale Rotor with On-Blade Elevons," Journal of the AHS, pp. 96-106, April, 2001.
[9] Nikhil, A. K., and Chopra, I., " Open-Loop Hover Testing of a Smart Rotor Model," AIAA Journal, Vol. 40, No. 8, pp. 1495-1502, August 2002.
[10] Kurdila, A. J., Li, J., Strganac, T., and Webb, G., " Nonlinear Control Methodologies for Hysteresis in PZT Actuated On-Blade Elevons," Journal of AerospaceEngineering, pp. 167-176, October 2003.
[11] Greenberg, J. M., " Airfoil in Sinusoidal Motion in a Pulsating Stream," NACA TN 1326, June 1947.
[12] Garrick, I. E., " Propulsion of Flapping and Oscillating Airfoil," NACA TR567, May 1936.
[13] Abourahma, Ahmed Ali Mohammed Hassan, " The Analysis of Wake-Induced Unsteady Aerodynamics as Related to Higher Harmonic Control," Ph.D. Dissertation, Naval Postgraduate School, December 1993.
[14] Peters, D. A., and Johnson, M. J., " Finite-State Airloads For Deformable Airfoils on Fixed and Rotating Wings," AD-Vol. 44, Aeroelasticity and Fluid Structure Interaction Problems, ASME, pp. 1~28, 1994.
[15] Dowell, E. H., and Tang, D., " Nonlinear Aeroelasticity and Unsteady Aerodynamics," AIAA Journal, Vol. 40, No. 9, pp. 1697-1707, September 2002.
[16] Hodges, D. H., " Nonlinear Equations of Motion for Cantilever Rotor Blades in Hover with Pitch-Link Flexibility, Twist, Precone, Droop, Sweep, Torque Offset and Blade Root Offset," NASA TM X-73, 112, May 1976.
[17] Hodges, D. H. and Dowell, E. H., " Nonlinear Equations of Motion for the Elastic Bending and Torsion of Twisted Nonuniform Rotor Blades," NASA TN D-7818, December 1974.
[18] Hodges, D. H., Crespo da Silva, M. R. M., and Peters, D. A., " Nonlinear Effects in the Static and Dynamic Behavior of Beams and Rotor Blades ," Vertica, Vol. 12, pp. 243-256, 1988.
[19] Hodges, D. H., Atilgan, A. R., Cesnik, C. E. S., and Fulton, M. V., " On A Simplified Strain Energy Function for Geometrically Nonlinear Behavior of Anisotropic Beams," Composite Engineering Vol. 2, pp. 513-526, 1992.
[20] Pai, P. F., and Nayfeh, A. H., " Three-Dimensional Nonlinear Vibrations of Composite Beams - I. Equations of Motion," Nonlinear Dynamics 1, Kluwer Academic Publishers, pp. 477-502, 1990. 
[21] Pai, P. F., and Nayfeh, A. H., " Three-Dimensional Nonlinear Vibrations of Composite Beams - II. Flapwise Excitations," Nonlinear Dynamics 2, Kluwer Academic Publishers, pp. 1-34, 1991.
[22] Pai, P. F., " Nonlinear Flexural-Flexural-Torsional Dynamics of Metallic and Composite Beams," Ph. D. Thesis, Dept. of Engineering Science and Mechanics, Virginia Polytechnic Institute and State University, Blacksburg, VA, April 1990.
[23] Pai, P. F., and Nayfeh, A. H., " A Nonlinear Composite Beam Theory," Nonlinear Dynamics 3,  Kluwer Academic Publishers, pp. 273-303, 1992
[24] Pai, P. F., and Palazotto, A. N., " Large-Deformation Analysis of Flexible Beams," International Journal of Solids and Structures, Vol. 33, No. 9, pp. 1335-1353, 1996.
[25] Pai, P. F., Anderson, T. J., and Weather, E. A., " Large-Deformation Tests and Total-Lagrangian Finite-Element Analysis of Flexible Beams," International Journal of Solids and Structures, Vol. 37, No. 21, pp. 2951-2980, 2000.
[26] Crespo da Silva, M. R. M., " A Comprehensive Analysis of The Dynamics of A Helicopter Rotor Blade," International Journal of Solids and Structures, Vol. 35, No. 7-8, pp. 619-635, 1998.
[27] de Andrade, D., and Peters, D. A., " Correlation of Experimental Flap-Lag-Torsion Damping - A Case Study," Mathematical and Computer Modeling - Rotorcraft Modeling, Part II, Pergamon, pp. 135-158, 1994.
[28] Donizeti de Andrade, " Application of Finite-State Inflow to Flap-Lag-Torsion Damping in Hover," Ph. D. Thesis, School of Aerospace Engineering, Georgia Institute of Technology, May 1992.
[29] Tang, D., and Dowell, E.H., " Damping Prediction for a Stalled Rotor in Flap-Lag with Experimental Correlation," Journal of the American Helicopter Society, Vol. 40, No. 4, pp. 79-89, October 1995.
[30] Tang, D., and Dowell, E. H., " Damping Prediction for Hingeless Rotor Aeroelastic Stability with Experimental Correlation," Journal of Aircraft, Vol. 33, No. 6, pp. 1071-1078, November-December 1996.
[31] Tang, D., and Dowell, E. H., " Nonlinear Rotor Aeroelastic Analysis with Stall and Advanced Wake Dynamics," Journal of Aircraft, Vol. 34, No. 5, pp. 679-687, September-October 1997.
[32] Shang, X., Hodges, D. H., and Peters, D. A., " Aeroelastic Stability of Composite Hingeless Rotors in Hover with Finite-State Unsteady Aerodynamics," Journal of the American Helicopter Society, Vol. 44, No. 3, pp. 206-221, July 1999.
[33] Pai, P. F., and Nayfeh, A. H., " Fully Nonlinear Theory of Curved and Twisted Composite Rotor Blades According for Warpings and Three-Dimensional Stress Effects," International Journal of Solids and Structures, Vol. 31, No. 9, pp. 1309-1340, 1994.
[34] 王怡仁，彭暉元，唐鵬駿，" 三維尾流及非線性旋翼葉片耦合模式之葉片弦向阻尼分析 "，中華民國航空太空學會第四十一屆學術研討會，桃園 龍潭 中正理工學院，民國八十八年十二月。
[35] 王怡仁，唐鵬駿，" 三維尾流及旋翼葉片拍撲-延遲-扭轉耦合模式之建立及應用- I. 均質穩態之模擬 "，中華民國航空太空學會第四十屆學術研討會，台中 逢甲大學，民國八十七年十二月。
[36] 李明峰，" 動態失速下主旋翼懸空及前飛時的平衡態分析 "， 淡江大學航太系碩士班碩士論文，民國八十九年六月。
[37] Jeon, S.-M., and Lee, I., " Aeroelastic Analysis of a Hingeless Rotor Blade in Forward Flight," AIAA Journal, Vol. 38, No. 5, pp. 843-850, May 2000.
[38] Tang, D., and Dowell, E. H., " Nonlinear Rotor Aeroelastic Analysis with Stall and Advanced Wake Dynamics," Journal of Aircraft, Vol. 34, No. 5, pp. 679-687, September-October 1997.
[39] Dessi, D., Mastroddi, F. and Morino, L., " Limit-Cycle Stability Reversal near A Hopf Bifurcation with Aeroelastic Applications," Journal of sound and Vibration, Vol. 256, No. 2, pp. 347-365, 2002.
[40] Millot, T. and Friedmann, P., " Vibration Reduction in Helicopter Rotors Using an Actively Controlled Partial Span Trailing Edge Flap Located on the Blades." Technical Report CR 4611, NASA, June 1994.
[41] Milgram, J., Chopra, I., and Straub, F., " Rotors with Trailing Edge Flaps: Analysis and Comparison with Experiment Data," Journal of the American Helicopter Society, Vol. 43, (4), October 1998.
[42] Koratkar, N. A. and Chopra, I., " Wind Tunnel Testing of a Mach-Scaled Rotor Model with Trailing-Edge Flaps". In Proceedings of the 56th Annual Forum of the American Helicopter Society, Virginia Beach, May 2000.
[43] Shen, J. and Chopra, I., " Aeroelastic Modeling of Trailing-Edge Flaps with Smart Material Actuators," In Proceedings of the 41st AIAA/ASME/ASCE/ AHS/ASC structure, structural dynamics, and materials conference, AIAA-2000-1622, page 14, Atlanta, GA. April 3-6, 2000.
[44] Shen, J. and Chopra, I., " Aeroelastic Stability of Smart Trailing-Edge Flap Helicopter Rotors," In Proceedings of the 42nd AIAA/ASME /ASCE/AHS/ASC structure, structural dynamics, and materials conference, AIAA-2001-1675, page 11, Seattle, WA, April 16-19, 2001.
[45] Shen, J. and Chopra, I., " Actuation Requirements for a Swashplateless Helicopter Control System With Trailing-Edge Flaps,". In Procedings of the 43rd AIAA/ASME /ASCE/AHS structures, structural dynamics, and materials conference and 10th AIAA/ASME/AHS adaptive structures conference, number AIAA-2002-1444, page 11, Denver, Colorado, April 2002. 
[46] Straub, F. K. and Charles, B. D., " Aeroelastic Analysis of Rotors with Trailing Edge Flaps Using Comprehensive Codes," Journal of the American Helicopter Society, Vol. 46, (3): 192-199, July 2001.
[47] Shen, J. and Chopra, I., " A Parametric Design Study for a Swashplateless Helicopter Rotor with Trailing-Edge Flaps," presented at the American Helicopter Society 58th Annual Forum. Montrel, Canada, Jun 11-13, 2002. 
[48] 林宏益，" 複材葉片之旋翼後緣襟翼分析 "， 淡江大學航空太空工程研究所碩士論文，民國九十三年六月。
[49] 蔡耀昌，" 旋轉複合材料梁之振動模態分析 "， 淡江大學航空太空工程研究所碩士論文，民國九十二年六月。