TWI220923B - A method for continuously measuring the damping ratio of structural systems - Google Patents

A method for continuously measuring the damping ratio of structural systems Download PDF

Info

Publication number
TWI220923B
TWI220923B TW92128786A TW92128786A TWI220923B TW I220923 B TWI220923 B TW I220923B TW 92128786 A TW92128786 A TW 92128786A TW 92128786 A TW92128786 A TW 92128786A TW I220923 B TWI220923 B TW I220923B
Authority
TW
Taiwan
Prior art keywords
signal
damping ratio
ratio
damping
phase
Prior art date
Application number
TW92128786A
Other languages
Chinese (zh)
Other versions
TW200514974A (en
Inventor
Wen-Lung Li
Original Assignee
Wen-Lung Li
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 Wen-Lung Li filed Critical Wen-Lung Li
Priority to TW92128786A priority Critical patent/TWI220923B/en
Application granted granted Critical
Publication of TWI220923B publication Critical patent/TWI220923B/en
Publication of TW200514974A publication Critical patent/TW200514974A/en

Links

Landscapes

  • Measuring Phase Differences (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

The present patent is a method for measuring the damping ratio of a structural system that is under base excitations. The method is accomplished by processing the response signal with respect to that of the input. In other words, one is to measure the responses of the system, which is excited by an input sinusoidal signal. Once the response signal is obtained, it is modulated by two different reference signals. The first one is the row input signal, while the second one is the phase-shifted input signal. The proposed phase given here is 90 degrees. The two modulated output signals are then put into the low-pass filter separately and their dc values are the ones the present method required. Dividing the dc values each other to get the tangent value of the phase lag, which is due the system damping ratio, one is able to use the equation provided to calculate the required damping ratio. Because the present method is based on the phase angle of the steady-state responses, it can be adopted and displayed in real time. In addition, the current method does not require the system damping to be small and valid for a wide range of damping ratios.

Description

1220923 五、發明說明(l) " ----- [技術領域] 本創作為一種適用於機械或土木結構系統之阻尼比量 檢測方法,特別適用於基座激振之系統,利用本創作可以 將實測之訊號,連續、即時轉換成系統之阻尼比,故 作可以增加結構動態量測儀器之相關功能。 ♦剧 [先前技術] 一般習知之系統阻尼量測方法,大都採取以暫離 方式量測,例如,對被量測系統施加一衝 二 其振幅縮小之程度,或以對數表現之對 等,即令時間t丨盥t夕屏;^ \ 以乐数方法 <ti,則對數遞減係數6為 1且t】 :Kl)51220923 V. Description of the Invention (l) " ----- [Technical Field] This creation is a damping ratio detection method suitable for mechanical or civil structure systems, especially for systems that are excited by a base. Use this creation The measured signal can be continuously and instantly converted into the damping ratio of the system, so it can increase the relevant functions of the structure dynamic measurement instrument. ♦ Drama [Prior art] Most commonly used system damping measurement methods are mostly measured in a temporary way, for example, applying a shock to the system being measured to reduce the amplitude of the amplitude, or a logarithmic equivalent, so that Time t 丨 wash t evening screen; ^ \ With the method of music number < ti, the logarithmic decreasing factor 6 is 1 and t]: Kl) 5

δ (1) 式中之Γ即為系統之阻尼比係數。即 度,然後假設系統阻尼甚小,忽略 2里振幅縮小之程 算求出。或者以勁度除以第一個 項=,再由U)式反 之比rA反推得系統阻尼比: 又田或最大振幅),所得 (2a) ·' 另外一類的方法,則將測 域,以下列方式表示: J系統之反應轉換到頻 (a)半功率或Q-值之頻寬八ω與系 比,或近似之系統之阪尼比係數表’-'、、、自…、'頻率ωη之 1220923 五'發明說明(2) (2b) (b)頻域曲線之適合(fitti 適合、(ϋ) Nyquist圖之圓適望(1)共振區段曲線 不管疋前述之哪一類的方法, 假設或適用於系統之阻尼比遠小於;t 缺點:(1) 多數之方法無法即時量測並顯示量測之糸統;或(2)絕大 時控制無法發择最大的作用。本創作=比’以致於即 方法,完全沒有習知各種方々 ^出之阻尼比量測 量測方法。 7缺點與限制,為一全新的 [發明内容及其理論推背景] 牲ιΛΛΓ夕審查6委員瞭解本創作之理論背景及其正確性, 、下列之内容說明本創作之主要推導過程: 參考圖一之系統,令一單自由度(SD〇F)之系統以基座 激振(54)u(t)作用,則其運動動方程式可以表示成: wz + c(i ~ //) + k(z ~ //) = 0 (3) 上式中之m代表系統質量(51 )、C代表系統阻尼(52)及k為 系統勁度(53)。故系統響應z(t) (55)可表示成 ζ(ί) l + j(2rc) a-r2) + /(2rC)δ (1) where Γ is the damping ratio coefficient of the system. That is, it is assumed that the system damping is very small, and it is calculated by ignoring the amplitude reduction process in 2 miles. Or divide the stiffness by the first term =, and then inversely derive the system damping ratio from the inverse ratio rA of U): (the field or the maximum amplitude), and obtain (2a) · 'Another type of method is to measure the domain, Expressed in the following way: The response of the J system is converted to the frequency (a) half-power or Q-value bandwidth eight ω and the ratio, or the approximate Hansby coefficient table of the system '-' ,,, ..., ' Frequency ωη of 1220923 Five 'invention description (2) (2b) (b) fit of frequency domain curve (fitti fit, (ϋ) circle of Nyquist diagram is desirable) (1) resonance section curve regardless of 疋It is assumed or applicable that the damping ratio of the system is much less than; t Disadvantages: (1) Most methods cannot measure and display the measurement system in real time; or (2) Control cannot select the maximum effect when it is extremely large. This creation = Ratio, so that it is the method, there is no familiarity with various methods of measuring the damping ratio measurement. 7 Disadvantages and limitations, a brand new [invention content and theoretical background] 6 members understand The theoretical background of this creation and its correctness, the following content explains the main derivation of this creation Process: Referring to the system of Figure 1, let a single-degree-of-freedom (SD0F) system act on the base (54) u (t), then its equation of motion can be expressed as: wz + c (i ~ / /) + k (z ~ //) = 0 (3) In the above formula, m represents the system mass (51), C represents the system damping (52), and k is the system stiffness (53). Therefore, the system response z (t ) (55) can be expressed as ζ (ί) l + j (2rc) a-r2) + / (2rC)

UeiQt = ZeK〇x^ (4)倘若基座激振u(t)(54)是頻率Ω諧波,或以複數式表示UeiQt = ZeK〇x ^ (4) If the base excitation u (t) (54) is a frequency Ω harmonic, or it is expressed in a complex formula

第9頁 1220923 五 發明說明(3) (5 ) u(t) = UeiQt ^式中…π 。另外,由於版尼存在,故造成輸入激振 (54)與系統響應間γ55)之相位差ρ正切,洽可表示成 (6) tM1^= (i-r2) + (2rc7 前述這些式子裡,已經使兩了下列之參數符號: (7 - 9) l4mk 、% =; 、 r"~ 其中(7 )式之ζ即為量測目標之阻尼比。 為更進一步說明本創作之背景理論,令輸入訊號(54)❿ (10) w(i) = i/sia(Qf) 則系統穩態響應(55)可表示成 (11) z(t)= Ζύη(Ωί-φ) 因此,若將之系統穩態響應(5 5 )與激振訊號(5 4 )做乘積調 變(60),或寫成 (12) g(t) = u(t) z(t) = Usin(pt) Ζήη(α-φ) 但, (13) 上式也可以利用三角等式寫成 厂.C/2 g(t) = —-[cos^ - cos(2Ωί - φ)[ 其中,rz由(4)式定義或Page 9 1220923 5 Description of the invention (3) (5) u (t) = UeiQt ^ In the formula ... π. In addition, due to the existence of the version, the phase difference ρ tangent between the input excitation (54) and the system response γ55) can be expressed as (6) tM1 ^ = (i-r2) + (2rc7) The following parameter symbols have been used: (7-9) l4mk,% = ;, r " ~ where ζ of the formula (7) is the damping ratio of the measurement target. To further explain the background theory of this creation, Let the input signal (54) ❿ (10) w (i) = i / sia (Qf) then the steady-state response (55) of the system can be expressed as (11) z (t) = ZOύ (Ωί-φ). Therefore, if The system steady-state response (5 5) and the excitation signal (5 4) are multiplied by (60), or written as (12) g (t) = u (t) z (t) = Usin (pt) ήη ( α-φ) However, (13) The above formula can also be written as a factory using a triangle equation. C / 2 g (t) = —- [cos ^-cos (2Ωί-φ) [where rz is defined by (4) or

r 〇 -λιιζ (14) Z |_(1 一 r2)2 +刚 \ 觀察(13)式,中括弧内之第一項為非時變項,故若以一低r 〇 -λιιζ (14) Z | _ (1 a r2) 2 + Gang \ Observe (13), the first term in the brackets is not a time-varying term, so if a low

第10頁 1220923 五、發明說明(4) 故 通滤波器(low〜Pass filter,LPF,61)將它過據出 (15)Page 10 1220923 V. Description of the invention (4) Therefore, the low-pass filter (LPF, 61) passes it out (15)

_IIlpF cosd 上式裡之n I I ·· I I "表示過濾後之非時變值。 另一方面,若將激振訊號(54)再乘積調變前,先做四 分之一週期或p/2相移(phase shift) (62),以作為系 統響應之第一個參考訊號,再依前述相同之積货變(^ 、 低通過濾(64)模式,可得 、 (16) 以及 g(〇 = u(i~r/4)«z(i)_IIlpF cosd n I I ·· I I in the above formula represents a time-invariant value after filtering. On the other hand, if the excitation signal (54) is subjected to product modulation, a quarter cycle or p / 2 phase shift (62) is performed as the first reference signal of the system response. According to the same product change (^, low-pass filtering (64) mode as above, we can obtain, (16), and g (〇 = u (i ~ r / 4) «z (i)

&_LPF 兩訊號相除 (17) ⑽ 2 因此,由(15)及(17)式,並比較(6)式後 (65)下式成立& _LPF Division of two signals (17) ⑽ 2 Therefore, from (15) and (17), and after comparing (6), (65) holds

G 2ζΓ (18) (i-n^(2rcy 同時,也可以由重新安排(18)式、並解出系統阻尼比( ,r2 ±^jr4 -4G2(l-r2)G 2ζΓ (18) (i-n ^ (2rcy At the same time, the equation (18) can also be rearranged and the system damping ratio (, r2 ± ^ jr4 -4G2 (l-r2)

(19) , 4rG 圖二說明上述之計算方法過程與創作構想。 此外,由於系統阻尼比為一正實數值,故依r之選 擇區域不同,下式必須成立: '(19), 4rG Figure 2 illustrates the above calculation process and creative concept. In addition, because the system's damping ratio is a positive real value, the following formula must be established according to the choice of r:

第11頁 1220923 五、發明說明(5) 2 if r <l.〇 (20) J^TG 〇>10 再者’觀察(18)式,當 r—1.0 (即 , #、 G〜G 从曰…咸a 、丨阻尼比Γ值很小時),則Page 11 1220923 V. Description of the invention (5) 2 if r < l.〇 (20) J ^ TG 〇 > 10 Furthermore, observe the formula (18), when r—1.0 (ie, #, G ~ G Since ... a, the damping ratio Γ is very small), then

來二,。,若當r 小值時,G是可以用Q ίϊ於:i計算式仍然適用。另外,依據傅利葉 :數’輸入之激振波可以視任意之週期波,並不一定是諧 [實施方式] 整理前述之理論背景,本創作之實施方式主要在將受 測糸統以一已知的週期諧波激振,同時以感測器(71)擷取 輸入訊號(54),參考圖一;另一方面,當受測系統受到激 振((54)後,其響應訊號(55),可以另一感測器(72)擷 取。接下來,可以利用軟體或電路硬體,依循圖二或本創 作之方法程序處理即可。圖三及圖四分別本創作發展時之 實驗訊號例放大示意圖,其中圖三為經由迴路54 6〇61 之訊號;圖四則為相同訊號但以迴路54_62_63處理後之 訊號。 另外,由於本創作之量測方法乃是在時域進行,故可 以即時量測、即時顯示,圖五說明本創作可即時量測之概 念。假定經時間L後、在某時間u進行一次量測,依據 實驗之取樣速度A t,只要在下一個取樣時間(12二七^ + ^ΊΙΗ 第12頁 1220923 五、發明說明(6) △ t ),即可依據本創作,進行下一循環之量測。即每隔 △ t即完成圖一之所有訊號迴路,並輸出顯示於螢幕,故 本創作可以連續、即時達成阻尼比之量測。 ,唯以上所述者,僅為本創作之較佳實施方式例而已 當不能以之限定本創作所實施之範圍,即大凡依本創作申 =專利範®所作之均#變化與修飾,皆應㈣於本創作專 :::之範圍内’謹請t審查委員明鑑,並祈惠准,是 1220923 圖式簡單說明 [圖示說明] 圖一本創作之阻尼比量測流程方塊圖。 圖二本創作之阻尼比量測方法之訊號處理方塊圖。 圖三本創作實施例之經由迴路5 4 - 6 0 - 61處理之訊號放大 示意圖。 圖四本創作實施例之經由迴路54- 62-63處理之訊號放大 示意圖。 圖五本創作實施例連續量測法之實施概念說明圖。 圖號說明: 50 系 統 之 參考 基 準 或 固定 面 51 系 統 質 量 52 系 統 阻 尼 53 系 統 勁 度 54 m 入 激振 訊 號 55 系 統 之 輸 出 響 應 訊 號 60 積 調 變 處 理 61 低 通 濾 波 處 理 62 訊 號 之 相 移 處理 63 積 調 變 處 理 64 低 通 濾 波 處 理 65 訊 號 乘 除 處 理 66 依 本創 作 提供 之 公 式計 算處理 71 感 測 器 72 感 測 器 Ω 輸 入 激 振 訊 號 之 頻 率 g(t) 經積調變處理後之訊號Come two. If G is small, G can be used for Q ίϊ: i is still applicable. In addition, according to the Fourier: number input shock wave can be viewed as any periodic wave, not necessarily harmonic [Embodiment] To organize the theoretical background, the implementation of this creation is mainly to test the system to a known The periodic harmonic is excited, and the sensor (71) is used to capture the input signal (54), refer to Figure 1. On the other hand, when the system under test is excited ((54), its response signal (55) , Can be captured by another sensor (72). Next, you can use software or circuit hardware to process according to Figure 2 or the method of this creation. Figure 3 and Figure 4 are the experimental signals during the development of this creation. The enlarged diagram of the example, in which Figure 3 is the signal via circuit 54 6061; Figure 4 is the same signal but processed by circuit 54_62_63. In addition, since the measurement method of this creation is performed in the time domain, it can be performed immediately Measurement, real-time display, Figure 5 illustrates the concept of real-time measurement in this creation. Assume that after time L, a measurement is performed at a certain time u, and according to the sampling speed A t of the experiment, as long as the next sampling time (12 27 ^ + ^ Ί ΙΗ Page 12 1220923 V. Description of the invention (6) △ t), according to this creation, the next cycle of measurement can be performed. That is, all the signal circuits of Figure 1 are completed every △ t, and the output is displayed on the screen. Therefore, this creation can achieve the measurement of the damping ratio continuously and in real time. However, the above mentioned are only examples of the best implementation of this creation and should not be used to limit the scope of this creation. = Patent 范 ® 作 之 均 # Variations and modifications should all be within the scope of this creative book ::: 'I would like to invite the reviewing committee to make a clear reference, and pray for your approval. It is a simple illustration of the 1220923 diagram. A block diagram of the damping ratio measurement process of this creation. Figure 2 A block diagram of signal processing of the damping ratio measurement method of this creation. Fig. 3 An enlarged schematic diagram of the signal processed by the circuit 5 4-6 0-61 in this embodiment of the creation. Figure 4 The enlarged schematic diagram of the signal processed by the circuit 54- 62-63 in this creative embodiment. Figure 5 The conceptual illustration of the implementation of the continuous measurement method in this creative embodiment. Figure No. Description: 50 Reference datum or fixed surface of the system 51 System mass 52 System damping 53 System stiffness 54 m Induced vibration signal 55 System output response signal 60 Product modulation processing 61 Low-pass filtering processing 62 Signal phase shift processing 63 Product modulation processing 64 Low-pass filtering 65 Signal multiplication and division Process 66 Calculate process 71 according to the formula provided by this creation 71 Sensor 72 Sensor Ω Input frequency g (t) of the excitation signal Signal after product modulation processing

第14頁Page 14

Claims (1)

1220923 六、申請專利範圍 [專利範圍] 1. 一種連續量測系統阻尼比之量測或檢測方法,係將系統 以弦波或週期波訊號激振,量測系統之響應訊號,再將響 應訊號分別與激振訊號、9 0度相移後之激振訊號做乘積調 變,分別做低通處理,擷取乘積調變後兩訊號之直流位 準,計算其比值,並代入下列公式計算所得者: r2±」r4 一 4G2a — r2) 公式中符號之代表意義分別為: · ζ* ··受測系統之阻尼比; Γ :激振訊號之頻率與系統自然頻率之比; G :乘積調變訊號之直流位準比值。 2.如申請專利範圍第1項所稱之連續量測系統阻尼比之量 測或檢測方法,得將結果立即顯示,即可達到具有即時量 測或監測之功能。1220923 VI. Scope of Patent Application [Patent Scope] 1. A method for measuring or detecting the damping ratio of a continuous measurement system, which excites the system with a sine wave or periodic wave signal, measures the response signal of the system, and then responds to the signal Perform product modulation with the excitation signal and the 90 ° phase-shifted excitation signal, respectively, and perform low-pass processing, respectively, extract the DC levels of the two signals after the product modulation, calculate the ratio, and substitute the following formula to calculate (R2 ± ”r4 4G2a — r2) The meanings of the symbols in the formula are: · ζ * · damping ratio of the system under test; Γ: ratio of the frequency of the excitation signal to the natural frequency of the system; G: product adjustment DC level ratio of variable signal. 2. The measurement or detection method of the damping ratio of the continuous measurement system, as described in item 1 of the scope of the patent application, can display the results immediately, which can achieve the function of instant measurement or monitoring. 第15頁Page 15
TW92128786A 2003-10-17 2003-10-17 A method for continuously measuring the damping ratio of structural systems TWI220923B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
TW92128786A TWI220923B (en) 2003-10-17 2003-10-17 A method for continuously measuring the damping ratio of structural systems

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
TW92128786A TWI220923B (en) 2003-10-17 2003-10-17 A method for continuously measuring the damping ratio of structural systems

Publications (2)

Publication Number Publication Date
TWI220923B true TWI220923B (en) 2004-09-11
TW200514974A TW200514974A (en) 2005-05-01

Family

ID=34132852

Family Applications (1)

Application Number Title Priority Date Filing Date
TW92128786A TWI220923B (en) 2003-10-17 2003-10-17 A method for continuously measuring the damping ratio of structural systems

Country Status (1)

Country Link
TW (1) TWI220923B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109632224A (en) * 2018-12-07 2019-04-16 武昌船舶重工集团有限公司 A kind of ship pedestal ratio of damping measuring device and measurement method

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109632224A (en) * 2018-12-07 2019-04-16 武昌船舶重工集团有限公司 A kind of ship pedestal ratio of damping measuring device and measurement method

Also Published As

Publication number Publication date
TW200514974A (en) 2005-05-01

Similar Documents

Publication Publication Date Title
JP3132628B2 (en) Coriolis mass flowmeter
JP2008046121A5 (en)
CN110907827B (en) Motor transient distortion measuring method and system
Li et al. Signal processing methods for Coriolis Mass Flow Metering in two-phase flow conditions
TWI220923B (en) A method for continuously measuring the damping ratio of structural systems
KR100839436B1 (en) The method of power frequency estimation using the difference between the gain and cosine and sine filter
Xu et al. Frequency-feature based antistrong-disturbance signal processing method and system for vortex flowmeter with single sensor
JP2013024577A (en) Angle detection method by resolver
CN104020350A (en) Voltage fundamental component detection method for overcoming frequency perturbation
Lu Lock-in frequency measurement with high precision and efficiency
CN110263482A (en) A kind of vortex impedance method for solving and device based on cross correlation algorithm
JP2787078B2 (en) measuring device
CN108957118A (en) A kind of reactive power calculating method
JP5394945B2 (en) Synchronous verification device
KR101538738B1 (en) Apparatus and method for frequency measurement using 3-level discreet fourier transform
JP7080757B2 (en) Impedance measuring device and impedance measuring method
JP2004294177A (en) Sampling type measuring device
JP2015122701A (en) Synchronization detection circuit and resistance measuring apparatus
CN110927504A (en) Analysis method, device and equipment of combined network based on LCR
WO2023139871A1 (en) Impedance measurement system and method
RU2577828C1 (en) Method of measuring phase shift of low frequency filter synchronous detector
TW201224469A (en) Method of measuring signal delay time
JP5022545B2 (en) Phase difference calculation method
WO2023139870A1 (en) Impedance measurement device and method
CN107421435A (en) Homodyne orthogonal laser interference signal generation method based on digit shifter

Legal Events

Date Code Title Description
MM4A Annulment or lapse of patent due to non-payment of fees