CN103942394A - Amplitude limiting optimizing method based on magneto-rheological damper - Google Patents
Amplitude limiting optimizing method based on magneto-rheological damper Download PDFInfo
- Publication number
- CN103942394A CN103942394A CN201410171379.2A CN201410171379A CN103942394A CN 103942394 A CN103942394 A CN 103942394A CN 201410171379 A CN201410171379 A CN 201410171379A CN 103942394 A CN103942394 A CN 103942394A
- Authority
- CN
- China
- Prior art keywords
- centerdot
- matrix
- damper
- method based
- magneto
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Vibration Prevention Devices (AREA)
- Buildings Adapted To Withstand Abnormal External Influences (AREA)
Abstract
The invention discloses an amplitude limiting optimizing method based on a magneto-rheological damper. The amplitude limiting optimizing method based on the magneto-rheological damper comprises the following steps that a control law is modified into uc=U(f)*H{(fc-f)f} by combining characteristics of the magneto-rheological damper and particularities of an aqueduct structure; the force corresponding to a force and relative speed hysteretic loop upper portion intersection is adopted as characterization, voltages are progressively increased by 0.1v to conduct simulation, and the disperse relation between the force and the voltages is obtained. According to the amplitude limiting optimizing method based on the magneto-rheological damper, the control law is modified by orienting to the characteristics of the magneto-rheological damper and combining the particularities of the aqueduct structure, the force corresponding to the force and relative speed hysteretic loop upper portion intersection is adopted as the characterization, the voltages are progressively increased by 0.1v to conduct the simulation, the disperse relation between the force and the voltages is obtained, the method is simple, operation is convenient, changes are conducted on an original basis, and the problems that in an existing amplitude limiting optimizing method, when certain low-amplitude vibration is caused, large control voltages are still be loaded, the changing amplitude of control force can be large, and local earthquake responses are possibly enlarged are well solved.
Description
Technical field
The invention belongs to MR damper technical field, relate in particular to a kind of method based on MR damper amplitude limit optimum.
Background technology
At present, South to Northern Water Diversion Project of China is the large hydraulic engineering that the present situation for alleviating China North China and Water Resource in Northwest shortage is built, analysis on aqueduct structure is as the important component part in south water to north water delivery engineering, on Route, to apply one of more main crossing, and the water delivery line of work centerline passes through China's Areas of High Earthquake Intensity region, the seismic fortification intensity of the many sections of circuit reaches 8 degree, 49 of the total aqueducts of whole engineering water delivery general main canal, cumulative length reaches 5520m, and flow reaches 500m
3/ s, once there is earthquake, water delivery completely will be interrupted, and it closes on important trunk railways such as capital is wide, capital nine, may cause serious secondary disaster, therefore the quake-resistant safety of water delivery engineering is extremely important.
The amplitude limit optimal algorithm that Dyke proposes for MR damper has at present obtained application comparatively widely, and its control law can be described as:
u
c=u
maxH{(f
c-f)f}
In formula, H{} is Heaviside step function, and this formula meaning is, if the damping force f that MR produces
cequal active controlling force f, keep current voltage constant; If f
cbe less than f, be adjusted to by applying voltage the control coupling expectation value that maximal value provides to impel MR; In other situations, voltage is made as to zero.
But, in the time that maximum voltage is set, there is certain subjectivity in this algorithm, and in working control, only have two grades of voltages adjustable, be null value and maximal value, to cause like this some low amplitude vibration time, still have larger control voltage-drop loading, like this change amplitude of control can be larger, likely can strengthen local seismic response.
Summary of the invention
The object of the embodiment of the present invention is to provide a kind of method based on MR damper amplitude limit optimum, when being intended to solve existing amplitude limit optimal algorithm and causing the vibration of some low amplitude, still there is larger control voltage-drop loading, the change amplitude of control can be larger, likely can strengthen the problem of local seismic response.
The embodiment of the present invention is achieved in that a kind of method based on MR damper amplitude limit optimum, should the method based on MR damper amplitude limit optimum comprise the following steps:
In conjunction with the singularity of MR damper characteristic and analysis on aqueduct structure, control law is revised as:
u
c=U(f)*H{(f
c-f)f}
The employing power masterpiece corresponding with relative velocity hysteretic loop top intersection point is to characterize, and simulates with 0.1v incremental voltage, obtains the discrete relationship of power and voltage.
Further, discrete point is carried out to matching, finally U (f) is expressed as:
U(f)=(af+b)*H{af+b}
In formula, a, b is the Relation Parameters of power and voltage.
Further, the method for this method amplitude limit optimum control based on MR damper amplitude limit optimum is:
In the structure of n degree of freedom, add r controller, the structure half active control system equation of motion under the consistent input of earthquake motion is:
X(t
0)=X
0
In formula, M, C and K are respectively quality, damping and the stiffness matrix of structure n × n dimension; X (t),
with
respectively displacement, speed and the vector acceleration of structure n dimension;
it is earthquake ground motion acceleration time-histories; U (t) is r × 1 dimension control vector; H
0it is n × 1 dimension earthquake ground motion acceleration location matrix; B
sthat n × r ties up half ACTIVE CONTROL damper position matrix;
Introduce state vector
The equation of motion of structure is expressed as state equation:
Z(t
0)=Z
0
Wherein,
C=[0 I], Y is observation output, for absolute acceleration is using as feedback, v (t) is for measuring noise;
The performance objective function of define system is:
Wherein, Q and R are weight matrix, make performance objective function get minimum value, can be in the hope of feedback gain G=R
-1b
tp is r × 2n dimension feedback of status gain matrix; P is that 2n × 2n ties up matrix, can be by Riccati matrix algebra equation solution;
Owing to being fed back to absolute acceleration, so need to introduce state estimation vector
replace state vector Z, make
keep optimum for output feedback problem, state estimation vector
produced by Kalman wave filter:
Formula median filter gain L can be expressed as:
L=P
1C
TV
-1
V=E[v (t) v
t(t)] for measuring noise covariance matrix, E[] expression mathematical expectation, P
1solution for following formula:
AP
1+P
1A
T+N-P
1C
TV
-1CP
1=0
In formula, N=E[w (t) w
t(t)] be input noise covariance matrix, Kalman wave filter makes the progressive covariance of state estimation error
reach minimum.
Method based on MR damper amplitude limit optimum provided by the invention, by for MR damper characteristic, in conjunction with the singularity of analysis on aqueduct structure, control law is made to improvement, the employing power masterpiece corresponding with relative velocity hysteretic loop top intersection point is sign, simulate with 0.1v incremental voltage, obtain the discrete relationship of power and voltage, method is simple, easy to operate, make a change on the original basis, when preferably resolving existing amplitude limit optimal algorithm and causing the vibration of some low amplitude, still there is larger control voltage-drop loading, the change amplitude of control can be larger, likely can strengthen the problem of local seismic response.
Brief description of the drawings
Fig. 1 is the method flow diagram based on MR damper amplitude limit optimum that the embodiment of the present invention provides;
Fig. 2 is the control that provides of the embodiment of the present invention and the discrete relationship schematic diagram of voltage.
It is lower across span centre and the right time-histories across span centre response schematic diagram in the input of ElcentroNS seismic event that Fig. 3, Fig. 4 are respectively analysis on aqueduct structure.
Embodiment
In order to make object of the present invention, technical scheme and advantage clearer, below in conjunction with embodiment, the present invention is further elaborated.Should be appreciated that specific embodiment described herein, only in order to explain the present invention, is not intended to limit the present invention.
Below in conjunction with drawings and the specific embodiments, application principle of the present invention is further described.
As shown in Figure 1, the method based on MR damper amplitude limit optimum of the embodiment of the present invention comprises the following steps:
S101: in conjunction with the singularity of MR damper characteristic and analysis on aqueduct structure, change control rule;
S102: the employing power masterpiece corresponding with relative velocity hysteretic loop top intersection point is to characterize, and simulates with 0.1v incremental voltage, obtains the discrete relationship of power and voltage.
Concrete steps of the present invention are:
The first step, in conjunction with the singularity of MR damper characteristic and analysis on aqueduct structure, changes control law into:
u
c=U(f)*H{(f
c-f)f}
Wherein, U (f) is for characterizing the continuous function of voltage and power; H{} is Heaviside step function.
Second step, the employing power masterpiece corresponding with relative velocity hysteretic loop top intersection point is to characterize, and simulates with 0.1v incremental voltage, obtains the discrete relationship of power and voltage, as shown in Figure 2.
Active controlling force, by semi-active control law u=k (f
c, calculating f), draws the result of semi-active control law, result outputs to MR damper, finally exports half active controlling force.
As can be seen from Figure 2, linear relationship is comparatively obvious, and discrete point is carried out to matching, finally U (f) is expressed as:
U(f)=(af+b)*H{af+b}
In formula, a, b is the Relation Parameters of power and voltage.
The method of amplitude limit optimum control of the present invention is:
In the structure of n degree of freedom, add r controller, the structure half active control system equation of motion under the consistent input of earthquake motion is:
X(t
0)=X
0
\*MERGEFORMAT (1)
In formula, M, C and K are respectively quality, damping and the stiffness matrix of structure n × n dimension; X (t),
with
respectively displacement, speed and the vector acceleration of structure n dimension;
it is earthquake ground motion acceleration time-histories; U (t) is r × 1 dimension control vector; H
0it is n × 1 dimension earthquake ground motion acceleration location matrix; B
sthat n × r ties up half ACTIVE CONTROL damper position matrix.
Introduce state vector
The equation of motion of structure can be expressed as state equation:
Z(t
0)=Z
0 \*MERGEFORMAT (2)
Wherein,
C=[0 I], Y is observation output, be absolute acceleration using as feedback herein, v (t) is measurement noise.
The performance objective function of define system is
Wherein, Q and R are weight matrix.Make performance objective function get minimum value, can be in the hope of feedback gain G=R
-1b
tp is r × 2n dimension feedback of status gain matrix; P is that 2n × 2n ties up matrix, can be by Riccati matrix algebra equation solution.
Owing to being fed back to absolute acceleration, so need to introduce state estimation vector
replace state vector Z, make
keep optimum for output feedback problem, state estimation vector
produced by Kalman wave filter:
\*MERGEFORMAT (4)
Formula median filter gain L can be expressed as
L=P
1C
TV
-1 \*MERGEFORMAT (5)
V=E[v (t) v
t(t)] for measuring noise covariance matrix, E[] expression mathematical expectation, P
1for the solution of following formula
AP
1+P
1A
T+N-P
1C
TV
-1CP
1=0 \*MERGEFORMAT (6)
In formula, N=E[w (t) w
t(t)] be input noise covariance matrix.Kalman wave filter makes the progressive covariance of state estimation error
reach minimum.
As shown in Figure 3 and Figure 4, be the damping effect schematic diagram of the embodiment of the present invention, can find out by Fig. 3 and Fig. 4, the present invention has good damping performance.
It is lower across span centre and the response of the right time-histories across span centre in the input of ElcentroNS seismic event that Fig. 3 and Fig. 4 are respectively analysis on aqueduct structure, in figure, directions X coordinate time is time shaft, unit is second, Y-direction coordinate axis Displacement is displacement response, unit is rice, Z direction coordinate Velocity is speed responsive, unit is meter per second, in figure, the representative of the dotted portion of curve is the time-histories response of aqueduct in the time not adding control device, in figure, the representative of the solid line part of curve is the time-histories response of analysis on aqueduct structure under this control program, from solid line and dotted line to recently, there is following characteristics:
1, within the scope of seismic event input maximum value (two time periods of 5-10s and 15-20s), solid line (control program) numerical value will significantly be less than dotted line (without control) numerical value, illustrates that control program has good control effect to the peak value of seismic response.
2, in the less scope of seismic event input value, the value envelope of solid line, within dotted line numerical value, illustrates that control program also has certain effect, controls effect steady, does not occur enlarge-effect.
It is lower across span centre and the response of the right time-histories across span centre in the input of ElcentroNS seismic event that Fig. 3 and Fig. 4 are respectively analysis on aqueduct structure, in figure, directions X coordinate time is time shaft, unit is second, Y-direction coordinate axis Displacement is displacement response, unit is rice, Z direction coordinate Velocity is speed responsive, unit is meter per second, in figure, the representative of the dotted portion of curve is the time-histories response of aqueduct in the time not adding control device, in figure, the representative of the solid line part of curve is the time-histories response of analysis on aqueduct structure under this control program, from solid line and dotted line to recently, there is following characteristics:
1. within the scope of seismic event input maximum value (two time periods of 5-10s and 15-20s), solid line (control program) numerical value will significantly be less than dotted line (without control) numerical value, illustrates that control program has good control effect to the peak value of seismic response.
In the less scope of seismic event input value, the value envelope of solid line, within dotted line numerical value, illustrates that control program also has certain effect, controls effect steady, does not occur enlarge-effect.
The foregoing is only preferred embodiment of the present invention, not in order to limit the present invention, all any amendments of doing within the spirit and principles in the present invention, be equal to and replace and improvement etc., within all should being included in protection scope of the present invention.
Claims (3)
1. the method based on MR damper amplitude limit optimum, is characterized in that, should the method based on MR damper amplitude limit optimum comprise the following steps:
In conjunction with the singularity of MR damper characteristic and analysis on aqueduct structure, control law is revised as:
u
c=U(f)*H{(f
c-f)f}
The employing power masterpiece corresponding with relative velocity hysteretic loop top intersection point is to characterize, and simulates with 0.1v incremental voltage, obtains the discrete relationship of power and voltage.
2. the method based on MR damper amplitude limit optimum as claimed in claim 1, is characterized in that, discrete point is carried out to matching, finally U (f) is expressed as:
U(f)=(af+b)*H{af+b}
In formula, a, b is the Relation Parameters of power and voltage.
3. the method based on MR damper amplitude limit optimum as claimed in claim 1, is characterized in that, the method for being somebody's turn to do the method amplitude limit optimum control based on MR damper amplitude limit optimum is:
In the structure of n degree of freedom, add r controller, the structure half active control system equation of motion under the consistent input of earthquake motion is:
X(t
0)=X
0
In formula, M, C and K are respectively quality, damping and the stiffness matrix of structure n × n dimension; X (t),
with
respectively displacement, speed and the vector acceleration of structure n dimension;
it is earthquake ground motion acceleration time-histories; U (t) is r × 1 dimension control vector; H
0it is n × 1 dimension earthquake ground motion acceleration location matrix; B
sthat n × r ties up half ACTIVE CONTROL damper position matrix;
Introduce state vector
The equation of motion of structure is expressed as state equation:
Z(t
0)=Z
0
Wherein,
C=[0 I], Y is observation output, for absolute acceleration is using as feedback, v (t) is for measuring noise;
The performance objective function of define system is:
Wherein, Q and R are weight matrix, make performance objective function get minimum value, can be in the hope of feedback gain G=R
-1b
tp is r × 2n dimension feedback of status gain matrix; P is that 2n × 2n ties up matrix, can be by Riccati matrix algebra equation solution;
Owing to being fed back to absolute acceleration, so need to introduce state estimation vector
replace state vector Z, make
keep optimum for output feedback problem, state estimation vector
produced by Kalman wave filter:
Formula median filter gain L can be expressed as:
L=P
1C
TV
-1
V=E[v (t) v
t(t)] for measuring noise covariance matrix, E[] expression mathematical expectation, P
1solution for following formula:
AP
1+P
1A
T+N-P
1C
TV
-1CP
1=0
In formula, N=E[w (t) w
t(t)] be input noise covariance matrix, Kalman wave filter makes the progressive covariance of state estimation error
reach minimum.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410171379.2A CN103942394B (en) | 2014-04-25 | 2014-04-25 | Amplitude limiting optimizing method based on magneto-rheological damper |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410171379.2A CN103942394B (en) | 2014-04-25 | 2014-04-25 | Amplitude limiting optimizing method based on magneto-rheological damper |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103942394A true CN103942394A (en) | 2014-07-23 |
CN103942394B CN103942394B (en) | 2017-02-15 |
Family
ID=51190062
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410171379.2A Expired - Fee Related CN103942394B (en) | 2014-04-25 | 2014-04-25 | Amplitude limiting optimizing method based on magneto-rheological damper |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103942394B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105260530A (en) * | 2015-09-29 | 2016-01-20 | 合肥工业大学 | Modeling method for magneto-rheological damper inverse model and application thereof |
CN105930938A (en) * | 2016-05-10 | 2016-09-07 | 重庆大学 | Method of load shedding of floating wind turbine generator system based on semi-active structure control of magneto rheological damper |
CN111338208A (en) * | 2020-03-02 | 2020-06-26 | 南京林业大学 | Vehicle side-tipping and smoothness coordination control method |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100417773C (en) * | 2006-06-07 | 2008-09-10 | 东南大学 | High lossing energy composite magnetorheopectic damping device |
CN102767589B (en) * | 2012-07-02 | 2014-07-16 | 江苏科技大学 | Magneto-rheological damper adjustable in damping channel height |
-
2014
- 2014-04-25 CN CN201410171379.2A patent/CN103942394B/en not_active Expired - Fee Related
Non-Patent Citations (4)
Title |
---|
DYKE S J, SPENCER JR B F, SAIN M K, ET AL.: "Modeling and control of magnetorheological dampers for seismic response reduction", 《SMART MATERIALS AND STRUCTURES》 * |
王博,任克彬,徐建国,黄亮: "大型渡槽结构半主动变阻尼控制研究", 《水利与建筑工程学报》 * |
郑明燕,杨世浩,李正农: "大型渡槽智能隔震数值模拟", 《武汉大学学报(工学版)》 * |
黄亮,侯玉洁,徐建国: "磁流变阻尼器在渡槽抗震中的应用研究", 《磁性材料及器件》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105260530A (en) * | 2015-09-29 | 2016-01-20 | 合肥工业大学 | Modeling method for magneto-rheological damper inverse model and application thereof |
CN105260530B (en) * | 2015-09-29 | 2018-03-09 | 合肥工业大学 | A kind of modeling method of MR damper inversion model and its application |
CN105930938A (en) * | 2016-05-10 | 2016-09-07 | 重庆大学 | Method of load shedding of floating wind turbine generator system based on semi-active structure control of magneto rheological damper |
CN105930938B (en) * | 2016-05-10 | 2019-06-21 | 重庆大学 | Floatation type blower off-load method based on MR damper semi-active structural co ntrol |
CN111338208A (en) * | 2020-03-02 | 2020-06-26 | 南京林业大学 | Vehicle side-tipping and smoothness coordination control method |
CN111338208B (en) * | 2020-03-02 | 2022-05-20 | 南京林业大学 | Vehicle side-tipping and smoothness coordination control method |
Also Published As
Publication number | Publication date |
---|---|
CN103942394B (en) | 2017-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106934185B (en) | A kind of multiple dimensioned flow simulating method of the fluid structurecoupling of elastic fluid | |
Silva et al. | Stochastic analysis of nonlinear wave energy converters via statistical linearization | |
CN107506564A (en) | Consider stake Longitudinal vibration analysis method in vertical fluctuation effect radial direction heterogeneous soil | |
CN104573274B (en) | Structural finite element model updating method based on displacement time-histories area under vehicular load | |
Domaneschi | Feasible control solutions of the ASCE benchmark cable‐stayed bridge | |
CN104268791A (en) | Health assessment method for 500 kV high-voltage power transmission line in mountain land occurrence environment | |
CN103942394A (en) | Amplitude limiting optimizing method based on magneto-rheological damper | |
CN107121932A (en) | Motor servo system error symbol integrates Robust Adaptive Control method | |
Domany et al. | Dynamic responses and mitigation of limit cycle oscillations in Van der Pol–Duffing oscillator with nonlinear energy sink | |
CN103399996A (en) | Multi-target topological optimization design method for flexible mechanism for fast servo tool rest | |
CN108446460A (en) | A kind of radial direction is heterogeneous, the layered soil body pile foundation Longitudinal vibration analysis method in longitudinal direction | |
Ke et al. | Stability and reinforcement analysis of superlarge exhaust cooling towers based on a wind tunnel test | |
CN106054670B (en) | A kind of super mangneto driver hysteresis modeling method based on time lag | |
Kim et al. | On the use of coupled wind, wave, and current fields in the simulation of loads on bottom-supported offshore wind turbines during hurricanes: March 2012-September 2015 | |
CN102176120B (en) | Method for optimally controlling magnetorheological damping of random seismic excitation system | |
Macchelli | Boundary energy shaping of linear distributed port-Hamiltonian systems | |
CN105701279B (en) | Non-classical structural dynamic response frequency domain method | |
da Silva et al. | Wave energy converter modelling, control, and power take-off design | |
CN104020669A (en) | Semi-active shock absorption controlling device and method for large aqueduct under longitudinal shock input | |
Gao et al. | Dynamic evolution characteristics of a fractional order hydropower station system | |
CN105162377B (en) | A method of establishing generator excitation controller energy function | |
CN104358327B (en) | Damping method of random-rigidity eccentric structure | |
Roy et al. | Small-signal stability assessment of active distribution networks with dynamic loads | |
CN103926836A (en) | Shock absorption control method of aqueduct structure under action of random loads | |
CN101853323B (en) | Modeling method for acting force of full-automatic pulsator washing machine suspension system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170215 Termination date: 20170425 |