CN105335610B - 一种计算实时混合模拟试验瞬时时滞和幅值误差的方法 - Google Patents
一种计算实时混合模拟试验瞬时时滞和幅值误差的方法 Download PDFInfo
- Publication number
- CN105335610B CN105335610B CN201510690506.4A CN201510690506A CN105335610B CN 105335610 B CN105335610 B CN 105335610B CN 201510690506 A CN201510690506 A CN 201510690506A CN 105335610 B CN105335610 B CN 105335610B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- window function
- displacement
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种计算实时混合模拟试验瞬时时滞和幅值误差的方法,通过预估结构周期和非线性程度,选择合适长度的窗函数;对实时混合模拟的计算位移和测量位移按照选定的窗函数进行窗变换,计算该窗函数上的时滞和幅值误差作为窗函数末时刻的瞬时时滞和幅值误差;将窗函数随时间移动,计算不同时刻的瞬时时滞和幅值误差,直至实验结束。本发明方法可以对实时混合模拟的实验效果进行在线评估,避免由于不合理补偿导致的实验构件破坏,并为补偿方法提供必要的信息参考。
Description
技术领域
本发明为一种能够计算实时混合模拟试验的瞬时时滞和幅值误差的方法,用于试验中在线评估实时混合模拟的实验效果,避免由于不合理补偿导致的实验构件破坏以及可能出现的人身伤害,并为实时混合模拟试验的补偿方法提供参考。
背景技术
实时混合模拟实验方法将混合模拟实验和振动台实验的优点结合在一起,能够在大尺度或者足尺的条件下对实验结构进行实时加载,从而克服有效的克服振动台试验尺寸效应和拟动力试验无法反映速度对结构的影响的问题。该实验方法,该方法将结构分为试验子结构和数值子结构两部分,其中,实验子结构在实验室中进行实时实验,而数值子结构在计算机中进行数值模拟,子结构之间通过伺服作动器进行信息交换。
由于伺服系统固有的特性,作动器到达指令位移会有一个延迟,该延迟称为时滞。同时,由于实验中的误差,作动器的测量位移与计算位移之间在幅值上同样存在误差,该误差称为幅值误差。时滞和幅值误差都会对实验结果造成很大影响。
时滞和幅值误差产生的原因复杂,且受到的影响因素较多,因此在实验过程中时滞和幅值误差是随时间变化而变化的。而目前对时滞和幅值误差的计算往往从整体的角度出发,计算整个实验的平均时滞和平均幅值误差,与实际情况不符。在实际的实验中,结构的最大响应往往持续时间非常短,约2s到8s之间,而整个实验的持时可以长达30s甚至更多,因此瞬时时滞和幅值误差的计算是非常有必要的。另外,从安全和经济角度考虑,在实验中应避免由于不合理补偿导致的实验构件破坏,也需要对瞬时时滞和幅值误差进行实时监控。本发明即针对上述问题展开。
发明内容
本发明的目的是为了克服传统评价方法无法计算实时混合模拟瞬时时滞和幅值误差的不足,提供一种计算实时混合模拟试验瞬时时滞和幅值误差的方法,采用移动窗函数的方法,计算实时混合模拟的瞬时时滞和幅值误差,可以对实时混合模拟的实验效果进行在线评估,避免由于不合理补偿导致的实验构件破坏以及可能出现的人身伤害,并为补偿方法提供必要的信息参考。
本发明采用的技术方案为:一种计算实时混合模拟试验瞬时时滞和幅值误差的方法,包括以下步骤:
第一步:估算结构的周期,并预测结构在实验中是否会出现非线性以及非线性程度;
第二步:确定窗函数的长度。对于线性结构,窗函数长度取结构周期的2倍;对于非线性结构,根据非线性的程度,窗函数长度取结构进入非线性结构前周期的2到4倍;
第三步:实验开始时,将窗函数的起始时刻分别与计算位移和测量位移的初始时刻重合。当实验中计算位移和测量位移数据长度与窗函数长度相同时,对计算位移和测量位移进行窗变换并计算该窗函数上的时滞和幅值误差,该时滞和幅值误差可以认为是窗函数末时刻的瞬时时滞和幅值误差。
第四步:将窗函数向后移动,移动的时间长度等于测量位移(或者计算位移)相邻两点的时间。按照第三步的方法,计算此窗函数上的时滞和幅值误差,作为窗函数末时刻的瞬时时滞和幅值误差;
第五步:继续移动窗函数,计算窗函数的时滞和幅值误差,直至实验结束。
作为优选,所述第一步中结构的周期可以通过能量法、等效质量法或顶点位移法进行估算。对于比较规则的N层结构,也可以采用以下近似公式计算:
框架结构:T=(0.08~0.10)N (1,a)
框剪结构:T=(0.06~0.08)N (1,b)
剪力墙结构:T=(0.06~0.08)N (1,c)
对于结构的非线性以及非线性程度,可以通过比较结构可能的最大位移与屈服位移的大小得到。
作为优选,所述第三步中采用实时混合模拟频域评价指标计算该窗函数上的时滞和幅值误差:
A=‖FEI‖ (2,c)
φ=arc tan[Im(FEI)/Re(FEI)] (2,d)
d=φ/(2π·feq) (2,e)
其中yI(f)和yo(f)分别代表计算位移和测量位移加窗后的傅里叶变换,p代表窗函数长度的一半,Im(·)和Re(·)分别代表FEI的虚部和实部,FEI代表频域评价指标,feq表示等效频率,A和d分别代表实时混合模拟的的幅值和等效时滞,φ表示实时混合模拟的相位差。
本发明的有益效果:(1)可以计算出实时混合模拟实验任意时刻的时滞和幅值误差,从而能够对实时混合模拟实验进行在线评估;(2)当实验中瞬时时滞和幅值误差过大时停止实验,避免由于不合理补偿导致的实验构件破坏以及可能出现的人身伤害;(3)瞬时时滞和幅值误差可以为补偿方法提供必要的依据,并对补偿方法进行修正,从而提高补偿效果。
附图说明
图1为本发明方法的示意图。
图中1:实时混合模拟测量位移或计算位移;2:第一个窗函数;3第二个窗函数;4;时间轴。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步的说明。
一种计算实时混合模拟试验瞬时时滞和幅值误差的方法,包括以下步骤:
第一步:估算结构的周期,并预测结构在实验中是否会出现非线性以及非线性程度。结构的周期可以通过能量法、等效质量法和顶点位移法等方法进行估算。对于比较规则的N层结构,也可以采用以下近似公式计算:
框架结构:T=(0.08~0.10)N (1,a)
框剪结构:T=(0.06~0.08)N (1,b)
剪力墙结构:T=(0.06~0.08)N (1,c)
对于结构的非线性以及非线性程度,可以通过比较实时混合模拟测量位移或计算位移1的最大值与屈服位移的大小得到。
第二步:确定窗函数的长度。窗函数的长度选择是计算瞬时时滞的关键问题,当窗函数的长度较大时,计算得到的幅值误差和时滞会更加准确。但由于移动窗函数采用整个窗函数的平均幅值误差和时滞对实验进行在线评估,过大的窗函数同样会降低计算结果的时域分辨率,导致计算结果存在偏差。
因此,为了同时保证计算结果的准确性和足够的时域分辨率,一般选择结构周期的2倍作为窗函数的长度。然而,当结构进入非线性后刚度下降,从而导致结构周期变大。非线性程度越大,结构周期变化越大。因此,对于非线性结构,窗函数的长度根据非线性程度窗函数长度取结构进入非线性结构前周期的2到4倍。
第三步:实验开始时,将窗函数的起始时刻分别与计算位移和测量位移的初始时刻重合。当计算位移和测量位移数据长度与窗函数长度相同时,对计算位移和测量位移施加第一个窗函数2,并采用实时混合模拟频域评价指标计算该窗函数上的时滞和幅值误差,该时滞和幅值误差可以认为是窗函数末时刻的瞬时时滞和幅值误差。
实时混合模拟频域评价指标计算该窗函数上的时滞和幅值误差方法如下:
A=‖FEI‖ (2,c)
φ=arc tan[Im(FEI)/Re(FEI)] (2,d)
d=φ/(2π·feq) (2,e)
其中yI(f)和yo(f)分别代表计算位移和测量位移加窗后的傅里叶变换,p代表窗函数长度的一半,Im(·)和Re(·)分别代表FEI的虚部和实部,FEI代表频域评价指标,feq表示等效频率,A和d分别代表实时混合模拟的的幅值和等效时滞,φ表示实时混合模拟的相位差。
第四步:将窗函数向后移动,移动的时间长度等于测量位移(或者计算位移)相邻两点的时间,施加第二个窗函数3。测量位移(或者计算位移)相邻两点的时间是指作动器可以记录的位移中,距离最短的两个位移之间的时间差,一般去1/1024s。按照第三步的方法,计算此窗函数上的时滞和幅值误差,作为末时刻的瞬时时滞和幅值误差;
第五步:继续在时间轴4上移动窗函数,计算窗函数的时滞和幅值误差,直至实验结束。
以上结合附图对本发明的实施方式做出详细说明,但本发明不局限于所描述的实施方式。对本领域的普通技术人员而言,在本发明的原理和技术思想的范围内,对这些实施方式进行实施方式进行多种变化、修改、替换和变形仍落入本发明的保护范围内。
Claims (1)
1.一种计算实时混合模拟试验瞬时时滞和幅值误差的方法,其特征在于:包括以下步骤:
第一步:估算结构的周期,并预测结构在实验中是否会出现非线性以及非线性程度;
第二步:确定窗函数的长度;对于线性结构,窗函数长度取结构周期的2倍;对于非线性结构,根据非线性的程度,窗函数长度取结构进入非线性结构前周期的2到4倍;
第三步:实验开始时,将窗函数的起始时刻分别与计算位移和测量位移的初始时刻重合;当实验中计算位移和测量位移数据长度与窗函数长度相同时,对计算位移和测量位移进行窗变换并计算该窗函数上的时滞和幅值误差,该时滞和幅值误差认为是窗函数末时刻的瞬时时滞和幅值误差;
第四步:将窗函数向后移动,移动的时间长度等于测量位移、或者计算位移相邻两点的时间;按照第三步的方法,计算此窗函数上的时滞和幅值误差,作为窗函数末时刻的瞬时时滞和幅值误差;
第五步:继续移动窗函数,计算窗函数的时滞和幅值误差,直至实验结束;
所述第一步中结构的周期通过能量法、等效质量法或顶点位移法进行估算;对于比较规则的N层结构,采用以下近似公式计算:
框架结构:T=(0.08~0.10)N (1,a)
框剪结构:T=(0.06~0.08)N (1,b)
剪力墙结构:T=(0.06~0.08)N (1,c)
对于结构的非线性以及非线性程度,通过比较结构可能的最大位移与屈服位移的大小得到;
所述第三步中采用实时混合模拟频域评价指标计算该窗函数上的时滞和幅值误差:
<mrow>
<mi>F</mi>
<mi>E</mi>
<mi>I</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</munderover>
<mo>{</mo>
<mfrac>
<mrow>
<msub>
<mi>y</mi>
<mi>o</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>y</mi>
<mi>I</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>y</mi>
<mi>I</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>p</mi>
</munderover>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>y</mi>
<mi>I</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>}</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>,</mo>
<mi>a</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mi>f</mi>
<mrow>
<mi>e</mi>
<mi>q</mi>
</mrow>
</msup>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>p</mi>
</munderover>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mo>|</mo>
<mo>(</mo>
<mrow>
<msub>
<mi>y</mi>
<mi>I</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>&CenterDot;</mo>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>p</mi>
</munderover>
<mo>|</mo>
<mo>|</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>y</mi>
<mi>I</mi>
</msub>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>|</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>,</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
</mrow>
A=||FEI|| (2,c)
φ=arc tan[Im(FEI)/Re(FEI)] (2,d)
d=φ/(2π·feq) (2,e)
其中yI(f)和yo(f)分别代表计算位移和测量位移加窗后的傅里叶变换,p代表窗函数长度的一半,Im(·)和Re(·)分别代表FEI的虚部和实部,FEI代表频域评价指标,feq表示等效频率,A和d分别代表实时混合模拟的的幅值和等效时滞,φ表示实时混合模拟的相位差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510690506.4A CN105335610B (zh) | 2015-10-22 | 2015-10-22 | 一种计算实时混合模拟试验瞬时时滞和幅值误差的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510690506.4A CN105335610B (zh) | 2015-10-22 | 2015-10-22 | 一种计算实时混合模拟试验瞬时时滞和幅值误差的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105335610A CN105335610A (zh) | 2016-02-17 |
CN105335610B true CN105335610B (zh) | 2017-12-19 |
Family
ID=55286130
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510690506.4A Active CN105335610B (zh) | 2015-10-22 | 2015-10-22 | 一种计算实时混合模拟试验瞬时时滞和幅值误差的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105335610B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107356523B (zh) * | 2017-06-20 | 2019-08-20 | 东南大学 | 基于离散切线刚度估计的实时混合模拟试验反馈力修正方法 |
CN109827763B (zh) * | 2019-01-21 | 2020-09-11 | 东南大学 | 粘滞阻尼器检测效果评价方法 |
CN109883666B (zh) * | 2019-01-21 | 2020-09-11 | 东南大学 | 一种粘滞阻尼器检测位移输入信号修正方法 |
CN110376894B (zh) * | 2019-07-26 | 2023-03-31 | 哈尔滨工业大学 | 一种适用于实时混合试验的两级时滞补偿方法 |
CN111323185B (zh) * | 2020-03-16 | 2021-11-12 | 西安建筑科技大学 | 一种空间框架模型自适应实时混合仿真试验加载方法 |
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 |
CN114459712B (zh) * | 2022-01-11 | 2023-12-29 | 东南大学 | 一种基于自回归模型的地震模拟振动台实验补偿方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103076194A (zh) * | 2012-12-31 | 2013-05-01 | 东南大学 | 实时混合模拟试验效果的频域评价方法 |
-
2015
- 2015-10-22 CN CN201510690506.4A patent/CN105335610B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103076194A (zh) * | 2012-12-31 | 2013-05-01 | 东南大学 | 实时混合模拟试验效果的频域评价方法 |
Non-Patent Citations (8)
Title |
---|
"A frequency response analysis approach for quantitative assessment of actuator tracking for real-time hybrid simulation";Tong Guo等;《Smart Materials and Structures》;20140430;第23卷(第4期);第1-13页 * |
"Analysis of decimation techniques to improve computational efficiency of a frequency-domain evaluation approach for real-time hybrid simulation";Tong Guo等;《Smart Materials and Structures》;20140630;第14卷(第6期);第1197-1220页 * |
"Application of Frequency-domain Evaluation Method in Real-time Hybrid Simulation";TONG Guo等;《5th International Conference on Advances in Experimental Structural Engineering》;20131109;第1-7页 * |
"FEI based Online Actuator Tracking Assessment for Real-Time Hybrid";Tong Guo等;《6th International Conference on Advances in Experimental Structural Engineering》;20150802;第1-7页 * |
"Hanning自卷积窗函数及其谐波分析应用";温和等;《中国科学(E 辑: 技术科学)》;20090630;第39卷(第6期);第1190-1198页 * |
"On the Use of Windows for Harmonic Analysis with the Discrete Fourier Transform";FREDRIC J. HARRIS等;《PROCEEDINGS OF THE IEEE》;19780131;第66卷(第1期);第51-83页 * |
"实时混合模拟试验作动器位移追踪的频域评价方法";郭彤等;《工程力学》;20140430;第31卷(第4期);第171-177页 * |
"实时混合模拟试验的频域评估方法研究";徐伟杰;《万方数据企业知识服务平台》;20140918;第2、4节 * |
Also Published As
Publication number | Publication date |
---|---|
CN105335610A (zh) | 2016-02-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105335610B (zh) | 一种计算实时混合模拟试验瞬时时滞和幅值误差的方法 | |
CN105588718B (zh) | 机床主轴综合性能检测/监测试验系统及方法 | |
CN103575523B (zh) | 基于FastICA-谱峭度-包络谱分析的旋转机械故障诊断方法 | |
CN103175679B (zh) | 四旋翼飞行器旋翼特性综合测试系统 | |
CN104034407B (zh) | 减小旋转机械扭振信号脉冲测量法中周期性误差的方法 | |
CN109323791B (zh) | 基于增量切割法的复合材料板材的残余应力分布测量方法 | |
CN102759448A (zh) | 基于柔性时域平均的齿轮箱故障检测方法 | |
CN104914846B (zh) | 基于自适应滑模观测器的电连接器间歇性失效检测方法 | |
CN105136423A (zh) | 考虑摩擦力的自由振动动导数试验的数据分析方法 | |
CN103345158B (zh) | 通风盘式制动器试验台及其电惯量模拟控制方法 | |
CN106768763A (zh) | 一种板损伤快速检测方法 | |
CN102253240A (zh) | 一种石英挠性加速度计安全使用边界确定方法 | |
CN103076194A (zh) | 实时混合模拟试验效果的频域评价方法 | |
CN103335831A (zh) | 一种制动器电惯量模拟试验台及其电惯量模拟控制方法 | |
CN105116003A (zh) | 一种利用热膨胀曲线计算两相转变比例的测量方法 | |
CN105571874A (zh) | 一种发动机振动不平衡相位的实时测量方法 | |
CN208606915U (zh) | 一种钢管构件双轴残余应力的测量装置 | |
CN102507055B (zh) | 一种小试样插销试验机裂纹试验断裂功的测量方法 | |
CN103344371B (zh) | 电惯量轮胎制动试验台及其电惯量控制方法 | |
CN103528775A (zh) | 一种基于响应灵敏度的结构健康检测方法 | |
CN101762422B (zh) | 拉力试验机拉伸测试方法 | |
CN107748817B (zh) | 一种考虑非比例附加强化的高温多轴本构关系确定方法 | |
CN106124196A (zh) | 一种风机齿轮箱传动链效率实验测试方法 | |
Wang et al. | Study on the propagation of coupling shock waves with phase transition under combined tension-torsion impact loading | |
CN104573197B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |