彎曲 (力學)

**彎曲**（**bending**）也稱為**屈曲**（**flexure**），為材料力學的名詞．是指一形狀狹長的結構件固體，受到和其長軸垂直的外力時，固體變形的情形。

若結構件在某一方向長度很長，另外二方向的尺寸是該方向尺寸的1/10或更小，即滿足上述形狀狹長的定義，當其長度顯著的大於其寬度及厚度，此結構件會稱為梁。例如衣櫉中的橫桿會因為衣架上衣服的總重量而變形，就是梁受力彎曲的例子。另外方面，是指另外一種結構，其長寬是相近的數量級，但其厚度較長寬要小很多。例如一個大直徑、薄壁，長度和直徑相當的管子橫放，一側固定，上方乘載重量，是殼層受力彎曲的例子之一。

若沒有具體說明物體的形狀，「彎曲」可以指任何外形物體的彎曲。在工程應用上，有時會說明彎曲物體的形狀，例如「桿的彎曲」（bending of rods）、「梁的彎曲」（bending of beams）、（bending of plates）、殼層彎曲（bending of shells）等。
## 梁的準靜彎曲（1919年）

梁在受到側向力時，其內部會變形，而且會產生應力。準靜（quasi-static）彎曲是假設彎曲產生的形變以及應力不隨時間而變化。若考慮一水平梁，兩側固定，中間受到往下的力，梁的上半部會受到壓縮力，而梁的下半部會受到拉伸力。由於其側向力，會產生以下兩種內應力：
* 和側向力平行的剪應力，以及和受力方向垂直的互補剪應力
* 梁-{}-上半部受到的直接壓應力，以及梁下半部受到的直接拉伸應力。

後二個力會形成力偶或是，兩個力大小相等，方向相反。會抵抗梁受力彎曲時的變形特性。若一些簡化的假設成立，可以精確估測梁-{}-上的應力分佈。
### 歐拉-伯努力彎曲理論

在細長梁的歐拉-伯努力棟樑方程中，重要的假設是「平面截面在受力後仍維持平面」，換句話說，不考慮剪力對截面的影響（無剪切形變）。而且線性的應力應力分佈只適用於最大應力小於材料屈服強度情形。（條目會探討應力超過屈服強度的情形。）在材料屈服時，截面受到的最大應力（在離梁的最遠的位置）定義為。

考慮以下條件成立的梁： 
* 梁一開始是細長筆直的，就算有錐度，錐度也很小可以省略。
* 材料是各向同性（或是）、線性彈性，在任一截面都同質（不過沿著軸不一定要同質）。
* 只考慮小的形變。

此情形下，描述梁形變（w）可以近似為：
:\cfrac{\mathrm{d}^2 w(x)}{\mathrm{d} x^2}=\frac{M(x)}{E(x) I(x)}

其形變相對x的二次導數可以視為曲率，E是杨氏模量，I是截面的截面二次轴矩，M是梁的內部彎矩。

若細長梁沿著其軸向也是同質，而且其截面不會沿著軸而變化，而且也因為側向負荷q(x)而變形，可以得到：
:
EI~\cfrac{\mathrm{d}^4 w(x)}{\mathrm{d} x^4} = q(x)

這就是梁彎曲的歐拉-伯努力方程（歐拉-伯努力棟樑方程）。

若找到了梁形變位移的解，可以用下式計算彎矩（M）及剪力（Q）
:
M(x) = -EI~\cfrac{\mathrm{d}^2 w}{\mathrm{d} x^2} ~;~~ Q(x) = \cfrac{\mathrm{d}M}{\mathrm{d}x}.


簡單的梁彎曲常用歐拉-伯努力棟樑方程來分計。可以適用簡單彎曲理論的條件如下：
# 梁受到的力是。表示剪力為零，也沒有軸向負載以及沿著軸旋轉的負載。
# 材料為各向同性（或是），而且同質。
# 材料遵守胡克定律（在線性彈性範圍內，不會有塑性形變）。
# 梁一開始是筆直的，截面不會隨著軸而變化。
# 至少一個梁有在彎曲平面上的對稱軸。
# 梁的組成成份適當，使其最後是因為彎曲而失效，不是因為破裂、皺摺或挫曲而失效。
# 梁的截面在彎曲過程中仍維持平面。

在承受彎曲負荷時，沿著梁-{}-的軸上會產生壓縮力及伸張力。這些力會造成梁內部的應力。最大壓縮應力出現在梁的最上方，最大拉伸應力出現在梁的最下方。在中間的應力會呈線性變化，因此在梁內部，有一些點的彎曲應力為零。這些點的轨迹的軌跡稱為中性軸（neutral axis）。因為中性軸不受應力，中性軸附近的應力也比較小，若截面中的各個位置都是均質的，在中性軸附近的材料利用率較低，比較不經濟。寬緣樑（如工字梁）及桁架在應力較少的部份減少材料使用，因此可以減少總材料使用。

在純彎曲下，決定彎曲應力的公式如下：
:\sigma_x = \frac{M_z y}{I_z} = \frac{M_z}{W_z}
其中
*{\sigma_x}為彎曲應力
*M_z – 相對中性軸*z*的矩
*y – 相對中性軸的垂直距離
*I_z – 相對中性軸的截面二次轴矩
*W_z - 相對中性軸的阻力矩（Resistance Moment）。W_z = I_z / y
### 歐拉-伯努力彎曲理論的延伸
#### 塑性彎曲

方程式\sigma = \tfrac{M y}{I_x}只在最大應力（離中性軸最遠的位置）小於材料降伏應力的情形下。若負荷更大，則應力分布就會是非線性分析，延展性材料最後會進入「塑性鉸鏈」（plastic hinge）的情形，也就是在梁的各處應力大小都等於降伏應力，在中性軸的位置出現應力的不連續，從壓應力轉變成拉伸壓力。塑性鉸鏈狀態一般會用在鋼結構設計時的极限状态。
#### 複合或非對稱的彎曲
上述推導只用在截面對稱的條件。針對非對稱截面的均質梁，梁的最大彎曲應力如下：

: \sigma_x(y,z) = \frac {M_z~I_y - M_y~I_{yz}} {I_y~I_z - I_{yz}^2}y + \frac {M_y~I_z - M_z~I_{yz}} {I_y~I_z - I_{yz}^2}z

其中y,z是截面上一點的坐標M_y及M_z是相對位於幾何中心y軸和z軸的彎軸。I_y和I_z是相對y軸和z軸的截面二次軸矩，I_{yz}是面積乘積矩（Product moment of area）。用這些公式可以計算任意截面、任意彎矩下，在截面任意點的彎曲應力。其中M_y, M_z, I_y, I_z, I_{yz}不會隨截面上的不同位置而改變。

#### 大彎曲形變

針對梁的大幅形變，可以用以下延伸版本的公式來計算。其假設如下：

# 假設截面維持平面，在彎曲前和彎曲後的截面都維持平面
# 和截面法向量垂直的剪應力和正向應力沒有影響平行截面的正向應力。

若彎曲半徑\rho不到截面高度h的十倍時，需以大變形來考慮：

:\rho 

配合上述假設，大彎曲形變的應力為：

:
\sigma = \frac {F} {A} + \frac {M} {\rho A} + {\frac {M} y{\frac {\rho}{\rho +y}}


其中
:F為正向力
:A為截面積
:M為彎矩
:\rho為局部的彎曲半徑（目前部份的彎曲半徑）
:{\mathrm{d}x^4} - m\omega^2\hat{w} = 0
 
上式的通解為
:
   \hat{w} = A_1\cosh(\beta x) + A_2\sinh(\beta x) + A_3\cos(\beta x) + A_4\sin(\beta x)
 
其中A_1,A_2,A_3,A_4是常數，且

  \beta := \left(\cfrac{m}{EI}~\omega^2\right)^{1/4}


### 鐵木辛柯﹣瑞利理論

瑞利在1877年修正了歐拉-伯努力彎曲理論，多考慮了梁的截面的轉動慣量造成的效應。鐵木辛柯在1922年進一步的修正，多考慮了剪力的影響。鐵木辛柯﹣瑞利理論中允許梁的中表面之法向的剪力形變。

針對細長、各向同性、均質、截面積不變的梁，在受到動態側向負載下，鐵木辛柯﹣瑞利理論的動態彎曲方程為

: 
\begin{align}
& EI~\frac{\partial^4 w}{\partial x^4} + m~\frac{\partial^2 w}{\partial t^2} - \left(J + \frac{E I m}{k A G}\right)\frac{\partial^4 w}{\partial x^2~\partial t^2} + \frac{J m}{k A G}~\frac{\partial^4 w}{\partial t^4} \\[6pt]
= {} & q(x,t) + \frac{J}{k A G}~\frac{\partial^2 q}{\partial t^2} - \frac{EI}{k A G}~\frac{\partial^2 q}{\partial x^2}
\end{align}


其中J = \tfrac{mI}{A}是截面的极惯性矩，m = \rho A 是梁單位長度下的質量，\rho是梁的密度，A是截面積，G是剪力模量，k是剪力修正係數（shear correction factor）。針對帕松比\nu接近0.3的材料，剪力修正係數近似為
:
    \begin{align}
    k &= \frac{5 + 5\nu}{6 + 5\nu} \quad \text{rectangular cross-section}\\[6pt]
      &= \frac{6 + 12\nu + 6\nu^2}{7 + 12\nu + 4\nu^2} \quad \text{circular cross-section}
    \end{align}
 
#### 無外力振動
若梁沒有側向力時，鐵木辛柯﹣瑞利理論的簡諧振動方程如下
:
   EI~\cfrac{\mathrm{d}^4 \hat{w}}{\mathrm{d} x^4} + m\omega^2\left(\cfrac{J}{m} + \cfrac{E I}{k A G}\right)\cfrac{\mathrm{d}^2 \hat{w}}{\mathrm{d} x^2} + m\omega^2\left(\cfrac{\omega^2 J}{k A G}-1\right)~\hat{w} = 0
 
此式求解時，考慮w的各階導數需要是相同形式（才能互相抵消），因此其解為e^{kx}的形式。因此可以導出特徵方程式
:
   \alpha~k^4 + \beta~k^2 + \gamma = 0 ~;~~ \alpha := EI ~,~~ \beta := m\omega^2\left(\cfrac{J}{m} + \cfrac{E I}{k A G}\right) ~,~~ \gamma := m\omega^2\left(\cfrac{\omega^2 J}{k A G}-1\right)
 
此四次方程的解為
:
   k_1 = +\sqrt{z_+} ~,~~ k_2 = -\sqrt{z_+} ~,~~ k_3 = +\sqrt{z_-} ~,~~ k_4 = -\sqrt{z_-}
 
其中
:
   z_+ := \cfrac{-\beta + \sqrt{\beta^2 - 4\alpha\gamma}}{2\alpha} ~,~~
   z_-:= \cfrac{-\beta - \sqrt{\beta^2 - 4\alpha\gamma}}{2\alpha}
 
無外力振動下，鐵木辛柯﹣瑞利梁方程的通解為
:
   \hat{w} = A_1~e^{k_1 x} + A_2~e^{-k_1 x} + A_3~e^{k_3 x} + A_4~e^{-k_3 x}
 

## 準靜平板彎曲

梁的特點是其中一個方向的尺寸遠大於另外二個方向的尺寸。一結構若其中一個方向的尺寸遠小於另外二個方向的尺寸，則稱為平板。有許多理論要描述平板在受力下的形變以及應力分佈（板理论），其中有二種理論比較常用，分別是
* 克希荷夫–勒夫平板理論（Kirchhoff–Love theory of plates，也稱為經典平板理論）
* 明德林–賴斯納平板理論（Mindlin–Reissner plate theory，也稱為一階平板理論）
### 克希荷夫–勒夫平板理論
的假設是
* 和中表面垂直的直線在形變後仍然是直線
* 和中表面垂直的直線在形變後仍然和中表面直線垂直
* 在形變前後，平板的厚度不會變化。

上述的假設意味著
:
   \begin{align}
     u_\alpha(\mathbf{x}) & = - x_3~\frac{\partial w^0}{\partial x_\alpha}
         = - x_3~w^0_{,\alpha} ~;~~\alpha=1,2 \\
     u_3(\mathbf{x}) & = w^0(x_1, x_2)
   \end{align}

其中\mathbf{u}是板上一點的形變，而w^0是中表面上的位移。

應變和位移的關係如下
:
  \begin{align}
    \varepsilon_{\alpha\beta} & =
      - x_3~w^0_{,\alpha\beta} \\
    \varepsilon_{\alpha 3} & = 0 \\
    \varepsilon_{33} & = 0
  \end{align}


平衡方程式為
:
     M_{\alpha\beta,\alpha\beta} + q(x) = 0 ~;~~ M_{\alpha\beta} := \int_{-h}^h x_3~\sigma_{\alpha\beta}~dx_3

其中q(x)是和平板表面垂直的力

若以位移來表示，在沒有外力下，各向同性、線彈性平板的平衡方程為
:
   w^0_{,1111} + 2~w^0_{,1212} + w^0_{,2222} = 0
 
若以直接張量表示法，可以表示如下
:
   \nabla^2\nabla^2 w = 0
 
### 明德林–賴斯納平板理論
的假設類似克希荷夫–勒夫平板理論，和中表面垂直的直線在形變後仍然是直線，而且不會延展。但是和克希荷夫–勒夫平板理論不同的是：表面垂直的直線在形變後不一定仍和中表面垂直。平板的位移為
:
   \begin{align}
     u_\alpha(\mathbf{x}) & = - x_3~\varphi_\alpha ~;~~\alpha=1,2 \\
     u_3(\mathbf{x}) & = w^0(x_1, x_2)
   \end{align}

其中\varphi_\alpha是法向的旋轉量。

依照以上假設，應變和位移之間的關係是
:
  \begin{align}
    \varepsilon_{\alpha\beta} & =
      - x_3~\varphi_{\alpha,\beta} \\
    \varepsilon_{\alpha 3} & = \cfrac{1}{2}~\kappa\left(w^0_{,\alpha}- \varphi_\alpha\right) \\
    \varepsilon_{33} & = 0
  \end{align}

其中\kappa是剪力修正係數。

平衡方程為
:
   \begin{align}
     & M_{\alpha\beta,\beta}-Q_\alpha = 0 \\
     & Q_{\alpha,\alpha}+q = 0
   \end{align}
 
其中
:
   Q_\alpha := \kappa~\int_{-h}^h \sigma_{\alpha 3}~dx_3

## 平板的動態彎曲

### 克希荷夫平板的動態彎曲
平板的動態理論會決定在平板上震波的傳播，以及其駐波以及振動模態。克希荷夫平板動態彎曲的統御方程如下
:
     M_{\alpha\beta,\alpha\beta} - q(x,t) = J_1~\ddot{w}^0 - J_3~\ddot{w}^0_{,\alpha\alpha}

其中（假設板的密度為\rho = \rho(x)）
:
   J_1 := \int_{-h}^h \rho~dx_3  ~;~~
   J_3 := \int_{-h}^h x_3^2~\rho~dx_3

以及
:
   \ddot{w}^0 = \frac{\partial^2 w^0}{\partial t^2} ~;~~
   \ddot{w}^0_{,\alpha\beta} = \frac{\partial^2 \ddot{w}^0}{\partial x_\alpha\, \partial x_\beta} 


以下是一些圓形平板的振動模態。

File:Drum vibration mode01.gif|mode *k* = 0, *p* = 1
File:Drum vibration mode02.gif|mode *k* = 0, *p* = 2
File:Drum vibration mode12.gif|mode *k* = 1, *p* = 2

## 參見
* 剪切強度
* 振動學
* 
* 
* 
* 连续介质力学
* 
* 

* 
* 剪切强度
* 
* 振动
* 
* 

## 參考資料

## 外部連結
* Flexure formulae
* Beam flexure, stress formulae and calculators