吳　穹，吳　振

(沈陽(yáng)航空航天大學(xué) 遼寧省飛行器復(fù)合材料結(jié)構(gòu)分析與仿真重點(diǎn)實(shí)驗(yàn)室，遼寧 沈陽(yáng) 110136)

?

Reddy修正偶應(yīng)力模型及層合板熱尺度效應(yīng)

吳穹，吳振

(沈陽(yáng)航空航天大學(xué) 遼寧省飛行器復(fù)合材料結(jié)構(gòu)分析與仿真重點(diǎn)實(shí)驗(yàn)室，遼寧 沈陽(yáng) 110136)

摘要：引入纖維和基體材料的材料尺度參數(shù)，提出了適于分析復(fù)合材料層合板熱尺度效應(yīng)的Reddy修正偶應(yīng)力理論。此理論中將轉(zhuǎn)角矢量看作獨(dú)立變量，曲率張量不再對(duì)稱，但可以退化為各向同性的修正偶應(yīng)力。使用Reddy位移場(chǎng)模式，考慮溫度在x軸、y軸方向引起的熱應(yīng)變，分析了簡(jiǎn)支板的熱彎曲問題。數(shù)值結(jié)果表明，對(duì)于細(xì)觀復(fù)合材料層合板熱彎曲問題，材料的尺度效應(yīng)不能忽略。

關(guān)鍵詞：Reddy型修正偶應(yīng)力理論；復(fù)合材料層合板；熱彎曲；尺度效應(yīng)

實(shí)驗(yàn)表明[1-2]，當(dāng)金屬材料尺寸進(jìn)入微米量級(jí)時(shí)，材料的強(qiáng)度和剛度有所增強(qiáng)，這種現(xiàn)象被稱為尺度效應(yīng)。經(jīng)典連續(xù)體力學(xué)無(wú)法解釋該現(xiàn)象，為了分析尺度效應(yīng)，各國(guó)學(xué)者相繼發(fā)展了偶應(yīng)力理論及應(yīng)變梯度理論。上世紀(jì)60年代，Mindlin等人[3-5]提出了簡(jiǎn)化的偶應(yīng)力模型。2002年，Yang等[6]重新定義曲率張量，引入偶應(yīng)力力矩平衡方程，使得應(yīng)變張量與應(yīng)力張量對(duì)稱，因而本構(gòu)關(guān)系中只需引入一個(gè)材料尺度參數(shù)。2006年P(guān)ark和Gao[7]研究了修正偶應(yīng)力的Bernoulli-Euler梁模型。2008年Ma等[8]建立了Timoshenko梁模型。Tsiatas[9]在2009年建立了偶應(yīng)力Kirchhoff板模型。2011年，Ma等[10]基于修正偶應(yīng)力建立了Mindlin板模型。2012年，Alireza Nateghi等[11]研究了基于修正偶應(yīng)力的功能梯度梁的熱尺度效應(yīng)。同年，M.Simsek等[12]對(duì)功能梯度Timoshenko梁作了靜彎曲分析。2012年，Asari等[13]基于修正偶應(yīng)力對(duì)功能梯度Timoshenko梁做了熱后屈曲的尺度效應(yīng)分析。2012年，陳萬(wàn)吉等[14]定義新的曲率張量并建立復(fù)合材料層合板的本構(gòu)關(guān)系，認(rèn)為新的曲率張量不再對(duì)稱，但可以與各向同性的曲率張量等價(jià)；建立了Reddy層合梁模型[15]，根據(jù)新修正偶應(yīng)力分析復(fù)合材料層合梁的自由振動(dòng)問題[16]；2013年基于該理論用整體局部位移模式建立層合梁模型[17]。綜上所述，各國(guó)學(xué)者基于修正偶應(yīng)力分析了復(fù)合材料彎曲和自由振動(dòng)問題，然而沒有涉及復(fù)合材料層合板熱尺度效應(yīng)?；诖?，本文推導(dǎo)了考慮熱應(yīng)變的修正偶應(yīng)力理論，并基于修正偶應(yīng)力分析復(fù)合材料層合板熱尺度效應(yīng)。

2理論模型

2.1位移模型

為了得到Reddy理論，首先使用如下初始位移場(chǎng)：

u(x,y,z)=u0(x,y)+zu1(x,y)+z2u2(x,y)+z3u3(x,y)

v(x,y,z)=v0(x,y)+zv1(x,y)+z2v2(x,y)+z3v3(x,y)

v(x,y,z)=w0(x,y)

(1)

其中，u0、v0和w0表示中面上點(diǎn)(x,y)的位移；u1和v1表示中面法線關(guān)于y軸和x軸的轉(zhuǎn)角；u2、u3、v2和v3是沿z坐標(biāo)的三次展開項(xiàng)。

通過使用橫向剪切自由表面條件，Reddy理論最終位移場(chǎng)可寫為

w(x,y,z)=w0(x,y)

(2)

溫度場(chǎng)如下：

ΔT=f(z)T(x,y)

(3)

其中，f(z)為沿厚度方向分布的溫度構(gòu)形，T(x,y)為面內(nèi)溫度分布函數(shù)。

ωz=0

(4)

2.2新修正偶應(yīng)力及本構(gòu)關(guān)系

Yang等[6]提出的修正偶應(yīng)力中，應(yīng)變張量和曲率張量是對(duì)稱的，定義如下：

(5)

應(yīng)力張量和偶應(yīng)力張量定義如下：

(6)

(7)

(8)

其中，

mxymyx}T

(9)

(10)

={εx-αxΔTεy-αyΔTγxyγxzγyzχx

χyxxyχyx}T

(11)

αx，αy分別表示沿x軸方向和y軸方向的熱膨脹系數(shù)；

(12)

經(jīng)坐標(biāo)變換，第k層在整體坐標(biāo)中的應(yīng)力應(yīng)變關(guān)系

(13)

其中，Qk=TkTCkTk,Tk為坐標(biāo)轉(zhuǎn)換矩陣。

(14)

m=cosφk,n=sinφk。

2.3平衡方程

根據(jù)最小勢(shì)能原理，可得到平衡方程。最小勢(shì)能原理可寫為，

δπp=δU-δW=0

(15)

(16)

δW=∫Ωf δudv+∫?ΩV δuds

(17)

其中，δπp表示彈性體總的勢(shì)能，δU和δW分別表示彈性體應(yīng)變能和外力所做的功，f為層合板上的體力，V為在板邊界處的邊界力。只考慮受熱載荷的影響，外載荷不再考慮，故f=V=0。

把公式(2)帶入幾何方程，得到宏觀部分的應(yīng)變，同時(shí)將公式(4)帶入(5)中得到偶應(yīng)力部分應(yīng)變。最終都帶入式子(16)中，利用變分得到平衡方程，如下：

(18)

對(duì)四邊簡(jiǎn)支板，位移邊界條件

在x=0和x=L處，v0=0，v1=0，w0=0

在y=0和y=L處，u0=0，u1=0，w0=0

2.4解析解

在橫截面(y,z)纖維斷面與基體相互作用，纖維相當(dāng)于是夾雜，其材料尺度參數(shù)遠(yuǎn)大于在截面(x,z)中以長(zhǎng)纖維和基體作為夾雜的材料尺度參數(shù)[14]，即b?m。因此對(duì)于復(fù)合材料層合板，可以假設(shè)m=0，以簡(jiǎn)化剛度矩陣。

由式子(12)和(14)，得到整體坐標(biāo)中的剛度系數(shù)，如下：

(19)

為滿足四邊簡(jiǎn)支方板的邊界條件，根據(jù)Navier方法，位移u0、u1、v0、v1和w0都以含未知系數(shù)的雙傅里葉級(jí)數(shù)展開。各位移函數(shù)表示如下：

(20)

將式子(19)中的剛度系數(shù)歸到式子(13)的矩陣Qk中，然后將式子(13)、(20)帶入到(18)中，最終可得到式子(20)中各項(xiàng)系數(shù)。

3算例

為了驗(yàn)證微觀尺度下復(fù)合材料層合板的尺度效應(yīng)，采用[0°/90°/0°]正交鋪設(shè)的四邊簡(jiǎn)支方板。長(zhǎng)度L=200μm,板厚h=25μm，材料的屬性參數(shù)：E1/E2=15,E2=10GPa，G12/E2=0.5,G22/E2=0.335 6,υ12=0.3,υ22=0.49,αx=0.015×10-6/K,αy=1×10-6/K。

取在截面(x=L/2，y=L/2)的主應(yīng)力，應(yīng)力表達(dá)式如下：

圖1　在板(x=L/2,y=L/2)處的主應(yīng)力δx和δy

圖2給出了橫向位移分布。在圖中可以看到偶應(yīng)力理論下板的變形要小于經(jīng)典彈性力學(xué)中板的變形。

圖3分別給出了層合板在x=L/2和y=L/2處板的轉(zhuǎn)角。在圖中可以看到隨著尺度參數(shù)的增大板的轉(zhuǎn)角在減小(逆時(shí)針為正)。

數(shù)值結(jié)果表明，復(fù)合材料進(jìn)入細(xì)觀尺度時(shí)，材料的尺度效應(yīng)不能忽略。

圖2　橫向位移(y=L/2)

圖3　轉(zhuǎn)角(x=L/2、y=L/2)

4總結(jié)

參考文獻(xiàn)(References)：

[1]FLECK N A,MULLER G M,ASHBY MF,et al.Strain gradient plasticity:theory and experiment[J].Acta Metall Mater,1994,42(2):475-487.

[2]LAM D C,YANG F,CHONG A C M,et al.Experiments and theory in strain gradient elasticity[J].Mech Phys Solids,2003,51(8):1477-508.

[3]MINDLIN R D,TIERSTEN H F.Effects of couple stresses in linear elasticity[J].Archive for Rational Mechanics and Analysis,1962,11(11):415-448.

[4]TIUPIN R A.Elastic materials with couple stresses[J].Arch Rational Mech Anal,1962,11(1):385-414.

[5]MINDLIN R D.Influence of couple-stresses on stress concentrations[J].Experimental Mechan-ics,1962,3(1):1-7.

[6]YANG F,CHONG A M,Lam D C C,et al.Couple stress based strain gradient theory for elasticity[J].International Journal of Solids and Structures,2002,39(10):2731-2743.

[7]PARK S K,GAO X L.Bernoulli-Euler beam model based on a modified couple stress theory[J].Micromech Microeng 2006;16(11):2355-2359.

[8]MA H M,GAO X L,Reddy J N.A microstructure dependent Timoshenko beam model basedon a modified couple stress theory[J].Journal of the Mechanics and Physics of Solids,2008,56(12):3379-3391.

[9]TSIATAS G C.A new Kirchhoff plate model based on a modified couple stress theory[J].International Journal of Solids and Structures,2009,46(13):2757-2764.

[10]MA HM,GAO X L,REDDY J N.A non-classical Mindlin plate model based on a modified coup-le stress theory[J],2011(220):217-235.

[11]ALIREZA N,MAZAHER S.Thermal effect on size dependent behavior of functionally graded microbeams based on modified couple stress theory[J].Composite Structure,2012，8:97-110.

[12]SIMSEK M,OCATURK T K,AKBAS S D.Static bending of a functionally graded microscale Timoshenko beam based on the modified couple stress theory.Composite Structures,2013，95(1):740-747.

[13]ANSARI R,FAGHIH M ,GHOLAMI S R,et al.Thermal postbuckling behavior of size-dependent functionally graded Timoshenko microbeams[J].International Journal of Non-Linear Mechani-cs,2013,50(50):127-135.

[14]CHEN W J,MA X,LI L.A model of composite laminated Reddy plate based on a new modified couple stress theory[J].Composite Structures,2012,94(7):2143-2156.

[15]CHEN W J,CHEN W W,SZE K Y.A model of composite laminate Reddy beam based on a modified couple stress theory[J].Composite Structures,2012,94(8):2599-2609.

[16]CHEN W J,LI X P.Size-dependent free vibration analysis of composite laminated Timoshenko beam based on new modified couple stress theory[J].Arch Appl Mech,2012,83(3):431-444.

[17]CHEN W J,SI J L.A model of composite laminated beam based on the global-local theory andnew modified couple stress theory[J].Composite Structures,2013,103(3):99-107.

(責(zé)任編輯：宋麗萍英文審校：隋華)

Thermal size effect of laminated composite plate based on Reddy-type modified couple stress theory

WU Qiong,WU Zhen

(Key Laboratory of Liaoning Province for Composite Structural Analysis of Aircraft and Simulation,Shenyang Aerospace University,Shenyang 110136,China)

Abstract：By introducing material length scale parameter of fiber and matrix,a Reddy-type modified couple stress theory is developed to analyze the thermal size effect of composite laminated plate.In the proposed model,the rotation vectors are considered as independent variables,and curvature tensor is asymmetric.However,this model is still adaptable to the isotropic material.By considering thermal strain induced by temperature in X-axis and Y-axis,the thermal bending of simply-supported plate is analyzed by using the Reddy-type modified couple stress model.Numerical results show that the size effect of material has to be considered for the thermal bending of laminated composite micro-plate.

Key words：Reddy-type modified couple stress theory;laminated composite plate;thermal bending;size effect

doi:10.3969/j.issn.2095-1248.2016.01.003

中圖分類號(hào)：TP391.7; V261.2

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：2095-1248(2016)01-0013-05

作者簡(jiǎn)介：吳穹(1990-)，男，遼寧沈陽(yáng)人，碩士研究生，主要研究方向：多層度層合板熱分析，E-mail:705180054@qq.com；吳振(1977-)男，黑龍江佳木斯人，博士，教授，博士生導(dǎo)師，主要研究方向：復(fù)合材料力學(xué)研究，E-mail: wuzhenhk@163.com。

基金項(xiàng)目：國(guó)家自然科學(xué)基金(項(xiàng)目編號(hào)：11272217)

收稿日期：2015-03-10