利用奇異函數法對輕型卡車縱梁彎曲剛度分析

左印波，孫江平，蔡欣

（安徽江淮汽車集團股份有限公司，安徽 合肥 230601）

利用奇異函數法對輕型卡車縱梁彎曲剛度分析

左印波，孫江平，蔡欣

（安徽江淮汽車集團股份有限公司，安徽 合肥 230601）

基于奇異函數法，并結合單位載荷法，推導出在一般受力狀態下的輕型卡車縱梁撓曲通用方程，用于求解車架縱梁最大撓度，進行彎曲剛度校核。

奇異函數；縱梁；撓曲方程；單位載荷；彎曲剛度

CLC NO.: U467.3 Document Code: A Article ID: 1671-7988 (2017)05-54-03

引言

輕型卡車的車架縱梁是車輛主要承載元件，承受安裝在其上面的各種裝置、人員和貨物的重量以及來自懸架系統的反作用力，因此須對縱梁進行彎曲剛度分析。

傳統的車架縱梁撓度的計算方法是在車架縱梁不同荷載、截面變參數等處分段，分別建立撓度方程，再聯立相鄰段的位移邊界條件來求解。當所受載荷復雜或縱梁為變截面、彎縱梁時，這種方法十分煩瑣整個過程煩瑣，邊界條件多。而利用奇異函數來求解，可以建立計算方程進行統一表達，從而能夠求解上得到很大的簡化，并且利于編程，利用計算機進行求解。

在進入20世紀90年代后，我國一些學者把奇異函數用于力學問題的求解[1-3]，但在汽車車架縱梁上的分析上的應用卻極少，有學者利用奇異函數僅對特定的變截面直縱梁強度分析進行了初步嘗試[4]。

本文基于奇異函數法，并結合單位載荷法，針對輕型載貨汽車常用的直縱梁、彎縱梁、變截面梁等車架縱梁，推導出統一撓曲方程，可精確計算出各點撓度及最大撓度，進行縱梁彎曲剛度分析。

1、奇異函數的定義

奇異函數又稱麥考利函數，是一個由多項奇異函數構成的函數族，由下式給出[1]：

當n≥0時，

當n＜0時，

當n=0時稱為單位階躍函數，n=-l時稱為δ函數，這里n為整數，它的基本微分為：

奇異函數的積分為：

2、利用奇異函數表達常見載荷及剛度

車架縱梁常見的受力情況及剛度變化，可用如下奇異函數表達[5]。

（1）作用在x=ai處的集中力偶Mi可表達為：（2）作用在x=bj處的集中力Pj可表達為：

（3）作用于（ck,dk）處的均布線載荷qk可表達為：

抗彎剛度EI沿整個構件長度方向發生階形變化,截面抗彎剛度EI發生t次，另t次的抗彎剛度的倒數(即柔度)的改變量βt稱為變柔度系數[6-7]，則有：

由式（8）可以寫出受彎構件的變參數（抗彎剛度EI發生改變）對構件的影響[5]：

3、基于奇異函數的彎曲剛度分析

本文以計算難度較大的變截面彎縱梁為基礎建立通用撓曲方程。

將輕型卡車前后懸架與車架縱梁相互作用力，簡化為支撐縱梁的支座，簡化后縱梁的受力情況如下：

圖1 變截面彎縱梁受力示意圖

其中：

A（0，0）—車架縱梁最前端；

B（xB，yB）—車架縱梁前拐點；

C（xc，yc）—車架縱梁后拐點；

l—車架縱梁沿x軸方向長度；

h—為車架上表面前后y方向高度落差；

D（l, h）—車架縱梁最后端；

R2—車架縱梁支座反力，坐標位置為(xR2,h)；

Pj、qk力的方向均為垂直于X軸；

?—車架縱梁折彎角度，00≤?＜900。

圖1表示變截面彎縱梁在x、y坐標下的受力情況，依據奇異函數的定義及性質，可以寫出作用在梁上沿x軸分布的等效載荷集度為：

利用截面法，發現在圖1所示的x、y坐標及載荷下（Pj、qk力的方向均為垂直于x軸），無論是在縱梁直斷還是折彎段所受的彎矩均與x軸坐標相關，對（11）式利用（4）公式，進行積分兩次有如下彎矩方程[4,7]：

如果梁上作用有L個集中力偶、M個集中載荷、N個均布載荷，則有如下彎矩方程：

其中，支座反力R1、R2可以通過力的平衡方程求出。

計算節點CPU支持AVX-512指令集，大幅度提升了處理器的浮點運算能力，GPU節點采用最新的NVIDIA Tesla V100型號GPU，并采用NVLink方式連接，單卡GPU雙精度浮點運算速度理論值可達7.8TFLOPS。計算網絡采用omni-path架構，安裝有intel、gcc等編譯器，作業管理采用slurm調度系統，集群管理采用計算中心自主開發的集群監控管理套件。

縱梁中心線實為雙折線，縱梁撓度計算需要沿縱梁中心線路徑積分，若用截面法分段處理，進行受力分析，建立撓度方程，邊界情況復雜，工作繁瑣。為簡化工作量，本文依據虛功原理，引入單位載荷法求解車架縱梁各點的撓度。

在縱梁上任意一點E點（xE,yE），施加垂直于縱梁軸線方向的單位力，依據單位載荷法[9]則有：

其中：

I(x)—車架縱梁慣性距變化函數。

單位載荷法計算時縱梁虛擬受力，存在如圖3、圖4兩種情況所示。

圖3 在AB或CD段施加單位載荷受力狀態

圖4 在BC段施加單位載荷受力狀態

利用奇異函數通用方程（13）及力的平衡方程，可以求出兩種狀態下的彎矩方程。

依據奇異函數可將上述兩種情況整合，則有：

常見的汽車縱梁慣性距變化往往為以下兩種慣性距變化的組合：

階梯梁：其變化函數表達見（8）～（9）公式；

線性變慣距梁：其變化函數見（10）公式。這里簡化起見令縱梁慣性距函數為I(x)，實際計算時分別帶入（8）～（10），利用奇異函數的定義及性質進行計算，不單獨展開。

依據微分原理，ds與dx存在以下微分關系：

將公式（14）進一步拓展，則有如下計算：

將（13）、（17）公式以及實際的I(x)函數組合公式帶入公式（19）中，即可求出任意一點E點（xE,yE）的撓度ωE。

當?=00時為特殊狀態下的直梁，積分路徑則為沿x軸方向積分，其計算公式簡化為：

當?=00時，積分路徑則為沿x軸方向積分，也可依據材料力學中彎矩與撓度的微分關系，對q(x)進行積分，可得：

其中：C、D為常數。

將（13）、（17）公式以及實際的I(x)函數組合公式帶入（21）公式中，再通過輕型卡車前后懸架與車架縱梁相互作用力簡化成支撐縱梁支座后的邊界條件。

以及奇異函數的積分公式（4）可形成一個奇異函數通用方程，求解車架縱梁上任意一點的撓度。

當?=00、EI為常量時，上述公式（21）可進一步簡化為：

其中，ω0、θ0可通過邊界條件（22）式求出。

4、結束語

本文利用奇異函數，結合單位載荷法，針對輕型載貨汽車常用的直縱梁、彎縱梁、變截面梁等車架縱梁，推導出通用的撓曲方程，克服了原來對變參數、變積分路徑問題要進行分段求解的不足，并且求解過程規整，有利于編制程序，求解出最大撓度ωmax進行彎曲剛度校核，極大解決了人工計算的復雜性。

[1] 王燮山.奇異函數在力學中的應用[M].科學出版社,1993.

[2] 劉念超.奇異函數解幾種超靜定結構[J]．長沙交通學院學報. 1998.14(2)：20—25.

[3] 荊振華.奇異函數在材料力學中的應川[J].遼寧省交通高等專科學校學報,1997,5(1)：52—58.

[4] 俞德津.利用奇異函數法對專用汽車縱梁的強度分析[J].專用汽車.2003(5)，6—8.

[5] 曾紀鵬.用奇異函數建立變參數連續梁撓度統一方程的新方法[J].四川建筑,2010,31(1)：128～130.

[6] 徐彬.高層建筑結構分析的奇異甬數方法[D].華南理工大學,2000.

[7] 徐彬.奇異函數建立高層結構分析模型的方法及應用[J].昆明理大學學報,2000,25(6)：106～109.

[8] 王燮山.用奇異函數求解某些變形梁的變形[J].力學與實踐,1984（4）：53～55.

[9] 李廉錕.結構力學[M].高等教育出版社，2004,97-99.

The bending stiffness analysis on the longitudinal beam of the light truck by using the singular function

Zuo Yingbo, Sun Jiangping, Cai Xin
( Anhui Jianghuai Automobile Co., Ltd, Anhui Hefei 230601）

By using the singular function and the unit load method,the deflection equation for the longitudinal beam of the light truck on the general stress state can be derived in this article.Then it can calculates the maximum deflection of the longitudinal beam by using the deflection equation for checking the bending stiffness.

Singular function; Longitudinal beam; Deflection curve equation; Unit load; Bending stiffness

U467.3

A

1671-7988 (2017)05-54-03

左印波（1986-），總布置工程師，就職于安徽江淮汽車集團股份有限公司，技術中心輕型商用車研究院，主要研究方向為整車總布置。

10.16638/j.cnki.1671-7988.2017.05.018