崔海軍

（揚州工業(yè)職業(yè)技術(shù)學(xué)院，江蘇 揚州225127）

0 引言

受壓桿件的穩(wěn)定性計算是工程結(jié)構(gòu)中常見的分析問題［1］，如橋梁結(jié)構(gòu)中的桁架和高橋墩、鋼結(jié)構(gòu)工程中的斜梁、立柱［2］和支撐體系［3］、巖土工程中的鉆桿［4］以及工程建設(shè)中的腳手架桿件［5］。材料力學(xué)臨界力歐拉公式是在中心受壓直桿理想壓桿模型［6-7］、剪切變形忽略假定下推導(dǎo)的［8］。但由于桿件的初始曲率、豎向作用的初始偏心以及水平作用，必然導(dǎo)致壓桿的剪切變形［8］。

在以往文獻(xiàn)中，壓桿的臨界力按不同的約束條件，從壓桿的撓曲線近似微分方程求解，方法繁瑣，也不能全部推導(dǎo)證明［9］。文獻(xiàn)［10-11］利用彎矩微分方程和力的邊界條件推導(dǎo)了理想壓桿臨界力的歐拉公式，形成統(tǒng)一公式，簡化了推導(dǎo)過程［10-11］，但壓桿的穩(wěn)定問題的實質(zhì)是整體變形效應(yīng)問題，用力的邊界條件求解容易出現(xiàn)不確定和解的不一致問題［12］。文獻(xiàn)［9］以一端固定，一端鉸支的細(xì)長壓桿撓曲線微分方程推導(dǎo)了理想壓桿臨界力的歐拉公式，體現(xiàn)了桿的整體變形效應(yīng)［9］，但壓桿的撓曲位移僅為彎曲位移，沒有考慮剪切變形引起的位移，且壓桿失穩(wěn)模型上端鉸支，與不同約束條件下壓桿實際屈曲狀態(tài)不完全吻合。與文獻(xiàn)［9-11］不同，本文以細(xì)長壓桿在彎曲變形和剪切變形下平衡條件建立其撓曲線微分方程，推導(dǎo)不同約束邊界條件下考慮剪切變形影響下壓桿臨界力公式。

1 基本算式

如圖1所示，一發(fā)生剪切變形的細(xì)長受壓桿，下端固定，上端在平面內(nèi)自由。細(xì)長桿件受壓力F，同時在自由端附加一水平作用力FR，代替實際桿件在桿件上端可能受到的水平支反力，以及附加一外力偶MR，代替實際桿件在桿件上端可能受到的約束外力偶。壓桿彈性模量為E，橫截面慣性矩為I，剪切模量為G，橫截面積為A，自由端在平面內(nèi)微彎形態(tài)下平衡時最大水平位移為δ，建立圖示坐標(biāo)系，則細(xì)長壓桿在垂直于軸線的總位移W（x）可表示為：

式中，wf（x）為彎曲變形引起的位移，ws（x）為剪切變形引起的位移。取距離細(xì)長壓桿下端x的截面，該截面的轉(zhuǎn)角為 φ（x），則 φ（x）為：

由材料力學(xué)知：

對于剪力Q（x），在文獻(xiàn)［12］按下式方法一選?。?/p>

在文獻(xiàn)［13-14］按下式方法二選取：

但是，在撓度變形問題中，sin［φ（x）］＝φ（x）偏差較大，一般不再成立。所以，方法二剪力取值既不適用于小撓度變形問題，也不適用于大撓度變形問題，僅適用于文獻(xiàn)［13-14］中所闡述的橡膠材料的臨界力計算。本文按照方法一進(jìn)行考慮剪切變形狀態(tài)下的壓桿臨界力分析。

對于方法一：

式中，V為截面形狀系數(shù)，矩形截面取6/5，圓形截面取10/9。對于V的取值，文獻(xiàn)［16］尚有異議，因此，也有文獻(xiàn)取V為1，即忽略截面形狀系數(shù)。

對式（1）兩邊進(jìn)行求導(dǎo)，并將式（2）和式（8）代入得：

由式（9）得：

對式（10）兩邊進(jìn)行求導(dǎo)得：

將式（3）和式（4）代入式（11）得：

解此微分方程，通解為：

則w（x）的一階導(dǎo)數(shù)為

由式（9）及式（15）得：

圖1 壓桿計算模型

圖2 兩端鉸支壓桿失穩(wěn)模型

2 幾種桿端約束形式下考慮剪切變形的細(xì)長壓桿的臨界力

2.1 兩端鉸支的壓桿

如圖2所示，對于兩端鉸支的壓桿，δ＝0，MR＝0；由于壓桿上端水平支反力FR與下端水平支反力要平衡，構(gòu)成力偶，但桿件不受外力偶作用，因此FR＝0。位移邊界條件為：x＝0，w＝0；x＝l，w＝0。將邊界條件 x＝0，w＝0代入式（14）得：

B＝0

將邊界條件 x＝l，w＝0代入式（14）得：

A sin kl＝0

因為B＝0，A與B不可能同時為0，即A不等于0，則

sin kl＝0

由此求得 kl＝nπ，n＝0，1，2…

根據(jù)工程實際問題，取n＝1，將k代入式（13）得壓桿臨界壓力

式中

在式（17）中，當(dāng)剪切剛度＝∞時

即為不考慮剪切變形時的歐拉公式。

2.2 一端固定，一端自由的壓桿

如圖3所示，對于一端固定，一端自由的細(xì)長壓桿，自由端不能提供水平支反力及約束力偶，因此FR＝0，MR＝0。位移邊界條件為：x＝0，w＝0，φ＝0；x＝l，w＝δ。將邊界條件 x＝0，φ＝0代入式（16）得：

故A＝0

將FR＝0及位移邊界條件x＝0，w＝0代入式（14）得：

B＝-δ

將 x＝l，w＝δ代入式（14）得：

δ（1-cos kl）＝δ

則能使細(xì)長壓桿撓曲線成立的條件是：

cos kl＝0

從而得到

kl＝nπ/2，n＝1，3，5…

其最小解為n＝1，將k代入式（13）得壓桿臨界壓力

式中

在式（20）中，當(dāng)剪切剛度＝∞時

即為不考慮剪切變形時的歐拉公式。

2.3 兩端固定的壓桿

如圖4所示，對于兩端固定的壓桿，δ＝0；由于兩端固定壓桿邊界約束、所受外力、變形對壓桿中點上下對稱，壓桿兩端所受約束外力偶相互平衡，而壓桿上端水平支反力FR與下端水平支反力要平衡，構(gòu)成力偶，但桿件除固定端約束外力偶作用，不受其它外力偶作用，因此FR＝0（也可以根據(jù)對稱條件，水平支反力對壓桿中點反對稱，直接判斷水平支反力 FR＝0）。位移邊界條件為：x＝0，w＝0，φ＝0；x＝l，w＝0，φ＝0。

圖3 一端固定、一端自由壓桿失穩(wěn)模型

圖4 兩端固定壓桿失穩(wěn)模型

將位移邊界條件x＝0，φ＝0代入式（16）得

故A＝0

將位移邊界條件x＝0，w＝0代入式（14）得：

將位移邊界條件 x＝l，w＝0，φ＝0代入式（14）、式（16）得：

由此求得 kl＝2nπ，n＝0，1，2…

根據(jù)工程實際問題，取n＝1，將k代入式（13）得壓桿臨界壓力

式中

在式（20）中，當(dāng)剪切剛度＝∞時

即為不考慮剪切變形時的歐拉公式。

2.4 一端固定、一端鉸支的壓桿

如圖5所示，對于一端固定、另一端鉸支的壓桿，δ＝0；上端鉸支座端不能提供約束力偶，因此MR＝0；下端固定端有約束外力偶，這就要求壓桿上端鉸支座提供水平支反力FR，以平衡固定端約束外力偶。位移邊界條件為：x＝0，w＝0，φ＝0；x＝l，w＝0。

將位移邊界條件x＝0，w＝0代入式（14）得：

由壓桿受力特征分析知，壓桿下端剪力為-FR，則由式（5）、式（15）可得：

將位移邊界條件x＝l，w＝0代入式（14）得：

A sin kl+B cos kl＝0

將A、B代入上式得

tan kl＝kl

上式為超越方程，利用三角函數(shù)表或采用圖解法求解，求得最小正數(shù)值解

kl＝4.49

將k代入式（13）得壓桿臨界壓力

式中

在式（26）中，當(dāng)剪切剛度＝∞時

即為不考慮剪切變形時的歐拉公式。

圖5 一端固定、一端鉸支壓桿失穩(wěn)模型

圖6 一端固定、一端定向支撐壓桿失穩(wěn)模型

2.5 一端固定、一端定向支撐的壓桿

如圖6所示，一端固定、一端定向支撐的壓桿，即為兩端固定但可沿桿件橫向相對移動，因此FR＝0。根據(jù)桿件受力及變形的對稱性，對稱截面w（x/2）＝δ/2，M（x/2）＝0，則由式（4）得：

對于一端固定、一端定向支撐的壓桿，位移邊界條件為 x＝0，w＝0，φ＝0；x＝l，w＝δ，φ＝0。

將位移邊界條件 x＝0，w＝0及式（29）代入式（14）得

將位移邊界條件x＝0，φ＝0及式（28）代入式（16）得 A＝0

將位移邊界條件x＝l，w＝δ代入式（14）得

coskl＝-1

從而得到

其最小解為n＝0，將k代入式（13）得壓桿臨界壓力

式中

在式（30）中，當(dāng)剪切剛度＝∞時

即為不考慮剪切變形時的歐拉公式。

3 考慮剪切變形狀態(tài)下壓桿臨界力統(tǒng)一公式

由2.1～2.5五種不同約束條件下考慮桿件剪切變形狀態(tài)下壓桿臨界力分析，其壓桿臨界力公式可以統(tǒng)一寫成

Fcr表達(dá)式與材料力學(xué)歐拉公式一致，為長度因素，兩端鉸支，μ＝1，一端固定、另一端自由μ＝2，兩端固定μ＝0.5，一端固定、另一端鉸支μ＝0.7，一端固定、另一端定向支撐μ＝1

4 結(jié)論

（1）提出了一端固定、另一端平面內(nèi)自由的細(xì)長壓桿考慮剪切變形下壓桿失穩(wěn)計算模型。

（2）對計算模型在剪切變形和彎曲變形平衡條件下推導(dǎo)細(xì)長壓桿統(tǒng)一的撓曲線方程w（x）和轉(zhuǎn)角方程φ（x），w（x）由彎曲位移和剪切位移兩部分組成，φ（x）為彎曲位移方程的一階導(dǎo)數(shù)。

（3）針對五種約束條件下的壓桿，利用統(tǒng)一的撓曲線方程w（x）和轉(zhuǎn)角方程φ（x）推導(dǎo)了考慮剪切變形狀態(tài)下壓桿穩(wěn)定的臨界力公式及統(tǒng)一公式。