拉壓不同模量矩形板的雙向彎曲問題

張良飛,姚文娟

(上海大學土木工程系,上海 200444)

拉壓不同模量矩形板的雙向彎曲問題

張良飛,姚文娟

(上海大學土木工程系,上海 200444)

拉壓不同模量矩形板的雙向彎曲的中性軸可以從兩個彎曲方向考慮.基于不同模量理論,利用靜力平衡方程推導了不同模量矩形板的中性軸位置,再利用Kantorovich變分法求解了不同模量矩形板的撓曲線方程,并將得到的數(shù)值解和有限元解進行比較,二者較為吻合.計算結果表明,當拉壓不同模量的差異較大時,不同模量彎曲矩形板的撓度不宜采用相同模量經(jīng)典板殼理論.該方法為分析不同模量矩形板和其他結構形式的板的彎曲問題提供了求解思路,并為其在工程中的應用提供了一定的理論參考.

Kantorovich變分法;不同模量;雙向彎曲

工程上廣泛應用的許多材料,如混凝土、金屬、石墨、塑料、生物材料等都具有拉壓不同模量特性,特別是近幾年發(fā)展起來的復合材料以及聚合材料的不同模量之比高達4,其雙模量特性所引發(fā)的力學性質(zhì)差異已經(jīng)到了不容忽視的程度.因此,由雙模量材料制成的結構力學研究備受關注.最初的雙模量問題的研究是由美國學者Jone[1]、Bert[2]及前蘇聯(lián)學者Ambartsumyan[3]提出,1986年Ambartsumyan[3]發(fā)表了拉壓不同模量彈性問題的專著,總結了雙模量的基本概念和基本假設,并提出了主應力的正負為拉壓判定準則.

早期,對于有限元的研究還主要基于有限元數(shù)值的方法.Rafaele等[4]和Patel等[5]基于Bert[2]提出的雙模量模型——以纖維的縱向應變的正負為拉壓判定準則,分析了雙模量層壓復合板應力計變形;針對有限元方法迭代速度慢的問題,He等[6]應用了剪切模量法提高了計算速度;Yang等[7-8]用光滑技術處理雙模量模型,避免了迭代中應力狀態(tài)的判斷從而簡化了計算,并建立了初應力有限元及神經(jīng)網(wǎng)絡預測模型[8].

根據(jù)雙模量材料特殊的非線性特點,Yao等利用流動坐標系統(tǒng)結合分段積分法推導了雙模量彎壓柱[9]、橫力彎曲梁[10]、擋土墻的中性軸[11]、應力和位移的解析解;何曉婷等[12-13]采用簡化等效剛度法對不同模量彎曲梁和柱進行了分析;Qu[14]通過引入Yao的等效彈性模量,推導了雙模量土工格室變形表達式.

對于薄板彎曲問題,高潮等[15]基于不同彈性理論,采用有限元數(shù)值方法分析了彎曲板;吳曉等[16-17]用變分法研究了不同模量圓板和矩形板的彎曲問題,在研究過程中簡化了對于中性軸位置確定的步驟,并假設雙向彎曲截面上的中性軸在同一平面.實際上不同模量板雙向彎曲時中性層不在同一平面上,且其位置的確定不僅依賴于結構材料,而且還與結構形狀、邊界條件、外荷載等有關,是諸多因素所致的具有非線性現(xiàn)象的力學問題,至今尚未有學者對此問題進行解析研究.本工作首先基于Ambartsumyan[3]的不同模量理論,利用靜力平衡方程推導出不同模量彎曲板的中性軸方程,計算得到中性軸位置的數(shù)值解;接著用Kantorovich變分法研究了不同模量彎曲板變形問題;最后分析了拉壓不同模量對四邊固支矩形板彎曲變形的影響.

1 拉壓不同模量彈性理論

1.1 基本概念和初始假設

在絕對值相同的拉應力和壓應力的作用下,材料會發(fā)生絕對值不同的拉應變和壓應變,即材料具有不同的拉伸彈性模量Et和壓縮彈性模量Ec.在軸向應力作用下,材料在拉伸和壓縮時的應力應變關系是非線性的.

假設所研究的物體是變形連續(xù)的固體,物質(zhì)是均勻的、各向同性的,但由于主應力符號的不同而產(chǎn)生不同的彈性性質(zhì),即一維拉伸時所研究材料的任何方向有彈性模量Et,一維壓縮時有Ec,材料在任意應力狀態(tài)下只發(fā)生彈性小變形,并服從連續(xù)彈性介質(zhì)的一般規(guī)律.經(jīng)典彈性理論的基本方程和關系式保持不變,不同之處僅反映在物理方程(本構關系)中.

研究對象為薄板,以如下兩個計算假定為基礎:①垂直于中面方向的正應變εz可以不計;②應力分量τzx,τzy和σz遠小于其他3個應力分量,其引起的形變可以不計.

1.2 應力應變關系

根據(jù)不同模量彈性理論,假定主應變和主應力分別為{σI}和{εI},在主方向上的本構關系為式中,[a]為柔度矩陣,是由主應力的符號確定;[D]為彈性矩陣,是由主應變的方向確定.且這些矩陣滿足如下關系:

由不同模量彈性理論本構關系可知,對于平面問題,[D]和[a]的矩陣單元可以表示如下.

2 中性軸位置的確定

圖1 計算模型Fig.1 Calculation model

由方程組(12)可知,不同模量薄板兩彎曲截面應力改變符號的位置(由η1,η2決定)與拉壓模量比、拉壓泊松比、板邊彎矩之比有關,方程組(12)是這些參量的非線性函數(shù),為二階非線性方程組,對其進行求解可以得到中性軸位置的數(shù)值解.

3不同模量彎曲薄板的Kantorovich變分解

利用微分符號和變分符號的可交換性,對式(19)分部積分得

4 算例和討論

4.1 算例與計算結果

薄板尺寸為a=2 000 mm,b=1 000 mm,h=200 mm,受均布荷載q=2 MPa,M1= 20×106N·m,M2=5×106N·m.為了具體分析拉壓彈性模量對矩形板彎曲變形的影響,取如下材料彈性模量計算:①保持E=(Ec+Et)/2=4 000 MPa不變,不同模量比Et/Ec在1.5～5.0范圍內(nèi)變化;②保持Ec=1 000 MPa不變,Et/Ec在1.5～5.0范圍內(nèi)變化.

分別采用本工作推求的不同模量薄板撓曲線公式以及不同模量有限元數(shù)值解對四邊固支矩形薄板進行求解,并用經(jīng)典板殼理論計算出同模量時板的撓度.有限元數(shù)值模擬用Ansys 12.0,采用8節(jié)點Solid185單元,利用小位移靜力分析進行求解(部分結果見表1).

表1 不同模量矩形板中性軸位置、位移計算結果Table 1 Neutral axis location and displacement of rectangular plates with diferent modulus

4.2分析與討論

運用本方法所推求的板中心處撓度與不同模量有限元模擬結果吻合較好(見表1及圖2),其誤差均在5%以內(nèi),驗證了本方法的可靠性.

計入不同模量后,隨著材料壓、拉模量比Et/Ec的變化,算例中薄板的板中心撓度呈有規(guī)律的變化(見圖2).當雙模量的平均值E=(Et+Ec)/2保持不變,隨著Et/Ec的增大,不同模量板中心撓度呈非線性增大:當Et/Ec從1.5增大到2.5時,相應的板中心撓度增大了4 mm;當Et/Ec從2.5增大到5.0時,板中心撓度增大了22 mm.如果Ec不變,當Et/Ec增大時,則板中心處的撓度呈非線性減小:當Et/Ec從1.5增大到3.5時,相應的板中心撓度減小21 mm;當Et/Ec從3.5增大到5.0時,板中心撓度減小3 mm.這說明當雙模量的平均值E=(Et+Ec)/2保持不變時,隨著Et的增大,則Ec減小,板整體剛度分配不均勻,導致板中心撓度增大;如果Ec不變,當Et/Ec增大時,板的整體模量增大,則板的撓度減小;但當Et/Ec繼續(xù)增大時,拉、壓模量的離散性在增大,即不均性在增大,因此板中心撓度曲線最終趨于平緩.

圖2 E=4 000 MPa和Ec=1 000 MPa時利用本方法和有限元求解的板中心處撓度Fig.2 Finite element method solution and the proposed solution of displacement at the center of the plate when E=4 000 MPa and Ec=1 000 MPa

運用本方法算得的不同模量板的撓曲線和運用經(jīng)典板殼理論算得的相同模量板撓曲線如圖3和4所示.從圖中可以看出,撓曲線沿板長邊方向均呈拋物線形狀.當E=(Et+Ec)/2不變時,隨著Et/Ec增大,不同模量板的撓度均大于相同模量板的撓度,且板撓度呈非線性增大趨勢,沿板邊向板中心方向板的撓度在增大.當Ec不變時,隨著Et/Ec增大,薄板的撓度非線性減小,不同模量板撓度均小于相同模量板的豎向位移,沿板中心向板邊方向板的撓度在減小.

圖3 E=4 000 MPa時板撓曲線Fig.3 Displacement curves when E=4 000 MPa

圖4 Ec=1 000 MPa時板撓曲線Fig.4 Displacement curves when Ec=1 000 MPa

在本算例中,當Et/Ec=2.5時,保持E不變的不同模量板中心撓度為?34 mm,運用經(jīng)典板殼理論算得相同模量板的中心撓度為?27.1 mm,二者之間相差達到25.5%,超過20%;當Ec不變,增大Et時亦如此,故此時不宜再按照相同模量經(jīng)典理論計算.

5 結論

(1)在矩形薄板x和y方向作用有不同的彎矩時,薄板雙向彎曲,中性層不在一個平面上.本工作將板沿板厚方向分3層(理論上分層越多,越接近真實受力狀態(tài)),每一層默認為同一類應力區(qū)域,引入不同模量的本構關系,得到不同方向截面上的中性軸,用Kantorovich變分法計算出薄板的撓度曲線表達式.結果表明,這種近似計算的方法和不同模量有限元模擬方法較為吻合,二者誤差小于5%.

(2)在保持E=(Et+Ec)/2和Ec不變的兩種情況下,隨著Et/Ec增大,板的撓度表現(xiàn)出不同的變化趨勢,前者呈非線性增大趨勢,后者減小后趨于平緩.

(3)在本算例中,當Et/Ec≥2.5時,運用經(jīng)典板殼理論算得的板位移值和不同模量數(shù)值解之間相差超過20%,此時不宜按照相同模量經(jīng)典理論計算,而應采用不同模量理論.

[1]JONE R M.Stress-strain relations for materials with diferent moduli in tension and compression[J].J AIAA,1977,15:16-23.

[2]BERTCW.Modelsforfbrouscompositeswithdiferentpropertiesintensionand compression[J].J Engg Math Tech(ASME),1977,99:344-349.

[3]AMBARTsuMYAN S A.Elasticity theory of diferent modulus[M].Beijing:China Railway Press, 1986.

[4]RAFFAELE Z,FABRIZIO G.Damage evolution in bimodular laminated composites under cyclic loading[J].Compo Struct,2001,53:381-402.

[5]PATEL B P,LELE A V,GANAPATHI M,et al.Thermo-fexural analysis of thick laminates of bimodulus composites materials[J].Compo Struct,2004,63:11-20.

[6]HE X T,ZHENG Z L,SuN J Y,et al.Convergence analysis of a fnite element method based on diferent moduli in tension and compression[J].Int J Solids Struct,2009,46:3734-3740.

[7]YANG H T,ZHu Y L.Solving elasticity problems with bi-modulus via a smoothing technique[J].Chinese J Comput Mech,2006,23:19-23.

[8]YANG H T,Xu M L.Solving inverse bimodular problems via artifcial neutral network[J]. Inverse Problems Sci Eng,2009,17:999-1017.

[9]YAO W J,YE Z M.Analytical solution of bending-compression column using diferent tensioncompression modulus[J].Applied Mathematics and Mechanics(English Edition),2004,25: 901-909.

[10]YAO M J,YE Z M.Analytical solution for bending beam subject to later for with diferent modulus[J].Applied Mathematics and Mechanics(English Edition),2004,25:1107-1117.

[11]YAO W J,YE Z M.The analytical and numerical solution of retaining wall based on elastic theory of diferent modulus[J].Journal of Shanghai Jiao Tong University,2004,38:1022-1027.

[12]何曉婷,陳山林,孫俊貽.不同模量簡支梁均布荷載下的彈性力學解[J].工程力學,2007,24(10): 51-56.

[13]何曉婷,鄭周練,陳山林.拉壓不同模量彎壓柱的近似彈性力學解[J].重慶大學學報,2008,31(3): 339-343.

[14]Qu C Z.Deformation of geocell with diferent tensile and compressive modulus[J].Elec J Geotech Eng,2009,14:1-14.

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

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

[17]吳曉,楊立軍,孫晉.雙模量圓板彎曲變形的計算分析[J].西安建筑科技大學學報,2009,41(4): 485-488.

Biaxial bending of rectangular plates with diferent modulus

ZHANG Liangfei,YAO Wenjuan
(Department of Civil Engineering,Shanghai University,Shanghai 200444,China)

An eutral axis can be considered from two bending directions to solve the biaxial bending problem of rectangular plates with diferent modulus.Based on the diferent modulus theories,an equation of the neutral axis location is derived using the static equilibrium equation of rectangular plates with diferent modulus.The defection curve equation is solved with the Kantorovich method.Its solution agrees well with the fnite element solution.The result shows that the classical shell theory is not applicable in calculating defection of a rectangular plate when tensile and compressive modulus are quite diferent.The proposed method provides an approach to analyze the bending problem of other structure forms of plate with diferent modulus,and theoretical reference for engineering applications.

Kantorovich variational method;diferent modulus;biaxial bending

TU 36

A

1007-2861(2017)01-0128-10

10.3969/j.issn.1007-2861.2015.02.003

2015-01-03

姚文娟(1957—),女,教授,博士生導師,博士,研究方向為結構工程.E-mail:wenjuan@mail.shu.edu.cn