李昕，沈斌（同濟(jì)大學(xué) 中德學(xué)院，上海 201804）

各向同性材料的強(qiáng)度校核在Matlab中的實(shí)現(xiàn)

李昕，沈斌
（同濟(jì)大學(xué) 中德學(xué)院，上海 201804）

在汽車設(shè)計(jì)生產(chǎn)中，為了能高效及時(shí)的分析車身材料的特性及其之間相互影響的關(guān)系，需要開發(fā)一款利用Matlab評(píng)估材料特性的軟件。論文介紹了由各向同性材料制成的具有不同截面的桿件在受到多個(gè)線性疊加外部靜載荷時(shí)強(qiáng)度校核的方法以及應(yīng)力集中現(xiàn)象對(duì)強(qiáng)度的影響；其次簡(jiǎn)述了該強(qiáng)度校核過程在Matlab中實(shí)現(xiàn)的算法思路以及主要的算法流程；最后對(duì)軟件計(jì)算結(jié)果和實(shí)際工程應(yīng)用中的計(jì)算結(jié)果進(jìn)行比較和誤差分析，保證計(jì)算結(jié)果的可靠性。該算法的實(shí)現(xiàn)保證了汽車在實(shí)現(xiàn)輕量化生產(chǎn)的同時(shí)并沒有損失其必要的安全性。

各向同性；強(qiáng)度校核；Matlab編程

0　引言

在汽車設(shè)計(jì)生產(chǎn)過程中，為了保證汽車安全性能和舒適度，需要對(duì)車身材料的性能，結(jié)構(gòu)及其之間相互影響的關(guān)系進(jìn)行研究分析。因此有必要開發(fā)一款利用計(jì)算機(jī)對(duì)材料性能以及各性能之間的關(guān)系進(jìn)行有效評(píng)估的軟件系統(tǒng)。這款軟件系統(tǒng)是基于Matlab的編程語(yǔ)言開發(fā)實(shí)現(xiàn)的，同時(shí)材料強(qiáng)度是評(píng)定材料性能時(shí)最重要的標(biāo)準(zhǔn)之一，在汽車車身制造中各向同性材料也占有重要的位置。因此本文著重研究了各向同性材料的強(qiáng)度校核在Matlab中的實(shí)現(xiàn)。

本文實(shí)現(xiàn)了利用Matlab對(duì)不同截面的桿件或梁在線性疊加的靜載荷下的強(qiáng)度計(jì)算，并分析了應(yīng)力集中現(xiàn)象對(duì)強(qiáng)度的影響，最后將Matlab計(jì)算出來(lái)的結(jié)果和實(shí)際計(jì)算結(jié)果進(jìn)行了比較和誤差分析。

1　理論基礎(chǔ)

1.1單一載荷下的應(yīng)力計(jì)算

強(qiáng)度是指某種材料或者由某種材料制成的構(gòu)件抵抗斷裂或永久變形的能力。桿件最常見的四種基本變形為：拉壓，剪切，彎曲，扭轉(zhuǎn)。桿件受拉壓和彎曲所產(chǎn)生的應(yīng)力為正應(yīng)力，受剪切和扭轉(zhuǎn)所產(chǎn)生的應(yīng)力為剪切應(yīng)力。

桿件受到拉壓應(yīng)力的大小等于其單位橫截面積上所受的軸向力的大小。受到切應(yīng)力大小等于其單位剪切面積上所受剪切力的大小。桿件上某位置的最大彎曲應(yīng)力為此位置引起形變的彎矩Mb與此截面的抗彎截面模量W的比值。桿上某位置的扭轉(zhuǎn)應(yīng)力等于該位置所受扭矩Mt與該截面的抗扭截面模量Wp的比值。表1中列出了桿件分別受到四種基本變形后產(chǎn)生的應(yīng)力計(jì)算公式。

表1　單一載荷下桿件應(yīng)力計(jì)算公式Tab.1 Stress calculation formula ofrods underthe single loads

1.2截面形狀及相關(guān)幾何量

根據(jù)表1中的公式，要計(jì)算桿件在外力作用下產(chǎn)生應(yīng)力的大小，不僅需要外部載荷，還要計(jì)算截面面積，抗彎和抗扭截面模量。這些根據(jù)桿件截面形狀及尺寸推導(dǎo)得出的量稱為與截面相關(guān)的幾何量。要確定這些幾何量的大小，就要先確定待研究的截面形狀。根據(jù)在汽車中常見桿件和梁的截面形狀分析，確定六種截面為主要研究對(duì)象，即橢圓，長(zhǎng)方形，I型，T型，Z型和U型截面。其他截面均可由此六種形狀推導(dǎo)得出?？箯澖孛婺Ａ縒和抗扭截面模量Wp是兩個(gè)分別用于計(jì)算截面上產(chǎn)生的最大彎曲應(yīng)力和扭轉(zhuǎn)應(yīng)力的幾何量。它們可以分別通過計(jì)算軸慣性矩I或極慣性矩Ip與在橫截面上離中性軸最遠(yuǎn)的點(diǎn)和中性軸間的距離amax的比值得出：

1.3載荷形式

桿件所受的彎矩是計(jì)算彎曲應(yīng)力時(shí)一個(gè)重要的物理量。它與梁的類型以及所受彎曲載荷的形式有關(guān)。本文中主要研究常見的簡(jiǎn)支梁和懸臂梁。對(duì)每種梁又分別研究了五種外部載荷類型，見表2。根據(jù)每種梁在不同載荷下的受力分析可以求出其關(guān)于梁長(zhǎng)x的彎矩方程，此彎矩方程的極值，即為梁所受的最大彎矩的位置及大小，由此可求出所產(chǎn)生的最大彎曲應(yīng)力。

表2　引起簡(jiǎn)支梁和懸臂梁發(fā)生彎曲變形的五種載荷形式Tab.2 Five load forms that cause the bending deformation of simply supported beam and cantilever

1.4胡克定律和強(qiáng)度理論

根據(jù)胡克定律，當(dāng)材料處于線彈性變形范圍內(nèi)時(shí)，其變形所產(chǎn)生的應(yīng)力σ，τ與應(yīng)變?chǔ)?，γ成線性比例關(guān)系，比例系數(shù)為彈性模量E或剪切模量G：

要校核桿件的強(qiáng)度還需合適的強(qiáng)度理論作為基礎(chǔ)，常用的強(qiáng)度理論有最大拉應(yīng)力，最大切應(yīng)力以及形狀改變比能理論。最大拉應(yīng)力理論認(rèn)為材料發(fā)生脆性斷裂的主要原因是最大拉應(yīng)力σ1達(dá)到極限值σb，它用于校核脆性材料發(fā)生斷裂破壞時(shí)的強(qiáng)度。其強(qiáng)度條件：

式中：S—安全系數(shù)；［σ］—許用拉應(yīng)力。最大切應(yīng)力理論用于校核塑性材料發(fā)生塑性流動(dòng)破壞或者滑移斷裂時(shí)的強(qiáng)度。理論上認(rèn)為無(wú)論構(gòu)件處于何種應(yīng)力狀態(tài)，構(gòu)件危險(xiǎn)點(diǎn)的最大切應(yīng)力τmax達(dá)到了材料在單向拉伸屈服時(shí)的極限剪應(yīng)力τ0，是發(fā)生塑性流動(dòng)的主要原因。由于：

用主應(yīng)力表達(dá)的最大切應(yīng)力理論的強(qiáng)度條件為：

式中：σ1，σ3—構(gòu)件危險(xiǎn)點(diǎn)的最大和最小主應(yīng)力；σ—材料屈服極限。形狀改變比能理論認(rèn)為：塑性材料發(fā)生屈服的主要原因，是構(gòu)件危險(xiǎn)點(diǎn)的最大形狀改變比能達(dá)到了一個(gè)極限值。該理論用主應(yīng)力表達(dá)的強(qiáng)度條件為：

式中：σ1，σ2，σ3—構(gòu)件危險(xiǎn)點(diǎn)的主應(yīng)力。將這三個(gè)強(qiáng)度理論分別用應(yīng)力轉(zhuǎn)換公式轉(zhuǎn)換為用拉壓應(yīng)力，彎曲應(yīng)力和扭轉(zhuǎn)應(yīng)力表達(dá)的平面應(yīng)力狀態(tài)下的強(qiáng)度條件公式如下：

式中：σzd—構(gòu)件危險(xiǎn)點(diǎn)的最大拉或壓應(yīng)力；σb—構(gòu)件危險(xiǎn)點(diǎn)的最大彎曲應(yīng)力；τt—構(gòu)件危險(xiǎn)點(diǎn)的最大扭轉(zhuǎn)應(yīng)力。通過計(jì)算桿件受到線性疊加的合成載荷下所產(chǎn)生最大拉壓應(yīng)力，彎曲應(yīng)力及扭轉(zhuǎn)應(yīng)力，判斷其可能的失效形式后通過相應(yīng)的強(qiáng)度理論，可以對(duì)桿件的強(qiáng)度利用公式（8）～（10）進(jìn)行校核。

1.5應(yīng)力集中現(xiàn)象

在實(shí)際工程中所使用的桿件往往是經(jīng)過開槽，挖孔，做軸肩后的桿件。這會(huì)使桿件的橫截面積尺寸發(fā)生局部突變。根據(jù)實(shí)驗(yàn)可得，在尺寸發(fā)生突變的橫截面應(yīng)力分布是不均勻的。從圖1可以看出，應(yīng)力在孔和凹槽的邊緣發(fā)生急劇升高，這種由于尺寸變化引起局部應(yīng)力增大的現(xiàn)象叫做應(yīng)力集中現(xiàn)象。

應(yīng)力集中系數(shù)αk可以反映出桿件在靜態(tài)外載荷作用下應(yīng)力集中的程度，它是應(yīng)力集中處的最大應(yīng)力 σmax與凈截面上的平均應(yīng)力σn的比值：

圖1　圓孔和凹槽處的應(yīng)力集中現(xiàn)象Fig.1 Stress concentration at the hole and the groove

應(yīng)力集中系數(shù)與構(gòu)件的材料無(wú)關(guān)，而與切槽的深度，孔的大小以及變截面的過渡圓弧大小有關(guān)。表1中所羅列的應(yīng)力計(jì)算公式均用于計(jì)算等直桿在靜載荷作用下凈截面上的平均應(yīng)力大小，在進(jìn)行強(qiáng)度校核時(shí)，還需要考慮到應(yīng)力集中現(xiàn)象，求得應(yīng)力集中處的最大應(yīng)力，再通過強(qiáng)度理論求得用于與許用應(yīng)力進(jìn)行比較的比較應(yīng)力。常見幾種構(gòu)件的應(yīng)力集中系數(shù)，可以在相關(guān)的工程手冊(cè)中查詢。

2　強(qiáng)度計(jì)算在Matlab中的實(shí)現(xiàn)

2.1計(jì)算總流程

利用Matlab實(shí)現(xiàn)強(qiáng)度計(jì)算時(shí)，首先要根據(jù)構(gòu)件所受外部載荷的種類及大小分別計(jì)算出最大拉壓，彎曲和扭轉(zhuǎn)應(yīng)力。但此應(yīng)力僅是構(gòu)件凈截面上所受的最大平均應(yīng)力。因此還必須要考慮到實(shí)際工程的應(yīng)力集中現(xiàn)象。各應(yīng)力分別乘以相應(yīng)的應(yīng)力集中系數(shù)，即可獲得應(yīng)力集中現(xiàn)象下的截面最大應(yīng)力。將計(jì)算得出的最大應(yīng)力帶入合適的強(qiáng)度條件公式中，即可得出比較應(yīng)力。如果比較應(yīng)力沒有超過材料的許用應(yīng)力，即可說(shuō)桿件在這種外部載荷下滿足強(qiáng)度條件，強(qiáng)度足夠。計(jì)算總流程見圖2。

2.2最大合成應(yīng)力的計(jì)算

為了求出構(gòu)件凈截面的最大平均應(yīng)力，要首先分別求出構(gòu)件危險(xiǎn)截面的所受外部載荷的大小。想要求出凈截面上最大拉壓應(yīng)力的大小，可根據(jù)胡克定律將可以引起拉壓應(yīng)力的外部載荷（拉/壓力）線性疊加，求出合成載荷，然后根據(jù)表1中的應(yīng)力計(jì)算公式計(jì)算構(gòu)件的合成應(yīng)力的大小。用同樣的方法對(duì)構(gòu)件所受的扭矩進(jìn)行疊加然后利用表1中的公式，也可以求得合成扭轉(zhuǎn)應(yīng)力的大小。

但要計(jì)算構(gòu)件所受的最大彎矩要先在程序中確認(rèn)梁的種類，其次要分析引起彎曲應(yīng)力的外部載荷。這些外部載荷可以根據(jù)胡克定律分解為單一的外部載荷，每種單一外部載荷都有對(duì)應(yīng)的關(guān)于梁長(zhǎng)x的撓度和彎矩方程，將這些單獨(dú)的彎矩方程線性疊加之后便得到了一個(gè)關(guān)于梁長(zhǎng)的總彎矩方程。利用 Matlab程序語(yǔ)言求出這個(gè)方程在梁長(zhǎng)范圍內(nèi)的極值，即為梁所受最大合成彎矩，然后利用彎曲應(yīng)力計(jì)算公式，可以求出凈截面在受多個(gè)線性疊加外載荷時(shí)產(chǎn)生的最大彎曲應(yīng)力的大小。計(jì)算合成彎曲應(yīng)力的算法流程如圖3所示。

圖2　在Matlab中實(shí)現(xiàn)強(qiáng)度計(jì)算的總流程圖Fig.2 Master Flowchart of the strength calculation with Matlab

圖3　計(jì)算合成彎曲應(yīng)力的算法流程圖Fig.3 Algorithmic Flowchart of the calculation of resultant bending stress

2.3應(yīng)力集中系數(shù)的計(jì)算

應(yīng)力集中系數(shù)一般可以根據(jù)構(gòu)件尺寸，應(yīng)力集中類型和受載荷形式在應(yīng)力集中系數(shù)曲線上讀出。但在Matlab中這種方式效率低下并且不易于代碼的實(shí)現(xiàn)。觀察常見的應(yīng)力集中系數(shù)曲線可以發(fā)現(xiàn)，這些曲線可以歸納為關(guān)于構(gòu)件尺寸和危險(xiǎn)截面尺寸的方程。利用這些方程，就可以在Matlab中計(jì)算出近似的應(yīng)力集中系數(shù)。

文獻(xiàn)［4］中給出了通過實(shí)驗(yàn)和數(shù)值分析方法進(jìn)行驗(yàn)證得出的應(yīng)力集中系數(shù)近似方程。通過近似方程計(jì)算得出的應(yīng)力集中系數(shù)和圖表讀取的應(yīng)力集中系數(shù)有大約10%的誤差，但這樣的誤差精度對(duì)于我們強(qiáng)度校核的影響是很小的，于是在Matlab實(shí)現(xiàn)過程中，使用了近似方程進(jìn)行應(yīng)力集中系數(shù)的計(jì)算。

表3是文獻(xiàn)［4］中給出的其中一個(gè)近似方程。它用于計(jì)算方形桿和圓形桿在具有切槽和軸肩時(shí)分別受拉壓、彎曲和扭轉(zhuǎn)時(shí)的應(yīng)力集中系數(shù)。

表3　圓形桿和對(duì)稱方形桿應(yīng)力集中系數(shù)計(jì)算近似方程Tab.3 Approximate equation for calculating the stress concentration factor of round and symmetrical square rods

3　計(jì)算結(jié)果驗(yàn)證及誤差分析

為保證Matlab計(jì)算結(jié)果的正確性，要確保每個(gè)子程序計(jì)算結(jié)果正確。在整個(gè)程序中，其他子程序均可通過簡(jiǎn)單的公式帶入來(lái)驗(yàn)證計(jì)算結(jié)果的準(zhǔn)確性，但是在計(jì)算最大彎矩時(shí)用到了 Matlab中的方程極值計(jì)算函數(shù)fminbnd（f（x），a，b），為了確定這個(gè)函數(shù)可以準(zhǔn)確的計(jì)算出合成彎矩方程的極值，我們進(jìn)行了如下的驗(yàn)證。

如圖4所示，一個(gè)懸臂梁同時(shí)受到了單力、均布和遞減載荷的作用發(fā)生彎曲，產(chǎn)生彎曲應(yīng)力?，F(xiàn)將其分解為分別單獨(dú)受到這三種載荷的三個(gè)相同的懸臂梁，可得到如下的三個(gè)單獨(dú)彎矩圖。

圖4　受到線性疊加載荷的懸臂梁線性分解為三個(gè)受單一載荷的懸臂梁Fig.4 Linear decomposition of cantilever by linear superposition load in three cantilevers by single loads

根據(jù)這三個(gè)彎矩圖可以看出，懸臂梁受到最大彎矩位置在梁首，即x=0處，最大彎矩的大小為分解后三個(gè)梁首所受彎矩之和，即：

將相同的數(shù)據(jù)帶入Matlab程序中，可得到：

xk=4.5085e-05

Bkmax=2.3333e+03

比較兩個(gè)計(jì)算結(jié)果可以看出，Matlab可以準(zhǔn)確計(jì)算出合成總彎矩方程的極值Bkmax。xk為極值出現(xiàn)的位置，由于Matlab程序本身的計(jì)算精度這個(gè)值產(chǎn)生了微小的誤差，可忽略不計(jì)。

其次在用Matlab求應(yīng)力集中系數(shù)時(shí)采用了近似方程，而非傳統(tǒng)的讀取應(yīng)力集中系數(shù)曲線，因此兩個(gè)結(jié)果間存在誤差。經(jīng)大量數(shù)據(jù)測(cè)試，利用近似方程求出的應(yīng)力集中系數(shù)普遍超過曲線中讀取結(jié)果的10%左右。這說(shuō)明，當(dāng)Matlab計(jì)算出的構(gòu)件強(qiáng)度滿足強(qiáng)度條件時(shí)，在實(shí)際工程中這個(gè)構(gòu)件一定滿足強(qiáng)度要求不會(huì)發(fā)生失效。因此，計(jì)算應(yīng)力集中系數(shù)而產(chǎn)生的誤差使強(qiáng)度得到了更為保守安全的校核結(jié)果，是可接受的誤差。

4　結(jié)論

本文利用胡克定律和強(qiáng)度理論，實(shí)現(xiàn)了利用Matlab編程計(jì)算具有不同截面的各向同性材料制成的桿件在靜態(tài)的復(fù)合外部載荷下的強(qiáng)度計(jì)算。對(duì)計(jì)算結(jié)果進(jìn)行誤差分，確保在適度誤差下桿件在實(shí)際工程中也不會(huì)發(fā)生失效，從而保證校核結(jié)果的可靠性。通過在Matlab中實(shí)現(xiàn)各向同性材料的強(qiáng)度計(jì)算，可以在汽車產(chǎn)品設(shè)計(jì)的早期階段就對(duì)相關(guān)構(gòu)件進(jìn)行強(qiáng)度校核，保證汽車在足夠的強(qiáng)度和安全性能下達(dá)到輕量化生產(chǎn)的目標(biāo)。

［1］宋子康，蔡文安.材料力學(xué)［M］.上海同濟(jì)大學(xué)出版社，1997.

［2］同濟(jì)大學(xué)航空航天與力學(xué)學(xué)院基礎(chǔ)力學(xué)教學(xué)研究部.材料力學(xué)［M］.上海同濟(jì)大學(xué)出版社，2011.

［3］Gross，D.；Hauger，W.；Schr?der，J.；Wall，W.A.：Technische Mechanik 1.Statik.10.Auflg.Berlin，Heidelberg：Springer-Verlag，2011.

［4］Haibach，E.：Betriebsfestigkeit-Verfahren und Daten zur Bauteilberechnung.3.Auflg.Berlin，Heidelberg：Springer-Verlag，2006.

Strength Calculation of Isotropic Materials with Matlab

LI Xin，SHEN Bin
（Sino-German Kollege，Tongji University Shanghai，Shanghai 201804，China）

In order to realize the efficient and fast analysis of the properties of various materials in carrosserie，it is of realistic significance to take use of Matlab to assess the material properties.This article describes the methods of strength calculation，when the rods with different cross-sections，which are made of isotropic materials，are subjected to multiple linear superposition of external static load.Besides，the effects of stress concentration are also discussed in this article.What's more，the main ideas of the algorithm，which are used to implement the strength calculation of isotropic materials with Matlab，are presented after the theoretical basis.Finally，to ensure the reliability of the results，the results from Software and practical calculation are compared and the errors are analyzed.The implementation of this algorithm proves that the security of vehicles can be ensured during the realization of lightweight production.

isotropic；strength calculation；Matlab programming

TQ9

Adoi：10.3969/j.issn.1002-6673.2015.05.026

1002-6673（2015）05-068-04

2015-07-14

李昕（1989-），陜西人，碩士研究生。研究方向：機(jī)械制造及其自動(dòng)化；沈斌（1955-），上海人，現(xiàn)任同濟(jì)大學(xué)中德學(xué)院機(jī)械信息技術(shù)基金教研室教授、主任、中國(guó)自動(dòng)化學(xué)會(huì)制造技術(shù)專業(yè)委員會(huì)委員、中國(guó)機(jī)電一體化應(yīng)用技術(shù)協(xié)會(huì)理事。研究領(lǐng)域：制造系統(tǒng)與自動(dòng)化技術(shù)、網(wǎng)絡(luò)化設(shè)計(jì)與制造及其相關(guān)技術(shù)、CAD/CAPP/CAM及PDM集成技術(shù)、智能制造技術(shù)。