日字形薄壁桿件解析解與數值解對比分析

陳彪來

（甘肅交通職業技術學院，甘肅 蘭州 730070）

0 引言

在工程實際中薄壁桿件受扭狀態普遍存在，薄壁桿件雖然有部分特性類似于實體桿件，但對其進行力學性能分析時，橫截面的翹曲變形不再是次要因素而忽略不計。因此由薄壁桿件組成的薄壁結構，無論從薄壁桿件自身還是薄壁結構整體來講，其受力反應和變形都十分突出，必須對其分析。薄壁桿件力學的解析理論主要有不考慮截面外形輪廓線變形的烏曼斯基的閉口薄壁桿件的約束扭轉理論和符拉索夫的廣義坐標法，可以通過這些理論對薄壁桿件進行分析計算得到解析解，但其公式煩瑣，推導復雜。而有限元方法由于其強大的實用性，可以較精確地考慮各種復雜形狀和各種復雜的邊界條件，因而成為研究最多、應用最廣的數值方法。因而可以通過解析解與數值解的比較分析誤差產生的原因，驗證薄壁桿件在荷載作用下受翹曲效應的影響，并且得到有限元單元類型選擇的相關原則。

1 多室閉合薄壁桿件解析解計算

1.1 豎向荷載作用下的彎曲正應力

如圖1所示為薄壁桿件自由端截面圖，坐標原點在對稱中心處。因為其關于x軸和y軸對稱，所以它的形心、主極點、彎曲中心、剪切中心都與對稱中心重合。將豎向力往形心簡化可以得到剪力Qy和扭矩T，如圖2所示。

圖1 截面中線圖（單位：cm）

圖2 豎向力向形心簡化

根據材料力學：

則

1.2 軸力作用下的彎曲正應力

偏心力向形心簡化可以得到一個軸向拉力和兩個彎矩Mx和My，如圖3和圖4所示。

圖3 模型軸向拉力作用位置圖

圖4 軸向拉力向形心簡化后等效荷載作用位置圖

由上述彎矩作用的方向可以得到截面右上角點處最大的正應力：

如圖5所示，根據薄壁桿件理論首先在各室分別設置一個切口，將兩個切口分別取在1點的左側和右側，在截面的各節點處給以標號，剪力流q0i表示節點i處的靜定剪力流。應用公式計算時，由各切口處開始，依次計算各節點處的剪力流，顯然在切口處q0i為零。

1室各點處的靜定剪力流：

圖5 多室閉合薄壁桿件彎曲剪應

根據對稱性，2室各節點處截面靜定剪力流與1室對稱節點處剪力流相等。

為了求解超靜定剪力流qi（i=1，2，3，4，…，n），需依公式建立二元一次方程組。利用以下公式計算各室各板上的積分值：

可通過下面的方程組求出兩個附加剪力流q1、q2：

解得：

剪應力分布如圖6～圖8所示。

1.3 多室閉合薄壁桿件的自由扭轉分析

解得：

扭轉慣性矩：It=2×（1.2×10-3×0.03×2）=1.44×10-4（m4）

圖6 q0分布圖（單位：kN/m）

圖7 qi分布圖（單位：kN/m）

圖8 q分布圖（單位：kN/m）

利用公式（6），可求得各室自由扭轉剪力流：

薄壁桿件懸臂梁在自由端豎向力作用下產生扭矩 Mt=-40×0.05=-2（kN·m），可得：

故截面剪應力順時針轉,其分布如圖9所示。

圖9 τ分布圖（單位：MPa）

1.4 多室閉合薄壁桿件約束扭轉

根據截面中線圖可得極慣性矩：

由此得修正系數:

特征常數：

根據閉口薄壁桿件加載情況，得到約束扭轉微分方程：

根據加載情況，將豎向力與軸力分開求解微分方程。

1.4.1 豎向力作用下邊界條件

固定端：

自由端：

由以上邊界條件可得桿件雙力矩函數和彎扭力矩函數：

由雙力矩和彎扭力矩函數可計算固定端處的正應力：

翹曲正應力：

自由端：

1.4.2 軸向力作用下邊界條件

固定端：

自由端：

由以上邊界條件可得桿件雙力矩函數和彎扭力矩函數：

翹曲剪應力：

由雙力矩和彎扭力矩函數可計算固定端處的正應力：

固定端：

2 多室閉合薄壁桿件數值解計算

2.1 豎向力彎曲正應力數值解

利用ANSYS建模進行數值解計算，在豎向力作用下彎曲正應力數值解與解析解分布相同，其正應力分布關于水平軸呈反對稱分布，上邊緣為正，下邊緣為負。在角點處由于集中效應的影響，其正應力遠比解析解數值要大。角點處的正應力值為284.09 MPa，模型中線角點處的數值解為204 MPa。因為剪切滯后效應的影響，ANSYS數值解分布圖上下板的正應力不是均勻分布，在三塊腹板上也不再是線性分布。剪切滯后是由于剪切變形造成應力分布不均勻的現象，此時靠近節點處的腹板不再符合平截面假定，應力增大。

2.2 豎向彎曲剪應力數值解

模型橫截面上在各板的中間遠離角點的地方，剪應力分布比較均勻，且應力值不是很大。而在角點處存在局部應力集中的現象，這些點處的剪應力值會急劇增大，比其他各點處的值要大9～10倍。

2.3 豎向力作用整體分析

在偏心的豎向力作用下，結構會受到橫力彎曲、豎向剪力及扭矩的作用，所以截面上的應力包括彎曲正應力、彎曲剪應力、扭轉剪應力、翹曲剪應力、翹曲正應力。總正應力和總剪應力為各項的線性疊加。下面來分析薄壁桿件解析解中的翹曲效應對薄壁桿件的影響。

固定端右上角點最大正應力：

σ=σz+σω=177.78-0.913 2=176.867（MPa）

固定端右腹板最大剪應力：

τ=τw+τk+τω=1.54+1.667+0.467=3.674（MPa）

從上述解析解結果可以看出，翹曲正應力占總應力的0.913 2/176.867=0.52%，扭轉剪應力占總應力的2.134/3.674=58.1%。

2.3.1 正應力數值解比較

單元結點22830處的總正應力值是159.92 MPa,與解析解數值176.867 MPa比較，其誤差為9.6%。從2.1節的ANSYS計算結點22838數值解中,彎曲正應力為284.09 MPa,總正應力為283.84 MPa，所以固定端右上角點處的翹曲正應力是σω=283.84-284.09=-0.25（MPa），其翹曲正應力數值解與解析解相比，誤差為72.6%，誤差比較大。數值解的翹曲正應力0.25 MPa占總正應力283.84 MPa的比例約為0.09%，這一數值與解析解相比差別不大。

2.3.2 剪應力數值解比較

單元結點91325處的總剪應力值是3.378 MPa，從2.1節ANSYS計算的結點91325處數值解中,彎曲剪應力為1.120 4MPa，所以右腹板中間處的扭轉剪應力是 τn=3.378-1.102 4=2.275 6（MPa），其自由扭轉剪應力數值解與解析解十分接近，誤差為13.7%。數值解的扭轉剪應力2.275 6 MPa占總剪應力3.378 MPa的比例約為60.7%，這也與解析解的比值接近。

3 兩種求解方法比較及原因分析

（1）閉口薄壁桿件在此假定基礎上，運用薄壁桿件理論計算公式，計算約束扭轉剪應力及約束扭轉正應力。但是在用ANSYS計算閉口薄壁桿件時，運算過程中不會考慮到平截面假定。由上述分析可知，薄壁桿件解析解理論不考慮截面外形輪廓線的變形，這與實際情況是不符的。而用有限元軟件ANSYS進行數值計算反映了截面輪廓線的變形，所以兩者產生了較大的誤差。

（2）剪力滯后效應會使ANSYS數值解在界面邊緣的值不均勻，使得角點處的最大值遠大于解析解。如在第2節ANSYS數值解與解析解的比較中，由于橫力彎曲所引起的固定端正應力在上下板上應該是均勻分布的，但是從ANSYS數值解計算圖中可以看出在上下頂板上應力分布是不均勻的，特別是靠近節點板處，在腹板和頂板的交界處最大，隨著離開節點距離的增大而逐漸減小，因此正應力的橫向分布呈曲線形狀，使得頂、底板的中間小而兩邊大的分布狀態，使得彎曲正應力呈現不均勻分布的現象。

（3）局部集中效應會導致最大值出現在集中力作用點處，以及在腹板與頂、底板的交界處出現最大值。由于集中效應所造成的應力最大值與按薄壁桿件理論求得的解析解最大值有極大的誤差，所以比較力的作用點處的值是沒有意義的。根據圣維南原理，在遠離力的作用點處的應力值可以應用等效荷載去求解，其與理論求解結果是相吻合的。在力的作用截面上，由于薄壁桿件的厚度較小，在截面上按照圣維南原理受到集中荷載的影響較大，其用ANSYS計算的數值解與理論解析解的誤差也是相當大的。

4 結語

（1）用實體對閉口薄壁桿件建立有限元模型，可以很好地得出正應力及剪應力在截面上的分布情況，而且與用薄壁桿件理論得到的解析解分布相當吻合，如圖10和圖11所示。

圖10 數值解與解析解吻合

圖11 ANSYS計算模型

在角點處有集中效應、剪切滯后效應現象出現，所以在角點處的應力比數值解要大得多，應力分布也不是線性分布。由于殼單元在處理上與實體單元存在差異，尤其是在角點處得到的數值解與解析解十分接近，誤差不會超過10%。以豎向力的橫力彎曲為例，解析解在有角點處的數值為44.4 MPa，用殼單元分析得到的數值解是45.698 MPa,兩者誤差僅為2.9%,是相當吻合的。

（2）對于兩種荷載（豎向力、軸力）下的閉口薄壁桿件的翹曲效應不一樣，在軸力作用下自由端的翹曲效應比較大。解析解得到的翹曲正應力占到了66.2%,用實體分析得到的結果是57.2%，用殼單元分析得到的結果是50.4%，兩種翹曲正應力數值解所占比例都小于解析解結果。由于在集中荷載作用面上會有應力集中,圣維南原理不適用，所得到的總應力比實際情況偏大，并不是簡單的線性疊加。