吳曉，黃翀，孫晉

(湖南文理學院 土木建筑工程學院，湖南 常德，415000)

在工程實際中，石墨、增強復(fù)合材料、金屬合金、陶瓷、玻璃、鑄鐵等許多材料都具有拉壓彈性模量不同的雙模量性質(zhì)，所以，用雙模量本構(gòu)關(guān)系對這些材料制成的結(jié)構(gòu)進行計算分析已備受關(guān)注[1-3]。對于拉壓彈性模量不同的雙模量材料，彈性系數(shù)不僅依賴于結(jié)構(gòu)材料，而且與結(jié)構(gòu)材料、形狀、邊界條件及外載荷有關(guān)[4-6]。在梁、彈性平面等問題的結(jié)構(gòu)中，人們考慮了材料的雙模量特性[7-10]，并采用Kantorovich法研究了柱形桿的扭轉(zhuǎn)問題，但未見采用Kantorovich法研究雙模量簡支梁的平面應(yīng)力問題的報道。經(jīng)典彈性理論研究線性分布荷載作用下簡支梁的平面應(yīng)力問題多采用半逆法或三角級數(shù)法，而采用半逆法時需要確定的待定常數(shù)較多，采用三角級數(shù)法時存在收斂慢且計算過程復(fù)雜繁瑣等缺陷。為此，本文作者采用Kantorovich法研究在線性分布荷載作用下雙模量簡支梁的平面應(yīng)力問題。

1 材料力學應(yīng)力表達式

由于雙模量梁在外荷載作用下彎曲時，會形成彈性模量不同的拉伸區(qū)和壓縮區(qū)，由彈性理論可知雙彈性模量梁彎曲時的應(yīng)力和應(yīng)變關(guān)系為：

式中：σ1為拉伸區(qū)應(yīng)力；σ2為壓縮區(qū)應(yīng)力；E1為拉伸區(qū)的彈性模量；E2為壓縮區(qū)的彈性模量；y為計算點至中性軸的距離；ρ為中性層的曲率半徑。

據(jù)彈性理論，雙模量梁彎曲時橫截面內(nèi)力應(yīng)滿足以下關(guān)系：

將式(1)代入(2)可得：

圖1中，雙模量簡支梁中點處的彎矩為

式中：q0為梁上線性分布載荷最小值；q1為梁上線性分布載荷最大值；l為梁的跨長。

圖1 雙模量簡支梁Fig.1 Bimodulous simply supported beam

因此，圖1 中雙模量簡支梁中點處的最大拉應(yīng)力σ1及最大壓應(yīng)力σ2分別為：

2 簡支梁的Kantorovich解

對于圖1所示線性分布荷載作用下的雙模量簡支梁，其線性分布荷載的集度表達式為

雙模量簡支梁上下面邊界條件及中性軸上的應(yīng)力條件分別為：

假設(shè)線性分布荷載作用下雙模量簡支梁受拉區(qū)應(yīng)力函數(shù)φ1(x,y)和受壓區(qū)應(yīng)力函數(shù)φ2(x,y)分別為：

式中：Y1(y)為梁拉伸區(qū)變量y的貢獻函數(shù)；Y2(y)為梁壓縮區(qū)變量y的貢獻函數(shù)。

則線性分布荷載作用下雙模量簡支梁的應(yīng)力表達式為：

式中：“+”表示拉伸區(qū)；“－”表示壓縮區(qū)。

對于平面應(yīng)力問題，線性分布荷載作用下雙模量簡支梁的余能為

將式(10)代入式(11)可得：

對式(12)進行一階變分可得：

利用微分符號與變分符號可交換性及邊界條件式(8)可得：

將函數(shù)g1(y,Y1,Y1′,Y1′′) 和g2(y,Y2,Y2′,′ )表達式代入式(8)得：

由式(15)可以求得Y1(y)和Y2(y)的表達式：

將式(16)代入式(10)可求得拉伸區(qū)及壓縮區(qū)應(yīng)力表達式為：

利用式(17)及邊界條件式(8)可得：

常數(shù)K1，K2，K3，K4，K5，K6，K7，K8，K9，K10，K11，K12，R1，R2，R3，R4和R5分別為：

3 算例分析與討論

首先采用有限元法檢驗本文 Kantorovich解的計算精度。假設(shè)某雙模量簡支梁的b=1，E1=93.2 GPa，E2=124.36 GPa，采用式(17)和(6)及有限元法計算雙模量簡支梁中點處的最大拉應(yīng)力及最大壓應(yīng)力。其中：采用有限元建立雙模量簡支梁實體模型，模型由2層梁組成，梁高為100 mm，其中上層梁高46.4 mm，材料模型為mat1，E2=124.36 GPa；下層梁高為53.6 mm，材料模型為 mat2，E1=93.2 GPa。單元為 8節(jié)點SOLID185單元。具體的計算結(jié)果見表1～4。其中：q0和q1見圖1；σ1為采用材料力學計算的最大拉應(yīng)力；σ2為采用材料力學計算的最大壓應(yīng)力。

表1 q1=q0時采用本文方法所得雙模量簡支梁中點處彎曲應(yīng)力Table 1 Bending stress of bimodulous simply supported beam at midpoint obtained by method in the paper at q1=q0 MPa

表2 q1=2q0時采用本文方法所得雙模量簡支梁中點處彎曲應(yīng)力Table 2 Bending Stress of bimodulous simply supported beam at midpoint obtained by method in the paper at q1=2q0 MPa

表3 q1=q0時采用材料力學方法所得雙模量簡支梁中點處彎曲應(yīng)力Table 3 Bending Stress of bimodulous simply supported beam at midpoint obtained by method of material mechanics at q1=q0 MPa

表4 q1=2q0時采用材料力學方法所得雙模量簡支梁中點處彎曲應(yīng)力Table 4 Bending Stress of bimodulous simply supported beam at midpoint obtained by method of material mechanics at q1=2q0 MPa

由表1～2中本文方法與ANSYS計算結(jié)果的比較可知：采用Kantorovich法研究雙模量簡支梁中點處彎曲應(yīng)力的計算結(jié)果與有限元法的計算結(jié)果很接近，這說明采用 Kantorovich法研究雙模量簡支梁彎曲應(yīng)力的計算精度高，所得雙模量簡支梁的應(yīng)力公式是可靠的。

由表1～2中采用相同彈性模量彈性理論計算結(jié)果與采用雙模量彈性理論計算結(jié)果的比較可知：采用相同彈性模量彈性理論得到雙模量簡支梁中點處彎曲應(yīng)力，與采用雙模量彈性理論得到雙模量簡支梁中點處彎曲應(yīng)力的誤差均在 5%以上。采用相同彈性模量彈性理論研究雙模量簡支梁的彎曲應(yīng)力，拉壓區(qū)的彎曲應(yīng)力絕對值相等。事實上，雙模量簡支梁的拉壓區(qū)的彎曲應(yīng)力絕對值不相等，隨雙模量簡支梁拉壓區(qū)的彈性模量變化而變化。在本文中，雙模量簡支梁的拉壓區(qū)的彎曲應(yīng)力相差均達10%以上，所以，對于雙模量簡支梁的平面應(yīng)力問題，不宜采用相同彈性模量彈性理論求解，而應(yīng)該采用雙模量彈性理論求解。

由表1～2可知：在外載荷作用下，隨著雙模量簡支梁長高比的增大，雙模量簡支梁的拉壓區(qū)的彎曲應(yīng)力也隨著增大。因為隨著雙模量簡支梁長高比的增大，將導致線性分布荷載對雙模量簡支梁產(chǎn)生的彎矩也相應(yīng)增大。

比較表1與表2中的計算結(jié)果可知：在線性分布荷載(q1=2q0)作用下，雙模量簡支梁中點處的彎曲應(yīng)力要大于在均布載荷(q1=q0)作用下雙模量簡支梁中點處的彎曲應(yīng)力。這是因為線性分布荷載對雙模量簡支梁產(chǎn)生的彎矩大于均布載荷對雙模量簡支梁產(chǎn)生的彎矩。

本文方法計算結(jié)果與材料力學方法計算結(jié)果見表3和表4。從表3和表4可以看出：采用材料力學方法研究雙模量簡支梁的彎曲應(yīng)力，得到的雙模量簡支梁拉伸區(qū)的彎曲應(yīng)力與本文 Kantorovich法得到的雙模量簡支梁拉伸區(qū)的彎曲應(yīng)力很接近；但是，采用材料力學方法得到的雙模量簡支梁的彎曲應(yīng)力與本文Kantorovich法得到的雙模量簡支梁壓縮區(qū)的彎曲應(yīng)力相差很大，相對誤差均在50%以上，并且有限元法的計算結(jié)果也證明了這一點。這說明采用材料力學方法研究雙模量簡支梁的彎曲應(yīng)力有其局限性。

4 結(jié)論

(1) 采用Kantorovich法研究雙模量簡支梁彎曲應(yīng)力的計算精度高，所得雙模量簡支梁的應(yīng)力公式是可靠的。

(2) 采用材料力學方法研究雙模量簡支梁的彎曲應(yīng)力有其局限性。對于雙模量簡支梁的平面應(yīng)力問題，不宜采用相同彈性模量彈性理論求解，而應(yīng)該采用雙模量彈性理論求解。

(3) 在外載荷作用下，隨著雙模量簡支梁長高比的增大，雙模量簡支梁的拉壓區(qū)的彎曲應(yīng)力也隨著增大；在線性分布荷載作用下，雙模量簡支梁的彎曲應(yīng)力大于在均布載荷作用下雙模量簡支梁的彎曲應(yīng)力。

[1] Medri G. A nonlinear elastic model for isotropic materials with different behavior in tension and compression[J]. Transactions of the ASME, 1982, 104(1): 26-28.

[2] Bert C W, Reddy J N, Chao W C, et al. Vibration of thick rectangular plates of bimodulous composite material[J]. Journal of Applied Mechanics, 1981, 48(2): 371-376.

[3] Srinivasan R S, Ramachandra L S. Axisymmetric nonlinear dynamic response of bimodulous annular plates[J]. Journal of Vibration and Acoustics, 1990, 112(2): 202-205.

[4] 李戰(zhàn)莉, 黃再興. 雙模量泡沫材料等效彈性模量的細觀力學估算方法[J]. 南京航空航天大學學報, 2006, 38(4): 464-468.

LI Zhanli, HUANG Zaixing. Meso-mechanical method for estimating equivalent elastic modulus of foam-solid with double-modulus[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2006, 38(4): 464-468.

[5] 曾紀杰. 對中柔度壓桿的雙模量理論的修正[J]. 機械強度,2006, 28(3): 462-464.

ZENG Jijie. Revision of the formula with bimodulusim intermediate column[J]. Journal of Mechanical Strength, 2006,28(3): 462-464.

[6] 蔡來生, 俞煥然. 拉壓模量不同彈性物質(zhì)的本構(gòu)[J]. 西安科技大學學報, 2009, 29(1): 17-21.

CAI Laisheng, YU Huanran. Constitutive relation of elastic materials with different elastic moduli in tension and compression[J]. Journal of Xi'an University of Science and Technology, 2009, 29(1): 17-21.

[7] 吳曉, 楊立軍, 孫晉. 雙模量圓板彎曲變形的計算分析[J]. 西安建筑科技大學學報: 自然科學版, 2009, 41(1): 88-92.

WU Xiao, YANG Lijun, SUN Jin. Bending deformation calculation of bimodulous circular plate[J]. Journal of Xi'an University of Architecture & Technology: Natural Science Edition, 2009, 41(1): 88-92.

[8] 羅戰(zhàn)友, 夏建中, 龔曉南. 不同拉壓模量及軟化特性材料的柱形孔擴張問題的統(tǒng)一解[J]. 工程力學, 2008, 25(9): 79-84.

LUO Zhanyou, XIA Jianzhong, GONG Xiaonan. Unified solution for expansion of cylindrical cavity in strain-softening materials with different elastic moduli in tension and compression[J]. Engineering Mechanics, 2008, 25(9): 79-84.

[9] 阿巴爾楚米揚. 不同模量彈性理論[M]. 鄔瑞鋒, 張允真, 譯.北京: 中國鐵道出版社, 1986: 11-22.

Ambartsumyan S A. Elasticity theory of different modulus[M].WU Ruifeng, ZHANG Yunzhen, transl. Beijing: China Railway Press, 1986: 11-22.

[10] 高潮, 劉相斌, 呂顯強. 用拉壓不同模量理論分析彎曲板[J].計算力學學報, 1998, 15(4): 448-455.

GAO Chao, LIU Xiangbin, Lü Xianqiang. Analysis for the plate with the theory of different extension compression elastic modulis[J]. Chinese Journal of Computaional Mechanics, 1998,15(4): 448-455.