考虑前屈曲耦合变形时功能梯度圆板的稳定性分析

Stability Analysis of Functionally Graded Circular Plates Considering Pre-Buckling Coupling Deformation

摘要: 功能梯度材料(FGM)结构的非均匀性来自于材料组份沿某一方向梯度变化，因此会存在拉弯耦合效应。在周边简支边界条件下，只要有面内载荷就会产生挠度。本文基于经典板理论，研究了机械载荷和热载荷作用下FGM简支圆板的非线性变形及稳定性问题。假设FGM圆板的材料性质只沿厚度方向进行变化，利用能量原理推导出FGM圆板的平衡方程，得到了包含前屈曲耦合变形影响的控制方程，并用打靶法进行求解。讨论了外载荷、前屈曲耦合变形以及材料的温度依赖性质等因素对FGM圆板非线性变形和稳定性的影响。

Abstract: The inhomogeneity of functionally graded materials (FGM) structure comes from the gradient change of material components along a certain direction, so there will be stretch bending coupling effect. As long as the in-plane load is applied, deflection will occur, in the case of simply supported boundary conditions. Based on the classical plate theory, the nonlinear deformation and stability of FGM simply supported circular plates under mechanical and thermal loads are studied in this paper. Assuming that the material properties of FGM circular plate only change along the thickness direction, the equilibrium equation of FGM circular plate is derived by using the energy principle, and the governing equation including the influence of pre-buckling coupling deformation is obtained, which is solved by shooting method. The effects of external loads, pre-buckling coupling deformation and temperature dependence of material properties on the nonlinear deformation and stability of FGM circular plate are discussed.

1. 引言

功能梯度材料是一种材料组分按梯度变化的新型材料，因其优异性能被广泛关注。关于FGM结构稳定性分析的研究结果已经十分丰富。基于Love-Kirchhoff假设和Sander非线性应变位移关系，Najafizadeh和Eslami [1] 研究了受均匀径向压力的简支和夹紧FGM圆板的屈曲问题。几年后，Najafizadeh和Heydari [2] 又给出了基于高阶剪切变形理论的封闭解。Ma和Wang [3] [4] 分析了FGM圆板的弯曲和过屈曲，发现FGM的过屈曲行为不同于均匀板。对于具有任意材料特性的FGM板的后屈曲行为，Lal等 [5] 得到了后屈曲响应的二阶统计量。而Li等 [6] 通过打靶法数值求解，直接得到了有几何缺陷的FGM圆板的非线性热–机械载荷后屈曲响应。Samsam等 [7] 和Mohammadi等 [8] 研究了FGM矩形薄板的屈曲问题。利用有限元方法，Ghomshei和Abbasi [9] 分析了变厚度FGM圆环板的热屈曲。后来，李清禄等 [10] 研究了变厚度FGM圆板的热后屈曲，结果表明径向厚度的变化并不会影响热临界载荷值，但会影响屈曲以后的平衡路径。而Zhao [11] 发现，当热载荷和机械载荷的值以及它们的加载顺序满足一定的组合时，固支圆板不可避免地会发生二次屈曲(跳跃屈曲)。最近，Trinh等 [12] 给出了多孔FGM板随机屈曲的封闭形式的临界载荷，也有一些研究人员 [13] [14] 分析了裂纹对FGM板稳定性的影响。

FGM结构由于具有非均匀性，在简支边界条件下会引起前屈曲耦合变形，即只要面内载荷存在，就产生耦合挠度。这种情形下，FGM简支圆板的屈曲问题变得十分复杂，与均匀材料圆板有着显著的区别 [3] [4] [15] [16]。本文中屈曲是指面内载荷有微小变化就会引起圆板挠度的大幅度增加或减小，在其载荷–挠度曲线上出现一个“平台” [15] [16] [17]。杨帆和马连生 [18] 的研究结果表明，前屈曲耦合变形对FGM圆板的稳定性确实有一定影响。而过去的一些研究成果没有考虑前屈曲耦合变形，掉进了Quta和Leissa [15] [16] 所称的“陷阱”，这个观点也被很多研究人员的成果 [3] [19] [20] [21] [22] [23] 所证实。

据作者所知，关于前屈曲耦合变形对FGM结构稳定性影响的研究成果很少见诸于报道。本文将在文献 [18] 的基础上，进一步研究面内载荷作用下FGM简支圆板的稳定性问题，讨论前屈曲耦合变形、外载荷和材料的温度依赖性质等因素对FGM圆板稳定性的影响。

2. 基本方程

研究一个半径为b，厚度为h的FGM圆板，设板承受均布径向压力p或受可沿板厚度方向变化的热场作用。假设材料组分只沿厚度方向变化，并且材料性质P(弹性模量E、热膨胀系数α及热传导率K等)可表示为：

$P\left(z\right)=\left(P_c-P_m\right){\left(\frac{h+2z}{2h}\right)}^n+P_m$ (1)

这里，下标c和m分别表示陶瓷和金属材料，上标n是梯度指数。

具有温度依赖的材料性质P写为：

$P=P_0\left(P_{-1}{T}^{-1}+1+P_{1}T+P_2{T}^2+P_3{T}^3\right)$ (2)

考虑以下一个稳态温度场，设该温度场仅沿板的厚度方向变化：

$-\frac{\text{d}}{\text{d}z}\left(K\left(z\right)\frac{\text{d}T\left(z\right)}{\text{d}z}\right)=0$ (3)

式中， $T\left(h/2\right)=T_c$， $T\left(-h/2\right)=T_m$，积分上式得到温度场：

$T\left(z\right)=T_m+\left(T_c-T_m\right){\displaystyle {\int }_{-h/2}^z\frac{\text{d}z}{K\left(z\right)}}/{\displaystyle {\int }_{-h/2}^{h/2}\frac{\text{d}z}{K\left(z\right)}}$ (4)

经典板理论下的位移场：

$U_r\left(r,z\right)=U\left(r\right)-z{W}_{,r},U_z\left(r,z\right)=W\left(r\right)$ (5)

这里，U和W分别是板中面的径向位移和横向位移。

板的应变与位移关系：

${\epsilon }_r=U_{,r}+\frac{{\left(W_{,r}\right)}^2}{2}-z{W}_{,rr}={\epsilon }_r^0+z{k}_r,{\epsilon }_{\theta }=\frac{U}{r}-\frac{z{W}_{,r}}{r}={\epsilon }_{\theta }^0+z{k}_{\theta }$ (6)

内力位移关系为：

$N_r=A_{11}\left(\frac{\partial U}{\partial r}+\nu \frac{U}{r}+\frac{1}{2}{\left(\frac{\partial W}{\partial r}\right)}^2\right)-B_{11}\left(\frac{{\partial }^{2}W}{\partial r^2}+\frac{\nu }{r}\frac{\partial W}{\partial r}\right)-N^T$ (7a)

$N_{\theta }=A_{11}\left(\nu \frac{\partial U}{\partial r}+\frac{U}{r}+\frac{\nu }{2}{\left(\frac{\partial W}{\partial r}\right)}^2\right)-B_{11}\left(\nu \frac{{\partial }^{2}W}{\partial r^2}+\frac{1}{r}\frac{\partial W}{\partial r}\right)-N^T$ (7b)

$M_r=B_{11}\left(\frac{\partial U}{\partial r}+\nu \frac{U}{r}+\frac{1}{2}{\left(\frac{\partial W}{\partial r}\right)}^2\right)-D_{11}\left(\frac{{\partial }^{2}W}{\partial r^2}+\frac{\nu }{r}\frac{\partial W}{\partial r}\right)-M^T$ (7c)

$M_{\theta }=B_{11}\left(\nu \frac{\partial U}{\partial r}+\frac{U}{r}+\frac{\nu }{2}{\left(\frac{\partial W}{\partial r}\right)}^2\right)-D_{11}\left(\nu \frac{{\partial }^{2}W}{\partial r^2}+\frac{1}{r}\frac{\partial W}{\partial r}\right)-M^T$ (7d)

其中， $A_{11}$， $B_{11}$ 和 $D_{11}$ 分别是拉压刚度、耦合刚度和抗弯刚度， $N^T$ 和 $M^T$ 分别是热薄膜力和热弯矩：

$\left(A_{11},B_{11},D_{11}\right)={\displaystyle {\int }_{-h/2}^{h/2}{Q}_{11}\left(1,z,z^2\right)\text{d}z},Q_{11}=\frac{E\left(z\right)}{1-{\nu }^2}$ (8a)

$\left(N^T,M^T\right)={\displaystyle {\int }_{-h/2}^{h/2}\frac{E\left(z\right)}{1-\nu }\alpha \left(z\right)T\left(z\right)\left(1,z\right)\text{d}z}$ (8b)

根据能量原理，可以推出FGM圆板的平衡方程与边界条件，其无量纲形式如下：

$\frac{\text{d}}{\text{d}x}\left[\frac{1}{x}\frac{\text{d}}{\text{d}x}\left(xu\right)\right]+\left[\frac{\text{d}w}{\text{d}x}\frac{{\text{d}}^{2}w}{\text{d}{x}^2}+\frac{1-\nu }{2x}{\left(\frac{\text{d}w}{\text{d}x}\right)}^2\right]-F_1\frac{\text{d}}{\text{d}x}{\nabla }^{2}w=0$ (9a)

$\begin{array}{l}{\nabla }^{4}w-F_2\left[\nu \frac{\text{d}u}{\text{d}x}+\frac{u}{x}+\frac{\nu }{2}{\left(\frac{\text{d}w}{\text{d}x}\right)}^2\right]\frac{1}{x}\frac{\text{d}w}{\text{d}x}-F_2\left[\frac{\text{d}u}{\text{d}x}+\nu \frac{u}{x}+\frac{1}{2}{\left(\frac{\text{d}w}{\text{d}x}\right)}^2\right]\frac{{\text{d}}^{2}w}{\text{d}{x}^2}\\ +F_3\left[{\left(\frac{{\text{d}}^{2}w}{\text{d}{x}^2}\right)}^2+\frac{2\nu }{x}\frac{\text{d}w}{\text{d}x}\frac{{\text{d}}^{2}w}{\text{d}{x}^2}+\frac{1}{x^2}{\left(\frac{\text{d}w}{\text{d}x}\right)}^2\right]+\bar{N}{\nabla }^{2}w=0\end{array}$ (9b)

这里， ${\nabla }^2=\frac{1}{x}\frac{\text{d}}{\text{d}x}\left(x\frac{\text{d}}{\text{d}x}\right)$， ${\nabla }^4=\frac{1}{x}\frac{\text{d}}{\text{d}x}\left[x\frac{\text{d}}{\text{d}x}\left(\frac{1}{x}\frac{\text{d}}{\text{d}x}\left(x\frac{\text{d}}{\text{d}x}\right)\right)\right]$。

板中心处条件：

$u=\frac{\text{d}w}{\text{d}x}=0,\underset{x\to 0}{\lim }\left(\frac{{\text{d}}^{3}w}{\text{d}{x}^3}+\frac{1}{x}\frac{{\text{d}}^{2}w}{\text{d}{x}^2}\right)=0$ (10)

周边简支边界条件：

$u=w=0,F_4\left[\frac{\text{d}u}{\text{d}x}+\frac{1}{2}{\left(\frac{\text{d}w}{\text{d}x}\right)}^2\right]-\left(\frac{{\text{d}}^{2}w}{\text{d}{x}^2}+\frac{\nu }{x}\frac{\text{d}w}{\text{d}x}\right)-\bar{M}=0$ (11)

无量纲量如下：

$\begin{array}{l}x=\frac{r}{b},w=\frac{W}{h},u=\frac{Ub}{h^2},F_1=\frac{F_3}{F_2},F_2=\frac{A_{11}{h}^2}{\Omega },F_3=\frac{B_{11}h}{\Omega },F_4=\frac{B_{11}h}{D_{11}},\\ \bar{N}=\frac{N^T{b}^2}{\Omega },\bar{M}=\frac{M^T{b}^2}{D_{11}h},\Omega =D_{11}-\frac{B_{11}^2}{A_{11}},F_p=\frac{p{b}^2}{A_{11}{h}^2}\end{array}$

求解上述问题，可以得到FGM简支圆板的载荷–挠度曲线(图1，图3，图5)。从这些图中可以发现，在某些情况下(如图5)，横向挠度随着热载荷的增加而增大，不出现“平台”。从图1和图3均可以看到，对于任一条载荷–挠度曲线，在某个载荷数值附近，只要载荷数值小幅增加，挠度数值就会大幅增大。即外载荷的微小增加导致了圆板横向挠度的迅速增大，这意味着相应FGM圆板的抗弯刚度急剧降低，近乎于丧失对变形的抵抗能力。因此，当这样一个“平台”出现时，可以认为圆板发生了屈曲 [17] [18]。

下面分析前屈曲耦合变形对FGM简支圆板屈曲的影响。

设在外因素F作用下，FGM圆板的位移为 $\left(\bar{u},\bar{w}\right)$，当外因素增加∆F后，位移变为：

$u=\bar{u}+\Delta u,w=\bar{w}+\Delta w$ (12)

将式(12)代入式(9)，并忽略掉高阶小量后，得到关于∆u，∆w的线性方程：

$\begin{array}{l}{\nabla }^4\Delta w-F_2\left[\nu \frac{\text{d}\Delta u}{\text{d}x}+\frac{\Delta u}{x}+\nu \frac{\text{d}\Delta w}{\text{d}x}\frac{\text{d}\bar{w}}{\text{d}x}\right]\frac{1}{x}\frac{\text{d}\bar{w}}{\text{d}x}-F_2\left[\nu \frac{\text{d}\bar{u}}{\text{d}x}+\frac{\bar{u}}{x}+\frac{\nu }{2}{\left(\frac{\text{d}\bar{w}}{\text{d}x}\right)}^2\right]\frac{1}{x}\frac{\text{d}\Delta w}{\text{d}x}\\ -F_2\left[\frac{\text{d}\Delta u}{\text{d}x}+\frac{\nu }{x}\Delta u+\frac{\text{d}\bar{w}}{\text{d}x}\frac{\text{d}\Delta w}{\text{d}x}\right]\frac{{\text{d}}^2\bar{w}}{\text{d}{x}^2}-F_2\left[\frac{\text{d}\bar{u}}{\text{d}x}+\frac{\nu }{x}\bar{u}+\frac{1}{2}{\left(\frac{\text{d}\bar{w}}{\text{d}x}\right)}^2\right]\frac{{\text{d}}^2\Delta w}{\text{d}{x}^2}\\ +F_3\left[\left(2\frac{{\text{d}}^2\bar{w}}{\text{d}{x}^2}\frac{{\text{d}}^2\Delta w}{\text{d}{x}^2}\right)+\frac{2\nu }{x}\left(\frac{\text{d}\Delta w}{\text{d}x}\frac{{\text{d}}^2\bar{w}}{\text{d}{x}^2}+\frac{\text{d}\bar{w}}{\text{d}x}\frac{{\text{d}}^2\Delta w}{\text{d}{x}^2}\right)+\frac{2}{x^2}\frac{\text{d}\bar{w}}{\text{d}x}\frac{\text{d}\Delta w}{\text{d}x}\right]+\bar{N}{\nabla }^2\Delta w=0\end{array}$ (13a)

$\frac{\text{d}}{\text{d}x}\left[\frac{1}{x}\frac{\text{d}}{\text{d}x}\left(x\Delta u\right)\right]+\left[\frac{{\text{d}}^2\bar{w}}{\text{d}{x}^2}\frac{\text{d}\Delta w}{\text{d}x}+\frac{\text{d}\bar{w}}{\text{d}x}\frac{{\text{d}}^2\Delta w}{\text{d}{x}^2}+\frac{1-\nu }{2x}\frac{\text{d}\bar{w}}{\text{d}x}\frac{\text{d}\Delta w}{\text{d}x}\right]-F_1\frac{\text{d}}{\text{d}x}{\nabla }^2\Delta w=0$ (13b)

结合具体的边界条件，求出前屈曲耦合变形 $\bar{u}$ 和 $\bar{u}$ (其中含有待定的外载荷)。然后将其代入式(13)，并在关于∆u和∆w的齐次边界条件下，利用∆u和∆w有非零解这一条件来确定临界屈曲载荷。显然，所求出的临界屈曲载荷包含了满足式(9)的前屈曲耦合变形的影响。如果不存在非零解，则说明该问题并无屈曲现象发生。

当仅考虑均匀径向压力时，上式中令 $\bar{N}=0$ 即可。

3. 数值结果与讨论

本文的分析以Si₃N₄和SUS304为例，取泊松比 $\nu =0.3$，材料性质参数值见表1。

Table 1. Temperature-dependent coefficients for ceramic (Si3N4) and metals (SUS304)

表1. 陶瓷和金属的材料常数

3.1. 受机械载荷作用

我们先分析受均匀径向压力p的FGM简支圆板，边界条件为：

$\begin{array}{l}w=0,F_4\left[\frac{\text{d}u}{\text{d}x}+\nu \frac{u}{x}+\frac{1}{2}{\left(\frac{\text{d}w}{\text{d}x}\right)}^2\right]-\left(\frac{{\text{d}}^{2}w}{\text{d}{x}^2}+\frac{\nu }{x}\frac{\text{d}w}{\text{d}x}\right)=0\\ \left[\frac{\text{d}u}{\text{d}x}+\nu \frac{u}{x}+\frac{1}{2}{\left(\frac{\text{d}w}{\text{d}x}\right)}^2\right]-F_1\left(\frac{{\text{d}}^{2}w}{\text{d}{x}^2}+\frac{\nu }{x}\frac{\text{d}w}{\text{d}x}\right)=-F_p\end{array}$ (14)

图1是均匀径向压力作用下的FGM简支可移圆板的载荷–挠度曲线。图中，1表示考虑前屈曲耦合变形，0表示未考虑前屈曲耦合变形。可以看出，从载荷作用之初，就有挠度产生，但在载荷–挠度曲线上存在quta所称的“平台” [15] [16] [17]。数值结果也表明，此种边界条件下FGM圆板的屈曲前变形是微小变形。因此，我们可以将式(9)进行线性化：

$\begin{array}{l}\frac{\text{d}}{\text{d}x}\left[\frac{1}{x}\frac{\text{d}}{\text{d}x}\left(x\bar{u}\right)\right]-F_1\frac{\text{d}}{\text{d}x}{\nabla }^2\bar{w}=0\\ {\nabla }^4\bar{w}=0\end{array}$ (15a-b)

结合边界条件式(14)解得：

$\bar{u}=-\lambda u_{0}x,\bar{w}=\lambda w_0\left(1-x^2\right)$ (16a-b)

将式(16)代入式(13)，可得包含前屈曲耦合变形影响的控制方程：

$\begin{array}{l}\frac{{\text{d}}^2\Delta u}{\text{d}{x}^2}+\frac{1}{x}\frac{\text{d}\Delta u}{\text{d}x}-\frac{\Delta u}{x^2}-c_{n1}\frac{{\text{d}}^3\Delta w}{\text{d}{x}^3}-c_{n1}\left(1+\frac{c_{n2}}{1+\nu }\lambda x^2\right)\frac{1}{x}\frac{{\text{d}}^2\Delta w}{\text{d}{x}^2}\\ +c_{n1}\left(1-c_{n2}\frac{2-\nu }{1+\nu }\lambda x^2\right)\frac{1}{x^2}\frac{\text{d}\Delta w}{\text{d}x}=0\end{array}$ (17a)

$\begin{array}{l}\frac{{\text{d}}^4\Delta w}{\text{d}{x}^4}+\frac{2}{x}\frac{{\text{d}}^3\Delta w}{\text{d}{x}^3}-\left[1-c_{n2}\left(c_{n3}-2c_{n4}\left(1+\nu \right)\right)\lambda x^2+\frac{c_{n2}^2{c}_{n4}}{2\left(1+\nu \right)}{\lambda }^2{x}^4\right]\frac{1}{x^2}\frac{{\text{d}}^2\Delta w}{\text{d}{x}^2}\\ +\left[1+c_{n2}\left(c_{n3}-2c_{n4}\left(1+\nu \right)\right)\lambda x^2-\frac{c_{n2}^2{c}_{n4}\left(2+3\nu \right)}{2\left(1+\nu \right)}{\lambda }^2{x}^4\right]\frac{1}{x^3}\frac{\text{d}\Delta w}{\text{d}x}+c_{n2}{c}_{n5}\lambda \left(\frac{\text{d}\Delta u}{\text{d}x}+\frac{\Delta u}{x}\right)=0\end{array}$ (17b)

这里，

$\begin{array}{l}{D}_c=\frac{E_c{h}^3}{12\left(1-{\nu }^2\right)},\lambda =\frac{p{b}^2}{D_c},c_{n1}=\frac{B_{11}}{A_{11}h},c_{n2}=\frac{D_c}{\Omega },u_0=\frac{D_{11}{c}_{n2}}{A_{11}{h}^2\left(1+\nu \right)},\\ w_0=\frac{F_1{c}_{n2}}{2\left(1+\nu \right)},c_{n3}=\frac{D_{11}}{\Omega },c_{n4}=\frac{B_{11}^2}{A_{11}\Omega \left(1+\nu \right)},c_{n5}=\frac{B_{11}h}{\Omega }\end{array}$

利用打靶法求解上式即可。

Figure 1. Load vs. deflection under mechanical load

图1. 机械载荷下载荷–挠度曲线

Figure 2. Critical load vs. gradient under mechanical load

图2. 机械载荷下临界载荷随n变化曲线

图2是临界载荷随梯度指数变化的曲线，虚线参见文献 [24] 中的公式。可以看出，临界载荷随着n的增加而减小。这是因为n越大，金属的组分越多，陶瓷的组分越少，从而导致抗弯刚度减小，这与文献 [25] 的结论一致。并且，前屈曲耦合变形使得FGM圆板的临界载荷增大。

3.2. 受热载荷作用

考虑两种边界条件：

简支不可移：

$u=w=0,F_4\left[\frac{\text{d}u}{\text{d}x}+\frac{1}{2}{\left(\frac{\text{d}w}{\text{d}x}\right)}^2\right]-\left(\frac{{\text{d}}^{2}w}{\text{d}{x}^2}+\frac{\nu }{x}\frac{\text{d}w}{\text{d}x}\right)-\bar{M}=0$ $x=1$ (18)

简支可移：

$\begin{array}{l}w=0,\left[\frac{\text{d}u}{\text{d}x}+\nu \frac{u}{x}+\frac{1}{2}{\left(\frac{\text{d}w}{\text{d}x}\right)}^2\right]-F_1\left(\frac{{\text{d}}^{2}w}{\text{d}{x}^2}+\frac{\nu }{x}\frac{\text{d}w}{\text{d}x}\right)-\overline{N_0}=0,\\ F_4\left[\frac{\text{d}u}{\text{d}x}+\nu \frac{u}{x}+\frac{1}{2}{\left(\frac{\text{d}w}{\text{d}x}\right)}^2\right]-\left(\frac{{\text{d}}^{2}w}{\text{d}{x}^2}+\frac{\nu }{x}\frac{\text{d}w}{\text{d}x}\right)-\bar{M}=0\end{array}$ $x=1$(19)

这里， $\overline{N_0}=\frac{N^T{b}^2}{A_{11}{h}^2}$。结合边界条件，利用打靶法可以得到热载荷作用下的载荷–挠度曲线。

图3是热载荷作用下，简支不可移FGM圆板的载荷挠度曲线，这里， $T_r=T_m/T_c$。实线表示考虑材料性质温度依赖性(TD)，虚线表示未考虑材料性质温度依赖性(TID)。从图可知，FGM圆板在同一载荷下可能存在三种构形。受热载荷的简支不可移圆板的载荷–挠度曲线与夹紧圆板的载荷–挠度曲线 [3] 有明显的不同，夹紧边界的圆板发生典型的分支屈曲，而简支不可移的圆板出现了前屈曲耦合变形，这与受均匀径向压力的简支可移圆板的变形很相似，都出现了“平台” [17]。但随着圆板两侧的温差增大，即 $T_r$ 越来越大或越来越小，“平台”越来越不明显。另外，在图3(c)中，n = 100的曲线与其它的曲线是不一样的。结合图3(a)，图3(b)，图3(d)，我们发现，这是随着温差 $T_r$ 变化，FGM圆板的构形在发生翻转，且均温时，n = 100的FGM圆板中陶瓷的材料组份非常大，很难发生翻转。再者，当均温时两个分支最接近，温差越大，两个分支距离越远。

(a) T_r = 0.1 (b) T_r = 0.5 (c) T_r = 1 (d) T_r = 5

Figure 3. Load deflection curve of FGM simply supported and immovable circular plate under non-uniform thermal loads for different values of temperature ratio T_r

图3. 非均匀热载荷FGM简支不可移圆板载荷–挠度曲线

在均温或温差较小时，前屈曲耦合挠度很明显，且属于小变形，因此，我们可以将式(9)线性化：

$\frac{\text{d}}{\text{d}x}\left[\frac{1}{x}\frac{\text{d}}{\text{d}x}\left(x\bar{u}\right)\right]-F_1\frac{\text{d}}{\text{d}x}{\nabla }^2\bar{w}=0$ (20a)

${\nabla }^4\bar{w}+\bar{N}{\nabla }^2\bar{w}=0$ (20b)

结合边界条件式(18)解得前屈曲耦合变形为：

$\bar{u}=F_1\beta A_0\left(J_1\left(\beta \right)x-J_1\left(\beta x\right)\right),\bar{w}=A_0\left(J_0\left(\beta x\right)-J_0\left(\beta \right)\right)$ (21)

其中， $\beta =\sqrt{\bar{N}},A_0=\frac{\bar{M}}{H},H=\left(1-F_1{F}_4\right)\beta J_0\left(\beta \right)+\beta J_1\left(\beta \right)\left(\nu -1+2F_1{F}_4\right)$，并且 $J_0$ 和 $J_1$ 分别为零阶第一类Bessel函数和一阶第一类Bessel函数。

将式(21)代入式(13)，可得热载荷作用下受前屈曲耦合变形影响的控制方程：

$\begin{array}{l}\frac{{\text{d}}^2\Delta u}{\text{d}{x}^2}+\frac{1}{x}\frac{\text{d}\Delta u}{\text{d}x}-\frac{\Delta u}{x^2}-F_1\frac{{\text{d}}^3\Delta w}{\text{d}{x}^3}-\left(F_1+x{A}_0\beta J_1\left(\beta x\right)\right)\frac{1}{x}\frac{{\text{d}}^2\Delta w}{\text{d}{x}^2}\\ +\left[F_1-x^2{A}_0\beta \left(J_0\left(\beta x\right)-\frac{\nu }{x}{J}_1\left(\beta x\right)\right)\right]\frac{1}{x^2}\frac{\text{d}\Delta w}{\text{d}x}=0\end{array}$ (22a)

$\begin{array}{l}\frac{{\text{d}}^4\Delta w}{\text{d}{x}^4}+\frac{2}{x}\frac{{\text{d}}^3\Delta w}{d{x}^3}-\left[1+x\beta A_0{J}_1\left(\beta x\right)\left(F_1{F}_2-2F_3\right)\left(1-\nu \right)+\frac{1}{2}{x}^2{A}_0^2{\beta }^2{\left(J_1\left(\beta x\right)\right)}^2\right]\frac{1}{x^2}\frac{{\text{d}}^2\Delta w}{\text{d}{x}^2}\\ -\left[-x^2\beta A_0\left(J_0\left(\beta x\right)\left(F_1{F}_2-2F_3\right)-\left(1+\nu \right)F_1{F}_2{J}_1\left(\beta \right)\right)\right]\frac{1}{x^2}\frac{{\text{d}}^2\Delta w}{\text{d}{x}^2}\\ +\left[1+x\beta A_0{J}_1\left(\beta x\right)\left(F_1{F}_2-2F_3\right)\left(1-\nu \right)+\frac{3\nu }{2}{F}_2{x}^2{A}_0^2{\beta }^2{\left(J_1\left(\beta x\right)\right)}^2\right]\frac{1}{x^3}\frac{\text{d}\Delta w}{\text{d}x}\\ +\left[-x^2\left(\left(1+\nu \right)\beta A_0{F}_1{F}_2{J}_1\left(\beta \right)-\nu \beta A_0{J}_0\left(\beta x\right)\left(F_1{F}_2-2F_3\right)\right)\right]\frac{1}{x^3}\frac{\text{d}\Delta w}{\text{d}x}\\ -\left[F_2{A}_0^2{\beta }^2{J}_1\left(\beta x\right)\left(J_0\left(\beta x\right)-\frac{1}{x}{J}_1\left(\beta x\right)\right)\right]\frac{\text{d}\Delta w}{\text{d}x}+\bar{N}\left(\frac{{\text{d}}^2\Delta w}{\text{d}{x}^2}+\frac{1}{x}\frac{\text{d}\Delta w}{\text{d}x}\right)\\ +\frac{F_2}{x}\beta A_0{J}_1\left(\beta x\right)\left(\nu \frac{\text{d}\Delta u}{\text{d}x}+\frac{\Delta u}{x}\right)+F_2\beta A_0\left[J_0\left(\beta x\right)-\frac{1}{x}{J}_1\left(\beta x\right)\right]\left(\frac{\text{d}\Delta u}{\text{d}x}+\frac{\nu }{x}\Delta u\right)=0\end{array}$ (22b)

利用打靶法求解即可得到热载荷作用下受前屈曲耦合变形影响的临界载荷。

图4是临界载荷随梯度指数变化曲线，不考虑前屈曲耦合变形的数据参见文献 [24] 中的公式。数据表明，前屈曲耦合变形对稳定性是有影响的。非均温时，前屈曲耦合变形使得临界载荷增大，而考虑材料性质温度依赖性会让临界载荷变小，但依然大于不考虑前屈曲耦合变形时的临界载荷。均温时，前屈曲耦合变形只在局部范围内对临界载荷有明显影响，当n > 20后，几乎没有影响。随着 $T_r$ 逐渐增大，即陶瓷侧温度越来越低，金属侧温度越来越高，FGM圆板的临界载荷也在逐渐减小。

图5是热载荷作用下简支可移FGM圆板的载荷–挠度曲线。这是热弯曲现象，相同热载作用下，考虑温度依赖性时产生的挠度更小，这可能是因为温度依赖性减弱了拉弯耦合效应对挠度的影响。

(a) T_r = 0.1 (b) T_r = 0.5 (c) T_r = 1 (d) T_r = 5

Figure 4. Critical load vs. gradient of FGM simply supported and immovable circular plate under non-uniform thermal load for different values of temperature ratio T_r

图4. 非均匀热载荷下临界载荷随梯度指数n变化曲线

Figure 5. Load deflection curve of FGM simply supported and movable circular plate under uniform thermal load

图5. 均温下FGM简支可移圆板载荷–挠度曲线

4. 结论

本文基于经典板理论，研究了FGM简支圆板的非线性变形以及稳定性问题。推导了机械载荷和热载荷作用下，FGM简支圆板受前屈曲耦合变形影响的控制方程，分析了前屈曲耦合变形对FGM圆板稳定性的影响，结果表明：

1) 机械载荷作用下，考虑前屈曲耦合变形时，FGM简支圆板的临界载荷增大。

2) 非均匀热载荷作用下，简支不可移边界下FGM圆板的前屈曲耦合变形使得临界载荷增大；而均温时，前屈曲耦合变形只在局部范围内对临界载荷有明显影响。

3) 材料的温度依赖性质降低了FGM圆板的临界载荷。

4) 均温下的简支可移FGM圆板的挠度随着热载荷的增加而增大，不出现“平台”，只是通常意义下的弯曲。

5) 陶瓷侧温度越来越低，金属侧温度越来越高时，FGM圆板的临界载荷越小。

6) 同一热载荷下，FGM圆板可能存在三种构形；随着上下面的温差变化，FGM圆板的构形会发生翻转。

基金项目

本文得到国家自然科学基金11862012和12062010的资助。