[go: up one dir, main page]

CN109726521B - Method for eliminating suspension support influence in frequency response function aiming at free mode test - Google Patents

Method for eliminating suspension support influence in frequency response function aiming at free mode test Download PDF

Info

Publication number
CN109726521B
CN109726521B CN201910119087.7A CN201910119087A CN109726521B CN 109726521 B CN109726521 B CN 109726521B CN 201910119087 A CN201910119087 A CN 201910119087A CN 109726521 B CN109726521 B CN 109726521B
Authority
CN
China
Prior art keywords
frequency response
response function
point
suspension
displacement
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.)
Active
Application number
CN201910119087.7A
Other languages
Chinese (zh)
Other versions
CN109726521A (en
Inventor
任军
张强豪
何文浩
吴正虎
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hubei University of Technology
Original Assignee
Hubei University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hubei University of Technology filed Critical Hubei University of Technology
Priority to CN201910119087.7A priority Critical patent/CN109726521B/en
Publication of CN109726521A publication Critical patent/CN109726521A/en
Application granted granted Critical
Publication of CN109726521B publication Critical patent/CN109726521B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Vehicle Body Suspensions (AREA)

Abstract

The invention discloses a method for eliminating suspension support influence in a frequency response function aiming at free mode test, wherein in the free mode test, a tested structure is usually suspended and supported by adopting a soft rope or a soft spring and the like, and the free boundary condition is approximately simulated; however, the suspension support can introduce additional rigidity and damping to the test structure, and further bring a certain influence on the measured frequency response function; the invention relates to a method for eliminating suspension support influence in a frequency response function aiming at free modal test, which can eliminate the influence of suspension support rigidity and damping on the measured frequency response function and accurately acquire the frequency response function of a free structure by measuring the frequency response function related to a suspension point under the condition of giving the support rigidity and the support damping.

Description

针对自由模态测试消除频响函数中悬挂支承影响的方法Method for Eliminating Suspension Support Effects in Frequency Response Functions for Free Mode Testing

技术领域technical field

本发明属于模态测试技术领域,涉及一种频响函数数据修正的方法,具体涉及一种针对自由模态测试消除频响函数中悬挂支承影响的方法,此处悬挂支承影响主要指悬挂支承刚度和阻尼对频响函数引入的附加影响。The invention belongs to the technical field of modal testing, and relates to a method for correcting frequency response function data, in particular to a method for eliminating the influence of suspension supports in frequency response functions for free modal tests, where the influence of suspension supports mainly refers to the stiffness of suspension supports and the additional effect of damping on the frequency response function.

背景技术Background technique

自由模态测试是比较常见的一类模态试验,以获取自由结构的动态特性为目的。这类测试中,被测结构应当处于类似于“悬浮”的自由状态。然而实践中,这种状态很难实现,通常会采用柔软的绳索或软弹簧等对被测结构进行悬挂支承,近似模拟自由边界条件。显然,悬挂支承会给测试结构引入附加刚度和阻尼,进而给测试结果带来一定程度的影响,最直接影响的就是测量的频响函数。Free modal test is a relatively common type of modal test, which aims to obtain the dynamic characteristics of free structures. In this type of test, the structure under test should be in a free state similar to "suspension". However, in practice, this state is difficult to achieve. Usually, soft ropes or soft springs are used to suspend and support the structure under test to approximate free boundary conditions. Obviously, the suspension support will introduce additional stiffness and damping to the test structure, which will affect the test results to a certain extent, and the most direct impact is the measured frequency response function.

发明内容Contents of the invention

为了解决上述技术问题,本发明公开了一种针对自由模态测试消除频响函数中悬挂支承影响的方法。In order to solve the above technical problems, the present invention discloses a method for eliminating the influence of suspension support in the frequency response function for free mode testing.

本发明所采用的技术方案是:一种针对自由模态测试消除频响函数中悬挂支承影响的方法,其特征在于:针对单点悬挂,假定结构体C的自身刚度足够小而不能忽略悬挂支承的影响,结构体C悬挂于s点处,悬挂支承效果等效为ks刚度的弹簧和大小为cs的阻尼,p点为激励点,l点为测量点,其中“p点”和“l点”为结构体上任意两点;为了消除悬挂支承的影响,在s点施加“虚拟支承”,“虚拟支承”中的“负刚度”k′s=-ks和“负阻尼”c′s=-cs抵消了悬挂支承的影响,使结构体等效为自由状态;The technical solution adopted in the present invention is: a method for eliminating the influence of the suspension support in the frequency response function for the free mode test, which is characterized in that: for the single-point suspension, it is assumed that the self-stiffness of the structure C is small enough to ignore the suspension support The structure C is suspended at point s, the suspension support effect is equivalent to a spring with a stiffness of k s and a damper with a size of c s , point p is the excitation point, and point l is the measurement point, where "point p" and " Point l" is any two points on the structure; in order to eliminate the influence of suspension support, a "virtual support" is applied at point s, and the "negative stiffness" k′ s = -k s and "negative damping" c in the "virtual support" ′ s = -c s offsets the influence of the suspension support, making the structure equivalent to a free state;

为便于分析,将s点拆分为两个点s点和s′点,在“虚拟支承”中,用s′点来表述s点;For the convenience of analysis, point s is split into two points, s and s′, and in the “virtual support”, point s is expressed by point s′;

视悬挂的结构体C为子结构Ⅰ,“虚拟支承”为子结构Ⅱ,则子结构Ⅰ和子结构Ⅱ结合即得到自由结构体;Assuming that the suspended structure C is substructure I, and the "virtual support" is substructure II, the combination of substructure I and substructure II results in a free structure;

对子结构Ⅰ和Ⅱ分别进行分析;Analyze substructures I and II separately;

对于子结构Ⅰ,建立频域下的输出位移响应和输入力的关系,For substructure Ⅰ, the relationship between the output displacement response and the input force in the frequency domain is established,

Figure BDA0001971219910000011
Figure BDA0001971219910000011

其中,Xl和Xs分别表示子结构Ⅰ在l点和s点处的位移响应,Fp和Rs分别表示子结构Ⅰ在p点和s点所受的作用力;

Figure BDA0001971219910000021
表示悬挂状态下l、s点之间的位移频响函数;
Figure BDA0001971219910000022
表示悬挂状态下s、p点之间的位移频响函数;
Figure BDA0001971219910000023
表示悬挂状态下s点处的原点位移频响函数;Among them, X l and X s represent the displacement response of substructure I at points l and s, respectively, and F p and R s represent the forces on substructure I at points p and s, respectively;
Figure BDA0001971219910000021
Indicates the displacement frequency response function between points l and s in the suspended state;
Figure BDA0001971219910000022
Indicates the displacement frequency response function between points s and p in the suspended state;
Figure BDA0001971219910000023
Indicates the origin displacement frequency response function at point s in the suspended state;

对于子结构Ⅱ,同样建立输出位移响应和输入力的关系,For substructure II, the relationship between the output displacement response and the input force is also established,

X′s=αs′s′Rs′ (2)X' s = α s's' R s' (2)

式中,In the formula,

Figure BDA0001971219910000024
Figure BDA0001971219910000024

其中,X′s′和Rs′分别表示子结构Ⅱ在s点处的位移响应和所受的作用力;αs′s′表示子结构Ⅱ在s点处的原点位移频响函数;jω中j是复数单位,ω为频响函数自变量;Among them, X′ s′ and R s′ represent the displacement response and the acting force of substructure II at point s respectively; α s′s′ represents the frequency response function of the origin displacement of substructure II at point s; jω where j is a complex unit, and ω is an independent variable of the frequency response function;

连接点s或s′处力和位移的约束条件分别为:The constraint conditions of force and displacement at connection point s or s′ are respectively:

Rs+Rs′=0 (4)R s +R s′ = 0 (4)

Xs=Xs′ (5)X s =X s′ (5)

联立式(1)~(5),得到Simultaneous formulas (1)~(5), get

Figure BDA0001971219910000025
Figure BDA0001971219910000025

式中,αlp表示自由状态下l、p点位移频响函数;

Figure BDA0001971219910000026
表示悬挂状态下l、p点位移频响函数;为了避免混淆,命名自由状态下的频响函数为“准确频响函数”,悬挂状态下的频响为“实测频响函数”,即αlp为准确频响函数,
Figure BDA0001971219910000027
均为实测频响函数;In the formula, α lp represents the displacement frequency response function of points l and p in the free state;
Figure BDA0001971219910000026
Indicates the displacement frequency response function of points l and p in the suspension state; in order to avoid confusion, the frequency response function in the free state is named "accurate frequency response function", and the frequency response in the suspension state is "measured frequency response function", that is, α lp is the exact frequency response function,
Figure BDA0001971219910000027
Both are measured frequency response functions;

根据式(6)可知,准确频响函数αlp和实测频响函数

Figure BDA0001971219910000028
有所差异,而这种差异正是由于支承的刚度和阻尼导致;而且在悬挂刚度ks和阻尼cs给定的情况下,通过测量与悬挂点s相关的位移频响函数可计算得到准确频响函数αlp,即达到消除位移频响函数
Figure BDA0001971219910000031
中支承影响的目的;According to formula (6), it can be seen that the exact frequency response function α lp and the measured frequency response function
Figure BDA0001971219910000028
There is a difference, and this difference is caused by the stiffness and damping of the support; and when the suspension stiffness k s and damping c s are given, the exact frequency response function can be calculated by measuring the displacement frequency response function related to the suspension point s Frequency response function α lp , that is, to achieve the elimination of displacement frequency response function
Figure BDA0001971219910000031
The purpose of supporting influence in the medium;

以上所分析均针对的位移频响函数,而实际测量中大多采用加速度传感器,所测量频响函数为加速度频响函数。对于加速度频响的情况,只需将式(6)左右乘以-ω2,即得到加速度频响和准确加速度频响之间的关系式(7)The above analyzes are all aimed at the displacement frequency response function, but in the actual measurement, acceleration sensors are mostly used, and the measured frequency response function is the acceleration frequency response function. For the case of acceleration frequency response, just multiply the left and right sides of formula (6) by -ω 2 to get the relationship between acceleration frequency response and accurate acceleration frequency response (7)

Figure BDA0001971219910000032
Figure BDA0001971219910000032

式中A均代表加速度频响函数,Alp表示自由结构体l、p点间的加速度频响函数;

Figure BDA0001971219910000033
表示悬挂状态下l、p点之间的加速度频响函数;
Figure BDA0001971219910000034
表示悬挂状态下s、p点之间的加速度频响函数;
Figure BDA0001971219910000035
表示悬挂状态下l、s点之间的加速度频响函数;
Figure BDA0001971219910000036
表示悬挂状态下s点加速度原点频响函数;只要支承刚度ks和支承阻尼cs大小已知,就能通过测量与悬挂点相关的频响函数来获取自由结构的频响函数Alp。In the formula, A represents the acceleration frequency response function, and A lp represents the acceleration frequency response function between points l and p of the free structure;
Figure BDA0001971219910000033
Indicates the acceleration frequency response function between points l and p in the suspended state;
Figure BDA0001971219910000034
Indicates the acceleration frequency response function between points s and p in the suspended state;
Figure BDA0001971219910000035
Indicates the acceleration frequency response function between points l and s in the suspended state;
Figure BDA0001971219910000036
Represents the frequency response function of the acceleration origin at point s in the suspension state; as long as the support stiffness k s and support damping c s are known, the frequency response function Alp of the free structure can be obtained by measuring the frequency response function related to the suspension point.

针对多点悬挂,假定结构体C的自身刚度足够小而不能忽略悬挂支承的影响,结构体C处于多点悬挂状态,悬挂点分别为1、2、…、n;每个悬挂点支承刚度为ki,阻尼大小为ci,i=1、2、…、n;p点为激励点,l点为测量点,其中“p点”和“l点”为结构体上任意两点;为了消除悬挂支承的影响,依次在每个悬挂点处施加“虚拟支承”,“虚拟支承”中的“负刚度”k′i=-ki和“负阻尼”c′i=-ci分别抵消各悬挂支承的影响,从而得到自由状态的结构体;For multi-point suspension, it is assumed that the self-stiffness of structure C is small enough to ignore the influence of suspension support, structure C is in a multi-point suspension state, and the suspension points are 1, 2, ..., n respectively; the support stiffness of each suspension point is k i , the damping size is c i , i=1, 2, ..., n; point p is the excitation point, point l is the measurement point, where "point p" and "point l" are any two points on the structure; for To eliminate the influence of the suspension support, apply a "virtual support" at each suspension point in turn, and the "negative stiffness" k′ i = -ki and "negative damping" c′ i = -ci in the "virtual support" are offset respectively The influence of each suspension support, so as to obtain the structure in a free state;

为便于分析,将悬挂点si点拆分为两个点si点和s′i点,在“虚拟支承”中,用s′i点来表述si点;For the convenience of analysis, the suspension point s i is split into two points s i and s′ i , and in the “virtual support”, s i is used to express s i ;

将原悬挂的结构体C视为子结构Ⅰ,将虚拟支承视为子结构Ⅱ、Ⅲ、……;对于子结构Ⅰ,建立频域下的输出位移响应和输入力的关系,The original suspended structure C is regarded as substructure I, and the virtual support is regarded as substructure II, III, ...; for substructure I, the relationship between the output displacement response and the input force in the frequency domain is established,

Figure BDA0001971219910000037
Figure BDA0001971219910000037

Figure BDA0001971219910000038
Figure BDA0001971219910000038

其中,k=1、2、…、n,Xl和Xk分别表示子结构Ⅰ在l点和k点处的位移Among them, k=1, 2, ..., n, X l and X k represent the displacement of substructure I at point l and point k respectively

响应,Fp

Figure BDA0001971219910000047
和Ri分别表示子结构Ⅰ在p点、Si点和i点所受的作用力;
Figure BDA0001971219910000041
表示悬挂状态下l、p点之间的位移频响函数;
Figure BDA0001971219910000042
表示悬挂状态下l、Si点之间的位移频响函数;
Figure BDA0001971219910000043
表示悬挂状态下k、p点之间的位移频响函数;
Figure BDA0001971219910000044
表示悬挂状态下k、i点之间的位移频响函数;Response, F p ,
Figure BDA0001971219910000047
and R i denote the acting force on substructure I at point p, S i and i, respectively;
Figure BDA0001971219910000041
Indicates the displacement frequency response function between points l and p in the suspended state;
Figure BDA0001971219910000042
Indicates the displacement frequency response function between points l and S i in the suspended state;
Figure BDA0001971219910000043
Indicates the displacement frequency response function between points k and p in the suspension state;
Figure BDA0001971219910000044
Indicates the displacement frequency response function between points k and i in the suspended state;

对于子结构Ⅱ、Ⅲ、……,同样建立输出位移响应和输入力的关系,For substructures II, III, ..., the relationship between the output displacement response and the input force is also established,

X′k=α′kkR′k (10)X′ k =α′ kk R′ k (10)

式中,In the formula,

Figure BDA0001971219910000045
Figure BDA0001971219910000045

其中,X′k和Rk分别表示子结构Ⅱ、Ⅲ、……在k点处的位移响应和所受的作用力;α′kk表示子子结构Ⅱ、Ⅲ、……在k点处的原点位移频响函数;jω中j是复数单位,ω为频响函数自变量;Among them, X′ k and R k represent the displacement responses and acting forces of substructures II, III, ... at point k, respectively; Origin displacement frequency response function; j in jω is a complex unit, and ω is an independent variable of frequency response function;

连接点s或s′处力和位移的约束条件分别为:The constraint conditions of force and displacement at connection point s or s′ are respectively:

Rk+R′k=0 (12)R k + R′ k = 0 (12)

Xk=X′k (13)X k = X′ k (13)

联立式(8)~(13),得到Simultaneous formula (8)~(13), get

α(s)·F=X (14)α (s) F=X (14)

式中,In the formula,

Figure BDA0001971219910000046
Figure BDA0001971219910000046

求解式(14),得Solving formula (14), we get

Figure BDA0001971219910000051
Figure BDA0001971219910000051

式中,αlp表示自由结构体l、p点间的位移频响函数;

Figure BDA0001971219910000052
表示悬挂结构体的位移频响函数;α′11,α′22……表示添加的虚拟支承分别在1、2……点处的原点位移频响函数;式(15)表明,对位移频响函数
Figure BDA0001971219910000053
修正需要额外测量与悬挂点1,2…n相关的多个位移频响函数,所需测量的位移频响函数数量随着悬挂点的增多而增加;In the formula, α lp represents the displacement frequency response function between points l and p of the free structure;
Figure BDA0001971219910000052
Indicates the displacement frequency response function of the suspended structure; α′ 11 , α′ 22 ... represent the origin displacement frequency response function of the added virtual support at points 1, 2... respectively; Equation (15) shows that the displacement frequency response function
Figure BDA0001971219910000053
The correction requires additional measurement of multiple displacement frequency response functions related to suspension points 1, 2...n, and the number of displacement frequency response functions required to be measured increases with the number of suspension points;

从式(15)推算出,对于N点悬挂的情况,修正

Figure BDA0001971219910000054
需要额外测量(N2+5N)/2组位移频响函数,分别为原点位移频响函数
Figure BDA0001971219910000055
和跨点位移频响函数
Figure BDA0001971219910000056
这些位移频响函数均与悬挂点相关;以上结论是针对激励点p、响应点l和悬挂点1,2…n均不重合的情况,当它们其中某些点重合时,式(15)会适当简化,修正工作中所需测量的位移频响函数也会适当减少;It is deduced from formula (15), for the case of N-point suspension, the correction
Figure BDA0001971219910000054
Additional measurements (N 2 +5N)/2 sets of displacement frequency response functions are required, which are the origin displacement frequency response functions
Figure BDA0001971219910000055
and across point displacement frequency response function
Figure BDA0001971219910000056
These displacement frequency response functions are all related to the suspension point; the above conclusion is for the case that the excitation point p, the response point l and the suspension point 1, 2...n do not overlap. When some of them overlap, the formula (15) will be Appropriate simplification, the displacement frequency response function that needs to be measured in the correction work will also be appropriately reduced;

对于加速度频响的情况,只需将式(15)两端分别乘以-ω2即得到加速度频响函数A的表达形式,For the case of acceleration frequency response, it is only necessary to multiply both ends of formula (15) by -ω 2 to obtain the expression form of the acceleration frequency response function A,

Figure BDA0001971219910000057
Figure BDA0001971219910000057

式中A均代表加速度频响函数,Alp表示自由结构体l、p点间的加速度频响函数;

Figure BDA0001971219910000058
表示悬挂状态下l、p点之间的加速度频响函数;
Figure BDA0001971219910000059
表示悬挂状态下n点加速度原点频响函数;A′11,A′22……A′nn分别表示添加的虚拟支承在1、2……n点处的加速度频响函数。In the formula, A represents the acceleration frequency response function, and A lp represents the acceleration frequency response function between points l and p of the free structure;
Figure BDA0001971219910000058
Indicates the acceleration frequency response function between points l and p in the suspended state;
Figure BDA0001971219910000059
Indicates the frequency response function of the acceleration origin of point n in the suspension state; A′ 11 , A′ 22 ... A' nn represent the acceleration frequency response function of the added virtual support at points 1, 2...n, respectively.

本发明涉及一种针对自由模态测试消除频响函数中悬挂支承影响的方法,该方法在给定支承刚度和支承阻尼的情况下,通过测量与悬挂点相关的频响函数,即可消除悬挂支承刚度和阻尼对所测量频响函数的影响,准确获取自由结构的频响函数。The invention relates to a method for eliminating the influence of the suspension support in the frequency response function for the free mode test. Under the condition of given support stiffness and support damping, the method can eliminate the suspension support by measuring the frequency response function related to the suspension point. The influence of support stiffness and damping on the measured frequency response function, to accurately obtain the frequency response function of the free structure.

附图说明Description of drawings

图1为本发明背景技术中自由模态测试悬挂支承示意图;Fig. 1 is the free mode test suspension support schematic diagram in the background technology of the present invention;

图2为本发明实施例的结构体自由模态测试示意图(单点悬挂);Fig. 2 is the structure free modal test schematic diagram (single-point suspension) of the embodiment of the present invention;

图3为本发明实施例的悬挂结构体转换为自由结构体示意图;Fig. 3 is a schematic diagram of the conversion of a suspended structure into a free structure according to an embodiment of the present invention;

图4为本发明实施例的结构体自由模态测试示意图(多点悬挂);Fig. 4 is the schematic diagram of structure free modal test (multi-point suspension) of the embodiment of the present invention;

图5为本发明实施例的悬挂结构体转变为自由结构体示意图;Fig. 5 is a schematic diagram of the transformation of a suspended structure into a free structure according to an embodiment of the present invention;

图6为本发明实施例的自由梁模态测试示意图;Fig. 6 is the free beam modal test schematic diagram of the embodiment of the present invention;

图7为本发明实施例的测量的频响函数

Figure BDA0001971219910000061
与准确频响函数A24对比示意图;Fig. 7 is the frequency response function of the measurement of the embodiment of the present invention
Figure BDA0001971219910000061
Schematic diagram of comparison with the exact frequency response function A 24 ;

图8为本发明实施例的修正后的频响函数

Figure BDA0001971219910000062
与准确频响函数A24对比示意图;Fig. 8 is the revised frequency response function of the embodiment of the present invention
Figure BDA0001971219910000062
Schematic diagram of comparison with the exact frequency response function A 24 ;

图9为本发明实施例的修正后的频响函数

Figure BDA0001971219910000063
(未消除支承阻尼)与准确频响函数A24对比示意图。Fig. 9 is the revised frequency response function of the embodiment of the present invention
Figure BDA0001971219910000063
(without eliminating the support damping) and the accurate frequency response function A 24 comparison diagram.

具体实施方式detailed description

为了便于本领域普通技术人员理解和实施本发明,下面结合附图及实施例对本发明作进一步的详细描述,应当理解,此处所描述的实施示例仅用于说明和解释本发明,并不用于限定本发明。In order to facilitate those of ordinary skill in the art to understand and implement the present invention, the present invention will be described in further detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the implementation examples described here are only used to illustrate and explain the present invention, and are not intended to limit this invention.

请见图2,本发明提供的一种针对自由模态测试消除频响函数中悬挂支承影响的方法,针对单点悬挂,结构体C悬挂于s点处,悬挂支承效果等效为ks刚度的弹簧和大小为cs的阻尼。p点为激励点,l点为测量点。假定结构体C的自身刚度足够小而不能忽略悬挂支承的影响。为了消除悬挂支承的影响,可在s点施加“虚拟支承”,如图3中左侧的虚线框内部分。“虚拟支承”中的“负刚度”(k′s=-ks)和“负阻尼”(c′s=-cs)抵消了悬挂支承的影响,使结构体等效为图3右侧的自由状态。图3左图中s点和s′点实际代表同一点,为便于分析,拆分为两个点。Please see Fig. 2, a method for eliminating the influence of suspension support in the frequency response function provided by the present invention for free mode testing, for single-point suspension, the structure C is suspended at point s, and the suspension support effect is equivalent to k s stiffness A spring of size c s and a damper of size c s. Point p is the excitation point, and point l is the measurement point. Assume that the self-stiffness of structure C is small enough to ignore the influence of suspension supports. In order to eliminate the influence of suspension support, a "virtual support" can be applied at point s, as shown in the dotted line box on the left in Figure 3. The “negative stiffness” (k′ s =-k s ) and “negative damping” (c′ s =-c s ) in the “virtual support” offset the influence of the suspension support, making the structure equivalent to the right side of Figure 3 free state. Point s and point s′ in the left figure of Figure 3 actually represent the same point, and are split into two points for the convenience of analysis.

本质上,对原悬挂结构体C施加“虚拟支承”属于结构修改的范畴,因此,可以用结构动力学修改的方法来进行建模。若视悬挂的结构体C为子结构Ⅰ,虚拟支承为子结构Ⅱ,则子结构Ⅰ和子结构Ⅱ结合即可得到如图3右侧的自由结构体。对子结构Ⅰ和Ⅱ分别进行分析。对于子结构Ⅰ,可建立频域下的输出位移响应和输入力的关系,Essentially, imposing "virtual support" on the original suspension structure C belongs to the category of structural modification, so the method of structural dynamic modification can be used for modeling. If the suspended structure C is regarded as substructure I, and the virtual support is regarded as substructure II, then the free structure on the right side of Figure 3 can be obtained by combining substructure I and substructure II. Substructures I and II were analyzed separately. For substructure I, the relationship between the output displacement response and the input force in the frequency domain can be established,

Figure BDA0001971219910000071
Figure BDA0001971219910000071

其中,Xl和Xs分别表示子结构Ⅰ在l点和s点处的位移响应,Fp和Rs分别表示子结构Ⅰ在p点和s点所受的作用力;

Figure BDA0001971219910000072
表示悬挂状态下l、s点之间的位移频响函数;
Figure BDA0001971219910000073
表示悬挂状态下s、p点之间的位移频响函数;
Figure BDA0001971219910000074
表示悬挂状态下s点处的原点位移频响函数;Among them, X l and X s represent the displacement response of substructure I at points l and s, respectively, and F p and R s represent the forces on substructure I at points p and s, respectively;
Figure BDA0001971219910000072
Indicates the displacement frequency response function between points l and s in the suspended state;
Figure BDA0001971219910000073
Indicates the displacement frequency response function between points s and p in the suspended state;
Figure BDA0001971219910000074
Indicates the origin displacement frequency response function at point s in the suspended state;

对于子结构Ⅱ,同样建立输出位移响应和输入力的关系,For substructure II, the relationship between the output displacement response and the input force is also established,

X′s=αs′s′Rs′ (2)X' s = α s's' R s' (2)

式中,In the formula,

Figure BDA0001971219910000075
Figure BDA0001971219910000075

其中,X′s′和Rs′分别表示子结构Ⅱ在s点处的位移响应和所受的作用力;αs′s′表示子结构Ⅱ在s点处的原点位移频响函数;jω中j是复数单位,ω为频响函数自变量;Among them, X′ s′ and R s′ represent the displacement response and the acting force of substructure II at point s respectively; α s′s′ represents the frequency response function of the origin displacement of substructure II at point s; jω where j is a complex unit, and ω is an independent variable of the frequency response function;

连接点s或s′处力和位移的约束条件分别为The constraint conditions of force and displacement at connection point s or s′ are respectively

Rs+Rs′=0 (4)R s +R s′ = 0 (4)

Xs=Xs′ (5)X s =X s′ (5)

联立式(1)~(5),得到Simultaneous formulas (1)~(5), get

Figure BDA0001971219910000076
Figure BDA0001971219910000076

式中,αlp表示自由状态下l、p点位移频响函数;

Figure BDA0001971219910000081
表示结构体在悬挂状态下l、p点之间的位移频响函数;为了避免混淆,命名自由状态下的频响函数为“准确频响函数”,悬挂状态下的频响函数为“实测频响函数”,即αlp为准确频响函数,
Figure BDA0001971219910000082
均为实测频响函数;In the formula, α lp represents the displacement frequency response function of points l and p in the free state;
Figure BDA0001971219910000081
Indicates the displacement frequency response function between points l and p of the structure in the suspension state; in order to avoid confusion, the frequency response function in the free state is named "accurate frequency response function", and the frequency response function in the suspension state is "measured frequency response function". Response function", that is, α lp is the exact frequency response function,
Figure BDA0001971219910000082
Both are measured frequency response functions;

根据式(6)可知,准确频响函数αlp和实测频响函数

Figure BDA0001971219910000083
有所差异,而这种差异正是由于支承的刚度和阻尼导致;而且在悬挂刚度ks和阻尼cs给定的情况下,通过测量与悬挂点s相关的频响函数可计算得到准确频响函数αlp,即达到消除频响函数
Figure BDA0001971219910000084
中支承影响的目的;According to formula (6), it can be seen that the exact frequency response function α lp and the measured frequency response function
Figure BDA0001971219910000083
There is a difference, and this difference is caused by the stiffness and damping of the support; and when the suspension stiffness k s and damping c s are given, the exact frequency can be calculated by measuring the frequency response function related to the suspension point s Response function α lp , that is, to eliminate the frequency response function
Figure BDA0001971219910000084
The purpose of supporting influence in the medium;

以上所分析均针对的位移频响函数,而实际测量中大多采用加速度传感器,所测量频响函数为加速度频响函数。对于加速度频响的情况,只需将式(6)左右乘以-ω2,即得到加速度频响和准确加速度频响之间的关系式(7)The above analyzes are all aimed at the displacement frequency response function, but in the actual measurement, acceleration sensors are mostly used, and the measured frequency response function is the acceleration frequency response function. For the case of acceleration frequency response, just multiply the left and right sides of formula (6) by -ω 2 to get the relationship between acceleration frequency response and accurate acceleration frequency response (7)

Figure BDA0001971219910000085
Figure BDA0001971219910000085

式中A均代表加速度频响函数,Alp表示自由结构体l、p点间的加速度频响函数;

Figure BDA0001971219910000086
表示悬挂状态下l、p点之间的加速度频响函数;
Figure BDA0001971219910000087
表示悬挂状态下s、p点之间的加速度频响函数;
Figure BDA0001971219910000088
表示悬挂状态下l、s点之间的加速度频响函数;
Figure BDA0001971219910000089
表示悬挂状态下s点加速度原点频响函数;只要支承刚度ks和支承阻尼cs大小已知,就能通过测量与悬挂点相关的频响函数来获取自由结构的频响函数Alp。In the formula, A represents the acceleration frequency response function, and A lp represents the acceleration frequency response function between points l and p of the free structure;
Figure BDA0001971219910000086
Indicates the acceleration frequency response function between points l and p in the suspended state;
Figure BDA0001971219910000087
Indicates the acceleration frequency response function between points s and p in the suspended state;
Figure BDA0001971219910000088
Indicates the acceleration frequency response function between points l and s in the suspended state;
Figure BDA0001971219910000089
Represents the frequency response function of the acceleration origin at point s in the suspension state; as long as the support stiffness k s and support damping c s are known, the frequency response function Alp of the free structure can be obtained by measuring the frequency response function related to the suspension point.

考虑到有些情况下需要对被测结构进行多点悬挂,以下进一步推导了多点悬挂情况下支承影响的消除方法,单点支承情况可视为此处多点支承的一种特例情况。对于多点悬挂的情况,同样可采取施加“负刚度”、“负阻尼”的方法来消除频响函数中的支承影响。Considering that in some cases it is necessary to carry out multi-point suspension on the structure under test, the method for eliminating the influence of support in the case of multi-point suspension is further deduced below, and the case of single-point support can be regarded as a special case of multi-point support here. For the case of multi-point suspension, the method of applying "negative stiffness" and "negative damping" can also be adopted to eliminate the support effect in the frequency response function.

请见图4,针对多点悬挂,假定结构体C的自身刚度足够小而不能忽略悬挂支承的影响,结构体C处于多点悬挂状态,悬挂点分别为1、2、…、n;每个悬挂点支承刚度为ki,阻尼大小为ci,i=1、2、…、n;p点为激励点,l点为测量点,其中“p点”和“l点”为结构体上任意两点;为了消除悬挂支承的影响,依次在每个悬挂点处施加“虚拟支承”(如图5左侧虚线框内的部分),“虚拟支承”中的“负刚度”k′i=-ki和“负阻尼”c′i=-ci分别抵消各悬挂支承的影响,从而得到自由状态的结构体;Please see Figure 4, for multi-point suspension, it is assumed that the self-stiffness of structure C is small enough to ignore the influence of suspension support, structure C is in a multi-point suspension state, and the suspension points are 1, 2, ..., n; each The support stiffness of the suspension point is k i , the damping size is c i , i=1, 2, ..., n; point p is the excitation point, point l is the measurement point, where "point p" and "point l" are the points on the structure Any two points; in order to eliminate the influence of the suspension support, a "virtual support" is applied at each suspension point in turn (as shown in the part in the dotted line box on the left side of Figure 5), and the "negative stiffness"k' i in the "virtual support" = -k i and "negative damping"c' i = -c i respectively counteract the influence of each suspension support, so as to obtain the free state structure;

为便于分析,将悬挂点si点拆分为两个点si点和s′i点,在“虚拟支承”中,用s′i点来表述si点;For the convenience of analysis, the suspension point s i is split into two points s i and s′ i , and in the “virtual support”, s i is used to express s i ;

对原悬挂结构体C添加虚拟支承的过程可视为结构动力学修改的过程。若将原悬挂的结构体C视为子结构Ⅰ,将虚线框内虚拟支承视为子结构Ⅱ,Ⅲ……。对于子结构Ⅰ,可建立频域下的输出位移响应和输入力的关系,The process of adding virtual supports to the original suspension structure C can be regarded as the process of structural dynamics modification. If the original suspended structure C is regarded as substructure I, the virtual support inside the dotted line frame is regarded as substructure II, III.... For substructure I, the relationship between the output displacement response and the input force in the frequency domain can be established,

Figure BDA0001971219910000091
Figure BDA0001971219910000091

Figure BDA0001971219910000092
Figure BDA0001971219910000092

其中,k=1、2、…、n,Xl和Xk分别表示子结构Ⅰ在l点和k点处的位移响应,Fp

Figure BDA0001971219910000093
和Ri分别表示子结构Ⅰ在p点、Si点和i点所受的作用力;
Figure BDA0001971219910000094
表示悬挂状态下l、p点之间的位移频响函数;
Figure BDA0001971219910000095
表示悬挂状态下l、Si点之间的位移频响函数;
Figure BDA0001971219910000096
表示悬挂状态下k、p点之间的位移频响函数;
Figure BDA0001971219910000097
表示悬挂状态下k、i点之间的位移频响函数;Among them, k=1, 2,..., n, X l and X k represent the displacement response of substructure I at point l and point k respectively, F p ,
Figure BDA0001971219910000093
and R i denote the acting force on substructure I at point p, S i and i, respectively;
Figure BDA0001971219910000094
Indicates the displacement frequency response function between points l and p in the suspended state;
Figure BDA0001971219910000095
Indicates the displacement frequency response function between points l and S i in the suspended state;
Figure BDA0001971219910000096
Indicates the displacement frequency response function between points k and p in the suspension state;
Figure BDA0001971219910000097
Indicates the displacement frequency response function between points k and i in the suspended state;

对于子结构Ⅱ、Ⅲ、……,同样建立输出位移响应和输入力的关系,For substructures II, III, ..., the relationship between the output displacement response and the input force is also established,

X′k=α′kkR′k (10)X′ k =α′ kk R′ k (10)

式中,In the formula,

Figure BDA0001971219910000098
Figure BDA0001971219910000098

其中,X′k和Rk分别表示子结构Ⅱ、Ⅲ、……在k点处的位移响应和所受的作用力;α′kk表示子子结构Ⅱ、Ⅲ、……在k点处的原点位移频响函数;jω中j是复数单位,ω为频响函数自变量;Among them, X′ k and R k represent the displacement responses and acting forces of substructures II, III, ... at point k, respectively; Origin displacement frequency response function; j in jω is a complex unit, and ω is an independent variable of frequency response function;

连接点s或s′处力和位移的约束条件分别为:The constraint conditions of force and displacement at connection point s or s′ are respectively:

Rk+R′k=0 (12)R k + R′ k = 0 (12)

Xk=X′k (13)X k = X′ k (13)

联立式(8)~(13),得到Simultaneous formula (8)~(13), get

α(s)·F=X (14)α (s) F=X (14)

式中,In the formula,

Figure BDA0001971219910000101
Figure BDA0001971219910000101

求解式(14),得Solving formula (14), we get

Figure BDA0001971219910000102
Figure BDA0001971219910000102

式中,αlp表示自由结构体l、p点间的位移频响函数;

Figure BDA0001971219910000103
表示悬挂结构体的位移频响函数;α′11,α′22……表示添加的虚拟支承分别在1、2……点处的原点位移频响函数;式(15)表明,对位移频响函数
Figure BDA0001971219910000104
修正需要额外测量与悬挂点1,2…n相关的多个位移频响函数,所需测量的位移频响函数数量随着悬挂点的增多而增加;In the formula, α lp represents the displacement frequency response function between points l and p of the free structure;
Figure BDA0001971219910000103
Indicates the displacement frequency response function of the suspended structure; α′ 11 , α′ 22 ... represent the origin displacement frequency response function of the added virtual support at points 1, 2... respectively; Equation (15) shows that the displacement frequency response function
Figure BDA0001971219910000104
The correction requires additional measurement of multiple displacement frequency response functions related to suspension points 1, 2...n, and the number of displacement frequency response functions required to be measured increases with the number of suspension points;

从式(15)推算出,对于N点悬挂的情况,修正

Figure BDA0001971219910000105
需要额外测量(N2+5N)/2组位移频响函数,分别为原点位移频响函数
Figure BDA0001971219910000106
和跨点位移频响函数
Figure BDA0001971219910000107
这些位移频响函数均与悬挂点相关;以上结论是针对激励点p、响应点l和悬挂点1,2…n均不重合的情况,当它们其中某些点重合时,式(15)会适当简化,修正工作中所需测量的位移频响函数也会适当减少;It is deduced from formula (15), for the case of N-point suspension, the correction
Figure BDA0001971219910000105
Additional measurements (N 2 +5N)/2 sets of displacement frequency response functions are required, which are the origin displacement frequency response functions
Figure BDA0001971219910000106
and across point displacement frequency response function
Figure BDA0001971219910000107
These displacement frequency response functions are all related to the suspension point; the above conclusion is for the case that the excitation point p, the response point l and the suspension point 1, 2...n do not overlap. When some of them overlap, the formula (15) will be Appropriate simplification, the displacement frequency response function that needs to be measured in the correction work will also be appropriately reduced;

对于加速度频响的情况,只需将式(15)两端分别乘以-ω2即得到加速度频响函数A的表达形式,For the case of acceleration frequency response, it is only necessary to multiply both ends of formula (15) by -ω 2 to obtain the expression form of the acceleration frequency response function A,

Figure BDA0001971219910000111
Figure BDA0001971219910000111

式中A均代表加速度频响函数,Alp表示自由结构体l、p点间的加速度频响函数;

Figure BDA0001971219910000112
表示悬挂状态下l、p点之间的加速度频响函数;
Figure BDA0001971219910000113
表示悬挂状态下n点加速度原点频响函数;A′11,A′22……A′nn分别表示添加的虚拟支承在1、2……n点处的加速度频响函数。In the formula, A represents the acceleration frequency response function, and A lp represents the acceleration frequency response function between points l and p of the free structure;
Figure BDA0001971219910000112
Indicates the acceleration frequency response function between points l and p in the suspended state;
Figure BDA0001971219910000113
Indicates the frequency response function of the acceleration origin of point n in the suspension state; A′ 11 , A′ 22 ... A' nn represent the acceleration frequency response function of the added virtual support at points 1, 2...n, respectively.

为了验证多点悬挂支承影响消除方法的可行性,选用两点悬挂的自由梁模态测试。如图6,自由梁沿长度方向均分五个测点,其中的1、5点为悬挂点,每个悬挂点处由一个支承弹簧和一个支承阻尼等效。梁的尺寸和物理参数分别为:1×0.04×0.003m3和7850kg/m3。系统参数选为:k1=600N/m,k5=300N/m,c1=0.3N·s/m,c5=0.5N·s/m。选择梁的第4点为激励点,第2点为响应测量点。本实验目标是消除频响函数

Figure BDA0001971219910000114
中的支承影响。In order to verify the feasibility of the multi-point suspension support effect elimination method, the free beam modal test with two-point suspension is selected. As shown in Figure 6, the free beam is equally divided into five measuring points along the length direction, of which points 1 and 5 are suspension points, and each suspension point is equivalent to a support spring and a support damper. The dimensions and physical parameters of the beam are: 1×0.04×0.003m 3 and 7850kg/m 3 , respectively. The system parameters are selected as: k 1 =600N/m, k 5 =300N/m, c 1 =0.3N·s/m, c 5 =0.5N·s/m. Select the fourth point of the beam as the excitation point, and the second point as the response measurement point. The goal of this experiment is to eliminate the frequency response function
Figure BDA0001971219910000114
support effects in .

可根据式(16)对

Figure BDA0001971219910000115
进行修正,修正后的频响函数为
Figure BDA0001971219910000116
According to formula (16) to
Figure BDA0001971219910000115
After correction, the frequency response function after correction is
Figure BDA0001971219910000116

Figure BDA0001971219910000117
Figure BDA0001971219910000117

上式右边

Figure BDA0001971219910000118
为需要“测量的”频响函数,这里可直接对悬挂梁有限元建模计算获取。为了检验修正效果,需要将修正后的频响函数
Figure BDA0001971219910000119
和准确频响函数A24进行对比,而A24可通过自由梁有限元建模计算获取。如图7,相比于准确频响函数A24,“测量的”频响函数
Figure BDA00019712199100001110
在悬挂支承的影响下,三阶共振峰均向右移动,并且在原第一阶共振峰左侧多出两个共振峰(约3.8Hz和7Hz处)。同时可看出,三阶共振峰的移动幅度,第一阶最大,第二节其次,第三阶最小。当消除悬挂支承影响后,修正后的频响函数
Figure BDA0001971219910000121
与准确的频响函数A24完全一致,如图8所示。但如果仅消除支承影响中的支承刚度影响,而忽略支承阻尼,修正后的效果如图9所示,总体修正效果较好,但在共振峰处
Figure BDA0001971219910000122
幅值略小于A24,反共振峰处
Figure BDA0001971219910000123
幅值略大于A24,这种误差将导致提取的模态阻尼比不准确。因此,对于高质量模态测试,既要消除支承刚度的影响,还需考虑支承阻尼的影响。The right side of the above formula
Figure BDA0001971219910000118
For the "measured" frequency response function, it can be directly obtained by modeling and calculating the finite element of the suspension beam. In order to test the correction effect, the corrected frequency response function needs to be
Figure BDA0001971219910000119
It is compared with the accurate frequency response function A 24 , and A 24 can be obtained through free beam finite element modeling calculation. As shown in Figure 7, compared to the exact FRF A 24 , the “measured” FRF
Figure BDA00019712199100001110
Under the influence of the suspension support, the third-order formant moves to the right, and there are two more formants (about 3.8Hz and 7Hz) on the left side of the original first-order formant. At the same time, it can be seen that the movement range of the third-order formant is the largest in the first order, followed by the second, and the smallest in the third. After eliminating the effect of the suspension support, the modified frequency response function
Figure BDA0001971219910000121
It is completely consistent with the accurate frequency response function A 24 , as shown in FIG. 8 . However, if only the support stiffness effect in the support effect is eliminated, and the support damping is ignored, the corrected effect is shown in Figure 9. The overall correction effect is good, but at the resonance peak
Figure BDA0001971219910000122
The amplitude is slightly smaller than A 24 , at the anti-resonant peak
Figure BDA0001971219910000123
The magnitude is slightly larger than A 24 , and this error will lead to inaccurate extracted modal damping ratios. Therefore, for high-quality modal testing, it is necessary to eliminate the influence of support stiffness and also consider the influence of support damping.

悬挂支承是模态测试实验中一种常用的支承方式,主要用于模拟自由边界条件。但这种支承方式引入的附加刚度和阻尼会给测量的频响函数带来影响,尤其当被测结构件刚度较小时,这种影响影响十分显著。本发明针对悬挂支承方式的自由模态测试,提出一种从测量的频响函数中消除支承影响的方法,支承条件中同时考虑了支承刚度和支承阻尼。结果表明:在给定支承刚度ks和支承阻尼cs的情况下,通过测量与悬挂点相关的频响函数,即可消除悬挂支承刚度和阻尼的影响,获取自由结构的频响函数。Suspension support is a commonly used support method in modal test experiments, and it is mainly used to simulate free boundary conditions. However, the additional stiffness and damping introduced by this support method will have an impact on the measured frequency response function, especially when the stiffness of the measured structural member is small, the impact is very significant. Aiming at the free mode test of the suspension support mode, the invention proposes a method for eliminating the influence of the support from the measured frequency response function, and the support stiffness and support damping are considered in the support condition. The results show that: given the support stiffness k s and support damping c s , by measuring the frequency response function related to the suspension point, the influence of the suspension support stiffness and damping can be eliminated, and the frequency response function of the free structure can be obtained.

应当理解的是,本说明书未详细阐述的部分均属于现有技术。It should be understood that the parts not described in detail in this specification belong to the prior art.

应当理解的是,上述针对较佳实施例的描述较为详细,并不能因此而认为是对本发明专利保护范围的限制,本领域的普通技术人员在本发明的启示下,在不脱离本发明权利要求所保护的范围情况下,还可以做出替换或变形,均落入本发明的保护范围之内,本发明的请求保护范围应以所附权利要求为准。It should be understood that the above-mentioned descriptions for the preferred embodiments are relatively detailed, and should not therefore be considered as limiting the scope of the patent protection of the present invention. Within the scope of protection, replacements or modifications can also be made, all of which fall within the protection scope of the present invention, and the scope of protection of the present invention should be based on the appended claims.

Claims (4)

1.一种针对自由模态测试消除频响函数中悬挂支承影响的方法,其特征在于:针对单点悬挂,假定结构体C的自身刚度足够小而不能忽略悬挂支承的影响,结构体C悬挂于s点处,悬挂支承效果等效为ks刚度的弹簧和大小为cs的阻尼,p点为激励点,l点为测量点,其中“p点”和“l点”为结构体上任意两点;为了消除悬挂支承的影响,在s点施加“虚拟支承”,“虚拟支承”中的“负刚度”k′s=-ks和“负阻尼”c′s=-cs抵消了悬挂支承的影响,使结构体等效为自由状态;1. A method for eliminating the influence of suspension support in the frequency response function for free mode testing, characterized in that: for single-point suspension, assuming that the self stiffness of structure C is small enough to ignore the influence of suspension support, structure C suspension At point s, the suspension support effect is equivalent to the spring with k s stiffness and the damping with c s , point p is the excitation point, and point l is the measurement point, where "p point" and "l point" are the Any two points; in order to eliminate the influence of the suspension support, a "virtual support" is applied at point s, and the "negative stiffness" k′ s = -k s and "negative damping" c′ s = -c s in the "virtual support" cancel The influence of the suspension support is eliminated, so that the structure is equivalent to a free state; 为便于分析,将s点拆分为两个点s点和s′点,在“虚拟支承”中,用s′点来表述s点;For the convenience of analysis, point s is split into two points, s and s′, and in the “virtual support”, point s is expressed by point s′; 视悬挂的结构体C为子结构Ⅰ,“虚拟支承”为子结构Ⅱ,则子结构Ⅰ和子结构Ⅱ结合即得到自由结构体;Assuming that the suspended structure C is substructure I, and the "virtual support" is substructure II, the combination of substructure I and substructure II results in a free structure; 对子结构Ⅰ和Ⅱ分别进行分析;Analyze substructures I and II separately; 对于子结构Ⅰ,建立频域下的输出位移响应和输入力的关系,For substructure Ⅰ, the relationship between the output displacement response and the input force in the frequency domain is established,
Figure FDA0001971219900000011
Figure FDA0001971219900000011
其中,Xl和Xs分别表示子结构Ⅰ在l点和s点处的位移响应,Fp和Rs分别表示子结构Ⅰ在p点和s点所受的作用力;
Figure FDA0001971219900000012
表示悬挂状态下l、s点之间的位移频响函数;
Figure FDA0001971219900000013
表示悬挂状态下s、p点之间的位移频响函数;
Figure FDA0001971219900000014
表示悬挂状态下s点处的原点位移频响函数;
Figure FDA0001971219900000015
表示悬挂状态下l、p点之间的位移频响函数;
Among them, X l and X s represent the displacement response of substructure I at points l and s, respectively, and F p and R s represent the forces on substructure I at points p and s, respectively;
Figure FDA0001971219900000012
Indicates the displacement frequency response function between points l and s in the suspended state;
Figure FDA0001971219900000013
Indicates the displacement frequency response function between points s and p in the suspended state;
Figure FDA0001971219900000014
Indicates the origin displacement frequency response function at point s in the suspended state;
Figure FDA0001971219900000015
Indicates the displacement frequency response function between points l and p in the suspended state;
对于子结构Ⅱ,同样建立输出位移响应和输入力的关系,For substructure II, the relationship between the output displacement response and the input force is also established, X′s′=αs′s′Rs′ (2)X's' = α s's' R s' (2) 式中,In the formula,
Figure FDA0001971219900000016
Figure FDA0001971219900000016
其中,X′s′和Rs′分别表示子结构Ⅱ在s点处的位移响应和所受的作用力;αs′s′表示子结构Ⅱ在s点处的原点位移频响函数;jω中j是复数单位,ω为频响函数自变量;Among them, X′ s′ and R s′ represent the displacement response and the acting force of substructure II at point s respectively; α s′s′ represents the frequency response function of the origin displacement of substructure II at point s; jω where j is a complex unit, and ω is an independent variable of the frequency response function; 连接点s或s′处力和位移的约束条件分别为The constraint conditions of force and displacement at connection point s or s′ are respectively Rs+Rs′=0 (4)R s +R s′ = 0 (4) Xs=Xs′ (5)X s =X s′ (5) 联立式(1)~(5),得到Simultaneous formulas (1)~(5), get
Figure FDA0001971219900000021
Figure FDA0001971219900000021
式中,αlp表示自由状态下l、p点之间的位移频响函数;
Figure FDA0001971219900000022
表示悬挂状态下l、p点之间的位移频响函数;为了避免混淆,命名自由状态下的频响函数为“准确频响函数”,悬挂状态下的频响为“实测频响函数”,即αlp为准确频响函数,
Figure FDA0001971219900000023
Figure FDA0001971219900000024
均为实测频响函数;
In the formula, α lp represents the displacement frequency response function between points l and p in free state;
Figure FDA0001971219900000022
Indicates the displacement frequency response function between points l and p in the suspension state; in order to avoid confusion, the frequency response function in the free state is named "accurate frequency response function", and the frequency response in the suspension state is "measured frequency response function", That is, α lp is the exact frequency response function,
Figure FDA0001971219900000023
Figure FDA0001971219900000024
Both are measured frequency response functions;
根据式(6)可知,准确频响函数αlp和实测频响函数
Figure FDA0001971219900000025
有所差异,而这种差异正是由于支承的刚度和阻尼导致;而且在悬挂刚度ks和阻尼cs给定的情况下,通过测量与悬挂点s相关的频响函数可计算得到准确频响函数αlp,即达到消除频响函数
Figure FDA0001971219900000026
中支承影响的目的。
According to formula (6), it can be seen that the exact frequency response function α lp and the measured frequency response function
Figure FDA0001971219900000025
There is a difference, and this difference is caused by the stiffness and damping of the support; and when the suspension stiffness k s and damping c s are given, the exact frequency can be calculated by measuring the frequency response function related to the suspension point s Response function α lp , that is, to eliminate the frequency response function
Figure FDA0001971219900000026
The purpose of supporting influence.
2.根据权利要求1所述的针对自由模态测试消除频响函数中悬挂支承影响的方法,其特征在于:对于加速度频响的情况,需将式(6)左右乘以-ω2,即得到加速度频响和准确加速度频响之间的关系式(7)2. the method for eliminating the influence of suspension support in the frequency response function for free mode testing according to claim 1 is characterized in that: for the situation of acceleration frequency response, formula (6) needs to be multiplied by -ω 2 left and right, namely Get the relationship between acceleration frequency response and accurate acceleration frequency response (7)
Figure FDA0001971219900000027
Figure FDA0001971219900000027
式中,A均代表加速度频响函数,Alp表示自由结构体l、p点间的加速度频响函数;
Figure FDA0001971219900000028
表示悬挂状态下l、p点之间的加速度频响函数;
Figure FDA0001971219900000029
表示悬挂状态下s、p点之间的加速度频响函数;
Figure FDA00019712199000000210
表示悬挂状态下l、s点之间的加速度频响函数;
Figure FDA00019712199000000211
表示悬挂状态下s点加速度原点频响函数;只要支承刚度ks和支承阻尼cs大小已知,就能通过测量与悬挂点相关的频响函数来获取自由结构的频响函数Alp
In the formula, A represents the acceleration frequency response function, and A lp represents the acceleration frequency response function between points l and p of the free structure;
Figure FDA0001971219900000028
Indicates the acceleration frequency response function between points l and p in the suspended state;
Figure FDA0001971219900000029
Indicates the acceleration frequency response function between points s and p in the suspended state;
Figure FDA00019712199000000210
Indicates the acceleration frequency response function between points l and s in the suspended state;
Figure FDA00019712199000000211
Represents the frequency response function of the acceleration origin at point s in the suspension state; as long as the support stiffness k s and support damping c s are known, the frequency response function Alp of the free structure can be obtained by measuring the frequency response function related to the suspension point.
3.一种针对自由模态测试消除频响函数中悬挂支承影响的方法,其特征在于:针对多点悬挂,假定结构体C的自身刚度足够小而不能忽略悬挂支承的影响,结构体C处于多点悬挂状态,悬挂点分别为1、2、…、n;每个悬挂点支承刚度为ki,阻尼大小为ci,i=1、2、…、n;p点为激励点,l点为测量点,其中“p点”和“l点”为结构体上任意两点;为了消除悬挂支承的影响,依次在每个悬挂点处施加“虚拟支承”,“虚拟支承”中的“负刚度”ki′=-ki和“负阻尼”c′i=-ci分别抵消各悬挂支承的影响,从而得到自由状态的结构体;3. A method for eliminating the influence of suspension supports in the frequency response function for free mode testing, characterized in that: for multi-point suspension, it is assumed that the self-stiffness of structure C is small enough to ignore the influence of suspension supports, and structure C is in In the multi-point suspension state, the suspension points are 1, 2, ..., n respectively; the support stiffness of each suspension point is ki , the damping size is ci, i =1, 2, ..., n; point p is the excitation point, l point is the measurement point, where "point p" and "point l" are any two points on the structure; in order to eliminate the influence of suspension support, a "virtual support" is applied at each suspension point in turn, and the "virtual support" in "virtual support" Negative stiffness " ki '=- ki and "negative damping"c' i = -ci respectively counteract the influence of each suspension support, so as to obtain a structure in a free state; 为便于分析,将悬挂点si点拆分为两个点si点和s′i点,在“虚拟支承”中,用s′i点来表述si点;For the convenience of analysis, the suspension point s i is split into two points s i and s′ i , and in the “virtual support”, s i is used to express s i ; 将原悬挂的结构体C视为子结构Ⅰ,将虚拟支承视为子结构Ⅱ、Ⅲ、……;对于子结构Ⅰ,建立频域下的输出位移响应和输入力的关系,The original suspended structure C is regarded as substructure I, and the virtual support is regarded as substructure II, III, ...; for substructure I, the relationship between the output displacement response and the input force in the frequency domain is established,
Figure FDA0001971219900000031
Figure FDA0001971219900000031
Figure FDA0001971219900000032
Figure FDA0001971219900000032
其中,k=1、2、…、n,Xl和Xk分别表示子结构Ⅰ在l点和k点处的位移响应,Fp
Figure FDA0001971219900000033
和Ri分别表示子结构Ⅰ在p点、Si点和i点所受的作用力;
Figure FDA0001971219900000034
表示悬挂状态下l、p点之间的位移频响函数;
Figure FDA0001971219900000035
表示悬挂状态下l、Si点之间的位移频响函数;
Figure FDA0001971219900000036
表示悬挂状态下k、p点之间的位移频响函数;
Figure FDA0001971219900000037
表示悬挂状态下k、i点之间的位移频响函数;
Among them, k=1, 2,..., n, X l and X k represent the displacement response of substructure I at point l and point k respectively, F p ,
Figure FDA0001971219900000033
and R i denote the acting force on substructure I at point p, S i and i, respectively;
Figure FDA0001971219900000034
Indicates the displacement frequency response function between points l and p in the suspended state;
Figure FDA0001971219900000035
Indicates the displacement frequency response function between points l and S i in the suspended state;
Figure FDA0001971219900000036
Indicates the displacement frequency response function between points k and p in the suspension state;
Figure FDA0001971219900000037
Indicates the displacement frequency response function between points k and i in the suspended state;
对于子结构Ⅱ、Ⅲ、……,同样建立输出位移响应和输入力的关系,For substructures II, III, ..., the relationship between the output displacement response and the input force is also established, X′k=α′kkR′k (10)X′ k =α′ kk R′ k (10) 式中,In the formula,
Figure FDA0001971219900000038
Figure FDA0001971219900000038
其中,X′k和Rk分别表示子结构Ⅱ、Ⅲ、……在k点处的位移响应和所受的作用力;α′kk表示子子结构Ⅱ、Ⅲ、……在k点处的原点位移频响函数;jω中j是复数单位,ω为频响函数自变量;Among them, X′ k and R k represent the displacement responses and acting forces of substructures II, III, ... at point k, respectively; Origin displacement frequency response function; j in jω is a complex unit, and ω is an independent variable of frequency response function; 连接点s或s′处力和位移的约束条件分别为:The constraint conditions of force and displacement at connection point s or s′ are respectively: Rk+R′k=0 (12)R k + R′ k = 0 (12) Xk=X′k (13)X k = X′ k (13) 联立式(8)~(13),得到Simultaneous formula (8)~(13), get α(s)·F=X (14)α (s) F=X (14) 式中,In the formula,
Figure FDA0001971219900000041
Figure FDA0001971219900000041
求解式(14),得Solving formula (14), we get
Figure FDA0001971219900000042
Figure FDA0001971219900000042
式中,αlp表示自由结构体l、p点间的位移频响函数;
Figure FDA0001971219900000043
……表示悬挂结构体的位移频响函数;α′11,α′22……表示添加的虚拟支承分别在1、2……点处的原点位移频响函数;式(15)表明,对位移频响函数
Figure FDA0001971219900000044
修正需要额外测量与悬挂点1,2…n相关的多个位移频响函数,所需测量的位移频响函数数量随着悬挂点的增多而增加;
In the formula, α lp represents the displacement frequency response function between points l and p of the free structure;
Figure FDA0001971219900000043
...represents the displacement frequency response function of the suspended structure; α′ 11 , α′ 22 ...represents the origin displacement frequency response function of the added virtual support at points 1, 2... respectively; Equation (15) shows that for displacement frequency response function
Figure FDA0001971219900000044
The correction requires additional measurement of multiple displacement frequency response functions related to suspension points 1, 2...n, and the number of displacement frequency response functions required to be measured increases with the number of suspension points;
从式(15)推算出,对于N点悬挂的情况,修正
Figure FDA0001971219900000045
需要额外测量(N2+5N)/2组位移频响函数,分别为原点位移频响函数
Figure FDA0001971219900000051
和跨点位移频响函数
Figure FDA0001971219900000052
这些位移频响函数均与悬挂点相关;以上结论是针对激励点p、响应点l和悬挂点1,2…n均不重合的情况,当它们其中某些点重合时,式(15)会适当简化,修正工作中所需测量的位移频响函数也会适当减少。
It is deduced from formula (15), for the case of N-point suspension, the correction
Figure FDA0001971219900000045
Additional measurements (N 2 +5N)/2 sets of displacement frequency response functions are required, which are the origin displacement frequency response functions
Figure FDA0001971219900000051
and across point displacement frequency response function
Figure FDA0001971219900000052
These displacement frequency response functions are all related to the suspension point; the above conclusion is for the case that the excitation point p, the response point l and the suspension point 1, 2...n do not overlap. When some of them overlap, the formula (15) will be With appropriate simplification, the frequency response function of displacement that needs to be measured in the correction work will also be appropriately reduced.
4.根据权利要求3所述的针对自由模态测试消除频响函数中悬挂支承影响的方法,其特征在于:对于加速度频响的情况,只需将式(15)两端分别乘以-ω2即得到加速度频响函数A的表达形式,4. the method for eliminating the influence of suspension support in the frequency response function according to claim 3, is characterized in that: for the situation of acceleration frequency response, only need to multiply the two ends of formula (15) by -ω respectively 2 That is, the expression form of the acceleration frequency response function A is obtained,
Figure FDA0001971219900000053
Figure FDA0001971219900000053
式中A均代表加速度频响函数,Alp表示自由结构体l、p点间的加速度频响函数;
Figure FDA0001971219900000054
表示悬挂状态下l、p点之间的加速度频响函数;
Figure FDA0001971219900000055
表示悬挂状态下n点处加速度原点频响函数;A′11,A′22……A′nn分别表示添加的虚拟支承在1、2……n点处的原点加速度频响函数。
In the formula, A represents the acceleration frequency response function, and A lp represents the acceleration frequency response function between points l and p of the free structure;
Figure FDA0001971219900000054
Indicates the acceleration frequency response function between points l and p in the suspended state;
Figure FDA0001971219900000055
Indicates the frequency response function of the acceleration origin at point n in the suspension state; A′ 11 , A′ 22 ... A′ nn represent the acceleration frequency response functions of the added virtual support at points 1, 2...n respectively.
CN201910119087.7A 2019-02-18 2019-02-18 Method for eliminating suspension support influence in frequency response function aiming at free mode test Active CN109726521B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910119087.7A CN109726521B (en) 2019-02-18 2019-02-18 Method for eliminating suspension support influence in frequency response function aiming at free mode test

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910119087.7A CN109726521B (en) 2019-02-18 2019-02-18 Method for eliminating suspension support influence in frequency response function aiming at free mode test

Publications (2)

Publication Number Publication Date
CN109726521A CN109726521A (en) 2019-05-07
CN109726521B true CN109726521B (en) 2022-12-16

Family

ID=66301523

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910119087.7A Active CN109726521B (en) 2019-02-18 2019-02-18 Method for eliminating suspension support influence in frequency response function aiming at free mode test

Country Status (1)

Country Link
CN (1) CN109726521B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110569568A (en) * 2019-08-21 2019-12-13 湖北工业大学 A Method to Eliminate the Effects of Suspension Stiffness and Sensor Mass in Free Mode Testing

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013134853A1 (en) * 2012-03-13 2013-09-19 The Royal Institution For The Advancement Of Learning/Mcgill University Methods and systems for optimizing frequency modulation atomic force microscopy
US20170297402A1 (en) * 2016-04-13 2017-10-19 GM Global Technology Operations LLC Detection and reconstruction of suspension height sensor faults
CN106871933B (en) * 2017-03-03 2019-01-29 湖北工业大学 A kind of method of multisensor additional mass in elimination frequency response function
CN108489696B (en) * 2018-03-26 2019-03-12 东南大学 A kind of method of spring additional stiffness in elimination frequency response function

Also Published As

Publication number Publication date
CN109726521A (en) 2019-05-07

Similar Documents

Publication Publication Date Title
Larsen et al. Modal analysis of wind turbine blades
Addari et al. Experimental and numerical investigation of coupled microvibration dynamics for satellite reaction wheels
CN106323528B (en) A kind of high-precision dynamic measurement of pile of rope tensility
Zhang et al. Nonlinear vibration characteristics of bridge deck section models in still air
Puel et al. Parameter identification of nonlinear time-dependent rubber bushings models towards their integration in multibody simulations of a vehicle chassis
KR102287994B1 (en) Methods and systems for testing coupled hybrid dynamic system
CN106052743B (en) A kind of sensor mass of assessing influences frequency response function the method for size
Dai et al. Robust control of vortex‐induced vibration in flexible bridges using an active tuned mass damper
CN109726521B (en) Method for eliminating suspension support influence in frequency response function aiming at free mode test
Zhu et al. Removing mass loading effects of multi-transducers using Sherman-Morrison-Woodbury formula in modal test
US11371923B2 (en) Automatic system compliance estimation and correction for mechanical testing systems
CN106871933A (en) A kind of method of multisensor additional mass in elimination frequency response function
CN106407599A (en) Semi-active suspension amplitude variation characteristic modeling method
JP5203851B2 (en) Rigidity evaluation support device, rigidity evaluation support method and program
Faverjon et al. Identification of an open crack in a beam using an a posteriori error estimator of the frequency response functions with noisy measurements
Helsen et al. Global static and dynamic car body stiffness based on a single experimental modal analysis test
JP2016529528A (en) Test system with compliant actuator assembly and repetitively obtained drive
Kranjc et al. An interface force measurements-based substructure identification and an analysis of the uncertainty propagation
CN110569568A (en) A Method to Eliminate the Effects of Suspension Stiffness and Sensor Mass in Free Mode Testing
Kashani et al. Identification of dynamic characteristics of nonlinear joint based on the optimum equivalent linear frequency response function
Maheux et al. Nonlinear Wind Tunnel Tests of Cable-Supported Bridges
CN106643847B (en) A Design Method of Universal Elastic Vibration Isolation Device for Sensor Mounting Measurement
Klaus et al. Model parameter identification from measurement data for dynamic torque calibration
Flore et al. Dynamic Mechanical Analysis of an Aircraft Wing with emphasis on vibration modes change with loading
JP6159954B2 (en) Method for evaluating characteristics of connecting members

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant