CN104008234B - 密集模态含阻尼结构模型修正方法 - Google Patents
密集模态含阻尼结构模型修正方法 Download PDFInfo
- Publication number
- CN104008234B CN104008234B CN201410199353.9A CN201410199353A CN104008234B CN 104008234 B CN104008234 B CN 104008234B CN 201410199353 A CN201410199353 A CN 201410199353A CN 104008234 B CN104008234 B CN 104008234B
- Authority
- CN
- China
- Prior art keywords
- lambda
- phi
- equation
- psi
- matrix
- 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.)
- Expired - Fee Related
Links
Landscapes
- Complex Calculations (AREA)
Abstract
本发明公开了一种密集模态含阻尼结构模型修正方法,其主要特点是:设定有限元各矩阵的误差是各单元相应子矩阵的误差累积而成,并假设初始的各单元子矩阵与准确矩阵只相差修正因子;在修正算法中首先将动力学方程转换为状态方程,并求特征解;对不同阶的初始状态方程与假设的准确的状态方程的特征方程分别前乘计算与试验得到的不同阶复振型,再相减并整理,得到关于修正因子的线性方程组;用最小二乘法计算修正因子,进而得到修正后的有限元矩阵。该方法实施规范,便于计算机编程实施,由于不要求获知试验与计算模态的一一对应情况,因此适用于密集模态含阻尼结构的模型修正。修正后的模型不仅可用于模态计算,还可用于动力学响应计算分析。
Description
技术领域
本发明涉及一种密集模态含阻尼结构模型修正方法,特别适用于飞机、汽车等具有密集模态特点的含阻尼结构的模型修正。
背景技术
在实际工程中,许多机械结构都是具有密集模态的含阻尼结构,利用试验结果修正初始建立的有限元模型是这类结构动力学准确建模的必由之路。但截至目前模型修正都是利用正则模态测试结果修正无阻尼模型;而且一般要求模态分离清楚,即不能出现密集模态现象,以保证测试模态与有限元计算模态的阶次对应关系明确。2007年Sau-Lon JamesHu基于经典无阻尼***提出了CMCM法,首次避开了模型修正中模态配对难题,但不适用于工程中最广泛存在的阻尼***,而且其基本方程假设也决定了无法考虑阻尼的影响。中国专利CN101794338A公开了一种基于结构模态试验的矩阵型动力学模型修正方法,该方法利用多次同时改变结构质量分布和刚度后产生的新结构的模态试验结果,运用矩阵计算和代数方程求解,计算出有限元模型的质量阵和刚度阵的修正量。虽然该方法提高了修正精度,但它既不能处理含阻尼结构的模型修正问题,也不能解决含密集模态结构的模型修正问题。上述方法修正得到的模型除精度问题外,由于不涉及阻尼矩阵的修正,因此只能用于进行结构模态计算,而不能有效用于动力学响应计算分析。
发明内容
本发明的目的在于克服现有技术的不足,提供一种密集模态含阻尼结构模型修正方法。该方法一方面能够修正含阻尼结构的动力学模型,通过该方法获得的模型能够有效地应用于动力学响应计算分析。另一方面在修正过程中不要求获知试验与计算模态的一一对应情况,因此特别适用于密集模态含阻尼结构的模型修正。而且本发明的计算方法实施规范,便于用计算机编程实施。
为实现上述目的,本发明采用了如下技术方案:
本发明密集模态含阻尼结构模型修正方法,其特征在于:对于初始建立的有限元模型,假设其刚度、质量、阻尼矩阵分别为K、M、C,各矩阵规模为m×m,而各矩阵分别为Ne个元素的单元矩阵叠加构成;
通过对含阻尼结构进行复模态测试,来得到结构的多个特征解,具体可假设包括2N个复特征值及2N个复特征向量其中,r=1,2,…2N;
利用K、M、C构建A,B矩阵,建立状态方程,实现微分方程的降阶,基于复模态理论与矩阵计算,构建出以各元素刚度、质量、阻尼的比例修正系数为未知量的线性方程组,再通过最小二乘法来获得这些修正系数的近似解;
粘性阻尼***的自由振动方程为
其中,u(t)表示位移;在振动理论中,该方程为运动微分方程。
其特征值方程
(λ2Μ+λC+Κ)Φ=0 (2)
在振动理论中,λ表示特征方程的特征根,Φ表示特征方程的特征向量。
引入由位移和速度所组成的2N维状态向量
改写为状态方程形式,此时,方程(1)可写作由状态向量描述的一阶线性微分方程组
在振动理论中,v(t)表示状态向量,由速度和位移组成,v(0)为零时刻的状态向量的值,值的大小为。
在振动理论中,u0表示零时刻的位移。
设***在状态空间中的运动为
v(t)=Ψeλt (6)
在振动理论中,Ψ表示状态空间方程的特征向量。
相应的特征值问题为
或
BΨ=λAΨ (8)
将(8)展开与(2)比较,不难看出它们具有相同的特征值λr(r=1,2,…2N)并且特征向量满足以下关系
在振动理论中,Ψr为状态空间描述的第r阶特征向量,Φr为物理空间描述的第r阶特征向量。
对于含有粘性阻尼结构的分析模型”,设Ψi,λi分别为分析模型的第i阶复振型和第i阶复频率,根据(8)式有
BΨi=λiAΨi (10)
同样地假设对应结构的真实的质量,刚度和阻尼矩阵分别为K*,M*,C*,并设分别为“真实模型”的第j阶复振型和第j阶复频率;它们自然也满足
在振动理论中,对于真实结构,A*、B*为的状态空间方程参数,的状态空间方程第j阶特征向量。
假定真实模型与分析模型具有以下修正关系
其中,Κn,Μn,Cn是全局坐标下第n个单元的刚度矩阵,质量矩阵和阻尼矩阵,αn,βn,γn是相应的修正系数,Ne是分析模型中有限单元的个数;
方程(10)两边左乘方程(11)两边左乘(Ψi)T得到
和
已知A,B是对称矩阵,将(17)式两边转置得
(18)式减(19)式,并用(14)(15)(16)整理得
在振动理论中,为第i阶特征向量的转置。
展开上式并令:
得到
写成矩阵形式就得到
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)可以获得修正后准确的的刚度、质量、阻尼矩阵K′,M′,C′。
附图说明
图1是本发明密集模态含阻尼结构模型修正方法的流程图。
具体实施方式
下面结合附图对本发明的技术方案进行详细说明。
本发明密集模态含阻尼结构模型修正方法的特征在于:
对于初始建立的有限元模型,假设其刚度、质量、阻尼矩阵分别为K、M、C,各矩阵规模为m×m,而各矩阵分别为Ne个元素的单元矩阵叠加构成;
通过对含阻尼结构进行复模态测试,来得到结构的多个特征解,具体可假设包括2N个复特征值及2N个复特征向量其中,r=1,2,…2N。该实测得到复模态不必连续按阶次排列,并允许一些阶次的复模态未被测量出来,更不要求知道测量复模态与计算模态的一一对应情况;
利用K、M、C构建A,B矩阵,建立状态方程,实现微分方程的降阶,基于复模态理论与矩阵计算,构建出以各元素刚度、质量、阻尼的比例修正系数为未知量的线性方程组,再通过最小二乘法来获得这些修正系数的近似解;
粘性阻尼***的自由振动方程为
在振动理论中,u(t)表示位移;在振动理论中,该方程为运动微分方程。具体可以参考北京航空航天大学出版社,胡海岩主编的《机械振动基础》。
其特征值方程
(λ2Μ+λC+Κ)Φ=0 (2)
在振动理论中,λ表示特征方程的特征根,Φ表示特征方程的特征向量。
引入由位移和速度所组成的2N维状态向量
改写为状态方程形式,此时,方程(1)可写作由状态向量描述的一阶线性微分方程组
在振动理论中,v(t)表示状态向量,由速度和位移组成,v(0)为零时刻的状态向量的值,值的大小为v0。
在振动理论中,u0表示零时刻的位移。
设***在状态空间中的运动为
v(t)=Ψeλt (6)
在振动理论中,Ψ表示状态空间方程的特征向量。
相应的特征值问题为
或
BΨ=λAΨ (8)
将(8)展开与(2)比较,不难看出它们具有相同的特征值λr(r=1,2,…2N)并且特征向量满足以下关系
在振动理论中,Ψr为状态空间描述的第r阶特征向量,Φr为物理空间描述的第r阶特征向量。
对于含有粘性阻尼结构的分析模型”,设Ψi,λi分别为分析模型的第i阶复振型和第i阶复频率,根据(8)式有
BΨi=λiAΨi (10)
同样地假设对应结构的真实的质量,刚度和阻尼矩阵分别为K*,M*,C*,并设分别为“真实模型”的第j阶复振型和第j阶复频率;它们自然也满足
在振动理论中,对于真实结构,A*、B*为的状态空间方程参数,的状态空间方程第j阶特征向量。
假定真实模型与分析模型具有以下修正关系
其中,Κn,Μn,Cn是全局坐标下第n个单元的刚度矩阵,质量矩阵和阻尼矩阵,αn,βn,γn是相应的修正系数,Ne是分析模型中有限单元的个数;
方程(10)两边左乘方程(11)两边左乘(Ψi)T得到
和
已知A,B是对称矩阵,将(17)式两边转置得
(18)式减(19)式,并用(14)(15)(16)整理得
在振动理论中,为第i阶特征向量的转置。
展开上式并令:
得到
写成矩阵形式就得到
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)可以获得修正后准确的的刚度、质量、阻尼矩阵K′,M′,C′。
本发明的优点是:实施规范,便于计算机编程实施,由于不要求获知试验与计算模态的一一对应情况,因此适用于密集模态含阻尼结构的模型修正。修正后的模型不仅可用于模态计算,还可用于动力学响应计算分析。
Claims (1)
1.一种密集模态含阻尼结构模型修正方法,其特征在于:
对于初始建立的有限元模型,假设其刚度、质量、阻尼矩阵分别为K、M、C,各矩阵规模为m×m,而各矩阵分别为Ne个元素的单元矩阵叠加构成;
通过对含阻尼结构进行复模态测试,来得到结构的多个特征解,具体可假设包括2N个复特征值及2N个复特征向量其中,r=1,2,…2N;
利用K、M、C构建A,B矩阵,建立状态方程,实现微分方程的降阶,基于复模态理论与矩阵计算,构建出以各元素刚度、质量、阻尼的比例修正系数为未知量的线性方程组,再通过最小二乘法来获得这些修正系数的近似解;
粘性阻尼***的自由振动方程为
在振动理论中,该方程为运动微分方程,u(t)表示位移;
其特征值方程
(λ2Μ+λC+Κ)Φ=0 (2)
在振动理论中,λ表示特征方程的特征根,Φ表示特征方程的特征向量;
引入由位移和速度所组成的2N维状态向量
改写为状态方程形式,此时,方程(1)可写作由状态向量描述的一阶线性微分方程组
在振动理论中,v(t)表示状态向量,由速度和位移组成,v(0)为零时刻的状态向量的值,值的大小为v0;
其中
在振动理论中,u0表示零时刻的位移;
设***在状态空间中的运动为
v(t)=Ψeλt (6)
在振动理论中,Ψ表示状态空间方程的特征向量;
相应的特征值问题为
或
BΨ=λAΨ (8)
将(8)展开与(2)比较,不难看出它们具有相同的特征值λr,并且特征向量满足以下关系
其中,r=1,2,…2N; (9)
在振动理论中,Ψr为状态空间描述的第r阶特征向量,Φr为物理空间描述的第r阶特征向量;
对于含有粘性阻尼结构的分析模型,设Ψi,λi分别为分析模型的第i阶复振型和第i阶复频率,根据(8)式有
BΨi=λiAΨi (10)
同样地假设对应结构的真实的刚度,质量和阻尼矩阵分别为K*,M*,C*,并设分别为“真实模型”的第j阶复振型和第j阶复频率;它们自然也满足
在振动理论中,对于真实结构,A*、B*为的状态空间方程参数,的状态空间方程第j阶特征向量;
其中
假定真实模型与分析模型具有以下修正关系
其中,Κn,Μn,Cn是全局坐标下第n个单元的刚度矩阵,质量矩阵和阻尼矩阵,αn,βn,γn是相应的修正系数,Ne是分析模型中有限单元的个数;
方程(10)两边左乘方程(11)两边左乘(Ψi)T得到
和
已知A,B是对称矩阵,将(17)式两边转置得
(18)式减(19)式,并用(14)(15)(16)整理得
在振动理论中,为第i阶特征向量的转置;
展开上式并令:
得到
写成矩阵形式就得到
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)可以获得修正后准确的的刚度、质量、阻尼矩阵K′,M′,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 CN104008234A (zh) | 2014-08-27 |
CN104008234B true 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) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105956213B (zh) * | 2016-03-07 | 2018-12-21 | 重庆大学 | 一种机敏约束层阻尼板结构模态状态观测器的设计方法 |
CN106682301A (zh) * | 2016-12-23 | 2017-05-17 | 西安电子科技大学 | 一种面向多点激励振动试验的***模型修正方法 |
CN107292045B (zh) * | 2017-06-30 | 2020-07-28 | 南京航空航天大学 | 一种含阻尼结构有限元模型的复模态评估方法 |
-
2014
- 2014-05-12 CN CN201410199353.9A patent/CN104008234B/zh not_active Expired - Fee Related
Non-Patent Citations (1)
Title |
---|
基于模态试验与优化的静动力学模型转换;王轲;《第十届全国振动理论及应用学术会议论文集》;20111231;第591-598页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104008234A (zh) | 2014-08-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103970964B (zh) | 一种挠性卫星模态参数在轨辨识方法 | |
CN105843073B (zh) | 一种基于气动力不确定降阶的机翼结构气动弹性稳定性分析方法 | |
CN102436176B (zh) | 基于神经网络的微陀螺仪控制*** | |
CN111368466B (zh) | 一种基于频响函数参数修正的机械振动预测方法 | |
CN107220403A (zh) | 飞行器弹性模态的控制关联建模方法 | |
CN103886160B (zh) | 一种基于基础激励响应数据的考虑阻尼的模型修正方法 | |
CN110631792B (zh) | 基于卷积神经网络的抗震混合试验模型更新方法 | |
Winter et al. | CFD-based aeroelastic reduced-order modeling robust to structural parameter variations | |
CN104008234B (zh) | 密集模态含阻尼结构模型修正方法 | |
Lind et al. | Robust flutter margin analysis that incorporates flight data | |
Stern et al. | Estimation of dynamic stability coefficients for aerodynamic decelerators using CFD | |
CN107633132A (zh) | 空间分布动态载荷的等效识别方法 | |
CN101794338B (zh) | 基于结构模态试验的矩阵型动力学模型修正方法 | |
CN111159934A (zh) | 一种铝基复合材料动力学模拟方法 | |
Mavriplis et al. | Results using NSU3D for the first aeroelastic prediction workshop | |
Kiviaho et al. | Application of a time-accurate aeroelastic coupling framework to flutter-constrained design optimization | |
CN103902764B (zh) | 基于Householder变换的无约束结构静力分析方法 | |
Thomas et al. | Static/dynamic correction approach for reduced-order modeling of unsteady aerodynamics | |
CN105093933B (zh) | 一种确定lpv变增益控制器的方法 | |
CN115146402A (zh) | 高温合金的氧化过程模拟方法及装置 | |
CN102567561B (zh) | 一种离散化板壳结构的建模方法 | |
CN106342305B (zh) | 一种面向多任务要求的测试性指标确定方法 | |
Callejo et al. | Vehicle suspension identification via algorithmic computation of state and design sensitivities | |
Righi et al. | Uncertainties Quantification of CFD-Based Flutter Prediction | |
CN110298073B (zh) | 集成神经网络与物理***模型的换挡负载动态模拟方法 |
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 |