螺栓联接梁结构结合部非线性特性参数辨识

本资料包含pdf文件1个，下载需要1积分

螺栓联接是机械结构系统最常见的联接形式，通常情况下，螺栓联接结构的结合部会对系统的动力学特性产生很大的影响。与刚性联接对照，结合部的力学特性参数会使结构的固有频率、刚度和阻尼等都会发生改变，而表现出非线性动力学特性。

因此，为了有效创建螺栓联接结构系统的动力学模型以及对振动响应进行精确预测，需要辨识结合部的动力学特性参数(主要是刚度和阻尼参数)。

近年来关于结合部参数辨识的研究有很多，通常采用实验数据与有限元模型相结合，即通过最小化实验测试数据和有限元分析模型预测结果之间误差的方法来获得最合适的结合部特性参数[1 ]。但是，由于有限元方法及实验测试都存在着误差 ，误差的相互叠加使这种参数辨识方法很难精确辨识出结合部的力学特性参数。

鉴于这种情况，很多学者从不同的角度提出建立各种解析模型来辨识螺栓结合部特性参数的方法，例如，Ma等提出通过对比无螺栓联接和有螺栓联接的梁的动力学特性 ]，来辨识结合部的非线性参数 。Ahmadian提出建立-种通用联接单元 的方式来表示螺栓联接处的非线性 ]，并将这种单元引入到结构控制方程中，最后通过实验测试获得的线收稿 日期 ：2012-03-30；修订 日期 ：2012-12-19基金项目：国家自然科学基金资助项 目(50905029)性固有频率、非线性频响函数来辨识结构的阻尼和刚度。Lee将非线性结构的响应表示为 Volterra级数[6]，通过对高阶频响函数的分析，辨识出结构的线性和非线性参数 。由于对螺栓联接结合部非线性产生机理的认识不清，上述对螺栓结合部参数辨识的研究还未得到公认，因而，围绕螺栓结合部参数辨识还需要深人开展研究 。

力状态映射法(Force-state mapping)是通过在某阶固有频率附近定频激励结构，测得其激励和响应的时域信号，根据力学模型中的函数关系对非线性参数进行辨识的-种方法l7]。按照模态理论，相互耦合的多自由动力学系统均可以解耦为若干个单自由度系统，因而力状态映射法可识别对应各阶次的非线性特性参数，有着较为广泛的应用。例如，Edward等利用 Force-state mapping方法识别了-些强 非线 性 结构 的非 线 性模 型参 数 ]，Hadid对Force-state mapping方法质量和等效质量的计算方法进行了分析[g]，提出-种基于线性拟合和迭代的方法来获得非线性模型中准确的等效质量。

本文以螺栓联接组合梁为例，研究采用力状态映射法识别螺栓组合梁非线性结合部特性参数的方法。在分析结合部非线性模型以及对非线性结合部动力学特性 表征测试 的基础 上，提 出采 用 简化 的1wan模型模拟接触面之间的非线性动力学现象，从而建立了螺栓联接组合梁非线性解析模型▲-g 担-学E程工V动m振。

振 动 工 程 学 报 第 26卷步，提出采用力状态映射法识别螺栓组合梁结合部的动力学特性参数的方法及相应的测试流程。最后，进行了实例研究，应用所提出的方法，辨识出了所研究的螺栓组合梁结合部的刚度及阻尼特性参数。

1 螺栓结合部非线性模型及非线性表征测试1.1 螺栓结合部非线性模型在研究螺栓联接结构的非线性模型时，许多学者对结合部的物理机理进行了分析，并提出多种模型来表示这种物理现象。例如，Ferri利用弹簧和粘性摩擦阻尼组成的滞回系统来表示联接部结合面之间的接触问题[1。。。Gaul和 Segalman概 括了-些可以描述联接结构的模型，如 Jenkins模型、1wan模型、LuGre模型和Valanis模型等[1 ]，可用以描述螺栓联接结合部非线性产生机理 。

其中并联 1wan模型由于物理概念清晰，便于建模得到了较为广泛的应用[1 ]。1wan模型将螺栓结合部的力学特性视为若干弹簧 k 和库伦摩擦阻尼力 F 的组合。对于每-个子系统，当外载荷较小时 ，整个系统可以认为是黏着运动 ，随着载荷的逐渐增加，子系统开始克服库伦摩擦，逐渐蜕化为完全的滑移运动，且临界状态会出现拍击现象。因此，并联 1wan模型可以很好地模拟联接单元的接触面在外载荷作用下产生的滑移和拍击等非线性行为。

1.2 螺栓结合部非线性表征测试用实验来描述 螺栓结合部 的非线性 动力 学行为 ，即螺栓结合部 的非线性表征测试。图 1为本文所研究 的螺栓组合梁表征测试实验系统 。

测试不同激振力幅下的频响函数是表征结构系统非线性特征的最有效方式之-l1 。激振力幅值由2 N变换到 4 N时，在-阶固有频率附近螺栓组合梁的非线性 FRF曲线见 图 2。从 图 2可以看出，螺栓联接结构的 1阶固有频率随着激振力幅值的增加减小，表现出明显的软化效应。

结合螺柱联接梁的非线性表征实验，1wan模型所描述的力-位移关系可进-步简化为[5]e(Ay)- kyAy- kyNAy。 (1)这里将线性结构中的刚度方程转化为引入库伦摩擦和滑移之后的包含立方项的刚度方程，其中e为非线性弹性力，k 和k 分别表示线性项和立方项的刚度系数。

(a)测试实物图(a)Test system(b)测试方案图(单位：ram)(b)Test schem e(Unit：ram)图 1 螺栓组合梁非线性表征测试系统Fig.1 Nonlinear characterization test system for the bol-ted built-up structure频率/Hz图 2 不同激振力幅下螺栓联接结构的频响函数Fig.2 The frequency response function of bolted built-upstructure at the different exciting levels2 螺栓联接梁结构非线性动力学模型的建立螺栓组合梁非线性动力学建模的思路是：首先建立起结构 的线 性模型 ，进 -步采用 简化 的 1wan模型模拟接触面之间的非线性动力学现象，最后将非线性项引入到动力学模型中，完成螺栓组合梁的非线性动力学模型的创建。

2.1 螺栓组合梁的线性动力学模型创建图3(a)为由两根材料及几何尺寸相同的梁通过单个螺栓联接组成螺栓联接梁结构。图 3(b)为螺栓联接梁的动力学模型，联接部处包含了联接刚第 2期 孙 伟，等：螺栓联接梁结构结合部非线性特性参数辨识度及阻尼系数。

(a)螺栓联接梁结构(a)The bolted built-up s仇lctIlre(b)螺栓联接梁的动力学模型(b)The dynamics model ofthe bolted beam s协lcture图 3 螺栓联接梁结构及 简化动力 学模 型Fig.3 The bolted beam structure and its simplified dy-namics model图 3中.D，A，E和 分别代表梁的密度、截面面积、杨氏模量和惯性矩，愚和C则表示刚度及阻尼系数。Y。( ， )和 Y (z，)分别表示两个梁的横向位移。

在外界激励载荷较小时，螺栓联接结合部处不发生相对位移 ，可视为刚性联接 。由此，按照欧拉伯努利梁理论，并引入边界条件及连续性条件，图3所示的梁结构在 1阶固有频率附近的模态坐标运动方程可表示为口( dz z dx)·q( ) (j。 (x)dx )dz) -I 1 (z)f1(z， )dxI 1 (z)f2(z，)dx (2)式中 (z)和 (z)分别表示左侧及右侧联接梁的 1阶模态振型函数，q( )为模态坐标，厂 (z， )和-厂2(z，)分别是作用在两根梁的激振力。

式(2)左侧两项分别是弹性力项和惯性力项，而式(2)右侧为激振力项，由此可将方程(1)进-步表示为m k，q( )-Q (3)这里 m 为等效模态质量，志 为等效模态刚度，Q(t)为等效模态力。相应项的具体表达式为m e - (J。 )dzJf/2 )dz)(4)肼 )。dz田 )。dz(5)rl/2 rlQ-j。91l(z)厂 (z，t)dx r( )厂z( ，t)dx(6)2.2 螺栓组合梁的非线性动力学模型创建在螺栓联接结合部处，当外界载荷的幅值超过某-临界值时 ，联接单元的变形变大 ，在联接部的结合面上将产生滑移和拍击的现象。此时结合部处的刚度及阻尼效应将显现出来，且刚度值不为常数，而变成相对位移 Ay的函数。相对于方程(3)，此时的方程将会出现非线性弹性力项 e。

结合试验中观察到的软式非线性现象，并考虑结构的阻尼，按照简化的1wan模型，最终获得螺栓联接梁在模态坐标下的非线性动力学方程为( ) c，q( ) k，q( )- k3q。( )-Q (7)式中 c 为等效阻尼，k。为立方刚度项系数♂合部的线性刚度也融入到等效刚度 k 中，进-步通过正则变换，获得在正则坐标下的螺栓联接梁方程的表达为( ) 口 。q- 是3q。( )-Q (8)方程(8)中的 愚。及 C 描述了螺栓组合梁结合部的非线性刚度及阻尼参数，接下来将研究这些参数的辨识方法 。

3 力状态映射法及面向组合梁的结合部动力学参数测试流程3.1 力状态映射法及时域数据处理方法方程(8)为典型的 Duffing系统，按照力状态映射法，可变换为rq( 1) -q。( 1) q(t1)lq(2) -q。(2) Q(t2)l i i ；lq( ) -q。(t ) 口( )是3 -l I C目 J(9)这里 t ，t。，，t 为时间序列。

进-步，可得螺栓连接结构非线性特性参数的辨识方程是3 -- q。( 1) q(t1)- q。(2) q(t2)- q。( ) q(t ))-)-)-(10)188 振 动 工 程 学 报 第 26卷这里 ，(·) 表示伪逆矩阵。

只要测得包含 1阶固有频率的任-时间序列的响应信号，其他两种响应信号可通过对信号进行时域积分或微分变换来获得 。记录并输入所获得的位移、速度和加速度时间序列响应信号以及激振力信号，则可辨识出螺栓联接结构结合部的非线性特性参数。

由式(10)的辨识公式可知，准确获得各种响应信号是精确辨识螺栓联接结构结合部非线性特性参数的前提。而响应信号时域微积分过程产生的误差很容易积累在最终数据上，产生辨识误差。因此，在实际测试中，首先测得的通常是加速度信号，为了减小误差，试验中采用单-频率的正弦波激励，在实际积分之前首先对测得的加速度信号做如下处理星( )-∑[-A sin(n∞ )B cos(nwt)](11)即将信号可以表示为多次谐波的叠加，通过对时域加速度信号进行各次谐波的最小二乘拟合，可以获得其中的各次系数 A 和B 。这样，加速度信号就表示为解析方程的形式 ，然后对方程再进行解析积分就可以获得较精确的速度和位移信号，解析表达式为士-∑ [-A COS(Yt∞ ) -i-1 ，MDB sin(nwt)]C1 (12))- -志 [ i B cos(ncot)]C1tC2 (13)振动时，速度和位移信号的平均值为 0，所以式中C1-0，C20。

这种方法通过时域信号的解析表达式进行积分 ，避免了数值积分方法 的不稳定性和容易产生相位偏差的缺点 ，可以较为准确地获得所有 的时域 响应信号，减小了误差。

3.2 测试流程力状态映射法测试螺栓组合梁结合部特性参数的流程可概括为以下 6个关键步骤：(1)组合梁线性固有频率测试可采用经典的实验模态测试方法，例如，通过锤击 、电磁激振器、振动台等激励设备激励组合梁，用加速度传感器拾振完成模态测试▲-步获得结构的频响函数，从频响函数中获得组合梁的固有频率值。

(2)测试包含某阶固有频率的非线性频响测试非线性频响的目的是为了确认能够激发出结合部非线性特性的频带及激励幅值。当在某-激励幅度下，组合梁的结构的固有频率发生明显偏移，则说明组合梁结合部已表现出非线性振动特性，可记录这个频带和激励幅值。为了实现结合部特性参数的精确辨识，通出选取固有频率附件-小段频带(小于 5 Hz)，另外从测试的方法上可以采用振动台或激振器进行扫频测试 。

(3)定频定幅激励组合梁结构件在所确定的频带内选择若干激励频率点(通常大于 5个)，用所确定 的激振力幅定频激励试件 ，记录在每个频率下指定拾振点的响应信号和对应的激振力时域信号。

(4)时域响应信号的处理按 3.1所述的方法 ，对在各频率点下获得 的时域信号进行变换获得与加速度信号相对应的速度及位移响应信号。

(5)模态坐标变换所记录的激振力及位移响应信号均是用物理坐标描述的，而式(10)是用模态坐标系表示的辨识公式，为此需对时域响应信号及激振力信号进行模态坐标变换 。变换公式可描述为q( )- y(xp， )/‰ (xp)，口( )- (zP， )/9。(zP)，( )-3)(xp， )/9 (zP)，Q( )F(zF， )· 。(zF) (14)式中 y(x， )， (z， )，y(x， )为测试获得的拾振点P处位移、速度、加速度响应信号，F为激励点 F处的激振力信号。‰为正则化模态振型函数，可采取有限元法或解析法 ，获得组合梁的正则化函数值 。

(6)非线性参数识别利用 Force-state mapping方法，将某 阶共振频率附近各频率点模态坐标系下的时域响应信号、激振力信号代入到式(10)，可以求解出对应状态下组合梁结合部处的刚度及阻尼参数∩取平均值作为所辨识的结合部处立方刚度值，而螺栓联接结合部的阻尼通常于组合梁的振动幅值相关，即阻尼具有振幅依赖性。参照文[17]，结合部处的等效阻尼可表达为C -Co c1z c2z2 (15)进-步，可根据获得的等效阻尼系数和响应幅值，拟合出式(15)中的各个系数 。

4 研究实例螺栓联接梁组合结构见图 1，按照 3.2所述流第 2期 孙 伟，等：螺栓联接梁结构结合部非线性特性参数辨识程，利用力状态映射法辨识螺栓组合梁结合部力学特性参数。

4.1 实验测试首先，使用力矩扳手将试件固定在夹具上，夹具上螺栓的预紧力为 35 N·m，联接梁上螺栓的预紧力为 10 N·m，实验过程中保持不变。采用随机激励扫频测试获得螺栓联接梁的线性频响函数，如图4蜃频率 /Hz图4 螺栓联接梁的线性频响函数Fig.4 Linear FRF of the bolted beam structure进行实验模态分析，可以获得联接梁的前四阶线性固有频率见表 1。

表 1 螺栓联接梁的线性固有频率Tab.1 Linear nature frequencies of the bolted beam structure1阶 2阶 3阶 4阶实验模态分析/Hz 27.34 155.47 454.69 861.72有限元计算/Hz 27.30 137.8O 509.9O 861.80在 1阶固有频率 附近选取-小段频带 (26.9～27.5 Hz)，从高频到低频进行扫频测试，不断改变激振力幅，可以发现当激振力幅大于2 N时，螺栓联接结合部处出现非线性现象。随着激振力增大，非线性现象逐步明显(见图 2)。

分别选取 2和4 N作为激振力幅，在所选频段内 5个频率点进行定幅定频激励试件，记录时域加速度及激振力信号。按照式(I1)至(13)对时域加速度信号进行处理，分别变换成速度及位移响应信号。

4.2 参数辨识在进行力状态映射法辨识螺栓组合梁结合部特性参数前，按照测试流程，应对测试获得的物理坐标下的时域响应信号及激振力信号进行坐标变换。按耋 塑 立 篁辱避戚蕊 蕊 吾垂 圆蕊婴盔蕊 时间/图5 激振力 2 N、激振频率 26.9 Hz时的响应及激励信号Fig.5 Response and excitation signals at 2 N exciting照式(14)，应获得组合梁的线性模态振型函数。

这里采用有限元法获得螺栓组合梁的线性模态振形函数，具体方法为：建立螺栓组合梁的有限元模型，不断修改螺栓组合梁的结合部刚度系数，直到有限元计算的固有频率与实验模态测试的固有频率相- 致。最终有限元计算结果见表 1，此时，输出组合梁的经正则化处理的 1阶模态振型，如图 6。

根据离散点的位置，可从中读趣励处和拾振点处的正则化模态振型函数值，分别为‰ (xp)- 2.714 734q(xF)- 0.142 283利用式(14)完成坐标变换，进-步利用力状态映射法，输入模态坐标下的时域响应信号及激振力信号，可辨识出螺栓组合梁结合部处的特性参数。

表 2列出了在两种激振力幅值下，结构识别出的固有频率和立方刚度项系数的平均参数值。图 7为识别的不同频率点处对应这两个激励幅值的结合部等效阻尼系数。

X/m图 6 1阶模态振型图Fig.6 First-order modal shape19O 振 动 工 程 学 报 第 26卷表 2 螺栓联接 梁固有频 率及 立方 刚度项识别结果Tab.2 Identification results of nature frequencies and cubicstiffness for the bolted beam structure固有频率(c， ) 立方刚度( 。)F-2 N F-4 N F-2 N F-4 N27.26 27.24 1.002× 10 6.524× 10。吕邑 频率 /Hz图 7 不同频率对应两个激励幅值的结合部等效阻尼系数Fig.7 Equivalent damping coefficient of the joint for thedifferent frequency at the tWO exciting levels进-步，根据测得的等效阻尼系数和各频率点的响应幅值，可拟合出式(15)描述阻尼特性的二次方程系数，具体值见表3。实验数据及拟合曲线见图8和9。

表3 螺栓组合梁结合部阻尼系数识别结果Tab.3 Identification results of damping coefficient for thejoint of the bolted beam structure。

、。吕邑 ‰ ，m ×10图 8 2 N时阻尼方程系数的拟合Fig.8 Coefficients fitting of the damping equation at 2 Nexciting level5 结 论本文研究采用力状态映射法识别螺栓组合梁非蜀 1.0 1.5 2.0 2.5 3.O 3.5 4.0 4.5 5.0 5.5g衄 /m X 10-图 9 4 N时阻尼方程系数的拟合Fig.9 Coefficients fitting of the damping equation at 4 Nexciting level线性结合部特性参数的方法，结论包括：(1)对所研究的螺栓组合梁表征测试，发现组合梁表现出软式非线性特性，因而可用三次立方刚度描述结构的非线性特征。

(2)采用简化的 1wan模型可以很好的模拟联接单元的接触面在外载荷作用下产生的刚度非线性行为。采用连续梁理论，并引入 1wan模型所描述的结合部非线性行为，可建立起螺栓联接组合梁非线性解析模型。

(3)对应解析模型，采用力状态映射法，具体包括：组合梁线性固有频率测试、测试包含某阶固有频率的非线性频响、定频定幅激励组合梁结构件、时域响应信号的处理、模态坐标变换、非线性参数识别等6个关键步骤，可有效辨识螺栓组合梁结合部的刚度、阻尼等动力学特性参数。

(4)辨识结果显示在同样的螺栓预紧力下，随着激振力的增加螺栓联接结构的整体阻尼会增加，表明能量耗散增加；立方刚度项系数会减小，表明结构的整体刚度减校