董繼蕾 繆廣紅 張宏學(xué) 盧小雨

【摘?? 要】? 以兩根等長等寬等高的矩形截面異質(zhì)雙材料疊合懸臂梁為對象，研究雙材料疊合懸臂梁受集中荷載作用時(shí)的解析解和有限元解。先利用彈性理論對雙材料疊合懸臂梁進(jìn)行內(nèi)力分析，然后利用有限元分析軟件對雙材料疊合懸臂梁進(jìn)行數(shù)值模擬分析，并將模擬結(jié)果和應(yīng)力解析解進(jìn)行比較。結(jié)果表明：疊合梁上、下兩層梁中，彎矩-剪力是按照其剛度來分布的，當(dāng)上層梁自由端上所承受的外力總和大于或等于某值時(shí)，上層梁與下層梁的彎矩才能維持相對穩(wěn)定的接觸；解析解和有限元解比較一致，驗(yàn)證了彈性理論應(yīng)力解析解的正確性。

【關(guān)鍵詞】?? 疊合梁；集中荷載；彈性理論；有限元分析

Analytical Solution and Finite Element Analysis of Bimaterial Composite Cantilever Beams Subjected to Concentrated Loads

Dong Jilei， Miao Guanghong， Zhang Hongxue， Lu Xiaoyu

（Anhui University of Science and Technology， Huainan 232001，China）

【Abstract】??? The analytical and finite element solutions of two heterogeneous bimaterial composite cantilever beams with rectangular cross-sections of equal length， width， and height under concentrated loads are studied. Firstly， the internal force analysis of the bimaterial composite cantilever beam is conducted by using elastic theory. Then， numerical simulation analysis of the bimaterial composite cantilever beam is carried out by using finite element analysis software， and the simulation results are compared with the stress analytical solution. The results show that in the upper and lower beams of the composite beam， the bending moment shear force is distributed according to its stiffness. When the total external force borne by the free end of the upper beam is greater than or equal to a certain value， the bending moment of the upper beam and the lower beam can maintain relatively stable contact; the analytical solution is consistent with the finite element solution， which verifies the correctness of the stress analytical solution in elastic theory.

【Key words】???? composite beams; concentrated loads; elastic theory; finite element analysis

〔中圖分類號〕? TB301???????????? ?? ????? ???〔文獻(xiàn)標(biāo)識碼〕? A????? ???????????? 〔文章編號〕 1674 - 3229（2024）02- 0052 - 05

[收稿日期]?? 2023-12-20

[基金項(xiàng)目]?? 國家自然科學(xué)基金“基于濕-熱-彈本構(gòu)和剛?cè)狁詈侠碚摰淖儤D(zhuǎn)風(fēng)力機(jī)葉片動力學(xué)特性研究”（119002002）；2021年度安徽省高等學(xué)校省級質(zhì)量工程項(xiàng)目（2021xxkc033）；省級研究生線下示范課程：彈塑性力學(xué)（2022xxsfkc023）

[作者簡介]?? 董繼蕾（1990- ），女，碩士，安徽理工大學(xué)力學(xué)與光電物理學(xué)院實(shí)驗(yàn)師，研究方向：力學(xué)實(shí)驗(yàn)教學(xué)與學(xué)科競賽。

0???? 引言

疊合梁是由兩種或兩種以上材料組合而成的梁構(gòu)件，多種材料疊合而成的疊合梁能充分發(fā)揮每一種材料的性能，在力學(xué)性能、經(jīng)濟(jì)效益和綠色環(huán)保等方面展示出單種材料梁無法比擬的優(yōu)勢［1］。

由于疊合梁的種種優(yōu)勢，其在工程應(yīng)用中有很高的實(shí)用價(jià)值，目前已有大量關(guān)于疊合梁的研究成果。易超［2］根據(jù)材料力學(xué)、彈性力學(xué)相關(guān)理論對受任意線性荷載下的懸臂梁進(jìn)行分析，提出了滿足相容方程和邊界條件的應(yīng)力函數(shù)表達(dá)式，采用半逆解法得到了懸臂梁的應(yīng)力和位移解。譚宸［3］基于MATLAB平臺，開發(fā)了一個(gè)結(jié)構(gòu)良好、易于使用的計(jì)算懸臂梁結(jié)構(gòu)彈性變形的有限單元分析程序。利用該程序，可以進(jìn)一步加深學(xué)生對采用有限元法分析問題編制程序的理解。萬澤青等［4］為分析由不同材料組成的任意多層合梁的彎曲性能，以一矩形截面層合梁為例，采用相當(dāng)截面法確定了組合截面中性軸的位置，利用靜力平衡條件建立任意多層復(fù)合材料層合梁橫截面上任一點(diǎn)處彎曲切應(yīng)力的一般公式。劉一華等［5］采用彈性力學(xué)的應(yīng)力函數(shù)法，求得了雙材料疊合懸臂梁在自由端受集中力作用時(shí)的理論解，所得的應(yīng)力理論解可退化為已有的單材料懸臂梁在自由端受集中力作用時(shí)的理論解。郝英奇等［6］研究了鋼纖維混凝土組合材料中應(yīng)力波演化機(jī)理及衰減規(guī)律。Bouazaoui.L等［7］對粘結(jié)鋼—混凝土組合梁進(jìn)行了一系列的力學(xué)性能試驗(yàn)，將薄板鋼梁和混凝土板材用膠粘劑進(jìn)行粘接組合，結(jié)果表明，該方法能夠有效地提高混凝土構(gòu)件的承載能力。Nasser-Eddine等［8］通過將兩個(gè)相同梁的兩端面接觸連接起來，驗(yàn)證通過邊界指數(shù)形式可以控制穩(wěn)定系統(tǒng)。

以上研究主要集中在兩種異質(zhì)材料的疊層結(jié)構(gòu)上，國內(nèi)外關(guān)于雙材料疊合梁受荷載作用時(shí)的解析解和有限元分析還很少，采用彈性理論進(jìn)行分析的研究更是鮮見。本文將在已有研究成果的基礎(chǔ)上，利用彈性理論對雙材料疊合梁受集中荷載作用時(shí)進(jìn)行解析分析和有限元分析，以檢驗(yàn)其彈性力學(xué)解析解的正確性，并將有限元分析與理論分析進(jìn)行對比。

1???? 疊合懸臂梁彈性理論解析解

如圖1所示，設(shè)由兩根等長等寬等高異質(zhì)材料的矩形截面梁疊合的疊合懸臂梁，疊合懸臂梁的長度為[L]，寬度為[b]，高度為[h]，在自由端受集中向下的集中力[P]作用。假設(shè)上梁的彈性模量為[E1]，泊松比為[μ1]；下梁的彈性模量為[E2]，泊松比為[μ2]。

1.1?? 疊合懸臂梁軸力

上、下層梁中的軸力分別為：

[FN1=-h0σx1bdy=0FN2=0hσx2bdy=0]??????? （1）

其中：[FN1]、[FN2]分別為上層梁和下層梁的軸力值；[σx1]、[σx2]分別為上層梁和下層梁的正應(yīng)力值。

1.2?? 疊合懸臂梁剪力

上、下層梁中的剪力分別為：

[FQ1=-h0τxy1bdy=E1I1E1I1+E2I2FFQ2=0hτxy2bdy=E2I2E1I1+E2I2F] （2）

其中：[FQ1]、[FQ2]分別為上層梁和下層梁的剪力值；[τxy1]、[τxy2]分別為上層梁和下層梁的剪應(yīng)力值。

1.3?? 疊合懸臂梁彎矩

上、下層梁中的彎矩分別為：

[M1=-h0σx1ybdy=E1I1E1I1+E2I2Fx-lM2=0hσx2ybdy=E2I2E1I1+E2I2Fx-l] （3）

其中：[M1]、[M2]分別為上層梁和下層梁的彎矩值；[I1=I2=bh312]分別為上層梁和下層梁的慣性矩。

通過對疊合梁橫截面受力情況的分析，可以看出在疊合梁上、下兩層中，彎矩-剪力是按照其剛度來分布的，也就是說，彎矩-剪力分布在疊合梁上、下兩層中。而且，在梁的端部，當(dāng)上層梁所受的外力之和[F1≥E1I1FE1I1+E2I2]時(shí)，疊合梁在荷載作用下發(fā)生變形，上層梁和下層梁結(jié)構(gòu)之間始終存在著一定的接觸。相反，在上面受到的外力的總和[F1]滿足[F1≤E1I1FE1I1+E2I2]時(shí)，上層梁和下層梁兩個(gè)層間存在著相互分離和獨(dú)立的變形過程，其變形過程由有限元軟件得到，如圖2所示。

2???? 有限單元法結(jié)果檢驗(yàn)

有限元法是一種基于有限元理論的數(shù)值計(jì)算方法，它的基本思想是：將一個(gè)連續(xù)的結(jié)構(gòu)離散為有限多個(gè)單元，在每一個(gè)單元內(nèi)設(shè)置有限多個(gè)結(jié)點(diǎn)，將一個(gè)連續(xù)體系視為一系列僅在結(jié)點(diǎn)上相連的單元所組成的整體。采用有限元方法，以場函數(shù)的結(jié)點(diǎn)值為基元，對每個(gè)單元引入一種近似的內(nèi)插函數(shù)，來描述各單元內(nèi)的場函數(shù)分布，并運(yùn)用變分原理，構(gòu)建有限元方程，將連續(xù)區(qū)域內(nèi)的無窮自由度問題轉(zhuǎn)化為離散區(qū)域內(nèi)的有限自由度問題。

在此基礎(chǔ)上，本文提出了一種新的數(shù)值分析方法。有限元法在固體力學(xué)、流體力學(xué)、熱傳導(dǎo)、電磁、聲學(xué)、生物力學(xué)等方面有著重要的應(yīng)用。有限元法有如下特點(diǎn)：

（1）整個(gè)系統(tǒng)離散為有限個(gè)元素。

（2）將能量最低原則和函數(shù)的數(shù)值性質(zhì)聯(lián)系起來，建立相應(yīng)的數(shù)學(xué)模型。

（3）處理過程非常簡明。

（4）線性、非線性均適用。

（5）離散處理需要在整個(gè)局部區(qū)域進(jìn)行，這個(gè)過程需要大量的輸出空間和計(jì)算機(jī)內(nèi)存，并且解決問題需要相當(dāng)長的時(shí)間。

（6）無限區(qū)域的問題較難仿真。

利用有限元分析對兩種材料的疊合懸臂梁進(jìn)行數(shù)值計(jì)算，以檢驗(yàn)其彈性力學(xué)解析解的正確性，并將結(jié)果與理論分析結(jié)果進(jìn)行對比。內(nèi)力分析結(jié)果表明，只有在上部自由端面上承受的荷載足夠大，才能保證上部自由端面上的應(yīng)力集中。由圣維南原理可知［9］，在梁上施加集中力或等效荷載時(shí)，其應(yīng)力變化微乎其微，為了分析方便，采用在右端加集中力方式進(jìn)行，按上、下層中不同的彎曲剛度（[E1I1≥E2I2]、[E1I1≤E2I2]）情況進(jìn)行數(shù)值求解，然后將數(shù)值解與彈性力學(xué)所得的解析解進(jìn)行比較［10］。

2.1[?? E1I1≥E2I2]

如圖3所示計(jì)算模型1，假設(shè)疊合梁的長度[L=160 mm]，上下層材料分別為鋼和鋁，其中上層鋼的彈性模量[E1=206 GPa]，泊松比[μ1=0.3]；下層鋁的彈性模量[E2=70 GPa]，泊松比[μ2=0.3]。上、下層梁的高度[h=23 mm]，且梁的寬度相同[b=24 mm]，自由端受向下的集中荷載[F=2 KN]。

有限元分析計(jì)算結(jié)果的位移云紋圖、應(yīng)力云紋圖如圖4-圖7所示。

將[yi]分別取-23、-18、-12、-6、-0、+0、6、12、18、23[ mm]處的[σxi]，[σyi] 和[τxyi]的理論值計(jì)算出來，然后將同一處的理論值和有限元值進(jìn)行比較，計(jì)算出相對誤差，相對誤差的計(jì)算公式為：（理論值-有限元值）/理論值×100%，計(jì)算結(jié)果如表1所示。

2.2?? [E1I1≤E2I2]

如圖8所示計(jì)算模型2，組合梁的長度[L=160 mm]，上、下層梁的材料分別為鋁和鋼，上層鋁的彈性模量[E1=70 GPa]，泊松比[μ1=0.3]；下層鋼的彈性模量[E2=206 GPa]，泊松比[μ2=0.3]。其中上、下層的高度[h=23 mm]，梁的寬度[b=24 mm]，自由端受向下集中荷載[F=2 KN]。

有限元模擬計(jì)算得出的應(yīng)力云紋圖、位移云紋圖如圖9-圖13所示。

將[yi]分別取-23、-18、-12、-6、-0、+0、6、12、18、23 [mm]處的[σxi]，[σyi] 和[τxyi]的理論值計(jì)算出來，然后將同一處的理論值和有限元值進(jìn)行比較，計(jì)算出相對誤差，相對誤差的計(jì)算公式為：（理論值-有限元值）/理論值×100%，計(jì)算結(jié)果如表2所示。

如上所述，雙材料疊合懸臂梁的彎曲應(yīng)力分布是線性的，在彎曲剛度越大的組合梁上，應(yīng)力值越大。在疊合結(jié)構(gòu)中，層間存在著相對滑移，并且在任意一段變形后，層間均不再是平直結(jié)構(gòu)。用有限元方法對其進(jìn)行了計(jì)算，結(jié)果表明，理論與數(shù)值研究結(jié)果之間誤差較小，證明理論分析結(jié)果是正確的。

3???? 結(jié)論

（1）根據(jù)雙材料疊合懸臂梁上層梁和下層梁結(jié)構(gòu)中的受力情況，發(fā)現(xiàn)上層梁和下層梁的剪力和彎矩均按剛度分布，且只有在上部自由端受到的外力之和大于或等于某一值的情況下，上下兩層才能始終保持接觸。

（2）利用 ANSYS軟件，對上層梁和下層梁的剛度比值進(jìn)行分析，對上層梁和下層梁不同的彈性模量進(jìn)行了詳細(xì)的數(shù)據(jù)分析，得到以下結(jié)論：①在不考慮上層梁和下層梁之間存在摩擦力的情況下，所得到的計(jì)算結(jié)果與理論計(jì)算的結(jié)果最為一致；②如果上層梁的彎曲剛度比下部的彎曲剛度小，則計(jì)算結(jié)果的誤差較大；③當(dāng)上部和下層梁的彈性模量增加時(shí)，計(jì)算結(jié)果和分析結(jié)果之間的誤差將增加。

［參考文獻(xiàn)］

［1］ 郭斌強(qiáng)，黨濤，晉洋濤，等.超寬鋼混疊合梁施工過程和力學(xué)性能研究［J］.建筑施工，2023，45（6）：1256-1259.

［2］ 易超.任意線性荷載下懸臂梁半逆解法研究［J］.四川建筑，2023，43（6）：165-168.

［3］ 譚宸.基于MATLAB平臺的懸臂梁靜力彈性分析［J］.山西建筑，2023，49（20）：42-45+50.

［4］ 萬澤青，陶陽，夏鈺.雙材料層合梁的彎曲切應(yīng)力分析［J］.青海大學(xué)學(xué)報(bào)，2020，38（2）：81-86.

［5］ 劉一華，朱立軍.雙材料疊合懸臂梁端部受集中力作用的理論解［J］.合肥工業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版），2010，33（3）：391-395.

［6］ 郝英奇，葉中豹，周新貴.鋼纖維混凝土材料中應(yīng)力波演化機(jī)理及衰減規(guī)律研究［J］.安徽科技學(xué)院學(xué)報(bào)，2023，37（1）：63-68.

［7］ Bouazaoui L ， Perrenot G ， Delmas Y ，et al.Experimental study of bonded steel concrete composite structures［J］.Journal of Constructional Steel Research， 2007， 63（9）：1268-1278.

［8］ Nasser-Eddine T. Stabilization of a laminated beam with interfacial slip by boundary controls［J］.Boundary Value Problems， 2015，2015（1）： 1-11.

［9］ 劉延彬.淺析平面運(yùn)動剛體的外力之功［J］.廊坊師范學(xué)院學(xué)報(bào)（自然科學(xué)版），2021，21（2）：43-45.

［10］羅吉安，周星越.有限單元法分析中應(yīng)力奇異問題的研究［J］.廊坊師范學(xué)院學(xué)報(bào)（自然科學(xué)版），2021，21（4）：42-45+50.