JP2017078591A - Vibration control device and vibration control method, and vibration control program - Google Patents

Vibration control device and vibration control method, and vibration control program Download PDF

Info

Publication number
JP2017078591A
JP2017078591A JP2015205582A JP2015205582A JP2017078591A JP 2017078591 A JP2017078591 A JP 2017078591A JP 2015205582 A JP2015205582 A JP 2015205582A JP 2015205582 A JP2015205582 A JP 2015205582A JP 2017078591 A JP2017078591 A JP 2017078591A
Authority
JP
Japan
Prior art keywords
signal
series signal
frequency
drive
target
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
Application number
JP2015205582A
Other languages
Japanese (ja)
Other versions
JP6480848B2 (en
Inventor
文男 相原
Fumio Aihara
文男 相原
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Emic Corp
Original Assignee
Emic Corp
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 Emic Corp filed Critical Emic Corp
Priority to JP2015205582A priority Critical patent/JP6480848B2/en
Publication of JP2017078591A publication Critical patent/JP2017078591A/en
Application granted granted Critical
Publication of JP6480848B2 publication Critical patent/JP6480848B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

PROBLEM TO BE SOLVED: To provide a vibration control device, a vibration control method, and a vibration control program which are high in convergence performance of an excitation signal to be outputted, strictly following the specified energy of a measurement reference and also following the performance variation of a vibration exciter device during vibration testing and capable of adjusting the variation range within a set range.SOLUTION: A vibration control device is inputted with a PSD profile to be tested which is a spectrum distribution of the vibration generated in a vibration testing device, and generates a target time series signal from the PSD profile to be tested. The vibration control device includes: target signal generation means for generating a target frequency series signal Tf from the target time series signal; drive signal generation means for generating a drive time series signal from the target frequency series signal Tf; and drive signal correction means for correcting a drive frequency series signal Df which is a frequency spectrum of the drive time series signal so that the target frequency series signal Tf and a response frequency series signal Rf of the frequency spectrum of the response time series signal of the vibration testing device are matched.SELECTED DRAWING: Figure 8

Description

本発明は、工業製品に使用環境下と同様のランダム振動を実験室内で与えることにより、工業製品の耐振性及び耐久性を評価するランダム振動試験を行うための振動制御装置および振動制御方法ならびに振動制御プログラムに関する。   The present invention relates to a vibration control device, a vibration control method, and a vibration for performing a random vibration test for evaluating the vibration resistance and durability of an industrial product by giving the industrial product a random vibration similar to that in an operating environment. It relates to the control program.

従来のランダム振動制御装置では、制御目標として専ら、加速度平均パワースペクトル密度(PSD:Power Spectrum Density)、二乗平均平方根(RMS:Root Mean Square)が用いられている。
具体的には、試験に供する被試験PSDプロファイルの目標平均PSDと、振動試験装置の応答信号より得られる応答平均PSDとの比較により、目標平均PSDと応答平均PSDが一致するように制御を行う。
In the conventional random vibration control device, the acceleration average power spectrum density (PSD) and the root mean square (RMS) are exclusively used as control targets.
Specifically, control is performed so that the target average PSD and the response average PSD coincide with each other by comparing the target average PSD of the PSD profile to be tested to be tested and the response average PSD obtained from the response signal of the vibration test apparatus. .

図12に基づいて、瞬時PSDと平均PSDについて説明する。
まず、長時間にわたる時系列信号から、フーリエ変換の処理単位に相当する一定時間区切りのフレーム単位に切り出した時系列信号をフーリエ変換することにより、フレーム単位のフーリエスペクトルを計算する。時系列信号から周波数系列信号への変換にはフーリエ変換が用いられるのが一般的であるが、ここでは、フーリエ変換を高速で行う高速フーリエ変換(FFT変換)を用いることとして説明をする。
Based on FIG. 12, the instantaneous PSD and the average PSD will be described.
First, a Fourier spectrum in a frame unit is calculated by Fourier transforming a time series signal cut out from a time series signal over a long period of time into a frame unit of a fixed time interval corresponding to a Fourier transform processing unit. In general, Fourier transform is used for the conversion from the time series signal to the frequency series signal, but here, it will be described as using fast Fourier transform (FFT transform) that performs Fourier transform at high speed.

次いで、前記フーリエ変換で計算されたフーリエ係数から下記式(4)を用いて瞬時PSDを求め、下記式(5)から二乗平均平方根(RMS)を求める。   Next, an instantaneous PSD is obtained from the Fourier coefficient calculated by the Fourier transform using the following equation (4), and a root mean square (RMS) is obtained from the following equation (5).

Figure 2017078591
なお、上記式(4)において、Nは所要周波数範囲の分割個数である。
Figure 2017078591
In the above formula (4), N is the number of divisions in the required frequency range.

Figure 2017078591
このとき、各瞬時PSDを瞬時PSD_0(f),瞬時PSD_1(f),・・・,瞬時PSD_i(f)と表し、これら瞬時PSD_0(f)〜瞬時PSD_i(f)を時間平均したものが平均PSD(f)である。
Figure 2017078591
At this time, each instantaneous PSD is expressed as an instantaneous PSD_0 (f), an instantaneous PSD_1 (f),..., An instantaneous PSD_i (f). PSD (f).

図13に基づいて、従来の振動制御装置100の構成と動作について説明する。
まず、振動試験装置150からの応答時系列信号Rtが、AD変換部107によりデジタル信号に変換される。なお、振動試験装置150は、被試験体に振動を与え、この被試験体からの応答信号を計測する応答検出器(図示せず)を有している。
Based on FIG. 13, the structure and operation | movement of the conventional vibration control apparatus 100 are demonstrated.
First, the response time series signal Rt from the vibration test apparatus 150 is converted into a digital signal by the AD conversion unit 107. The vibration test apparatus 150 includes a response detector (not shown) that applies vibration to the device under test and measures a response signal from the device under test.

デジタル信号に変換された応答時系列信号RtをFFT変換部108においてフーリエ変換することにより瞬時PSDを計算し、瞬時PSDを所定時間における時間平均化処理することにより平均PSDを生成する。   The response time series signal Rt converted to the digital signal is Fourier-transformed by the FFT conversion unit 108 to calculate an instantaneous PSD, and the instantaneous PSD is time-averaged at a predetermined time to generate an average PSD.

そして、PSD比較部101において、被試験PSDプロファイルから得た目標PSDと、平均PSDとを比較することによって、PSD補正操作量を生成する。
次いで、目標PSDとPSD補正操作量とから修正PSDを生成し、振動試験装置の伝達特性の振幅特性に基づき逆振幅特性乗算処理部102において逆振幅特性乗算処理を施し、駆動PSDを生成する。
Then, the PSD comparison unit 101 generates a PSD correction manipulated variable by comparing the target PSD obtained from the PSD profile under test with the average PSD.
Next, a modified PSD is generated from the target PSD and the PSD correction manipulated variable, and a reverse amplitude characteristic multiplication process is performed in the reverse amplitude characteristic multiplication processing unit 102 based on the amplitude characteristic of the transfer characteristic of the vibration test apparatus to generate a drive PSD.

また、フーリエ逆変換部103において、乱数生成手段104で生成された白色ランダム信号を駆動PSDで正規化するとともに、フーリエ逆変換処理により時系列波形に変換、展開し、駆動波形合成手段105により時系列波形の重ね合わせ合成処理を行い、DA変換部106によって、振動試験装置150へ入力するアナログ信号である駆動時系列信号Dtが生成される。   Further, in the Fourier inverse transform unit 103, the white random signal generated by the random number generation unit 104 is normalized by the driving PSD, converted into a time series waveform by the Fourier inverse transformation process, and developed, and the drive waveform synthesizing unit 105 converts the time. A series waveform overlapping process is performed, and the DA converter 106 generates a drive time series signal Dt that is an analog signal input to the vibration test apparatus 150.

一般的に、振動制御装置では、振動試験装置の破損防止や安全性確保のため、駆動時系列信号により振動試験装置に与えられる変位が加振限界領域を超えると、自動的に停止し、試験が終了してしまうため、出来る限り安全許容内でのプロファイルで試験条件を設定することが必要とされる。   Generally, a vibration control device automatically stops when the displacement applied to the vibration test device by the drive time series signal exceeds the excitation limit region to prevent damage and ensure safety of the vibration test device. Therefore, it is necessary to set the test conditions with a profile within safety tolerance as much as possible.

通常の加速度振動試験においては、加振周波数範囲の被試験プロファイルとしてのPSDプロファイルと上記式(5)で示されるRMSとが与えられ、加速度時系列信号の確率密度分布は、RMSの2乗が分散σ2の正規分布となることが必要となる。 In a normal acceleration vibration test, a PSD profile as a profile to be tested in the excitation frequency range and an RMS represented by the above equation (5) are given, and the probability density distribution of the acceleration time series signal is the square of the RMS. A normal distribution with variance σ 2 is required.

ランダム信号は確率統計的に生成される信号である。そのため、このような確率密度分布を有する加速度時系列信号は、振動試験装置の許容範囲を超えた信号を含んでいることが予測される。   A random signal is a signal generated stochastically. Therefore, it is predicted that the acceleration time series signal having such a probability density distribution includes a signal that exceeds the allowable range of the vibration test apparatus.

例え、確率統計的には稀な現象であったとしても、振動試験装置の許容範囲を超えた信号が発生すると、振動試験装置に損傷を与える恐れがあるため、通常は変位センサーが振動試験装置の許容範囲の変位を超えた場合に異常を検出して、振動試験を自動停止させる機能が働くため、試験目的は完遂されない。このため、許容範囲外の加速度信号を除去するためにクリッピングと呼ばれる手法が用いられている。   For example, even if this phenomenon is rare in terms of probability statistics, if a signal exceeding the allowable range of the vibration test equipment is generated, the vibration test equipment may be damaged. Since the function of automatically detecting the abnormality when the displacement of the allowable range is exceeded and stopping the vibration test automatically works, the test purpose is not completed. For this reason, a technique called clipping is used to remove acceleration signals outside the allowable range.

クリッピングは、加速度ランダム時系列を分散の平方根であるσに適正倍率(クリッピング倍率)を乗じた値で規制して、過剰な加速度信号が発生しないように操作する処理である。クリッピング倍率が適正に設定されるならば、確率統計的には稀に発生する過剰加速度だけが抑制されるので、全体の加速度平均スペクトルには大きな影響を与えない。多くは3.0から5.0の倍率値が用いられている。   Clipping is a process of controlling an acceleration random time series with a value obtained by multiplying σ, which is the square root of variance, by an appropriate magnification (clipping magnification) so as not to generate an excessive acceleration signal. If the clipping magnification is set appropriately, only excessive acceleration that rarely occurs in terms of probability statistics is suppressed, so that the entire acceleration average spectrum is not greatly affected. In many cases, a magnification value of 3.0 to 5.0 is used.

クリッピング処理は、ランダム加速度振動試験において過剰な加速度信号の発生を直接抑止することが目的として使われるとともに、振動信号の振動物理範囲(変位)を抑止する効果を期待して用いられている。   The clipping process is used for the purpose of directly suppressing the generation of an excessive acceleration signal in a random acceleration vibration test, and is used in anticipation of the effect of suppressing the vibration physical range (displacement) of the vibration signal.

図14は、σクリッピングの有無による加速度振幅の確率分布の違いを表すグラフ、図15は、σクリッピングの有無による変位振幅の確率分布の違いを表すグラフである。図14に示すσクリッピングによって加速度を抑制したとしても、図15に示すように変位の抑制はほとんど出来ていない。即ち、クリッピングによる顕著な変位抑制効果は観測されない。   FIG. 14 is a graph showing the difference in probability distribution of acceleration amplitude depending on the presence or absence of σ clipping, and FIG. 15 is a graph showing the difference in probability distribution of displacement amplitude depending on the presence or absence of σ clipping. Even if the acceleration is suppressed by the σ clipping shown in FIG. 14, the displacement is hardly suppressed as shown in FIG. That is, no significant displacement suppression effect due to clipping is observed.

このように、加速度信号に対するσクリッピングは、加振力を抑制する手法としては有効であるが、変位(振動物理範囲)を抑制する手法としては、極めて限定的な効果しか発揮しない。   As described above, σ clipping with respect to the acceleration signal is effective as a technique for suppressing the excitation force, but has a very limited effect as a technique for suppressing the displacement (vibration physical range).

さらに、従来の振動制御装置においては、応答時系列信号はFFT処理によりフーリエスペクトルに変換されるが、その後すぐに瞬間PSDへと変換され、さらに瞬間PSDの時間平均としての応答平均PSDが計算されることになる。このとき、目標平均PSDとの比較において振動制御が行われるが、振動制御における重要な情報である位相成分は無視されていた。さらに、加速度信号からもたらされる振動信号の振動物理範囲(変位)は振動試験装置に係る重要な限界物理量であるにも関わらず制御の対象外となっていた。   Furthermore, in the conventional vibration control apparatus, the response time series signal is converted into a Fourier spectrum by FFT processing, but immediately after that, it is converted into an instantaneous PSD, and a response average PSD as a time average of the instantaneous PSD is calculated. Will be. At this time, vibration control is performed in comparison with the target average PSD, but phase components that are important information in vibration control are ignored. Furthermore, the vibration physical range (displacement) of the vibration signal resulting from the acceleration signal is outside the scope of control despite being an important limit physical quantity related to the vibration test apparatus.

また、ランダム振動試験において直接生成されるのは加速度ランダム信号であり、それに基づく変位信号は、ランダム加速度の無限二階積分となることから、誤差蓄積を免れず、理論的に特定することが困難であった。   In addition, the random signal generated directly in the random vibration test is an acceleration random signal, and the displacement signal based on it is an infinite second-order integral of random acceleration. there were.

また、振動試験装置に変位ドリフトを抑制・防止するための機構が設けられている場合、振動試験装置において発生する変位は、その機構の影響も受けてしまうことになる。   Further, when a mechanism for suppressing and preventing displacement drift is provided in the vibration test apparatus, the displacement generated in the vibration test apparatus is also affected by the mechanism.

もし、PSDプロファイルに基づいて生成されるランダム加速度信号が所定のPSDプロファイルとRMS値とを維持しつつ、変位を抑制できれば、振動試験装置の適用範囲を拡大することが可能となる。   If the random acceleration signal generated based on the PSD profile can suppress the displacement while maintaining the predetermined PSD profile and the RMS value, the application range of the vibration test apparatus can be expanded.

このように変位抑制を行う方法としては、従来、被試験PSDプロファイルから予め最大変位を推定し、許容性能範囲を超えている場合には、被試験PSDプロファイルの修正をするような予防的対策や、上述するようなクリッピングなどにより過大加速度信号のカットなどが用いられている。   As a method for suppressing displacement in this way, conventionally, a precautionary measure such as estimating the maximum displacement from the PSD profile under test in advance and correcting the PSD profile under test when the allowable performance range is exceeded, The over-acceleration signal is cut by clipping as described above.

しかしながら、加速度に対するクリッピング方式は加速度ランダムに対して確率的な規制を発揮するものであり、変位自体を規制するものではなく、図15に示すように、変位信号の抑制効果としては限定的なものであった。   However, the clipping method for acceleration exhibits probabilistic regulation for random acceleration, and does not regulate displacement itself. As shown in FIG. 15, the effect of suppressing the displacement signal is limited. Met.

本発明では、このような現状に鑑み、出力される加振信号の収束性が高く、測定基準の規定エネルギーに厳密に従い、かつ、振動試験中の加振装置や被試験体の特性変化にも追従し、変位範囲を設定された範囲内に収めることができる振動制御装置および振動制御方法ならびに振動制御プログラムを提供することを目的とする。   In the present invention, in view of such a current situation, the output excitation signal has high convergence, strictly follows the prescribed energy of the measurement standard, and changes in characteristics of the excitation apparatus and the DUT during the vibration test. An object of the present invention is to provide a vibration control device, a vibration control method, and a vibration control program that can follow and keep a displacement range within a set range.

本発明は、前述するような加速度ランダム振動試験に関する従来技術における課題を解決するために発明されたものであって、本発明の振動制御装置は、
被試験体に振動を与える振動試験装置の振動を制御するための振動制御装置であって、
前記振動試験装置において発生させる振動のスペクトル分布である被試験PSD(Power Spectrum Density)プロファイルが入力され、該被試験PSDプロファイルから目標周波数系列信号Tfを生成する目標信号生成手段と、
目標周波数系列信号Tfから駆動時系列信号Dtを生成する駆動信号生成手段と、
前記目標周波数系列信号Tfと、前記振動試験装置の応答時系列信号Rtの周波数スペクトルである応答周波数系列信号Rfとが一致するように、前記駆動時系列信号Dtの周波数スペクトルである駆動周波数系列信号Dfを補正する駆動信号補正手段と、
を備え、
前記目標信号生成手段が、
前記被試験PSDプロファイルを持つランダム時系列信号を生成するランダム信号生成部と、
前記ランダム時系列信号から、目標時系列信号Ttを生成する目標時系列信号生成部と、
前記目標時系列信号Ttから、目標周波数系列信号Tfを生成する目標周波数系列信号生成部と、
を含むことを特徴とする。
The present invention was invented in order to solve the problems in the prior art relating to the acceleration random vibration test as described above.
A vibration control device for controlling the vibration of a vibration test device that applies vibration to a device under test,
A target signal generating means for receiving a PSD (Power Spectrum Density) profile, which is a spectrum distribution of vibration generated in the vibration test apparatus, and generating a target frequency series signal Tf from the PSD profile under test;
Drive signal generating means for generating a drive time series signal Dt from the target frequency series signal Tf;
The drive frequency series signal that is the frequency spectrum of the drive time series signal Dt so that the target frequency series signal Tf matches the response frequency series signal Rf that is the frequency spectrum of the response time series signal Rt of the vibration test apparatus. Drive signal correction means for correcting Df;
With
The target signal generating means is
A random signal generator for generating a random time series signal having the PSD profile under test;
A target time series signal generating unit for generating a target time series signal Tt from the random time series signal;
A target frequency sequence signal generator for generating a target frequency sequence signal Tf from the target time series signal Tt;
It is characterized by including.

この場合、前記駆動信号生成手段が、
前記目標周波数系列信号Tfから、駆動周波数系列信号Dfを生成する周波数系列信号処理部と、
前記駆動周波数系列信号Dfから、前記駆動時系列信号Dtを生成する時系列信号への再変換処理部と、
を含むことが好ましい。
In this case, the drive signal generating means is
A frequency sequence signal processing unit that generates a drive frequency sequence signal Df from the target frequency sequence signal Tf;
A re-conversion processing unit from the driving frequency series signal Df to a time series signal for generating the driving time series signal Dt;
It is preferable to contain.

また、前記駆動信号補正手段が、前記周波数系列信号処理部を含み、
前記周波数系列信号処理部が、
前記目標周波数系列信号Tfと、前記振動試験装置の応答検出器により検出された応答時系列信号Rtから生成される応答スペクトルである応答周波数系列信号Rfとが一致するように、前記駆動時系列信号Dtから生成される駆動スペクトルである前記駆動周波数系列信号Dfを補正するように構成されていることが好ましい。
Further, the drive signal correction means includes the frequency series signal processing unit,
The frequency series signal processing unit,
The drive time series signal so that the target frequency series signal Tf and the response frequency series signal Rf, which is a response spectrum generated from the response time series signal Rt detected by the response detector of the vibration test apparatus, coincide with each other. It is preferable that the driving frequency series signal Df which is a driving spectrum generated from Dt is corrected.

また、本発明の振動制御方法は、
被試験体に振動を与える振動試験装置の振動を制御するための振動制御方法であって、
前記振動試験装置において発生させる振動のスペクトル分布である被試験PSD(Power Spectrum Density)プロファイルから、目標時系列信号Ttを生成するとともに、該目標時系列信号Ttから目標周波数系列信号Tfを生成し、
前記目標周波数系列信号Tfから駆動時系列信号Dtを生成し、
前記目標周波数系列信号Tfと、前記振動試験装置の応答時系列信号Rtの周波数スペクトルである応答周波数系列信号Rfとが一致するように、前記駆動時系列信号Dtの周波数スペクトルである駆動周波数系列信号Dfを補正するとともに、
前記目標時系列信号Ttに対する変位信号を常に生成し、変位領域の設定限界値を常に判定して、該変位限界値を超えずに、かつ、ランダム振動、目標PSDプロファイル、RMS値を一定に維持し続ける加速度ランダム信号を生成することを特徴とする。
Further, the vibration control method of the present invention includes:
A vibration control method for controlling vibration of a vibration test apparatus that applies vibration to a device under test,
A target time series signal Tt is generated from a PSD (Power Spectrum Density) profile that is a spectrum distribution of vibrations generated in the vibration test apparatus, and a target frequency series signal Tf is generated from the target time series signal Tt.
Generating a driving time series signal Dt from the target frequency series signal Tf;
The drive frequency series signal that is the frequency spectrum of the drive time series signal Dt so that the target frequency series signal Tf matches the response frequency series signal Rf that is the frequency spectrum of the response time series signal Rt of the vibration test apparatus. While correcting Df,
A displacement signal for the target time series signal Tt is always generated, the set limit value of the displacement region is always determined, and the random vibration, the target PSD profile, and the RMS value are maintained constant without exceeding the displacement limit value. A random acceleration signal that continues to be generated is generated.

この場合、前記目標周波数系列信号Tfから、駆動周波数系列信号Dfを生成するとともに、
前記駆動周波数系列信号Dfから、前記駆動時系列信号Dtを生成することが好ましい。
In this case, while generating the drive frequency sequence signal Df from the target frequency sequence signal Tf,
The drive time series signal Dt is preferably generated from the drive frequency series signal Df.

また、前記目標周波数系列信号Tfと、前記振動試験装置の応答検出器により検出された応答時系列信号Rtから生成される応答スペクトルである応答周波数系列信号Rfとが一致するように、前記駆動時系列信号Dtから生成される駆動スペクトルである前記駆動周波数系列信号Dfを補正することが好ましい。   The target frequency series signal Tf and the response frequency series signal Rf, which is a response spectrum generated from the response time series signal Rt detected by the response detector of the vibration test apparatus, coincide with each other during the driving. It is preferable to correct the drive frequency series signal Df which is a drive spectrum generated from the series signal Dt.

また、本発明の振動制御プログラムは、
被試験体に振動を与える振動試験装置の振動を制御するための振動制御プログラムであって、
前記振動試験装置において発生させる振動のスペクトル分布である被試験PSD(Power Spectrum Density)プロファイルから、目標時系列信号Ttを生成するとともに、該目標時系列信号Ttから目標周波数系列信号Tfを生成し、さらに、前記目標時系列信号Ttに対する変位信号を常に生成し、変位領域の設定限界値を常に判定して、該変位限界値を超えずに、かつ、ランダム振動、目標PSDプロファイル、RMS値を一定に維持し続ける加速度ランダム信号を生成する目標信号生成手順と、
目標周波数系列信号Tfから駆動時系列信号Dtを生成する駆動信号生成手順と、
前記目標周波数系列信号Tfと、前記振動試験装置の応答時系列信号Rtの周波数スペクトルである応答周波数系列信号Rfとが一致するように、前記駆動時系列信号Dtの周波数スペクトルである駆動周波数系列信号Dfを補正する駆動信号補正手順と、
をコンピュータに実行させることを特徴とする。
The vibration control program of the present invention is
A vibration control program for controlling vibration of a vibration test apparatus that applies vibration to a device under test,
A target time series signal Tt is generated from a PSD (Power Spectrum Density) profile that is a spectrum distribution of vibrations generated in the vibration test apparatus, and a target frequency series signal Tf is generated from the target time series signal Tt. Further, a displacement signal for the target time series signal Tt is always generated, the set limit value of the displacement region is always determined, and the random vibration, the target PSD profile, and the RMS value are kept constant without exceeding the displacement limit value. A target signal generation procedure for generating an acceleration random signal that continues to be maintained,
A drive signal generation procedure for generating a drive time series signal Dt from the target frequency series signal Tf;
The drive frequency series signal that is the frequency spectrum of the drive time series signal Dt so that the target frequency series signal Tf matches the response frequency series signal Rf that is the frequency spectrum of the response time series signal Rt of the vibration test apparatus. A drive signal correction procedure for correcting Df;
Is executed by a computer.

この場合、前記駆動信号生成手順が、
前記目標周波数系列信号Tfから、駆動周波数系列信号Dfを生成する周波数系列信号処理手順と、
前記駆動周波数系列信号Dfから、前記駆動時系列信号Dtを生成する時系列信号復号処理手順と、
を含むことが好ましい。
In this case, the drive signal generation procedure is:
A frequency sequence signal processing procedure for generating a drive frequency sequence signal Df from the target frequency sequence signal Tf;
A time series signal decoding processing procedure for generating the drive time series signal Dt from the drive frequency series signal Df;
It is preferable to contain.

また、前記駆動信号補正手順が、前記周波数系列信号処理手順を含み、
前記周波数系列信号処理手順が、
前記目標周波数系列信号Tfと、前記振動試験装置の応答検出器により検出された応答時系列信号Rtから生成される応答スペクトルである応答周波数系列信号Rfとが一致するように、前記駆動時系列信号Dtから生成される駆動スペクトルである前記駆動周波数系列信号Dfを補正することが好ましい。
The drive signal correction procedure includes the frequency sequence signal processing procedure,
The frequency sequence signal processing procedure includes:
The drive time series signal so that the target frequency series signal Tf and the response frequency series signal Rf, which is a response spectrum generated from the response time series signal Rt detected by the response detector of the vibration test apparatus, coincide with each other. It is preferable to correct the drive frequency series signal Df which is a drive spectrum generated from Dt.

本発明によれば、目標周波数系列信号と応答周波数系列信号とが一致するように、駆動周波数系列信号を補正することで、振動試験装置の加速度頻度分布を抑制することなく、変位頻度分布のみを抑制でき、振動試験装置の振動物理範囲(変位)を超えずに安定かつ連続的な振動試験を行うことができる。   According to the present invention, by correcting the drive frequency sequence signal so that the target frequency sequence signal matches the response frequency sequence signal, only the displacement frequency distribution can be obtained without suppressing the acceleration frequency distribution of the vibration test apparatus. Stable and continuous vibration test can be performed without exceeding the vibration physical range (displacement) of the vibration test apparatus.

図1は、本発明の一実施例における振動制御装置の制御系を説明するための制御プロセス図である。FIG. 1 is a control process diagram for explaining a control system of a vibration control apparatus according to an embodiment of the present invention. 図2は、被試験PSDプロファイルの一例を示すグラフである。FIG. 2 is a graph showing an example of the PSD profile under test. 図3は、ランダム信号生成部及び目標時系列信号生成部において、被試験PSDプロファイルからピュアランダム信号を生成する流れを説明するためのフロー図である。FIG. 3 is a flowchart for explaining a flow of generating a pure random signal from the PSD profile under test in the random signal generation unit and the target time series signal generation unit. 図4は、フーリエスペクトル成分と複素平面の偏角θの関係を示すグラフである。FIG. 4 is a graph showing the relationship between the Fourier spectrum component and the deviation angle θ of the complex plane. 図5は、ピュアランダム時系列信号の一例を示すグラフである。FIG. 5 is a graph showing an example of a pure random time series signal. 図6は、変位抑止された加速度ランダム信号を生成する流れを説明するためのフロー図である。FIG. 6 is a flowchart for explaining a flow of generating an acceleration random signal whose displacement is suppressed. 図7は、RMS=12.5m/sec2で加速度ランダム試験を実行した際のPSDプロファイルと、応答信号の平均PSDとを示すグラフである。FIG. 7 is a graph showing the PSD profile when the acceleration random test is executed at RMS = 12.5 m / sec 2 and the average PSD of the response signal. 図8は、周波数系列信号処理部の構成を説明するための処理ブロック図である。FIG. 8 is a processing block diagram for explaining the configuration of the frequency sequence signal processing unit. 図9は、時系列信号への再変換処理部の構成を説明するための処理ブロック図である。FIG. 9 is a processing block diagram for explaining the configuration of the time-series signal reconversion processing unit. 図10は、本実施例の振動制御装置を用いて振動試験を行った際の加速度振幅の確率分布を示すグラフである。FIG. 10 is a graph showing a probability distribution of acceleration amplitude when a vibration test is performed using the vibration control apparatus of the present embodiment. 図11は、本実施例の振動制御装置を用いて振動試験を行った際の変位振幅の確率分布を示すグラフであり、本実施例により35%もの変位抑止効果が確認された。FIG. 11 is a graph showing a probability distribution of displacement amplitude when a vibration test is performed using the vibration control apparatus of the present embodiment. As a result, a displacement suppression effect of 35% was confirmed by the present embodiment. 図12は、瞬時PSDと平均PSDについて説明するための概念図である。FIG. 12 is a conceptual diagram for explaining the instantaneous PSD and the average PSD. 図13は、従来の振動制御装置の構成と動作を説明するための制御プロセス図である。FIG. 13 is a control process diagram for explaining the configuration and operation of a conventional vibration control apparatus. 図14は、クリッピングの有無による加速度振幅の確率分布の違いを表すグラフ。FIG. 14 is a graph showing a difference in probability distribution of acceleration amplitude depending on the presence or absence of clipping. 図15は、クリッピングの有無による変位振幅の確率分布の違いを表すグラフである。FIG. 15 is a graph showing the difference in probability distribution of displacement amplitude depending on the presence or absence of clipping.

以下、本発明の実施の形態(実施例)を図面に基づいてより詳細に説明する。
図1は、本実施例における振動制御装置10の制御系を説明するための制御プロセス図である。
Hereinafter, embodiments (examples) of the present invention will be described in more detail with reference to the drawings.
FIG. 1 is a control process diagram for explaining a control system of the vibration control device 10 in the present embodiment.

図1に示すように、振動制御装置10は、ランダム信号生成部12と、目標時系列信号生成部14と、目標スペクトル生成部16と、応答スペクトル生成部18と、駆動スペクトル生成部20と、周波数系列信号処理部22と、時系列信号への再変換処理部24と、DA変換部26と、AD変換部28と、を備えている。   As shown in FIG. 1, the vibration control device 10 includes a random signal generation unit 12, a target time series signal generation unit 14, a target spectrum generation unit 16, a response spectrum generation unit 18, a drive spectrum generation unit 20, A frequency series signal processing unit 22, a time series signal re-conversion processing unit 24, a DA conversion unit 26, and an AD conversion unit 28 are provided.

また、振動制御装置10は、DA変換部26及びAD変換部28を介して、被試験体に振動を与える振動試験装置50と接続されている。なお、振動試験装置50は、パワー増幅器及び応答検出器(図示せず)を含んでおり、DA変換部26から出力される加振信号がパワー増幅器を介して振動試験装置に入力されて、振動試験装置は加振信号に基づいて振動を発生させる。そして、その振動に応じて得られる被試験体からの応答振動を、応答検出器によって検出し、応答アナログ信号としてAD変換部28へ入力する。   Further, the vibration control device 10 is connected to a vibration test device 50 that applies vibration to the device under test via the DA conversion unit 26 and the AD conversion unit 28. The vibration test apparatus 50 includes a power amplifier and a response detector (not shown), and the vibration signal output from the DA converter 26 is input to the vibration test apparatus via the power amplifier to generate vibration. The test apparatus generates vibration based on the excitation signal. Then, the response vibration from the DUT obtained according to the vibration is detected by the response detector and input to the AD converter 28 as a response analog signal.

以下、各プロセスについて各々説明する。
(1)ランダム信号生成部12〜目標時系列信号生成部14
ランダム信号生成部12には、例えば、図2に示すような、振動試験装置50において発生させる振動のスペクトル分布である被試験PSDプロファイルが入力される。本実施例の振動制御装置10は、この被試験入力プロファイルに基づいて、加振信号を生成する。
Hereinafter, each process will be described.
(1) Random signal generator 12 to target time series signal generator 14
For example, a random signal generator 12 receives a PSD profile to be tested, which is a spectral distribution of vibrations generated in the vibration test apparatus 50 as shown in FIG. The vibration control apparatus 10 according to the present embodiment generates an excitation signal based on the input profile under test.

ランダム信号生成部12では、フレーム毎に、被試験PSDプロファイルを持つランダム信号を生成し、これを順次繰り返すことにより、被試験PSDプロファイルが持続するランダム時系列信号が生成される。   The random signal generator 12 generates a random signal having a PSD profile under test for each frame, and sequentially repeats this to generate a random time-series signal that maintains the PSD profile under test.

なお、ランダム時系列信号が、完全な白色雑音であるための必要条件は、時間領域での出現頻度が一様ランダムで、かつ、周波数領域での振幅分布がランダムで、かつ、周波数領域での位相分布が一様ランダムであることである。   The necessary condition for the random time-series signal to be completely white noise is that the frequency of appearance in the time domain is uniformly random, the amplitude distribution in the frequency domain is random, and the frequency domain The phase distribution is uniformly random.

本実施例のように、図2で示す被試験PSDプロファイルで定義されたランダム振動試験では、周波数領域での振幅分布が被試験PSDプロファイルに置き換えられ、時間領域での出現頻度がランダムで尚且つ周波数領域での位相分布がランダムであるランダム信号が生成されることになる。   As in this example, in the random vibration test defined by the PSD profile under test shown in FIG. 2, the amplitude distribution in the frequency domain is replaced with the PSD profile under test, the frequency of appearance in the time domain is random, and A random signal having a random phase distribution in the frequency domain is generated.

図3は、ランダム信号生成部12及び目標時系列信号生成部14において、被試験PSDプロファイルからピュアランダム信号を生成する流れを説明するためのフロー図である。
まず、被試験PSDプロファイルから位相がランダムなフーリエスペクトルを生成する。
FIG. 3 is a flowchart for explaining a flow of generating a pure random signal from the PSD profile under test in the random signal generation unit 12 and the target time series signal generation unit 14.
First, a Fourier spectrum having a random phase is generated from the PSD profile under test.

時系列信号f(t)と、そのフーリエ変換関数F(ω)は、下記式(1)で示す関係を有する。   The time series signal f (t) and its Fourier transform function F (ω) have a relationship represented by the following formula (1).

Figure 2017078591
上記式(1)において、n=・・・,−3,−2,−1,0,1,2,3,・・・である。
Figure 2017078591
In the above formula (1), n =..., −3, −2, −1, 0, 1, 2, 3,.

また、離散的フーリエ変換であるFFT処理後の関数における実数成分aと虚数成分bは、下記式(2)、(3)で示す関係を有する。   Moreover, the real component a and the imaginary component b in the function after the FFT process, which is a discrete Fourier transform, have a relationship represented by the following equations (2) and (3).

Figure 2017078591
Figure 2017078591

Figure 2017078591
とすれば、
Figure 2017078591
given that,

Figure 2017078591
Figure 2017078591

また、フーリエスペクトルの絶対値は式(3)、(4)を用いて、下記式(6)のように表すことができる。   Further, the absolute value of the Fourier spectrum can be expressed by the following equation (6) using equations (3) and (4).

Figure 2017078591
なお、上記式(4)において、Nは所要周波数範囲の分割個数である。
Figure 2017078591
In the above formula (4), N is the number of divisions in the required frequency range.

Figure 2017078591
Figure 2017078591

この絶対値を基準として、フーリエスペクトルの偏角をθとした時、PSD[n]を示すフーリエスペクトル成分は、下記式(7)と表すことができる。   With this absolute value as a reference, when the declination of the Fourier spectrum is θ, the Fourier spectrum component indicating PSD [n] can be expressed by the following formula (7).

Figure 2017078591
Figure 2017078591

すなわち、フーリエスペクトル成分と複素平面の偏角θの関係は、図4に示すようになる。
PSDは、図4に示す絶対値(PSD値)により決定されるが、位相値であるθには自由度が存在するため、位相値θをランダムに変動させることで、PSD分布が閉架しないランダム位相フーリエスペクトルを生成することができる。
That is, the relationship between the Fourier spectrum component and the deviation angle θ of the complex plane is as shown in FIG.
The PSD is determined by the absolute value (PSD value) shown in FIG. 4, but since there is a degree of freedom in the phase value θ, the PSD distribution is not randomized by randomly varying the phase value θ. A phase Fourier spectrum can be generated.

ランダム位相フーリエスペクトルを逆FFT変換することで、位相がバラバラな正弦波が重なり合った、被試験PSDプロファイルのスペクトル分布を有する正規分布ランダム時系列信号が生成される。   By subjecting the random phase Fourier spectrum to inverse FFT transformation, a normally distributed random time series signal having a spectrum distribution of the PSD profile under test in which sine waves with different phases are overlapped is generated.

しかしながら、デジタル信号処理で使われる逆FFT変換は離散的な処理であるため、全ての周波数成分を漏れなく含むわけではなく、フレーム位相が不規則に変化することによって欠落した周波数成分が埋め合わされる効果がある。   However, since the inverse FFT transform used in the digital signal processing is a discrete process, it does not include all frequency components without omission, and the missing frequency components are compensated for by irregularly changing the frame phase. effective.

それでも、ある程度の周波数成分は欠落してしまうため、正規分布ランダム時系列信号は擬似ランダム信号と呼ばれる。そこで、擬似ランダム信号からピュアランダム信号を生成するために、次の処理を施す。   Nevertheless, since a certain amount of frequency components are lost, the normally distributed random time series signal is called a pseudo-random signal. Therefore, the following processing is performed in order to generate a pure random signal from the pseudo random signal.

まず、擬似ランダム信号のフレーム毎の時系列信号に、ハン窓(Hann WindowまたはHanning Window)を適用する。
そして、ハン窓が適用されたランダム時系列信号を1/2フレーム分ずつオーバーラップ合成処理する。
First, a Hann window (Hanning Window) is applied to the time-series signal of each pseudo-random signal.
The random time series signal to which the Hann window is applied is overlap-combined by 1/2 frame.

ハン窓を適用する処理と、1/2フレーム分ずつオーバーラップ合成処理とにより、時系列が混ぜ合わさり、図5に示すような、被試験PSDプロファイルを維持したピュアランダム時系列信号が生成される。   The time series is mixed by the process of applying the Hann window and the overlap synthesis process by 1/2 frame, and a pure random time series signal maintaining the PSD profile under test as shown in FIG. 5 is generated. .

目標時系列信号生成部14は、ピュアランダム信号を目標時系列信号として、目標スペクトル生成部16に送り出す。   The target time series signal generation unit 14 sends the pure random signal as a target time series signal to the target spectrum generation unit 16.

(2)目標スペクトル生成部16
目標スペクトル生成部では、目標時系列信号をFFT処理することにより、目標瞬時フーリエスペクトル、目標瞬時PSD、目標平均PSDが求められ、これら全ての情報を有する目標周波数系列信号を生成する。
(2) Target spectrum generator 16
In the target spectrum generation unit, a target instantaneous Fourier spectrum, a target instantaneous PSD, and a target average PSD are obtained by performing FFT processing on the target time series signal, and a target frequency series signal having all these information is generated.

n+ibnは位相成分を含み、また、フレーム毎に求められるものであるため、これを瞬時フーリエスペクトルと呼ぶ。また、瞬時フーリエスペクトルから上記式(4)に基づいて瞬時PSDが算出され、瞬時PSDに対する時間平均を取ることで平均PSDが算出される。 Since a n + ib n includes a phase component and is obtained for each frame, this is called an instantaneous Fourier spectrum. Also, an instantaneous PSD is calculated from the instantaneous Fourier spectrum based on the above equation (4), and an average PSD is calculated by taking a time average with respect to the instantaneous PSD.

すなわち、目標スペクトル生成部16では、目標時系列信号Ttが入力されることにより、式(1)〜(7)に基づき、目標瞬時フーリエスペクトル、目標瞬時PSD、目標平均PSDが算出され、これら全ての情報を有する目標周波数系列信号Tfが生成される。   That is, the target spectrum generation unit 16 calculates the target instantaneous Fourier spectrum, the target instantaneous PSD, and the target average PSD based on the equations (1) to (7) by inputting the target time series signal Tt. A target frequency sequence signal Tf having the following information is generated.

以上の方法で生成された目標加速度ランダムの変位越えは、確率統計的にも稀にしか発生しない。しかし、実際の装置を用いた加振試験においては、一回でも変位越えが発生すると機械に損傷を与える恐れが生じるため、その時点で試験を停止せざるを得なくなる。このため、変位越えが発生しない加速度ランダムな目標時系列信号Ttを生成し続けるための目標信号生成手段は、下記の処理を行う加速度ランダム信号生成部をさらに設けることもできる。   The displacement of the target acceleration random generated by the above method rarely occurs in terms of probability statistics. However, in the vibration test using an actual apparatus, if the displacement is exceeded even once, the machine may be damaged, and the test must be stopped at that time. For this reason, the target signal generation means for continuously generating the acceleration random target time-series signal Tt that does not cause displacement excess can further include an acceleration random signal generation unit that performs the following processing.

図6は、変位抑止された加速度ランダム目標信号を生成する流れを説明するためのフロー図である。
加速度ランダム信号は、まず、フレーム毎に生成され、順次接続され周期的繰り返しの無い連続信号となる。
FIG. 6 is a flowchart for explaining a flow of generating a random acceleration target signal in which displacement is suppressed.
The acceleration random signal is first generated for each frame, and is sequentially connected to be a continuous signal without periodic repetition.

ここでは、フレームiが生成された次の段階として、変位越えが発生しないフレームi+1の加速度ランダム信号の生成について説明する。
フレームi+1の加速度ランダム信号は、上述するランダム信号生成部12〜目標時系列信号生成部14において生成される。
Here, generation of an acceleration random signal of frame i + 1 in which displacement over does not occur will be described as the next stage after frame i is generated.
The acceleration random signal of the frame i + 1 is generated in the random signal generation unit 12 to the target time series signal generation unit 14 described above.

同時に、フレームiの最終変位データを初期値として数値2階積分を行い正確な変位予測を行う。もし、この仮定での変位予測が予め設定された変位許容範囲を超えたならば、フレームi+1の加速度ランダム信号を破棄して、再度同じ手順で加速度ランダム信号を生成し、同様の変位予測を行う。   At the same time, accurate displacement prediction is performed by performing second-order numerical integration using the final displacement data of frame i as an initial value. If the displacement prediction under this assumption exceeds the preset allowable displacement range, the acceleration random signal of frame i + 1 is discarded, and the acceleration random signal is generated again by the same procedure, and the same displacement prediction is performed. .

変位予測値が予め設定された許容範囲に収まったならば、ここで生成された加速度ランダムのフレームデータをフレームi+1として採用する。この処理により所望の変位抑制がなされた加速度ランダム信号の生成がおこなわれる。   If the predicted displacement value falls within the preset allowable range, the acceleration random frame data generated here is adopted as the frame i + 1. This process generates an acceleration random signal in which a desired displacement is suppressed.

もし、求められる変位抑止量が40%を超えると、上述する方法を数回連続で試みても変位を抑止できない場合が発生する。この場合には、5Hz以下の低周波数領域のPSD値を適宜下げることにより変位の抑止を行うことができる。現実的にこのようなケースは大きな抑止量を設定した場合で稀にしか発生せず、平均PSDに大きな影響を与えるものではない。   If the required displacement suppression amount exceeds 40%, the displacement may not be suppressed even if the above-described method is tried several times in succession. In this case, displacement can be suppressed by appropriately reducing the PSD value in the low frequency region of 5 Hz or less. In reality, such a case rarely occurs when a large amount of deterrence is set, and does not greatly affect the average PSD.

このようにして得られた加速度ランダム信号を、次の元データとして順次処理を進める事により、変位抑止された加速度ランダム信号を連続的に生成することができる。   By sequentially processing the acceleration random signal thus obtained as the next original data, it is possible to continuously generate the acceleration random signal in which the displacement is suppressed.

以上の変位抑止アルゴリズムにより、PSDプロファイルで規定された加速度ランダム加振における、以下の成立条件を満たしつつ、実際の変位範囲を抑止することが可能となる。   With the above displacement suppression algorithm, it is possible to suppress the actual displacement range while satisfying the following conditions for the random acceleration of acceleration defined by the PSD profile.

このようにして生成された変位抑止の加速度ランダム目標信号は、図7に示すように、
(A)計測された応答信号の平均PSDが指定されたPSDプロファイルと一致する。
(B)RMS値が指定された値と一致する。
(C)加速度時系列の確率密度分布が正規分布となる。
(D)同じ信号の繰り返しが存在しない。
という加速度ランダム加振信号に求められる基本要件を完全に満たしている。
図11に示される通り、本実施例の振動制御装置を用いて加速度振動試験を行った結果では35%もの変位抑止効果が確認された。
The acceleration random target signal for displacement suppression generated in this way is as shown in FIG.
(A) The average PSD of the measured response signal matches the designated PSD profile.
(B) The RMS value matches the specified value.
(C) The probability density distribution of the acceleration time series is a normal distribution.
(D) There is no repetition of the same signal.
The basic requirements for an acceleration random vibration signal are completely satisfied.
As shown in FIG. 11, as a result of an acceleration vibration test using the vibration control apparatus of this example, a displacement suppression effect of 35% was confirmed.

(3)周波数系列信号処理部22
図8は、周波数系列信号処理部22の構成を説明するための処理ブロック図である。
図8に示すように、周波数系列信号処理部22は、目標平均PSDと応答平均PSDとを比較する平均PSD比較部221と、目標瞬時PSDを補正して補正目標周波数系列信号Tcを生成するスペクトル補正制御部222と、補正目標周波数系列信号Tcに逆伝達関数IGfを乗算することで駆動信号Drを生成する逆伝達関数乗算部223と、応答周波数系列信号Rfと駆動周波数系列信号Dfの比である伝達関数Gfを算出する伝達関数計算部224と、伝達関数Gfの逆特性である逆伝達関数IGfを算出する逆伝達関数計算部225と、を備えている。
(3) Frequency sequence signal processing unit 22
FIG. 8 is a processing block diagram for explaining the configuration of the frequency sequence signal processing unit 22.
As shown in FIG. 8, the frequency sequence signal processing unit 22 includes an average PSD comparison unit 221 that compares the target average PSD and the response average PSD, and a spectrum that corrects the target instantaneous PSD to generate a corrected target frequency sequence signal Tc. The correction control unit 222, the inverse transfer function multiplier 223 that generates the drive signal Dr by multiplying the correction target frequency sequence signal Tc by the inverse transfer function IGf, and the ratio of the response frequency sequence signal Rf and the drive frequency sequence signal Df A transfer function calculation unit 224 that calculates a certain transfer function Gf and an inverse transfer function calculation unit 225 that calculates an inverse transfer function IGf that is an inverse characteristic of the transfer function Gf are provided.

(4)駆動周波数信号生成部
平均PSD比較部221には、目標周波数系列信号Tfと応答周波数系列信号Rfが入力されることで、目標平均PSDと応答平均PSDとが比較され、それらの誤差である平均PSD補正ベクターTCvがフレーム毎に計算される。そして、目標周波数系列信号Tfの目標瞬時PSDに対する補正量として、スペクトル補正制御部222へと送られる。
(4) Drive frequency signal generation unit The average PSD comparison unit 221 receives the target frequency sequence signal Tf and the response frequency sequence signal Rf, and compares the target average PSD with the response average PSD. An average PSD correction vector TCv is calculated for each frame. Then, it is sent to the spectrum correction control unit 222 as a correction amount for the target instantaneous PSD of the target frequency series signal Tf.

スペクトル補正制御部222では、目標周波数系列信号Tfの目標瞬時PSDに平均PSD補正ベクターTCvを加算することで、フレーム毎に目標瞬時PSDが補正された補正目標周波数系列信号Tcが生成される。   The spectrum correction control unit 222 adds the average PSD correction vector TCv to the target instantaneous PSD of the target frequency series signal Tf, thereby generating a corrected target frequency series signal Tc in which the target instantaneous PSD is corrected for each frame.

このように構成することにより、スペクトル補正制御部222は、位相成分の補正を行ってはいないが、目標瞬時PSDの振幅補正として振動試験装置の特性変動に対して機敏な制御動作を行うことになる。   By configuring in this way, the spectrum correction control unit 222 does not correct the phase component, but performs an agile control operation with respect to the characteristic variation of the vibration test apparatus as the amplitude correction of the target instantaneous PSD. Become.

また、伝達関数計算部224では、応答周波数系列信号Rfと駆動周波数系列信号Dfの比である伝達関数Gfとその平均値を常時算出している。そして、伝達関数Gfの変化が所定量を超えた場合、最新の伝達関数Gfが逆伝達関数計算部225へと送られる。   In addition, the transfer function calculation unit 224 constantly calculates a transfer function Gf that is a ratio of the response frequency sequence signal Rf and the drive frequency sequence signal Df and an average value thereof. When the change in the transfer function Gf exceeds a predetermined amount, the latest transfer function Gf is sent to the inverse transfer function calculation unit 225.

逆伝達関数計算部225では、入力された伝達関数Gfに基づき、逆伝達関数IGfが計算され、逆伝達関数乗算部223へと送られ、上述するように、補正目標周波数系列信号Tcに対して逆伝達関数IGfが乗算され、駆動のためのフーリエスペクトル信号である駆動信号Drが生成される。   In the inverse transfer function calculation unit 225, the inverse transfer function IGf is calculated based on the input transfer function Gf, and is sent to the inverse transfer function multiplier 223. As described above, the inverse transfer function IGf is applied to the corrected target frequency sequence signal Tc. The inverse transfer function IGf is multiplied to generate a drive signal Dr that is a Fourier spectrum signal for driving.

(5)駆動時系列信号への再変換処理部24
図9は、時系列信号への再変換処理部24の構成を説明するための処理ブロック図である。
時系列信号への再変換処理部24では、IFFT変換部241において、駆動信号DrをIFFT変換処理により周波数系列信号から時系列信号に復号化する。そして、復号化された時系列信号は、フレーム間が不連続であるため、オーバーラップ処理部242において、オーバーラップ合成処理を施すことにより、連続的な一次駆動時系列信号Dtcを生成する。
(5) Re-conversion processing unit 24 for driving time series signals
FIG. 9 is a processing block diagram for explaining the configuration of the time-series signal reconversion processing unit 24.
In the time series signal reconversion processing unit 24, the IFFT conversion unit 241 decodes the drive signal Dr from the frequency series signal to the time series signal by IFFT conversion processing. Since the decoded time series signal is discontinuous between frames, the overlap processing unit 242 performs overlap synthesis processing to generate a continuous primary drive time series signal Dtc.

時系列信号への再変換処理部24の制御サイクルは、IFFT変換処理とオーバーラップ合成処理を除いて、制御サンプリングクロック(例えば、制御帯域2kHzの場合5.12×103sample/s)に同期して動作する。 The control cycle of the re-conversion processing unit 24 to the time series signal is synchronized with a control sampling clock (for example, 5.12 × 10 3 sample / s when the control band is 2 kHz) except for the IFFT conversion process and the overlap synthesis process. Works.

IFFT変換処理とオーバーラップ合成処理の演算は、時系列信号への再変換処理部24の動作のバックグラウンドにおいて処理され、この演算は1フレームの処理時間内に完了する。
そして、生成された駆動時系列信号Dtは、バッファーに蓄積され、次の制御サイクルにおいて使用される。
The operations of the IFFT conversion processing and the overlap synthesis processing are processed in the background of the operation of the time-series signal reconversion processing unit 24, and this calculation is completed within the processing time of one frame.
The generated drive time series signal Dt is accumulated in a buffer and used in the next control cycle.

周波数系列での信号演算では常に、論理的に求められる適正長(1フレーム)の時系列信号に対して一括してFFT演算処理が行われる。このため、周波数系列で演算された結果をIFFT変換処理した時系列信号では、フレーム間が不連続となってしまう。   In the signal calculation in the frequency series, the FFT calculation process is always performed on the time series signal of the appropriate length (one frame) that is logically obtained. For this reason, in a time-series signal obtained by performing IFFT conversion processing on the result calculated in the frequency series, the frames are discontinuous.

時系列信号は、通常、連続的であることを前提として制御演算が行われるため、フレーム間の不連続を解消し、周波数系列での信号演算で得られた時系列信号が連続性を保つようにオーバーラップ合成処理を行う。なお、本実施例では、1/2フレーム毎のオーバーラップ合成処理を行っている。   Since time-series signals are usually controlled on the assumption that they are continuous, discontinuity between frames is eliminated, and the time-series signals obtained by signal operations in the frequency sequence are kept continuous. The overlap synthesis process is performed. In this embodiment, the overlap synthesis process is performed every 1/2 frame.

また、目標・応答時系列信号比較部244では、目標時系列信号Ttと応答時系列信号Rtとを比較し、例えば、PID制御などに基づき、時系列信号補正量Dcを生成する。   The target / response time series signal comparison unit 244 compares the target time series signal Tt and the response time series signal Rt, and generates a time series signal correction amount Dc based on, for example, PID control.

時系列信号補正部243では、一次駆動時系列信号Dtcと時系列信号補正量Dcとに基づき補正処理が行われ、補正された一次駆動時系列信号Dtcは、駆動時系列信号Dtとして、図1に示すように、DA変換部26を介して振動試験装置50へ出力される。   The time series signal correction unit 243 performs correction processing based on the primary drive time series signal Dtc and the time series signal correction amount Dc, and the corrected primary drive time series signal Dtc is shown as a drive time series signal Dt in FIG. As shown in FIG. 6, the vibration is output to the vibration test apparatus 50 via the DA converter 26.

以上、このように構成された本実施例の振動制御装置10は、被試験PSDプロファイルが入力されることにより、ランダム信号生成部12、目標時系列信号生成部14、目標スペクトル生成部16を含む目標信号生成手段により、変位抑止の目的にかなった目標周波数系列信号Tfが生成される。   As described above, the vibration control apparatus 10 of the present embodiment configured as described above includes the random signal generation unit 12, the target time series signal generation unit 14, and the target spectrum generation unit 16 when the PSD profile to be tested is input. A target frequency series signal Tf suitable for the purpose of displacement suppression is generated by the target signal generation means.

次いで、周波数系列信号処理部22、時系列信号への再変換処理部24を含む駆動信号生成手段によりランダム時系列信号である駆動時系列信号Dtが生成される。
また、周波数系列信号処理部22を含む駆動信号補正手段により、目標周波数系列信号Tfと、振動試験装置50の応答検出器により検出された応答時系列信号Rtから生成される応答スペクトルである応答周波数系列信号Rfとが一致するように、駆動時系列信号Dtから生成される駆動スペクトルである駆動周波数系列信号Dfを補正することによって、新たな駆動信号Drが生成され、この駆動信号Drから駆動時系列信号生成手段において、新たな駆動時系列信号Dtが生成される。
Next, a drive time series signal Dt, which is a random time series signal, is generated by a drive signal generation unit including a frequency series signal processing unit 22 and a time series signal reconversion processing unit 24.
In addition, a response frequency that is a response spectrum generated by the drive signal correction means including the frequency sequence signal processing unit 22 from the target frequency sequence signal Tf and the response time series signal Rt detected by the response detector of the vibration test apparatus 50. A new drive signal Dr is generated by correcting the drive frequency series signal Df, which is a drive spectrum generated from the drive time series signal Dt, so that the series signal Rf matches, and the drive signal Dr is driven from this drive signal Dr. In the series signal generation means, a new drive time series signal Dt is generated.

このようなフィードバック制御を行うことによって、振動試験装置50における変位を抑制することができる。   By performing such feedback control, the displacement in the vibration test apparatus 50 can be suppressed.

(6)応答スペクトル生成部18
振動試験装置から入力される応答アナログ信号は、AD変換部により応答時系列信号Rtに変換され、応答スペクトル生成部へ送られる。
応答スペクトル生成部では、応答時系列信号Rtが入力されることにより、上記式(1)〜(6)に基づき、応答瞬時フーリエスペクトル、応答瞬時PSD、応答平均PSDが算出され、これら全ての情報を有する応答周波数系列信号Rfが生成される。
(6) Response spectrum generator 18
The response analog signal input from the vibration test apparatus is converted into a response time series signal Rt by the AD conversion unit and sent to the response spectrum generation unit.
The response spectrum generation unit calculates the response instantaneous Fourier spectrum, response instantaneous PSD, and response average PSD based on the above formulas (1) to (6) by inputting the response time series signal Rt, and all the information Is generated as a response frequency sequence signal Rf.

(7)駆動スペクトル生成部20
振動試験装置50の駆動信号として生成された駆動時系列信号Dtを、DA変換部により駆動アナログ信号に変換する前段で分岐し、駆動スペクトル生成部20に入力する。
すなわち、駆動スペクトル生成部20では、駆動時系列信号Dtが入力されることにより、上記式(1)〜(7)に基づき、駆動瞬時フーリエスペクトル、駆動瞬時PSD、駆動平均PSDが算出され、これら全ての情報を有する駆動周波数系列信号Dfが生成される。
(7) Drive spectrum generation unit 20
The drive time series signal Dt generated as the drive signal of the vibration test apparatus 50 is branched at the previous stage where it is converted into a drive analog signal by the DA converter and input to the drive spectrum generator 20.
That is, the drive spectrum generation unit 20 receives the drive time series signal Dt, and calculates the drive instantaneous Fourier spectrum, the drive instantaneous PSD, and the drive average PSD based on the above formulas (1) to (7). A drive frequency series signal Df having all information is generated.

以上、本発明の好ましい実施例を説明したが、本発明はこれに限定されることはなく、例えば、上述した各処理部については、装置構成に応じて、ハードウェアでも、ソフトウェアでも構わず、コンピュータにおいて上述する処理内容を実行するプログラムとしてもよいなど、本発明の目的を逸脱しない範囲で種々の変更が可能である。   The preferred embodiments of the present invention have been described above, but the present invention is not limited to this. For example, each processing unit described above may be hardware or software depending on the device configuration. Various modifications can be made without departing from the object of the present invention, such as a program for executing the above-described processing contents in a computer.

本実施例の振動制御装置10及び従来の振動制御装置を用いた振動試験を行った際の加速度振幅の確率分布と、変位振幅の確率分布を図10、11に示す。   FIGS. 10 and 11 show the probability distribution of the acceleration amplitude and the probability distribution of the displacement amplitude when the vibration test is performed using the vibration control device 10 of the present embodiment and the conventional vibration control device.

従来の振動制御装置を用いた場合、振動試験装置50に設けた変位計による計測では、振動試験装置50の最大変位は±115mmであった。   When a conventional vibration control apparatus is used, the maximum displacement of the vibration test apparatus 50 is ± 115 mm in measurement using a displacement meter provided in the vibration test apparatus 50.

一方、本実施例の振動制御装置10を用いた場合、振動試験装置50に設けた変位計による計測では、振動試験装置50の最大変位は±75mmに抑制することができ、100時間以上の連続動作でも35%の安定した変位抑止が確認された。   On the other hand, when the vibration control apparatus 10 of the present embodiment is used, the maximum displacement of the vibration test apparatus 50 can be suppressed to ± 75 mm in the measurement by the displacement meter provided in the vibration test apparatus 50, and it is continuously performed for 100 hours or more. In operation, stable displacement suppression of 35% was confirmed.

10 振動制御装置
12 ランダム信号生成部
14 目標時系列信号生成部
16 目標スペクトル生成部
18 応答スペクトル生成部
20 駆動スペクトル生成部
22 周波数系列信号処理部
221 平均PSD比較部
222 スペクトル補正制御部
223 逆伝達関数乗算部
224 伝達関数計算部
225 逆伝達関数計算部
24 時系列信号への再変換処理部
241 IFFT変換部
242 オーバーラップ処理部
243 時系列信号補正部
244 目標・応答時系列信号比較部
26 DA変換部
28 AD変換部
50 振動試験装置
100 振動制御装置
101 PSD比較部
102 逆振幅特性乗算処理部
103 フーリエ逆変換部
104 乱数生成手段
105 駆動波形合成手段
106 DA変換部
107 AD変換部
108 FFT変換部
150 振動試験装置
DESCRIPTION OF SYMBOLS 10 Vibration control apparatus 12 Random signal generation part 14 Target time series signal generation part 16 Target spectrum generation part 18 Response spectrum generation part 20 Drive spectrum generation part 22 Frequency series signal processing part 221 Average PSD comparison part 222 Spectral correction control part 223 Reverse transmission Function multiplication unit 224 Transfer function calculation unit 225 Inverse transfer function calculation unit 24 Re-conversion processing unit 241 IFFT conversion unit 242 Overlap processing unit 243 Time-series signal correction unit 244 Target / response time-series signal comparison unit 26 DA Conversion unit 28 AD conversion unit 50 Vibration test device 100 Vibration control device 101 PSD comparison unit 102 Inverse amplitude characteristic multiplication processing unit 103 Inverse Fourier transform unit 104 Random number generation unit 105 Drive waveform synthesis unit 106 DA conversion unit 107 AD conversion unit 108 FFT conversion 150 vibration test equipment

Claims (9)

被試験体に振動を与える振動試験装置の振動を制御するための振動制御装置であって、
前記振動試験装置において発生させる振動のスペクトル分布である被試験PSD(Power Spectrum Density)プロファイルが入力され、該被試験PSDプロファイルから目標時系列信号Ttを生成するとともに、該目標時系列信号Ttから目標周波数系列信号Tfを生成する目標信号生成手段と、
目標周波数系列信号Tfから駆動時系列信号Dtを生成する駆動信号生成手段と、
前記目標周波数系列信号Tfと、前記振動試験装置の応答時系列信号Rtの周波数スペクトルである応答周波数系列信号Rfとが一致するように、前記駆動時系列信号Dtの周波数スペクトルである駆動周波数系列信号Dfを補正する駆動信号補正手段と、
を備え、
前記目標信号生成手段が、
前記目標時系列信号Ttに対する変位信号を常に生成し、変位領域の設定限界値を常に判定して、該変位限界値を超えずに、かつ、時系列の確率密度の正規分布性、目標PSDプロファイル、RMS値を一定に維持し続ける加速度ランダム信号を生成する加速度ランダム信号生成部をさらに含むことを特徴とする振動制御装置。
A vibration control device for controlling the vibration of a vibration test device that applies vibration to a device under test,
A PSD (Power Spectrum Density) profile, which is a spectrum distribution of vibrations generated in the vibration test apparatus, is input, a target time series signal Tt is generated from the PSD profile under test, and a target is obtained from the target time series signal Tt. Target signal generating means for generating a frequency series signal Tf;
Drive signal generating means for generating a drive time series signal Dt from the target frequency series signal Tf;
The drive frequency series signal that is the frequency spectrum of the drive time series signal Dt so that the target frequency series signal Tf matches the response frequency series signal Rf that is the frequency spectrum of the response time series signal Rt of the vibration test apparatus. Drive signal correction means for correcting Df;
With
The target signal generating means is
A displacement signal for the target time series signal Tt is always generated, a set limit value of the displacement region is always determined, and the normal distribution of the time series probability density without exceeding the displacement limit value, the target PSD profile The vibration control device further includes an acceleration random signal generation unit that generates an acceleration random signal that keeps the RMS value constant.
前記駆動信号生成手段が、
前記目標周波数系列信号Tfから、駆動周波数系列信号Dfを生成する周波数系列信号処理部と、
前記駆動周波数系列信号Dfから、前記駆動時系列信号Dtを生成する時系列信号への再変換処理部と、
を含むことを特徴とする請求項1に記載の振動制御装置。
The drive signal generating means is
A frequency sequence signal processing unit that generates a drive frequency sequence signal Df from the target frequency sequence signal Tf;
A re-conversion processing unit from the driving frequency series signal Df to a time series signal for generating the driving time series signal Dt;
The vibration control device according to claim 1, comprising:
前記駆動信号補正手段が、前記周波数系列信号処理部を含み、
前記周波数系列信号処理部が、
前記目標周波数系列信号Tfと、前記振動試験装置の応答検出器により検出された応答時系列信号Rtから生成される応答スペクトルである応答周波数系列信号Rfとが一致するように、前記駆動時系列信号Dtから生成される駆動スペクトルである前記駆動周波数系列信号Dfを補正するように構成されていることを特徴とする請求項2に記載の振動制御装置。
The drive signal correction means includes the frequency series signal processing unit,
The frequency series signal processing unit,
The drive time series signal so that the target frequency series signal Tf and the response frequency series signal Rf, which is a response spectrum generated from the response time series signal Rt detected by the response detector of the vibration test apparatus, coincide with each other. The vibration control apparatus according to claim 2, wherein the vibration control apparatus is configured to correct the drive frequency series signal Df that is a drive spectrum generated from Dt.
被試験体に振動を与える振動試験装置の振動を制御するための振動制御方法であって、
前記振動試験装置において発生させる振動のスペクトル分布である被試験PSD(Power Spectrum Density)プロファイルから、目標時系列信号Ttを生成するとともに、該目標時系列信号Ttから目標周波数系列信号Tfを生成し、
前記目標周波数系列信号Tfから駆動時系列信号Dtを生成し、
前記目標周波数系列信号Tfと、前記振動試験装置の応答時系列信号Rtの周波数スペクトルである応答周波数系列信号Rfとが一致するように、前記駆動時系列信号Dtの周波数スペクトルである駆動周波数系列信号Dfを補正するとともに、
前記目標時系列信号Ttに対する変位信号を常に生成し、変位領域の設定限界値を常に判定して、該変位限界値を超えずに、かつ、ランダム振動、目標PSDプロファイル、RMS値を一定に維持し続ける加速度ランダム信号を生成することを特徴とする振動制御方法。
A vibration control method for controlling vibration of a vibration test apparatus that applies vibration to a device under test,
A target time series signal Tt is generated from a PSD (Power Spectrum Density) profile that is a spectrum distribution of vibrations generated in the vibration test apparatus, and a target frequency series signal Tf is generated from the target time series signal Tt.
Generating a driving time series signal Dt from the target frequency series signal Tf;
The drive frequency series signal that is the frequency spectrum of the drive time series signal Dt so that the target frequency series signal Tf matches the response frequency series signal Rf that is the frequency spectrum of the response time series signal Rt of the vibration test apparatus. While correcting Df,
A displacement signal for the target time series signal Tt is always generated, the set limit value of the displacement region is always determined, and the random vibration, the target PSD profile, and the RMS value are maintained constant without exceeding the displacement limit value. A vibration control method characterized by generating an acceleration random signal that continues.
前記目標周波数系列信号Tfから、駆動周波数系列信号Dfを生成するとともに、
前記駆動周波数系列信号Dfから、前記駆動時系列信号Dtを生成することを特徴とする請求項4に記載の振動制御方法。
A drive frequency series signal Df is generated from the target frequency series signal Tf,
The vibration control method according to claim 4, wherein the drive time series signal Dt is generated from the drive frequency series signal Df.
前記目標周波数系列信号Tfと、前記振動試験装置の応答検出器により検出された応答時系列信号Rtから生成される応答スペクトルである応答周波数系列信号Rfとが一致するように、前記駆動時系列信号Dtから生成される駆動スペクトルである前記駆動周波数系列信号Dfを補正することを特徴とする請求項5に記載の振動制御方法。   The drive time series signal so that the target frequency series signal Tf and the response frequency series signal Rf, which is a response spectrum generated from the response time series signal Rt detected by the response detector of the vibration test apparatus, coincide with each other. The vibration control method according to claim 5, wherein the drive frequency series signal Df, which is a drive spectrum generated from Dt, is corrected. 被試験体に振動を与える振動試験装置の振動を制御するための振動制御プログラムであって、
前記振動試験装置において発生させる振動のスペクトル分布である被試験PSD(Power Spectrum Density)プロファイルから、目標時系列信号Ttを生成するとともに、該目標時系列信号Ttから目標周波数系列信号Tfを生成し、さらに、前記目標時系列信号Ttに対する変位信号を常に生成し、変位領域の設定限界値を常に判定して、該変位限界値を超えずに、かつ、ランダム振動、目標PSDプロファイル、RMS値を一定に維持し続ける加速度ランダム信号を生成する目標信号生成手順と、
目標周波数系列信号Tfから駆動時系列信号Dtを生成する駆動信号生成手順と、
前記目標周波数系列信号Tfと、前記振動試験装置の応答時系列信号Rtの周波数スペクトルである応答周波数系列信号Rfとが一致するように、前記駆動時系列信号Dtの周波数スペクトルである駆動周波数系列信号Dfを補正する駆動信号補正手順と、
をコンピュータに実行させるための振動制御プログラム。
A vibration control program for controlling vibration of a vibration test apparatus that applies vibration to a device under test,
A target time series signal Tt is generated from a PSD (Power Spectrum Density) profile that is a spectrum distribution of vibrations generated in the vibration test apparatus, and a target frequency series signal Tf is generated from the target time series signal Tt. Further, a displacement signal for the target time series signal Tt is always generated, the set limit value of the displacement region is always determined, and the random vibration, the target PSD profile, and the RMS value are kept constant without exceeding the displacement limit value. A target signal generation procedure for generating an acceleration random signal that continues to be maintained,
A drive signal generation procedure for generating a drive time series signal Dt from the target frequency series signal Tf;
The drive frequency series signal that is the frequency spectrum of the drive time series signal Dt so that the target frequency series signal Tf matches the response frequency series signal Rf that is the frequency spectrum of the response time series signal Rt of the vibration test apparatus. A drive signal correction procedure for correcting Df;
Vibration control program for causing a computer to execute.
前記駆動信号生成手順が、
前記目標周波数系列信号Tfから、駆動周波数系列信号Dfを生成する周波数系列信号処理手順と、
前記駆動周波数系列信号Dfから、前記駆動時系列信号Dtを生成する時系列信号復号処理手順と、
を含むことを特徴とする請求項7に記載の振動制御プログラム。
The drive signal generation procedure includes:
A frequency sequence signal processing procedure for generating a drive frequency sequence signal Df from the target frequency sequence signal Tf;
A time series signal decoding processing procedure for generating the drive time series signal Dt from the drive frequency series signal Df;
The vibration control program according to claim 7, comprising:
前記駆動信号補正手順が、前記周波数系列信号処理手順を含み、
前記周波数系列信号処理手順が、
前記目標周波数系列信号Tfと、前記振動試験装置の応答検出器により検出された応答時系列信号Rtから生成される応答スペクトルである応答周波数系列信号Rfとが一致するように、前記駆動時系列信号Dtから生成される駆動スペクトルである前記駆動周波数系列信号Dfを補正することを特徴とする請求項8に記載の振動制御プログラム。
The drive signal correction procedure includes the frequency sequence signal processing procedure,
The frequency sequence signal processing procedure includes:
The drive time series signal so that the target frequency series signal Tf and the response frequency series signal Rf, which is a response spectrum generated from the response time series signal Rt detected by the response detector of the vibration test apparatus, coincide with each other. 9. The vibration control program according to claim 8, wherein the drive frequency series signal Df, which is a drive spectrum generated from Dt, is corrected.
JP2015205582A 2015-10-19 2015-10-19 Vibration control device, vibration control method, and vibration control program Active JP6480848B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2015205582A JP6480848B2 (en) 2015-10-19 2015-10-19 Vibration control device, vibration control method, and vibration control program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2015205582A JP6480848B2 (en) 2015-10-19 2015-10-19 Vibration control device, vibration control method, and vibration control program

Publications (2)

Publication Number Publication Date
JP2017078591A true JP2017078591A (en) 2017-04-27
JP6480848B2 JP6480848B2 (en) 2019-03-13

Family

ID=58666049

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2015205582A Active JP6480848B2 (en) 2015-10-19 2015-10-19 Vibration control device, vibration control method, and vibration control program

Country Status (1)

Country Link
JP (1) JP6480848B2 (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2018021781A (en) * 2016-08-02 2018-02-08 地方独立行政法人大阪産業技術研究所 Non-gaussian vibration control device
CN108195537A (en) * 2018-02-06 2018-06-22 东南大学 A kind of blade of aviation engine vibration fatigue test method based on amplitude controlling
WO2020158068A1 (en) * 2019-01-31 2020-08-06 Imv株式会社 Vibration control device
JP2021043033A (en) * 2019-09-10 2021-03-18 エミック株式会社 Vibration test device and failure prediction method thereof
JP2021096238A (en) * 2019-12-17 2021-06-24 エミック株式会社 Excitation capability prediction/evaluation device, excitation capability prediction/evaluation method, and excitation capability prediction/evaluation program for vibration test

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011027669A (en) * 2009-07-29 2011-02-10 Hitachi Plant Technologies Ltd Device and method for testing vibration
JP2012103240A (en) * 2010-10-14 2012-05-31 Microsignal Kk Random vibration test control apparatus
JP2012233814A (en) * 2011-05-06 2012-11-29 Microsignal Kk Vibration test control device
US20130231878A1 (en) * 2012-03-05 2013-09-05 Crystal Instruments Corporation System and method for generation and control of mehanical vibration

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011027669A (en) * 2009-07-29 2011-02-10 Hitachi Plant Technologies Ltd Device and method for testing vibration
JP2012103240A (en) * 2010-10-14 2012-05-31 Microsignal Kk Random vibration test control apparatus
JP2012233814A (en) * 2011-05-06 2012-11-29 Microsignal Kk Vibration test control device
US20130231878A1 (en) * 2012-03-05 2013-09-05 Crystal Instruments Corporation System and method for generation and control of mehanical vibration

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2018021781A (en) * 2016-08-02 2018-02-08 地方独立行政法人大阪産業技術研究所 Non-gaussian vibration control device
CN108195537A (en) * 2018-02-06 2018-06-22 东南大学 A kind of blade of aviation engine vibration fatigue test method based on amplitude controlling
WO2020158068A1 (en) * 2019-01-31 2020-08-06 Imv株式会社 Vibration control device
JP2020122734A (en) * 2019-01-31 2020-08-13 Imv株式会社 Vibration controller
US11879816B2 (en) 2019-01-31 2024-01-23 Imv Corporation Vibration control system
JP2021043033A (en) * 2019-09-10 2021-03-18 エミック株式会社 Vibration test device and failure prediction method thereof
JP7405393B2 (en) 2019-09-10 2023-12-26 エミック株式会社 Vibration test equipment and its failure prediction method
JP2021096238A (en) * 2019-12-17 2021-06-24 エミック株式会社 Excitation capability prediction/evaluation device, excitation capability prediction/evaluation method, and excitation capability prediction/evaluation program for vibration test
JP7086411B2 (en) 2019-12-17 2022-06-20 エミック株式会社 Vibration capacity prediction evaluation device, vibration capacity prediction evaluation method and vibration capacity prediction evaluation program in vibration test
JP7430433B2 (en) 2019-12-17 2024-02-13 エミック株式会社 Excitation ability prediction evaluation device, excitation ability prediction evaluation method, and excitation ability prediction evaluation program in vibration testing
JP7475076B2 (en) 2019-12-17 2024-04-26 エミック株式会社 Vibration capability prediction and evaluation device, vibration capability prediction and evaluation method, and vibration capability prediction and evaluation program for vibration testing
JP7490288B2 (en) 2019-12-17 2024-05-27 エミック株式会社 System for predicting and evaluating vibration capability in vibration testing, method for predicting and evaluating vibration capability, and program for predicting and evaluating vibration capability

Also Published As

Publication number Publication date
JP6480848B2 (en) 2019-03-13

Similar Documents

Publication Publication Date Title
JP6480848B2 (en) Vibration control device, vibration control method, and vibration control program
JP4429358B2 (en) System and method for simultaneously controlling the spectrum and kurtosis of random vibrations
Orallo et al. Harmonics measurement with a modulated sliding discrete Fourier transform algorithm
JP5490335B1 (en) Frequency response measuring device
US8942930B2 (en) System and method for generation and control of mechanical vibration
JP2012233814A (en) Vibration test control device
JPH0552219A (en) Magnetic bearing control method and its device
JP6518631B2 (en) Non-Gaussian Vibration Controller
EP1772739B1 (en) Method and apparatus for monitoring the condition of a machine
JP2012103240A (en) Random vibration test control apparatus
EP3571761B1 (en) Controlling a wind turbine generator
JP2008515059A (en) Frequency control in closed loop control systems.
JP4788656B2 (en) Power test system
CN110535091B (en) Microcomputer protection method and device adapting to frequency conversion characteristic of electrical equipment
JP4955196B2 (en) AC signal measuring device
Jiao et al. An approach for electrical harmonic analysis based on interpolation DFT
EP3199956A1 (en) Apparatus for determination of the frequency of an electrical signal and associated method
US20090261768A1 (en) System and method for adjusting pulse width modulation signals according to component temperatures in a computer
JP5363900B2 (en) Excitation control method, excitation control device, and vibration test device
JP4591333B2 (en) Measuring device for engine generated torque
KR101714480B1 (en) Phase estimating apparatus based on discrete fourier transform using variable sampling point
JP6214480B2 (en) Frequency response measuring device
Dunn et al. Sinewave parameter estimation using the fast fan-chirp transform
JP2011117794A (en) Device for measuring effective value
Kunac et al. A Laboratory Setup for Testing Frequency Estimation Algorithms Implemented on Real-Time Target

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20170919

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20170926

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20171122

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20180125

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20180626

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20180827

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20190129

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20190208

R150 Certificate of patent or registration of utility model

Ref document number: 6480848

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350