[发明专利]一种复合材料层合结构三维振动分析方法

摘要

本发明涉及的是一种应用于工程力学和振动工程领域的复合材料层合结构三维振动分析方法。本发明包括：提取结构的几何、材料和边界条件参数，对结构每一层沿厚度方向设置J个非均匀分布的计算平面；配置结构在计算平面上的位移；对每个计算平面上的结构位移进行面内分量展开并通过边界特征函数施加边界条件，得到计算平面上结构位移；计算第l层第j计算平面上结构的面内方向应变；求得结构在第l层任意位置的应变和应力；建立结构能量方程；求得结构特征方程。本发明通过空间分层取面把结构分解成多个计算平面，一方面降低结构维度，从而提高计算速度，节约计算成本，另一方面把结构化整为零，便于并行计算，从而提高计算效率。

主权项

1.一种复合材料层合结构三维振动分析方法，其特征是：(1)提取结构的几何、材料和边界条件参数，对结构每一层沿厚度方向设置J个非均匀分布的计算平面；其中第1个和第J个计算平面分别选取为该层的下表面和上表面；第2到第J‑1个计算平面在厚度方向上的位置选取为移位勒让德多项式PJ‑2(z)在区间[0,hl]的实零点；PJ‑2(z)表达式如下：其中，l指的是层合结构的第l层，hl为该层厚度；(2)配置结构在计算平面上的位移，即结构在第l层第j个计算平面上的结构位移设置为其中α，β为结构空间坐标系面内坐标，i＝1,2,3分别为结构位移在α，β和z方向上的分量；利用拉格朗日插值将结构任意位置的结构位移设定为其中，为步骤(1)中得到的多项式P_J‑2(z)在区间[0,h_l]的实零点；(3)利用Chebyshev多项式对每个计算平面上的结构位移进行面内分量展开并通过边界特征函数施加边界条件，得到计算平面上结构位移：其中Fi(α,β)＝(1‑α)a(1+α)b(1‑β)c(1+β)d为边界特征函数；a,b,c,d为边界特征函数系数，这些系数由边界条件确定；M，N为截断级数；(4)由步骤(2)和(3)计算第l层第j计算平面上结构的面内方向应变和横向应变和横向剪切应变和分别为：当j≠i；否则(5)由步骤(2)和(4)求得结构在第l层任意位置的应变和应力表达式为σl＝Cεl；σl＝[σlα,σlβ,σlz,τlαβ,τlαz,τlβz]T；εl＝[εlα,εlβ,εlz,γlαβ,γlαz,γlβz]T其中C为结构材料系数矩阵；(6)根据步骤(5)建立结构能量方程(U，T)：(7)在步骤(6)基础上建立结构拉格朗日能量泛函L＝U‑T，然后利用Ritz法求得结构特征方程：(K‑ω2M)＝0；其中ω为圆频率；(8)应用Arnoldi算法建立MATLAB求解器输出结构的振动特征数据并判定计算精度，若满足精度要求则输出振动特征数据的固有频率、模态，不满足则继续优化空间分层数量和增加面内位移展开式级数截取量。