CN104361206A - 一种结合面动态参数的测量方法 - Google Patents
一种结合面动态参数的测量方法 Download PDFInfo
- Publication number
- CN104361206A CN104361206A CN201410568646.XA CN201410568646A CN104361206A CN 104361206 A CN104361206 A CN 104361206A CN 201410568646 A CN201410568646 A CN 201410568646A CN 104361206 A CN104361206 A CN 104361206A
- Authority
- CN
- China
- Prior art keywords
- prime
- lambda
- theta
- minor structure
- transport function
- 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.)
- Pending
Links
Abstract
本发明公开了一种结合面动态参数的测量方法,提出对于由第一子结构、第三子结构经结合面第二子结构联接形成的整体联接结构,通过测试整体联接结构和第一子结构、第三子结构的传递函数,能分析得到中间第二子结构(结合面)的传递函数,从而形成了一种方便的结合面动态参数的测量方法。该结合面动态参数的测量方法建立了结合面子结构与整体联接结构和非结合面子结构的传递函数的关系式,并采用横向位移传递函数替换带转角的传递函数的方法,从而可通过测试横向平动位移传递函数,分析得到结合面动态参数。本发明用于结合面动态参数的测量,方法简单,操作方便,并解决了带转角的传递函数难测的问题。
Description
技术领域
本发明属于结构动力学参数测试技术领域,尤其涉及一种结构动力学建模中结合面动态参数的测量方法。
背景技术
参数识别问题是结构动力学建模中的一个极其重要的问题,许多研究者对于结合面的参数识别进行了大量研究。
但是,目前的一些方法通常只考虑结构纵向振动的形式,即结合面只有一个方向的自由度,没有考虑带转角自由度参数的识别。
发明内容
本发明实施例的目的在于提供一种结构动力学建模中结合面动态参数的测量方法,旨在解决目前的参数识别方法通常只考虑结构纵向振动的形式,即结合面只有一个方向的自由度,没有考虑带转角自由度参数的识别的问题。
本发明实施例是这样实现的,一种结合面动态参数的测量方法,该结构动力学建模中结合面动态参数的测量方法包括:通过测试整体联接结构A和第一子结构B、第三子结构D的传递函数,分析得到中间第二子结构C(结合面)的传递函数;建立结合面子结构与整体联接结构和非结合面子结构的传递函数的关系式;并采用横向位移传递函数替换带转角的传递函数。
进一步,建立结合面子结构与整体联接结构和非结合面子结构的传递函数的关系式的步骤如下:
建立梁组合结构的模型,系统可视为由第一子结构B、第二子结构C和第三子结构D通过四个坐标V1(挠度)、θ1(转角)、V2(挠度)、θ2(转角)联结而成的组合系统来研究,其中F和M分别代表图中的力和力偶矩:
对于第一子结构B:
对于第二子结构C:
对于第三子结构D:
约束方程,对于截面I:
对于截面I:
根据上述运动方程及约束方程,整体联接结构A的传递函数推导如下:
对于3个子结构,有:
根据约束方程得:
则式(4-11)、(4-12)变为:
将式(4-13)、(4-14)扩展并合并,得:
由上述推导得到,整体联结结构的传递函数矩阵为:
同理,由式(4-16)可知:
则
根据上述推导公式,通过测试整体联接结构A和第一子结构B、D的传递函数,能得到中间子结构的传递函数。
进一步,采用横向位移传递函数替换带转角的传递函数包括以下步骤:
第一子结构B传递函数为 则 则传递函数矩阵为:
第三子结构D传递函数为 则 则传递函数矩阵为:
中间第二子结构C,传递函数
f6=cosλl·shλl+sinλl·chλl,f7=sinλl+shλl,f8=sinλl-shλl;
f9=cosλl+chλl,f10=cosλl-chλl;
根据式(4-17),得:
式中,
而中间第二子结构C的传递函数为:
通过采用不含转角变量的传递函数来代换含转角变量的传递函数。
本发明提供的结构动力学建模中结合面动态参数的测量方法,提出对于由第一子结构B、第三子结构D经结合面(中间第二子结构C)联接形成的整体联接结构A,通过测试整体联接结构A和第一子结构B、第三子结构D的传递函数,能分析得到中间第二子结构C(结合面)的传递函数,从而形成了一种方便的结合面动态参数的测量方法。该结合面动态参数的测量方法建立了结合面子结构与整体联接结构和非结合面子结构的传递函数的关系式,并采用横向位移传递函数替换带转角的传递函数的方法,从而可通过测试横向平动位移传递函数,分析得到结合面动态参数。本发明用于结合面动态参数的测量,方法简单,操作方便,并解决了带转角的传递函数难测的问题。
附图说明
图1是本发明实施例提供的结构动力学建模中结合面动态参数的测量方法流程图;
图2是本发明实施例提供的梁组合结构联结图;
图3是本发明实施例提供的采用MATLAB语言编程画出传递函数的示意图;
图4是本发明实施例提供的得到整体结构A在联结界面处的传递函数值示意图;
图5是本发明实施例提供的传递函数幅值图;
图6是本发明实施例提供的传递函数实部图;
图7是本发明实施例提供的传递函数虚部图;
图8是本发明实施例提供的采用MATLAB画出横坐标在200000与2000000rad/s之间的幅值函数示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
下面结合附图及具体实施例对本发明的应用原理作进一步描述。
本发明实施例的结构动力学建模中带转角的传递的测量方法包括以下步骤:
S101:制作第一子结构B、第三子结构D,进行模态试验,分别测试第一子结构B、第三子结构D的传递函数;
s102:进行模态试验,分别测试第一子结构B、第三子结构D的传递函数;
s103:第一子结构B、第三子结构D经结合面连接形成整体联接结构A;
S104:进行模态试验,测试整体联接结构A的传递函数;
s105:采用横向位移传递函数替换带转角的传递函数;
s106:计算中间第二子结构C的传递函数,分析得到结合面动态参数。
结合以下的实施例对本发明的原理做进一步的叙述:
1、本发明的梁组合结构传递函数矩阵的理论推导:
梁组合结构的模型如图2所示,该系统可视为由第一子结构B、第二子结构C和第三子结构D通过四个坐标V1(挠度)、θ1(转角)、V2(挠度)、θ2(转角)联结而成的组合系统来研究,其中F和M分别代表图中的力和力偶矩:
对于第一子结构B:
对于第二子结构C:
对于第三子结构D:
约束方程,对于截面I:
对于截面I:
根据上述运动方程及约束方程,整体联接结构A的传递函数推导如下:
对于3个子结构,有:
根据约束方程得:
则式(4-11)、(4-12)变为:
将式(4-13)、(4-14)扩展并合并,得:
由上述推导得到,整体联结结构的传递函数矩阵为
同理,由式(4-16)可知:
则
根据上述推导公式,通过测试整体联接结构A和第一子结构B、第三子结构D的传递函数,能得到中间第二子结构C的传递函数,由于试验条件的限制,带转角的传递函数无法测试,故考虑用不含转角变量的传递函数来代换含转角变量的传递函数。
2、由于带转角的传递函数在试验中很难测得,根据左鹤声《机械阻抗方法与应用》推导出的两端自由梁振动传递函数的表达式,由带转角的传递函数和横向位移传递函数之间的关系,将带转角的传递函数用横向位移传递函数来替换;
具体包括以下步骤:
设第一子结构B传递函数为 则 则传递函数矩阵为:
第三子结构D传递函数为 则 则传递函数矩阵为:
中间第二子结构C,传递函数
f5=cosλl·shλl-sinλl·chλl;
f6=cosλl·shλl+sinλl·chλl,f7=sinλl+shλl,f8=sinλl-shλl;
f9=cosλl+chλl,f10=cosλl-chλl;
根据式(4-17),得:
式中,
而中间子结构的传递函数为:
通过采用不含转角变量的传递函数来代换含转角变量的传递函数,解决了带转角的传递函数不能测量的问题。
3、下面构造传递函数,对于无阻尼梁和阻尼梁进行验算:
3.1、无阻尼梁验算:
设一均匀弹性梁A,梁长lA=0.6m,弹性模量E=5×104MPa,截面尺寸a=0.02m,b=0.03m,密度ρ=2600kg/m3,根据公式写出梁A作自由振动时的传递函数值为:
令上式分母为0,即:
cosλl·chλl-1=0
解此超越方程即可得到各阶固有频率,可得到:
λ0l=0,λ1l=4.730,λ2l=7.853,λ3l=10.996;
由式(4-28b)得到相应的各阶固有频率为
代入λnl数值得:
将该梁均分为B、C、D三个子结构,lB=lC=lD=0.2m,其他参数与整体梁A均一样,根据上述推导过程,由公式(4-28b)可得到整体结构A在联结界面处的传递函数值;
本发明采用MATLAB语言编程画出传递函数的图象,从图3和表1中可看出,极值点频率与前面采用整体梁计算的结果基本符合,说明推导公式的正确性;
表1组合结构频率计算值
再将该梁分为B、C、D三个不等的子结构,lB=lD=0.25m,lC=0.1m,其他参数与整体梁A均一样,根据上述推导过程,由公式(4-28b)可得到整体结构A在联结界面处的传递函数值,如表2和图4,画出该传递函数图象为图4;
表2组合结构频率计算值
由图4可看出,用不等长三个子结构表示的整体梁的固有频率计算结果与理论结果一致,说明了该推导公式对无阻尼梁适用。
3.2有阻尼梁组合结构传递函数矩阵的验算:
设一均匀弹性梁A,梁长lA=0.6m,弹性模量E=5×104MPa,截面尺寸a=0.02m,b=0.03m,密度ρ=2600kg/m3。将梁A分为三个第一子结构B、C、D,每个子结构为划分为10个单元。
采用有限元方法验证,根据梁单元:
根据传递函数模态表达式:
来构造传递函数,写出传递函数矩阵。将传递函数代入到推导公式中,画出的传递函数图像为图5-图7;
表3 组合结构频率阻尼比计算值
由表3可看出,用不等长三个子结构表示的有阻尼梁的固有频率计算结果与理论结果一致,说明了该推导公式对阻尼梁也适用;
3.3、中间子结构的传递函数识别验算:
由于整体梁带转角的传递函数同样无法通过试验测试得到,所以根据第一子结构B、C、D的组合方式得到的整体梁A两个截面的传递函数矩阵表达式,采用无转角传递函数代替带转角传递函数的方法,即由式(4-18)可算出矩阵各元素的代换系数:
根据代换系数值,在试验过程中,只需要测试整体梁A不带转角的传递函数一个值,根据代换系数即可得到矩阵中其他的元素。
验算过程如下:
设一均匀弹性梁A,由三个第一子结构B、C、D组成。梁纵长lA=0.6m,弹性模量E=5×104MPa,截面尺寸a=0.02m,b=0.03m,阻尼系数为c=100N·s/m,密度ρ=2600kg/m3。三个第一子结构B、C、D的长度依次为lB=0.28m,lC=0.04m,lA=0.28m,其他参数与梁A一致。
同样采用有限元方法构造整体梁A的和梁B、D的传递函数,根据式(4-27)得到识别出的梁C的传递函数采用MATLAB画出横坐标在200000与2000000rad/s之间的幅值函数图像如图8;
表4 组合结构频率计算值
4、结合面模态参数识别:
4.1 基于分量分析法的模态参数识别:
根据对试件的动态性能参数的测试,获得了固有频率、模态阻尼比、振型和传递函数值。根据传递函数的表达式:
式中,当ω趋近于某阶模态的固有频率时,该模态起主导作用,称为主模态;在主模态附近,其他模态影响较小;如果系统的阻尼较小且模态偶合较轻,则r阶模态以外的其余模态导纳可近似地用复常数HC称为剩余导纳,和分别为HC的实部和虚部;
Hlp(ω)的实部与虚部可表示如下:
实验证明,实频特性曲线的峰值频率即为系统的固有频率,而模态阻尼比为式中ωb-ωa为r阶模态导纳的半功率带宽;
对主模态而言,当时,分别测出结构上各点的值(L=1、2,…,L)则当第r阶模态的振型系数列阵:
对r阶模态,当采用单点激励时,为常数。因此即可代表模态振型。因为振型只是反应振动形态,与振动大小无关,故常取归一化振型列阵。
对p点激励,l点响应的传递函数,当ω=ωr时,第r阶的值可表示如下:
若取原点导纳元素且对原点归一化,令φpr=1,得各阶模态刚度为:
必须指出,模态刚度数值与归一化有关,不同归一化有不同的模态刚度。而在模态刚度求出后,也可求出模态质量mr:
由此全部模态参数均可求得。
4.2实测结合面试件的模态参数识别:
由测试得到的粗骨料及其相应结合面试件的传递函数,读出各频率下的传递函数值,根据无转角和带转角传递函数之间的关系,对无法测试的带转角的传递函数,采用带转角传递函数代换无转角传递函数的方法,得到带转角的传递函数值。由此可得到整体结构试件A传递函数矩阵HA和两端子结构试件B、D的传递函数矩阵HB和HD。因为结合面的厚度很小,响应频率很高(理论计算大概为105数量级),试验条件无法测试,故本发明采用估算中间第二子结构C的模态阻尼比的方式,根据前面推导的组合传递函数公式,得到由实测B、D和估算的C传递函数组合而成的整体结构A传递函数,再与实测A的传递函数比较来确定第二子结构C的传递函数;
例如对卵石2-卵石5(用卵石混凝土C40湿筛砂浆粘接)结合面试件,分别测试得到卵石2和卵石5粗骨料试件和联结以后结合面试件对应位置的传递函数图像,读出相应频率下对应的传递函数实部A和虚部B值,写成A+Bi的形式,得到卵石2-卵石5结合面试件的实测传递函数值,如表5所示。
表5卵石2-卵石5结合面试件传递函数实测值
频率(Hz) | 卵石2 | 卵石5 | 整体结构 |
195.31 | 9.66e-1+1.47i | 1.47+1.40i | -9.62+9.89i |
214.84 | 1.04+1.62i | 1.53+1.73i | -4.68+13.9i |
234.38 | 1.11+1.73i | 1.79+1.81i | -4.87+15.5i |
273.44 | 1.20+1.89i | 2.18+1.89i | -24.8+9.41i |
292.97 | 1.07+2.05i | 2.5+1.92i | -1.81-6.56i |
312.50 | 1.55+2.26i | 2.43+1.71i | -14.0-6.33 |
332.03 | 1.61+2.20i | 2.41+1.84i | -11.6-5.71i |
488.28 | 1.92+2.82i | 4.07+2.66i | -8.38+6.49i |
507.81 | 2.10+3.17i | 3.62+2.17i | -6.88+9.33i |
527.34 | 2.26+3.15i | 3.72+2.67i | -2.43+10.6i |
546.88 | 2.24+3.28i | 3.95+3.04i | 1.43+8.04i |
566.41 | 2.27+3.52i | 4.26+3.35i | -1.19e-23.89i |
585.94 | 2.21+3.81i | 4.69+3.70i | -2.67+3.75i |
634.77 | 2.84+4.57i | 6.41+4.19i | -7.88+5.07i |
722.66 | 4.81+7.17i | 12.1+7.39i | -13.0-28.5i |
810.55 | 16.6-11.3i | 2.63-23.6i | -11.6+1.10i |
908.20 | 1.09-1.40i | 2.32-5.22i | -13.6+9.50i |
976.56 | 4.99e-1+3.02e-1i | 2.62-3.65i | -18.3+22.2i |
1015.63 | -8.89e-2+1.08i | 3.55-2.57i | -22.4+38.1i |
1064.45 | -8.02e-1+1.86i | 2.8-2.26i | -11.9+49.3i |
1103.52 | -2.0+2.53i | 1.09-2.72i | 14.6+29.0i |
1123.05 | -2.51+3.14i | 8.81-2.62i | 10.4+13.0i |
对卵石2和卵石5传递函数值采用多项式拟合的方式,得到各自的传递函数表达式 构造中间子结构传递函数表达式根据转角代换关系,得到三个结构的传递函数矩阵。逐步尝试变化中间子结构的阻尼系数c,使组合得到的整体结构传递函数与实测传递函数大致相等,以此确定中间粘接结构的阻尼系数c。分别估算出每个结合面试件的阻尼系数值,列如下表6-表8。
表6 400×20×10mm结合面试件识别结果
表7 400×20×20mm结合面试件识别结果
表8 400×20×30mm结合面试件识别结果
从表6-表8中列出数据可看出,识别出的结合面的阻尼系数值有一定的规律性。结合面材料强度大的,阻尼系数偏小。由于试验的误差,同种材料试件识别出的阻尼系数有所不同,但基本上差别不大。识别结果验证了材料的阻尼比随着强度的增大有减小的趋势。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种结合面动态参数的测量方法,其特征在于,该结合面动态参数的测量方法包括:制作第一子结构、第三子结构,进行模态试验,分别测试第一子结构、第三子结构的传递函数;进行模态试验,分别测试第一子结构、第三子结构的传递函数;第一子结构、第三子结构经结合面连接形成整体联接结构;进行模态试验,测试整体联接结构的传递函数;采用横向位移传递函数替换带转角的传递函数;计算中间第二子结构的传递函数,分析得到结合面动态参数。
2.如权利要求1所述的结构动力学建模中结合面动态参数的测量方法,其特征在于,建立结合面子结构与整体联接结构和非结合面子结构的传递函数的关系式的步骤如下:
建立梁组合结构的模型,系统视为由第一子结构、第二子结构和第三子结构通过四个坐标V1(挠度)、θ1(转角)、V2(挠度)、θ2(转角)联结而成的组合系统来研究,其中F和M分别代表图中的力和力偶矩:
对于第一子结构:
对于第二子结构:
对于第三子结构:
约束方程,对于截面I:
对于截面I:
根据上述运动方程及约束方程,整体联接结构的传递函数推导如下:
对于3个子结构,有:
根据约束方程得:
则式(4-11)、(4-12)变为:
将式(4-13)、(4-14)扩展并合并,得:
由上述推导得到,整体联结结构的传递函数矩阵为:
同理,由式(4-16)可知:
则
根据上述推导公式,通过测试整体联接结构和第一子结构、第三子结构的传递函数,能得到中间第二子结构的传递函数。
3.如权利要求1或2所述的结构动力学建模中结合面动态参数的测量方法,其特征在于,采用横向位移传递函数替换带转角的传递函数包括以下步骤:
第一子结构传递函数为 则 则传递函数矩阵为:
第三子结构传递函数为 则 则传递函数矩阵为:
中间第二子结构,传递函数
其中,f1=sinλl·shλl,f2=cosλl·chλl,f3=cosλl·chλl-1,f4=cosλl·chλl+1;f5=cosλl·shλl-sinλl·chλl;
f6=cosλl·shλl+sinλl·chλl,f7=sinλl+shλl,f8=sinλl-shλl;
f9=cosλl+chλl,f10=cosλl-chλl;
根据式(4-17),得:
式中,
而中间第二子结构的传递函数为:
通过采用不含转角变量的传递函数来代换含转角变量的传递函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410568646.XA CN104361206A (zh) | 2014-10-23 | 2014-10-23 | 一种结合面动态参数的测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410568646.XA CN104361206A (zh) | 2014-10-23 | 2014-10-23 | 一种结合面动态参数的测量方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104361206A true CN104361206A (zh) | 2015-02-18 |
Family
ID=52528466
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410568646.XA Pending CN104361206A (zh) | 2014-10-23 | 2014-10-23 | 一种结合面动态参数的测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104361206A (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102095562A (zh) * | 2010-12-31 | 2011-06-15 | 北京工业大学 | 滑动结合面动态特性测试装置及测试方法 |
CN102779200A (zh) * | 2011-12-14 | 2012-11-14 | 西安交通大学 | 一种包含微观表面形貌的结合面接触性能分析方法 |
-
2014
- 2014-10-23 CN CN201410568646.XA patent/CN104361206A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102095562A (zh) * | 2010-12-31 | 2011-06-15 | 北京工业大学 | 滑动结合面动态特性测试装置及测试方法 |
CN102779200A (zh) * | 2011-12-14 | 2012-11-14 | 西安交通大学 | 一种包含微观表面形貌的结合面接触性能分析方法 |
Non-Patent Citations (1)
Title |
---|
危小灵: "结合面动力参数的识别及其对混凝土阻尼影响的研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Liu et al. | Vibration response of multi-span fluid-conveying pipe with multiple accessories under complex boundary conditions | |
Mu et al. | Decoupled schemes for a non-stationary mixed Stokes-Darcy model | |
Bayesteh et al. | XFEM fracture analysis of shells: the effect of crack tip enrichments | |
Li et al. | Bandgap calculations of two-dimensional solid–fluid phononic crystals with the boundary element method | |
Benmeddour et al. | Numerical modeling of guided wave interaction with non-axisymmetric cracks in elastic cylinders | |
CN103049670B (zh) | 管道激励源识别及其振动响应预测方法 | |
Mirsaidov et al. | Spatial stress state and dynamic characteristics of earth dams | |
Intarit et al. | Penny-shaped crack in elastic medium with surface energy effects | |
Hild et al. | Residual a posteriori error estimators for contact problems in elasticity | |
Yong-Hui et al. | Practical formula to calculate tension of vertical cable with hinged-fixed conditions based on vibration method | |
Benveniste | Models of thin interphases with variable moduli in plane-strain elasticity | |
Wang et al. | A POSTERIORI ERROR ESTIMATE FOR STABILIZED FINITE ELEMENT METHODS FOR THE STOKES EQUATIONS. | |
CN104361206A (zh) | 一种结合面动态参数的测量方法 | |
Xie et al. | New mixed finite elements for plane elasticity and Stokes equations | |
Hamza et al. | Pre-designing of a mechatronic system using an analytical approach with dymola | |
Bermúdez et al. | Finite element analysis of sloshing and hydroelastic vibrations under gravity | |
Beretta et al. | An asymptotic formula for the displacement field in the presence of thin elastic inhomogeneities | |
JP4182111B2 (ja) | 繊維集合体の音響特性予測方法 | |
Semenyuk et al. | Stability and initial postbuckling behavior of anisotropic cylindrical shells under external pressure | |
Winter et al. | Frequency dependent modeling for the prediction of the sound transmission in timber constructions | |
Noguchi et al. | Topological derivative for an acoustic-elastic coupled system based on two-phase material model | |
Gatica et al. | Numerical analysis of a transmission problem with Signorini contact using mixed-FEM and BEM | |
Bogolyubov et al. | Mathematical simulation of an irregular waveguide with reentering edges | |
Pu et al. | The hydrostatic approximation of the Boussinesq equations with rotation in a thin domain | |
CN116754091B (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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150218 |