王志鹏，韩志军，王龙飞

（1. 太原理工大学机械与运载工程学院，山西 太原 030024；2. 太原理工大学生物医学工程学院，山西 太原 030024）

复合材料板具有较高的强度质量比、良好的耐腐蚀性和优异的可设计性，被广泛应用于航空航天和工业制造等领域[1]。在实际使用中复合材料板经常受到不同形式的冲击荷载，从而产生振动和屈曲问题，因此冲击载荷作用下复合材料板的动力稳定性问题备受关注。

近年来关于复合材料板的研究越来越多，尤其是冲击荷载作用下复合材料板的动力稳定性问题研究[2]，对工程部件结构设计和使用具有重要的意义。Sun 等[3]研究了在加热环境中应力波对功能梯度圆柱壳轴向冲击屈曲的影响；毛柳伟等[4]对弹性直杆在应力波作用下的动力分叉屈曲进行了分析与探讨，提出了求解应力波作用下直杆动力屈曲的数值方法；Lepik[5]讨论了在应力波影响下轴向压缩的弹塑性梁的屈曲；Zhang 等[6]分析了不确定初始几何缺陷对薄板屈曲的影响；Abdelaziz 等[7]利用双曲线剪切变形理论，分析了在各种边界条件下复合材料板的弯曲变形和屈曲；Kouchakzadeh 等[8]采用线性和旋转弹簧的均匀分布来模拟边界条件，对矩形层压复合板的屈曲进行了分析；Czapski 等[9]通过数值和实验方法，研究了残余应力对压缩至破坏期间薄壁层压板屈曲性能的影响。

在实际工程中，复合材料板多应用于振动环境，其在动力响应下的动态特性和振动分析是必不可少的，因而对该类材料的振动屈曲研究至关重要。Kuo[10]研究了两种非均匀分布纤维复合材料板的振动屈曲问题，Villarreal 等[11]对典型正交异性板的本征频率和振动屈曲进行了理论分析，Eftekhari 等[12]提出了通过组合应用有限元方法和微分正交方法求解矩形板的振动屈曲问题，Rehman 等[13]探讨了壳体结构的缺陷和损坏对结构振动屈曲的影响，Sayyad 等[14]将三角剪切变形理论应用于复合板的变形和振动屈曲研究。

关于复合材料板的振动屈曲问题已开展了较多的研究，但大多未考虑应力波效应对振动屈曲的影响，而动力屈曲一般与应力波相联系且具有局部发生的特点，研究含初始缺陷的复合材料板能更好地揭示实际工程中复合材料板在不同工况下发生动力屈曲的机理。基于此，本研究利用Kirchhoff 薄板理论和Hamilton 原理，考虑应力波效应，建立含初始几何缺陷的四边简支复合材料板的振动控制方程，得到板的屈曲临界荷载表达式，在此基础上通过数值计算讨论初始几何缺陷、振型函数初相位、铺层角度、屈曲模态阶数和铺层层数对复合材料板振动屈曲临界荷载的影响，为工程实际提供理论依据。

1 复合材料板的控制方程

复合材料板在x=La处为固定边界条件，其余3 边为简支边界条件，在z= 0 的中性面上受x方向的面内阶跃荷载N作用，如图1 所示。板在z方向上含初始几何缺陷w1，且x、y、z方向的位移分别为u、v、w。复合材料板的长度、宽度和厚度分别为La、Lb、h，由n层单层板组成， θ为单层板的铺层角度，即纤维材料铺设方向与x方向的夹角（见图1）。

图1 复合材料板结构示意图Fig. 1 Schematic diagram of composite plate structure

根据Kirchhoff 薄板理论及经典弹性理论，复合材料板的位移与应变、弹性曲面的曲率和扭率的表达式为

复合材料板在左端受面内冲击荷载N作用（见图1）时，应力波沿x方向在板内传播，其应力变化如图2 所示。

图2 应力波传播示意图Fig. 2 Schematic diagram of stress wave propagation

当应力波传播至波阵面位置Lcr（临界长度）时，板发生振动屈曲，板的内力Nt和应力波波速c分别表示为

板发生振动屈曲时的变形能可以表示为

根据 Kirchhoff 薄板理论及经典弹性理论，薄板中面在变形过程中没有伸长变形，将板的本构关系代入式(16)～式(18)中，略去含u0和v0的项，得到复合材料板在面向阶跃荷载激励下的控制方程

给予对照组的120例孕妇常规孕期护理,其主要护理内容为饮食、行为、常规检查、一般治疗等。而对照组的120例孕妇,则是以常规孕期护理为基础,增添孕期健康教育,其主要内容如下。

由屈曲模态确定的缺陷分布形式是板结构最有可能发生的屈曲形式，能够很好地确定结构的缺陷敏感性[18]。对于复合材料板在制造过程中出现的初始几何缺陷，引入屈曲模态的 ε倍变形作为初始几何缺陷[19]，可以表示为

式中：i、j为屈曲模态阶数，i,j=1,2,3,···；Rij为板的第(i,j)阶模态幅值； ε表示初始几何缺陷系数。

根据式(21)和式(22)，利用棣莫弗公式对控制方程式(20)中的各项进行求导并化简，得到

2 算例分析

利用MATLAB 数值分析应力波未反射时初始几何缺陷、初相位、铺层角度、屈曲模态阶数、铺设厚度以及铺层层数对复合材料板振动屈曲临界荷载的影响，使用的材料参数见表1[21]。

表1 复合材料板参数[21]Table 1 Material parameters of composite plate[21]

图3 显示了复合材料板的模态阶数i为1、2、3，j为1 时的屈曲模态。当板的x方向模态增大时，x方向的屈曲模态第一峰值增大且波数增加，而y方向的屈曲模态呈正对称分布。模态阶数的增加使板振动的屈曲模态变得更复杂。

图3 x 方向模态取值增大时板的屈曲模态Fig. 3 Buckling mode of composite plate with increasing mode value in x direction

设置7 组算例，分别以初始几何缺陷、初相位、铺层角度、x和y两个方向屈曲模态阶数、铺层层数及铺设厚度为变量进行算例分析，研究以上因素对板振动屈曲临界荷载的影响，算例参数见表2。

表2 算例分析参数表Table 2 Example analysis parameter table

将A 组数据代入式(25)中，可以得到不同初始缺陷系数对复合材料板振动屈曲的影响，如图4所示。由Lcr-Ncr曲线可知：在应力波传播过程中，Ncr呈指数下降，分为敏感区和非敏感区。应力波在Lcr＜ 0.4 m 区域传播时，Ncr的变化较陡峭，该区域为敏感区；应力波在Lcr＞ 0.4 m 区域传播时，Ncr的变化趋于平缓，该区域为非敏感区，因此敏感分界点为0.4 m。当选取的初始缺陷系数增大时，临界荷载Ncr也随之增大。在敏感区，初始缺陷系数对临界荷载Ncr的影响较大，且随应力波传播呈减小趋势。此外，初始缺陷系数对非敏感区的影响较小。图4 表明，初始几何缺陷系数越大，板越容易发生屈曲。

图4 不同初始缺陷系数条件下Ncr 与Lcr 的关系曲线Fig. 4 Relationship between Ncr and Lcr under different initial defect coefficients

图5 x 方向模态阶数不同时Ncr 与Lcr 的关系曲线Fig. 5 Relationship between Ncr and Lcr with different order of modes in x direction

将C 组数据代入式(25)中，可以得到不同的y方向模态阶数对复合材料板振动屈曲的影响，如图6 所示。由Lcr-Ncr曲线可知：当选取的y方向模态阶数增大时，临界荷载Ncr也随之变大。在应力波传播过程中，在敏感区y方向模态阶数对临界载荷基本没有影响，而在非敏感区有极小的影响。图6 表明，y方向模态阶数的变化对板屈曲临界荷载基本没有影响。

图6 y 方向模态阶数不同时Ncr 与Lcr 的关系曲线Fig. 6 Relationship between Ncr and Lcr with different order of modes in y direction

将D 组数据代入式(25)中，得到不同铺层角度对复合材料板振动屈曲的影响，如图7 所示。由Lcr-Ncr曲线可知：在敏感区，不同的铺设角度对临界荷载Ncr的影响较大，且随应力波的传播不断减小，到达非敏感区后趋于稳定。图7 表明，铺层角度小的单层板的层数越多，板的临界屈曲荷载越大，说明复合材料板的铺设角度直接影响板的屈曲临界荷载。

图7 不同铺层角度条件下Ncr 与Lcr 的关系曲线Fig. 7 Relationship between Ncr and Lcr under different laying angles

将E 组数据代入式(25)中，得到不同初相位对复合材料板振动屈曲的影响，如图8 所示。由Lcr-Ncr曲线可知：振型函数的初相位越大，对应的临界荷载越大。在敏感区，振型函数的初相位对临界荷载Ncr的影响较小，且随应力波的传播不断减小；到达非敏感区之后，影响趋于平缓。图8 表明，振型函数的初相位越大，板的屈曲临界荷载越大。

图8 不同初相位条件下Ncr 与Lcr 的关系曲线Fig. 8 Relationship between Ncr and Lcr under the condition of initial phase of different mode functions

将F 组数据代入式(25)中，得到不同铺层层数对复合材料板振动屈曲的影响，如图9 所示。由Lcr-Ncr曲线可知：当按照不同层数铺设时，敏感区的临界荷载Ncr的变化较大，且随应力波的传播不断减小；到达非敏感区之后变化较小并趋于平缓。图9 表明，对于厚度固定且对称铺设的板，当铺设层数达到7 时，其屈曲荷载随层数增加趋于稳定。

图9 不同铺层层数下Ncr 与Lcr 的关系曲线Fig. 9 Relationship between Ncr and Lcr under different laying modes

将G 组数据代入式(25)中，得到不同铺设厚度对复合材料板振动屈曲的影响，如图10 所示。由Lcr-Ncr曲线可知：板的铺设厚度越大，对应的临界荷载越大。在敏感区，不同的板厚对临界荷载Ncr的影响很大，且随应力波的传播不断减小，到达非敏感区后趋于稳定。图10 表明，复合材料板的铺设厚度将直接决定板的屈曲临界荷载。

图10 不同铺设厚度下Ncr 与Lcr 的关系曲线Fig. 10 Relationship between Ncr and Lcr under different thicknesses

3 结 论

基于Kirchhoff 薄板理论和Hamilton 变分原理，建立了具有初始几何缺陷的四边简支复合材料板的振动控制方程。采用伽辽金法，选取符合边界条件的振型函数求解控制方程，得到屈曲临界载荷表达式。数值计算结果表明：应力波在未发生反射前的传播过程中，复合材料板的振动屈曲临界载荷随着临界长度的增大、铺设厚度的减小、初始几何缺陷系数的增大、振型函数初相位的减小而减小；复合材料板的各层铺层角度与荷载作用方向的夹角越小，屈曲临界载荷越大，当对称铺设层数达7 层时，临界荷载趋于稳定。研究结果可为工程中复合材料板的结构设计与应用提供一定的参考。