多孔FGM矩形板的自由振动与临界屈曲载荷分析

滕兆春, 席鹏飞

(兰州理工大学 理学院, 甘肃 兰州 730050)

功能梯度材料(functionally gradient materials,FGM)是指材料的化学组成、微观结构、原子排列由一侧向另一侧呈连续梯度变化，从而使材料的性质和功能连续地呈梯度变化[1]，与传统的复合材料[2-3]相比，FGM具有能够确保应力分布平稳过度、增加不同材料界面的结合强度等特点。FGM最初被设计用作热障材料，后来被拓展到机械工程、光学、核能、土木工程、电子信息及航空航天等诸多领域用来提高产品的实际性能。作为功能和安全设计的前提，FGM结构的具体力学性能，一直倍受人们的关注和重视。Wang等[4]研究了全组分W/Cu功能梯度材料在稳态热载荷下的性能，研究结果表明，梯度过渡层确实可以改善热失配，从而延长疲劳寿命。Shishesaz等[5]基于应变梯度理论，研究了功能梯度纳米圆盘的热弹性行为，主要分析了温度变化对纳米圆盘应力和径向位移的影响。Stathopoulos等[6]研究了用于燃气轮机的功能梯度多层热障涂料的设计，发现功能梯度材料涂层有效地提高了燃气轮机的耐高温性能和热效率。

功能梯度材料在各领域的需求，促进了其制备技术的快速发展，各种制备方法先后出现，为制备具有功能梯度特性的材料提供理论依据和实验支撑。Wan等[7]用等离子产生高温高速热源，把熔融状态的原材料喷涂到基材上形成涂层，原料成分比例不同时可形成不同的功能涂层。但制备出的功能梯度材料有孔隙率高、层间结合力低、强度低的缺点，严重影响了功能梯度材料的性能。其他如SHS法[8]、薄膜积层法[9]等制备出的功能梯度材料也有孔隙率高、易脱落的缺陷。孔隙率对功能梯度材料的各种特性有着很大的影响。对于单相材料中孔隙率与弹性模量关系的理论研究[10-12]，研究者们提出了许多能够表示孔隙率与弹性常数关系的方程，这些方程主要适用于小孔隙率时的情况。Phani等[13]根据以往经验推导出了新的方程，该方程能较好地满足具有理想和非理想填充几何模型系统在不同孔隙率时弹性模量的精确理论解。而对于功能梯度材料这类多相材料，Ramakrishana等[14]通过研究孔隙率对体积模量、密度的影响，推导出了适用于大范围孔隙率的有效弹性模量和有效泊松比。Kiran等[15]用有限元法研究了多孔FG磁电弹性板的自由振动特性和静态特性，主要利用修正的混合率公式分析了孔隙率分布对板结构性能的影响以及对板的自由振动和静态特性的影响。Behravan等[16]基于径向和厚度方向上的微分求积和状态空间矢量技术的求解方法，研究了非对称Kerr弹性基础上受非轴对称和非均匀面力作用下变厚度多孔FGM圆板的三维磁弹性解。洪轲等[17]利用本征函数展开法，分析了孔隙率对FGM中Lamb波传播特性的影响。在最近的研究成果中，Fahsi等[18]提出了一种新的准三维剪切变形理论，并在此基础上分析了弹性地基上多孔FGM梁的弯曲、屈曲和自由振动，讨论了2种不同孔隙率分布和长细比对梁的跨中挠度、频率和临界屈曲载荷的影响。诸多研究显示，对含有孔隙FGM的力学性能进行研究具有重要意义，虽然已有一些关于多孔功能梯度结构的静动态问题的研究结果，但研究深度和广度仍显不足。

近些年来，研究者采用一种基于函数的Taylor公式展开来求解微分方程的新方法——微分变换法(differential transformation method,DTM)。该方法的特点是避开直接求未知函数的高阶导数值，通过对微分方程的微分变换得到相邻阶的高阶导数值之间的代数关系，再将初值代入代数关系递推得到不同阶导数的值。由于DTM较高的精度和计算效率，被逐渐引入到许多问题的求解中[19]。目前，对于四边受压多孔功能梯度材料矩形板的DTM求解还鲜有文献报道。

本文基于经典薄板理论，考虑了孔隙和梯度指数对FGM弹性常数的影响，采用DTM研究四边受压多孔FGM矩形板的自由振动和屈曲特性，并在退化后与现有文献值进行对比验证。在计算结果的基础上，进一步分析孔隙率、梯度指数，载荷等参数对多孔FGM矩形板振动特性和临界屈曲载荷的影响。

1 多孔FGM矩形板的有效物性参数

图1所示的笛卡尔坐标系中,多孔FGM矩形板厚度为h,长度为a,宽度为b,四边受到单位压载荷Nx和Ny的作用,xy平面为中面,z轴垂直于中面且与xy轴形成右手螺旋。在y=0和y=b处为简支边界(S),而在x=0和x=a处可为简支边界(S)、固定边界(C)或自由边界(F)。在本文的问题描述和讨论中,边界条件按照x=0,y=b,x=a,y=0的顺序给出。

图1 四面受压多孔FGM矩形板的几何模型

FGM在无孔隙状态下,其弹性模量E、密度ρ等沿着厚度方向的变化如下[20]

(1)

式中:n为梯度指数;Pc,Pm分别表示无孔隙时陶瓷和金属材料的物性参数。

如果FGM中存在孔隙时,可以假设孔隙的大小和分布是随机的,并且介质是连续致密的,在此基础上得到孔隙率θ与弹性模量E、泊松比ν和密度ρ的关系[12]。

孔隙率

(2)

式中:ρ(z)为致密密度；ρ0(z)为自然密度。

有效泊松比

(3)

有效弹性模量

(4)

式中:E(z)和E*(z)分别为θ=0,0<θ<1时的弹性模量；bc=2-3ν。将(1)～(2)式代入(4)式即可得到不同孔隙率下FGM矩形板的有效弹性模量。

2 控制微分方程的无量纲化及DTM变换

2.1 控制微分方程及无量纲化

本文中物理中面[21]定义如下

(5)

根据物理中面和经典薄板理论所给出的本构关系,引入Hamilton原理[22]推导四边受压多孔FGM矩形板运动的控制微分方程。已知Hamilton原理如下

(6)

式中:T,U和V分别为系统的动能、应变能和外力势能;δ为变分符号;t1和t2分别表示系统运动的初始时刻和终止时刻。这里只考虑横向位移,令w=w(x,y,t),可得

(7)～(9)式中

分别积分后得到

式中

式中:D11,D12是抗弯刚度;D33为抗扭刚度,当n趋近于0时,抗弯刚度变为Dc,当n趋近于无穷时,抗弯刚度变为Dm,I1为多孔FGM矩形板的惯性系数,I3为转动惯量。将(7)～(9)式代入(6)式,并忽略面内和转动惯性力,可得四边受压多孔FGM矩形板横向运动的控制微分方程

(10)

(11)

式中

由弹性稳定性理论可知,当临界载荷使得系统发生屈曲时,其固有频率变为零,故当(11)式中Ω=0时,其也可以表示四边受压多孔FGM矩形板屈曲时的无量纲控制微分方程。接下来考虑多孔FGM矩形板在X=0和X=1处的无量纲边界条件,其形式如(12)～(14)式所示：

(12)

(13)

(14)

2.2 无量纲控制微分方程及边界条件的DTM变换

解决微分方程的相关问题,通常采用数值或解析的方法,这里用一种半解析法——微分变换法(DTM)对问题进行求解,其原理可参考文献[19]。

依据DTM求解原理,将四边受压多孔FGM矩形板自由振动与屈曲的无量纲控制微分方程转换为如下的代数方程

(15)

边界条件的DTM变换如下

在X=0处:

简支(S):

(16)

固定(C):

(17)

自由(F):

(18)

在X=1处:

简支(S):

(19)

固定(C):

(20)

自由(F):

(21)

(22)

结合(15)至(21)式,边界条件为两边固定两边简支(CSCS)、一边固定三边简支(CSSS)和一边固定一边自由两对边简支(CSFS)的关于无量纲固有频率的特征方程如下

(23)

(24)

如果系统发生屈曲时,无量纲固有频率Ω=0,此时关于无量纲临界载荷的求解过程类似于(23)至(24)式,可得

(25)

式中，Ncr为矩形板发生屈曲时的无量纲临界载荷。

对于边界条件为四边简支(SSSS)和一边自由三边简支(SSFS)的关于无量纲固有频率的特征方程如下

(26)

同理,要求得无量纲的固有频率和临界屈曲载荷,必须满足

(27)

(28)

对于边界条件为两边自由两边简支(FSFS)的情况下,关于无量纲固有频率的特征方程如下

(29)

同理,要求得无量纲的固有频率和临界屈曲载荷,必须满足

(30)

(31)

式中，ξ1,ξ2为迭代误差限,取ξ1=ξ2=0.000 001。

3 计算结果及分析

表1 不同边界条件下长宽比对矩形板无量纲固有频率的影响(n=0.5)

下面计算并分析四边受压多孔FGM矩形板自由振动的无量纲固有频率与临界屈曲载荷。多孔FGM矩形板材料由Si3N4/SUS304组成,材料的物性参数分别为ρc=2 370 kg/m3,ρm=7 930 kg/m3,

Ec=322 GPa,Em=196 GPa,ν=0.3。图2表示n=1,H=0.1,λ=1时,孔隙率与刚度的变化关系曲线。结果表明,随着孔隙率的增大FGM矩形板的刚度减小。

图2 孔隙率θ对多孔FGM矩形板刚度D的影响

图3 梯度指数n对多种边界条件下多孔FGM矩形板Ω的影响

图4 梯度指数n对多种边界条件下多孔FGM矩形板Ncr的影响

在CSCS、CSSS、SSSS、CSFS、SSFS和FSFS边界条件下梯度指数n与无量纲临界屈曲载荷Ncr的关系曲线。结果显示,随着n的增大,无量纲临界屈曲载荷在减小。n的取值范围从0到100,当n在较小值范围内取值时,无量纲临界载荷减小趋势比较剧烈,当n在较大值范围取值时,无量纲临界载荷变化较为平缓。当n的值一定时,边界约束越强,无量纲临界屈曲载荷也越大。

图5 孔隙率θ对多种边界条件下多孔FGM矩形板Ω的影响

图6 孔隙率θ对多种边界条件下多孔 图7 长宽比λ对多种边界条件下多孔FGM矩形板Ω的影响FGM矩形板Ncr的影响

图8 无量纲载荷N对多种边界条件下多孔FGM矩形板Ω的影响

4 结 论

本文基于经典薄板理论和Hamilton原理建立了四边受压多孔FGM矩形板的自由振动与屈曲问题的控制微分方程并无量纲化,并运用微分变化法对无量纲控制微分方程及其边界条件进行变换,通过编写MATLAB计算程序求解出无量纲固有频率和临界屈曲载荷,进一步分析了参数值的变化对无量纲固有频率和临界屈曲载荷的影响。主要结论如下:

1) 在孔隙率θ一定时,多孔FGM矩形板无量纲的固有频率Ω与临界屈曲载荷Ncr随着梯度指数n的增大而减小,随着梯度指数n趋向于无穷大,多孔矩形板的特性更加趋于多孔金属矩形板的特性。

2) 在梯度指数n一定时,多孔FGM矩形板无量纲的固有频率Ω与临界屈曲载荷Ncr随着孔隙率θ的增大而减小。

3) 在孔隙率θ与梯度指数n一定时,多孔FGM矩形板无量纲固有频率Ω随着长宽比λ的增大而增大,并且随着长宽比λ的增大,无量纲固有频率Ω增加的速度也在增大。

4) 当无量纲载荷达到无量纲临界屈曲载荷Ncr时,无量纲固有频率Ω变为零,此时多孔FGM矩形板的自由振动具有无限大的振动周期,系统发生失稳。

5) 在相同的条件下,系统的边界约束越强,其无量纲固有频率Ω和临界屈曲载荷Ncr就越大。同时,高阶无量纲固有频率Ω随孔隙率θ和梯度指数n的增大而减小的幅度大于低阶无量纲固有频率Ω。