CN102880803A - 一种复杂机械结构的转动自由度频率响应函数计算方法 - Google Patents
一种复杂机械结构的转动自由度频率响应函数计算方法 Download PDFInfo
- Publication number
- CN102880803A CN102880803A CN2012103716461A CN201210371646A CN102880803A CN 102880803 A CN102880803 A CN 102880803A CN 2012103716461 A CN2012103716461 A CN 2012103716461A CN 201210371646 A CN201210371646 A CN 201210371646A CN 102880803 A CN102880803 A CN 102880803A
- Authority
- CN
- China
- Prior art keywords
- frequency response
- measuring point
- response function
- centerdot
- freedom
- 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
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明涉及一种基于响应耦合技术的转动自由度频率响应函数估计方法,主要适用于复杂机械结构中转动自由度频率响应函数的估计。本发明利用响应耦合(Receptance coupling)技术,首先根据需要估计的转动自由度频率响应函数,将复杂机械结构分解为子结构A和子结构B,并在子结构A上便于测量的位置选择第一测点,在子结构A和B的结合面处选择第二测点;然后利用锤击激励法测量三个平动自由度的频率响应函数,求解子结构A和子结构B结合面处的频率响应函数矩阵;最后可计算得到两个测点处与转动自由度有关的所有频率响应函数。该发明基于响应耦合技术,实施方便、计算结果准确,为转动自由度频率响应函数的估计提供了又一有效的技术。
Description
技术领域
本发明属于复杂机械结构的分析技术领域,涉及一种复杂机械结构的转动自由度频率响应函数计算方法。
背景技术
广泛应用的有限元技术,分析简单机械结构时可以达到很高的精度,然而用其对复杂机械结构进行建模和动力学分析时,很难得到令人满意的结果。原因之一就是人们对复杂结构中各子***间的耦合关系认识不清,在建模时常常对其进行不合理的简化,导致建立的模型精度不高,难以匹配实验结果。因此对于复杂机械结构的有限元模型,需要进行修正以提高精度。
高质量的频率响应函数是对机械结构的有限元模型进行成功修正的基础。工程中通常利用试验模态测试法测量结构的频率响应函数,由于实验条件限制,一般只能测量机械结构平动自由度的频率响应函数。对于转动自由度,由于角位移难以测量而且代价高昂,其频率响应函数难以利用试验直接获得。因此,迫切需要一种能准确估计转动自由度频率响应函数的方法。
目前估计转动自由度频率响应函数的方法主要有Yoshimura(YoshimuraT,Hosoya N.FRF estimation on rotational degrees of freedom of structures[J].Proceedings of the International Modal Analysis Conference-IMAC,2000,2:1667-1671.Yoshimura T,Hosoya N.结构转动自由度频率响应函数的估计[J].国际模态分析会议论文集,2000,2:1667-1671.)提出的T型块法(T-blockapproach)。该方法对T型块的安装要求很高,实施不便,计算过程复杂,精度较低。
响应耦合(Receptance Coupling)技术是对复杂***或结构的动态特性进行求解的一种方法(RenY,Beards CF.On Substructure Synthesis with FRFData[J].Journal of Sound and Vibration,1995,185(5):845-866.Ren Y,BeardsCF.利用FRF数据进行子结构综合.声振学报,1995,185(5):845-866)。该理论中,将复杂的***分解为若干个简单子结构,分别用解析法或实验法得到各个子***的频率响应函数,再根据子***之间的耦合关系,即共同边界上的平衡条件(Equilibrium condition)及相容条件(Compatibility condition)对子***进行合成,最终求得总***的动态响应。
发明内容
本发明解决的问题在于提供一种复杂机械结构的转动自由度频率响应函数计算方法,该方法是基于响应耦合技术,实施方便、计算结果准确,为复杂机械结构提供准确可靠的转动自由度频率响应函数估计方法。
本发明是通过以下技术方案来实现:
一种复杂机械结构的转动自由度频率响应函数计算方法,包括以下步骤:
1)根据被测量转动自由度频率响应函数的激励点和响应点,将待分析的复杂机械结构分解为子结构A和子结构B,子结构A能够用有限元方法准确建模,在子结构A上选择出第一测点,在子结构A和子结构B的结合面处选择出第二测点;
2)利用锤击激励法测量三个平动自由度的频率响应函数:
第一测点处平动自由度的原点频率响应函数g11,ff,第二测点处平动自由度的原点频率响应函数g22,ff,第一测点与第二测点之间平动自由度的频率响应函数g12,ff;
其中,g11,ff为激励点为第一测点,响应点为第一测点时测得的平动自由度频率响应函数;g12,ff为激励点为第二测点,响应点为第一测点时测得的平动自由度频率响应函数;g22,ff为激励点为第二测点,响应点为第二测点时测得的平动自由度频率响应函数;
3)子结构A在自由状态下的所有频率响应函数利用有限元模型得到数值解,然后利用响应耦合技术求解子结构A和子结构B结合面处的耦合频率响应函数矩阵H2;
4)由结合面处耦合频率响应函数矩阵H2计算得到第一测点处与转动自由度有关的频率响应函数,第二测点处与转动自由度有关的频率响应函数,第一测点与第二测点之间与转动自由度有关的频率响应函数。
所述的转动自由度频率响应函数为待分析的复杂机械结构上任意两个激励点和响应点之间与转动自由度有关的频率响应函数;根据激励点和响应点,选定第一测点和第二测点。
所述的子结构A和子结构B通过阻尼、转动刚度和平动刚度结合。
若转动自由度频率响应函数的激励点和响应点为不同测点,则在子结构A上选择出第一测点,在子结构A和子结构B的结合面处选择出第二测点;如果子结构A上选的第一测点是激励点,则子结构A和子结构B的结合面处选择的第二测点为响应点;如果子结构A上选的第一测点是响应点,则子结构A和子结构B的结合面处选择的第二测点为激励点;
若转动自由度频率响应函数的激励点和响应点为同一测点,则该测点即为子结构A上的第一测点,仍在子结构A和子结构B的结合面处选择出第二测点。
所述步骤2)利用锤击激励法测量三个平动自由度的频率响应函数时,激励点采用激振力锤来进行锤击,响应点利用加速度传感器来检测加速度振动响应信号,由信号采集***计算分析。
所述求解子结构A和子结构B结合面处的频率响应函数矩阵H2为:
假设测点在一个平面内的运动由平动和转动自由度组成,输入力F是由力f和力矩M组成的向量,输出响应X由平动位移x和转动位移θ组成,输入力与输出响应的关系为
(1)式中, 其中hij,ff为测点j的平动自由度与测点i的平动自由度之间的频率响应函数,hij,fM为测点j的转动自由度与测点i的平动自由度之间的频率响应函数,hij,Mf为测点j的平动自由度与测点i的转动自由度之间的频率响应函数,hij,MM为测点j的转动自由度与测点i的转动自由度之间的频率响应函数;
在复杂机械结构中第一测点处施加外力F1,只考虑第一测点和第二测点处的响应X1和X2,得到复杂机械结构的频率响应函数矩阵G11和G21如下:
(2)式中,HA,11为子结构A中第一测点的原点频率响应函数矩阵,HA,12为子结构A中第一测点和第二测点之间的频率响应函数矩阵,HA,21为子结构A中第二测点和第一测点之间的频率响应函数矩阵,HA,22为子结构A中第二测点的原点频率响应函数矩阵,H2为结合面处耦合频率响应函数矩阵;
HA,11,HA,12,HA,21,HA,22和H2的矩阵形式表达式为:
G11和G21的矩阵表达式为:
仅在第二测点处施加外力F2,可得频率响应函数矩阵G12和G22:
上式中,G12和G22的矩阵表达式为:
将G11,G21,G12和G22用频率响应函数矩阵形式表示如(4)式:
分别取频率响应函数矩阵G11,G12和G22中的第一个元素,得方程组式(5):
式(5)为一个含有4个未知数h2,ff,h2,fM,h2,Mf和h2,MM的方程组,子结构A在自由状态下的所有频率响应函数HA,11,HA,12,HA,21和HA,22利用有限元模型得到数值解;g11,ff,g12,ff和g22,ff为复杂机械结构上平动自由度的频率响应函数,由步骤2)得到;利用式(5)解出结合面处耦合频率响应函数矩阵H2。
所述将结合面处耦合频率响应函数矩阵H2代入式(4)中,计算得到第一测点处与转动自由度有关的频率响应函数g11,fM、g11,Mf和g11,MM,其中:g11,fM为第一测点处转动自由度与平动自由度之间的频率响应函数,g11,Mf为第一测点处平动自由度与转动自由度之间的频率响应函数,g11,MM为第一测点处转动自由度的原点频率响应函数;
第二测点处与转动自由度有关的频率响应函数g22,fM、g22,Mf和g22,MM,其中:g22,fM为第二测点处转动自由度与平动自由度之间的频率响应函数,g22,Mf为第二测点处平动自由度与转动自由度之间的频率响应函数,g22,MM为第二测点处转动自由度的原点频率响应函数;
第一测点和第二测点之间与转动自由度有关的频率响应函数g21,fM、g21,Mf、g21,MM、g12,fM、g12,Mf和g12,MM,其中:g21,fM为第一测点处转动自由度与第二测点处平动自由度之间的频率响应函数,g21,Mf为第一测点处平动自由度与第二测点处转动自由度之间的频率响应函数,g21,MM为第一测点处转动自由度与第二测点处转动自由度之间的频率响应函数,g12,fM为第二测点处转动自由度与第一测点处平动自由度之间的频率响应函数,g12,Mf为第二测点处平动自由度与第一测点处转动自由度之间的频率响应函数,g12,MM为第二测点处转动自由度与第一测点处转动自由度之间的频率响应函数。
所述的第一测点处转动自由度的原点频率响应函数g11,MM、第二测点处转动自由度的原点频率响应函数g22,MM,第一测点与第二测点之间转动自由度频率响应函数g21,MM,分别表示如下:
与现有技术相比,本发明具有以下有益的技术效果:
本发明提供的复杂机械结构的转动自由度频率响应函数计算方法,将给定机械结构分解为两个子结构,即子结构A和子结构B,从而利用子结构A的限元方法准确建模,以及子结构A和子结构B的响应耦合技术进行复杂机械结构的转动自由度频率响应函数的计算。
本发明提供的复杂机械结构的转动自由度频率响应函数计算方法,其实施方便,代价低廉:该方法只需利用锤击激励法测量三个平动自由度的频率响应函数,即第一测点处平动自由度的原点频率响应函数g11,ff,第二测点处平动自由度的原点频率响应函数g22,ff,第一测点与第二测点之间平动自由度的频率响应函数g12,ff,就可以对这两个测点间与转动自由度有关的所有频率响应函数进行计算,避免了对角位移的测量,因此实施方便,代价低廉。
本发明提供的复杂机械结构的转动自由度频率响应函数计算方法,由于采用了先进的响应耦合(Receptance Coupling)技术,计算结果准确。
本发明提供的复杂机械结构的转动自由度频率响应函数计算方法,可以应用于复杂机械结构动态特性辨识、动力学建模与仿真、有限元模型修正等方向。
附图说明
图1是复杂机械结构分解为子结构A和子结构B的示意图;
图2是悬臂梁结构拆分示意图;
图3是悬臂梁平动自由度频率响应函数测试示意图;
其中1为第一测量点,2为第二测量点;
图4-1是测量得到的平动自由度频率响应函数g11,ff,图4-2是测量得到的平动自由度频率响应函数g12,ff,图4-3是测量得到的平动自由度频率响应函数g22,ff;
图5是子结构A的有限元模型;
图6-1是估计得到的转动自由度频率响应函数g11,MM,图6-2是是估计得到的转动自由度频率响应函数g21,MM,图6-3是估计得到的转动自由度频率响应函数g22,MM。
具体实施方式
下面结合具体的实施例对本发明做进一步的详细说明,所述是对本发明的解释而不是限定。
参见图1,将待分析的复杂机械结构分解为子结构A和子结构B,子结构A能够用有限元方法准确建模,在子结构A上选择出第一测点,在子结构A和B的结合面处选择出第二测点;
进一步的,所述的子结构A和子结构B之间的通过阻尼、转动刚度和平动刚度结合。
谈到频率响应函数,一定要指出该频率响应函数对应的激励点和响应点。可以根据激励点和响应点,去选定第一测点和第二测点。若转动自由度频率响应函数不同,选择的测点也不同。这样一来,第一测点、第二测点与转动自由度有关的频率响应函数代表复杂机械结构的转动自由度频率响应函数了。
所述的转动自由度频率响应函数为待分析的复杂机械结构上任意两个激励点和响应点之间与转动自由度有关的频率响应函数;根据激励点和响应点,选定第一测点和第二测点。
具体的若转动自由度频率响应函数的激励点和响应点为不同测点,则在子结构A上选择出第一测点(可为激励点或响应点),在子结构A和B的结合面处选择出第二测点(可为激励点或响应点);如果子结构A上选的第一测点是激励点,则子结构A和子结构B的结合面处选择的第二测点为响应点;如果子结构A上选的第一测点是响应点,则子结构A和子结构B的结合面处选择的第二测点为激励点;
若转动自由度频率响应函数的激励点和响应点为同一测点,则该测点即为子结构A上的第一测点,仍在子结构A和B的结合面处选择出第二测点。
具体对悬臂梁结构的转动自由度频率响应函数进行估计,包括以下步骤:
1)根据被测量转动自由度频率响应函数的激励点和响应点,将待分析的复杂机械结构分解为子结构A和子结构B,子结构A能够用有限元方法准确建模,在子结构A上选择出第一测点,在子结构A和B的结合面处选择出第二测点;
参见图2悬臂梁结构示意图,将其分为如图所示的子结构A和子结构B,并选择第一测点1和第二测点2;在结合点处,将总结构分解为子结构A和子结构B;在子结构A的端部,选择第一测点,在结合面处选择第二测点;
2)利用锤击激励法测量三个平动自由度的频率响应函数:
第一测点处转动自由度的原点频率响应函数g11,MM,第二测点处转动自由度的原点频率响应函数g22,MM,第一测点与第二测点之间转动自由度的频率响应函数g21,MM;
参见图3,利用锤击激励法测量三个平动自由度的频率响应函数,具体步骤如下:
第一测点处平动自由度的原点频率响应函数g11,ff,利用力锤敲击第一测点,施加脉冲激励力,在第一测点处利用振动传感器拾取振动响应信号;
第二测点处平动自由度的原点频率响应函数g22,ff,利用力锤敲击第二测点,施加脉冲激励力,在第二测点处利用振动传感器拾取振动响应信号;
第一测点与第二测点之间平动自由度频率响应函数g12,ff,利用力锤敲击第二测点,施加脉冲激励力,在第一测点处利用振动传感器拾取振动响应信号
测量过程所使用的仪器型号为:激振力锤是美国PCB公司生产的086C03型ICP激振力锤,振动传感器是美国PCB公司生产的333B32型ICP加速度传感器,信号采集***软件杭州亿恒公司生产的AVANT数据采集***。
根据采集的激励力信号和振动响应信号,利用AVANT数据采集***的模态测试模块,计算频率响应函数。
参照图4-1~4-3所示,试验得到第一测点处平动自由度的原点频率响应函数g11,ff,第二测点处平动自由度的原点频率响应函数g22,ff,第一测点与第二测点之间平动自由度频率响应函数g12,ff
3)求解子结构A和子结构B结合面处的耦合频率响应函数矩阵H2。
假设测点在一个平面内的运动由平动和转动自由度组成,输入力F是由力f和力矩M组成的向量,输出响应X由平动位移x和转动位移θ组成,输入力与输出响应的关系为
式中, 其中hij,ff为测点j的平动自由度与测点i的平动自由度之间的频率响应函数,hij,fM为测点j的转动自由度与测点i的平动自由度之间的频率响应函数,hij,Mf为测点j的平动自由度与测点i的转动自由度之间的频率响应函数,hij,MM为测点j的转动自由度与测点i的转动自由度之间的频率响应函数;
在复杂机械结构中第一测点处施加外力F1,只考虑第一测点和第二测点处的响应X1和X2,得到复杂机械结构的频率响应函数矩阵G11和G21如下:
式中,HA,11为子结构A中第一测点的原点频率响应函数矩阵,HA,12为子结构A中第一测点和第二测点之间的频率响应函数矩阵,HA,21为子结构A中第二测点和第一测点之间的频率响应函数矩阵,HA,22为子结构A中第二测点的原点频率响应函数矩阵,H2为结合面处耦合频率响应函数矩阵;
HA,11,HA,12,HA,21,HA,22和H2的矩阵形式表达式为:
G11和G21的矩阵表达式为:
仅在第二测点处施加外力F2,可得频率响应函数G12和G22:
上式中,G12和G22的矩阵表达式为:
将G11,G21,G12和G22用频率响应函数矩阵形式表示如(4)式:
分别取频率响应函数矩阵G11,G12和G22中的第一个元素,得方程组(5):
h2,fM=h2,Mf
式(5)为一个含有4个未知数h2,ff,h2,fM,h2,Mf和h2,MM的方程组,子结构A在自由状态下的所有频率响应函数HA,11,HA,12,HA,21和HA,22利用有限元模型得到数值解;g11,ff,g12,ff和g22,ff为复杂机械结构上平动自由度的频率响应函数,由步骤2)得到;利用式(5)解出结合面处耦合频率响应函数矩阵H2。
将结合面处耦合频率响应函数矩阵H2代入式(4)中,计算得到第一测点处与转动自由度有关的频率响应函数g11,fM、g11,Mf和g11,MM,其中:g11,fM为第一测点处转动自由度与平动自由度之间的频率响应函数,g11,Mf为第一测点处平动自由度与转动自由度之间的频率响应函数,g11,MM为第一测点处转动自由度的原点频率响应函数。
第二测点处和转动自由度有关的频率响应函数g22,fM、g22,Mf和g22,MM,其中:g22,fM为第二测点处转动自由度与平动自由度之间的频率响应函数,g22,Mf为第二测点处平动自由度与转动自由度之间的频率响应函数,g22,MM为第二测点处转动自由度的原点频率响应函数。
第一测点与第二测点之间与转动自由度有关的频率响应函数g21,fM、g21,Mf、g21,MM、g12,fM、g12,Mf和g12,MM,其中:g21,fM为第一测点处转动自由度与第二测点处平动自由度之间的频率响应函数,g21,Mf为第一测点处平动自由度与第二测点处转动自由度之间的频率响应函数,g21,MM为第一测点处转动自由度与第二测点处转动自由度之间的频率响应函数,g12,fM为第二测点处转动自由度与第一测点处平动自由度之间的频率响应函数,g12,Mf为第二测点处平动自由度与第一测点处转动自由度之间的频率响应函数,g12,MM为第二测点处转动自由度与第一测点处转动自由度之间的频率响应函数。
具体的,参考图5所示,对于子结构A,可以用Timoshenko梁单元进行有限元建模,将其划分为10个单元,每个节点分别包含6个自由度,即3个平动(δx,δy,δz)和3个转动(γx,γy,γz)自由度。子结构A的边界条件为自由状态。
利用子结构A的有限元模型,仿真频率响应函数hA11,ff,hA21,ff,hA21,Mf,hA22,ff,hA22,Mf,hA22,fM,hA12,ff,hA12,fM,hA12,fM;然后将上述频率响应函数代入下面方程组中,
h2,fM=h2,Mf
求解得到h2,ff,h2,fM,h2,Mf和h2,MM,即可得到频率响应函数矩阵H2。
将h2,ff,h2,fM,h2,Mf和h2,MM代入下面方程中,
参照图6-1~6-3所示,计算得到第一测点处转动自由度的原点频率响应函数g11,MM、第二测点处转动自由度的原点频率响应函数g22,MM,第一测点与第二测点之间转动自由度频率响应函数g21,MM。
Claims (8)
1.一种复杂机械结构的转动自由度频率响应函数计算方法,其特征在于,包括以下步骤:
1)根据被测量转动自由度频率响应函数的激励点和响应点,将待分析的复杂机械结构分解为子结构A和子结构B,子结构A能够用有限元方法准确建模,在子结构A上选择出第一测点,在子结构A和子结构B的结合面处选择出第二测点;
2)利用锤击激励法测量三个平动自由度的频率响应函数:
第一测点处平动自由度的原点频率响应函数g11,ff,第二测点处平动自由度的原点频率响应函数g22,ff,第一测点与第二测点之间平动自由度的频率响应函数g12,ff;
其中,g11,ff为激励点为第一测点,响应点为第一测点时测得的平动自由度频率响应函数;g12,ff为激励点为第二测点,响应点为第一测点时测得的平动自由度频率响应函数;g22,ff为激励点为第二测点,响应点为第二测点时测得的平动自由度频率响应函数;
3)子结构A在自由状态下的所有频率响应函数利用有限元模型得到数值解,然后利用响应耦合技术求解子结构A和子结构B结合面处的耦合频率响应函数矩阵H2;
4)由结合面处耦合频率响应函数矩阵H2计算得到第一测点处与转动自由度有关的频率响应函数,第二测点处与转动自由度有关的频率响应函数,第一测点与第二测点之间与转动自由度有关的频率响应函数。
2.如权利要求1所述的复杂机械结构的转动自由度频率响应函数计算方法,其特征在于,所述的转动自由度频率响应函数为待分析的复杂机械结构上任意两个激励点和响应点之间与转动自由度有关的频率响应函数;根据激励点和响应点,选定第一测点和第二测点。
3.如权利要求1所述的复杂机械结构的转动自由度频率响应函数计算方法,其特征在于,所述的子结构A和子结构B通过阻尼、转动刚度和平动刚度结合。
4.如权利要求1所述的复杂机械结构的转动自由度频率响应函数计算方法,其特征在于,若转动自由度频率响应函数的激励点和响应点为不同测点,则在子结构A上选择出第一测点,在子结构A和子结构B的结合面处选择出第二测点;如果子结构A上选的第一测点是激励点,则子结构A和子结构B的结合面处选择的第二测点为响应点;如果子结构A上选的第一测点是响应点,则子结构A和子结构B的结合面处选择的第二测点为激励点;
若转动自由度频率响应函数的激励点和响应点为同一测点,则该测点即为子结构A上的第一测点,仍在子结构A和子结构B的结合面处选择出第二测点。
5.如权利要求1所述的复杂机械结构的转动自由度频率响应函数计算方法,其特征在于,步骤2)利用锤击激励法测量三个平动自由度的频率响应函数时,激励点采用激振力锤来进行锤击,响应点利用加速度传感器来检测加速度振动响应信号,由信号采集***计算分析。
6.如权利要求1所述的复杂机械结构的转动自由度频率响应函数计算方法,其特征在于,求解子结构A和子结构B结合面处的频率响应函数矩阵H2为:
假设测点在一个平面内的运动由平动和转动自由度组成,输入力F是由力f和力矩M组成的向量,输出响应X由平动位移x和转动位移θ组成,输入力与输出响应的关系为
(1)式中, 其中hij,ff为测点j的平动自由度与测点i的平动自由度之间的频率响应函数,hij,fM为测点j的转动自由度与测点i的平动自由度之间的频率响应函数,hij,Mf为测点j的平动自由度与测点i的转动自由度之间的频率响应函数,hij,MM为测点j的转动自由度与测点i的转动自由度之间的频率响应函数;
在复杂机械结构中第一测点处施加外力F1,只考虑第一测点和第二测点处的响应X1和X2,得到复杂机械结构的频率响应函数矩阵G11和G21如下:
(2)式中,HA,11为子结构A中第一测点的原点频率响应函数矩阵,HA,12为子结构A中第一测点和第二测点之间的频率响应函数矩阵,HA,21为子结构A中第二测点和第一测点之间的频率响应函数矩阵,HA,22为子结构A中第二测点的原点频率响应函数矩阵,H2为结合面处耦合频率响应函数矩阵;
HA,11,HA,12,HA,21,HA,22和H2的矩阵形式表达式为:
G11和G21的矩阵表达式为:
仅在第二测点处施加外力F2,可得频率响应函数矩阵G12和G22:
上式中,G12和G22的矩阵表达式为:
将G11,G21,G12和G22用频率响应函数矩阵形式表示如(4)式:
分别取频率响应函数矩阵G11,G12和G22中的第一个元素,得方程组式(5):
式(5)为一个含有4个未知数h2,ff,h2,fM,h2,Mf和h2,MM的方程组,子结构A在自由状态下的所有频率响应函数HA,11,HA,12,HA,21和HA,22利用有限元模型得到数值解;g11,ff,g12,ff和g22,ff为复杂机械结构上平动自由度的频率响应函数,由步骤2)得到;利用式(5)解出结合面处耦合频率响应函数矩阵H2。
7.如权利要求6所述的复杂机械结构的转动自由度频率响应函数计算方法,其特征在于,将结合面处耦合频率响应函数矩阵H2代入式(4)中,计算得到第一测点处与转动自由度有关的频率响应函数g11,fM、g11,Mf和g11,MM,其中:g11,fM为第一测点处转动自由度与平动自由度之间的频率响应函数,g11,Mf为第一测点处平动自由度与转动自由度之间的频率响应函数,g11,MM为第一测点处转动自由度的原点频率响应函数;
第二测点处与转动自由度有关的频率响应函数g22,fM、g22,Mf和g22,MM,其中:g22,fM为第二测点处转动自由度与平动自由度之间的频率响应函数,g22,Mf为第二测点处平动自由度与转动自由度之间的频率响应函数,g22,MM为第二测点处转动自由度的原点频率响应函数;
第一测点和第二测点之间与转动自由度有关的频率响应函数g21,fM、g21,Mf、g21,MM、g12,fM、g12,Mf和g12,MM,其中:g21,fM为第一测点处转动自由度与第二测点处平动自由度之间的频率响应函数,g21,Mf为第一测点处平动自由度与第二测点处转动自由度之间的频率响应函数,g21,MM为第一测点处转动自由度与第二测点处转动自由度之间的频率响应函数,g12,fM为第二测点处转动自由度与第一测点处平动自由度之间的频率响应函数,g12,Mf为第二测点处平动自由度与第一测点处转动自由度之间的频率响应函数,g12,MM为第二测点处转动自由度与第一测点处转动自由度之间的频率响应函数。
8.如权利要求7所述的复杂机械结构的转动自由度频率响应函数计算方法,其特征在于,所述的第一测点处转动自由度的原点频率响应函数g11,MM、第二测点处转动自由度的原点频率响应函数g22,MM,第一测点与第二测点之间转动自由度频率响应函数g21,MM,分别表示如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210371646.1A CN102880803B (zh) | 2012-09-29 | 2012-09-29 | 一种复杂机械结构的转动自由度频率响应函数计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210371646.1A CN102880803B (zh) | 2012-09-29 | 2012-09-29 | 一种复杂机械结构的转动自由度频率响应函数计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102880803A true CN102880803A (zh) | 2013-01-16 |
CN102880803B CN102880803B (zh) | 2014-04-23 |
Family
ID=47482125
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210371646.1A Active CN102880803B (zh) | 2012-09-29 | 2012-09-29 | 一种复杂机械结构的转动自由度频率响应函数计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102880803B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104198143A (zh) * | 2014-09-11 | 2014-12-10 | 江南大学 | 一种测定极脆部件频率响应函数的方法 |
CN104236834A (zh) * | 2014-09-24 | 2014-12-24 | 江南大学 | 一种在线测定极脆部件频率响应函数的高精度方法 |
CN104236833A (zh) * | 2014-09-24 | 2014-12-24 | 江南大学 | 一种测定极脆部件频率响应函数的高精度方法 |
CN104236829A (zh) * | 2014-09-11 | 2014-12-24 | 江南大学 | 一种测定极脆部件非耦合界面频率响应函数的方法 |
CN106197564A (zh) * | 2016-09-18 | 2016-12-07 | 湖南工学院 | 一种适用于识别连接结合面参数的方法与*** |
CN106503318A (zh) * | 2016-10-12 | 2017-03-15 | 华中科技大学 | 一种机床主轴工作状态下的刀具端频响函数辨识方法 |
CN107000151A (zh) * | 2014-12-08 | 2017-08-01 | 山特维克知识产权股份有限公司 | 测量装置和选择去屑加工刀具的运行参数的方法 |
CN108170643A (zh) * | 2017-12-26 | 2018-06-15 | 上海交通大学 | 基于原位测量频响函数的免拆分子结构频响函数识别方法 |
CN108489736A (zh) * | 2018-03-28 | 2018-09-04 | 中汽研(天津)汽车工程研究院有限公司 | 一种车身联接结构件耦联损耗因子测试***及方法 |
CN109557186A (zh) * | 2018-12-24 | 2019-04-02 | 暨南大学 | 一种机械装配动态质量检测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101458205A (zh) * | 2008-08-29 | 2009-06-17 | 华中科技大学 | 一种机床固定结合部动力学参数的识别方法 |
CN101504323A (zh) * | 2009-03-19 | 2009-08-12 | 上海交通大学 | 转角动柔度的间接测量方法 |
CN102135466A (zh) * | 2010-11-18 | 2011-07-27 | 西安航空动力股份有限公司 | 一种振动模态测试方法 |
US8113050B2 (en) * | 2006-01-25 | 2012-02-14 | The Regents Of The University Of California | Robust six degree-of-freedom micromachined gyroscope with anti-phase drive scheme and method of operation of the same |
-
2012
- 2012-09-29 CN CN201210371646.1A patent/CN102880803B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8113050B2 (en) * | 2006-01-25 | 2012-02-14 | The Regents Of The University Of California | Robust six degree-of-freedom micromachined gyroscope with anti-phase drive scheme and method of operation of the same |
CN101458205A (zh) * | 2008-08-29 | 2009-06-17 | 华中科技大学 | 一种机床固定结合部动力学参数的识别方法 |
CN101504323A (zh) * | 2009-03-19 | 2009-08-12 | 上海交通大学 | 转角动柔度的间接测量方法 |
CN102135466A (zh) * | 2010-11-18 | 2011-07-27 | 西安航空动力股份有限公司 | 一种振动模态测试方法 |
Non-Patent Citations (3)
Title |
---|
《机械工程学报》 20120502 曹宏瑞 等 机床-主轴耦合***动力学建模与模型修正 第48卷, 第3期 * |
宋攀 等: "实验模态综合法若干问题的研究", 《振动与冲击》 * |
曹宏瑞 等: "机床—主轴耦合***动力学建模与模型修正", 《机械工程学报》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104198143A (zh) * | 2014-09-11 | 2014-12-10 | 江南大学 | 一种测定极脆部件频率响应函数的方法 |
CN104236829A (zh) * | 2014-09-11 | 2014-12-24 | 江南大学 | 一种测定极脆部件非耦合界面频率响应函数的方法 |
CN104236834A (zh) * | 2014-09-24 | 2014-12-24 | 江南大学 | 一种在线测定极脆部件频率响应函数的高精度方法 |
CN104236833A (zh) * | 2014-09-24 | 2014-12-24 | 江南大学 | 一种测定极脆部件频率响应函数的高精度方法 |
CN107000151A (zh) * | 2014-12-08 | 2017-08-01 | 山特维克知识产权股份有限公司 | 测量装置和选择去屑加工刀具的运行参数的方法 |
CN107000151B (zh) * | 2014-12-08 | 2019-07-05 | 山特维克知识产权股份有限公司 | 测量装置和选择去屑加工刀具的运行参数的方法 |
CN106197564A (zh) * | 2016-09-18 | 2016-12-07 | 湖南工学院 | 一种适用于识别连接结合面参数的方法与*** |
CN106197564B (zh) * | 2016-09-18 | 2019-05-14 | 湖南工学院 | 一种适用于识别连接结合面参数的方法与*** |
CN106503318A (zh) * | 2016-10-12 | 2017-03-15 | 华中科技大学 | 一种机床主轴工作状态下的刀具端频响函数辨识方法 |
CN106503318B (zh) * | 2016-10-12 | 2019-05-21 | 华中科技大学 | 一种机床主轴工作状态下的刀具端频响函数辨识方法 |
CN108170643A (zh) * | 2017-12-26 | 2018-06-15 | 上海交通大学 | 基于原位测量频响函数的免拆分子结构频响函数识别方法 |
CN108170643B (zh) * | 2017-12-26 | 2021-02-09 | 上海交通大学 | 一种免拆分的子结构频响函数识别方法 |
CN108489736A (zh) * | 2018-03-28 | 2018-09-04 | 中汽研(天津)汽车工程研究院有限公司 | 一种车身联接结构件耦联损耗因子测试***及方法 |
CN109557186A (zh) * | 2018-12-24 | 2019-04-02 | 暨南大学 | 一种机械装配动态质量检测方法 |
CN109557186B (zh) * | 2018-12-24 | 2021-05-28 | 暨南大学 | 一种机械装配动态质量检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102880803B (zh) | 2014-04-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102880803B (zh) | 一种复杂机械结构的转动自由度频率响应函数计算方法 | |
Li et al. | Identification of damage locations for plate-like structures using damage sensitive indices: strain modal approach | |
Ramos | Damage identification on masonry structures based on vibration signatures | |
Cho et al. | Forced vibration analysis of arbitrarily constrained rectangular plates and stiffened panels using the assumed mode method | |
CN104457681A (zh) | 一种基于应变模态的梁结构动挠度监测方法 | |
Ozer et al. | Hybrid motion sensing and experimental modal analysis using collocated smartphone camera and accelerometers | |
CN102621225A (zh) | 一种路面及桥面铺装材料阻尼特性参数的测试方法 | |
Ramos et al. | Operational modal analysis of historical constructions using commercial wireless platforms | |
Le et al. | Computationally efficient stochastic approach for the fragility analysis of vertical structures subjected to thunderstorm downburst winds | |
CN104133950A (zh) | 一种悬臂梁运行模态分析实验方法及装置 | |
Sun et al. | Modal identification of vibrating structures using singular value decomposition and nonlinear iteration based on high-speed digital image correlation | |
CN103217212A (zh) | 一种加工中心的刀具刀尖点位移导纳的软测量方法及*** | |
CN103900826A (zh) | 实时监测汽车底盘结构疲劳损伤的方法 | |
CN112507585B (zh) | 基于emd和模型缩聚的动力响应重构方法及***、存储介质 | |
CN106202735A (zh) | 一种局部非线性地基土‑结构相互作用子结构试验方法 | |
Gupta | Inverse methods for load identification augmented by optimal sensor placement and model order reduction | |
Zhu et al. | A modal test method using sound pressure transducers based on vibro-acoustic reciprocity | |
CN104280457A (zh) | 一种自升式平台损伤识别方法及装置 | |
Wang et al. | Damage identification for beam structures using the laplace transform-based spectral element method and strain statistical moment | |
Pan | System identification of constructed civil engineering structures and uncertainty | |
Helfrick | An investigation of 3D digital image correlation for structural health monitoring and vibration measurement | |
CN108304653B (zh) | 晶体培植工艺环境振动影响定量快速评估技术 | |
Gibanica | Experimental-analytical dynamic substructuring: a state-space approach | |
CN100430731C (zh) | 微型惯性传感器件的芯核建模方法及芯核库 | |
CN104897276A (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 |