理論發展背景
傳統彈性線理論無法描述細長桿進入塑性階段後的力學回響。2000年代Johnson與余同希通過引入塑性彎矩-曲率本構關係,將彈性控制方程擴展為包含塑性變形的大撓度方程組,標誌著塑性線模型的誕生。該突破使得細長桿在發生顯著塑性變形時的位形預測成為可能。
本構關係特性
塑性線模型的核心在於建立修正的彎矩-曲率關係:當桿件處於彈性階段時滿足$M=EI\kappa$的線性關係,進入塑性階段後轉變為分段函式形式。對於理想彈塑性材料,塑性區的曲率增量$\Delta\kappa$與彎矩$M$呈非線性關係。2023年研究表明,考慮線性強化效應時需引入硬化係數$m$,其取值範圍為0.02-0.15以反映不同材料的強化特性。
控制方程推導
在笛卡爾坐標系下,平衡微分方程可表示為:$$\frac{d^2\theta}{ds^2} + \frac{P}{EI}\sin\theta = 0$$其中$\theta$為截面轉角,$s$為弧長坐標。塑性階段需將彎曲剛度$EI$替換為與當前曲率相關的等效剛度$\overline{EI}(\kappa)$,這使得方程呈現強非線性特徵。數值求解時通常採用有限差分法配合牛頓-拉夫遜疊代,收斂步長控制在$10^{-4}$量級。
工程套用實例
在桁架結構分析中,塑性線模型可計算三桿桁架的極限載荷。當第2桿屈服後,其餘桿件通過彈性約束提供附加剛度,最終得到的塑性極限載荷$P_S$較彈性極限提高約4%。2024年兵器工程領域套用該模型分析身管殘餘應力,通過應變協調方程導出塑性半徑$r_p$與反向屈服半徑$r_y$的解析表達式。
幾何非線性處理
大變形分析需採用真應力-對數應變度量:這使得控制方程中幾何參數(如曲率半徑)與材料參數(如屈服強度)產生耦合效應。研究表明,當細長桿撓度超過長度20%時,幾何非線性引起的誤差可達傳統小變形假設的8倍以上。
失效判據建立
基於Drucker公設推導出塑性線模型的極限載荷曲面:$$f(\sigma_{eq}, \varepsilonp) = \sigma{eq} - Y_0 - H\varepsilonp = 0$$式中$\sigma{eq}$為等效應力,$Y_0$為初始屈服強度,$H$為塑性模量,$\varepsilon_p$為累積塑性應變。該判據可預測結構發生過量塑性變形或失穩塌陷的臨界狀態。

塑性線的圖片