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 PDF

Info

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
Application number
CN 201010537154
Other languages
Chinese (zh)
Other versions
CN102004822B (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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN2010105371546A priority Critical patent/CN102004822B/en
Publication of CN102004822A publication Critical patent/CN102004822A/en
Application granted granted Critical
Publication of CN102004822B publication Critical patent/CN102004822B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

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

Space six-freedom parallel kinematic system frequency analysis method based on mode
(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 x · · ( t ) + C x · ( t ) + Kx ( t ) = f ( t )
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 x · · ( t ) + C x · ( t ) + Kx ( t ) = f ( t )
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:
M = 1839 0 0 0 428.9 0 0 1839 0 - 428.9 0 0 0 0 1638 0 0 0 0 - 428.9 0 1470 0 0 428 . 9 0 0 0 1470 0 0 0 0 0 0 2708
K = 4.341 × 10 7 0 0 0 - 7.085 × 10 6 0 0 4.341 × 10 7 0 7.085 × 10 6 0 0 0 0 9.286 × 10 7 0 0 0 0 7.085 × 10 6 0 8.379 × 10 7 0 0 - 7.085 × 10 6 0 0 0 8.379 × 10 7 0 0 0 0 0 0 1.545 × 10 8
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:
D = 2.653 × 10 4 0 0 0 - 1.839 × 10 4 0 0 2.653 × 10 4 0 1.839 × 10 4 0 0 0 0 5.668 × 10 4 0 0 0 0 1.256 × 10 4 0 6.235 × 10 4 0 0 - 1.256 × 10 4 0 0 0 6.235 × 10 4 0 0 0 0 0 0 5.707 × 10 4
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.
U = 0.3793 0 0 0 0 - 0.9253 0 - 0.3793 0 0 - 0.9253 0 0 0 0 - 1 0 0 0 - 0.9253 0 0 0.3793 0 - 0.9253 0 0 0 0 - 0.3793 0 0 - 1 0 0 0
E freq = 682808.8 682808.8 57070.11 56678.48 20866.53 20866.53
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:
ω n = 41.6 41.6 38.0 37.9 23 23
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 x · · ( t ) + C x · ( t ) + Kx ( t ) = f ( t )
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.
CN2010105371546A 2010-11-10 2010-11-10 Modality-based frequency analysis method for space six-freedom-degree parallel motion system Expired - Fee Related CN102004822B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (4)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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