[发明专利]密集模态含阻尼结构模型修正方法

摘要

本发明公开了一种密集模态含阻尼结构模型修正方法，其主要特点是设定有限元各矩阵的误差是各单元相应子矩阵的误差累积而成，并假设初始的各单元子矩阵与准确矩阵只相差修正因子；在修正算法中首先将动力学方程转换为状态方程，并求特征解；对不同阶的初始状态方程与假设的准确的状态方程的特征方程分别前乘计算与试验得到的不同阶复振型，再相减并整理，得到关于修正因子的线性方程组；用最小二乘法计算修正因子，进而得到修正后的有限元矩阵。该方法实施规范，便于计算机编程实施，由于不要求获知试验与计算模态的一一对应情况，因此适用于密集模态含阻尼结构的模型修正。修正后的模型不仅可用于模态计算，还可用于动力学响应计算分析。

主权项

一种密集模态含阻尼结构模型修正方法，其特征在于：对于初始建立的有限元模型，假设其刚度、质量、阻尼矩阵分别为K、M、C，各矩阵规模为m×m，而各矩阵分别为Ne个元素的单元矩阵叠加构成；通过对含阻尼结构进行复模态测试，来得到结构的多个特征解，具体可假设包括2N个复特征值及2N个复特征向量其中，r＝1,2,…2N；利用K、M、C构建A，B矩阵，建立状态方程，实现微分方程的降阶，基于复模态理论与矩阵计算，构建出以各元素刚度、质量、阻尼的比例修正系数为未知量的线性方程组，再通过最小二乘法来获得这些修正系数的近似解；粘性阻尼系统的自由振动方程为Mu··(t)+Cu·(t)+Ku(t)=0---(1)]]>在振动理论中，该方程为运动微分方程，u(t)表示位移；其特征值方程(λ2Μ+λC+Κ)Φ＝0    (2)在振动理论中，λ表示特征方程的特征根，Φ表示特征方程的特征向量；引入由位移和速度所组成的2N维状态向量v(t)=def[u(t)u·(t)]T---(3)]]>改写为状态方程形式，此时，方程(1)可写作由状态向量描述的一阶线性微分方程组Av·(t)+Bv(t)=0v(0)=v0---(4)]]>在振动理论中，v(t)表示状态向量，由速度和位移组成，v(0)为零时刻的状态向量的值，值的大小为v0；其中A=defCMM0,B=def-K00M.v0=defu0u·0---(5)]]>在振动理论中，u0表示零时刻的位移；设系统在状态空间中的运动为v(t)＝Ψeλt  (6)在振动理论中，Ψ表示状态空间方程的特征向量；相应的特征值问题为(B-λA)Ψ=(-K00M-λCMM0)Ψ~Ψ^=0---(7)]]>或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阶复频率；它们自然也满足B*Ψj*=λj*A*Ψj*---(11)]]>在振动理论中，对于真实结构，A*、B*为的状态空间方程参数，的状态空间方程第j阶特征向量；其中A*=defC*M*M*0---(12)]]>B*=def-K*00M*---(13)]]>假定真实模型与分析模型具有以下修正关系K*=K+Σn=1NeαnKn,---(14)]]>M*=M+Σn=1NeβnMn---(15)]]>C*=C+Σn=1NeγnCn---(16)]]>其中，Κn，Μn，Cn是全局坐标下第n个单元的刚度矩阵，质量矩阵和阻尼矩阵，αn,βn,γn是相应的修正系数，Ne是分析模型中有限单元的个数；方程(10)两边左乘方程(11)两边左乘(Ψi)T得到(Ψj*)TBΨi=λi(Ψj*)TAΨi---(17)]]>和(Ψi)TB*Ψj*=λj*(Ψi)TA*Ψj*.---(18)]]>已知A,B是对称矩阵，将(17)式两边转置得(Ψi)TBΨj*=λi(Ψi)TAΨj*---(19)]]>(18)式减(19)式，并用(14)(15)(16)整理得[ΦiTλiΦiT]-Σn=1NeαnKn00Σn=1NeβnMnΦj*λj*Φj*=[ΦiTλiΦiT](λj*-λi)C+λj*Σn=1NeγnCn(λj*-λi)M+λj*Σn=1NeβnMn(λj*-λi)M+λj*Σn=1NeβnMn0Φj*λj*Φj*---(20)]]>在振动理论中，为第i阶特征向量的转置；展开上式并令：Uij,n=(Φi)TKnΦj*Vij,n=(λj*)2(Φi)TMnΦj*Wij,n=λj*(Φi)TCnΦj*fij=(λi)2(Φi)TMΦj*-(λj*)2(Φi)TMΦj*+λi(Φi)TCΦj*-λj*(Φi)TCΦj*---(21)]]>得到Σn=1NeαnUij,n+Σn=1NeβnVij,n+Σn=1NeγnWij,n=fij---(22)]]>写成矩阵形式就得到Uα+Vβ+Wγ＝f    (23)其中U,V和W是Nm×Ne矩阵；α,β和γ是Ne维列向量，f是Nm维列向量；进一步地，(23)式可以写成Gτ＝f   (24)式中G＝[U V W]   (25)且τ=αβγ.---(26)]]>(24)式为超定方程组，用最小二乘求解就得到修正参数的估计值，然后通过式(14)‑(16)可以获得修正后准确的的刚度、质量、阻尼矩阵K′，M′，C′。