第十章模态综合方法Word文件下载.docx
- 文档编号:22429161
- 上传时间:2023-02-04
- 格式:DOCX
- 页数:73
- 大小:133.20KB
第十章模态综合方法Word文件下载.docx
《第十章模态综合方法Word文件下载.docx》由会员分享,可在线阅读,更多相关《第十章模态综合方法Word文件下载.docx(73页珍藏版)》请在冰豆网上搜索。
:
[K:
]
(10—6)
从而,
[M:
]=[:
T[m:
][〉][M■]丁T[m][J
(10—7)
(10—8)
T寸p}T[M]{p}V弓p}T[K]{p}
当应用拉格朗日方程来建立振动方程时,由于拉格朗日方程要求各独立,而{p}中有不独立的坐标。
u«
l
uH
訥訥uF
由对接位移条件(界面位移连续条件):
{Uj}={u/},有[Y]{p}=[j]{p}
写成约束方程的形式:
[C]{p}=0[C]二[[j]-[訂]
下面进行第二次坐标变换
将{p}分块写成
rI」Da
{p}'
.piJ
则
[[Cdd][Cdi]]」{0}
.Pl‘
{pd}一[Cdd门Cdi]{p}
{口=Pi}=[S]{q}
〔[I]一
心|-[Cdd门Cdi]I
[S]=
[I]一
[S]称为独立坐标变换矩阵。
从而
T占①T[M]{q}V记{q}T[K]{q}
[M]=[S]T[M][S][K]=[S]T[K][S]
由拉格朗日方程可得整个梁结构通过模态综合后的自由振动方程为:
[M]{q}[K]{q}={0}
(10-9)
Pi相互
(10—10)
(10—11)
(10—12)
(10—13)
(10—14)
(10—15)
(10—16)
(10—17)
(10—18)
(10—19)
(10—20)
相应的广义特征值问题为:
([K]一,2[M]){-}={0}(10-21)
其阶数为所有子结构分支模态总数减去界面对接坐标数。
对其进行求解,就可以得到整个梁结构的动力学特性。
对于一般动力学方程,也可以进行上述的变换过程,得到缩减了自由度的动
力学方程:
[M]{q}[C]{q}[K]{q}珂F(t)}
(10-22)
其中:
[C]“s]T[C][S]
(10-23)
{F(t)}“S]T{P(t)}
(10-24)
在模态综合法中,为了描述结构在空间的运动和变形状态,采用两类广义坐标来描述,分别为“物理(几何)坐标”和“模态坐标”,物理坐标描述结构各节点的几何坐标位置,而模态坐标则表示物理坐标响应中各个模态成份大小的量。
对于模态综合法中的“模态”一词,它比“振型”具有更加广义的内涵,它不仅指结构做主振动时的振型,而且还包括了结构在一些特定的外力或者结点位移作用下产生的静变形形态,这些静变形形态被认为是在整个结构振动时,各子结构可能产生的变形形态。
而“振型”则是一个狭义的概念,表示结构作主振动时的变形形式。
【模态综合法的基本步骤】
由上例可以看到,模态综合法的基本步骤可以分成如下六个步骤:
(1)按结构特点划分子结构
(2)计算并选择分支模态进行第一次模态坐标变换
(3)在全部模态坐标中,选择不独立的广义坐标
(4)由位移对接条件,形成广义坐标的约束方程,得到独立坐标变换阵[S]
(5)对组集得到的质量矩阵、刚度矩阵进行合同变换,得到独立坐标下的质量
矩阵,刚度矩阵,形成整个系统的振动方程
(6)根据坐标变换关系,再现子结构物理参数
由上可知,模态综合法的关键技术是如何选择子结构的分支模态。
0.2各种形式的分支模态
如前所述,分支模态就是在结构系统振动时,其子结构(分支结构)可能出现的变形形态。
在模态综合法中,分支结构分为两类:
受约束分支结构、有刚体运动的分支结构。
有刚体运动的分支结构又称为自由悬浮分支。
一、受约束子结构的分支模态
它的可能变形形态包括:
在各种附加约束或无附加约束下自由振动模态,在各种外力作用下的位移形
态,在各种给定的边界条件下的内部位移形态。
在进行模态综合时,只需要选其中一部分构成其分支模态,且各有其相应的名称。
【主模态】
分支主模态由下列子结构的特征方程决定:
([約-co2[ma]){©
a}={0}(10—25)
在确定分支主模态时,需要首先确定子结构的界面坐标处理状态,按照对界
面位移的处理方法,有三种分支主模态
固定界面主模态:
子结构的全部界面加上附加约束
自由界面主模态:
子结构的全部界面都没有附加约束,但子结构本身原有的约束(称为自然约束)仍然存在
混合界面主模态:
子结构的部分界面加上附加约束
在模态综合法中,假定主模态阵都已按质量归一化。
即:
[:
]丁[何][:
]二[Inn](10—26)
[:
]丁[約[:
]=diag[/](10—27)
如果模态综合法所使用的不是子结构的完全主模态矩阵,而是保留主模态集,即经过高阶模态截断后的部分低阶主模态,模态综合法的误差就由此而产生。
【约束模态】
约束模态是指对界面坐标的约束模态,它定义为:
在子结构的全部界面自由度上引入附加约束,然后让这些界面自由度依次产生单位位移,其它约束(包括自然约束和附加约束)则保持不变(即这些界面坐标都强制为零)。
由此产生的一系列子结构静变形位移,称为子结构对于界面坐标的约
1
自由界面
朿模态,简称约朿模态。
约朿模态的数目,等于界面自由度的数目,全部约朿模态就组成子结构的约束模态阵N;
],从约束模态的生成过程看到,它有点类
似于有限元法中的形函数。
显然,约束模态可以写成:
(10—28)
下标V表示子结构不受约束的自由度,C表示附加约束的自由度,[I)为单位阵,
表示界面坐标依次产生单位位移。
卜vc]为子结构内部坐标由于界面坐标依次有
单位位移时所产生的静态位移。
要让界面坐标依次产生单位位移,必须对界面
坐标施加一定的界面力,记界面力矩阵为[尺訂,则应有:
■[kvv][kvc]『vcl=[[Ovc]【
1kcv][kcc讣IcdJRcc].
(10—29)
分块展开第一行有:
「vcH-[kvv]j[kvc]
(10—30)
从而约束模态为:
rvc]〕=]-[kvv]巾kvc]IJIcc]j:
[Icc]一
(10—31)
如图所示为悬臂梁的约束模态示意图。
Vj
有=0
Vi-1
【附着模态】
模态综合法中的附着模态是对界面坐标的附着模态。
定义为:
对子结构的界
面不附加任何约束,而是在一个界面自由度上沿此自由度方向施加单位力,而其它自由度上无外力作用,由此得到的子结构静态位移向量,就是子结构对该界面自由度的附着模态。
显然这个定义只适合于受约束子结构。
依次在每个界面自由度上作用单位力,就可以得到一系列静态位移,也就构成子结构对其界面坐标的附着模态矩阵口Jo
根据附着模态的定义,附着模态矩阵口;
]由下式来确定:
[kww][kwa][-wa][°
wa]
」kaw][kaa][aa]][laa]
(10—32)
从而有:
F:
wa]Lk]Ji[0wa]L【g]严a]l[[屮aa]」][IaaL][GL
(10—33)
子结构的柔度矩阵为:
[gwa]
[gaa]
(10—34)
所以有:
卜wa]=gww]
「aa]=_[gaw]
[gwa]「[0wa]1
[gaa]』[laal
(10—35)
从而附着模态为:
(10—36)
『wa]]=]gwa]]恠aa/haa]]
=1
Pj"
【剩余附着模态】
在假设模态法建立系统的运动方程,求解其特征值问题时,要求所用到的假
设模态应该是线性无关的,但是如果用子结构的主模态和附着模态作为假设模
41
态集,会出现主模态与附着模态线性相关的问题。
另一方面,在使用子结构的主模态组成模态综合时的假设模态集,采用的是
经过高阶截断的主模态。
显而易见,如果在假设模态集中加入高阶主模态的信息,则可以提高模态综合的精度。
前面提到,附着模态与保留主模态线性相关。
如果从附着模态中减去与之不独立的低阶保留主模态,则可以得到高阶主模态的近似表达。
这种模态称为“剩余附着模态”
设受约束子结构的全部归一化主模态为
(10-37)
[d]为高阶主模态,即剩余主模
g]=Wk][%]]=[:
:
其中,「k]为低阶主模态,即保留主模态,态。
下标i,j分别表示内部自由度和界面自由度。
显然有
(10—38)
[:
」]T[K][:
」]二diag[2]
子结构柔度矩阵为:
(10—39)
[G]珂K]」=[:
」](diag[・2])」冷]T
=[%](diag®
2])」[%]T+[%](diag^d])」此]T
定义剩余柔度矩阵:
G]珂G]-[k](diag[初]『珂d](diag[•;
])[】]丁(10-40)
仿照(10-32)的定义,将[G]换成[Gd],得到剩余附着模态阵的定义式为:
(10—41)
【5]=[Gd]
F面讨论一下受约束子结构剩余附着模态的物理意义
由定义:
化戶[Gdl.ftZdKdiag[研]心d]Tf叮(10—42)
]一]一
而
[灯『]]=[%「[%d]T罟]1=叫]丁(10—43)
」1j]-」1j]-
由此,剩余附加约束模态可以写成:
[誓d]=[%](diag固訂円编]丁=[札]((diaQa訂)」[怀]丁)(10—44)
显然,[丫d]与[d]线性相关。
也就是说,剩余附着模态实际是进行主模态截断时,略去的高阶模态的一种线性组合。
因此用剩余附着模态作为子结构分支模态集的一个子集,是对保留主模态集的一个合理补集。
由于[d]与保留主模态[\]具有正交性,因而剩余附着模态的[讥]与[k]也是关于质量阵正交的。
[$[皿][讥]=[0](10—45)
二•有刚体运动子结构的分支模态
对于有刚体运动的子结构,其模态集中包含有全部的刚体模态,由此可知,得到的有刚体运动子结构的刚度矩阵是奇异的。
有刚体运动子结构的主模态的定义与受约束子结构的主模态定义相同,只是
还应包括相应与刚体自由度的刚体模态。
但是对于固定界面和混合界面的分支主模态中,如果附加界面约束消除了刚体自由度,则这时的分支主模态集中将没有刚体模态。
有刚体运动子结构的约束模态的定义与受约束子结构的约束模态定义相同,但是,只有在附加界面约束全部约束了子结构的全部运动时,才能求出约束模^态0
胃1「叽广[心丨(10—46)
-[I』一
【刚体模态】
对受不完全约束的子结构,即有刚体自由度的子结构,描述其无变形运动位移的模态称为刚体模态。
对于空间自由悬浮结构,最多只有六个刚体模态坐标,因此刚体模态数r满足0空r乞6。
对于有刚体位移的子结构,其刚体模态是十分重要的,其分支模态集中必须包含这些刚体模态。
刚体模态可以通过求解自由界面子结构主模态的特征方程得到。
也可以作为子结构约束模态的一种特殊情况求出,即当附加界面约束刚好约束住结构的全部刚体自由度时,这是求出的
约束模态,就是有刚体运动自由度的子结构的刚体模态
叫]1=]%]」["
|[lrr]「[[Irr]_
当界面附加约束超过了子结构的刚体自由度时,约束模态就是弹性位移和全部刚体位移的线性组合。
即这时的约束模态实际上包含了刚体模态。
【例】求图示系统的约束模态
一4一乞——虫
mkmkmx
显然系统有一个刚体模态,其刚度矩阵为:
故%]二k%]呱]T
k[-1
0][kccHk21
-n
-
'
-1
01
[K]=k
2
1一
取约束坐标集C二[U2U3],由方程(10—30)得到:
Vc]1
_|-[kw]J[kvcfi
「1
」lcc]-
[[Ice]-
[‘讥]二
01
如果取界面坐标集为:
C=(u3),此时约束模态数和刚体模态数相等,约束模态退化为刚体模态。
即
C[111]T
对有刚体运动的子结构,应该在子结构中引入适当的附加静定约束,并刚好能约束住子结构的刚体运动,然后才能求其附着模态。
将有刚体运动子结构的物理坐标分为三个子集,即U=R•A,W。
R为附加约束坐标集,它刚好约束住子结构的刚体位移,A二C-R,C为界面坐标集,W为内部坐标集。
根据附着模态的定义,在A集中的物理坐标上依次作用单位力,得到子结构的位移就是附着模态。
[kwr][<
wa]
%][屮aa]=
[krr]_[Ora]
[Owa]||
[Iaa]
[Rra]
[Rra]是在附加约束坐标集
R中产生的附加约束反力。
将(
10—48)分块展
开,得到:
.Kww]
[kaw]
(10—49)
•bw]
[kwa]
[gww]
[gwa]
[kaw]
[kaa]
|_[gaw]
(10—50)
砂wall
[gww]
贋aa]」Jgaw]
[Owal_[
[gaa]_||[laa]「[[gaa]」
gwa]I
(10—51)
从而有刚体自由度子结构的附着模态为:
『wa]1序a]=haa]
_[Ora]
•gwa]I
.gaa]
[Ora]J
(10—52)
图中[子结构是一个有刚体自由度的子结构,引入附加约束后,可以求出其
附着模态,如图所示
pm
【惯性释放附着模态】
女口前面所述,我们看到,对有刚体自由度子结构,弓I入附加静定约束后,在界面坐标方向施加单位作用力,得到附着模态。
而惯性释放附着模态是在释放这些附加约束后,在界面单位力作用下求得的附着模态。
但这时子结构会产生刚体位移,因此我们将一组相应的惯性力作用在子结构上,这些惯性力与界面力组成一个自平衡力系。
惯性释放附着模态就是该平衡力系下,子结构的不包含刚体模态的附着模态。
在界面作用力作用下,子结构的位移矢量{u(t)}是刚体位移分量{Ur(t)}和变
形位移分量{uE(t)}之和。
{U(t)}二{Ur(t)}{UE(t)}(10-53)
代入子结构运动方程
[M]{u}[K]{u}二{F}(10-54)
由于[K]{Ur(t)}=0,故得到:
[M]{Ue}[K]{Ue}二{F}-[M]{Ur}={Fe}(10-55)
{Fe}就是上面所说的由界面力与惯性力组成的自平衡力系。
假定已经求出了子结构的主模态
(10—56)
[「]=[[「][e]]
[r]与[e]分别为子结构的刚体模态矩阵和完全弹性变形模态矩阵。
那么有:
(10—57)
{Ur(t)}<
r]{qr(t)}{UE(t)}二[E]{qE(t)}
于是,
(10—58)
对方程(10-55)作坐标变换,并注意到[■]与[e]的正交性,得到:
[mr]{qr}=[r]T{F}
(10—59)
[mE]{qE}[kE]{qE}珂e]T{F}(10-60)
其中,
[mr]<
r]T[M][r]He]二[e]T[M][e][KeH[e]T[K][e]
(10—61)
由(10-59)得到:
{qr}=[mr]Tr]T{F}
(10—62)
所以
{Ur}珂r][mr]'
[r]T{F}(10-63)
{Fe}二{F}-[M]{Ur}=([I]-[M][r][mr]TF){F}=[R]{F}(10-64)
显然,用[R]矩阵对界面力矩阵进行线性变换,就可以得到子结构在界面力
和惯性力组成的自平衡力系。
如果主模态已经对质量阵归一化,则
[R]<
I]_[M][r][r]T
自平衡力系作用下的附着模态{~a}由下列方程确定:
上式中,
[K]{~J=[Fe]=[R][F]
[0wa]
[F]=l[Iaa]
[0raI
为界面作用力矩阵
[kww]
[krw]
[kwa]
[kaa]
[kra]
[kwr][~wa]
[kar]l[~aa]
[krr]」~ral
0wa]【
[Iaa]
J0ral
由此得到:
[~wa11
[®
a]=|[~aa]=[G][R][F]
[0raj
.gww][gwa]
[gaw]
J0rw]
[0ra]
[0wr]1
[0ar]
[0rr]
般情况下,
[~a]中含有子结构的刚体模态分量,为了消除刚体模态
[a]珂~a「[r][G]
[Cr]由下式决定:
[a]珂~a]-[川殆
从而:
[Cr]"
mrlTJ[M][~a]=[J[M][~a]
(10—67)
(10—68)
(10—69)
(10—70)
(10—71)
.令:
(10—72)
(10—73)
(10—74)
所以,惯性释放附着模态的计算公式为:
[a]=([l]-[r][J[M])[~a]珂R]T[~a]珂R]T[G][R][f]
【例】对图示系统,设U3为约束坐标,
Ui为内部坐标,
(1)确定系统的附着模态
(2)确定系统的自平衡力系
(3)确定系统的惯性释放附着模态
U2是附着坐标,
系统的刚度阵和质量阵为:
■2
[K]=k-1
■-1
■1
[M]=m0
00
(1)根据刚度矩阵,可以求得系统的柔度矩阵,
附着模态为:
Lwa]
gww
_gaw
[gwa]I
gwa2k
d1=ll,
gaa
{^a}=『aa]=paa]=化
=L_.
][0丄][0ra]_
k;
12'
(2)系统有一个刚体模态
{小[1
11]T
归一化后得到:
[111]T
自平衡力系为:
[0wa][
{Fe}=[R][F]=([I]-[M][®
]®
]T)[laa]
J°
raL
=(0
■0
_2
0m0
-0
T「01一1
2]L0J
自平衡力系示意图如图
m
■-11
l_T_|
(3)系统的惯性释放附着模态为:
[V;
]“R]T[G][R][F]
/2
-1-1〕
0〕
_r
Li
2,
3
2一
lk-
0一
9k
;
1|
l一化
在模态综合法中,常常使用的是分支的惯性释放附着模态而不是分支的
附着模态。
(整体结构振动时,惯性释放附着模态是子结构更加可能的变形
形态)
【剩余惯性释放附着模态】
定义系统的弹性柔度阵[Ge]为
[Ge]二[E][kE「[E]T(10-76)
[kE]为系统对模态坐标的刚度矩阵。
[kEp[E]T[K][E]=diag[2](10-77)
[Ge]二[e][〒[e]T(10-78)
将完全弹性变形主模态[e]分成保留主模态[\]和剩余主模态[d],即
[eP[k][d]】(10-79)
[Ge]=[k][2]'
[k]T[d][2]'
[d]T(10-80)
定义剩余柔度阵为:
[Gd]=[d][京Td]T"
Ge]-[k][$]Tk]T(10-81)
根据(10-65)式与(10-70)式可得:
T〜
[R][G][R]二[Ge](10-82)
故
[Gd円R]t[(〜][R]-[加:
]平『(10-83)
类似惯性释放附着模态,定义剩余惯性释放附着模态:
[一a]珂Gd][F]
以上介绍了有关各种分支模态的概念,在以后的各种模态综合法中会看
到各种模态综合方法的差异,就在于模态集的不同选取。
0・3固定界面模态综合法
子结构模态综合法的关键,是对子结构进行模态坐标变换
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 第十 章模态 综合 方法
![提示](https://static.bdocx.com/images/bang_tan.gif)