CN105320057A - Synchronous movement control method of double vibration tables on the basis of coordinate-transformation matrix - Google Patents
Synchronous movement control method of double vibration tables on the basis of coordinate-transformation matrix Download PDFInfo
- Publication number
- CN105320057A CN105320057A CN201510282045.7A CN201510282045A CN105320057A CN 105320057 A CN105320057 A CN 105320057A CN 201510282045 A CN201510282045 A CN 201510282045A CN 105320057 A CN105320057 A CN 105320057A
- Authority
- CN
- China
- Prior art keywords
- transition matrix
- freedom
- output
- input
- coordinate
- 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
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/18—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
- G05B19/19—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by positioning or contouring control systems, e.g. to control position from one programmed point to another or to control movement along a programmed continuous path
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M7/00—Vibration-testing of structures; Shock-testing of structures
- G01M7/02—Vibration-testing by means of a shake table
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Human Computer Interaction (AREA)
- Manufacturing & Machinery (AREA)
- Automation & Control Theory (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention discloses a synchronous movement control method of double vibration tables on the basis of a coordinate-transformation matrix. The synchronous movement control method comprises: converting linear measurement of two control sensors on the double vibration tables to corresponding two degrees of freedom such as translation and rotation, and inputting the coordinate-transformation matrix; converting signals outputted by a multiple-input-multiple-output system to driving signals actually inputted to the power amplifiers of the double vibration tables, and outputting the coordinate-transformation matrix; and synchronously controlling the two-degree of freedom motion, specifically, performing closed-loop control according to the translation degree of freedom reference test magnitude, the rotation degree of freedom reference test magnitude, the transfer function of the system and the measured two-degree of freedom response through the multiple-input-multiple-output system, and realizing the synchronous movement of the double vibration tables. The synchronous movement control method of double vibration tables on the basis of a coordinate-transformation matrix has high control precision, and may preferably solve the technical question of synchronous movement control of the double vibration tables.
Description
Technical field
The invention belongs to spacecraft dynamic environmental test field, be specifically related to a kind of Double earthquakes Synchronous motion control method based on coordinate conversion matrix.
Background technology
Double earthquakes is synchronized with the movement and encourages pilot system to mainly contain two kinds of control modes at present.A kind of is power synchronisation control means based on Double earthquakes moving-coil current feedback; A kind of is standard square formation control method based on two reference mark on Double earthquakes.Double earthquakes utilizes phase controller to regulate the drive singal inputing to power amplifier according to the phase differential of two power amplifier feedback signals based on the power synchronisation control means of moving-coil current feedback usually, coordinates two shaking table synchronized movement; And utilize multi-input multi-output control system directly to carry out the control of standard square formation based on the reference test condition that the standard square formation controlling party rule at two reference mark on Double earthquakes is identical according to two reference mark, make two shaking table synchronized movement.The current power synchronisation control means based on Double earthquakes moving-coil current feedback needs single shaft vibrating controller to carry out two platform synchronized movement control of shaking in conjunction with phase controller, system complex; And all there is the not high problem of control accuracy with the standard square formation control method based on two reference mark on Double earthquakes in the method.
Summary of the invention
The technical problem to be solved in the present invention is to provide a kind of Double earthquakes Synchronous motion control method based on coordinate conversion matrix to make control accuracy high, and can solve Double earthquakes Synchronous motion control technical matters preferably.
In order to solve the problems of the technologies described above, the invention provides a kind of Double earthquakes synchronisation control means based on coordinate conversion matrix, comprising the following steps:
(1) input coordinate transition matrix is built, described structure input coordinate transition matrix step comprises according to the distance determination input coordinate transition matrix on Double earthquakes between two control sensors, and linear measurement Double earthquakes controlling sensor is converted to corresponding translation and rotation two degree of freedom by acting as of input coordinate transition matrix;
(2) output coordinate transition matrix is built, described structure output coordinate transition matrix step comprise according to Double earthquakes excitation between distance determination output coordinate transition matrix, output coordinate transition matrix act as by multi-input multi-output control system export signal be converted to the actual drive singal inputing to Double earthquakes power amplifier;
(3) two degree-of-freedom motion synchro control, described two degree-of-freedom motion synchro control step is specially the white noise random signal that multi-input multi-output control system first sends the low magnitude of two-way, and this white noise random signal enters power amplifier and drives Double earthquakes to move after output coordinate transition matrix; On Double earthquakes, the response of two measurement points enters multi-input multi-output control system after input coordinate transition matrix, and control system calculates the transport function of whole system according to input and output signal; Multi-input multi-output control system calculates initial two-way drive singal according to the translational degree of freedom reference test magnitude arranged, rotational freedom reference test magnitude, system transter subsequently, and this two-way drive singal enters power amplifier and drives Double earthquakes to move after output coordinate transition matrix; Multi-input multi-output control system is according to two degrees of freedom response level and the translation of measuring and rotates two degrees of freedom reference test magnitude and compare and the error of calculation, revises the two-way drive singal of transmission in real time, thus realizes the control of Double earthquakes synchronized movement.
Double earthquakes Synchronous motion control method based on coordinate conversion matrix provided by the invention is converted to corresponding translation according to the control linear sensor measurement on two shaking tables and rotates two degree of freedom, test magnitude and the standard side of translational degree of freedom Control experiment magnitude of shaking is identical, the test magnitude of rotational freedom is then set to the lower numerical value of a relative translational movement degree of freedom, translational degree of freedom is realized by multi-input multi-output control system, and rotational freedom is suppressed by multi-input multi-output control system, thus realize Double earthquakes synchronized movement.Control method of the present invention controls to compare with direct standard square formation with two platform current feedback synchro control of shaking, and control accuracy is higher, solves Double earthquakes Synchronous motion control technical matters preferably.
Accompanying drawing explanation
Fig. 1 is Double earthquakes synchro control of the present invention response geometric layout figure;
Wherein, 1-first acceleration transducer, 2-second acceleration transducer, the two platform geometric center of shaking of 3-.
Fig. 2 is Double earthquakes synchro control of the present invention excitation geometric layout figure;
Wherein, 4-first shaking table excitation, the excitation of 5-second shaking table.
Fig. 3 is input coordinate transition matrix design of graphics of the present invention.
Fig. 4 is output coordinate transition matrix design of graphics of the present invention.
Fig. 5 is Double earthquakes two degree-of-freedom motion synchro control process flow diagram of the present invention.
Embodiment
Below in conjunction with accompanying drawing, a kind of Double earthquakes Synchronous motion control method based on coordinate conversion matrix of the present invention is further described.
Fig. 1 is Double earthquakes synchro control of the present invention response geometric layout figure.As shown in Figure 1, the distance between the first acceleration transducer 1 shaking table placed and the second acceleration transducer 2 is 2r.First acceleration transducer 1 and the second acceleration transducer 2 are r from the distance of Double earthquakes geometric center 3.How the linear acceleration of two Vertical dimension is measured the translational degree of freedom T and the rotational freedom R that are converted to test specimen by input coordinate transition matrix by describing below.
For translation, the first acceleration transducer 1, second acceleration transducer 2 need move in same direction, mathematic(al) mean method can be adopted to calculate, that is:
T
control=0.5A
1+0.5A
2..................................(1)
Wherein A1 represents the response of the first acceleration transducer; A2 represents the response of the second acceleration transducer; T
controlrepresent that translational degree of freedom controls magnitude.
For rotation, the first acceleration transducer 1, second acceleration transducer 2 needs direction of motion contrary, also can adopt mathematic(al) mean method to calculate, need consider the impact of acceleration of gravity constant G when linear acceleration being converted to rotary acceleration.That is:
Convert upper two formula to matrix form, then have:
Input coordinate transition matrix is expressed as follows:
Fig. 2 is Double earthquakes synchro control of the present invention excitation geometric layout figure.As shown in Figure 2, for output coordinate transition matrix, by the method for observing, output degree of freedom can be converted to linear incentive signal.Suppose the first shaking table excitation 4, second shaking table excitation 5 apart from the distance of two platform geometric center 3 of shaking be r ~.Utilize input transition matrix derivation method, can show that the relation that excitation degree of freedom and linear incentive export is as follows:
Wherein D1 represents the drive singal of the first shaking table excitation 4; D2 represents the drive singal of the second shaking table excitation 5; T
driverepresent that translational degree of freedom drives magnitude; R
driverepresent that rotational freedom drives magnitude.
Therefore the pumping signal of reality is by obtaining above formula matrix inversion, that is:
Finally show that output linearity pumping signal is as follows with the relation of excitation degree of freedom:
Output coordinate expression matrix is as follows:
Double earthquakes Synchronous motion control method based on coordinate conversion matrix comprises the steps:
(1) input coordinate transition matrix is built.Particularly, build input coordinate transition matrix to comprise according to the distance determination input coordinate transition matrix (see formula 4) (as shown in Figure 3) on Double earthquakes between two control sensors.Linear measurement Double earthquakes controlling sensor is converted to corresponding translation and rotation two degree of freedom by acting as of input coordinate transition matrix.
(2) output coordinate transition matrix is built.Particularly, build output coordinate transition matrix to comprise according to Double earthquakes excitation between distance determination output coordinate transition matrix (see formula 8) (as shown in Figure 4).The signal act as multi-input multi-output control system exports of output coordinate transition matrix is converted to the actual drive singal inputing to Double earthquakes power amplifier.
(3) two degree-of-freedom motion synchro control.Particularly, two degree-of-freedom motion synchro control is the white noise random signal that multi-input multi-output control system first sends the low magnitude of two-way, and this white noise random signal enters power amplifier and drives Double earthquakes to move after output coordinate transition matrix; On Double earthquakes, the response of two measurement points enters multi-input multi-output control system after input coordinate transition matrix, and control system calculates the transport function of whole system according to input and output signal; According to the translational degree of freedom reference test magnitude arranged, rotational freedom reference test magnitude, (rotational freedom reference test magnitude requires for translational degree of freedom test magnitude 10 to multi-input multi-output control system subsequently
-3doubly following), system transter calculates initial two-way output signal, this two-way output signal enters power amplifier and drives Double earthquakes to move after output coordinate transition matrix; Multi-input multi-output control system is according to two degrees of freedom response level and the translation of measuring and rotate two degrees of freedom reference test magnitude and compare the also error of calculation, the real-time two-way drive singal revising transmission, thus realize Double earthquakes synchronized movement control (as shown in Figure 5).
Although give detailed description and explanation to the specific embodiment of the present invention above; but what should indicate is; we can carry out various equivalence according to conception of the present invention to above-mentioned embodiment and change and amendment; its function produced do not exceed that instructions and accompanying drawing contain yet spiritual time, all should within protection scope of the present invention.
Claims (5)
1., based on a Double earthquakes Synchronous motion control method for coordinate conversion matrix, comprise following three steps:
(1) input coordinate transition matrix is built, described structure input coordinate transition matrix step comprises according to the distance determination input coordinate transition matrix on Double earthquakes between two control sensors, and linear measurement Double earthquakes controlling sensor is converted to corresponding translation and rotation two degree of freedom by acting as of input coordinate transition matrix;
(2) output coordinate transition matrix is built, described structure output coordinate transition matrix step comprise according to Double earthquakes excitation between distance determination output coordinate transition matrix, output coordinate transition matrix act as by multi-input multi-output control system export signal be converted to the actual drive singal inputing to Double earthquakes power amplifier;
(3) two degree-of-freedom motion synchro control, described two degree-of-freedom motion synchro control step is specially the white noise random signal that multi-input multi-output control system first sends the low magnitude of two-way, and this white noise random signal enters power amplifier and drives Double earthquakes to move after output coordinate transition matrix; On Double earthquakes, the response of two measurement points enters multi-input multi-output control system after input coordinate transition matrix, and control system calculates the transport function of whole system according to input and output signal; Multi-input multi-output control system calculates initial two-way output signal according to the translational degree of freedom reference test magnitude arranged, rotational freedom reference test magnitude, system transter subsequently, and this two-way output signal enters power amplifier and drives Double earthquakes to move after output coordinate transition matrix; Multi-input multi-output control system is according to two degrees of freedom response level and the translation of measuring and rotates two degrees of freedom reference test magnitude and compare and the error of calculation, revises the two-way drive singal of transmission in real time, thus realizes the control of Double earthquakes synchronized movement.
2., as claimed in claim 1 based on the Double earthquakes Synchronous motion control method of coordinate conversion matrix, it is characterized in that, described input coordinate transition matrix and described output coordinate transition matrix are the square formation of 2 × 2.
3. as claimed in claim 1 based on the Double earthquakes Synchronous motion control method of coordinate conversion matrix, it is characterized in that, rotational freedom test magnitude requires for translational degree of freedom test magnitude 10
-3doubly.
4., as claimed in claim 1 based on the Double earthquakes Synchronous motion control method of coordinate conversion matrix, it is characterized in that, described input coordinate transition matrix is:
5., as claimed in claim 1 based on the Double earthquakes Synchronous motion control method of coordinate conversion matrix, it is characterized in that, described output coordinate transition matrix is:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510282045.7A CN105320057B (en) | 2015-05-28 | 2015-05-28 | Double earthquakes Synchronous motion control method based on coordinate conversion matrix |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510282045.7A CN105320057B (en) | 2015-05-28 | 2015-05-28 | Double earthquakes Synchronous motion control method based on coordinate conversion matrix |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105320057A true CN105320057A (en) | 2016-02-10 |
CN105320057B CN105320057B (en) | 2018-04-06 |
Family
ID=55247657
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510282045.7A Expired - Fee Related CN105320057B (en) | 2015-05-28 | 2015-05-28 | Double earthquakes Synchronous motion control method based on coordinate conversion matrix |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105320057B (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106248327A (en) * | 2016-07-11 | 2016-12-21 | 北京强度环境研究所 | A kind of firing box Transport Vibration assay device and method |
CN107543671A (en) * | 2016-07-15 | 2018-01-05 | 北京卫星环境工程研究所 | More and exciting movement controller and digit phase isochronous controller matching process |
CN108254145A (en) * | 2017-12-29 | 2018-07-06 | 苏州东菱智能减振降噪技术有限公司 | A kind of control method for realizing multi shaker synchronous vibration |
CN112444367A (en) * | 2020-12-18 | 2021-03-05 | 中国工程物理研究院总体工程研究所 | Multi-vibration-table parallel-pushing single-shaft vibration test control method |
CN112629833A (en) * | 2019-09-25 | 2021-04-09 | 上海汽车集团股份有限公司 | Load acquisition method and device |
CN114646440A (en) * | 2022-03-16 | 2022-06-21 | 北京卫星环境工程研究所 | Three-axis six-degree-of-freedom vibration control input and output coordinate transformation matrix construction method |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005121371A (en) * | 2003-10-14 | 2005-05-12 | National Institute Of Advanced Industrial & Technology | Dynamic matrix sensitivity measuring method for force sensor, and instrument therefor |
CN101173876A (en) * | 2007-11-30 | 2008-05-07 | 哈尔滨工业大学 | Control method for three-axis six-freedom hydraulic vibration table based on kinematics |
CN101196429A (en) * | 2007-12-07 | 2008-06-11 | 哈尔滨工业大学 | Internal force decoupling control method for driving vibration table with redundant 6 freedom |
CN104133479A (en) * | 2014-08-08 | 2014-11-05 | 上海新跃仪表厂 | Test system and method for simulating flexible satellite three-axis attitude coupling movement with single-axis air bearing table |
-
2015
- 2015-05-28 CN CN201510282045.7A patent/CN105320057B/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005121371A (en) * | 2003-10-14 | 2005-05-12 | National Institute Of Advanced Industrial & Technology | Dynamic matrix sensitivity measuring method for force sensor, and instrument therefor |
CN101173876A (en) * | 2007-11-30 | 2008-05-07 | 哈尔滨工业大学 | Control method for three-axis six-freedom hydraulic vibration table based on kinematics |
CN101196429A (en) * | 2007-12-07 | 2008-06-11 | 哈尔滨工业大学 | Internal force decoupling control method for driving vibration table with redundant 6 freedom |
CN104133479A (en) * | 2014-08-08 | 2014-11-05 | 上海新跃仪表厂 | Test system and method for simulating flexible satellite three-axis attitude coupling movement with single-axis air bearing table |
Non-Patent Citations (3)
Title |
---|
CHEN G,ZHUANG H,DU X等: "Analysis of large-scale shaking table test of dynamic soil-subway station interaction", 《EARTHQUAKE ENGINEERING AND ENGINEERING VIBRATION》 * |
樊世超等: "多维动力学环境模拟试验技术研究", 《航天器环境工程》 * |
邱汉平等: "双振动台正弦振动同步控制新方法研究", 《航天器环境工程》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106248327A (en) * | 2016-07-11 | 2016-12-21 | 北京强度环境研究所 | A kind of firing box Transport Vibration assay device and method |
CN107543671A (en) * | 2016-07-15 | 2018-01-05 | 北京卫星环境工程研究所 | More and exciting movement controller and digit phase isochronous controller matching process |
CN107543671B (en) * | 2016-07-15 | 2019-02-19 | 北京卫星环境工程研究所 | More and exciting movement controller and digit phase isochronous controller matching process |
CN108254145A (en) * | 2017-12-29 | 2018-07-06 | 苏州东菱智能减振降噪技术有限公司 | A kind of control method for realizing multi shaker synchronous vibration |
CN108254145B (en) * | 2017-12-29 | 2020-06-26 | 苏州东菱智能减振降噪技术有限公司 | Control method for realizing synchronous vibration of multiple vibration tables |
CN112629833A (en) * | 2019-09-25 | 2021-04-09 | 上海汽车集团股份有限公司 | Load acquisition method and device |
CN112629833B (en) * | 2019-09-25 | 2024-04-19 | 上海汽车集团股份有限公司 | Load acquisition method and device |
CN112444367A (en) * | 2020-12-18 | 2021-03-05 | 中国工程物理研究院总体工程研究所 | Multi-vibration-table parallel-pushing single-shaft vibration test control method |
CN114646440A (en) * | 2022-03-16 | 2022-06-21 | 北京卫星环境工程研究所 | Three-axis six-degree-of-freedom vibration control input and output coordinate transformation matrix construction method |
CN114646440B (en) * | 2022-03-16 | 2023-06-20 | 北京卫星环境工程研究所 | Three-axis six-degree-of-freedom vibration control input/output coordinate transformation matrix construction method |
Also Published As
Publication number | Publication date |
---|---|
CN105320057B (en) | 2018-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105320057A (en) | Synchronous movement control method of double vibration tables on the basis of coordinate-transformation matrix | |
CN105196294B (en) | Reconfigurable mechanical arm decentralized control system and control method adopting position measuring | |
CN103558079B (en) | Based on the multiple degrees of freedom loading method of parallel institution driving force closed loop | |
US20180031457A1 (en) | Shear test device and test method of rock mass discontinuities under constant normal stiffness condition | |
Kang | Closed-form force sensing of a 6-axis force transducer based on the Stewart platform | |
CN104006964B (en) | Harmonic drive tooth mesh and deformation of flexible wheel dynamic characterization measurement system | |
JPS61501413A (en) | Load dynamics compensation circuit for servo hydraulic control system | |
CN101763038A (en) | Method for controlling structural modal vibration of dual-frame magnetic levitation control moment gyroscope | |
CN105015805B (en) | A kind of secondary plateform system | |
CN104407639A (en) | Servo control method and servo control device for high-frequency hydraulic shock excitation system | |
CN104807606B (en) | A kind of MATLAB STM32 hybrid powers pilot system and its test method | |
CN106052986A (en) | Three-dimensional six-degree-of-freedom force-displacement hybrid control method | |
CN108241289B (en) | Force-displacement hybrid control method for three-degree-of-freedom loading system | |
CN206348142U (en) | Double-wall pipe elastic support dynamic fatigue test device | |
CN202520846U (en) | Low-frequency vibration isolator aiming at foundation excitation | |
CN109668704B (en) | Separation type hybrid test system and test method | |
CN106202735A (en) | A kind of local nonlinearity foundation soil structural interaction sub-structural test method | |
CN105508823A (en) | Motion control method for flexible six-degree-of-freedom parallel precision adjustment device | |
CN103344506B (en) | A kind of structural face shear strength size effect Combined Trials method | |
CN109709822B (en) | Whole-satellite ground semi-physical test device for large-scale flexible spacecraft | |
CN102540903A (en) | Modelica-language-based simulation modeling method for pump truck boom system | |
CN209820722U (en) | Separated mixing test system | |
Jiang et al. | Analysis of dynamic characteristics of water hydraulic rotating angle self-servo robot joint actuator | |
Feng et al. | A novel adaptive balance-drive mechanism for industrial robots using a series elastic actuator | |
CN107830980A (en) | A kind of unidirectional shear model box of hollow variable cross-section |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180406 Termination date: 20210528 |
|
CF01 | Termination of patent right due to non-payment of annual fee |