CN107292045A - 一种含阻尼结构有限元模型的复模态评估方法 - Google Patents
一种含阻尼结构有限元模型的复模态评估方法 Download PDFInfo
- Publication number
- CN107292045A CN107292045A CN201710521580.2A CN201710521580A CN107292045A CN 107292045 A CN107292045 A CN 107292045A CN 201710521580 A CN201710521580 A CN 201710521580A CN 107292045 A CN107292045 A CN 107292045A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- complex
- vibration type
- mtd
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000013016 damping Methods 0.000 title claims abstract description 97
- 238000000034 method Methods 0.000 title claims abstract description 60
- 239000011159 matrix material Substances 0.000 claims abstract description 37
- 238000011156 evaluation Methods 0.000 claims description 75
- 238000012360 testing method Methods 0.000 claims description 46
- 238000010219 correlation analysis Methods 0.000 claims description 29
- 238000004088 simulation Methods 0.000 claims description 26
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000004458 analytical method Methods 0.000 claims description 7
- 238000005316 response function Methods 0.000 claims description 4
- 230000005484 gravity Effects 0.000 claims description 3
- 238000000205 computational method Methods 0.000 abstract 1
- 230000007547 defect Effects 0.000 description 5
- 238000012854 evaluation process Methods 0.000 description 5
- 238000013210 evaluation model Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Complex Calculations (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种含阻尼结构有限元模型的复模态评估方法，首先根据试验复频率和仿真复频率的相对误差求出复频率相关性系数，并通过模态置信度准则求出试验和仿真的复振型幅值相关性系数、通过定义的矩阵范数求出复振型相位相关性系数；然后使用梯形隶属函数，将复频率、复振型幅值及复振型相位相关性系数模糊化，得到该三类影响因素的模糊关系矩阵，将模糊关系矩阵及各阶模态权重集加权运算得到各影响因素的评估结果；最后将各评估结果及三类影响因素的权重集进行加权运算得到一个综合评估值。本发明考虑了阻尼对结构响应的影响，适用于含阻尼结构模型的修正、误差位置识别及预估模型的动响应准确性，其计算方法实施规范，便于用计算机编程实施。
Description
技术领域
本发明属于结构动力学建模技术领域，涉及含阻尼模型的修正、误差位置识别及预估模型响应技术，尤其是一种含阻尼结构有限元模型的复模态评估方法。
背景技术
模型评估是衡量有限元建模准确度的必要手段，准确的有限元模型可为结构的动特性分析提供极大的便利。目前常用的评估方法是从实频率和实振型两个单方面来评价有限元模型，但评估结果不能有效用于阻尼结构模型响应计算精度的预估。一般常用的结构动力学有限元模型评估方法是求频率相对误差和模态置信因子(MAC)，具体的方法是从试验所测得复模态中提取实模态的相关信息，然后与所建立的有限元模型计算结果进行频率和振型的相关性分析，频率相对误差用于评价模态频率的准确性，而模态置信因子用于评价模态振型的准确性。
但现有技术常用的经典方法即使评估的结果良好，模型用于响应计算可能仍具有较大误差。其产生误差的原因主要在于三个方面，一是经典方法评估结果多，难以整体把握模型；二是振型MAC评估不能反映阻尼引起的各个自由度上相位角的差异；三是在评估的过程中需要将试验测得的复模态实数化，会产生二次误差。“动力学有限元模型的模糊评估”(张安平,王轲.动力学有限元模型的模糊评估[J].南京航空航天大学学报,2006,(03):367-372.)提出了一种有限元模型模糊评估方法，该方法通过取不同的权重将频率和振型评估结果进行综合，用一个参数来衡量有限元模型，该方法虽然避免了多个参数评估的弊端，但由于从实模态评估有限元模型，忽略了阻尼的影响，依然不能有效用于模型响应计算精度的预估。
总之，上述现有技术由于未考虑阻尼因素的影响，仅适用于对有限元模型的模态准确性的评估，而不能有效用于准确地预估有限元模型的动响应。
发明内容
本发明的目的在于克服现有技术的不足，提供一种含阻尼结构有限元模型的复模态评估方法。该方法一方面考虑了阻尼对结构响应的影响，且避免了多个参数评估模型的弊端，最终仅通过一个参数整体衡量有限元模型，另一方面在评估过程中直接利用试验测得的复模态进行相关性分析，避免了实数化过程中产生的二次误差，而且可预估模型响应计算的准确性，其计算方法实施规范，便于用计算机编程实施。
为实现上述目的，本发明采用了如下技术方案。
本发明的一种含阻尼结构有限元模型的复模态评估方法，其特征在于，包括以下步骤：
步骤一、根据阻尼结构试验复振型ΨT，求出阻尼结构试验复振型幅值和阻尼结构试验复振型相位θT；并且根据阻尼结构仿真复振型ΨA，求出阻尼结构仿真复振型幅值和阻尼结构仿真复振型相位θA；
步骤二、阻尼结构试验复模态和阻尼结构仿真复模态的相关性分析：由阻尼结构试验复频率λT和阻尼结构仿真复频率λA求出复频率相关性系数Rf，由阻尼结构试验复振型幅值和阻尼结构仿真复振型幅值求出复振型幅值相关性系数MACF，由阻尼结构试验复振型相位θT和阻尼结构仿真复振型相位θA求出复振型相位相关性系数Rθ；
步骤三、确定各阶模态权重集AI和复频率、复振型幅值、复振型相位三类影响因素权重集AII；
步骤四、将所述的复频率相关性系数Rf、复振型幅值相关性系数MACF及复振型相位相关性系数Rθ模糊化，得到复频率模糊关系矩阵R1、复振型幅值模糊关系矩阵R2、复振型相位模糊关系矩阵R3；
步骤五、第一阶段评估：对所述的复频率、复振型幅值、复振型相位三类影响因素分别进行单独评估，将三类影响因素的模糊关系矩阵分别乘以各阶模态权重集AI，得到复频率评估结果B1、复振型幅值评估结果B2、复振型相位评估结果B3，由此可得到第一阶段评估结果BI；
步骤六、第二阶段评估：对所述的复频率、复振型幅值、复振型相位三类影响因素进行综合评估，将第一阶段评估结果BI乘以三类影响因素权重集AII，并将评价等级量化，得到一个综合评估值P。
所述步骤一包括以下过程：
设阻尼结构试验复振型为ΨT，阻尼结构仿真复振型为ΨA，则阻尼结构试验复振型相位θT和阻尼结构仿真复振型相位θA分别为：
阻尼结构试验复振型幅值和阻尼结构仿真复振型幅值分别为：
所述步骤二包括以下过程：
2-1)复频率相关性分析：
设阻尼结构试验复频率为λT，阻尼结构仿真复频率为λA，则复频率相对误差δf定义为：
复频率相关性系数Rf则可定为：
Rf＝1-δf (6)
复频率相关性系数Rf为m×1维的列向量，m为评估模态阶数；
2-2)复振型幅值相关性分析
所述的复振型幅值相关性系数MACF可通过模态置信度准则求出：
其中，分别为阻尼结构仿真第i阶、试验第j阶的复振型幅值，MACF反映复振型幅值的相关性大小，复振型幅值相关性系数MACF为m×m维矩阵，m为评估模态阶数，在进行综合评估时，取矩阵对角线元素作为相关性系数；
2-3)复振型相位相关性分析，包括三部分：各自由度相位相对误差分析、各阶模态相位相对误差分析及总体相位相对误差分析；
各自由度相位相对误差Δθ可定义为：
通过各自由度相位相对误差Δθ，可判断该模型的相位误差在各个自由度上的分布情况；
各阶模态相位相对误差δθ及总体相位相对误差ΔF可通过求解相关矩阵范数来分析，首先定义一个矩阵范数：
则各阶模态相位相对误差δθ定义为：
则复振型相位相关性系数Rθ定义为：
Rθ＝1-δθ (11)
其中，θTi与θAi分别为阻尼结构试验和仿真第i阶复振型相位，复振型相位相关性系数Rθ为m×1维矩阵，m为评估模态阶数；
总体相位相对误差ΔF定义为：
所述步骤三包括以下过程：
3-1)确定各阶模态权重集AI：
在所述的第一阶段评估中，由于复频率相关性分析、复振型幅值相关性分析及复振型相位相关性分析都是对m阶模态进行的，因此三类影响因素都包含m阶模态的影响，各阶模态权重集AI可根据各阶模态对响应的影响来判定，对结构响应影响较大的模态其权重系数高；
所述的确定各阶模态权重集AI的一种方法为：由研究人员选取结构上比较重要或具有代表性的点，设结构上该点的频响函数固有频率处的幅值分别为a1、a2、…、an，其中m为模态阶数，设：
a＝a1+a2+...+am (13)
则各阶模态权重集AI为：
3-2)确定所述的复频率、复振型幅值、复振型相位三类影响因素权重集AII：
在所述的第二阶段评估中，认为复振型幅值及复振型相位对响应的影响略大于复频率的影响，则根据层次分析法，应用AHP理论，采用求和法求得三类影响因素权重集AII如下：
AII＝[0.20 0.40 0.40] (15)。
所述的确定各阶模态权重集AI的另一种方法为：
将不同载荷引起的动响应做模态分解，即
u(t)＝Ψq(t)＝q1(t)ψ1+q2(t)ψ2+...+q2n(t)ψ2n (16)
在模态坐标下，根据各阶模态的系数q1(t)、q2(t)、...、q2n(t)幅值所占比重来确定该特定载荷下各阶模态权重集AI，比重计算与第一种方法相同；
AI为第一阶段评估的权重集，三类影响因素中各阶模态权重集均为AI。
所述步骤四包括以下过程：
将评判标准设为5个等级：很好、好、一般、较差、很差；根据模糊评估理论，首先确定隶属度函数，隶属度函数是各个影响因素从属于各种不同评判标准的程度，在此取梯形函数，如下式：
根据梯形隶属度函数，将上述所求得的复频率相关性系数Rf、复振型幅值相关性系数MACF及复振型相位相关性系数Rθ模糊化，可分别确定复频率模糊关系矩阵R1、复振型幅值模糊关系矩阵R2、复振型相位模糊关系矩阵R3。
所述步骤五包括以下过程：
将所述的复频率、复振型幅值、复振型相位三类影响因素的模糊关系矩阵分别乘以各阶模态权重集AI，则第一阶段的评估结果为：
所述步骤六包括以下过程：
将第一阶段的评估结果BI乘以所述的复频率、复振型幅值、复振型相位三类影响因素权重集AII，得到：
BII＝AII×BI (19)
将5个评价等级量化：
C＝[0.1 0.2 0.3 0.4 0.5]T (20)
最终评估结果为：
P＝2×BII×C (21)。
与现有技术相比，本发明包括以下优点和有益效果：
1.本发明考虑了阻尼对结构响应的影响，且避免了多个参数评估模型的弊端，最终仅通过一个参数评估模型来整体衡量有限元模型。
2.本发明在评估过程中直接利用试验测得的复模态进行相关性分析，避免了实数化过程中产生的二次误差。
2.本发明既能进行有限元模型的模态准确性评估，又能预估模型响应计算的准确性，而且计算方法实施规范，便于用计算机编程实施。
附图说明
图1是本发明的一种实施例的含阻尼结构有限元模型的复模态评估方法的流程图。
具体实施方式
下面结合附图对本发明做进一步详细说明。
复模态下的频率和振型评估方法主要从复频率、复振型幅值及复振型相位三个方面进行试验与数值计算结果的相关性评估及综合评估，该方法将复频率、复振型幅值及复振型相位视做三个影响因素，首先分别求出试验和仿真结果在复频率、复振型幅值、复振型相位三方面的相关性系数，而每个影响因素的相关性系数包含m阶模态的结果，然后通过各阶模态取不同权重形成每一个影响因素的评估结果，最后通过三个影响因素取不同权重形成最终的一个评估参数。其过程可分为六个步骤。
步骤1、根据阻尼结构试验复振型ΨT求出阻尼结构试验复振型幅值和阻尼结构试验复振型相位θT，根据阻尼结构仿真复振型ΨA求出阻尼结构仿真复振型幅值和阻尼结构仿真复振型相位θA。
设阻尼结构试验复振型为ΨT，阻尼结构仿真复振型为ΨA，则阻尼结构试验复振型相位θT和阻尼结构仿真复振型相位θA分别为：
阻尼结构试验复振型幅值和阻尼结构仿真复振型幅值分别为：
该步骤一过程不需要将复振型转化为实振型，且复振型幅值及相位值具有明确的物理意义。
步骤2、阻尼结构试验复模态和阻尼结构仿真复模态的相关性分析：由阻尼结构试验复频率λT和阻尼结构仿真复频率λA求出复频率相关性系数Rf，由阻尼结构试验复振型幅值和阻尼结构仿真复振型幅值求出复振型幅值相关性系数MACF，由阻尼结构试验复振型相位θT和阻尼结构仿真复振型相位θA求出复振型相位相关性系数Rθ。
2.1、复频率相关性分析
设阻尼结构试验复频率为λT，阻尼结构仿真复频率为λA，则复频率相对误差δf定义为：
复频率相关性系数Rf则可定为：
Rf＝1-δf (6)
复频率相关性系数Rf为m×1维的列向量，m为评估模态阶数，该过程不需要将复频率转化为实频率，即可直接进行复频率的相关性分析。
2.2、复振型幅值相关性分析
复振型幅值相关性系数MACF可通过求模态置信度准则求出：
其中，分别为阻尼结构仿真第i阶、试验第j阶的复振型幅值，MACF反映复振型幅值的相关性大小，复振型幅值相关性系数MACF为m×m维矩阵，m为评估模态阶数，在进行综合评估时，取矩阵对角线元素作为相关性系数。
2.3、复振型相位相关性分析
复振型相位相关性分析可分为三部分，各自由度相位相对误差分析、各阶模态相位相对误差分析及总体相位相对误差分析。
各自由度相位相对误差Δθ可定义为：
通过各自由度相位相对误差Δθ，可判断该模型的相位误差在各个自由度上的分布情况。
各阶模态相位相对误差δθ及总体相位相对误差ΔF可通过求解相关矩阵范数来分析，首先定义一个矩阵范数：
则各阶模态相位相对误差δθ定义为：
则复振型相位相关性系数Rθ定义为：
Rθ＝1-δθ (11)
其中，θTi与θAi分别为阻尼结构试验和仿真第i阶复振型相位，复振型相位相关性系数Rθ为m×1维矩阵，m为评估模态阶数。
总体相位相对误差ΔF定义为：
总体相位相对误差ΔF可定量评价总体振型相位值的误差大小，总体相对误差值越小，模型越好，该值可为比较不同模型之间的优劣提供数据支持。
振型相位值相关性分析考虑了阻尼引起的各自由度上相位值的差异，而实模态评估中并未有反映该差异的具体参数。
步骤3、确定各阶模态权重集AI和复频率、复振型幅值、复振型相位三类影响因素权重集AII。
复模相关性分析主要从三个单方面来评价该模型的好坏，现利用模糊理论和层次分析法，以评估响应计算精度为目标，提出一种能综合衡量该模型的参数。
根据复频率、复振型幅值、复振型相位三个影响因素可分为两个阶段的评估，第一阶段为复频率、复振型幅值、复振型相位值三个方面的模糊评估，第二阶段为三个影响因素的综合评估。
3.1、确定各阶模态权重集AI
在第一阶段评估中，由于复频率相关性分析、复振型幅值相关性分析及复振型相位相关性分析都是对m阶模态进行的，因此三类影响因素都包含m阶模态的影响，各阶模态权重集AI可根据各阶模态对响应的影响来判定，对结构响应影响较大的模态其权重系数高。
各阶模态权重集AI选取有两种方法，一种方法可根据试验频响函数固有频率处幅值来确定，由研究人员选取结构上比较重要或具有代表性的点，设结构上该点的频响函数固有频率处的幅值分别为a1、a2、…、an，其中m为模态阶数，设：
a＝a1+a2+...+am (13)
则各阶模态权重集AI为：
各阶模态权重集AI选取的另一种方法是将不同载荷引起的动响应做模态分解，即
u(t)＝Ψq(t)＝q1(t)ψ1+q2(t)ψ2+...+q2n(t)ψ2n (15)
在模态坐标下，根据各阶模态的系数q1(t)、q2(t)、...、q2n(t)幅值所占比重来确定该特定载荷下各阶模态权重集AI，比重计算与第一种方法相同。
AI为第一阶段评估的权重集，三类影响因素中各阶模态权重集均为AI。
3.2、确定三类影响因素权重集AII
在第二阶段评估中，认为复振型幅值及复振型相位对响应的影响略大于复频率的影响，则根据层次分析法，应用AHP理论，采用求和法求得三类影响因素权重集AII如下：
AII＝[0.20 0.40 0.40] (16)
步骤4、将复频率相关性系数Rf、复振型幅值相关性系数MACF及复振型相位相关性系数Rθ模糊化，得到复频率模糊关系矩阵R1、复振型幅值模糊关系矩阵R2、复振型相位模糊关系矩阵R3。
本发明将评判标准设为5个等级(很好、好、一般、较差、很差)，根据模糊评估理论，首先确定隶属度函数，隶属度函数是各个影响因素从属于各种不同评判标准的程度，本发明取梯形函数，如下式：
根据梯形隶属度函数，将上述所求得的复频率相关性系数Rf、复振型幅值相关性系数MACF及复振型相位相关性系数Rθ模糊化，可分别确定复频率模糊关系矩阵R1、复振型幅值模糊关系矩阵R2、复振型相位模糊关系矩阵R3。
步骤5、上述第一阶段评估：对所述三类影响因素分别进行单独评估，将三类影响因素的模糊关系矩阵分别乘以各阶模态权重集AI，得到复频率评估结果B1、复振型幅值评估结果B2、复振型相位评估结果B3，由此可得到第一阶段评估结果BI。
将三类影响因素的模糊关系矩阵分别乘以各阶模态权重集AI，则第一阶段的评估结果为：
步骤6、上述第二阶段评估：对所述三类影响因素进行综合评估，将第一阶段评估结果BI乘以三类影响因素权重集AII，并将评价等级量化，得到一个综合评估值P。
将第一阶段的评估结果BI乘以三类影响因素权重集AII，得到：
BII＝AII×BI (19)
将5个评价等级量化：
C＝[0.1 0.2 0.3 0.4 0.5]T (20)
最终评估结果为：
P＝2×BII×C (21)
最终评估结果P为0～1范围内的某一值，P值越接近1，模型越好，在整个评估过程中，各阶模态权重集的选择与对响应的贡献有关，并且复振型幅值和复振型相位权重系数较高，故该综合评估数值主要是反映模型响应求解精度的总体情况。该方法最终评估参数仅含有一个P值，避免了多个参数难以总体把握模型的弊端，为多个模型的排队优选提供了数据支持。
图1是本发明的一种实施例的含阻尼结构有限元模型的复模态评估方法的流程图。
如图1所示，本发明的步骤1、2、4、5、6均对应流程图中每层的主框图，其步骤3对应最底两层框图的分支框。
总之，本发明的优点是：考虑了阻尼对结构响应的影响，且避免了多个参数评估模型的弊端，最终仅通过一个参数整体衡量有限元模型；在评估过程中直接利用试验测得的复模态进行相关性分析，避免了实数化过程中产生的二次误差；本发明既能进行有限元模型的模态准确性评估，又能预估模型响应计算的准确性，而且计算方法实施规范，便于用计算机编程实施。
Claims (8)
1.一种含阻尼结构有限元模型的复模态评估方法，其特征在于，包括以下步骤：
步骤一、根据阻尼结构试验复振型ΨT，求出阻尼结构试验复振型幅值和阻尼结构试验复振型相位θT；并且根据阻尼结构仿真复振型ΨA，求出阻尼结构仿真复振型幅值和阻尼结构仿真复振型相位θA；
步骤二、阻尼结构试验复模态和阻尼结构仿真复模态的相关性分析：由阻尼结构试验复频率λT和阻尼结构仿真复频率λA求出复频率相关性系数Rf，由阻尼结构试验复振型幅值和阻尼结构仿真复振型幅值求出复振型幅值相关性系数MACF，由阻尼结构试验复振型相位θT和阻尼结构仿真复振型相位θA求出复振型相位相关性系数Rθ；
步骤三、确定各阶模态权重集AI和复频率、复振型幅值、复振型相位三类影响因素权重集AII；
步骤四、将所述的复频率相关性系数Rf、复振型幅值相关性系数MACF及复振型相位相关性系数Rθ模糊化，得到复频率模糊关系矩阵R1、复振型幅值模糊关系矩阵R2、复振型相位模糊关系矩阵R3；
步骤五、第一阶段评估：对所述的复频率、复振型幅值、复振型相位三类影响因素分别进行单独评估，将三类影响因素的模糊关系矩阵分别乘以各阶模态权重集AI，得到复频率评估结果B1、复振型幅值评估结果B2、复振型相位评估结果B3，由此可得到第一阶段评估结果BI；
步骤六、第二阶段评估：对所述的复频率、复振型幅值、复振型相位三类影响因素进行综合评估，将第一阶段评估结果BI乘以三类影响因素权重集AII，并将评价等级量化，得到一个综合评估值P。
2.根据权利要求1所述的一种含阻尼结构有限元模型的复模态评估方法，其特征在于：所述步骤一包括以下过程：
设阻尼结构试验复振型为ΨT，阻尼结构仿真复振型为ΨA，则阻尼结构试验复振型相位θT和阻尼结构仿真复振型相位θA分别为：
<mrow>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>T</mi>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>arctan</mi>
<mfrac>
<mrow>
<mi>Im</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&Psi;</mi>
<mrow>
<mi>T</mi>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>Re</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&Psi;</mi>
<mrow>
<mi>T</mi>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>A</mi>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>arctan</mi>
<mfrac>
<mrow>
<mi>Im</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&Psi;</mi>
<mrow>
<mi>A</mi>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>Re</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&Psi;</mi>
<mrow>
<mi>A</mi>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
阻尼结构试验复振型幅值和阻尼结构仿真复振型幅值分别为：
3.根据权利要求1所述的一种含阻尼结构有限元模型的复模态评估方法，其特征在于：所述步骤二包括以下过程：
2-1)复频率相关性分析：
设阻尼结构试验复频率为λT，阻尼结构仿真复频率为λA，则复频率相对误差δf定义为：
<mrow>
<msub>
<mi>&delta;</mi>
<mi>f</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>|</mo>
<mrow>
<msub>
<mi>&lambda;</mi>
<mi>A</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&lambda;</mi>
<mi>T</mi>
</msub>
</mrow>
<mo>|</mo>
</mrow>
<mrow>
<mo>|</mo>
<msub>
<mi>&lambda;</mi>
<mi>T</mi>
</msub>
<mo>|</mo>
</mrow>
</mfrac>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
1
复频率相关性系数Rf则可定为：
Rf＝1-δf (6)
复频率相关性系数Rf为m×1维的列向量，m为评估模态阶数；
2-2)复振型幅值相关性分析
所述的复振型幅值相关性系数MACF可通过模态置信度准则求出：
其中，分别为阻尼结构仿真第i阶、试验第j阶的复振型幅值，MACF反映复振型幅值的相关性大小，复振型幅值相关性系数MACF为m×m维矩阵，m为评估模态阶数，在进行综合评估时，取矩阵对角线元素作为相关性系数；
2-3)复振型相位相关性分析，包括三部分：各自由度相位相对误差分析、各阶模态相位相对误差分析及总体相位相对误差分析；
各自由度相位相对误差Δθ可定义为：
<mrow>
<msub>
<mi>&Delta;</mi>
<mrow>
<mi>&theta;</mi>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>|</mo>
<mrow>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>A</mi>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>T</mi>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>|</mo>
</mrow>
<mrow>
<mo>|</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>T</mi>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>|</mo>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
通过各自由度相位相对误差Δθ，可判断该模型的相位误差在各个自由度上的分布情况；
各阶模态相位相对误差δθ及总体相位相对误差ΔF可通过求解相关矩阵范数来分析，首先定义一个矩阵范数：
<mrow>
<mo>|</mo>
<mo>|</mo>
<mi>A</mi>
<mo>|</mo>
<mo>|</mo>
<mo>=</mo>
<msqrt>
<mrow>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<munder>
<mo>&Sigma;</mo>
<mi>j</mi>
</munder>
<msubsup>
<mi>a</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
则各阶模态相位相对误差δθ定义为：
<mrow>
<msub>
<mi>&delta;</mi>
<mi>&theta;</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>T</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>A</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>|</mo>
<mo>|</mo>
</mrow>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>T</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
则复振型相位相关性系数Rθ定义为：
Rθ＝1-δθ (11)
其中，θTi与θAi分别为阻尼结构试验和仿真第i阶复振型相位，复振型相位相关性系数Rθ为m×1维矩阵，m为评估模态阶数；
总体相位相对误差ΔF定义为：
<mrow>
<msub>
<mi>&Delta;</mi>
<mi>F</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>&theta;</mi>
<mi>T</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mi>A</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
</mrow>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>&theta;</mi>
<mi>T</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
4.根据权利要求1所述的一种含阻尼结构有限元模型的复模态评估方法，其特征在于：所述步骤三包括以下过程：
3-1)确定各阶模态权重集AI：
在所述的第一阶段评估中，由于复频率相关性分析、复振型幅值相关性分析及复振型相位相关性分析都是对m阶模态进行的，因此三类影响因素都包含m阶模态的影响，各阶模态权重集AI可根据各阶模态对响应的影响来判定，对结构响应影响较大的模态其权重系数高；
所述的确定各阶模态权重集AI的方法为：由研究人员选取结构上比较重要或具有代表性的点，设结构上该点的频响函数固有频率处的幅值分别为a1、a2、…、an，其中m为模态阶数，设：
a＝a1+a2+...+am (13)
则各阶模态权重集AI为：
<mrow>
<msup>
<mi>A</mi>
<mi>I</mi>
</msup>
<mo>=</mo>
<mo>&lsqb;</mo>
<mtable>
<mtr>
<mtd>
<mfrac>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mi>a</mi>
</mfrac>
</mtd>
<mtd>
<mfrac>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mi>a</mi>
</mfrac>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mfrac>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<mi>a</mi>
</mfrac>
</mtd>
</mtr>
</mtable>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
3-2)确定所述的复频率、复振型幅值、复振型相位三类影响因素权重集AII：
在所述的第二阶段评估中，认为复振型幅值及复振型相位对响应的影响略大于复频率的影响，则根据层次分析法，应用AHP理论，采用求和法求得三类影响因素权重集AII如下：
AII＝[0.20 0.40 0.40] (15)。
5.根据权利要求4所述的一种含阻尼结构有限元模型的复模态评估方法，其特征在于：所述的确定各阶模态权重集AI的方法为：
将不同载荷引起的动响应做模态分解，即
u(t)＝Ψq(t)＝q1(t)ψ1+q2(t)ψ2+...+q2n(t)ψ2n (16)
在模态坐标下，根据各阶模态的系数q1(t)、q2(t)、...、q2n(t)幅值所占比重来确定该特定载荷下各阶模态权重集AI，比重计算与第一种方法相同；
AI为第一阶段评估的权重集，三类影响因素中各阶模态权重集均为AI。
6.根据权利要求1所述的一种含阻尼结构有限元模型的复模态评估方法，其特征在于：所述步骤四包括以下过程：
将评判标准设为5个等级：很好、好、一般、较差、很差；根据模糊评估理论，首先确定隶属度函数，隶属度函数是各个影响因素从属于各种不同评判标准的程度，在此取梯形函数，如下式：
<mrow>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>x</mi>
<mo>&le;</mo>
<mi>a</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<mi>a</mi>
<mo>)</mo>
<mo>/</mo>
<mo>(</mo>
<mi>b</mi>
<mo>-</mo>
<mi>a</mi>
<mo>)</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>a</mi>
<mo><</mo>
<mi>x</mi>
<mo><</mo>
<mi>b</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mrow>
<mi>b</mi>
<mo>&le;</mo>
<mi>x</mi>
<mo>&le;</mo>
<mi>c</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<mi>c</mi>
<mo>-</mo>
<mi>x</mi>
<mo>)</mo>
<mo>/</mo>
<mo>(</mo>
<mi>d</mi>
<mo>-</mo>
<mi>c</mi>
<mo>)</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>c</mi>
<mo><</mo>
<mi>x</mi>
<mo>&le;</mo>
<mi>d</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>d</mi>
<mo><</mo>
<mi>x</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>17</mn>
<mo>)</mo>
</mrow>
</mrow>
根据梯形隶属度函数，将上述所求得的复频率相关性系数Rf、复振型幅值相关性系数MACF及复振型相位相关性系数Rθ模糊化，可分别确定复频率模糊关系矩阵R1、复振型幅值模糊关系矩阵R2、复振型相位模糊关系矩阵R3。
7.根据权利要求1所述的一种含阻尼结构有限元模型的复模态评估方法，其特征在于：所述步骤五包括以下过程：
将所述的复频率、复振型幅值、复振型相位三类影响因素的模糊关系矩阵分别乘以各阶模态权重集AI，则第一阶段的评估结果为：
<mrow>
<msup>
<mi>B</mi>
<mi>I</mi>
</msup>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>B</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>B</mi>
<mn>3</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>A</mi>
<mi>I</mi>
</msup>
<mo>&times;</mo>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>A</mi>
<mi>I</mi>
</msup>
<mo>&times;</mo>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>A</mi>
<mi>I</mi>
</msup>
<mo>&times;</mo>
<msub>
<mi>R</mi>
<mn>3</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>18</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
8.根据权利要求1所述的一种含阻尼结构有限元模型的复模态评估方法，其特征在于：所述步骤六包括以下过程：
将第一阶段的评估结果BI乘以所述的复频率、复振型幅值、复振型相位三类影响因素权重集AII，得到：
BII＝AII×BI (19)
将5个评价等级量化：
C＝[0.1 0.2 0.3 0.4 0.5]T (20)
最终评估结果为：
P＝2×BII×C (21)。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710521580.2A CN107292045B (zh) | 2017-06-30 | 2017-06-30 | 一种含阻尼结构有限元模型的复模态评估方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710521580.2A CN107292045B (zh) | 2017-06-30 | 2017-06-30 | 一种含阻尼结构有限元模型的复模态评估方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107292045A true CN107292045A (zh) | 2017-10-24 |
| CN107292045B CN107292045B (zh) | 2020-07-28 |
Family
ID=60099390
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710521580.2A Active CN107292045B (zh) | 2017-06-30 | 2017-06-30 | 一种含阻尼结构有限元模型的复模态评估方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107292045B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112541293A (zh) * | 2020-12-22 | 2021-03-23 | 中国人民解放军陆军装甲兵学院 | 一种含锰铜基阻尼合金的悬臂梁结构仿真分析方法 |
| CN115563844A (zh) * | 2022-09-27 | 2023-01-03 | 中国海洋大学 | 一种基于频响函数形状相关系数的模型修正方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7206709B2 (en) * | 2003-05-29 | 2007-04-17 | Carnegie Mellon University | Determination of damping in bladed disk systems using the fundamental mistuning model |
| CN102982202A (zh) * | 2012-11-15 | 2013-03-20 | 中国电子科技集团公司第二十八研究所 | 基于缺损模态的结构模型修正方法 |
| CN104008234A (zh) * | 2014-05-12 | 2014-08-27 | 南京航空航天大学 | 密集模态含阻尼结构模型修正方法 |
| JP2016103101A (ja) * | 2014-11-27 | 2016-06-02 | 株式会社東芝 | 耐震解析装置、方法及びプログラム |
| CN105701333A (zh) * | 2015-12-31 | 2016-06-22 | 河海大学 | 一种空间预应力索杆体系整体安全性评价方法 |
-
2017
- 2017-06-30 CN CN201710521580.2A patent/CN107292045B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7206709B2 (en) * | 2003-05-29 | 2007-04-17 | Carnegie Mellon University | Determination of damping in bladed disk systems using the fundamental mistuning model |
| CN102982202A (zh) * | 2012-11-15 | 2013-03-20 | 中国电子科技集团公司第二十八研究所 | 基于缺损模态的结构模型修正方法 |
| CN104008234A (zh) * | 2014-05-12 | 2014-08-27 | 南京航空航天大学 | 密集模态含阻尼结构模型修正方法 |
| JP2016103101A (ja) * | 2014-11-27 | 2016-06-02 | 株式会社東芝 | 耐震解析装置、方法及びプログラム |
| CN105701333A (zh) * | 2015-12-31 | 2016-06-22 | 河海大学 | 一种空间预应力索杆体系整体安全性评价方法 |
Non-Patent Citations (2)
| Title |
|---|
| 张安平等: "动力学有限元模型的模糊评估", 《南京航空航天大学学报》 * |
| 沈敏: "基于飞行测试模态的物理模型的建立及动载荷位置识别", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112541293A (zh) * | 2020-12-22 | 2021-03-23 | 中国人民解放军陆军装甲兵学院 | 一种含锰铜基阻尼合金的悬臂梁结构仿真分析方法 |
| CN112541293B (zh) * | 2020-12-22 | 2024-05-31 | 中国人民解放军陆军装甲兵学院 | 一种含锰铜基阻尼合金的悬臂梁结构仿真分析方法 |
| CN115563844A (zh) * | 2022-09-27 | 2023-01-03 | 中国海洋大学 | 一种基于频响函数形状相关系数的模型修正方法 |
| CN115563844B (zh) * | 2022-09-27 | 2023-08-22 | 中国海洋大学 | 一种基于频响函数形状相关系数的模型修正方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107292045B (zh) | 2020-07-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103065042B (zh) | 基于情景的多目标综合决策评价方法 | |
| CN105608263B (zh) | 一种面向涡轮叶盘结构寿命概率分析的自适应处理方法 | |
| US20100058257A1 (en) | Topology optimization method using equivalent static loads | |
| CN111143981B (zh) | 虚拟试验模型验证系统及方法 | |
| CN107085633B (zh) | 基于支持向量机的多点振动响应频域预测的装置及方法 | |
| CN105572572B (zh) | 基于wknn-lssvm的模拟电路故障诊断方法 | |
| CN106372278A (zh) | 一种联合考虑输入参数不确定性和代理模型不确定性的灵敏度分析方法 | |
| CN110457750B (zh) | 基于神经网络响应面的有载分接开关弹簧储能不足故障识别方法 | |
| JP2003521015A (ja) | 非線形モデルの発生及びそれを用いたシミュレーション試験用駆動信号の発生 | |
| CN113111547A (zh) | 基于缩减基的频域有限元模型修正方法 | |
| Chen | FE model validation for structural dynamics | |
| CN107341332A (zh) | 一种车用电机驱动系统的评价指标权重的确定方法 | |
| CN114583767B (zh) | 一种数据驱动的风电场调频响应特性建模方法及系统 | |
| CN113360983A (zh) | 一种边坡可靠度分析与风险评估方法 | |
| CN116127695A (zh) | 一种基于综合性能评价的生产线构建方法及系统 | |
| CN107292045B (zh) | 一种含阻尼结构有限元模型的复模态评估方法 | |
| CN110110406B (zh) | 一种基于Excel计算平台实现LS-SVM模型的边坡稳定性预测方法 | |
| CN109359388A (zh) | 一种复杂仿真系统可信度评估方法 | |
| CN111445025A (zh) | 确定业务模型超参数的方法和装置 | |
| CN110825583A (zh) | 一种针对云数据中心多指标融合的能效定性评估技术 | |
| CN114218819B (zh) | 基于自适应元模型的复杂结构有限元参数区间的修正方法 | |
| CN114970072A (zh) | 一种应用于配电网仿真的分区方案的评价方法及系统 | |
| CN112182739A (zh) | 一种飞行器结构非概率可信可靠性拓扑优化设计方法 | |
| CN111507003A (zh) | 基于质量相关动态特性提取的脱丁烷塔关键变量预测方法 | |
| Cudney et al. | Predicting vehicle cost using the T-method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |