CN102004822A - Modality-based frequency analysis method for space six-freedom-degree parallel motion system - Google Patents
Modality-based frequency analysis method for space six-freedom-degree parallel motion system Download PDFInfo
- Publication number
- CN102004822A CN102004822A CN 201010537154 CN201010537154A CN102004822A CN 102004822 A CN102004822 A CN 102004822A CN 201010537154 CN201010537154 CN 201010537154 CN 201010537154 A CN201010537154 A CN 201010537154A CN 102004822 A CN102004822 A CN 102004822A
- Authority
- CN
- China
- Prior art keywords
- space
- frequency
- matrix
- freedom
- parallel kinematic
- 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
Images
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The invention provides a modality-based frequency analysis method for a space six-freedom-degree parallel motion system. The invention is realized by comprising the following steps of: creating a modality equation of a parallel motion system according to a second-order vibration differential equation of the space docking six-freedom-degree parallel motion system; constructing an equation system matrix; decomposing the system matrix by utilizing a numerical analysis method; calculating the low-order modality set of the six-freedom-degree parallel motion system and a corresponding characteristic frequency; and analyzing the low-order modality frequency of the space six-freedom-degree parallel motion system. The invention breaks through the traditional method based on physical space analysis and design, can be applied to the design and the analysis of all space six-freedom-degree parallel motion systems to ensure that the designed mechanism satisfies the requirement on a designed bandwidth, provides effective basis for changing the system modality, provides a new method and concept for more effective six-freedom-degree structural design and optimization and has important significance for practical application and theoretical analysis of a project.
Description
(1) technical field
The present invention relates to mechanical technique, is a kind of space six-freedom parallel kinematic system frequency analysis method based on mode specifically.
(2) background technology
Space articulation mechanism six-freedom parallel kinematic system adopts hydraulically powered typical Stewart platform form, be with the different of parallel kinematic system of present existence, space articulation parallel kinematic system requirements kinematic system has super large space ability, overlength supporting leg impulse stroke and high system's frequency range.At present, the Stewart platform that satisfies the index request of space articulation kinematic system does not still have precedent at home, and also only the U.S. and Russia have the type docking mechanism kinematic system in the world.The parallel kinematic system works frequency range of this large-scale over-long supporting leg form can reach 10Hz, and current hydraulic-driven parallel kinematic simulation system frequency range 3-5Hz only, and work space is little.Require and the big parallel kinematic system of work space for high frequency range, ignore its structure and can't satisfy the actual techniques index with the system that driving mode will cause designing.Design the parallel kinematic structure that satisfies space articulation mechanism kinematic performance requirement, must carry out frequency-response analysis system.
(3) summary of the invention
The object of the present invention is to provide a kind of space six-freedom parallel kinematic system frequency analysis method based on mode.
The object of the present invention is achieved like this: according to the second order oscillatory differential equation of space articulation six-freedom parallel kinematic system, set up the modal equations of parallel kinematic system, the equationof structure system matrix, adopt numerical analysis method that system matrix is decomposed, calculate the lower mode collection and the characteristic of correspondence frequency thereof of six-freedom parallel kinematic system, analyze the lower mode frequency of space six-freedom parallel kinematic system; Space six-freedom parallel kinematic system works near the meta usually, and its lower mode frequency analysis mainly concentrates on meta on engineering is used, and this moment the second order model of vibration of parallel kinematic system is represented with following form:
M is the wide mass matrix of n * n in the formula; C is n * n broad sense damping matrix; K is n * n broad sense Stiffness Matrix; F (t) is the outside input stimulus of system;
Be analytic system meta model frequency, damping term and external drive item are removed, so, obtain undamped second order vibration equation, suppose that its model frequency is ω
n, Mode Shape is ψ, then obtains space six-freedom parallel kinematic system modal equations,
(K-ω
nM)ψ=0
The structure matrix of coefficients, matrix of coefficients equals the matrix product of the contrary and n dimension broad sense Stiffness Matrix of n * n dimension generalized mass battle array, adopt methods of numerical then, characteristic value decomposition, svd calculate mode and its characteristic of correspondence frequency of the needs investigation of space six-freedom parallel kinematic system, try to achieve the corresponding model frequency of system according to characteristic frequency, reflect according to the model frequency of analyzing whether designed system meets the demands.
A kind of space six-freedom parallel kinematic system frequency analysis method based on mode of the present invention is broken through the method for tradition based on physical space analysis and design, has especially overcome the universal method based on the design of single channel equivalent frequency.The present invention can be applied in the design and analysis of the six-freedom parallel system that has living space, make the mechanism of design satisfy the requirement of design frequency range, for a change system mode provides effective foundation, for more effective six degree of freedom structural design provides new method and thinking with optimizing, for engineering practical application and theoretical analysis, all significant.
(4) description of drawings
Fig. 1 is a space articulation mechanism six-freedom parallel kinematic system structural representation.
(5) embodiment
The invention will be further described for example below in conjunction with accompanying drawing.
Embodiment 1: in conjunction with Fig. 1, a kind of space six-freedom parallel kinematic system frequency analysis method of the present invention based on mode, second order oscillatory differential equation according to space articulation six-freedom parallel kinematic system, set up the modal equations of parallel kinematic system, the equationof structure system matrix, adopt numerical analysis method that system matrix is decomposed, calculate the lower mode collection and the characteristic of correspondence frequency thereof of six-freedom parallel kinematic system, analyze the lower mode frequency of space six-freedom parallel kinematic system; Space six-freedom parallel kinematic system works near the meta usually, and its lower mode frequency analysis mainly concentrates on meta on engineering is used, and this moment the second order model of vibration of parallel kinematic system is represented with following form:
M is the wide mass matrix of n * n in the formula; C is n * n broad sense damping matrix; K is n * n broad sense Stiffness Matrix; F (t) is the outside input stimulus of system;
Be analytic system meta model frequency, damping term and external drive item are removed, so, obtain undamped second order vibration equation, suppose that its model frequency is ω
n, Mode Shape is ψ, then obtains space six-freedom parallel kinematic system modal equations,
(K-ω
nM)ψ=0
The structure matrix of coefficients, matrix of coefficients equals the matrix product of the contrary and n dimension broad sense Stiffness Matrix of n * n dimension generalized mass battle array, adopt methods of numerical then, characteristic value decomposition, svd calculate mode and its characteristic of correspondence frequency of the needs investigation of space six-freedom parallel kinematic system, try to achieve the corresponding model frequency of system according to characteristic frequency, reflect according to the model frequency of analyzing whether designed system meets the demands.
Embodiment 2: a kind of space six-freedom parallel kinematic system frequency analysis method of the present invention based on mode, parallel kinematic system of described space articulation mechanism by motion platform (1), go up and connect hinge (2), hydraulic-driven supporting leg (3), connect hinge (4) and fixed base platform (5) down and form.It is 4.0m that meta props up leg length, stroke 1.850m.Last hinge garden diameter 3.0m, following hinge garden diameter 7.3m.Platform and load gross mass are 563Kg.At the space six-freedom parallel kinematic system under this structure, when keeping its meta in the second order oscillatory differential equation 6 * 6 of the broad sense Stiffness Matrix dimension Stiffness Matrix and 6 * 6 the dimension mass matrixes, then generalized mass matrix and Stiffness Matrix are:
Damping term in the equation and external drive item are removed, construct its matrix of coefficients by space six-freedom parallel kinematic system ground modal equations, this 6 * 6 maintains matrix number and is:
Adopt the method for numerical value, the factor arrays that is obtained by modal equations is carried out characteristic value decomposition, obtained the modal matrix and the characteristic frequency vector of this parallel institution.
The model frequency of the space articulation mechanism six-freedom parallel kinematic system that can try to achieve by the characteristic frequency that calculates, this model frequency is expressed as:
Six-freedom parallel system of space articulation mechanism to this structure carries out model frequency identification experiment, obtain its corresponding to the minimum model frequency of hydraulic pressure rigidity, this frequency is about 22Hz, and the first rank model frequency that identification obtains is close with the second rank model frequency, its mode mainly is mapped in X and Y translation direction, the 3rd rank mode is mapped in the heave direction, and fourth, fifth rank mode mainly is mapped in inclination and pitch orientation, and the 6th rank mode mainly is mapped in yaw direction.The consistance of test findings and analysis result of the present invention has well illustrated feasibility and the space six degree of freedom model frequency that proposes of the present invention validity of analyzing thought in engineering reality of the present invention.
Claims (1)
1. space six-freedom parallel kinematic system frequency analysis method based on mode, it is characterized in that: according to the second order oscillatory differential equation of space articulation six-freedom parallel kinematic system, set up the modal equations of parallel kinematic system, the equationof structure system matrix, adopt numerical analysis method that system matrix is decomposed, calculate the lower mode collection and the characteristic of correspondence frequency thereof of six-freedom parallel kinematic system, analyze the lower mode frequency of space six-freedom parallel kinematic system; Space six-freedom parallel kinematic system works near the meta usually, and its lower mode frequency analysis mainly concentrates on meta on engineering is used, and this moment the second order model of vibration of parallel kinematic system is represented with following form:
M is the wide mass matrix of n * n in the formula; C is n * n broad sense damping matrix; K is n * n broad sense Stiffness Matrix; F (t) is the outside input stimulus of system;
Be analytic system meta model frequency, damping term and external drive item are removed, so, obtain undamped second order vibration equation, suppose that its model frequency is ω
n, Mode Shape is ψ, then obtains space six-freedom parallel kinematic system modal equations,
(K-ω
nM)ψ=0
The structure matrix of coefficients, matrix of coefficients equals the matrix product of the contrary and n dimension broad sense Stiffness Matrix of n * n dimension generalized mass battle array, adopt methods of numerical then, characteristic value decomposition, svd calculate mode and its characteristic of correspondence frequency of the needs investigation of space six-freedom parallel kinematic system, try to achieve the corresponding model frequency of system according to characteristic frequency, reflect according to the model frequency of analyzing whether designed system meets the demands.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010105371546A CN102004822B (en) | 2010-11-10 | 2010-11-10 | Modality-based frequency analysis method for space six-freedom-degree parallel motion system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010105371546A CN102004822B (en) | 2010-11-10 | 2010-11-10 | Modality-based frequency analysis method for space six-freedom-degree parallel motion system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102004822A true CN102004822A (en) | 2011-04-06 |
CN102004822B CN102004822B (en) | 2012-10-31 |
Family
ID=43812182
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010105371546A Expired - Fee Related CN102004822B (en) | 2010-11-10 | 2010-11-10 | Modality-based frequency analysis method for space six-freedom-degree parallel motion system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102004822B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102622936A (en) * | 2012-01-29 | 2012-08-01 | 上海交通大学 | Parallel experimental device for six-degree-of-freedom gait simulation |
CN103558079A (en) * | 2013-10-21 | 2014-02-05 | 哈尔滨工业大学 | Multi-degree-of-freedom loading method based on parallel mechanism driving force closed loop |
CN105030481A (en) * | 2015-07-10 | 2015-11-11 | 哈尔滨天愈康复医疗机器人有限公司 | Redundant drive three-degree-of-freedom ankle rehabilitation medical robot |
CN106076825A (en) * | 2016-06-16 | 2016-11-09 | 安徽旺达铜业发展有限公司 | A kind of copper powder screening intelligent machine |
CN109033025A (en) * | 2018-07-18 | 2018-12-18 | 中国海洋大学 | Floating structure time domain response analysis method based on state-space model |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030111420A1 (en) * | 2001-12-14 | 2003-06-19 | Dresty John Edward | Process of forming an oxidizing agent in liquid by use of ringing magnetic flux |
WO2007102155A2 (en) * | 2006-03-08 | 2007-09-13 | Expert Dynamics Ltd. | Pcb design reliability simulation method and system |
CN101051328A (en) * | 2007-01-26 | 2007-10-10 | 西北工业大学 | Reusable parameter module model building method for space contineous deformation MEMS |
CN101143443A (en) * | 2007-09-20 | 2008-03-19 | 哈尔滨工程大学 | Six-freedom parallel robot decoupling method |
-
2010
- 2010-11-10 CN CN2010105371546A patent/CN102004822B/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030111420A1 (en) * | 2001-12-14 | 2003-06-19 | Dresty John Edward | Process of forming an oxidizing agent in liquid by use of ringing magnetic flux |
WO2007102155A2 (en) * | 2006-03-08 | 2007-09-13 | Expert Dynamics Ltd. | Pcb design reliability simulation method and system |
CN101051328A (en) * | 2007-01-26 | 2007-10-10 | 西北工业大学 | Reusable parameter module model building method for space contineous deformation MEMS |
CN101143443A (en) * | 2007-09-20 | 2008-03-19 | 哈尔滨工程大学 | Six-freedom parallel robot decoupling method |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102622936A (en) * | 2012-01-29 | 2012-08-01 | 上海交通大学 | Parallel experimental device for six-degree-of-freedom gait simulation |
CN102622936B (en) * | 2012-01-29 | 2014-08-20 | 上海交通大学 | Parallel experimental device for six-degree-of-freedom gait simulation |
CN103558079A (en) * | 2013-10-21 | 2014-02-05 | 哈尔滨工业大学 | Multi-degree-of-freedom loading method based on parallel mechanism driving force closed loop |
CN103558079B (en) * | 2013-10-21 | 2015-11-25 | 哈尔滨工业大学 | Based on the multiple degrees of freedom loading method of parallel institution driving force closed loop |
CN105030481A (en) * | 2015-07-10 | 2015-11-11 | 哈尔滨天愈康复医疗机器人有限公司 | Redundant drive three-degree-of-freedom ankle rehabilitation medical robot |
CN106076825A (en) * | 2016-06-16 | 2016-11-09 | 安徽旺达铜业发展有限公司 | A kind of copper powder screening intelligent machine |
CN109033025A (en) * | 2018-07-18 | 2018-12-18 | 中国海洋大学 | Floating structure time domain response analysis method based on state-space model |
Also Published As
Publication number | Publication date |
---|---|
CN102004822B (en) | 2012-10-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102004822B (en) | Modality-based frequency analysis method for space six-freedom-degree parallel motion system | |
Liu et al. | An efficient method to reduce ill-posedness for structural dynamic load identification | |
CN102567567B (en) | Finite element analysis based pipeline random-vibration fatigue life analyzing method | |
CN110298105B (en) | CCPDI-IMPM method for large deformation analysis of saturated porous medium | |
CN101982822B (en) | Modal modeling method of kinematic system with spatial six degrees of freedom | |
Luo et al. | Model order reduction for dynamic simulation of a flexible multibody system via absolute nodal coordinate formulation | |
CN105069241A (en) | Step-by-step analysis and prediction method for dynamic performances of rubber material structure | |
CN103353911B (en) | Modal method introduces initial imperfection subsea pipeline complete buckling two dimensional dynamic analogy method | |
CN111125963B (en) | Numerical simulation system and method based on Lagrange integral point finite element | |
CN102787677B (en) | Engineering structure seismic intensity prediction method | |
CN104200052A (en) | System and method for simulating hydraulic-drive hexapod robot | |
CN104281730A (en) | Great-rotating-deformation plate shell structure dynamic response finite element analysis method | |
CN104268322A (en) | Boundary processing technology of WENO difference method | |
CN110261141B (en) | Suspension tunnel vehicle-tunnel power coupling hybrid simulation test method and device | |
Dindorf et al. | Review of Hydro-Pneumatic Accumulator Models for the Study of the Energy Efficiency of Hydraulic Systems | |
CN104807606A (en) | MATLAB-STM32 hybrid power test system and test method thereof | |
CN103714206A (en) | Subsea pipeline two-dimensional static numerical simulation method for introducing initial defects based on modes | |
CN106202735A (en) | A kind of local nonlinearity foundation soil structural interaction sub-structural test method | |
CN117316358A (en) | Multi-scale analysis method for composite material structure and application thereof | |
CN102768493B (en) | Intelligent vibration control device of ocean platform | |
CN116611295A (en) | Offshore wind turbine blade dynamic response calculation method and system | |
CN111859738B (en) | Demonstration method for disaster-causing mechanism of multiple physical fields of pipeline in storm environment | |
CN102567561B (en) | Modeling method for discretizing plate shell structure | |
CN101853323B (en) | Modeling method for acting force of full-automatic pulsator washing machine suspension system | |
CN103729504A (en) | Method for simulating three-dimensional static numerical values of submarine pipelines with introduced initial defects on basis of modes |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121031 Termination date: 20131110 |