CN105387981B - Test specimen vibration test coordinate transformation method - Google Patents

Test specimen vibration test coordinate transformation method Download PDF

Info

Publication number
CN105387981B
CN105387981B CN201510708051.4A CN201510708051A CN105387981B CN 105387981 B CN105387981 B CN 105387981B CN 201510708051 A CN201510708051 A CN 201510708051A CN 105387981 B CN105387981 B CN 105387981B
Authority
CN
China
Prior art keywords
matrix
row
control point
point
test specimen
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510708051.4A
Other languages
Chinese (zh)
Other versions
CN105387981A (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.)
Shanghai Space Precision Machinery Research Institute
Original Assignee
Shanghai Space Precision Machinery Research Institute
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 Shanghai Space Precision Machinery Research Institute filed Critical Shanghai Space Precision Machinery Research Institute
Priority to CN201510708051.4A priority Critical patent/CN105387981B/en
Publication of CN105387981A publication Critical patent/CN105387981A/en
Application granted granted Critical
Publication of CN105387981B publication Critical patent/CN105387981B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/02—Details

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

The invention provides a kind of test specimen vibration test coordinate transformation method, including (1) to determine m vibration measurement point coordinates on test specimen, as control point；(2) m control point is divided into p groups (3) to each control point weights assigned coefficient；(4) matrix Q is establishedp×m, (5) establish matrix Ym×p, (6) are by matrix Qp×mWith matrix Ym×pMultiplication obtains diagonal matrix Y﹡ p×p；(7) by matrix Y﹡ p×pDiagonal values test specimen is controlled according to the virtual point as virtual point, the shake table.The present invention makes test specimen vibratory response distribution more uniform, has been extremely improved the control effect of vibration test.

Description

Test specimen vibration test coordinate transformation method
Technical field
The present invention relates to sample vibration test technology field, in particular it relates to a kind of test specimen vibration test Coordinate Conversion side Method.
Background technology
At present, in multiaxis multiple activation vibration experiment, the control effect of control algolithm is with controlling the dimension of matrix to have Direct relation, many finite element stimulation results and experiment also directly demonstrate the quantity and the quantity of excitation point when control point Identical, that is, when to control matrix be square formation, control effect is optimal.And in specific vibration environment, shake table limited amount, to protect Card being capable of more real imitation specimen vibration environment, it is desirable to which the control points of test specimen are also more and more, conventional square formation control Experiment demand can not be met.In this case, it is therefore necessary to existing experimental control method is improved, new control method Multiple control points can be control effectively, its vibratory response spectrum shape is composed shape close proximity to given control,
The dimension of the control effect of control algolithm and control matrix has direct relation, and square formation control has with rectangular matrix control Limit first simulation result also directly to demonstrate when the quantity at control point is identical with the quantity of excitation point, that is, it is square formation to control matrix When, control effect is optimal.When the quantity at control point exceedes excitation point quantity, a part of control point of system can not obtain effectively Control, experiment are difficult to obtain preferable result of the test.In this case, vibration test security can not be protected.
The content of the invention
For in the prior art the defects of, it is an object of the invention to provide a kind of test specimen vibration test coordinate transformation method.
According to an aspect of the present invention, there is provided a kind of test specimen vibration test coordinate transformation method, it is characterized in that, including such as Lower step：
(1) m vibration measurement point coordinates is determined on test specimen, as control point；
(2) m control point is divided into p groups, p is the shake table quantity of vibration test, and every group of control point number is respectively m1,m2…mp；
(3) to each control point weights assigned coefficient；
(4) matrix Q is establishedp×m, the first first group of behavior control point weight coefficient, preceding m1 is classified as weight coefficient, and remaining is Zero；Second second group of behavior control point weight coefficient, m1+1 are arranged to m1+m2 and are classified as corresponding weight coefficient, region zero；According to this Analogize；
(5) matrix Y is establishedm×p, m1 behavior Control points before first row, remaining is zero；Secondary series m1+1 rows are to m1+ M2 behavior Control points, remaining is zero；The rest may be inferred；
(6) by matrix Qp×mWith matrix Ym×pMultiplication obtains diagonal matrix Y﹡ p×p；
(7) by matrix Y﹡ p×pDiagonal values test specimen is controlled according to the virtual point as virtual point, the shake table System.
Preferably, the m control point is grouped according to the position positioned at the p shake table.
Preferably, the m is 4, and control point is respectively y1, y2 and y3, y4, and above-mentioned 4 control points are divided into two groups, y1 with Mono- group of y2, y3 and mono- group of y4, the weight coefficient at 4 control points is 1/2, first row and the second column element in matrix Q the first rows For 1/2, remaining is zero, the second row the 3rd row and the 4th is classified as 1/2, and remaining is zero, the row of first row the first row second point in matrix Y Not Wei y1 and y2, secondary series the third line and fourth line are respectively y3 and y4, and matrix Q is multiplied by matrix Y and obtains the column matrix Y of 2 row two﹡, Element is respectively y1/2+y2/2 on its diagonal, y3/2+y4/2, and two elements on diagonal are virtual controlling point, chooses two Coherence factor is 0.5 between virtual point, and phase difference is 90 °, determines control of the shake table to test specimen.
Compared with prior art, the present invention has following beneficial effect：
(1) more vibration test control points are changed into vibrating number of units identical virtual controlling using Coordinate Conversion Point, it can still provide for effectively controlling when the measurement control point of test specimen is more than shake table quantity；
(2) make test specimen vibratory response distribution more uniform；
(3) control method is not limited to the structure type of test specimen, while can be selected according to the vibratory response situation at each control point Weight coefficient is taken, there is stronger versatility to different tests environment, be common in various multiple activation vibration tests；
(4) control method of all kinds of multiple activation vibration tests is simplified, by the conversion to control point, it is difficult to reduce experiment Degree, improve test efficiency.
Brief description of the drawings
The detailed description made by reading with reference to the following drawings to non-limiting example, further feature of the invention, Objects and advantages will become more apparent upon：
Fig. 1 is that slender rod piece dual stage joins the pilot system figure that shakes；
Fig. 2 is the vibrating effect figure that slender rod piece dual stage joins experiment of shaking.
Embodiment
With reference to specific embodiment, the present invention is described in detail.Following examples will be helpful to the technology of this area Personnel further understand the present invention, but the invention is not limited in any way.It should be pointed out that the ordinary skill to this area For personnel, without departing from the inventive concept of the premise, various modifications and improvements can be made.These belong to the present invention Protection domain.
Test specimen vibration test is according to control point position structure and vibration level, and we are by m control point on test specimen It is divided into p groups, and respectively to control point weights assigned coefficient.Wherein p is the quantity of shake table.It is related to shake table during packet The control point of connection point is in same group.
Utilize each control point weight coefficient composition coordinate conversion matrix Qp×m, m control point is divided into p groups (every group of points point Not Wei m1, m2 ...), the first first group of behavior control point weight coefficient, preceding m1 is classified as weight coefficient, and remaining is zero；Second behavior Second group of control point weight coefficient, m1+1 are arranged to m1+m2 and are classified as corresponding weight coefficient, and remaining is zero；The rest may be inferred；
Matrix Ym×pBe made up of m control point, equally by m control point be divided into p groups (every group of points be respectively m1, M2 ...), m1 behavior Control points before first row, remaining is zero；Secondary series m1+1 row value m1+m2 behavior Control points, Remaining is zero；The rest may be inferred；
By matrix Qp×mWith matrix Ym×pMultiplication obtains diagonal matrix Y﹡ p×p, it is concerned between according to Vibration Condition, simulation point is set The parameters such as coefficient, phase difference, vibration virtual point is controlled using shake table.
Being shaken below by one with slender rod piece dual stage connection, the invention will be further described for the embodiment of pilot system.
Referring to accompanying drawing 1,4 control points in experiment be present：Y1, y2 and y3, y4, can be by above-mentioned 4 according to measuring point distribution situation Individual point is divided into two groups, y1 and mono- group of mono- group of y2, y3 and y4.And it is 1/2 to make 4 dominating pair of vertices answer the Coordinate Conversion factor, it is known that In this test, coordinate conversion matrix Q is the column matrix of 2 row 4, and matrix Y is the column matrix of 4 row 2, wherein coordinate conversion matrix Q first First row and the second column element are 1/2 in row, and remaining is zero, and the second row the 3rd arranges and the 4th is classified as 1/2, and remaining is zero.Matrix Y The middle row of first row the first row second is respectively y1 and y2, and secondary series the third line and fourth line are respectively y3 and y4.Multiplied using matrix Q The column matrix Y ﹡ of 2 row two are obtained with matrix Y, it is diagonal matrix, and element is respectively y1/2+y2/2, y3/2+y4/2 on diagonal. Two elements are the target spectrum of virtual controlling point on diagonal.According to experimental condition, coherence factor is between choosing two virtual points 0.5, phase difference is 90 °.
Referring to accompanying drawing 2, after the vibrating effect of test specimen is simulated by software, the effect of oscillating curve can be clearly seen that.
If one times of the control points expansion for pilot system that dual stage connection is shaken, also will not cause shadow to system control accuracy Ring.
The specific embodiment of the present invention is described above.It is to be appreciated that the invention is not limited in above-mentioned Particular implementation, those skilled in the art can make various deformations or amendments within the scope of the claims, this not shadow Ring the substantive content of the present invention.

Claims (3)

1. a kind of test specimen vibration test coordinate transformation method, it is characterised in that comprise the following steps：
(1) m vibration measurement point coordinates is determined on test specimen, as control point；
(2) m control point being divided into p groups, p is the shake table quantity of vibration test, and every group of control point number is respectively m1, m2…mp；
(3) to each control point weights assigned coefficient；
(4) matrix Q is establishedp×m, the first first group of behavior control point weight coefficient, preceding m1 is classified as weight coefficient, and remaining is zero；The Two second group of behavior control point weight coefficients, m1+1 are arranged to m1+m2 and are classified as corresponding weight coefficient, and remaining is zero；The rest may be inferred；
(5) matrix Y is establishedm×p, first group of control point of m1 behaviors before first row, remaining is zero；Secondary series m1+1 rows are to m1+m2 Second group of control point of behavior, remaining is zero；The rest may be inferred；
(6) by matrix Qp×mWith matrix Ym×pMultiplication obtains diagonal matrix Y﹡ p×p；
(7) by matrix Y﹡ p×pDiagonal values test specimen is controlled according to the virtual point as virtual point, the shake table.
2. test specimen vibration test coordinate transformation method according to claim 1, it is characterised in that m control point root It is grouped according to positioned at the position of the p shake table.
3. test specimen vibration test coordinate transformation method according to claim 1, it is characterised in that the m is 4, control point Respectively y1, y2 and y3, y4, above-mentioned 4 control points are divided into two groups, y1 and mono- group of mono- group of y2, y3 and y4,4 control points Weight coefficient is 1/2, and first row and the second column element are 1/2 in matrix Q the first rows, and remaining is zero, the second row the 3rd row and 4th is classified as 1/2, and remaining is zero, and the row of first row the first row second is respectively y1 and y2 in matrix Y, secondary series the third line and the 4th Row is respectively y3 and y4, and matrix Q is multiplied by matrix Y and obtains the column matrix Y of two row two﹡, element is respectively y1/2+y2/ on its diagonal 2, y3/2+y4/2, two elements on diagonal are virtual controlling point, and coherence factor is 0.5 between choosing two virtual points, phase difference For 90 °, control of the shake table to test specimen is determined.
CN201510708051.4A 2015-10-27 2015-10-27 Test specimen vibration test coordinate transformation method Active CN105387981B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510708051.4A CN105387981B (en) 2015-10-27 2015-10-27 Test specimen vibration test coordinate transformation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510708051.4A CN105387981B (en) 2015-10-27 2015-10-27 Test specimen vibration test coordinate transformation method

Publications (2)

Publication Number Publication Date
CN105387981A CN105387981A (en) 2016-03-09
CN105387981B true CN105387981B (en) 2017-12-01

Family

ID=55420523

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510708051.4A Active CN105387981B (en) 2015-10-27 2015-10-27 Test specimen vibration test coordinate transformation method

Country Status (1)

Country Link
CN (1) CN105387981B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109765106A (en) * 2019-01-08 2019-05-17 上海航天精密机械研究所 Three-dimensional loading system and method applied to cylindrical product
CN110068436B (en) * 2019-04-04 2021-03-30 航天科工防御技术研究试验中心 Vibration control implementation method and device for thermal vibration detection
CN114216637B (en) * 2021-12-24 2023-09-15 苏州东菱振动试验仪器有限公司 Multi-synchronous vibration test device and test method thereof
CN114646440B (en) * 2022-03-16 2023-06-20 北京卫星环境工程研究所 Three-axis six-degree-of-freedom vibration control input/output coordinate transformation matrix construction method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001133357A (en) * 1999-11-02 2001-05-18 Imv Kk Vibration test system
CN102122151A (en) * 2010-12-10 2011-07-13 南京航空航天大学 Control device and control method used for multi-dimensional random vibration test
CN202350998U (en) * 2011-10-28 2012-07-25 北京航空航天大学 Multi-excitation-source combined control test device

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3074358B2 (en) * 1991-12-05 2000-08-07 株式会社日立製作所 Vibration test apparatus, vibration test method and vibration response analysis method for structures

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001133357A (en) * 1999-11-02 2001-05-18 Imv Kk Vibration test system
CN102122151A (en) * 2010-12-10 2011-07-13 南京航空航天大学 Control device and control method used for multi-dimensional random vibration test
CN202350998U (en) * 2011-10-28 2012-07-25 北京航空航天大学 Multi-excitation-source combined control test device

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
"Applying Coordinate Transformations to Multi-DOF Shaker Control";Marcos A.Underwood 等;《SOUND AND VIBRATION》;20060131;第22-27页 *
"Matrix Power Control Algorithm for Multi-input Multi-output Random Vibration Test";CUI Xuli 等;《Chinese Journal of Aeronautics》;20111231(第24（2011）期);第741-748页 *
"Recent System Developments for Multi-Actuator Vibration Control";Marcos A.Underwood 等;《SOUND AND VIBRATION》;20011031;第1-6页 *
"多台振动试验控制方法研究";王新;《强度与环境》;20080430;第35卷(第2期);第7-11页 *
"多输入多输出振动试验应用综述";赵保平 等;《装备环境工程》;20060630;第3卷(第3期);第25-32页 *
"某记载雷达多自由度振动试验控制方法比较";司中柱 等;《环境试验》;20150228;第5-8页 *

Also Published As

Publication number Publication date
CN105387981A (en) 2016-03-09

Similar Documents

Publication Publication Date Title
CN105387981B (en) Test specimen vibration test coordinate transformation method
CN110083920B (en) Analysis method for random response of non-proportional damping structure under earthquake action
CN103778292A (en) Method for predicting fatigue life of BGA (Ball Grid Array) welding spot under heat-vibration combined loads
US20100145668A1 (en) Method for dynamic repartitioning in adaptive mesh processing
CN103983412A (en) Avionic device operating modal measuring method for vibration finite element model correction
CN105571866A (en) Engine vibration spectrum extraction method
Mitchell et al. A method for enforcement of Dirichlet boundary conditions in isogeometric analysis
Fuellekrug et al. Identification of weak non-linearities within complex aerospace structures
Mandal et al. Pressure based Eulerian approach for investigation of sloshing in rectangular water tank
CN109783935A (en) A kind of implementation method improving splash fluid stability based on ISPH
JP7135976B2 (en) CONTROLLER SYSTEM, SUPPORT DEVICE AND EVALUATION METHOD
CN107727350B (en) Micro-nano satellite vector vibration test method
ITFI970210A1 (en) METHOD FOR THE REPRESENTATION USING A SEMI-STRUCTURED CALCULATION GRID OF PHYSICAL PHENOMENA EXTENDING IN A SPATIAL DOMAIN
CN109117519B (en) Method for simulating transient temperature field of oil tank structure
CN103440391A (en) Semiconductor process corner scanning and simulating method based on numerical value selection function
Tondini et al. Hybrid simulation applied to fire testing: a newly conceived numerical framework
US10915678B1 (en) Surface mapping apparatuses and methods
Casini et al. Effect of damping on the nonlinear modal characteristics of a piecewice-smooth system through harmonic forced response
CA2028807C (en) Decentralized simulator
CN108510575B (en) Method for determining shielding relation of large object on 2D game oblique 45-degree map
Szilagyi et al. A systematic analysis of the N‐person chicken game
Berman Inherently incomplete finite element model and its effects on model updating
CN114077763A (en) Nuclear power plant containment vessel structure determining method and device
Uehara et al. Numerical simulation of homogeneous polycrystalline grain formation using multi-phase-field model
Bhothikhun et al. Adaptive DKT finite element for plate bending analysis of built-up structures

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