CN104008234A - 密集模态含阻尼结构模型修正方法 - Google Patents
密集模态含阻尼结构模型修正方法 Download PDFInfo
- Publication number
- CN104008234A CN104008234A CN201410199353.9A CN201410199353A CN104008234A CN 104008234 A CN104008234 A CN 104008234A CN 201410199353 A CN201410199353 A CN 201410199353A CN 104008234 A CN104008234 A CN 104008234A
- Authority
- CN
- China
- Prior art keywords
- lambda
- phi
- psi
- matrix
- sigma
- 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
Landscapes
- Complex Calculations (AREA)
Abstract
本发明公开了一种密集模态含阻尼结构模型修正方法,其主要特点是:设定有限元各矩阵的误差是各单元相应子矩阵的误差累积而成,并假设初始的各单元子矩阵与准确矩阵只相差修正因子;在修正算法中首先将动力学方程转换为状态方程,并求特征解;对不同阶的初始状态方程与假设的准确的状态方程的特征方程分别前乘计算与试验得到的不同阶复振型,再相减并整理,得到关于修正因子的线性方程组;用最小二乘法计算修正因子,进而得到修正后的有限元矩阵。该方法实施规范,便于计算机编程实施,由于不要求获知试验与计算模态的一一对应情况,因此适用于密集模态含阻尼结构的模型修正。修正后的模型不仅可用于模态计算,还可用于动力学响应计算分析。
Description
技术领域
本发明涉及一种密集模态含阻尼结构模型修正方法,特别适用于飞机、汽车等具有密集模态特点的含阻尼结构的模型修正。
背景技术
在实际工程中,许多机械结构都是具有密集模态的含阻尼结构,利用试验结果修正初始建立的有限元模型是这类结构动力学准确建模的必由之路。但截至目前模型修正都是利用正则模态测试结果修正无阻尼模型;而且一般要求模态分离清楚,即不能出现密集模态现象,以保证测试模态与有限元计算模态的阶次对应关系明确。2007年Sau-LonJamesHu基于经典无阻尼系统提出了CMCM法,首次避开了模型修正中模态配对难题,但不适用于工程中最广泛存在的阻尼系统,而且其基本方程假设也决定了无法考虑阻尼的影响。中国专利CN101794338A公开了一种基于结构模态试验的矩阵型动力学模型修正方法,该方法利用多次同时改变结构质量分布和刚度后产生的新结构的模态试验结果,运用矩阵计算和代数方程求解,计算出有限元模型的质量阵和刚度阵的修正量。虽然该方法提高了修正精度,但它既不能处理含阻尼结构的模型修正问题,也不能解决含密集模态结构的模型修正问题。上述方法修正得到的模型除精度问题外,由于不涉及阻尼矩阵的修正,因此只能用于进行结构模态计算,而不能有效用于动力学响应计算分析。
发明内容
本发明的目的在于克服现有技术的不足,提供一种密集模态含阻尼结构模型修正方法。该方法一方面能够修正含阻尼结构的动力学模型,通过该方法获得的模型能够有效地应用于动力学响应计算分析。另一方面在修正过程中不要求获知试验与计算模态的一一对应情况,因此特别适用于密集模态含阻尼结构的模型修正。而且本发明的计算方法实施规范,便于用计算机编程实施。
为实现上述目的,本发明采用了如下技术方案:
本发明密集模态含阻尼结构模型修正方法,其特征在于:对于初始建立的有限元模型,假设其刚度、质量、阻尼矩阵分别为K、M、C,各矩阵规模为m×m,而各矩阵分别为Ne个元素的单元矩阵叠加构成;
通过对含阻尼结构进行复模态测试,来得到结构的若干个特征解,具体可假设包括2N个复特征值(r=1,2,…2N)及2N个复特征向量(r=1,2,…2N);
利用K、M、C构建A,B矩阵,建立状态方程,实现微分方程的降阶,基于复模态理论与矩阵计算,构建出以各元素刚度、质量、阻尼的比例修正系数为未知量的线性方程组,再通过最小二乘法来获得这些修正系数的近似解;
粘性阻尼系统的自由振动方程为
其特征值方程
(λ2Μ+λC+Κ)Φ=0 (2)
引入由位移和速度所组成的2N维状态向量
改写为状态方程形式,此时,方程(1)可写作由状态向量描述的一阶线性微分方程组
其中
设系统在状态空间中的运动为
v(t)=Ψeλt (6)
相应的特征值问题为
或
BΨ=λAΨ (8)
将(8)展开与(2)比较,不难看出它们具有相同的特征值λr(r=1,2,…2N)并且特征向量满足以下关系
对于含有粘性阻尼结构的有限元模型(或“分析模型”),设Ψi,λi分别为分析模型的第i阶复振型和第i阶复频率,根据(3.34)式有
BΨi=λiAΨi (10)
同样地假设对应结构的真实的质量,刚度和阻尼矩阵分别为K*,M*,C*,并设分别为“真实模型”的第j阶复振型和第j阶复频率;它们自然也满足
其中
假定真实模型与分析模型具有以下修正关系
其中,Κn,Μn,Cn是全局坐标下第n个单元的刚度矩阵,质量矩阵和阻尼矩阵,αn,βn,γn是相应的修正系数,Ne是分析模型中有限单元的个数;
方程(10)两边左乘方程(11)两边左乘(Ψi)T得到
和
已知A,B是对称矩阵,将(17)式两边转置得
(18)式减(19)式,并用(14)(15)(16)整理得
展开上式并令:
得到
写成矩阵形式就得到
Uα+Vβ+Wγ=f (23)
其中U,V和W是Nm×Ne矩阵;α,β和γ是Ne维列向量,f是Nm维列向量;进一步地,(23)式可以写成
Gτ=f (24)
式中
G=[U V W] (25)
且
(24)式为超定方程组,用最小二乘求解就得到修正参数的估计值,然后通过式(14)-(16)可以获得修正后准确的的刚度、质量、阻尼矩阵Κ*,Μ*,C*。
附图说明
图1是本发明密集模态含阻尼结构模型修正方法的流程图。
具体实施方式
下面结合附图对本发明的技术方案进行详细说明。
本发明密集模态含阻尼结构模型修正方法的特征在于:
对于初始建立的有限元模型,假设其刚度、质量、阻尼矩阵分别为K、M、C,各矩阵规模为m×m,而各矩阵分别为Ne个元素的单元矩阵叠加构成;
通过对含阻尼结构进行复模态测试,来得到结构的若干个特征解,具体可假设包括2N个复特征值(r=1,2,…2N)及2N个复特征向量(r=1,2,…2N),该实测得到复模态不必连续按阶次排列,并允许一些阶次的复模态未被测量出来,更不要求知道测量复模态与计算模态的一一对应情况;
利用K、M、C构建A,B矩阵,建立状态方程,实现微分方程的降阶,基于复模态理论与矩阵计算,构建出以各元素刚度、质量、阻尼的比例修正系数为未知量的线性方程组,再通过最小二乘法来获得这些修正系数的近似解;
粘性阻尼系统的自由振动方程为
其特征值方程
(λ2Μ+λC+Κ)Φ=0 (2)
引入由位移和速度所组成的2N维状态向量
改写为状态方程形式,此时,方程(1)可写作由状态向量描述的一阶线性微分方程组
其中
设系统在状态空间中的运动为
v(t)=Ψeλt (6)
相应的特征值问题为
或
BΨ=λAΨ (8)
将(8)展开与(2)比较,不难看出它们具有相同的特征值λr(r=1,2,…2N)并且特征向量满足以下关系
对于含有粘性阻尼结构的有限元模型(或“分析模型”),设Ψi,λi分别为分析模型的第i阶复振型和第i阶复频率,根据(3.34)式有
BΨi=λiAΨi (10)
同样地假设对应结构的真实的质量,刚度和阻尼矩阵分别为K*,M*,C*,并设分别为“真实模型”的第j阶复振型和第j阶复频率;它们自然也满足
其中
假定真实模型与分析模型具有以下修正关系
其中,Κn,Μn,Cn是全局坐标下第n个单元的刚度矩阵,质量矩阵和阻尼矩阵,αn,βn,γn是相应的修正系数,Ne是分析模型中有限单元的个数;
方程(10)两边左乘方程(11)两边左乘(Ψi)T得到
和
已知A,B是对称矩阵,将(17)式两边转置得
(18)式减(19)式,并用(14)(15)(16)整理得
展开上式并令:
得到
写成矩阵形式就得到
Uα+Vβ+Wγ=f (23)
其中U,V和W是Nm×Ne矩阵;α,β和γ是Ne维列向量,f是Nm维列向量;进一步地,(23)式可以写成
Gτ=f (24)
式中
G=[U V W] (25)
且
(24)式为超定方程组,用最小二乘求解就得到修正参数的估计值,然后通过式(14)-(16)可以获得修正后准确的的刚度、质量、阻尼矩阵Κ*,Μ*,C*。
本发明的优点是:实施规范,便于计算机编程实施,由于不要求获知试验与计算模态的一一对应情况,因此适用于密集模态含阻尼结构的模型修正。修正后的模型不仅可用于模态计算,还可用于动力学响应计算分析。
Claims (1)
1.一种密集模态含阻尼结构模型修正方法,其特征在于:
对于初始建立的有限元模型,假设其刚度、质量、阻尼矩阵分别为K、M、C,各矩阵规模为m×m,而各矩阵分别为Ne个元素的单元矩阵叠加构成;
通过对含阻尼结构进行复模态测试,来得到结构的若干个特征解,具体可假设包括2N个复特征值(r=1,2,…2N)及2N个复特征向量(r=1,2,…2N);
利用K、M、C构建A,B矩阵,建立状态方程,实现微分方程的降阶,基于复模态理论与矩阵计算,构建出以各元素刚度、质量、阻尼的比例修正系数为未知量的线性方程组,再通过最小二乘法来获得这些修正系数的近似解;
粘性阻尼系统的自由振动方程为
其特征值方程
(λ2Μ+λC+Κ)Φ=0 (2)
引入由位移和速度所组成的2N维状态向量
改写为状态方程形式,此时,方程(1)可写作由状态向量描述的一阶线性微分方程组
其中
设系统在状态空间中的运动为
v(t)=Ψeλt (6)
相应的特征值问题为
或
BΨ=λAΨ (8)
将(8)展开与(2)比较,不难看出它们具有相同的特征值λr(r=1,2,…2N)并且特征向量满足以下关系
对于含有粘性阻尼结构的有限元模型(或“分析模型”),设Ψi,λi分别为分析模型的第i阶复振型和第i阶复频率,根据(3.34)式有
BΨi=λiAΨi (10)
同样地假设对应结构的真实的质量,刚度和阻尼矩阵分别为K*,M*,C*,并设分别为“真实模型”的第j阶复振型和第j阶复频率;它们自然也满足
其中
假定真实模型与分析模型具有以下修正关系
其中,Κn,Μn,Cn是全局坐标下第n个单元的刚度矩阵,质量矩阵和阻尼矩阵,αn,βn,γn是相应的修正系数,Ne是分析模型中有限单元的个数;
方程(10)两边左乘方程(11)两边左乘(Ψi)T得到
和
已知A,B是对称矩阵,将(17)式两边转置得
(18)式减(19)式,并用(14)(15)(16)整理得
展开上式并令:
得到
写成矩阵形式就得到
Uα+Vβ+Wγ=f (23)
其中U,V和W是Nm×Ne矩阵;α,β和γ是Ne维列向量,f是Nm维列向量;进一步地,(23)式可以写成
Gτ=f (24)
式中
G=[U V W] (25)
且
(24)式为超定方程组,用最小二乘求解就得到修正参数的估计值,然后通过式(14)-(16)可以获得修正后准确的的刚度、质量、阻尼矩阵Κ*,Μ*,C*。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410199353.9A CN104008234B (zh) | 2014-05-12 | 2014-05-12 | 密集模态含阻尼结构模型修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410199353.9A CN104008234B (zh) | 2014-05-12 | 2014-05-12 | 密集模态含阻尼结构模型修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104008234A true CN104008234A (zh) | 2014-08-27 |
CN104008234B CN104008234B (zh) | 2017-02-22 |
Family
ID=51368890
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410199353.9A Expired - Fee Related CN104008234B (zh) | 2014-05-12 | 2014-05-12 | 密集模态含阻尼结构模型修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104008234B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105956213A (zh) * | 2016-03-07 | 2016-09-21 | 重庆大学 | 一种机敏约束层阻尼板结构模态状态观测器的设计方法 |
CN106682301A (zh) * | 2016-12-23 | 2017-05-17 | 西安电子科技大学 | 一种面向多点激励振动试验的系统模型修正方法 |
CN107292045A (zh) * | 2017-06-30 | 2017-10-24 | 南京航空航天大学 | 一种含阻尼结构有限元模型的复模态评估方法 |
-
2014
- 2014-05-12 CN CN201410199353.9A patent/CN104008234B/zh not_active Expired - Fee Related
Non-Patent Citations (1)
Title |
---|
王轲: "基于模态试验与优化的静动力学模型转换", 《第十届全国振动理论及应用学术会议论文集》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105956213A (zh) * | 2016-03-07 | 2016-09-21 | 重庆大学 | 一种机敏约束层阻尼板结构模态状态观测器的设计方法 |
CN105956213B (zh) * | 2016-03-07 | 2018-12-21 | 重庆大学 | 一种机敏约束层阻尼板结构模态状态观测器的设计方法 |
CN106682301A (zh) * | 2016-12-23 | 2017-05-17 | 西安电子科技大学 | 一种面向多点激励振动试验的系统模型修正方法 |
CN107292045A (zh) * | 2017-06-30 | 2017-10-24 | 南京航空航天大学 | 一种含阻尼结构有限元模型的复模态评估方法 |
CN107292045B (zh) * | 2017-06-30 | 2020-07-28 | 南京航空航天大学 | 一种含阻尼结构有限元模型的复模态评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104008234B (zh) | 2017-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Taherizadeh et al. | A non-associated constitutive model with mixed iso-kinematic hardening for finite element simulation of sheet metal forming | |
CN103299308A (zh) | 利用实体建模计算任意形状密闭容器中的液位 | |
Chwalowski et al. | Preliminary computational analysis of the (hirenasd) configuration in preparation for the aeroelastic prediction workshop | |
CN102759573B (zh) | 基于频率变化的结构损伤定位与损伤程度的评估方法 | |
CN102778240B (zh) | 一种基于物理模型的陀螺贮存寿命加速试验方案确定方法 | |
CN103324798B (zh) | 基于区间响应面模型的随机模型修正方法 | |
CN105468851A (zh) | 一种确定飞机动态重量特性的方法 | |
CN104296924B (zh) | 动力总成惯性参数测试方法 | |
CN104992002A (zh) | 一种面向智能蒙皮天线的应变传感器布局方法 | |
CN103198215A (zh) | 一种基于差异进化支持向量机的坑外土体沉降预测方法 | |
CN104008234A (zh) | 密集模态含阻尼结构模型修正方法 | |
CN104281730A (zh) | 一种大转动变形的板壳结构动响应的有限元分析方法 | |
CN104424380A (zh) | 机械变形测量仪的模拟 | |
CN101794338A (zh) | 基于结构模态试验的矩阵型动力学模型修正方法 | |
CN104462022A (zh) | 飞行器动力学系统参数可辨识性分析方法 | |
CN104834795A (zh) | 包带连接结构接触摩擦非线性特性模拟方法及系统 | |
CN105136150A (zh) | 一种基于多次星敏感器测量信息融合的姿态确定方法 | |
CN104408326B (zh) | 一种对深空探测自主导航滤波算法的评估方法 | |
NO20111686A1 (no) | Konvertering av reservoirinndeling for skalering av belastning | |
Thomas et al. | Static/dynamic correction approach for reduced-order modeling of unsteady aerodynamics | |
CN106342305B (zh) | 一种面向多任务要求的测试性指标确定方法 | |
CN102930161A (zh) | 基于多舱段航天器混合模态预示方法 | |
CN105093933B (zh) | 一种确定lpv变增益控制器的方法 | |
Roos | Generation of flexible aircraft LFT models for robustness analysis | |
Giorgi et al. | NLFK4ALL: An open-source demonstration toolbox for computationally efficient nonlinear Froude-Krylov force calculations |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170222 |
|
CF01 | Termination of patent right due to non-payment of annual fee |