CN108829924A - The calculation method of the mechanics parameter of connection structure - Google Patents
The calculation method of the mechanics parameter of connection structure Download PDFInfo
- Publication number
- CN108829924A CN108829924A CN201810421619.8A CN201810421619A CN108829924A CN 108829924 A CN108829924 A CN 108829924A CN 201810421619 A CN201810421619 A CN 201810421619A CN 108829924 A CN108829924 A CN 108829924A
- Authority
- CN
- China
- Prior art keywords
- finite element
- model
- singular value
- element model
- vibration signal
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The present invention provides a kind of calculation method of the mechanics parameter of connection structure, including:According to the physical model of hanging equipment and connection structure, finite element model relevant to mechanics parameter is constructed;Identical exciting force is applied to the physical model and finite element model respectively, obtains the singular value vector of the singular value vector of the vibration signal of the physical model and the vibration signal of the finite element model;Using the error between physical model and the corresponding singular value vector of finite element model as objective function, mechanics parameter when making the objective function stable convergence is determined using optimization algorithm, as the hanging equipment and the actual mechanics parameter of connection structure.The present invention is using the vibration signal for surveying physical model, and calculate the singular value vector of vibration signal, strong technical support is provided to the identification of mechanics parameter, and since the vibration signal of physical model is derived from the test of actual physical structure, the mechanics parameter finally identified also more has convincingness.
Description
Technical field
The present invention relates to Structural Dynamics technical fields, more particularly, to the calculating side of the mechanics parameter of connection structure
Method.
Background technique
In field of track traffic, motor-car, which is hung oneself, to be hung for mostly using connection structure to be hung below car body.Carrying out power
It is most important to the accurate vibration equivalence modeling of the mechanics parameter progress of connection structure in the model of credit analysis, be also simultaneously
Urgent technical problem to be solved in the field.
Academic research is calculated around the analysis work based on finite element model, such as modal calculation, noise calculation, dynamic response
Achievement it is commonly reported that.In addition, equipment and its connection structure will be hung separately as research object, its finite element model is established, is led to
It crosses inquiry pertinent literature, empirical equation or dynamic test and obtains the parameter of connection structure, then study mechanical characteristic also often to have
Report.
Existing hanging equipment and its connection structure are in modeling process, if obtaining connection structure by inquiry document
Mechanics parameter causes the finite element model accuracy established limited due to batches of materials, working condition, the otherness of manufacturer;
And if obtaining the mechanics parameter of connection structure by quasi-static mechanical test, it can not be well reflected using the model of foundation
The kinetic characteristics of practical hanging connection structure.
Summary of the invention
The present invention provides a kind of mechanics of connection structure for overcoming the above problem or at least being partially solved the above problem
The calculation method of parameter.
According to an aspect of the present invention, a kind of calculation method of the mechanics parameter of connection structure is provided, including:
According to the physical model of hanging equipment and connection structure, finite element model relevant to mechanics parameter is constructed;
Identical exciting force is applied to the physical model and finite element model respectively, obtains the vibration of the physical model
The singular value vector of the vibration signal of the singular value vector of signal and the finite element model;
Using the error between physical model and the corresponding singular value vector of finite element model as objective function, calculated using optimization
Method determines mechanics parameter when making the objective function stable convergence, as the hanging equipment and the actual mechanics of connection structure
Parameter.
Preferably, described that identical exciting force is applied respectively to the physical model and finite element model, obtain the reality
The singular value vector of the vibration signal of the singular value vector and finite element model of the vibration signal of body Model, specially:
Exciting force is applied to the physical model, while acquiring pumping signal and vibration signal, the physical model is produced
Raw vibration signal carries out singular value decomposition, obtains the singular value vector of the vibration signal of physical model;
The pumping signal is loaded onto finite element model, the vibration signal of the finite element model is calculated, to described
The vibration signal of finite element model carries out singular value decomposition, obtains the singular value vector of the vibration signal of finite element model.
Preferably, the error using between physical model and the corresponding singular value vector of finite element model is as target letter
Number, specially:
Using following formula as the objective function:
Wherein, saIndicate the corresponding singular value vector of finite element model, seThe corresponding singular value vector of presentation-entity model,
Error between Δ s presentation-entity model and the corresponding singular value vector of finite element model.
Preferably, described that finite element relevant to mechanics parameter is constructed according to the physical model of hanging equipment and connection structure
Model, specially:
Grid dividing is carried out to the hanging equipment using the lumped mass dot element in finite element analysis software, using having
Connection unit in finite element analysis software carries out grid dividing to the connection structure, establishes the finite element model;
The quality of the physical model and rotary inertia are input in the finite element model, as finite element model
Quality and rotary inertia;The value range of mechanics parameter in finite element model is set.
Preferably, described that exciting force is applied to physical model, while pumping signal and vibration signal are acquired, specially:
Setting at least one measured point passes through acceleration sensing using equipment is hung described in power hammer excitation in hanging equipment
Device acquires the acceleration responsive at least one measured point in 20Hz, as vibration signal;
Wherein, the measurement direction of acceleration transducer is translation mode direction, and test frequency range is estimate modal frequency 2
Times, signal sampling resolution ratio is 2.56 times for estimating the mode frequency range upper limit.
Preferably, the optimization algorithm is one in ant group algorithm, genetic algorithm, simulated annealing and particle swarm algorithm
Kind.
Preferably, the mechanics parameter is stiffness parameters and damping parameter.
According to another aspect of the present invention, a kind of computing system of the mechanics parameter of connection structure is also provided, including:
Finite element model generation module, for the physical model according to hanging equipment and connection structure, building is joined with mechanics
The relevant finite element model of number;
Singular value vector obtains module, for applying identical excitation respectively to the physical model and finite element model
Power, obtain the singular value of the singular value vector of the vibration signal of the physical model and the vibration signal of the finite element model to
Amount;
Optimization module, for using the error between physical model and the corresponding singular value vector of finite element model as target letter
Number determines mechanics parameter when making the objective function stable convergence using optimization algorithm, as the hanging equipment and connection
The actual mechanics parameter of structure.
According to another aspect of the present invention, a kind of electronic equipment of the mechanics parameter of connection structure is also provided, including:
At least one processor;And
At least one processor being connect with the processor communication, wherein:
The memory is stored with the program instruction that can be executed by the processor, and the processor calls described program to refer to
Enable the method for being able to carry out the calculation method and its any alternative embodiment of the mechanics parameter of connection structure of the embodiment of the present invention.
According to another aspect of the present invention, a kind of non-transient computer readable storage medium is also provided, it is described non-transient
Computer-readable recording medium storage computer instruction, the computer instruction make the computer execute the embodiment of the present invention
The method of the calculation method of the mechanics parameter of connection structure and its any alternative embodiment.
The calculation method of the mechanics parameter of connection structure proposed by the present invention, it is consistent limited with physical model by constructing
Meta-model is adopted using the error between physical model and the singular value vector of the corresponding vibration signal of finite model as objective function
With the vibration signal of actual measurement physical model, and the singular value vector of vibration signal is calculated, with the power to finite element dynamics
It learns parameter to be identified, provides strong technical support to the identification of mechanics parameter, and since the vibration of physical model is believed
Number be derived from actual physical structure test, the mechanics parameter finally identified also more have convincingness.It is determined using optimization algorithm
Make the mechanics parameter when objective function stable convergence, as the hanging equipment and the actual mechanics parameter of connection structure,
Since the singular value vector of vibration signal is divided into sensitivity to the variation of mechanics parameter, using the error between singular value vector as target
Function, so that objective function can also occur significantly to change even if small variation has occurred in mechanics parameter.It is calculated by optimization
The value of mechanics parameter in method optimized FEMs model, so that the corresponding singular value vector of finite element model and entity mould after optimization
The corresponding singular value vector of type is very close, and the finite element model after achieving the effect that optimization being capable of actual response physical model.
Detailed description of the invention
Fig. 1 is the flow diagram according to the calculation method of the mechanics parameter of the connection structure of the embodiment of the present invention;
Fig. 2 is the mechanics determined when making the objective function stable convergence using optimization algorithm according to the embodiment of the present invention
The flow diagram of the method for parameter;
Fig. 3 is the functional block diagram according to the computing system of the mechanics parameter of the connection structure of the embodiment of the present invention;
Fig. 4 is to make the objective function stable convergence using optimization algorithm determination according to the computing system of the embodiment of the present invention
When mechanics parameter method flow diagram;
Fig. 5 is the block schematic illustration according to the electronic equipment of the embodiment of the present invention;
Fig. 6 is the convergence curve figure of the mechanics parameter of a specific embodiment of the invention;
Fig. 7 is the convergence curve figure of the objective function of a specific embodiment of the invention.
Specific embodiment
With reference to the accompanying drawings and examples, specific embodiments of the present invention will be described in further detail.Implement below
Example is not intended to limit the scope of the invention for illustrating the present invention.
In order to overcome the drawbacks described above of the prior art, the embodiment of the present invention provides a kind of meter of the mechanics parameter of connection structure
Calculation method, Fig. 1 are the flow diagram of the calculation method of the mechanics parameter of the connection structure of the embodiment of the present invention, as shown in Figure 1,
This method includes:
101, according to the physical model of hanging equipment and connection structure, finite element model relevant to mechanics parameter is constructed.
It should be noted that the physical model of hanging equipment and connection structure hangs equipment and connection structure itself.It hangs
The standby equipment mounted by connection structure is hung, for example crane hangs shot by hook, shot is to hang equipment, hook
Son is bindiny mechanism.The embodiment of the present invention is not intended to limit the type and calculation method of specific hanging equipment and connection structure
Applied field.Finite element model is made of node and element.The finite element model of the embodiment of the present invention is according to hanging equipment
And the complexity of connection structure, it is constructed using two methods:If hanging equipment and connection structure being relatively simple, direct basis
The geometry external form of mechanical structure establishes node and unit;If hanging equipment and connection structure it is complex, by point, line, surface,
Volume first establishes finite element model, then carries out physical grid division, to complete the foundation of finite element model.
102, identical exciting force is applied to the physical model and finite element model respectively, obtains the entity mould respectively
The singular value vector of the vibration signal of the singular value vector and finite element model of the vibration signal of type.
It should be noted that the acceleration responsive of physical model can be collected by applying exciting force to physical model,
That is vibration signal can get the singular value vector of vibration signal by further carrying out singular value decomposition to vibration signal.Institute
It is because singular value can effectively inhibit vibration signal also to need further to obtain singular value vector after obtaining vibration signal
In interference and background random noise, while enhancing the energy value of vibration signal, and by the verifying of inventor, vibration signal
Singular value vector it is very sensitive for the variation of mechanics parameter, i.e., can be more accurately using the singular value vector of vibration signal
Demarcate the mechanics parameter in finite element model.Using the vibration signal of actual measurement physical model, and calculate the singular value of vibration signal
Vector is identified with the mechanics parameter to finite element dynamics, provides strong technology to the identification of mechanics parameter
Support, and since the vibration signal of physical model is derived from the test of actual physical structure, the mechanics parameter finally identified is also more
Add with convincingness.It is understood that identical exciting force refers to the size in same position excitation in the embodiment of the present invention
Exciting force identical with direction.
103, using the error between physical model and the corresponding singular value vector of finite element model as objective function, utilization is excellent
Change algorithm determines mechanics parameter when making the objective function stable convergence, actual as the hanging equipment and connection structure
Mechanics parameter.
It should be noted that the singular value vector due to vibration signal is divided into sensitivity to the variation of mechanics parameter, it will be unusual
It is worth the error between vector as objective function, so that objective function can also be sent out even if small variation has occurred in mechanics parameter
Raw apparent variation.By the value of mechanics parameter in optimization algorithm optimized FEMs model, so that the finite element model after optimization
Corresponding singular value vector singular value vector corresponding with physical model is very close, and the finite element model after reaching optimization can
The effect of actual response physical model.It should be noted that found when by the calculation method practical application of the embodiment of the present invention,
It only needs Optimized Iterative 3-5 times, objective function stable convergence can be realized, the statics test of the prior art is substantially better than in efficiency
Test method or dynamic testing machine test method.
Mechanics ginseng when making the objective function stable convergence is determined using optimization algorithm Fig. 2 shows the embodiment of the present invention
The flow diagram of several methods, as shown, this method includes:
Step 201 applies exciting force to physical model, obtains the singular value vector of the vibration signal of physical model, executes
Step 202;
Step 202 generates new mechanics parameter using optimization algorithm, executes step 203;
Step 203 updates finite element model using new mechanics parameter, and exciting force is input to updated finite element mould
In type, the singular value vector of new vibration signal is obtained, executes step 204;
Step 204, by the singular value vector generation of the singular value vector of new vibration signal and the vibration signal of physical model
Enter objective function, executes step 205;
Step 205, judge objective function whether stable convergence, if so, 206 are thened follow the steps, if it is not, then return step
202;
Step 206 terminates optimization.
On the basis of the above embodiments, the same position to the physical model and finite element model applies excitation
Power, to obtain the singular value of the singular value vector of the vibration signal of the physical model and the vibration signal of finite element model respectively
Vector, specially:
Exciting force is applied to the physical model, while acquiring pumping signal and vibration signal, the physical model is produced
Raw vibration signal carries out singular value decomposition, obtains the singular value vector of the vibration signal of physical model;
The pumping signal is loaded onto finite element model, the vibration signal of the finite element model is calculated, to described
The vibration signal of finite element model carries out singular value decomposition, obtains the singular value vector of the vibration signal of finite element model.
As it will be appreciated by persons skilled in the art that in finite element model calculate vibration signal node position with
The position consistency that the measured point of vibration signal is acquired in physical model, ensure that calculating in this way is the vibration letter of the same position
Number.The specific algorithm of singular value decomposition belongs to common knowledge, and this will not be repeated here.
On the basis of the above embodiments, the mistake by between physical model and the corresponding singular value vector of finite element model
Difference is used as objective function, specially:
Using following formula as the objective function:
Wherein, saIndicate the corresponding singular value vector of finite element model, seThe corresponding singular value vector of presentation-entity model,
Error between Δ s presentation-entity model and the corresponding singular value vector of finite element model.As shown from the above formula, the present invention is implemented
For example using the relative error between two singular value vectors as objective function, relative error can more reflect the credible journey of measurement
Degree.
On the basis of the various embodiments described above, according to the building of the physical model of hanging equipment and connection structure and mechanics parameter
Relevant finite element model, specially:
Grid dividing is carried out to the hanging equipment using the lumped mass dot element in finite element analysis software, using even
Order member carries out grid dividing to the connection structure, establishes the finite element model;
It should be noted that the foundation of finite element model is that mechanical structure is converted to multinode to be connected with unit, section
Point is the coordinate of a point in mechanical structure, specifies a number and coordinate position.When node foundation after the completion of, need using
Mechanical structure is connected into unit according to node, and completes finite element model by unit appropriate.Unit selection correctness, will
Determine last analysis result.Finite element analysis software (such as ANSYS) provides the list of certain amount heterogeneity and type
Member, each unit have fixed number, such as LINK is No. 1 unit, and SOLID is No. 45 unit.Cell type is
The unit type of mechanical structure system, such as desk can be made of table-top unit and feet unit, in embodiments of the present invention, will
It hangs equipment and grid dividing is carried out by lumped mass dot element (MASS), regard a lumped mass as equipment is hung using MASS
Point can play simplified model, reduce calculation amount, accelerate the effect of calculating speed;By connection structure by connection unit (BUSH)
Grid dividing is carried out, the rigidity parameters and damping parameter of bindiny mechanism can be defined using BUSH.
The quality of the physical model and rotary inertia are input in the finite element model, as finite element model
Quality and rotary inertia;The value range of mechanics parameter in finite element model is set.
It should be noted that the quality of physical model and rotary inertia are input to finite element model, be in order to construct with
Physical model finite element model as consistent as possible, and the value range of mechanics parameter in finite element model is set, then be in order to
The value bound of mechanics parameter is set, the value of mechanics parameter is constrained when optimizing algorithm.
On the basis of the above embodiments, described that exciting force is applied to physical model, while acquiring pumping signal and vibration
Signal, specially:
Setting at least one measured point passes through acceleration sensing using equipment is hung described in power hammer excitation in hanging equipment
Device acquires the acceleration responsive at least one measured point in 20Hz, and as vibration signal, the measurement direction of acceleration transducer is
Be translatable mode direction, and test frequency range is 2 times for estimating modal frequency, and signal sampling resolution ratio is to estimate the mode frequency range upper limit
2.56 again.
It should be noted that using low frequency since the embodiment of the present invention is concerned with the intrinsic frequency within 20Hz
Acceleration transducer acquires the acceleration responsive in 0~20Hz.Why selection be translatable mode direction, be because the direction be convenient for
It is same convenient for placement acceleration transducer with power hammer excitation.In view of being likely to occur error between calculating and estimation, therefore surveying
Test scope can be amplified when examination, so test frequency range is set as estimating twice of modal frequency.According to aromatic sampling
Theorem, signal sampling resolution ratio be estimate 2 times of mode frequency range, but in practical application, in view of in the frequency band be concerned about without mixed
It is folded, while computer disposal for convenience, it is 2.56 times for estimating the mode frequency range upper limit by signal sampling resolution ratio.
On the basis of the above embodiments, the optimization algorithm is specially ant group algorithm, genetic algorithm, simulated annealing
One of with particle swarm algorithm.It is verified, it, can be in the variation model of mechanics parameter using the optimization efficiency highest of genetic algorithm
Interior carry out global optimizing is enclosed, the precision of parameter identification had both been ensure that, and had in turn ensured the speed of iterative calculation.
On the basis of the above embodiments, the mechanics parameter is stiffness parameters and damping parameter.
Fig. 3 shows a kind of functional block diagram of the computing system of the mechanics parameter of connection structure of the embodiment of the present invention, such as
Shown in Fig. 3, computing system includes:
Finite element model generation module 301, for the physical model according to hanging equipment and connection structure, building and mechanics
The relevant finite element model of parameter.
It should be noted that the physical model of hanging equipment and connection structure hangs equipment and connection structure itself.It hangs
The standby equipment mounted by connection structure is hung, for example crane hangs shot by hook, shot is to hang equipment, hook
Son is bindiny mechanism.The embodiment of the present invention is not intended to limit the type and calculation method of specific hanging equipment and connection structure
Applied field.Finite element model is made of node and element.The finite element model of the embodiment of the present invention is according to hanging equipment
And the complexity of connection structure, it is constructed using two methods:If hanging equipment and connection structure being relatively simple, direct basis
The geometry external form of mechanical structure establishes node and unit;If hanging equipment and connection structure it is complex, by point, line, surface,
Volume first establishes finite element model, then carries out physical grid division, to complete the foundation of finite element model.
Singular value vector obtains module 302, identical sharp for applying respectively to the physical model and finite element model
Power is encouraged, the singular value of the singular value vector of the vibration signal of the physical model and the vibration signal of the finite element model is obtained
Vector.
It should be noted that the acceleration responsive of physical model can be collected by applying exciting force to physical model,
That is vibration signal can get the singular value vector of vibration signal by further carrying out singular value decomposition to vibration signal.Institute
It is because singular value can effectively inhibit vibration signal also to need further to obtain singular value vector after obtaining vibration signal
In interference and background random noise, simultaneously effective enhance vibration signal energy value, and by inventor verifying, vibration
The singular value vector of dynamic signal is very sensitive for the variation of mechanics parameter, i.e., can be more using the singular value vector of vibration signal
Accurately demarcate the mechanics parameter in finite element model.Using the vibration signal of actual measurement physical model, and calculate vibration signal
Singular value vector is identified with the mechanics parameter to finite element dynamics, is provided effectively to the identification of mechanics parameter
Technical support, and due to the vibration signal of physical model be derived from actual physical structure test, finally identify mechanics ginseng
Number also more has convincingness.
Optimization module 303, for using the error between physical model and the corresponding singular value vector of finite element model as mesh
Scalar functions determine mechanics parameter when making the objective function stable convergence using optimization algorithm, as the hanging equipment and
The actual mechanics parameter of connection structure.
It should be noted that the singular value vector due to vibration signal is divided into sensitivity to the variation of mechanics parameter, it will be unusual
It is worth the error between vector as objective function, so that objective function can also be sent out even if small variation has occurred in mechanics parameter
Raw apparent variation.By the value of mechanics parameter in optimization algorithm optimized FEMs model, so that the finite element model after optimization
Corresponding singular value vector singular value vector corresponding with physical model is very close, and the finite element model after reaching optimization can
The effect of actual response physical model.
The computing system that Fig. 4 shows the embodiment of the present invention makes the objective function stable convergence using optimization algorithm determination
When mechanics parameter method flow diagram, as shown,
Step 401, singular value vector obtain module and apply exciting force to physical model, obtain the vibration signal of physical model
Singular value vector, execute step 402;
Step 402, optimization module generate new mechanics parameter using optimization algorithm, execute step 403;
Step 403, finite element model generation module update finite element model using new mechanics parameter, and singular value vector obtains
Exciting force is input in updated finite element model by modulus block, is obtained the singular value vector of new vibration signal, is executed step
Rapid 404;
Step 404, optimization module are unusual by the vibration signal of the singular value vector of new vibration signal and physical model
It is worth vector and substitutes into objective function, executes step 405;
Step 405, optimization module judge objective function whether stable convergence, if so, 406 are thened follow the steps, if it is not, then returning
Return step 402;
Step 406 terminates optimization.
On the basis of the above embodiments, singular value vector obtains module and specifically includes:
Physical model singular value vector acquiring unit for applying exciting force to the physical model, while acquiring excitation
Signal and vibration signal carry out singular value decomposition to the vibration signal that the physical model generates, obtain the vibration of physical model
The singular value vector of signal;
Finite element model singular value vector acquiring unit is counted for the pumping signal to be loaded onto finite element model
The vibration signal for calculating the finite element model carries out singular value decomposition to the vibration signal of the finite element model, obtains limited
The singular value vector of the vibration signal of meta-model.
As it will be appreciated by persons skilled in the art that in finite element model calculate vibration signal node position with
The position consistency that the measured point of vibration signal is acquired in physical model, ensure that calculating in this way is the vibration letter of the same position
Number.The specific algorithm of singular value decomposition belongs to common knowledge, and this will not be repeated here.
On the basis of the above embodiments, the mistake by between physical model and the corresponding singular value vector of finite element model
Difference is used as objective function, specially:
Using following formula as the objective function:
Wherein, saIndicate the corresponding singular value vector of finite element model, seThe corresponding singular value vector of presentation-entity model,
Error between Δ s presentation-entity model and the corresponding singular value vector of finite element model.As shown from the above formula, the present invention is implemented
For example using the relative error between two singular value vectors as objective function, relative error can more reflect the credible journey of measurement
Degree.
On the basis of the above embodiments, finite element model generation module specifically includes:
Finite element generation unit, for utilizing the lumped mass dot element in finite element analysis software to the hanging equipment
Grid dividing is carried out, grid dividing is carried out to the connection structure using connection unit, establishes the finite element model;
It should be noted that the foundation of finite element model is that mechanical structure is converted to multinode to be connected with unit, section
Point is the coordinate of a point in mechanical structure, specifies a number and coordinate position.When node foundation after the completion of, need using
Mechanical structure is connected into unit according to node, and completes finite element model by unit appropriate.Unit selection correctness, will
Determine last analysis result.Finite element analysis software (such as ANSYS) provides the list of certain amount heterogeneity and type
Member, each unit have fixed number, such as LINK is No. 1 unit, and SOLID is No. 45 unit.Cell type is
The unit type of mechanical structure system, such as desk can be made of table-top unit and feet unit, in embodiments of the present invention, will
It hangs equipment and grid dividing is carried out by lumped mass dot element (MASS), regard a lumped mass as equipment is hung using MASS
Point can play simplified model, reduce calculation amount, accelerate the effect of calculating speed;By connection structure by connection unit (BUSH)
Grid dividing is carried out, the rigidity parameters and damping parameter of bindiny mechanism can be defined using BUSH.
Parameter setting unit, for the quality of the physical model and rotary inertia to be input to the finite element model
In, quality and rotary inertia as finite element model;The value range of mechanics parameter in finite element model is set.
It should be noted that the quality of physical model and rotary inertia are input to finite element model, be in order to construct with
Physical model finite element model as consistent as possible, and the value range of mechanics parameter in finite element model is set, then be in order to
The value bound of mechanics parameter is set, the value of mechanics parameter is constrained when optimizing algorithm.
On the basis of the above embodiments, singular value vector, which obtains, applies exciting force to physical model described in module, simultaneously
Pumping signal and vibration signal are acquired, specially:
Setting at least one measured point passes through acceleration sensing using equipment is hung described in power hammer excitation in hanging equipment
Device acquires the acceleration responsive at least one measured point in 20Hz, and as vibration signal, the measurement direction of acceleration transducer is
Be translatable mode direction, and test frequency range is 2 times for estimating modal frequency, and signal sampling resolution ratio is to estimate the mode frequency range upper limit
2.56 again.
It should be noted that using low frequency since the embodiment of the present invention is concerned with the intrinsic frequency within 20Hz
Acceleration transducer acquires the acceleration responsive in 0~20Hz.Why selection be translatable mode direction, be because the direction be convenient for
It is same convenient for placement acceleration transducer with power hammer excitation.In view of being likely to occur error between calculating and estimation, therefore surveying
Test scope can be amplified when examination, so test frequency range is set as estimating twice of modal frequency.According to aromatic sampling
Theorem, signal sampling resolution ratio be estimate 2 times of mode frequency range, but in practical application, in view of in the frequency band be concerned about without mixed
It is folded, while computer disposal for convenience, it is 2.56 times for estimating the mode frequency range upper limit by signal sampling resolution ratio.
Fig. 5 shows the block schematic illustration of electronic equipment of the embodiment of the present invention.
Reference Fig. 5, the electronic equipment, including:Processor (processor) 501, memory (memory) 502 and total
Line 503;
Wherein, the processor 501 and memory 502 complete mutual communication by the bus 503;
The processor 501 is used to call the program instruction in the memory 502, to execute above-mentioned each method embodiment
Provided method, for example including:According to the physical model of hanging equipment and connection structure, building is relevant to mechanics parameter to be had
Limit meta-model;Identical exciting force is applied to the physical model and finite element model respectively, obtains the vibration of the physical model
The singular value vector of the vibration signal of the singular value vector and finite element model of dynamic signal;By physical model and finite element mould
Error between the corresponding singular value vector of type makes the objective function stable convergence using optimization algorithm determination as objective function
When mechanics parameter, as the hanging equipment and the actual mechanics parameter of connection structure.
Another embodiment of the present invention provides a kind of non-transient computer readable storage medium, the non-transient computer is readable
Storage medium stores computer instruction, and the computer instruction executes the computer provided by above-mentioned each method embodiment
Method, for example including:According to the physical model of hanging equipment and connection structure, finite element mould relevant to mechanics parameter is constructed
Type;Identical exciting force is applied to the physical model and finite element model respectively, obtains the vibration signal of the physical model
Singular value vector and the finite element model vibration signal singular value vector;Physical model and finite element model is corresponding
Singular value vector between error as objective function, power when making the objective function stable convergence is determined using optimization algorithm
Parameter is learned, as the hanging equipment and the actual mechanics parameter of connection structure.
In order to more intuitively show the embodiment of the present invention connection structure mechanics parameter calculation method technical effect,
Now an example is illustrated.In this example, it with the mass block of a 2kg, is hung using 4 rubber ropes, simulates motor-car
Equipment and its hanging connection, finite element model are 1 lumped mass dot element, 1 Bush unit, 2 nodes.In this example
Implementation process in, in initial finite element model, stiffness parameters are set as 2190.0N/mm, and damping parameter is set as 10.0Ns/
Mm, recognition result are rigidity 3109.0N/mm, damping parameter 11.7Ns/mm.Fig. 6 shows the receipts of mechanics parameter in this example
Curve graph is held back, Fig. 7 shows the convergence curve figure of objective function in this example, it can be seen from figures 6 and 7 that damping parameter
It tends towards stability with stiffness parameters after 3 iteration, expression has optimized completion, ultrahigh in efficiency.
The apparatus embodiments described above are merely exemplary, wherein described, unit can as illustrated by the separation member
It is physically separated with being or may not be, component shown as a unit may or may not be physics list
Member, it can it is in one place, or may be distributed over multiple network units.It can be selected according to the actual needs
In some or all of the modules achieve the purpose of the solution of this embodiment.Those of ordinary skill in the art are not paying creativeness
Labour in the case where, it can understand and implement.
Through the above description of the embodiments, those skilled in the art can be understood that each embodiment can
It realizes by means of software and necessary general hardware platform, naturally it is also possible to pass through hardware.Based on this understanding, on
Stating technical solution, substantially the part that contributes to existing technology can be embodied in the form of software products in other words, should
Computer software product may be stored in a computer readable storage medium, such as ROM/RAM, magnetic disk, CD, including several fingers
It enables and using so that a computer equipment (can be personal computer, server or the network equipment etc.) executes each implementation
Method described in certain parts of example or embodiment.
Finally it should be noted that:The above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations;Although
Present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that:It still may be used
To modify the technical solutions described in the foregoing embodiments or equivalent replacement of some of the technical features;
And these are modified or replaceed, technical solution of various embodiments of the present invention that it does not separate the essence of the corresponding technical solution spirit and
Range.
Claims (10)
1. a kind of calculation method of the mechanics parameter of connection structure, which is characterized in that including:
According to the physical model of hanging equipment and connection structure, finite element model relevant to mechanics parameter is constructed;
Identical exciting force is applied to the physical model and finite element model respectively, obtains the vibration signal of the physical model
Singular value vector and the finite element model vibration signal singular value vector;
It is true using optimization algorithm using the error between physical model and the corresponding singular value vector of finite element model as objective function
Surely make the mechanics parameter when objective function stable convergence, as the hanging equipment and the actual mechanics ginseng of connection structure
Number.
2. calculation method according to claim 1, which is characterized in that described to the physical model and finite element model point
Do not apply identical exciting force, obtains the singular value vector of the vibration signal of the physical model and the vibration letter of finite element model
Number singular value vector, specially:
Exciting force is applied to the physical model, while acquiring pumping signal and vibration signal, the physical model is generated
Vibration signal carries out singular value decomposition, obtains the singular value vector of the vibration signal of physical model;
The pumping signal is loaded onto finite element model, the vibration signal of the finite element model is calculated, to described limited
The vibration signal of meta-model carries out singular value decomposition, obtains the singular value vector of the vibration signal of finite element model.
3. calculation method according to claim 1, which is characterized in that described that physical model and finite element model is corresponding
Error between singular value vector is as objective function, specially:
Using following formula as the objective function:
Wherein, saIndicate the corresponding singular value vector of finite element model, seThe corresponding singular value vector of presentation-entity model, Δ s table
Show the error between physical model and the corresponding singular value vector of finite element model.
4. calculation method according to claim 1, which is characterized in that the entity according to hanging equipment and connection structure
Model construction finite element model relevant to mechanics parameter, specially:
Grid dividing is carried out to the hanging equipment using the lumped mass dot element in finite element analysis software, utilizes finite element
The connection unit analyzed in software carries out grid dividing to the connection structure, establishes the finite element model;
The quality of the physical model and rotary inertia are input in the finite element model, the quality as finite element model
And rotary inertia;The value range of mechanics parameter in finite element model is set.
5. calculation method according to claim 2, which is characterized in that apply exciting force to the physical model, adopt simultaneously
Collect pumping signal and vibration signal, specially:
Setting at least one measured point is adopted using equipment is hung described in power hammer excitation by acceleration transducer in hanging equipment
Collect the acceleration responsive at least one measured point in 20Hz, as vibration signal;
Wherein, the measurement direction of acceleration transducer is translation mode direction, and test frequency range is 2 times for estimating modal frequency, letter
Number sampling resolution is 2.56 times for estimating the mode frequency range upper limit.
6. calculation method according to claim 1, which is characterized in that the optimization algorithm be ant group algorithm, genetic algorithm,
One of simulated annealing and particle swarm algorithm.
7. calculation method according to claim 1, which is characterized in that the mechanics parameter is stiffness parameters and damping ginseng
Number.
8. a kind of computing system of the mechanics parameter of connection structure, which is characterized in that including:
Finite element model generation module, for the physical model according to hanging equipment and connection structure, building and mechanics parameter phase
The finite element model of pass;
Singular value vector obtains module and obtains for applying identical exciting force respectively to the physical model and finite element model
Obtain the singular value vector of the singular value vector of the vibration signal of the physical model and the vibration signal of the finite element model;
Optimization module, for using the error between physical model and the corresponding singular value vector of finite element model as objective function,
Mechanics parameter when making the objective function stable convergence is determined using optimization algorithm, as the hanging equipment and connection structure
Actual mechanics parameter.
9. a kind of electronic equipment of the mechanics parameter of connection structure, which is characterized in that including:
At least one processor;And
At least one processor being connect with the processor communication, wherein:
The memory is stored with the program instruction that can be executed by the processor, and the processor calls described program to instruct energy
Enough execute method according to any one of claims 1 to 7.
10. a kind of non-transient computer readable storage medium, which is characterized in that the non-transient computer readable storage medium is deposited
Computer instruction is stored up, the computer instruction makes the computer execute method according to any one of claims 1 to 7.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810421619.8A CN108829924B (en) | 2018-05-04 | 2018-05-04 | Method for calculating mechanical parameters of connection structure |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810421619.8A CN108829924B (en) | 2018-05-04 | 2018-05-04 | Method for calculating mechanical parameters of connection structure |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108829924A true CN108829924A (en) | 2018-11-16 |
CN108829924B CN108829924B (en) | 2022-07-08 |
Family
ID=64147447
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810421619.8A Active CN108829924B (en) | 2018-05-04 | 2018-05-04 | Method for calculating mechanical parameters of connection structure |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108829924B (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105424797A (en) * | 2015-11-05 | 2016-03-23 | 北京航空航天大学 | Device and method for performing modal testing on inflatable flexible film structure based on hammering excitation method |
CN106529055A (en) * | 2016-11-18 | 2017-03-22 | 南京航空航天大学 | Model updating method based on strain modal shape correlation |
US20170177769A1 (en) * | 2015-12-22 | 2017-06-22 | Dassault Systemes Simulia Corp. | Effectively Solving Structural Dynamics Problems With Modal Damping In Physical Coordinates |
CN107885908A (en) * | 2017-10-18 | 2018-04-06 | 中车青岛四方机车车辆股份有限公司 | A kind of method for building up of the laminate dynamically equivalent model based on mode of oscillation |
-
2018
- 2018-05-04 CN CN201810421619.8A patent/CN108829924B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105424797A (en) * | 2015-11-05 | 2016-03-23 | 北京航空航天大学 | Device and method for performing modal testing on inflatable flexible film structure based on hammering excitation method |
US20170177769A1 (en) * | 2015-12-22 | 2017-06-22 | Dassault Systemes Simulia Corp. | Effectively Solving Structural Dynamics Problems With Modal Damping In Physical Coordinates |
CN106529055A (en) * | 2016-11-18 | 2017-03-22 | 南京航空航天大学 | Model updating method based on strain modal shape correlation |
CN107885908A (en) * | 2017-10-18 | 2018-04-06 | 中车青岛四方机车车辆股份有限公司 | A kind of method for building up of the laminate dynamically equivalent model based on mode of oscillation |
Non-Patent Citations (3)
Title |
---|
展铭: "基于代理模型的 L 型双加筋板结构", 《中国优秀硕士学位论文全文数据库 (工程科技Ⅱ辑)》 * |
庆光蔚: "基于复模态和频响修正方法的", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
黎定仕 等: "基于结构阻尼的机械结合部动力学模型研究", 《振动与冲击》 * |
Also Published As
Publication number | Publication date |
---|---|
CN108829924B (en) | 2022-07-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106525226B (en) | Evaluation method and system based on-site vibration load recognition | |
CN110146810A (en) | A kind of determination method and device of linear motor test parameter | |
CN107885908A (en) | A kind of method for building up of the laminate dynamically equivalent model based on mode of oscillation | |
CN105631090A (en) | Finite element model optimization device and method | |
CN110287633B (en) | Fatigue analysis method and device for engine aftertreatment bracket | |
CN109362020A (en) | A kind of Numerical Simulation Analysis method of speaker frame dynamic stiffness | |
CN104156549B (en) | The seat cushion equivalent stiffness of seat model of vibration and the discrimination method of equivalent damping | |
CN102810127A (en) | Virtual vibration test system for spacecraft | |
JP2015032295A (en) | Damped vibration analysis method | |
CN104075846B (en) | A kind of based on the rotor unbalance degree recognition methods calculating reverse technology | |
CN108804766A (en) | The recognition methods of damping parameter | |
CN108573084B (en) | Environmental vibration test method and system | |
CN104182592B (en) | The discrimination method of equivalent three mass of human body of seat human body vibrating model | |
CN112149221B (en) | Oil tank fluid-solid coupling transient analysis method and system considering liquid additional mass | |
CN108829924A (en) | The calculation method of the mechanics parameter of connection structure | |
CN102270249B (en) | Method for identifying characteristic frequency of parts | |
CN112129546B (en) | Method and device for testing modal parameters of suspension system, computer equipment and medium | |
Gaonkar | Modal analysis of exhaust system to optimize mounting hanger location | |
CN106599395A (en) | Numerical simulation calculation method for noise of oil immersed transformer | |
CN114912329A (en) | Modeling method and device of battery pack model, electronic equipment and storage medium | |
CN108491644A (en) | A kind of blade of wind-driven generator aerofoil profile critical wind speed of flutter computational methods and equipment | |
CN104408224B (en) | The human body equivalent stiffness of seat human body vibrating model and the discrimination method of damping | |
CN108304653B (en) | Quantitative and rapid evaluation technology for environmental vibration influence of crystal cultivation process | |
CN108090846B (en) | Method and device for constructing power grid low-frequency oscillation case library | |
CN106548137B (en) | Two degree freedom system Identification of Structural Parameters method based on vibration response signal |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |