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 PDFInfo
- 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
Links
Landscapes
- Measuring Phase Differences (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
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)
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)
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 |
-
2003
- 2003-10-17 TW TW92128786A patent/TWI220923B/en not_active IP Right Cessation
Cited By (1)
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 |