CN103076194A - 实时混合模拟试验效果的频域评价方法 - Google Patents
实时混合模拟试验效果的频域评价方法 Download PDFInfo
- Publication number
- CN103076194A CN103076194A CN2012105931903A CN201210593190A CN103076194A CN 103076194 A CN103076194 A CN 103076194A CN 2012105931903 A CN2012105931903 A CN 2012105931903A CN 201210593190 A CN201210593190 A CN 201210593190A CN 103076194 A CN103076194 A CN 103076194A
- Authority
- CN
- China
- Prior art keywords
- frequency
- input signal
- real
- output signal
- time
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种实时混合模拟试验效果的频域评价方法,对实时混合模拟的输入和输出进行傅里叶变化,用变换后输入的频域结果除以变换后输出的频域结果,将结果按照特定的加权方法加权,可以得到评价指标FEI,计算评价指标FEI的幅值A和相位φ,若值A越接近于1和相位φ越接近于0,说明试验的效果越好。同时,采用与计算评价指标FEI相同的加权方法计算输入信号的等效频率feq,用相位φ除以等效频率feq可以得到实时混合模拟的时滞。该方法可以通过选取计算位移或者指令位移作为输入,测量位移作为输出,对实时混合模拟的效果进行评价,或者对不同的补偿方法进行评价。
Description
技术领域
本发明为一种实时混合模拟试验基于频域分析的评价方法,用于计算实时混合模拟的时滞并评价实时混合模拟的试验效果,同时可以比较实时混合模拟的不同算法和补偿方法的优劣。
背景技术
传统的土木工程结构试验方法包括拟静力试验、拟动力试验及地震模拟振动台试验等。拟静力试验是按照一定的荷载或变形控制方式,对试件进行反复的低周加载,使试件逐渐从弹性阶段过渡至破坏,其优点是经济、实用,但不能全面地反映结构的动力特性。拟动力试验可以将计算机的计算控制与结构试验结合在一起,同时在拟动力试验中可以采用子结构技术以达到节约成本的目的。由于试验在较大的时间尺度内进行,可以很容易的得到结构在不同荷载下的反映,但缺点是大时间尺度内进行试验不仅不能反映速度对结构的影响,同时还会导致应变速率效应。地震模拟振动台试验是目前最准确的结构抗震性能试验方法,但由于振动台的高昂造价和承载能力的限制,在进行大型结构试验时往往需要对模型进行缩尺,尺寸效应会对试验结果造成不同程度的影响,同时振动台只能处理地面振动对结构的作用,而无法反映其他类型,如风荷载的影响。
实时混合模拟是在拟动力试验的基础上发展起来的,该试验方法采用作动器实时加载,从而在保留拟动力试验优点的同时,克服拟动力试验不能反映速度对结构的影响的问题。积分算法和时滞补偿是沟通数值子结构和试验子结构的两个重要环节,在实时混合模拟过程中,积分算法在前一步的恢复力和外力的基础上求得作动器下一步的计算位移。由于***误差,作动器总是不能实时地达到预定的位移,而是存在一个“时滞”,该时滞相当于对结构增加一个负阻尼,如果不能进行合理的补偿,会对试验的准确性有很大的影响,时滞补偿的目的是对计算位移做某种修正得到指令位移,指令位移经过作动器的时滞延迟之后得到的测量位移与计算位移相同。
目前对实时混合模拟的评价方法都是在时域上进行的,时域评价指标具有计算简单、物理意义明确的优点,但存在只能从某一个方面(如位移误差或能量误差)来描述试验的效果的缺点,很难反映输入与输出的时间误差,从而不能全面地反映试验效果。同时,不同的算法和补偿方法的优劣性也难以比较。本发明即针对上述问题展开。
发明内容
技术问题:为了克服传统的时域评价方法在评价实时混合模拟效果的不足,本发明从频域的角度出发,提出了一种实时混合模拟试验效果的频域评价方法,可以对实时混合模拟的试验误差和时滞进行准确的评价,从而指导后续的试验数据处理工作。
技术方案:本发明解决其技术问题所采用的技术方案包含以下步骤:
第一步:对实时混合模拟的输入信号和输出信号分别乘以一个窗函数,得到加窗后的输入信号和输出信号;
第二步:对第一步加窗后的输入信号和输出信号进行快速傅立叶变换,得到输入信号和输出信号的频谱;
第三步:将第二步得到的输入信号的频谱对应于每个频率的值,分别除以输出信号的频谱在相同频率下对应的值,得到各频率下输入信号与输出信号的频谱比值;
第四步:将第二步得到的输入信号的频谱对应于每个频率的值取模,进一步得到模的n次幂,并对各频率下模的n次幂求和,得到各频率下模的n次幂之和;
第五步:将第四步得到的各频率下模的n次幂,除以各频率下模的n次幂之和,得到输入信号中该频率的权重;
第六步:在每一个频率下,将第三步得到的输入信号与输出信号的频谱比值,分别乘以对应频率下第五步得到的权重,然后将各个频率下的乘积结果进行累加,得到频域分析指标FEI;
第七步:计算第六步得到的频域分析指标FEI的幅值A和相位φ,幅值A用来判断试验的幅值误差大小,而相位φ用于计算时滞;
第八步:对第二步得到的输入信号的频谱,将输入信号的各个频率分别乘以其在第五步中的权重,并进行加权求和,得到输入的等效频率feq;
第九步:采用第七步得到的相位φ除以第八步得到的到输入信号的等效频率feq,得到实时混合模拟的时滞d。
其中,第七步得到的幅值A、相位φ和第九步得到的时滞d即为评价实时混合模拟的指标,若幅值A越接近于1、相位φ和时滞d越接近于0,说明试验的效果越好。
有益效果:采用本发明的“实时混合模拟试验效果的频域评价方法”后,能够准确的计算实时混合模拟的时滞和试验误差。对数据加窗的方法,可以消除信号在傅里叶变换过程中的频谱泄露现象,从而提高分析的准确性。在实时混合模拟的分析过程中,若将计算位移作为输入信号,测量位移作为输出信号,可以得到实时混合模拟的时滞,从而可以对补偿方法进行评价;若将指令位移作为输入信号,测量位移作为输出信号,可以得到实时混合模拟的试验误差和试验效果。因此,通过实施混合模拟频域评价方法,可以对试验效果进行评价的同时,还能比较不同补偿方法的优劣。
具体实施方式
第一步:对实时混合模拟的输入信号和输出信号加窗,减小频谱泄露对频域分析的影响。对于窗函数的选取,以汉宁窗(Hanning Window)为例,该步骤可以写为:
I(t)=I(t)′·H(L)(1,a)
O(t)=O(t)′·H(L)(1,b)
其中I(t)’和O(t)’分别代表输入信号和输出信号,I(t)和O(t)分别代表加窗输入信号和输出信号,H(L)代表长度为L的汉宁窗;
第二步:对第一步加窗后的输入信号和输出信号进行快速傅立叶变换,得到输入信号和输出信号的频谱,该频谱包含两个变量:一个是频率,一个是对应频率下频谱的值;
第三步:将第二步得到的输入信号的频谱对应于每个频率的值,分别除以输出信号的频谱在相同频率下对应的值,得到各频率下输入信号与输出信号的频谱比值;
第四步:将第二步得到的输入信号的频谱对应于每个频率的值取模,进一步得到模的n次幂,并对各频率下模的n次幂求和,得到各频率下模的n次幂之和;
第五步:将第四步得到的各频率下模的n次幂,除以各频率下模的n次幂之和,得到输入信号中该频率的权重;
第六步:在每一个频率下,将第三步得到的输入信号与输出信号的频谱比值,乘以第五步得到的权重,再将所有频率下的乘积进行累加,得到频域分析指标FEI:
其中,fft[I(t)]j和fft[O(t)]j分别表示对输入、输出信号做FFT变换所得到的对应于第j个频率的频域结果(含实部和虚部),k表示分析数据长度的一半;
第七步:计算第六步得到的频域分析指标FEI的幅值A和相位φ,幅值A用来判断试验的幅值误差大小,而相位φ用于计算时滞:
A=||FEI||(2,a)
φ=arctan[i(FEI)/r(FEI)](3,b)
其中,||□||表示求模。i()和r()分别表示虚部和实部,在相同情况下幅值A越接近于1和相位φ越接近于0,说明试验的效果越好;
第八步:对第二步得到的输入信号的频谱,将输入信号的各个频率分别乘以其在第五步中的权重,并进行加权求和,得到输入的等效频率feq:
其中,n的取值必须与第四步中的n相同,从计算得到的等效频率要接近线性结构的自振频率角度考虑,建议n的值取2;
第九步:采用第七步得到的相位φ除以第八步得到的等效频率feq,可以得到实时混合模拟的时滞d:
d=φ/(2π·feq)(5)
若d越接近0,试验的效果越好,当d大于0,说明补偿不足,当d小于0,说明补偿出现了过补偿。
Claims (1)
1.一种实时混合模拟试验效果的频域评价方法,其特征在于该方法包含以下步骤:
第一步:对实时混合模拟的输入信号和输出信号分别乘以一个窗函数,得到加窗后的输入信号和输出信号;
第二步:对第一步加窗后的输入信号和输出信号进行快速傅立叶变换,得到输入信号和输出信号的频谱;
第三步:将第二步得到的输入信号的频谱对应于每个频率的值,分别除以输出信号的频谱在相同频率下对应的值,得到各频率下输入信号与输出信号的频谱比值;
第四步:将第二步得到的输入信号的频谱对应于每个频率的值取模,进一步得到模的n次幂,并对各频率下模的n次幂求和,得到各频率下模的n次幂之和;
第五步:将第四步得到的各频率下模的n次幂,除以各频率下模的n次幂之和,得到输入信号中该频率的权重;
第六步:在每一个频率下,将第三步得到的输入信号与输出信号的频谱比值,分别乘以对应频率下第五步得到的权重,然后将各个频率下的乘积结果进行累加,得到频域分析指标FEI;
第七步:计算第六步得到的频域分析指标FEI的幅值A和相位φ,幅值A用来判断试验的幅值误差大小,而相位φ用于计算时滞;
第八步:对第二步得到的输入信号的频谱,将输入信号的各个频率分别乘以其在第五步中的权重,并进行加权求和,得到输入的等效频率feq;
第九步:采用第七步得到的相位φ除以第八步得到的到输入信号的等效频率feq,得到实时混合模拟的时滞d。
其中,第七步得到的幅值A、相位φ和第九步得到的时滞d即为评价实时混合模拟的指标,若幅值A越接近于1、相位φ和时滞d越接近于0,说明试验的效果越好。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210593190.3A CN103076194B (zh) | 2012-12-31 | 2012-12-31 | 实时混合模拟试验效果的频域评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210593190.3A CN103076194B (zh) | 2012-12-31 | 2012-12-31 | 实时混合模拟试验效果的频域评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103076194A true CN103076194A (zh) | 2013-05-01 |
CN103076194B CN103076194B (zh) | 2014-12-17 |
Family
ID=48152799
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210593190.3A Active CN103076194B (zh) | 2012-12-31 | 2012-12-31 | 实时混合模拟试验效果的频域评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103076194B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105335610A (zh) * | 2015-10-22 | 2016-02-17 | 东南大学 | 一种计算实时混合模拟试验瞬时时滞和幅值误差的方法 |
CN107356523A (zh) * | 2017-06-20 | 2017-11-17 | 东南大学 | 基于离散切线刚度估计的实时混合模拟试验反馈力修正方法 |
CN109827763A (zh) * | 2019-01-21 | 2019-05-31 | 东南大学 | 粘滞阻尼器检测效果评价方法 |
CN109883666A (zh) * | 2019-01-21 | 2019-06-14 | 东南大学 | 一种粘滞阻尼器检测位移输入信号修正方法 |
WO2021184753A1 (zh) * | 2020-03-16 | 2021-09-23 | 西安建筑科技大学 | 一种空间框架模型自适应实时混合仿真试验加载方法 |
US11150159B2 (en) | 2020-03-16 | 2021-10-19 | Xi'an University Of Architecture And Technology | Adaptive loading method for real-time hybrid simulation testing of space frame model |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001005854A (ja) * | 1999-06-24 | 2001-01-12 | Osaka Gas Co Ltd | 操作パネルシミュレータ及び操作パネルシミュレーション方法 |
CN1826632A (zh) * | 2003-03-31 | 2006-08-30 | 国际商业机器公司 | 用于语音信号的组合频域和时域音高提取的***和方法 |
CN101813497A (zh) * | 2010-04-22 | 2010-08-25 | 南京大学 | 一种布里渊散射谱实时频谱分析装置及其数据处理方法 |
CN102121991A (zh) * | 2010-01-08 | 2011-07-13 | 郑州威科姆科技股份有限公司 | 基于北斗一代卫星信号接收的干扰抑制方法及装置 |
-
2012
- 2012-12-31 CN CN201210593190.3A patent/CN103076194B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001005854A (ja) * | 1999-06-24 | 2001-01-12 | Osaka Gas Co Ltd | 操作パネルシミュレータ及び操作パネルシミュレーション方法 |
CN1826632A (zh) * | 2003-03-31 | 2006-08-30 | 国际商业机器公司 | 用于语音信号的组合频域和时域音高提取的***和方法 |
CN100589178C (zh) * | 2003-03-31 | 2010-02-10 | 国际商业机器公司 | 用于语音信号的组合频域和时域音高提取的***和方法 |
CN102121991A (zh) * | 2010-01-08 | 2011-07-13 | 郑州威科姆科技股份有限公司 | 基于北斗一代卫星信号接收的干扰抑制方法及装置 |
CN101813497A (zh) * | 2010-04-22 | 2010-08-25 | 南京大学 | 一种布里渊散射谱实时频谱分析装置及其数据处理方法 |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105335610A (zh) * | 2015-10-22 | 2016-02-17 | 东南大学 | 一种计算实时混合模拟试验瞬时时滞和幅值误差的方法 |
CN105335610B (zh) * | 2015-10-22 | 2017-12-19 | 东南大学 | 一种计算实时混合模拟试验瞬时时滞和幅值误差的方法 |
CN107356523A (zh) * | 2017-06-20 | 2017-11-17 | 东南大学 | 基于离散切线刚度估计的实时混合模拟试验反馈力修正方法 |
CN109827763A (zh) * | 2019-01-21 | 2019-05-31 | 东南大学 | 粘滞阻尼器检测效果评价方法 |
CN109883666A (zh) * | 2019-01-21 | 2019-06-14 | 东南大学 | 一种粘滞阻尼器检测位移输入信号修正方法 |
CN109883666B (zh) * | 2019-01-21 | 2020-09-11 | 东南大学 | 一种粘滞阻尼器检测位移输入信号修正方法 |
WO2021184753A1 (zh) * | 2020-03-16 | 2021-09-23 | 西安建筑科技大学 | 一种空间框架模型自适应实时混合仿真试验加载方法 |
US11150159B2 (en) | 2020-03-16 | 2021-10-19 | Xi'an University Of Architecture And Technology | Adaptive loading method for real-time hybrid simulation testing of space frame model |
Also Published As
Publication number | Publication date |
---|---|
CN103076194B (zh) | 2014-12-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103076194B (zh) | 实时混合模拟试验效果的频域评价方法 | |
CN102129520B (zh) | 动载荷的自适应时域识别方法 | |
CN101697084B (zh) | 一种基于rls滤波器的电液伺服***随机振动控制方法 | |
CN102122322B (zh) | 动载荷的自适应时域识别方法 | |
CN105353789B (zh) | 连续振动信号时间历程复现控制方法 | |
CN101216512A (zh) | 一种非正弦周期信号实时高精度检测方法 | |
Mao et al. | The construction and comparison of damage detection index based on the nonlinear output frequency response function and experimental analysis | |
CN104792450A (zh) | 一种设备对安装基座激励力的间接测量方法 | |
CN109299539A (zh) | 一种利用曲线转角指标评估悬臂梁刚度薄弱环节的方法 | |
CN107356523B (zh) | 基于离散切线刚度估计的实时混合模拟试验反馈力修正方法 | |
Islam et al. | DATA ANALYSIS METHODOLOGIES FOR HYDRODYNAMIC EXPERIMENTS IN WAVES. | |
CN106324340A (zh) | 一种同步相量和频率测量动态性能的方法 | |
CN105335610A (zh) | 一种计算实时混合模拟试验瞬时时滞和幅值误差的方法 | |
CN104462862A (zh) | 一种基于三次b样条尺度函数的机械结构动载荷识别方法 | |
CN101702092A (zh) | 一种基于卡尔曼滤波器的电液伺服***随机振动控制方法 | |
CN101917163B (zh) | 一种改善非正弦周期信号电液激振控制波形的方法 | |
CN109446552B (zh) | 多轴相关随机激励下结构疲劳寿命时域计算方法 | |
CN103399489B (zh) | 一种电液加速度伺服***的波形再现控制方法 | |
CN104834271B (zh) | 基于动刚度评价的直接进给轴伺服参数优化方法 | |
CN104931777A (zh) | 一种基于两条dft复数谱线的信号频率测量方法 | |
Damadipour et al. | Seismic modal identification using a new approach based on weighted transmissibility | |
CN115728513A (zh) | 一种新型冲击加速度测量装置及其测量方法 | |
CN103398832B (zh) | 一种电液加速度伺服***波形再现控制方法 | |
CN102567561B (zh) | 一种离散化板壳结构的建模方法 | |
CN102788662A (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 | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20130501 Assignee: ARCHITECTS & ENGINEERS CO., LTD. OF SOUTHEAST University Assignor: SOUTHEAST University Contract record no.: X2020980010481 Denomination of invention: Frequency domain evaluation method for real time hybrid simulation test Granted publication date: 20141217 License type: Common License Record date: 20201230 |