JP4747336B2 - Variable structure semi-active seismic isolation system - Google Patents

Variable structure semi-active seismic isolation system Download PDF

Info

Publication number
JP4747336B2
JP4747336B2 JP2005102494A JP2005102494A JP4747336B2 JP 4747336 B2 JP4747336 B2 JP 4747336B2 JP 2005102494 A JP2005102494 A JP 2005102494A JP 2005102494 A JP2005102494 A JP 2005102494A JP 4747336 B2 JP4747336 B2 JP 4747336B2
Authority
JP
Japan
Prior art keywords
variable
value
damping
damper
control
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2005102494A
Other languages
Japanese (ja)
Other versions
JP2006283826A (en
Inventor
和夫 吉田
良和 北川
一郎 長島
洋三 篠崎
隆英 小林
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.)
Taisei Corp
Keio University
Hitachi Astemo Ltd
Original Assignee
Taisei Corp
Keio University
Hitachi Automotive Systems Ltd
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 Taisei Corp, Keio University, Hitachi Automotive Systems Ltd filed Critical Taisei Corp
Priority to JP2005102494A priority Critical patent/JP4747336B2/en
Publication of JP2006283826A publication Critical patent/JP2006283826A/en
Application granted granted Critical
Publication of JP4747336B2 publication Critical patent/JP4747336B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Vibration Prevention Devices (AREA)
  • Fluid-Damping Devices (AREA)
  • Buildings Adapted To Withstand Abnormal External Influences (AREA)

Description

本発明は、免震構造物に設けられた可変減衰ダンパの減衰力を制御することにより、地震時における上部構造の応答を低減化させる可変構造セミアクティブ免震システムに関するものである。   The present invention relates to a variable structure semi-active seismic isolation system that reduces the response of an upper structure during an earthquake by controlling the damping force of a variable damping damper provided in the base isolation structure.

近年、地震に対して構造物の安全性を確保するために、当該構造物の基礎部分や中間階に積層ゴムや滑り支承による免震装置を介装して、地震等によって地盤から構造物に伝播しようとする振動を緩和させて、構造物の躯体に生じる応力や変形を少なくするとともに、さらに上下部構造間に、粘弾性ダンパやオイルダンパー等の減衰ダンパを介装して、上記振動を積極的に減衰させる各種のパッシブ免震システムが適用されている。   In recent years, in order to ensure the safety of a structure against an earthquake, the base part and intermediate floor of the structure are provided with a base rubber and a base isolation device using a sliding bearing, and the structure can be changed from the ground to the structure by an earthquake. The vibration to be propagated is alleviated to reduce the stress and deformation generated in the frame of the structure. Further, a damping damper such as a viscoelastic damper or an oil damper is interposed between the upper and lower structures to suppress the vibration. Various passive seismic isolation systems that actively attenuate are applied.

しかしながら、このようなパッシブ免震システムにおいては、減衰ダンパにおける減衰力が一定であるために、様々なレベルの地震動に対して上記構造物の下部あるいは上部に生じる変位、速度あるいは加速度といった応答を最適に制御することができない。   However, in such a passive seismic isolation system, since the damping force in the damping damper is constant, the response such as displacement, velocity, or acceleration generated at the lower or upper part of the structure is optimal for various levels of ground motion. Can not be controlled.

そこで、上記減衰ダンパとして、その減衰力が可変な可変減衰ダンパを用いるとともに、速度センサによって上記可変減衰ダンパーの速度を計測し、得られた検出信号からコンピュータにより最適制御理論などに基づいて最適化した制御力指令値Uを求めて、上記可変減衰ダンパで実現可能な減衰係数値C等を選択して当該可変減衰ダンパを制御するセミアクティブ免震システムが開発されている。   Therefore, as the damping damper, a variable damping damper whose damping force is variable is used, and the speed of the variable damping damper is measured by a speed sensor, and the obtained detection signal is optimized based on the optimal control theory by a computer. A semi-active seismic isolation system has been developed in which the control force command value U is obtained, the damping coefficient value C that can be realized by the variable damping damper is selected, and the variable damping damper is controlled.

ところで、このようなセミアクティブ免震システムにおいては、地震規模が小さい場合には、可変減衰ダンパの減衰力を小さく設定することにより、高い免震効果を発揮させることができる反面、免震構造物における変位は大きくなる。他方、地震規模が大きい場合に、可変減衰ダンパの減衰力を大きく設定すれば、免震効果は小さくなるものの、上記変位を小さくすることができる。   By the way, in such a semi-active seismic isolation system, if the magnitude of the earthquake is small, it is possible to exert a high seismic isolation effect by setting the damping force of the variable damping damper small. The displacement at becomes larger. On the other hand, if the magnitude of the earthquake is large and the damping force of the variable damping damper is set large, the above-mentioned displacement can be reduced although the seismic isolation effect is reduced.

さらに、特に高層免震建物においては、振動数が高い領域において、その2次以上の高次振動モードによって上部の加速度が大きくなる場合がある。
このため、上記セミアクティブ免震システムにおいては、可変減衰ダンパ自体が奏する制振効果とともに、上記特性を勘案して、構造物に生じる加速度が全体として小さくなるようなきめ細かな制御を行うことが重要となる。
Further, particularly in a high-rise base-isolated building, the upper acceleration may be increased by the higher-order vibration mode of the second or higher order in the high-frequency region.
For this reason, in the semi-active seismic isolation system, it is important to perform fine control so that the acceleration generated in the structure is reduced as a whole in consideration of the above characteristics in addition to the vibration damping effect produced by the variable damping damper itself. It becomes.

しかしながら、上述した従来のセミアクティブ免震システムにあっては、先ず最適制御理論に基づいて制御力指令値Uを求め、当該指令値Uに基づいて可変減衰ダンパにおける実現可能な減衰係数値C等を選択する方法を採っているために、制御性能のロバスト性が悪いという問題点があった。   However, in the above-described conventional semi-active seismic isolation system, first, the control force command value U is obtained based on the optimal control theory, and the damping coefficient value C that can be realized in the variable damping damper based on the command value U. Since the method of selecting is adopted, there is a problem that the robustness of the control performance is poor.

このような問題点を解決する可変減衰ダンパ等の制御理論として、下記非特許文献1に開示されているような、双線形最適制御理論が定式化されている。
この双線形最適制御理論によれば、上述した構造物の応答検出信号から、直接可変減衰ダンパにおける減衰係数を最適化することができ、よって上記最適制御理論を用いたシステムの制御性能を大幅に改善することができるという利点がある。
As a control theory for a variable damping damper or the like that solves such a problem, a bilinear optimal control theory as disclosed in Non-Patent Document 1 below is formulated.
According to this bilinear optimal control theory, it is possible to optimize the damping coefficient in the variable damping damper directly from the response detection signal of the structure described above, and thus greatly improve the control performance of the system using the optimal control theory. There is an advantage that it can be improved.

しかしながら、上記双線形最適制御理論を用いた免震システムにおいても、構造物および可変減衰ダンパをモデル化する際に生じる誤差により、制御性能のロバスト性が低下してしまうという問題点を有している。   However, even in the seismic isolation system using the bilinear optimal control theory, there is a problem that the robustness of the control performance is reduced due to an error generated when modeling the structure and the variable damping damper. Yes.

また、いずれの制御理論によっても、可変減衰ダンパに対する減衰力設定値と、当該可変減衰ダンパで実際に発生する減衰力との間に生じる遅れが無視されることになるために、高い精度で上記可変減衰ダンパの最適な減衰力を制御することが難しいという問題点がある。
「吉田和夫、藤尾忠洋、双線形最適制御理論とセミアクティブ免震への応用」、日本機械学会論文集(C編)、67巻656号、2001.4.第96〜102頁
In addition, in any control theory, the delay between the damping force set value for the variable damping damper and the damping force actually generated by the variable damping damper is ignored, so the above-mentioned with high accuracy. There is a problem that it is difficult to control the optimum damping force of the variable damping damper.
“Kazuo Yoshida, Tadahiro Fujio, Bilinear Optimal Control Theory and Application to Semi-active Seismic Isolation”, Transactions of the Japan Society of Mechanical Engineers (C), Vol. 67, No. 656, 2001.4. Pages 96-102

この発明は、かかる事情に鑑みてなされたもので、可変減衰ダンパにおける遅れも考慮することができるとともに、さらに制御対象となる構造物や可変減衰ダンパのモデル化誤差に対する制御性能のロバスト性も一段と高めることができる可変構造セミアクティブ免震システムを提供することを課題とするものである。   The present invention has been made in view of such circumstances, and can take into account the delay in the variable damping damper, and further enhance the robustness of the control performance against the modeling error of the structure to be controlled and the variable damping damper. It is an object to provide a variable structure semi-active seismic isolation system that can be enhanced.

上記課題を解決するために、請求項1に記載の発明は、上下部構造物間に免震装置が介装された免震構造物の上記上部構造と上記下部構造との間に設けられた可変減衰ダンパと、地震時における上記上部構造および下部構造の応答を検出するセンサと、このセンサからの検出信号に基づいて上記可変減衰ダンパへの指令信号を演算・出力するコンピュータを備えた制御手段と、この制御手段から出力された上記指令信号により上記可変減衰ダンパによる減衰力を切り換える切換手段とを備えた可変構造セミアクティブ免震システムにおいて、上記コンピュータが、上記免震構造物自体の状態方程式に、上記可変減衰ダンパの減衰係数設定値とダンパ速度との積で規定される減衰力設定値と上記可変減衰ダンパで実際に発生する減衰力との一次遅れ系の関係式であって、かつ時定数を一定と近似した関係式を連成させた、全体システムの状態方程式に対して双線形最適制御理論を適用して求めた制御則を切換え超平面として設定し、上記切換え超平面に基づいて上記可変減衰ダンパの減衰係数を切り換えて上記上部構造物の応答変形および/または応答加速度を制御する可変構造制御法を採用して上記指令信号を出力することを特徴とするものである。 In order to solve the above problem, the invention according to claim 1 is provided between the upper structure and the lower structure of the seismic isolation structure in which the seismic isolation device is interposed between the upper and lower structures. Control means comprising a variable damping damper, a sensor for detecting the response of the upper structure and the lower structure during an earthquake, and a computer for calculating and outputting a command signal to the variable damping damper based on a detection signal from the sensor And a switching means for switching the damping force by the variable damping damper according to the command signal output from the control means, wherein the computer has a state equation of the seismic isolation structure itself. , the primary and the damping force actually generated by the variable damping attenuation coefficient setting value of the damper and the damping force setting value defined by the product of the damper speed and the variable attenuation damper It is based a relational expression, and the time constant was formed continuously the relational expression approximate constant hyperplane switched control law determined by applying the bilinear optimal control theory with respect to the equation of state of the entire system The command signal is output by adopting a variable structure control method in which the damping coefficient of the variable damping damper is switched based on the switching hyperplane to control response deformation and / or acceleration of the upper structure. It is characterized by this.

また、請求項に記載の発明は、請求項に記載の上記コンピュータが、予め上記免震構造物の低応答値および高応答値並びに上記低応答値よりも低い制御開始を判断する応答値が設定されており、上記センサからの上記検出信号が上記制御開始を判断する応答値未満である場合に上記減衰力を最大値に設定し、上記制御開始を判断する応答値を超えて上記低応答値未満である場合に上記減衰力を最小値に設定し、上記検出信号が上記低応答値を超えた以降は、上記高応答値以下の範囲である場合に、上記演算を実行し、上記高応答値を超えた場合に上記減衰力を最大値に設定することを特徴とするものである。 In the invention described in claim 2 , the computer according to claim 1 determines in advance a low response value and a high response value of the seismic isolation structure, and a control value lower than the low response value. Is set, the damping force is set to the maximum value when the detection signal from the sensor is less than the response value for determining the start of control, and the low value exceeds the response value for determining the start of control. When the damping force is less than the response value, the damping force is set to the minimum value, and after the detection signal exceeds the low response value, the calculation is performed when the range is equal to or less than the high response value. When the high response value is exceeded, the damping force is set to a maximum value.

さらに、請求項に記載の発明は、請求項1または2に記載の上記可変減衰ダンパが、上記減衰力が低減衰力と高減衰力との2段階に可変であることを特徴とするものである。 Furthermore, the invention according to claim 3 is characterized in that the variable damping damper according to claim 1 or 2 is characterized in that the damping force is variable in two stages of a low damping force and a high damping force. It is.

また、請求項に記載の発明は、請求項1〜のいずれかに記載の上記コンピュータが、上記上部構造と下部構造との相対速度または上記可変減衰ダンパの応答速度が規定値以下である場合に、上記指令信号を更新しない、すなわち減衰力を切り換えないことを特徴とするものである。 According to a fourth aspect of the present invention, in the computer according to any one of the first to third aspects, a relative speed between the upper structure and the lower structure or a response speed of the variable damping damper is less than a specified value. In this case, the command signal is not updated, that is, the damping force is not switched.

請求項に記載の発明は、請求項1〜のいずれかに記載の発明において、上記可変減衰ダンパが、内部に流体が充填された圧力室を有するシリンダと、このシリンダ内に移動自在に設けられて上記圧力室を2室に仕切るピストンと、このピストンと一体化されて上記シリンダから延出するロッドと、上記2室間を連通させる連通経路とを備えた流体圧ダンパであり、かつ上記切換手段が、上記連通経路に配設されて当該連通経路を開閉させることにより上記減衰力を切り換える電磁弁であるとともに、上記制御手段が、上記電磁弁に供給される駆動電流値を検出する検出手段と、この検出手段により検出された電流値が異常か否かを判定する判定手段を備えることを特徴とするものである。 According to a fifth aspect of the present invention, in the invention according to any one of the first to fourth aspects, the variable damping damper includes a cylinder having a pressure chamber filled with a fluid therein, and is movable in the cylinder. A fluid pressure damper provided with a piston provided to partition the pressure chamber into two chambers, a rod integrated with the piston and extending from the cylinder, and a communication path communicating between the two chambers; The switching means is an electromagnetic valve that is disposed in the communication path and switches the damping force by opening and closing the communication path, and the control means detects a drive current value supplied to the electromagnetic valve. It is characterized by comprising detection means and determination means for determining whether or not the current value detected by the detection means is abnormal.

また、請求項に記載の発明は、請求項に記載の判定手段が、上記電磁弁を駆動する上記指令信号に対する上記駆動電流値がしきい未満の場合、および上記電磁弁の駆動を停止する上記指令信号に対する上記駆動電流値がしきい値以上の場合に上記異常と判断する判断手段と、この判断手段が異常と判断した上記電磁弁の駆動電流を遮断する電流遮断手段とを備えることを特徴とするものである。 Further, the invention according to claim 6, the determination means according to claim 5, the drive current value with respect to the command signal for driving the electromagnetic valve is of less than the threshold value, and the driving of the solenoid valve A determination unit configured to determine that the abnormality is present when the drive current value with respect to the command signal to be stopped is equal to or greater than a threshold; and a current blocking unit configured to block the drive current of the electromagnetic valve that is determined to be abnormal by the determination unit. It is characterized by this.

請求項1〜のいずれかに記載の発明においては、制御対象となる免震構造物の状態方程式に、さらに可変減衰ダンパの減衰力設定値と可変減衰ダンパで実際に発生する減衰力との関係式を連成させた全体システムの状態方程式に対して、双線形最適制御理論を適用しているので、減衰係数の切り換えに対する実際の可変減衰ダンパにおける遅れも考慮した制御を行うことができる。 In the invention according to any one of claims 1 to 6 , the state equation of the seismic isolation structure to be controlled includes the damping force setting value of the variable damping damper and the damping force actually generated by the variable damping damper. Since the bilinear optimal control theory is applied to the state equation of the whole system in which the relational expressions are coupled, it is possible to perform control in consideration of the delay in the actual variable damping damper with respect to switching of the damping coefficient.

加えて、上記双線形最適制御理論を適用して設定した切換え超平面に基づいて、上記可変減衰ダンパの減衰係数を切り換えるスライディングモード制御法等の可変構造制御法を採用しているために、モデル化に際して生じる誤差を補償してロバスト性を一段と高めることが可能になる。   In addition, a variable structure control method such as a sliding mode control method for switching the damping coefficient of the variable damping damper based on the switching hyperplane set by applying the bilinear optimal control theory is used. It is possible to further enhance the robustness by compensating for an error that occurs during the conversion.

この際に、上記減衰力設定値と可変減衰ダンパで実際に発生する減衰力の関係式として、一次遅れ系を用いるとともに、当該関係式における時定数を一定と近似すれば、双線形最適制御理論を適用して切換え超平面を設定する演算を行うに際して、上記関係式を連成した全体システムの状態方程式の取り扱いが容易となる。 At this time, as relation of the damping force actually generated by the damping force setting value and the variable attenuation damper with use of first-order lag system, lever approximating a constant time constant in the equation, bilinear optimal control When the calculation is performed to set the switching hyperplane by applying theory, it becomes easy to handle the state equation of the entire system coupled with the above relational expressions.

また、通常の風や周辺環境から伝わる一時的な小振動等に起因して免震構造物に小さな振動が生じた場合には、上記可変減衰ダンパに対する細かな制御を行う必要が無く、他方大地震時に極端に大きな振動が作用した場合には、上記制御を行っても、得られた所望の減衰係数が可変減衰ダンパにおいて調整可能な最大値を超えてしまうことになる。   If small vibrations occur in the seismic isolation structure due to normal wind or temporary small vibrations transmitted from the surrounding environment, it is not necessary to perform fine control on the variable damping damper. If extremely large vibrations are applied during an earthquake, the desired damping coefficient obtained will exceed the maximum value that can be adjusted by the variable damping damper even if the above control is performed.

そこで、請求項に記載の発明においては、予め上記免震構造物の低応答値と高応答値並びに制御開始を判断する応答値を設定しておき、上記制御開始を判断する応答値未満である場合に上記減衰力を最大値に固定し、上記制御開始を判断する応答値を超えて上記低応答値未満である場合に上記減衰力を最小値に固定し、上記低応答値を超えた以後は、上記高応答値以下の範囲である場合にのみ制御を行い、上記高応答値を超えた場合は上記減衰力を最大値に固定して、無用な制御を回避するようにしたものである。 Therefore, in the invention described in claim 2 , a low response value and a high response value of the seismic isolation structure and a response value for determining the start of control are set in advance, and less than the response value for determining the start of control. In some cases, the damping force is fixed to the maximum value, and when the response value for determining the start of control exceeds the low response value, the damping force is fixed to the minimum value, and the low response value is exceeded. Thereafter, control is performed only when the value is within the range of the high response value, and when the high response value is exceeded, the damping force is fixed at the maximum value to avoid unnecessary control. is there.

さらに、本発明者等の解析および実験によれば、可変減衰ダンパの減衰力を2段階に切り換え可能とするのみで、十分に細かな制御を実現することが可能であり、それ以上の複数段に切り換えても制御性能に大幅な向上が見られないことが判明している。そこで、請求項に記載の発明のように、可変減衰ダンパとして、上記減衰力が2段階切り換えのものを用いれば、構造が簡易になるとともに、コンピュータにおける演算・制御が容易になる。 Furthermore, according to the analysis and experiment by the present inventors, it is possible to realize sufficiently fine control only by making it possible to switch the damping force of the variable damping damper to two stages, and a plurality of more stages. It has been found that there is no significant improvement in control performance even when switching to. Therefore, as in the invention described in claim 3 , if a variable damping damper having a two-step switching damping force is used, the structure becomes simple and the calculation and control in the computer becomes easy.

また、上記免震構造物においては、平常時においても、周囲から伝搬する微振動によって、可変減衰ダンパに小さな速度の動きが生じる。例えば、上記可変減衰ダンパとして、オイルダンパーを用いた場合には、上記微振動により例えば0.1kine程度の速度の動きが常時生じてしまう。このような場合には、上記免震システムとしての制御を行うことが必要でないために、請求項に記載の発明のように、上下部構造間の相対速度あるいは可変減衰ダンパの応答速度が規定値以下、例えば1kine以下であるような場合には、当該可変減衰ダンパに対する上記制御を行わないことが好ましい。 Further, in the above seismic isolation structure, even at normal times, a small speed movement occurs in the variable damping damper due to fine vibrations propagating from the surroundings. For example, when an oil damper is used as the variable damping damper, a movement with a speed of, for example, about 0.1 kine always occurs due to the fine vibration. In such a case, since it is not necessary to perform control as the seismic isolation system, the relative speed between the upper and lower structures or the response speed of the variable damping damper is defined as in the invention described in claim 4. When the value is less than the value, for example, less than 1 kine, it is preferable not to perform the above control on the variable damping damper.

さらに、請求項に記載の発明によれば、上記可変減衰ダンパとして、オイルダンパ等の流体圧ダンパを用いた場合に、減衰力の切換手段となる電磁弁の駆動電流値の判定手段を設ければ、どの電磁弁が異常か否かを監視することができ、よって制御時のシステムの信頼性を高められるとともに、メンテナンス作業の手間を大幅に緩和することができる。 Further, according to the invention described in claim 5, when a fluid pressure damper such as an oil damper is used as the variable damping damper, a means for determining a driving current value of an electromagnetic valve serving as a damping force switching means is provided. As a result, it is possible to monitor which solenoid valve is abnormal. Therefore, the reliability of the system at the time of control can be improved, and the labor of maintenance work can be greatly reduced.

この際に、請求項に記載の発明のように、上記判定手段として、予め電磁弁に供給される駆動電流値のしきい値を設定しておき、当該しきい値と指令信号に対する実際の駆動電流値とを対比することにより、容易かつ正確に上記電磁弁の異常を判断することが可能になる。 At this time, as in the invention described in claim 6 , as the determination means, a threshold value of the drive current value supplied to the solenoid valve is set in advance, and an actual value for the threshold value and the command signal is set. By comparing with the drive current value, it is possible to easily and accurately determine the abnormality of the electromagnetic valve.

図1〜図18は、本発明の可変構造セミアクティブ免震システム(以下、単に免震システムとも称す。)の一実施形態を示すものである。
図1〜図3に示すように、この免震システムは、基礎(下部構造)1と建物(上部構造)2との間に、積層ゴム支承による免震装置3が介装された免震構造物を対象とするものであり、基礎1と建物2との間に設けられた可変減衰オイルダンパ(可変減衰ダンパ)4および減衰係数固定型オイルダンパ4´と、基礎1、免震層の建物2側および建物2の頂部に設けられて地震時における加速度を検出するセンサ5a、5b、5cおよび変位センサ5dと、これらセンサ5a、5b、5c、5dからの検出信号に基づいてオイルダンパ4への指令信号を演算・出力する制御手段6と、この制御手段6から出力された指令信号によりオイルダンパ4の減衰力を切り換える電磁弁7とから概略構成されている。
FIGS. 1-18 shows one Embodiment of the variable structure semi-active seismic isolation system (henceforth only a seismic isolation system) of this invention.
As shown in FIG. 1 to FIG. 3, this seismic isolation system has a base isolation structure in which a base rubber (lower structure) 1 and a building (upper structure) 2 are provided with a base rubber isolation device 3 with a laminated rubber support. A variable damping oil damper (variable damping damper) 4 and a fixed damping coefficient type oil damper 4 'provided between a foundation 1 and a building 2, a foundation 1, a base isolation layer building Sensors 5a, 5b, 5c and a displacement sensor 5d that are provided on the second side and on the top of the building 2 to detect acceleration during an earthquake, and to the oil damper 4 based on detection signals from these sensors 5a, 5b, 5c, 5d The control means 6 for calculating / outputting the command signal and the electromagnetic valve 7 for switching the damping force of the oil damper 4 by the command signal output from the control means 6 are generally configured.

ここで、制御手段6は、上記演算を行うとともにその結果を指令信号として出力するコンピュータ8を有する制御回路12と、この制御回路12からの指令信号に基づいて電磁弁7へ駆動電流を供給する電源回路9とを備えている。   Here, the control means 6 supplies the drive current to the electromagnetic valve 7 on the basis of the control circuit 12 having the computer 8 that performs the above calculation and outputs the result as a command signal, and the command signal from the control circuit 12. And a power supply circuit 9.

さらに、この電源回路9には、制御回路12から出力された指令信号により励磁されてオン状態に切り替わるリレー10と、リレー10から供給される駆動電流を測定して電流の変化を検出する交流電流センサ(検出手段)11とを有している。ここで、リレー10としては、例えば半導体を用いた無接点なソリッドステートリレー(SSR)などが用いられている。   Further, the power supply circuit 9 includes a relay 10 that is excited by a command signal output from the control circuit 12 and switches to an on state, and an alternating current that detects a change in current by measuring a drive current supplied from the relay 10. And a sensor (detection means) 11. Here, as the relay 10, for example, a non-contact solid state relay (SSR) using a semiconductor is used.

交流電流センサ11は、電磁弁に供給される駆動電流値を検出するためのものであり、その検出信号が制御手段6に出力されるようになっている。他方、制御手段6のコンピュータ8には、交流電流センサ11によって検出された電磁弁7に供給される電流値が異常値か否かを判定する制御プログラム(判定手段)が格納されている。   The alternating current sensor 11 is for detecting the drive current value supplied to the electromagnetic valve, and the detection signal is output to the control means 6. On the other hand, the computer 8 of the control means 6 stores a control program (determination means) for determining whether or not the current value supplied to the electromagnetic valve 7 detected by the alternating current sensor 11 is an abnormal value.

また、可変減衰オイルダンパ4は、図4〜図6に示すように、内部の圧力室15に作動油が充填されたシリンダ16と、シリンダ16の他端から圧力室15に挿入されたピストン17と、一端がピストン17に結合されたピストンロッド18と、ピストン17を貫通する液通路19に設けられ圧力室15の右室15bから左室15aへの作動油の流れのみを許容する逆止弁20と、作動油が充填された補助タンク21と、補助タンク21と圧力室15との間を連通する連通経路22とから概略構成されており、連通経路22に交流電流で駆動する上記電磁弁7が配設されている。この電磁弁7は、電源回路9から供給される交流の駆動電流がオンになると開弁し、駆動電流がオフになると閉弁する構成である。   As shown in FIGS. 4 to 6, the variable damping oil damper 4 includes a cylinder 16 in which the internal pressure chamber 15 is filled with hydraulic oil, and a piston 17 inserted into the pressure chamber 15 from the other end of the cylinder 16. A check rod which is provided in a piston rod 18 having one end coupled to the piston 17 and a liquid passage 19 penetrating the piston 17 and allows only the flow of hydraulic oil from the right chamber 15b of the pressure chamber 15 to the left chamber 15a. 20, an auxiliary tank 21 filled with hydraulic oil, and a communication path 22 that communicates between the auxiliary tank 21 and the pressure chamber 15, and the solenoid valve that is driven by an alternating current in the communication path 22. 7 is disposed. The electromagnetic valve 7 is configured to open when an alternating drive current supplied from the power supply circuit 9 is turned on and close when the drive current is turned off.

そして、ピストンロッド18のXa方向の端部は、回動可能な継ぎ手23を介して建物2に連結されている。また、シリンダ16のXb方向の端部は、回動可能な継ぎ手24を介して基礎1に連結されている。この結果、オイルダンパ4は、基礎1に地震による振動が伝播されると、シリンダ16とピストン17との間で伸縮動作が行われて振動を減衰させる減衰力を発生させることにより、建物2に伝播される振動エネルギを減少させるようになっている。   The end of the piston rod 18 in the Xa direction is connected to the building 2 via a rotatable joint 23. Further, the end of the cylinder 16 in the Xb direction is connected to the foundation 1 via a rotatable joint 24. As a result, when the vibration due to the earthquake is propagated to the foundation 1, the oil damper 4 is expanded and contracted between the cylinder 16 and the piston 17 to generate a damping force that attenuates the vibration. The transmitted vibration energy is reduced.

また、補助タンク21の内部は、作動油が貯留される油層と空気層が形成されたリザーバ室21aになっており、連通経路22は、シリンダ16内に形成された筒状空間からなり、圧力室15の外周に形成されている。この圧力室15の右室15bは、吸込弁25を介して連通経路22に連通されており、ピストン17がXa方向に移動する過程で吸込弁25が開弁することにより、シリンダ15から押し出されたピストンロッド18の体積分の作動油が補助タンク21から連通経路22を介して右室15bに供給される。また、ピストン17がXb方向に移動する過程では、吸込弁25が閉弁し、逆止弁20が開弁することにより、右室15bの作動油がXb方向に押し込まれたピストンロッド18の体積分液通路19を通過して左室15aへ流入するようになっている。   The auxiliary tank 21 has a reservoir chamber 21a in which an oil layer and an air layer for storing hydraulic oil are formed, and the communication path 22 is formed of a cylindrical space formed in the cylinder 16 and has a pressure. It is formed on the outer periphery of the chamber 15. The right chamber 15b of the pressure chamber 15 is connected to the communication path 22 via the suction valve 25, and is pushed out of the cylinder 15 by opening the suction valve 25 in the process of moving the piston 17 in the Xa direction. The hydraulic oil corresponding to the volume of the piston rod 18 is supplied from the auxiliary tank 21 to the right chamber 15 b via the communication path 22. Further, in the process of moving the piston 17 in the Xb direction, the volume of the piston rod 18 in which the hydraulic oil in the right chamber 15b is pushed in the Xb direction by closing the suction valve 25 and opening the check valve 20. It passes through the liquid separation passage 19 and flows into the left chamber 15a.

さらに、圧力室15の図中左方には、左室15aに連通された通路26を開閉する第1の調圧弁27と、左室15aに連通された通路28を開閉する複数(本実施例では、4個)の上記電磁弁7および複数(本実施例では、2個)の第2の調圧弁29と、連通経路22に連通された複数(本実施例では、5個)のリリーフ弁30とが配設されている。
ここで、第1の調圧弁27、第2の調圧弁29およびリリーフ弁30は、例えば、弁体をコイルバネにより付勢する構成であり、コイルバネのバネ定数を任意の値に選択することにより、開弁圧力を任意の値に設定することができる。リリーフ弁30は、第1の調圧弁27よりも開弁圧力が高くなるようにバネ定数が設定されており、ピストン速度が所定以上に高くなったときに開弁するように設定されている。
Further, on the left side of the pressure chamber 15 in the figure, a first pressure regulating valve 27 that opens and closes a passage 26 that communicates with the left chamber 15a, and a plurality of (this embodiment) that opens and closes a passage 28 that communicates with the left chamber 15a. Then, four solenoid valves 7 and a plurality (two in this embodiment) second pressure regulating valves 29 and a plurality (five in this embodiment) relief valves communicated with the communication path 22. 30.
Here, the 1st pressure regulation valve 27, the 2nd pressure regulation valve 29, and the relief valve 30 are the structures which energize a valve body with a coil spring, for example, By selecting the spring constant of a coil spring to arbitrary values, The valve opening pressure can be set to an arbitrary value. The relief valve 30 has a spring constant set so that the valve opening pressure is higher than that of the first pressure regulating valve 27, and is set to open when the piston speed becomes higher than a predetermined value.

また、電磁弁7と第2の調圧弁29は、通路28に直列に配置されており、電磁弁7が開弁したときに第2の調圧弁29に作動油の圧力が作用するように構成されている。
ピストン17がXa方向に移動する動作過程では、逆止弁20が閉弁しており、シリンダ16から押し出されたピストンロッド18の体積分の作動油が左室15aから第1の調圧弁27および連通経路22を介して補助タンク21へ吐出される。この場合、作動油が第1の調圧弁27を通過する過程で抵抗が付与されて減衰力が発生する。さらに、ピストン17の速度が加速された場合には、複数のリリーフ弁30が開弁して左室15aから吐出された作動油を連通経路22へ逃がすようになっている。
Further, the solenoid valve 7 and the second pressure regulating valve 29 are arranged in series with the passage 28, and are configured such that the hydraulic pressure acts on the second pressure regulating valve 29 when the solenoid valve 7 is opened. Has been.
In the operation process in which the piston 17 moves in the Xa direction, the check valve 20 is closed, and the hydraulic oil corresponding to the volume of the piston rod 18 pushed out from the cylinder 16 flows from the left chamber 15a to the first pressure regulating valve 27 and It is discharged to the auxiliary tank 21 through the communication path 22. In this case, resistance is applied in the process of the hydraulic oil passing through the first pressure regulating valve 27, and a damping force is generated. Further, when the speed of the piston 17 is accelerated, the plurality of relief valves 30 are opened so that the hydraulic oil discharged from the left chamber 15 a is released to the communication path 22.

また、制御回路12からの指示により電磁弁7に駆動電流が供給されて電磁弁7が開弁しているときは、第2の調圧弁29および連通経路22を介して補助タンク21へ吐出されるようになっている。   Further, when a drive current is supplied to the electromagnetic valve 7 according to an instruction from the control circuit 12 and the electromagnetic valve 7 is opened, the electromagnetic valve 7 is discharged to the auxiliary tank 21 via the second pressure regulating valve 29 and the communication path 22. It has become so.

ここで、以上の構成からなるオイルダンパ4におけるピストン動作に伴う減衰力の減衰特性について説明する。
図7はピストン速度に対する減衰係数の変化を示すグラフである。図7中、グラフIは、ピストン速度が比較的小さい段階で減衰力が発生するように減衰係数を変化させた場合を示しており、グラフIIはピストン速度が比較的大きい段階で減衰力が発生するように減衰係数を変化させた場合を示している。すなわち、このオイルダンパ4では、後述するように基礎1から入力された振動の加速度が比較的小さく、振幅が大きい場合には、グラフIで示すような減衰特性を有し、基礎1から入力された振動の加速度及び振幅が大きい場合には、グラフIIで示すような減衰特性を有する。
Here, the damping characteristic of the damping force accompanying the piston operation in the oil damper 4 having the above configuration will be described.
FIG. 7 is a graph showing changes in the damping coefficient with respect to the piston speed. In FIG. 7, graph I shows the case where the damping coefficient is changed so that the damping force is generated when the piston speed is relatively low, and graph II shows that the damping force is generated when the piston speed is relatively high. The case where the attenuation coefficient is changed is shown. That is, the oil damper 4 has a damping characteristic as shown in the graph I when the acceleration of vibration input from the foundation 1 is relatively small and the amplitude is large as will be described later. When the acceleration and amplitude of the vibration are large, it has a damping characteristic as shown in graph II.

そして、このオイルダンパ4においては、上記グラフI、IIの勾配の差違から減衰係数がC1H、C2H、C1L、C2Lの4種類に設定されている。そして、制御回路12のコンピュータ8は、後述するように加速度センサ5a、5b、5cから出力された検出信号に基づいて建物2の振動を抑制するための減衰力を求め、図7に示すグラフI、IIに示す減衰力に対応する減衰係数が得られるように電磁弁7を開弁または閉弁させる。   And in this oil damper 4, the attenuation coefficient is set to four types, C1H, C2H, C1L, and C2L, from the difference of the gradient of the said graphs I and II. Then, the computer 8 of the control circuit 12 obtains a damping force for suppressing the vibration of the building 2 based on the detection signals output from the acceleration sensors 5a, 5b, and 5c as will be described later, and the graph I shown in FIG. The electromagnetic valve 7 is opened or closed so that a damping coefficient corresponding to the damping force shown in FIG.

すなわち、電磁弁7を閉弁させた状態では、第1の調圧弁27が開弁し、電磁弁7を開弁させた状態では、第1の調圧弁27および第2の調圧弁29が開弁して減衰係数を切り替える。
図8に示されるように、減衰係数C1Hの状態では、第1の調圧弁27が開弁しており、図中実線矢印で示すように、左室15aからリザーバ室21aへ流出される作動油が、第1の調圧弁27の開弁によって絞られることになり、図7に示されるC1H領域の減衰力が発生される。この際に、作動油の圧力が上昇すると、さらにリリーフ弁30が開弁することにより、図中点線矢印で示すように作動油が流れ、減衰係数C2H領域の減衰力が発生される。
That is, when the solenoid valve 7 is closed, the first pressure regulating valve 27 is opened, and when the solenoid valve 7 is opened, the first pressure regulating valve 27 and the second pressure regulating valve 29 are opened. To change the damping coefficient.
As shown in FIG. 8, in the state of the damping coefficient C1H, the first pressure regulating valve 27 is open, and the hydraulic fluid that flows out from the left chamber 15a to the reservoir chamber 21a as shown by the solid line arrow in the figure. However, the first pressure regulating valve 27 is throttled, and the damping force in the C1H region shown in FIG. 7 is generated. At this time, when the pressure of the hydraulic oil rises, the relief valve 30 is further opened, so that the hydraulic oil flows as shown by a dotted arrow in the figure, and a damping force in the damping coefficient C2H region is generated.

これに対して、図9に示されるように、電磁弁7を開弁させると、第1の調圧弁27および電磁弁7、第2の調圧弁29が開弁し、図中実線矢印で示すように、左室15aからリザーバ室21aへ流出される作動油が、第1の調圧弁27および第2の調圧弁29の開弁によって絞られることになり、図7に示されるC1L領域の減衰力が発生される。また、同様にこの状態から、さらに作動油の圧力が上昇すると、リリーフ弁30が開弁することにより、図中点線矢印で示すように作動油が流れ、減衰係数C2L領域の減衰力が発生される。   On the other hand, as shown in FIG. 9, when the solenoid valve 7 is opened, the first pressure regulating valve 27, the solenoid valve 7, and the second pressure regulating valve 29 are opened, which are indicated by solid arrows in the figure. As described above, the hydraulic oil flowing out from the left chamber 15a to the reservoir chamber 21a is throttled by opening the first pressure regulating valve 27 and the second pressure regulating valve 29, and the damping of the C1L region shown in FIG. Force is generated. Similarly, when the hydraulic oil pressure further increases from this state, the relief valve 30 opens, and the hydraulic oil flows as shown by the dotted arrow in the figure, and a damping force in the damping coefficient C2L region is generated. The

そこで次に、このようなオイルダンパ4における減衰力の切換を行わせるために、上記コンピュータ8において実行する制御方法について説明する。
(1)オイルダンパ4の減衰係数設定値cとダンパ速度の積で規定される減衰力設定値と可変減衰ダンパで実際に発生する減衰力fの下記関係式の設定

Figure 0004747336
上記関係式は、実験的な方法またはMaxwellモデルに基づく近似により設定することができる。 Then, next, the control method performed in the said computer 8 in order to perform switching of the damping force in such an oil damper 4 is demonstrated.
(1) Setting of the following relational expression of the damping force setting value defined by the product of the damping coefficient setting value c of the oil damper 4 and the damper speed and the damping force f d actually generated by the variable damping damper
Figure 0004747336
The above relational expression can be set by an experimental method or approximation based on the Maxwell model.

先ず、図10に示すMaxwellモデルに基づいて、1次遅れ系に基づく減衰力設定値とダンパ減衰力との関係式を評価すると、

Figure 0004747336

この(3)式は,ダンパのピストン全体の速度と減衰係数の積からなる減衰力設定値に対して,実際に発生するダンパ減衰力fが1次遅れとなる関係を表している。また、時定数Tは、減衰係数に依存している。 First, based on the Maxwell model shown in FIG. 10, when evaluating the relational expression between the damping force set value based on the first-order lag system and the damper damping force,
Figure 0004747336

The equation (3), with respect to the damping force setting value that is the product of velocity and attenuation coefficient of the entire piston damper, the damper damping force f d that actually generated represents the relationship: the first-order lag. In addition, the time constant T d depends on the attenuation coefficient.

この時定数Tを、CH=36.75kN・s/cmとCL=12.25kN・s/cmの2値に減衰係数が切換え可能な最大減衰力980kNの実大可変減衰ダンパでは、Maxwell部ばね定数の実測結果(k=294kN/cm)に対応した時定数は,以下の通りである。
CH=36.75kN・s/cm : T=0.125sec
CL=12.25kN・s/cm : T=0.042sec
なお、上記オイルダンパ4においては、変位制御三角波加振により時定数が図11に示す通り実測されている。
In the full-size variable damping damper having a maximum damping force of 980 kN, the time constant Td is changed to a binary value of CH = 36.75 kN · s / cm and CL = 12.25 kN · s / cm. The time constant corresponding to the actual measurement result (k = 294 kN / cm) of the spring constant is as follows.
CH = 36.75 kN · s / cm: T d = 0.125 sec
CL = 12.25 kN · s / cm: T d = 0.042 sec
In the oil damper 4, the time constant is actually measured as shown in FIG. 11 by displacement control triangular wave excitation.

上記Maxwellモデルに基づいて評価した時定数は、図11に示した実測結果に比べて若干小さめであるが、減衰係数が大きいほど時定数が大きい実測結果の傾向を良く表している。
そこで、本実施形態においては、時定数Tが一定と近似して、

Figure 0004747336

により双線形の状態方程式を作成する方法について定式化する。ここで、上記時定数Tとしては、CHとCLとに対応した時定数の平均値を使う方法や、最大値を使う方法を採用することができる。 Although the time constant evaluated based on the Maxwell model is slightly smaller than the actual measurement result shown in FIG. 11, the larger the damping coefficient, the better the tendency of the actual measurement result having a larger time constant.
Therefore, in this embodiment, the time constant T d is approximated to be constant,
Figure 0004747336

Formulate how to create a bilinear equation of state. Here, as the time constant Td , a method using an average value of time constants corresponding to CH and CL or a method using a maximum value can be employed.

(2)全体システムの状態方程式の作成
次に、(1)において得られた関係式を免震構造物自体の状態方程式に連成させた全体システムの状態方程式を作成する。
先ず、Maxwellモデルでモデル化したm台の減衰係数固定型オイルダンパ4´を備えた免震構造物を、x、y各方向を1台の可変減衰オイルダンパ4で制御する場合の状態方程式は、

Figure 0004747336
(2) Creation of State Equation of Entire System Next, the state equation of the entire system is created by coupling the relational expression obtained in (1) with the state equation of the seismic isolation structure itself.
First, the equation of state when a seismic isolation structure having m damping coefficient fixed type oil dampers 4 ′ modeled by the Maxwell model is controlled by one variable damping oil damper 4 in each of the x and y directions is as follows. ,
Figure 0004747336

ここで、

Figure 0004747336
here,
Figure 0004747336

なお、以下の定式化は、可変減衰オイルダンパ4を各方向1台であると想定しているが、可変減衰オイルダンパ4の台数については、Fの自由度を増やすことによって全く同様に定式化することができる。 Note that the following formulation is a variable damping oil damper 4 is assumed to be one in each direction, for the number of variable damping oil damper 4, just as in formulation by increasing the degree of freedom F d Can be

次いで、可変減衰オイルダンパ部を,減衰力の時間遅れを一定と近似する(4)式でモデル化すると、

Figure 0004747336
Next, when the variable damping oil damper part is modeled by the equation (4) that approximates the time delay of the damping force as follows,
Figure 0004747336

そこで、(5)式と(6)式とから、オイルダンパ4における実際の遅れを考慮した全体システムの状態方程式を作成すると、下記(7)式になる。

Figure 0004747336
Therefore, when the state equation of the entire system in consideration of the actual delay in the oil damper 4 is created from the equations (5) and (6), the following equation (7) is obtained.
Figure 0004747336

さらに、地震動の加速度は、白色雑音wを入力とする成形フィルタによって生成される有色雑音であると仮定する。ここで、成形フィルタとして、下記のバンドパス特性を有する伝達関数を想定すると、

Figure 0004747336
Further, it is assumed that the acceleration of the ground motion is colored noise generated by a shaping filter having white noise w as an input. Here, assuming a transfer function having the following bandpass characteristics as a shaping filter,
Figure 0004747336

ちなみに、図12は、成形フィルタの諸元としてω=100、ζ=0.707、ω=0.1、ζ=0.6と設定した場合の伝達関数を示すものである。
次いで、(8)式の伝達関数を状態方程式に書き直すと、

Figure 0004747336
Incidentally, FIG. 12 shows the transfer function when ω d = 100, ζ d = 0.707, ω h = 0.1, and ζ h = 0.6 as the dimensions of the shaping filter.
Next, when rewriting the transfer function of equation (8) into the state equation,
Figure 0004747336

x、y両方向に対して、例えば同じ成形フィルタを設定すると、

Figure 0004747336
For example, if the same shaping filter is set for both the x and y directions,
Figure 0004747336

ここで、w(t)は、スカラーの白色雑音として以下のような性質を持つ。ただし、E(*)は数学的期待値を、δ(*)はDiracのDelta関数を表す。
E[w(t)]=0、 E[w(t)w(τ)]=w(t)δ(t−τ)
Here, w (t) has the following properties as scalar white noise. Here, E (*) represents a mathematical expectation value, and δ (*) represents Dirac's Delta function.
E [w (t)] = 0, E [w (t) w (τ) T ] = w (t) δ (t−τ)

そして、(7)式および(9)式から、外乱を包含した拡大系の全体システムの状態方程式を設定すると以下の通りとなる。

Figure 0004747336
From the equations (7) and (9), the state equation of the entire system of the expanded system including the disturbance is set as follows.
Figure 0004747336

(10)式が、可変減衰オイルダンパ4の減衰係数設定値とダンパ速度の積で規定される減衰力設定値と当該可変減衰オイルダンパ4で実際に発生する減衰力の関係式と、制御対象となる上記免震構造物自体の状態方程式とを連成させた全体システムの状態方程式である。   (10) is a relational expression between the damping force setting value defined by the product of the damping coefficient setting value of the variable damping oil damper 4 and the damper speed, and the damping force actually generated in the variable damping oil damper 4, and the control target It is a state equation of the whole system that is coupled with the state equation of the seismic isolation structure itself.

(3)可変構造制御の切り換え超平面の設定
次に、(11)式で示す全体システムの状態方程式に対して、先ず双線形最適制御理論を適用して双線形最適制御則を計算する。
例えば前記非特許文献1に記載されているように、上記双線形最適制御理論は、(11)式の右辺第二項[X]uに見られるように,制御入力uと状態量[X]とが積の形で存在する、いわゆる双線形システムに分類されるシステムに対する制御理論である。
(3) Setting of switching hyperplane of variable structure control Next, a bilinear optimal control law is calculated by first applying a bilinear optimal control theory to the state equation of the whole system expressed by equation (11).
For example, as described in Non-Patent Document 1, the bilinear optimal control theory is obtained by the control input u and the state quantity [X] as seen in the second term [X] u on the right side of the equation (11). Is a control theory for systems classified as so-called bilinear systems in which and exist in the form of products.

したがって、可変減衰オイルダンパ4を制御する場合、可変減衰オイルダンパ4の減衰係数と速度が積の形で状態方程式中に現れる双線形システムになるため、本理論は効果を発揮する。
なお、通常の最適制御理論は,下式に示す線形システムに適用される。

Figure 0004747336
Accordingly, when the variable damping oil damper 4 is controlled, the present theory is effective because the damping coefficient and the speed of the variable damping oil damper 4 become a bilinear system appearing in the state equation in the form of a product.
The normal optimal control theory is applied to the linear system shown in the following equation.
Figure 0004747336

したがって、最適制御理論を適用して可変減衰ダンパの減衰係数を決定する場合は、上述したように、先ず可変減衰ダンパが発生する減衰力として制御入力uの最適制御則を計算し、次にダンパ速度で制御入力を割り算して減衰係数の指令値を評価する必要がある。
これに対して、双線形最適制御理論では、(11)式から直接制御入力、すなわち減衰係数の最適制御則が計算できる利点がある。これにより、制御性能の向上は勿論のこと、減衰係数の切換りが滑らかになるなどの効果が期待できる。
Accordingly, when the optimum control theory is applied to determine the damping coefficient of the variable damping damper, as described above, first, the optimum control law of the control input u is calculated as the damping force generated by the variable damping damper, and then the damper. It is necessary to divide the control input by the speed and evaluate the command value of the damping coefficient.
On the other hand, the bilinear optimal control theory has an advantage that the control input directly, that is, the optimal control law of the damping coefficient can be calculated from the equation (11). As a result, not only can the control performance be improved, but an effect such as smooth switching of the attenuation coefficient can be expected.

上記双線形最適制御理論では、下式に示す評価関数を導入する。

Figure 0004747336

ここで,Qは状態量に関する重み係数マトリックス、Rは[X]u、すなわち可変減衰オイルダンパ4の減衰力に対する重み係数行列であり、E(*)は数学的期待値を表すものである。 In the bilinear optimal control theory, the evaluation function shown in the following equation is introduced.
Figure 0004747336

Here, Q is a weighting coefficient matrix relating to the state quantity, R is [X] u, that is, a weighting coefficient matrix for the damping force of the variable damping oil damper 4, and E (*) represents a mathematical expectation value.

(11)式のシステムに対して、制御区間[0、T]において、

Figure 0004747336

となるような双線形最適制御則Sを求める。 In the control section [0, T] for the system of (11),
Figure 0004747336

A bilinear optimal control law S E such that

このような最適制御問題は、例えばダイナミックプログラミングの手法を用いて解くことができる。境界条件を考慮に入れHamilton-Jacobi-Bellman方程式に代入すると、求める超平面は、

Figure 0004747336

ただし、Pは以下のRicatti方程式の解である。
Figure 0004747336
Such an optimal control problem can be solved using, for example, a dynamic programming technique. Substituting into the Hamilton-Jacobi-Bellman equation taking into account the boundary conditions, the hyperplane to be obtained is
Figure 0004747336

Where P is the solution of the following Ricatti equation.
Figure 0004747336

(4)可変構造制御の設定
(11)式について、(14)式のSを得て制御を行うには、上記免震システム全体の状態量を使用して計算する必要がある。このため、免震構造物の全個所にわたってセンサを設置して、多次元にわたる膨大な量の演算を行う必要があり、コンピュータ8の負荷の問題等から現実的ではない。
このため、先ず、低次元化モデルを用いた可変構造制御を行う。
免震構造物部分のみを(5)式から複素モードを考慮して低次元化し、可変減衰オイルダンパ4の遅れを考慮したシステム全体の状態方程式を作成する。
(4) Setting of variable structure control Regarding the expression (11), in order to perform the control by obtaining the SE of the expression (14), it is necessary to calculate using the state quantity of the entire seismic isolation system. For this reason, it is necessary to install sensors over all parts of the seismic isolation structure and perform a huge amount of calculations over multiple dimensions, which is not realistic due to the load problem of the computer 8 and the like.
For this reason, first, variable structure control using a reduced-order model is performed.
Only the seismic isolation structure is reduced in order from the equation (5) in consideration of the complex mode, and the state equation of the entire system is created in consideration of the delay of the variable damping oil damper 4.

(5)式から次式の複素固有値問題を解く。

Figure 0004747336

Φは複素モード行列、Rは共役な複素固有値を持つ2行2列の対角行列であり、rは実固有値を持つm行m列の対角行列である。 The complex eigenvalue problem of the following equation is solved from the equation (5).
Figure 0004747336

Φ is a complex mode matrix, R i is a diagonal matrix of two rows and two columns with conjugate complex eigenvalues, r m is the diagonal matrix of m rows and m columns with a real eigenvalues.

ここで、複素モード行列の要素は,以下の通りである。

Figure 0004747336

ここで、σとτは実数ベクトル、i=(―1)1/2である。 Here, the elements of the complex mode matrix are as follows.
Figure 0004747336

Here, σ i and τ i are real vectors, i = (− 1) 1/2 .

そして、次式の変換行列を定義し、(5)式に対して座標変換を行うと、複素モードを実数化することができる。

Figure 0004747336

ここで、A=T−1AT、B=T−1B、E=T−1Eである。 If a transformation matrix of the following equation is defined and coordinate transformation is performed on the equation (5), the complex mode can be converted to a real number.
Figure 0004747336

Here, A q = T −1 AT, B q = T −1 B, and E q = T −1 E.

次に、実数モードを除外して実数化された複素モードをr次元まで抽出して低次元化すると、

Figure 0004747336
Next, when the real mode is excluded and the realized complex mode is extracted down to r-dimension and reduced in dimension,
Figure 0004747336

上式から、可変減衰オイルダンパ4の遅れを考慮した低次元化したシステム全体の状態方程式を設定すると、

Figure 0004747336
From the above equation, when setting the state equation of the entire reduced system considering the delay of the variable damping oil damper 4,
Figure 0004747336

次いで、(14)式で評価した全体モデルに対する双線形最適制御則に対して(20)式の変換行列から制御対象モードに対する変換行列Tを抽出して作成した下記の変換行列Hを用いて、(24)式で表される低次元化モデルに対する双線形最適制御則を計算する。

Figure 0004747336

はTから制御対象モードのモードベクトルを抽出した変換行列であり、IとIは単位行列である。 Then, using the transformation matrix H below created by extracting the transformation matrix T R for the control target mode from the bilinear optimal control law with respect to (20) of the transformation matrix for the entire model was evaluated in (14) The bilinear optimal control law for the reduced-order model expressed by the equation (24) is calculated.
Figure 0004747336

T R is a transformation matrix obtained by extracting mode vector of the controlled object mode from T, I F and I d is the identity matrix.

なお、低次元化モデルの状態量qは、カルマンフィルターなどのオブザーバを用いて推定することができる。例えば、外乱包含拡大系の全体状態方程式である(11)式の観測状態量をYとすると、出力方程式は以下の通り記述することが出来る。
=C (27)
ここで、(27)式の出力方程式に、(26)式で用いた変換行列Hを適用すると、低次元化した状態方程式に対する出力方程式が求められる。
=C (28)
ただし、C=CHである。
The state quantity q R of reduced-order model can be estimated using the observer, such as a Kalman filter. For example, if an overall state equation of the disturbance include expansion system (11) an observation state of formulas and Y E, the output equation can be described as follows.
Y E = C E X E (27)
Here, when the transformation matrix H used in the equation (26) is applied to the output equation of the equation (27), an output equation for the reduced state equation is obtained.
Y R = C R q R (28)
However, C R = C E H.

したがって、(23)式の状態方程式および(28)式の出力方程式の組み合わせに対して、免震システムの検知個所やその数を設定してオブザーバーを構成すればよいことになる。
このような低次元化モデルでの可変構造制御は、(26)式を切換え超平面として、等価制御入力を設定する。

Figure 0004747336
Therefore, an observer may be configured by setting the detection locations and the number of seismic isolation systems for the combination of the equation of state (23) and the output equation (28).
In the variable structure control in such a reduced-dimensional model, the equivalent control input is set using the equation (26) as a switching hyperplane.
Figure 0004747336

非線形制御入力は、切換面からの偏差を抑制する働きがあるが、制御対象の状態が切換面近傍に近付いた時に制御が高速で切換えられる、いわゆるチャタリングを起こす可能性がある。そのため、(30)式では、δとδという定数を導入して、チャタリングを除去するために制御入力の平滑化を行っている。 The non-linear control input has a function of suppressing deviation from the switching surface, but may cause so-called chattering in which the control is switched at high speed when the state of the controlled object approaches the vicinity of the switching surface. Therefore, in equation (30), constants δ 1 and δ 2 are introduced to smooth the control input in order to eliminate chattering.

(29)式と(30)式とにより、建物2の加速度を低減化させるべくシステムに与えるべき制御入力uは、

Figure 0004747336

により計算される。 According to the equations (29) and (30), the control input u to be given to the system to reduce the acceleration of the building 2 is
Figure 0004747336

Is calculated by

この結果、CHとCLとの2値切換え式の可変減衰オイルダンパ4を用いた場合の可変構造制御の切換え制御則は、(31)式から以下のように決定される。

Figure 0004747336

コンピュータ8は、可変減衰オイルダンパ4におけるx方向とy方向の減衰係数切換え基準値が、それぞれcSX、cSYであるとすると、X方向についての減衰係数cがこれ以下となった場合には、制御すべき減衰係数cとしてCLを設定し、切り換え指令信号を発する。また、逆に上記減衰係数cがcSXを超えた場合には、制御すべき減衰係数cとしてCHを設定し、切り換え指令信号を発する。 As a result, the switching control law of the variable structure control in the case of using the variable damping oil damper 4 of the binary switching type between CH and CL is determined as follows from the formula (31).
Figure 0004747336

If the damping coefficient switching reference values in the x direction and y direction in the variable damping oil damper 4 are c SX and c SY , respectively, the computer 8 determines that the damping coefficient c in the X direction is less than this. Then, CL is set as the attenuation coefficient c X to be controlled, and a switching command signal is issued. Conversely, when the attenuation coefficient c exceeds c SX , CH is set as the attenuation coefficient c X to be controlled, and a switching command signal is issued.

また、Y方向についても同様に、減衰係数cがcSY以下となった場合には、制御すべき減衰係数cとしてCLを設定し、切り換え指令信号を発する。他方、上記減衰係数cがcSYを超えた場合には、制御すべき減衰係数cとしてCHを設定して切り換え指令信号を発する。 Similarly, the Y direction, when the damping coefficient c is equal to or less than c SY sets CL as the attenuation coefficient c Y to be controlled, issues a switching instruction signal. On the other hand, when the damping coefficient c exceeds c SY, set the CH issues a switching command signal as a damping coefficient c Y to be controlled.

さらに、この免震システムにおいては、予めコンピュータ8に、当該免震構造物における低応答値Aminおよび高応答値Amax並びに低応答値Aminよりも低い制御開始を判断する応答値Ainiが設定されており、センサ5a、5b、5cからの検出信号が制御開始を判断する応答値Aini未満である場合に減衰力c、cを最大値CHに設定し、制御開始を判断する応答値Ainiを超えて低応答値Amin未満である場合に減衰力c、cを最小値CLに設定するようになっている。 Further, in this seismic isolation system, advance in the computer 8, the response value A ini for determining lower control start than the low response value A min and a high response value A max and low response value A min in the seismic isolation structure When it is set and the detection signals from the sensors 5a, 5b, and 5c are less than the response value A ini for determining the start of control, the damping forces c X and c Y are set to the maximum value CH to determine the start of control. When the response value A ini is exceeded and less than the low response value A min , the damping forces c X and c Y are set to the minimum value CL.

そして、上記検出信号が、低応答値Aminを初めて超えた以降は、高応答値Amax以下の範囲である場合に、上記演算を実行し、高応答値Amaxを超えた場合に減衰力c、cを最大値CHに設定するようになっている。
加えて、平常時のノイズに対して制御手段12が応答しないように、コンピュータ8は、基礎1のセンサ5aと建物2下部のセンサ5bとによって得られる上下部構造間の相対速度または可変減衰オイルダンパ4の応答速度が規定値以下である場合に、上記減衰力の切換えを行わないように設定されている。
Then, after the detection signal exceeds the low response value Amin for the first time, the above calculation is performed when the detection signal is in the range of the high response value Amax or less, and when the detection signal exceeds the high response value Amax , the damping force c X and c Y are set to the maximum value CH.
In addition, the relative speed between the upper and lower structures obtained by the sensor 5a of the foundation 1 and the sensor 5b of the lower part of the building 2 or a variable damping oil is obtained so that the control means 12 does not respond to noise during normal times. When the response speed of the damper 4 is equal to or less than a specified value, the damping force is set not to be switched.

このようにして、コンピュータ8から可変減衰オイルダンパ4に対する減衰係数切換の指令信号が発せられ、電源回路9から供給される交流の駆動電流がオンになると電磁弁7が開弁し、駆動電流がオフになると電磁弁7が閉弁することになる。
この際に、この免震システムにおいては、交流電流センサ11によって電磁弁7に供給される交流電流を測定し、検出された電流値が異常値か否かにより電磁弁7の作動状態が正常か異常かがコンピュータ8に格納された制御プログラム(判定手段)によって判定されるようになっている。
In this way, when the computer 8 issues a command signal for switching the damping coefficient to the variable damping oil damper 4, and when the AC driving current supplied from the power supply circuit 9 is turned on, the electromagnetic valve 7 is opened and the driving current is reduced. When turned off, the solenoid valve 7 is closed.
At this time, in this seismic isolation system, the alternating current supplied to the electromagnetic valve 7 by the alternating current sensor 11 is measured, and whether the operating state of the electromagnetic valve 7 is normal or not is determined by whether or not the detected current value is an abnormal value. Abnormality is determined by a control program (determination means) stored in the computer 8.

図13は、駆動電流の波形の一例を示す波形図である。図13中、グラフIaは電磁弁7に供給される交流電流の波形を示し、グラフIbは交流電流を平均化した駆動電流の実効値を示すものである。
図13に示すように、交流電流センサ11により測定された電磁弁7の駆動電流Iaは、制御回路12において、波形整形及び平均化される。制御回路12は、平均化した波形に基づいて電磁弁の異常の有無を判定する。本実施形態における平均化処理では、任意の周期の1サイクルに亘る二乗平均平方根の実効値(RMS:Root mean square value)Ibを求める。
FIG. 13 is a waveform diagram showing an example of the waveform of the drive current. In FIG. 13, the graph Ia shows the waveform of the alternating current supplied to the electromagnetic valve 7, and the graph Ib shows the effective value of the drive current obtained by averaging the alternating current.
As shown in FIG. 13, the drive current Ia of the electromagnetic valve 7 measured by the alternating current sensor 11 is waveform-shaped and averaged in the control circuit 12. The control circuit 12 determines whether there is an abnormality in the electromagnetic valve based on the averaged waveform. In the averaging process in the present embodiment, an effective value (RMS: Root mean square value) Ib over one cycle of an arbitrary period is obtained.

図14は制御回路12が電磁弁7を開弁させて減衰係数を切り替えたときの駆動電流の波形を示す図である。
図14に示すように、制御回路12が減衰係数指令をCH(C1HまたはC2H)からCL(C1LまたはC2L)に切り替えたときは、瞬間的に大きな起動電流aが流れ、電磁弁7が開弁動作し、その後安定し、規定保持電流bが流れ続ける。例えば、電磁弁7で異物の挟み込みなどにより動作が途中で止まった場合は、起動電流と保持電流との間の電流値が流れ続ける。そして、制御回路12では、減数係数指令をCH(C1HまたはC2H)からCL(C1LまたはC2L)に切り替えたときから保持電流bになるまでの時間をdTとし、微小電流をIとし、dTとIをしきい値として異常判定処理を行う。
FIG. 14 is a diagram showing the waveform of the drive current when the control circuit 12 opens the electromagnetic valve 7 and switches the attenuation coefficient.
As shown in FIG. 14, when the control circuit 12 switches the damping coefficient command from CH (C1H or C2H) to CL (C1L or C2L), a large starting current a instantaneously flows and the solenoid valve 7 opens. It operates and then stabilizes, and the specified holding current b continues to flow. For example, when the operation stops midway due to foreign matter being caught by the electromagnetic valve 7, the current value between the starting current and the holding current continues to flow. Then, the control circuit 12, the time of the subtrahend coefficients command from CH (C1H or C2H) until the holding current b from the time of switching to CL (C1L or C2L) and dT, the minute current and I 0, and dT Abnormality determination processing is performed using I0 as a threshold value.

図15は制御回路12が電磁弁7を閉弁させて減衰係数を切り替えたときの駆動電流の波形を示す図である。
図15に示すように、制御回路12が減衰係数指令をCL(C1LまたはC2L)からCH(C1HまたはC2H)に切り替えたときは、電流が0Aに低下することになる。そして、制御回路12では、減衰係数指令をCL(C1LまたはC2L)からCH(C1HまたはC2H)に切り替えたときから、保持電流bから微小電流をIになるまでの時間をdTとし、微小電流をIとし、dTとIをしきい値として異常判定処理を行う。
FIG. 15 is a diagram showing the waveform of the drive current when the control circuit 12 closes the electromagnetic valve 7 and switches the attenuation coefficient.
As shown in FIG. 15, when the control circuit 12 switches the attenuation coefficient command from CL (C1L or C2L) to CH (C1H or C2H), the current decreases to 0A. Then, in the control circuit 12, the time from when the attenuation coefficient command is switched from CL (C1L or C2L) to CH (C1H or C2H) until the minute current is changed from the holding current b to I 0 is defined as dT. It was a I 0, the abnormality determination process of dT and I 0 as a threshold.

なお、上記dTとIのしきい値を小さくした場合、誤検出となる可能性が高くなり、dTとIのしきい値を大きくした場合、電磁弁7の異常を検出できないおそれが生じる。そのため、本実施形態では、dTを切り替えたときから規定保持電流になるまでの時間の2倍程度に設定し、Iを規定保持電流の1/5程度に設定する。
ここで、制御回路12が実行する制御処理について図16〜図18を参照して説明する。
If the threshold values of dT and I 0 are reduced, the possibility of erroneous detection increases. If the threshold values of dT and I 0 are increased, there is a possibility that an abnormality of the electromagnetic valve 7 cannot be detected. . Therefore, in this embodiment, set to about twice the time until the prescribed holding current from the time of switching the dT, set to about 1/5 of the prescribed holding current I 0.
Here, the control processing executed by the control circuit 12 will be described with reference to FIGS.

図16に示されるように、制御回路12は、S11において、制御中か否かを確認する。S11で制御中でないときは、S12に進み、定期点検時刻かどうかをチェックする。
S12で定期点検時刻になると、S13に進み、指定された電磁弁7に開指令を出力する。これにより、当該電磁弁7は開弁動作を行う。次のS14では、電磁弁開異常判定処理を実行する。
As shown in FIG. 16, the control circuit 12 confirms whether or not the control is being performed in S11. When the control is not in S11, the process proceeds to S12 to check whether it is a periodic inspection time.
When the periodic inspection time comes in S12, the process proceeds to S13, and an open command is output to the designated solenoid valve 7. Thereby, the said electromagnetic valve 7 performs valve opening operation | movement. In next S14, an electromagnetic valve opening abnormality determination process is executed.

S15では、電磁弁開異常判定処理の判定結果が正常かどうかをチェックする。S15において、電磁弁開異常判定処理の判定結果が正常である場合は、S16に進み、電磁弁7に閉指令を出力する。これにより、当該電磁弁7は閉弁動作を行う。次のS17では、電磁弁閉異常判定処理を実行する。
S18では、電磁弁閉異常判定処理の判定結果が正常かどうかをチェックする。S18において、電磁弁開異常判定処理の判定結果が正常である場合は、S19に進み、次に異常判定を行う電磁弁7を指定する。そして、上記S11に戻る。
In S15, it is checked whether the determination result of the electromagnetic valve open abnormality determination process is normal. If the determination result of the electromagnetic valve open abnormality determination process is normal in S15, the process proceeds to S16, and a close command is output to the electromagnetic valve 7. Thereby, the said electromagnetic valve 7 performs valve closing operation | movement. In the next S17, an electromagnetic valve closing abnormality determination process is executed.
In S18, it is checked whether the determination result of the electromagnetic valve closing abnormality determination process is normal. In S18, when the determination result of the electromagnetic valve opening abnormality determination process is normal, the process proceeds to S19, and the electromagnetic valve 7 that performs abnormality determination next is designated. Then, the process returns to S11.

また、上記S15及びS18で判定結果が異常である場合には、S20に進み、当該電磁弁7の異常内容を記憶する。これにより、多数配置された電磁弁7のうちどの電磁弁7でどのような異常(開弁不良あるいは閉弁不良か)が容易に分かる。
上記S11において、制御中であるときは、S21に進み、減衰係数指令値がCHからCLに切り替わったかどうかをチェックする。S21において、減衰係数指令値がCHからCLに切り替わった場合には、S22に進み、上記S14と同様に電磁弁開異常判定処理を実行する。
If the determination result is abnormal in S15 and S18, the process proceeds to S20, and the abnormality content of the electromagnetic valve 7 is stored. Thereby, it can be easily understood which abnormality (open valve failure or closed valve failure) is in any one of the many solenoid valves 7 arranged.
When the control is being performed in S11, the process proceeds to S21 to check whether or not the attenuation coefficient command value has been switched from CH to CL. In S21, when the damping coefficient command value is switched from CH to CL, the process proceeds to S22, and the electromagnetic valve opening abnormality determination process is executed as in S14.

また、上記S21において、減衰係数指令値がCHからCLに切り替わっていない場合には、S23に進み、減衰係数指令値がCLからCHに切り替わったかどうかをチェックする。S23において、減衰係数指令値がCLからCHに切り替わった場合には、S24に進み、上記S17と同様に電磁弁閉異常判定処理を実行する。
そして、S25に進み、S22またはS24において、電磁弁閉異常判定処理の判定結果が正常かどうかをチェックする。S25において、電磁弁開異常判定処理の判定結果が正常である場合は、S26に進み、次に異常判定を行う電磁弁7を指定する。そして、上記S11に戻る。
If the attenuation coefficient command value has not been switched from CH to CL in S21, the process proceeds to S23 to check whether the attenuation coefficient command value has been switched from CL to CH. If the damping coefficient command value is switched from CL to CH in S23, the process proceeds to S24, and the electromagnetic valve closing abnormality determination process is executed in the same manner as in S17.
Then, the process proceeds to S25, and in S22 or S24, it is checked whether or not the determination result of the electromagnetic valve closing abnormality determination process is normal. In S25, when the determination result of the electromagnetic valve opening abnormality determination process is normal, the process proceeds to S26, and the electromagnetic valve 7 that performs abnormality determination next is designated. Then, the process returns to S11.

しかし、S25において、電磁弁開異常判定処理の判定結果が異常である場合は、S20に進み、当該電磁弁7の異常内容を記憶する。これにより、多数配置された電磁弁7のうちどの電磁弁7でどのような異常(開弁不良あるいは閉弁不良か)が容易に分かる。よって、制御回路12は、どの電磁弁7で異常があるか否かを監視することが可能になる。そのため、メンテナンス作業の手間を大幅に緩和することが可能になると共に、電磁弁7の異常発生を判定した場合には、直ちに異常のある電磁弁7を修理することが可能になる。   However, if the determination result of the electromagnetic valve opening abnormality determination process is abnormal in S25, the process proceeds to S20, and the abnormality content of the electromagnetic valve 7 is stored. Thereby, it can be easily understood which abnormality (open valve failure or closed valve failure) is in any one of the many solenoid valves 7 arranged. Therefore, the control circuit 12 can monitor which electromagnetic valve 7 is abnormal. Therefore, it is possible to greatly reduce the labor of maintenance work, and when it is determined that an abnormality has occurred in the electromagnetic valve 7, it is possible to immediately repair the abnormal electromagnetic valve 7.

ここで、上記S14及びS22で実行された電磁弁開異常判定処理について説明する。
図17に示されるように、制御回路50は、電磁弁開異常判定処理が開始されると、S31において、交流電流センサ11により電磁弁7に供給される駆動電流を計測する。次のS32では、dT時間(しきい値)以内に電磁弁7の駆動電流がIアンペア以上となるかどうかをチェックする。S32において、dT時間以内に電磁弁7の駆動電流がIアンペア以上となったときは、S33に進み、判定結果を正常とする。また、S32において、dT時間以内に電磁弁7の駆動電流がIアンペア以上とならなかったときは、S34に進み、判定結果を異常とする。
Here, the electromagnetic valve opening abnormality determination process executed in S14 and S22 will be described.
As shown in FIG. 17, when the electromagnetic valve open abnormality determination process is started, the control circuit 50 measures the drive current supplied to the electromagnetic valve 7 by the AC current sensor 11 in S31. In next S32, it is checked whether or not the drive current of the solenoid valve 7 is equal to or more than 10 amperes within the dT time (threshold). In S32, when the drive current of the solenoid valve 7 becomes equal to or greater than 10 amperes within dT time, the process proceeds to S33 and the determination result is normal. In S32, when the drive current of the solenoid valve 7 does not become 10 amperes or more within dT time, the process proceeds to S34, and the determination result is abnormal.

このように、電磁弁開異常判定処理では、図14に示すように、減数係数指令をCH(C1HまたはC2H)からCL(C1LまたはC2L)に切り替えたときから保持電流bになるまでの時間をdTとし、微小電流をIとし、dTとIをしきい値として異常判定処理を行うため、電磁弁7の開弁動作不良の有無を正確に判定することができる。 Thus, in the electromagnetic valve opening abnormality determination process, as shown in FIG. 14, the time from when the reduction coefficient command is switched from CH (C1H or C2H) to CL (C1L or C2L) until the holding current b is reached. Since abnormality determination processing is performed using dT, a minute current as I 0, and dT and I 0 as threshold values, it is possible to accurately determine the presence or absence of the valve opening malfunction of the electromagnetic valve 7.

ここで、上記S17及びS24で実行された電磁弁閉異常判定処理について説明する。
図18に示されるように、制御回路12は、電磁弁閉異常判定処理が開始されると、S41において、交流電流センサ11により電磁弁7に供給される駆動電流を計測する。次のS42では、dT時間(しきい値)以内に電磁弁7の駆動電流がIアンペア未満となるかどうかをチェックする。S42において、dT時間以内に電磁弁7の駆動電流がIアンペア未満となったときは、S43に進み、判定結果を正常とする。また、S42において、dT時間以内に電流がIアンペア未満とならなかったときは、S44に進み、判定結果を異常とする。
Here, the electromagnetic valve closing abnormality determination process executed in S17 and S24 will be described.
As shown in FIG. 18, when the electromagnetic valve closing abnormality determination process is started, the control circuit 12 measures the drive current supplied to the electromagnetic valve 7 by the AC current sensor 11 in S41. In the next S42, it is checked whether or not the drive current of the solenoid valve 7 becomes less than 10 amperes within the dT time (threshold value). In S42, when the drive current of the solenoid valve 7 becomes less than 10 amperes within dT time, the process proceeds to S43 and the determination result is normal. If the current does not become less than 10 amperes within dT time in S42, the process proceeds to S44, and the determination result is abnormal.

このように、電磁弁閉異常判定処理では、図15に示すように、減数係数指令をCL(C1LまたはC2L)からCH(C1HまたはC2H)に切り替えたときから、保持電流bから微小電流をIになるまでの時間をdTとし、微小電流をIとし、dTとIをしきい値として異常判定処理を行うため、電磁弁7の閉弁動作不良の有無を正確に判定することができる。よって、制御回路12は、多数の電磁弁7を作動させてどの電磁弁7で異常があるか否かを監視することが可能になる。そのため、メンテナンス作業の手間を大幅に緩和することが可能になると共に、電磁弁7の異常発生を判定した場合には、直ちに制御を停止させてシステムの安全性を高めることができ、また異常のある電磁弁7を修理することが可能になる。 In this way, in the electromagnetic valve closing abnormality determination process, as shown in FIG. 15, the minute current is reduced from the holding current b to I after the reduction coefficient command is switched from CL (C1L or C2L) to CH (C1H or C2H). Since the time to reach 0 is dT, the minute current is I 0 and dT and I 0 are threshold values, abnormality determination processing is performed. it can. Therefore, the control circuit 12 can monitor which electromagnetic valve 7 is abnormal by operating a large number of electromagnetic valves 7. Therefore, it is possible to greatly reduce the labor of maintenance work, and when it is determined that an abnormality has occurred in the solenoid valve 7, it is possible to immediately stop the control and increase the safety of the system. A certain solenoid valve 7 can be repaired.

(実施例)
本発明の可変構造セミアクティブ免震システムの効果を実証するために、図19および図20に示すような、地上部が13階建ての高層棟と、5階建ての低層棟からなる免震構造物に、本発明に係る可変構造セミアクティブ免震システムを適用した場合についての解析を行った。
この可変構造セミアクティブ免震システムにおいては、図20に示すように、x方向とy方向の両方向に対して、可変減衰オイルダンパを5台、減衰係数固定型オイルダンパ(パッシブオイルダンパ)を5台設置することとした。
(Example)
In order to demonstrate the effect of the variable structure semi-active seismic isolation system of the present invention, as shown in FIGS. 19 and 20, the ground-isolated structure is composed of a 13-story high-rise building and a 5-story low-rise building. We analyzed the case where the variable structure semi-active seismic isolation system according to the present invention was applied to the object.
In this variable structure semi-active seismic isolation system, as shown in FIG. 20, five variable damping oil dampers and five fixed damping coefficient oil dampers (passive oil dampers) are provided in both the x and y directions. It was decided to install a stand.

また、高層棟頂部のx、y方向加速度、低層棟頂部のx、y方向加速度、免震層下部のx、y方向入力加速度、免震層のx,y方向相対変位を観測状態量として使用し、カルマンフィルターにより前記(24)式の低次元化モデルの状態量を推定した。なお、本実施例においては、免震構造物の低次5モードの応答を制御対象とした。   In addition, x and y direction acceleration at the top of the high-rise building, x and y direction acceleration at the top of the low-rise building, x and y direction input acceleration at the bottom of the base isolation layer, and x and y direction relative displacement of the base isolation layer are used as observation state quantities. Then, the state quantity of the reduced-order model of the equation (24) was estimated by the Kalman filter. In this embodiment, the response of the low-order five modes of the seismic isolation structure is set as a control target.

図22は、上記パッシブダンパと可変減衰オイルダンパのそれぞれの減衰力特性を示すものである。図22に示すように、可変減衰オイルダンパとしては、減衰係数が2段階に切換え可能な可変減衰ダンパを使用した。
また、図23に示すように、本実施例においても、上記実施形態において図2示したシステム構成と同様の構成を使用した。
FIG. 22 shows the damping force characteristics of the passive damper and the variable damping oil damper. As shown in FIG. 22, as the variable damping oil damper, a variable damping damper whose damping coefficient can be switched in two stages was used.
Further, as shown in FIG. 23, the same configuration as the system configuration shown in FIG.

図24は、免震構造物の解析モデルを示すものである。低層棟各質点については層間のせん断ばねを考慮してx、y、θの3自由度を、免震層についてもx、y、θの3自由度を、高層棟各質点については層間のせん断ばねと曲げまねの効果を考慮して、x、y、θ、θY、θの5自由度を設定して解析モデルを作成した。解析モデルの固有値解析結果を表1に示す。 FIG. 24 shows an analysis model of a seismic isolation structure. For low-rise buildings the mass points in consideration of the shear spring interlayer x, y, three degrees of freedom of the theta Z, x also the base isolation layer, y, three degrees of freedom of the theta Z, layers for high-rise buildings the mass points In consideration of the effects of shear spring and bending imitation, an analysis model was created by setting five degrees of freedom of x, y, θ X , θ Y, and θ Z. Table 1 shows the eigenvalue analysis results of the analysis model.

Figure 0004747336
Figure 0004747336

次いで、静的縮合を行い、最終的に、低層棟、免震層、高層棟の全ての質点についてx、y、θの3自由度を持つ制御設計用解析モデルを作成した。
次いで、0.1Hzから10Hzまでフラットな特性を持つホワイトノイズを地震動の入力加速度として用いて、可変減衰オイルダンパを制御する場合、全てのオイルダンパをパッシブダンパにした場合、可変減衰オイルダンパの減衰係数をCL(cmin)に固定した場合、およびCH(cmax)に固定した場合について、制御性能を比較した。ここで、可変減衰オイルダンパの減衰係数の減衰係数切換え基準値は、
SX=cSY=25.0kN・s/cmとした。
Then carried out statically condensation, ultimately, low-rise buildings, created isolation layer, x for all mass points of high-rise buildings, y, the analysis model for control design with three degrees of freedom of the theta Z.
Next, when controlling variable damping oil dampers using white noise with flat characteristics from 0.1 Hz to 10 Hz as the input acceleration of seismic motion, if all oil dampers are passive dampers, damping of variable damping oil dampers The control performance was compared when the coefficient was fixed to CL (c min ) and when the coefficient was fixed to CH (c max ). Here, the damping coefficient switching reference value of the damping coefficient of the variable damping oil damper is
c SX = c SY = 25.0 kN · s / cm.

図25の(a)、(c)は、それぞれ免震層におけるx方向およびy方向の変形と地震動の入力加速度に対する伝達関数を示すものであり、図25(b)、(d)は、それぞれ高層棟の頂部におけるx方向およびy方向の加速度と地震動の入力加速度に対する伝達関数を示すものである。   25 (a) and 25 (c) show the transfer functions for the x- and y-direction deformations in the seismic isolation layer and the input acceleration of the ground motion. FIGS. 25 (b) and 25 (d) The transfer function with respect to the acceleration of x direction and y direction in the top part of a high-rise building and the input acceleration of a ground motion is shown.

図25に示す伝達関数の比較から、免震層の変形については、制御時の応答はパッシブの応答とほぼ同等であることがわかる。他方、高層棟頂部の加速度伝達率については、1次モードに対応する0.2Hz付近では、制御時はパッシブとほぼ同等であるが、2次モードのx方向における2Hz付近、y方向における1Hzを越えた付近の加速度伝達率は、それぞれ制御時の結果が一番小さくなっている。
この結果、制御によって免震層の変形をパッシブとほぼ同等に保ちつつ、建物頂部の応答加速度を低減する効果が得られることがわかる。
From the comparison of the transfer functions shown in FIG. 25, it is understood that the response at the time of control is almost equivalent to the passive response for the deformation of the seismic isolation layer. On the other hand, the acceleration transmissibility at the top of the high-rise building is approximately the same as passive at the time of control around 0.2 Hz corresponding to the primary mode, but around 2 Hz in the x direction of the secondary mode and 1 Hz in the y direction. As for the acceleration transmissibility in the vicinity, the result at the time of control is the smallest.
As a result, it can be seen that the effect of reducing the response acceleration at the top of the building can be obtained while maintaining the deformation of the seismic isolation layer substantially the same as the passive by the control.

次に、地震動として、1968年八戸波NS方向加速度記録を、最大速度25kineに基準化して作用させた場合の制御時の応答解析結果と、全てのオイルダンパをパッシブダンパにした場合の制御性能を時刻歴応答解析により比較した。   Next, as seismic motion, the 1968 Hachinohe NS direction acceleration record was standardized to a maximum speed of 25 kine, and the response analysis result at the time of control and the control performance when all oil dampers were made passive dampers Comparison was made by time history response analysis.

図26は、上記解析結果を示すものであり、(a)は免震層の応答変形を、(b)は高層棟頂部の応答加速度を、(c)は可変減衰ダンパの減衰力をパッシブと比較して示したものである。なお、同図(d)は、この際の可変減衰オイルダンパの減衰係数の切換え状況を示すものである。   FIG. 26 shows the results of the above analysis, where (a) shows the response deformation of the seismic isolation layer, (b) shows the response acceleration at the top of the high-rise building, and (c) shows the damping force of the variable damping damper as passive. It is shown in comparison. FIG. 6D shows the switching state of the damping coefficient of the variable damping oil damper at this time.

図26に示す時刻歴応答解析の結果から、制御時の免震層の応答変形は、パッシブとほぼ同等であるのに対して、高層棟頂部の応答加速度は、パッシブに比べて30%程度低減されていることから、地震時における免震構造物全体としての加速度を効果的に低減化し得ることが判る。   From the results of the time history response analysis shown in FIG. 26, the response deformation of the seismic isolation layer at the time of control is almost the same as that of the passive, whereas the response acceleration at the top of the high-rise building is reduced by about 30% compared to the passive. Therefore, it can be seen that the acceleration of the seismic isolation structure as a whole during an earthquake can be effectively reduced.

本発明に係る可変構造セミアクティブ免震システムの一実施形態を示す全体の概略構成図である。1 is an overall schematic configuration diagram showing an embodiment of a variable structure semi-active seismic isolation system according to the present invention. 上記免震システムのブロック図である。It is a block diagram of the said seismic isolation system. 図2の電源回路部分を詳細に示すブロック図である。It is a block diagram which shows the power supply circuit part of FIG. 2 in detail. 図1の可変減衰オイルダンパの構成を示す縦断面図である。It is a longitudinal cross-sectional view which shows the structure of the variable damping oil damper of FIG. 図4のV―V線視断面図である。FIG. 5 is a cross-sectional view taken along line VV in FIG. 4. 図4のVI―VI線視断面図である。FIG. 6 is a sectional view taken along line VI-VI in FIG. 4. 上記可変減衰オイルダンパのピストン速度に対する減衰係数の変化を示すグラフである。It is a graph which shows the change of the damping coefficient with respect to the piston speed of the said variable damping oil damper. 上記減衰係数をCHとする場合の各弁の開閉状態を示す系統図である。It is a systematic diagram which shows the opening-and-closing state of each valve when the said attenuation coefficient is set to CH. 上記減衰係数をCLとする場合の各弁の開閉状態を示す系統図である。It is a systematic diagram which shows the open / closed state of each valve when the said attenuation coefficient is set to CL. 可変減衰オイルダンパのMaxwellモデルを示す図である。It is a figure which shows the Maxwell model of a variable damping oil damper. 加減減衰オイルダンパの減衰係数の切換え時間遅れの実測値を示すグラフである。It is a graph which shows the actual value of the switching time delay of the damping coefficient of an adjustable damping oil damper. 成形フィルタのバンドパス特性を示すグラフである。It is a graph which shows the band pass characteristic of a shaping | molding filter. 上記可変減衰オイルダンパの電磁弁への駆動電流の変化を示す波形図である。It is a wave form diagram which shows the change of the drive current to the solenoid valve of the said variable damping oil damper. 減衰係数をCHからCLに切り換えた際の駆動電流の変化を示す波形図である。It is a wave form diagram which shows the change of the drive current at the time of switching an attenuation coefficient from CH to CL. 減衰係数をCLからCHに切り換えた際の駆動電流の変化を示す波形図である。It is a wave form diagram which shows the change of the drive current at the time of switching an attenuation coefficient from CL to CH. 制御回路が実行する制御処理を示すフローチャートである。It is a flowchart which shows the control processing which a control circuit performs. 電磁弁開の異常判定処理を示すフローチャートである。It is a flowchart which shows the abnormality determination process of a solenoid valve opening. 電磁弁閉の異常判定処理を示すフローチャートである。It is a flowchart which shows the abnormality determination process of a solenoid valve closing. 本実施例において解析対象とした免震構造物を示す縦断面図である。It is a longitudinal cross-sectional view which shows the seismic isolation structure made into the analysis object in a present Example. 図19のダンパ配置を示す平面図である。FIG. 20 is a plan view showing the damper arrangement of FIG. 19. 本実施例に使用した可変減衰オイルダンパのピストン速度と減衰力との関係を示すグラフである。It is a graph which shows the relationship between the piston speed and damping force of the variable damping oil damper used for the present Example. 本実施例に使用したオイルダンパの特性を示す図表である。It is a graph which shows the characteristic of the oil damper used for the present Example. 本実施例のシステム構成図である。It is a system configuration diagram of the present embodiment. 本実施例の免震構造物の解析モデルを示す図である。It is a figure which shows the analysis model of the seismic isolation structure of a present Example. 本実施例におけるホワイトノイズ加振時の伝達関数の解析結果を示すグラフである。It is a graph which shows the analysis result of the transfer function at the time of white noise excitation in a present Example. 本実施例における時刻歴応答解析の結果を示すグラフである。It is a graph which shows the result of the time history response analysis in a present Example.

符号の説明Explanation of symbols

1 基礎(下部構造)
2 建物(上部構造)
3 免震装置
4 可変減衰オイルダンパ(可変減衰ダンパ)
5a、5b、5c、5d センサ
6 制御手段
7 電磁弁(切換手段)
8 コンピュータ
9 電源回路
11 交流電流センサ
12 制御回路
1 Foundation (substructure)
2 Building (superstructure)
3 Seismic isolation device 4 Variable damping oil damper (variable damping damper)
5a, 5b, 5c, 5d Sensor 6 Control means 7 Solenoid valve (switching means)
8 Computer 9 Power supply circuit 11 AC current sensor 12 Control circuit

Claims (6)

上下部構造物間に免震装置が介装された免震構造物の上記上部構造と上記下部構造との間に設けられた可変減衰ダンパと、地震時における上記上部構造および下部構造の応答を検出するセンサと、このセンサからの検出信号に基づいて上記可変減衰ダンパへの指令信号を演算・出力するコンピュータを備えた制御手段と、この制御手段から出力された上記指令信号により上記可変減衰ダンパによる減衰力を切り換える切換手段とを備えた可変構造セミアクティブ免震システムにおいて、
上記コンピュータは、上記免震構造物自体の状態方程式に、上記可変減衰ダンパの減衰係数設定値とダンパ速度との積で規定される減衰力設定値と上記可変減衰ダンパで実際に発生する減衰力との一次遅れ系の関係式であって、かつ時定数を一定と近似した関係式を連成させた、全体システムの状態方程式に対して双線形最適制御理論を適用して求めた制御則を切換え超平面として設定し、上記切換え超平面に基づいて上記可変減衰ダンパの減衰係数を切り換えて上記上部構造物の応答変形および/または応答加速度を制御する可変構造制御法を採用して上記指令信号を出力することを特徴とする可変構造セミアクティブ免震システム。
The variable damping damper provided between the upper structure and the lower structure of the base isolation structure in which the base isolation device is interposed between the upper and lower structures, and the response of the upper structure and the lower structure at the time of the earthquake. A sensor for detecting, a control means having a computer for calculating and outputting a command signal to the variable damping damper based on a detection signal from the sensor, and the variable damping damper by the command signal output from the control means. In the variable structure semi-active seismic isolation system with switching means for switching the damping force by
The computer includes a damping force setting value defined by a product of a damping coefficient setting value and a damper speed of the variable damping damper and a damping force actually generated by the variable damping damper in the state equation of the base isolation structure itself. The control law obtained by applying the bilinear optimal control theory to the state equation of the entire system , which is a relational expression of the first-order lag system and the relational expression approximating that the time constant is constant. The command signal is set using a variable structure control method which is set as a switching hyperplane and switches the damping coefficient of the variable damping damper based on the switching hyperplane to control response deformation and / or response acceleration of the upper structure. A variable structure semi-active seismic isolation system characterized by
上記コンピュータは、予め上記免震構造物の低応答値および高応答値並びに上記低応答値よりも低い制御開始を判断する応答値が設定されており、上記センサからの上記検出信号が上記制御開始を判断する応答値未満である場合に上記減衰力を最大値に設定し、上記制御開始を判断する応答値を超えて上記低応答値未満である場合に上記減衰力を最小値に設定し、上記検出信号が上記低応答値を超えた以降は、上記高応答値以下の範囲である場合に、上記演算を実行し、上記高応答値を超えた場合に上記減衰力を最大値に設定することを特徴とする請求項1に記載の可変構造セミアクティブ免震システム。 The computer is set in advance with a low response value and a high response value of the seismic isolation structure, and a response value for determining a control start lower than the low response value, and the detection signal from the sensor starts the control start. The damping force is set to the maximum value when the response value is less than the low response value, and the damping force is set to the minimum value when the response value is less than the low response value. After the detection signal exceeds the low response value, the calculation is performed when the detection signal is within the range of the high response value. When the detection signal exceeds the high response value, the damping force is set to the maximum value. The variable structure semi-active seismic isolation system according to claim 1. 可変減衰ダンパは、上記減衰力が低減衰力と高減衰力との2段階に可変であることを特徴とする請求項1または2に記載の可変構造セミアクティブ免震システム。 The variable damping damper is a variable structure semi-active seismic isolation system according to claim 1 or 2, wherein the damping force is variable in two stages of a low damping force and a high damping force . 上記コンピュータは、上記上部構造と下部構造との相対速度または上記可変減衰ダンパの応答速度が規定値以下である場合に、上記指令信号を更新しないことを特徴とする請求項1ないし3のいずれかに記載の可変構造セミアクティブ免震システム。 4. The computer according to claim 1 , wherein the command signal is not updated when a relative speed between the upper structure and the lower structure or a response speed of the variable damping damper is equal to or lower than a specified value . Variable structure semi-active seismic isolation system described in 1. 上記可変減衰ダンパは、内部に流体が充填された圧力室を有するシリンダと、このシリンダ内に移動自在に設けられて上記圧力室を2室に仕切るピストンと、このピストンと一体化されて上記シリンダから延出するロッドと、上記2室間を連通させる連通経路とを備えた流体圧ダンパであり、かつ上記切換手段は、上記連通経路に配設されて当該連通経路を開閉させることにより上記減衰力を切り換える電磁弁であるとともに、
上記制御手段は、上記電磁弁に供給される駆動電流値を検出する検出手段と、この検出手段により検出された電流値が異常か否かを判定する判定手段を備えることを特徴とする請求項1ないし4のいずれかに記載の可変構造セミアクティブ免震システム。
The variable damping damper includes a cylinder having a pressure chamber filled with a fluid therein, a piston that is movably provided in the cylinder and divides the pressure chamber into two chambers, and the piston integrated with the piston. A fluid pressure damper including a rod extending from the two chambers and a communication path for communicating between the two chambers, and the switching means is disposed in the communication path and opens and closes the communication path. A solenoid valve that switches the force,
The said control means is provided with the detection means which detects the drive current value supplied to the said solenoid valve, and the determination means which determines whether the electric current value detected by this detection means is abnormal. The variable structure semi-active seismic isolation system according to any one of 1 to 4.
上記判定手段は、上記電磁弁を駆動する上記指令信号に対する上記駆動電流値がしきい値未満の場合、および上記電磁弁の駆動を停止する上記指令信号に対する上記駆動電流値がしきい値以上の場合に上記異常と判断する判断手段と、この判断手段が異常と判断した上記電磁弁の駆動電流を遮断する電流遮断手段とを備えることを特徴とする請求項に記載の可変構造セミアクティブ免震システム。 The determination means is configured such that the drive current value for the command signal for driving the electromagnetic valve is less than a threshold value, and the drive current value for the command signal for stopping the drive of the electromagnetic valve is greater than or equal to a threshold value. 6. The variable structure semi-active immunity according to claim 5 , further comprising: a determination unit that determines that the abnormality is abnormal, and a current interruption unit that interrupts a drive current of the solenoid valve that is determined to be abnormal by the determination unit. Seismic system.
JP2005102494A 2005-03-31 2005-03-31 Variable structure semi-active seismic isolation system Active JP4747336B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2005102494A JP4747336B2 (en) 2005-03-31 2005-03-31 Variable structure semi-active seismic isolation system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2005102494A JP4747336B2 (en) 2005-03-31 2005-03-31 Variable structure semi-active seismic isolation system

Publications (2)

Publication Number Publication Date
JP2006283826A JP2006283826A (en) 2006-10-19
JP4747336B2 true JP4747336B2 (en) 2011-08-17

Family

ID=37405965

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2005102494A Active JP4747336B2 (en) 2005-03-31 2005-03-31 Variable structure semi-active seismic isolation system

Country Status (1)

Country Link
JP (1) JP4747336B2 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101892874B1 (en) * 2017-11-28 2018-08-28 권정수 seismic retrofitting structure of concrete structure
KR101921837B1 (en) * 2017-11-30 2019-02-13 신도식 seismic retrofitting structure of concrete structure

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5104238B2 (en) * 2007-11-12 2012-12-19 株式会社大林組 Isolation device
JP5653600B2 (en) * 2009-08-05 2015-01-14 清水建設株式会社 Seismic isolation system and control method of seismic isolation system
JP5431185B2 (en) * 2010-01-25 2014-03-05 大成建設株式会社 Control system for variable damping damper in damping structure
JP6515483B2 (en) * 2014-10-20 2019-05-22 住友電気工業株式会社 Optical fiber drawing apparatus and optical fiber drawing method

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2685200B2 (en) * 1988-02-03 1997-12-03 株式会社デンソー Piezo actuator drive
JPH0324340A (en) * 1989-06-20 1991-02-01 Ishikawajima Harima Heavy Ind Co Ltd Control method for active system dynamic vibration reducer
JP2835943B2 (en) * 1995-10-27 1998-12-14 セイコー精機株式会社 Magnetic bearing control device
JP3606044B2 (en) * 1998-04-23 2005-01-05 鹿島建設株式会社 Damping method for structures and damper for damping
JP3467642B2 (en) * 1999-06-29 2003-11-17 学校法人慶應義塾 Seismic isolation system
JP4103540B2 (en) * 2001-10-30 2008-06-18 トヨタ自動車株式会社 Sliding mode controller
JP4109891B2 (en) * 2002-04-19 2008-07-02 キヤノン株式会社 Active vibration control apparatus, exposure apparatus, and device manufacturing method

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101892874B1 (en) * 2017-11-28 2018-08-28 권정수 seismic retrofitting structure of concrete structure
KR101921837B1 (en) * 2017-11-30 2019-02-13 신도식 seismic retrofitting structure of concrete structure

Also Published As

Publication number Publication date
JP2006283826A (en) 2006-10-19

Similar Documents

Publication Publication Date Title
JP4747336B2 (en) Variable structure semi-active seismic isolation system
CN106078750B (en) A kind of legged type robot hydraulic drive unit Active Compliance Control method
Feng et al. Identification and compensation of non-linear friction for a electro-hydraulic system
Guo et al. Parametric adaptive estimation and backstepping control of electro-hydraulic actuator with decayed memory filter
Ding et al. Self-tuning pressure-feedback control by pole placement for vibration reduction of excavator with independent metering fluid power system
Guo et al. Coupled-disturbance-observer-based position tracking control for a cascade electro-hydraulic system
He et al. Adaptive robust dead-zone compensation control of electro-hydraulic servo systems with load disturbance rejection
Dahunsi et al. System identification and neural network based PID control of servo-hydraulic vehicle suspension system
Song et al. RBF neural network based sliding mode control of a lower limb exoskeleton suit
Kong et al. Force control compensation method with variable load stiffness and damping of the hydraulic drive unit force control system
Maneetham et al. Modeling, simulation and control of high speed nonlinear hydraulic servo system
Friedman et al. Development and experimental validation of a new control strategy considering device dynamics for large-scale MR dampers using real-time hybrid simulation
Pagano et al. Modelling and control of a hydraulically actuated shaking table employed for vibration absorber testing
Vazquez et al. Second order sliding mode control of a mobile hydraulic crane
Na et al. Robust aeroelastic control of flapped wing systems using a sliding mode observer
Sørensen et al. Numerical and experimental study of a novel concept for hydraulically controlled negative loads
Jooei et al. Semi-active control of building structures using variable stiffness device and fuzzy logic
Lu et al. Seismic test of least-input-energy control with ground velocity feedback for variable-stiffness isolation systems
Khaje-Karamodin et al. Semi-active control of structures using a neuro-inverse model of MR dampers
Durbha et al. A nonlinear spring model of hydraulic actuator for passive controller design in bilateral tele-operation
Pfeiffer et al. New hydraulic system modelling
Henikl et al. Modeling and simulation of large-scale manipulators with hydraulic actuation
Mutlu et al. System level performance and reliability investigation of hydraulic circuits using physics based models
Conrad et al. It-tools concept for simulation and design of water hydraulic mechatronic test facilities for motion control and operation in environmentally sensitive application areas
Singh et al. Active aeroelastic control using the receptance method

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20080130

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20080313

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A712

Effective date: 20090827

RD03 Notification of appointment of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7423

Effective date: 20090827

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20090827

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20100125

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20100928

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20101124

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: 20110329

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20110427

R150 Certificate of patent or registration of utility model

Ref document number: 4747336

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140527

Year of fee payment: 3

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250