[发明专利]一种基于子结构的复合材料弹性参数识别方法

摘要

本发明提供了一种基于子结构的复合材料弹性参数识别方法，建立复合材料子结构有限元模型，根据子结构理论对复合材料子结构模型进行动力学缩聚；缩聚后子结构特征矩阵装配到残余结构上，计算得到复合材料全模型模态信息；提取全模型模态数据，计算模态频率对残余结构弹性参数的相对灵敏度；将试验和有限元模拟的模态频率残差的二范数作为目标函数，利用迭代优化算法最小化目标函数。本发明通过考虑了子结构的复合材料建模，将模型待识别部分定义为残余结构，通过全模型模态频率对残余结构弹性参数的相对灵敏度分析和模态振型匹配，采用优化迭代算法识别复合材料待识别参数，节省计算资源，提高计算效率，具有十分重要的工程意义。

主权项

1.一种基于子结构的复合材料弹性参数识别方法，其特征在于：包括以下步骤：(1)建立复合材料子结构有限元模型，根据子结构理论对复合材料子结构模型进行动力学缩聚；(2)将缩聚后子结构特征矩阵装配到残余结构上，计算得到复合材料全模型模态信息；(3)提取全模型模态数据，计算模态频率对残余结构弹性参数的相对灵敏度；(4)将试验和有限元模拟的模态频率残差的二范数作为目标函数，利用迭代优化算法最小化目标函数；步骤(1)包括以下步骤：(1.1)建立复合材料精细化有限元模型，施加分析边界条件与载荷；将复合材料结构中材料参数已知部分定义为子结构，将材料参数未知部分定义为残余结构，分别对子结构和残余结构赋予确定性材料参数p和初始参数并将两部分有限元模型分别输出,其中为残余结构待识别的材料参数p_r的初始值；(1.2)根据所给出的材料参数p和利用有限元方法分别求解子结构和残余结构特征矩阵，并得到子结构在物理坐标u下的动力学方程：式中，M、C、K分别表示子结构的质量矩阵、阻尼矩阵和刚度矩阵,F为子结构在物理坐标下的外载荷；子结构分为内部自由度o集和边界自由度b集，将子结构动力学方程写成分块矩阵形式：式中，下角标oo和bb分别表示特征矩阵的内部分块和外部分块，ob和bo表示耦合分块；根据子结构动力学方程的分块形式，求解子结构的约束模态矩阵ψ和固定界面模态矩阵Φ_oo：其中，I_bb为与Φ_ob维度相同的单位矩阵，ω_o为固定界面特征方程对应的主模态特征频率；将约束模态矩阵和固定界面模态矩阵组装为转换矩阵H：(1.3)复合材料子结构模型动力学缩聚，得到模态坐标下的特征矩阵和动力学方程；由子结构在物理坐标u下运动方程转换到缩减的模态坐标p下的运动方程为：其中，式中，分别表示子结构在模态坐标下的质量矩阵，阻尼矩阵和刚度矩阵；p为模态坐标；为子结构在模态坐标下的外部载荷。