CN108570553B - A kind of determination method of the oscillating aging excited frequency based on the strain vibration shape - Google Patents
A kind of determination method of the oscillating aging excited frequency based on the strain vibration shape Download PDFInfo
- Publication number
- CN108570553B CN108570553B CN201810282496.4A CN201810282496A CN108570553B CN 108570553 B CN108570553 B CN 108570553B CN 201810282496 A CN201810282496 A CN 201810282496A CN 108570553 B CN108570553 B CN 108570553B
- Authority
- CN
- China
- Prior art keywords
- frequency
- vibration
- component
- strain
- vibration shape
- 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.)
- Expired - Fee Related
Links
Classifications
-
- C—CHEMISTRY; METALLURGY
- C21—METALLURGY OF IRON
- C21D—MODIFYING THE PHYSICAL STRUCTURE OF FERROUS METALS; GENERAL DEVICES FOR HEAT TREATMENT OF FERROUS OR NON-FERROUS METALS OR ALLOYS; MAKING METAL MALLEABLE, e.g. BY DECARBURISATION OR TEMPERING
- C21D10/00—Modifying the physical properties by methods other than heat treatment or deformation
-
- C—CHEMISTRY; METALLURGY
- C21—METALLURGY OF IRON
- C21D—MODIFYING THE PHYSICAL STRUCTURE OF FERROUS METALS; GENERAL DEVICES FOR HEAT TREATMENT OF FERROUS OR NON-FERROUS METALS OR ALLOYS; MAKING METAL MALLEABLE, e.g. BY DECARBURISATION OR TEMPERING
- C21D11/00—Process control or regulation for heat treatments
-
- 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]
-
- 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
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P10/00—Technologies related to metal processing
- Y02P10/20—Recycling
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Organic Chemistry (AREA)
- Crystallography & Structural Chemistry (AREA)
- Mechanical Engineering (AREA)
- Materials Engineering (AREA)
- Metallurgy (AREA)
- Thermal Sciences (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
Based on the determination method of the oscillating aging excited frequency of the strain vibration shape, the surface residual stress distribution including obtaining component；It derives the strain vibration shape and is displaced the relationship of the vibration shape, establish the governing equation of the strain vibration shape；The finite element model that component is established using finite element software carries out numerical value model analysis to component；Obtain each rank strain vibration shape and intrinsic frequency of component；Determine the excited frequency f of oscillating aging0；With (f0‑f0× 20%) Hz is initial excited frequency, carries out frequency sweep vibration to component, obtains excited frequency f of the frequency as oscillating aging when vibration amplitude maximum1.The excited frequency that the determination method of oscillating aging excited frequency proposed by the present invention based on the strain vibration shape determines is using the residual stress distribution of component surface as foundation, and resonance function of strain is criterion, help to obtain ideal effect of vibration stress relief.
Description
Technical field
The present invention relates to vibration aging technology field, a kind of oscillating aging excited frequency based on the strain vibration shape is refered in particular to really
Determine method.
Technical background
Eliminating the residual stress during component fabricates is an important research class in machine building industry field
Topic.Vibration aging technology because have the characteristics that equipment is simple, short processing time, high treating effect and it is residual as being concerned
Residue stress technology for eliminating.Domestic and international experts and scholars have carried out in-depth study to Vibration Aging Process and theory, and obtain
Theoretical result and application example abundant.By it is domestic and international experts and scholars' the study found that influence effect of vibration stress relief pass
Key technological parameter specifically includes that excited frequency, exciting amplitude and exciting time.The selection of excited frequency is mainly used at present
Frequency sweep method carries out frequency sweep vibration to component, find the biggish resonant frequency of component vibration amplitude, then in its resonant frequency or
Vibration stress relief treatment is carried out under person's subresonance frequency.However when determining the excited frequency of oscillating aging by frequency sweep method, usually simultaneously
The residual stress distribution state of material and the regularity of distribution of the vibration shape are not considered, are significantly limited vibration aging on eliminating remnants and are answered
The effect of power.Determine that deficiency existing for oscillating aging excited frequency, the present invention propose a kind of based on the strain vibration shape for frequency sweep method
Oscillating aging excited frequency determination method.
Summary of the invention
Deficiency existing for oscillating aging excited frequency is determined for frequency sweep method, and the present invention proposes a kind of based on the strain vibration shape
The determination method of oscillating aging excited frequency.
Based on the determination method of the oscillating aging excited frequency of the strain vibration shape, comprising the following steps:
(1), the surface residual stress distribution of component is obtained；
(2), it derives the strain vibration shape and is displaced the relationship of the vibration shape, establish the governing equation of the strain vibration shape；
(3), the finite element model that component is established using finite element software carries out numerical value model analysis to component；
(4), each rank strain vibration shape and intrinsic frequency of component are obtained；
(5), on the basis of step (1) and (4), the excited frequency f of oscillating aging is determined0；
(6), with (f0-f0× 20%) Hz is initial excited frequency, carries out frequency sweep vibration to component, obtains vibration amplitude most
Excited frequency f of the frequency as oscillating aging when big1。
The governing equation process of the strain vibration shape is established in step (2) are as follows: based on being displaced the vibration shape, derive the strain vibration shape with
It is displaced the relationship of the vibration shape, establishes the governing equation of the strain vibration shape, is laid the foundation for subsequent numerical value model analysis.One typical
Finite elements are by node j, k, and m and corresponding connection straight line define, then in the unit any point displacement components uiIt can be with
It is expressed as
ui=Niai (1)
In formula: NiIt is the function given in advance based on node location, referred to as shape function matrix, aiIndicate i-th of unit institute
The modal displacement vector being made of modal displacement.
When the displacement of any position in known units, then the strain stress of any position in unit can be determined accordinglyi,
Its expression formula is
εi=Siui (2)
In formula: SiFor Differential operator Matrix.If system is divided into n unit, equation as n formula (2) is shared,
It can be expressed as
It can be abbreviated as
ε=Su (4)
In formula: ε is the strain vector of n unit, and u is the modal displacement vector of n unit.Formula (4) establish strain with
The relationship of displacement is established in unit coordinate system, it is contemplated that it is consecutive variations that each unit is displaced in connecting node, then
Displacement of elemental node vector can be transformed into global coordinate by the coordinate transform of formula (5)
U=Tug (5)
In formula: T is transformation matrix of coordinates, ugFor the modal displacement vector in global coordinate.Formula (5) is updated to
(4) in, the strain Equation in available global coordinate is
Using finite element model, in global coordinate, the undamped motion equation of system is
In formula: ügFor the node acceleration vector in global coordinate, FgFor load vectors.Enable Fg=Fmaxejωt, then node
Motion vector be ug=Umaxejωt, it is available to substitute into equation (7)
(Kg-ω2M)Umax=Fmax (8)
Using mode superposition method, the response of modal displacement can be expressed as
In formula:To be inherently displaced vibration shape matrix, can be expressed asMatrix Y can be expressed as Y
=diag [Y1Y2…Yn], wherein Yr=(kr-ω2mr)-1, krFor r rank modal stiffness, can be expressed asmr
For r rank modal mass, can be expressed as
Equation (9) is updated in equation (6), it is available
It enablesThe available strain vibration shape of differentiating is carried out to the displacement vibration shape, and matrix S is differential operator,
Therefore ψ is known as strain mode vibration shape matrix, ψrFor the r rank strain mode vibration shape, with the r rank displacement modes vibration shapeIt is corresponding.
From the angle of energy, the r rank strain mode vibration shape reaction be component r rank self-energy distribution, with intrinsic frequency and
The displacement modes vibration shape is equally also a parameter for reflecting component inherent characteristic.
In step (5), on the basis of step (1) and (4), the excited frequency f of oscillating aging is determined0Including following step
It is rapid:
(5.1), on the basis of step (1), the region where the larger residual stress of component is determined；
(5.2), on the basis of step (4), the region where strain vibration shape large strain is determined；
(5.3), when strain the region where vibration shape large strain it is consistent with the region where the larger residual stress of component when,
Intrinsic frequency corresponding to the strain vibration shape is recorded, and is denoted as f0。
In step (6), with (f0-f0× 20%) Hz is initial excited frequency, carries out frequency sweep vibration to component, obtains vibration
Excited frequency f of the frequency as oscillating aging when amplitude maximum1The following steps are included:
(6.1), the initial excited frequency f of frequency sweep vibration2It is set as (f0-f0× 20%) Hz, then using 10Hz as step-length by
Step increases the frequency of frequency sweep vibration；The vibration amplitude on component is acted on when recording each excited frequency respectively, is vibrated
Frequency when amplitude maximum, and it is denoted as f21；
(6.2), the initial excited frequency of frequency sweep vibration is set as (f21- 10) Hz is stepped up frequency sweep vibration by step-length of 1Hz
Dynamic frequency；Then the process for repeating step (6.1), obtains frequency when vibration amplitude maximum, and be denoted as f22, when as vibrating
The excited frequency f of effect1。
Further, the finite element software is ANSYS finite element software.
Technical concept of the invention is: the excited frequency that the method proposed through the invention determines is with the residual of component surface
Residue stress is distributed as foundation, and resonate function of strain ψrFor criterion, it help to obtain ideal effect of vibration stress relief.
Beneficial effects of the present invention are as follows:
1, traditional frequency sweep vibratory drilling method determines that the excited frequency of oscillating aging is with 0Hz for initial excited frequency, and this hair
The bright mode by using finite element numerical simulation primarily determines the excited frequency of oscillating aging, then based on this frequency,
The excited frequency that frequency sweep vibrates determining oscillating aging is carried out to component, the time of frequency sweep vibration can be reduced, improve oscillating aging
Working efficiency.
2, it is not intended that the residual stress distribution of component when traditional frequency sweep vibratory drilling method determines the excited frequency of oscillating aging
State, and the determination method of the oscillating aging excited frequency proposed by the present invention based on the strain vibration shape, are with the residual of component surface
Based on residue stress distribution, and the strain vibration shape of combination member, so determining excited frequency can be improved oscillating aging
Effect.
3, the strain vibration shape of component is more sensitive to the localized variation of component compared to the displacement vibration shape, is able to reflect component office
The energy variation situation of portion position reflects the dynamic stress distribution situation of component partial position, therefore based on straining the vibration shape
Determining component oscillating aging excited frequency can be improved the effect of oscillating aging.
Detailed description of the invention
Flow diagram of the Fig. 1 based on the determination method of the oscillating aging excited frequency of the strain vibration shape.
Fig. 2 model analysis component schematic diagram.
The schematic diagram of Fig. 3 Laser Surface Treatment technique.
Fig. 4 residual stress measuring point schematic diagram.
Component 1# residual stress test result before and after Fig. 5 oscillating aging: (a) x-axis is to residual stress；(b) y-axis is answered to remnants
Power.
Component 2# residual stress test result before and after Fig. 6 oscillating aging: (a) x-axis is to residual stress；(b) y-axis is answered to remnants
Power.
Fig. 7 is non-porous and has the low order of pole and high order mode displacement.
Fig. 8 is non-porous and has the low order of pole and high order mode strain.
Specific embodiment
Referring to attached drawing, the present invention is further illustrated:
Based on the determination method of the oscillating aging excited frequency of the strain vibration shape, comprising the following steps:
(1), the surface residual stress distribution of component is obtained；
Specific implementation details are as follows:
45 are selected in experiment#Carbon constructional quality steel prepares experimental member.Experimental member (the i.e. structure of numerical value model analysis
Part) size be 230mm × 210mm × 6mm, as shown in Figure 2.It is answered in order to enable the surface of component generates biggish stretching remnants
Power generates residual stress using Laser Surface Treatment technique in experiment, as shown in Figure 3.This technique is by laser and structure
Part interaction zone quickly heats up to change the microstructure form of component, and the surface of component generates thermal damage and generates just
Begin to stretch residual stress.Laser Surface Treatment technique is unfolded on SLC-2030D planer-type CO2 Multifunction laser processing machine, swashs
The model Rofin SLAB DC030 of light device.
It is as shown in Figure 4 that component residual stress measuring point and strain rosette paste direction.In experiment to two regions of sample into
The edge of row Laser Surface Treatment, the centre distance sample in each laser thermal damage region is 55mm.14 are chosen in experiment
Measuring point is used to evaluate the residual stress of vibration front and back sample, and wherein measuring point 1#~7# is divided into one group, for trying before test vibration
The residual stress of sample, measuring point 1'#~7'# are divided into another group, the residual stress for sample after test vibration.Each group of measuring point
Symmetrically be distributed on the straight line perpendicular to vertical center line, the distance apart from vertical center line be respectively 0mm, 40mm, 60mm,
The distance of 80mm, every group of measuring point test coupon edge are 80mm.By taking measuring point 4 as an example, stickup direction such as Fig. 4 institute of strain rosette
Show, the stickup direction of strain rosette is consistent with the stickup direction of strain rosette at measuring point 4 at remaining measuring point.Pass through strain-gauge test sample
The release strain stress of drill hole1、ε2、ε3, and will release strain substitute into following equation can solve sample drill hole x-axis to
With y-axis to residual stress σxAnd σy
In formula: E is elasticity modulus, ν is Poisson's ratio；For equal two axial stresses calibration factor,For the school of shear stress
Quasi- coefficient.The device of orifice method measurement residual stress includes: that ZDL-II type orifice method surveys residual stress drilling equipment, YE2538 journey
Control statical strain indicator, 1.5 strain rosette of TJ120-1.5- φ.Bore diameter is 2mm, and drilling depth is also 2mm；It is surveyed according to orifice method
Test-object is quasi-, and selection a is 0.192, b 0.482.
Vibration stress relief treatment front part 1# and 2# residual stress test result is as shown in Figure 5 and Figure 6.Component 1# and component 2#
The test result of residual stress shows that the biggish residual stress region of component is laser surface burn processing center before oscillating aging
Region, with the increase far from laser surface Burn Center region distance, residual stress is then gradually decreased.
(2), it derives the strain vibration shape and is displaced the relationship of the vibration shape, establish the governing equation of the strain vibration shape；
(3), the finite element model that component is established using finite element software carries out numerical value model analysis to component；
Specific implementation details are as follows:
Establish physical model: the component of numerical value model analysis is exactly 45 used in experiment#Steel tabular component, long 230mm,
Wide 210mm, thickness 6mm, as shown in Figure 2.In order to study the shadow that the variation of component partial structure is displaced the vibration shape and the strain vibration shape to it
It rings, numerical value model analysis is carried out to the tabular component of two kinds of structures, wherein the center of tabular component 1# is without through-hole, referred to as non-porous structure
Part, in addition there is the through-hole of φ 6mm at the center of a tabular component 2#, is known as having pole.Due to the component three of numerical value model analysis
It is relatively easy to tie up geometrical model, selects the 3-D geometric model for directly establishing component in ANSYS classics gui interface herein.
Material setting: it needs to use 45 when numerical value model analysis#The density p of steel, elastic modulus E, Poisson's ratio ν, specific number
Value is respectively ρ=7850kg/m3, E=200GPa, ν=0.28.
Finite element model: using the SOLID95 unit of 20 nodes to component grid division, its three-dimensional finite element mould is established
Type.In order to improve the precision of solution, in the number difference that three length of test specimen, width and thickness direction setting units divide
It is 55,45 and 3, as shown in Figure 2.Then component grid division can be obtained by the way of scanning (Volume Sweep)
To the grid of rule, the precision of solution can also be improved.In addition, SOLID95 unit has intermediate node, this will make solution
As a result more accurate.
Method for solving: it is solved using the Block Lanczos method (piecemeal Block Lanczos) in ANSYS model analysis module
The displacement modes vibration shape and the strain mode vibration shape of component.Block Lanczos method characteristic value solver is ANSYS model analysis mould
Default solver in block, using Lanczos algorithm.This algorithm can not only guarantee the accuracy solved, can more mention
The speed that height solves.
The result of numerical value model analysis are as follows: the no pole lower mode vibration shape (modal frequency 572Hz)；There is pole low
Rank Mode Shape (modal frequency 571Hz).The no pole high order mode vibration shape (modal frequency 3231Hz)；There is pole height
Rank Mode Shape (modal frequency 3229Hz).Low order that is non-porous and having pole and high order mode are displaced along shown in Fig. 2
The regularity of distribution of mapping path is as shown in Figure 7；Low order that is non-porous and having pole and high order mode strain are reflected along shown in Fig. 2
The regularity of distribution of rays diameter is as shown in Figure 8.From fig.7, it can be seen that non-porous and have the low order of pole and high order mode shift value basic
Overlap, even if drill hole it is non-porous and have pole low order and high order mode shift value do not become significantly
Change, shows to be displaced the vibration shape to the localized variation of structure and insensitive.As seen from Figure 8, low order and high-order non-porous and that have pole
What place of the modal strain value other than drilling substantially overlaped, but low order that is non-porous in drill hole and having pole
Apparent mutation but has occurred with high order mode strain value, shows that straining the vibration shape is sensitive to the localized variation of structure, and
High order mode strain value is more obvious in the mutation of drill hole, this also illustrates why from strain the vibration shape angle determine vibration
The excited frequency of timeliness.
(4), each rank strain vibration shape and intrinsic frequency of component are obtained；
(5), on the basis of step (1) and (4), the excited frequency f of oscillating aging is determined0；
(6), with (f0-f0× 20%) Hz is initial excited frequency, carries out frequency sweep vibration to component, obtains vibration amplitude most
Excited frequency f of the frequency as oscillating aging when big1。
The governing equation process of the strain vibration shape is established in step (2) are as follows: based on being displaced the vibration shape, derive the strain vibration shape with
It is displaced the relationship of the vibration shape, establishes the governing equation of the strain vibration shape, is laid the foundation for subsequent numerical value model analysis.One typical
Finite elements are by node j, k, and m and corresponding connection straight line define, then in the unit any point displacement components uiIt can be with
It is expressed as
ui=Niai (1)
In formula: NiIt is the function given in advance based on node location, referred to as shape function matrix, aiIndicate i-th of unit institute
The modal displacement vector being made of modal displacement.
When the displacement of any position in known units, then the strain stress of any position in unit can be determined accordinglyi,
Its expression formula is
εi=Siui (2)
In formula: SiFor Differential operator Matrix.If system is divided into n unit, equation as n formula (2) is shared,
It can be expressed as
It can be abbreviated as
ε=Su (4)
In formula: ε is the strain vector of n unit, and u is the modal displacement vector of n unit.Formula (4) establish strain with
The relationship of displacement is established in unit coordinate system, it is contemplated that it is consecutive variations that each unit is displaced in connecting node, then
Displacement of elemental node vector can be transformed into global coordinate by the coordinate transform of formula (5)
U=Tug (5)
In formula: T is transformation matrix of coordinates, ugFor the modal displacement vector in global coordinate.Formula (5) is updated to
(4) in, the strain Equation in available global coordinate is
Using finite element model, in global coordinate, the undamped motion equation of system is
In formula: ügFor the node acceleration vector in global coordinate, FgFor load vectors.Enable Fg=Fmaxejωt, then node
Motion vector be ug=Umaxejωt, it is available to substitute into equation (7)
(Kg-ω2M)Umax=Fmax (8)
Using mode superposition method, the response of modal displacement can be expressed as
In formula:To be inherently displaced vibration shape matrix, can be expressed asMatrix Y can be expressed as Y
=diag [Y1Y2…Yn], wherein Yr=(kr-ω2mr)-1, krFor r rank modal stiffness, can be expressed asmrFor
R rank modal mass can be expressed as
Equation (9) is updated in equation (6), it is available
It enablesThe available strain vibration shape of differentiating is carried out to the displacement vibration shape, and matrix S is differential operator,
Therefore ψ is known as strain mode vibration shape matrix, ψrFor the r rank strain mode vibration shape, with the r rank displacement modes vibration shapeIt is corresponding.
From the angle of energy, the r rank strain mode vibration shape reaction be component r rank self-energy distribution, with intrinsic frequency and
The displacement modes vibration shape is equally also a parameter for reflecting component inherent characteristic.
In step (5), on the basis of step (1) and (4), the excited frequency f of oscillating aging is determined0Including following step
It is rapid:
(5.1), on the basis of step (1), the region where the larger residual stress of component is determined；
(5.2), on the basis of step (4), the region where strain vibration shape large strain is determined；
(5.3), when strain the region where vibration shape large strain it is consistent with the region where the larger residual stress of component when,
Intrinsic frequency corresponding to the strain vibration shape is recorded, and is denoted as f0。
Specific implementation details are as follows:
The modal strain of component residual stress test result and Fig. 8 before complex chart 5 and Fig. 6 vibration stress relief treatment is analyzed
As a result it is found that region where high-order strain vibration shape large strain can region where the larger residual stress of covering member, because
The excited frequency of this oscillating aging is 3231Hz, i.e. the oscillating aging that component 2# is determined by the method for numerical value model analysis swashs
Vibration frequency is 3231Hz.
It can be improved the validity of effect of vibration stress relief to verify method proposed by the present invention, while choosing oscillating aging
Excited frequency be 572Hz, i.e. component 1# is by the excited frequency of the oscillating aging of the method for numerical value model analysis determination
572Hz。
In step (6), with (f0-f0× 20%) Hz is initial excited frequency, carries out frequency sweep vibration to component, obtains vibration
Excited frequency f of the frequency as oscillating aging when amplitude maximum1The following steps are included:
(6.1), the initial excited frequency f of frequency sweep vibration2It is set as (f0-f0× 20%) Hz, then using 10Hz as step-length by
Step increases the frequency of frequency sweep vibration；The vibration amplitude on component is acted on when recording each excited frequency respectively, is vibrated
Frequency when amplitude maximum, and it is denoted as f21；
(6.2), the initial excited frequency of frequency sweep vibration is set as (f21- 10) Hz is stepped up frequency sweep vibration by step-length of 1Hz
Dynamic frequency；Then the process for repeating step (6.1), obtains frequency when vibration amplitude maximum, and be denoted as f22, when as vibrating
The excited frequency f of effect1。
Specific implementation details:
It is (572-572 × 20%) Hz that initial excited frequency, which is arranged, for component 1#, then according to step (6.1) and step
(6.2) it is 717Hz that method, which searches out the maximum frequency of vibration amplitude,.
It is (3231-3231 × 20%) that initial excited frequency, which is arranged, for component 2#, then according to step (6.1) and step
(6.2) it is 3311Hz that method, which searches out the maximum frequency of vibration amplitude,.
When the oscillating aging excited frequency that the method proposed through the invention determines vibrates component 1# and component 2#
Effect processing.Residual stress test result after component 1# and component 2# vibration stress relief treatment is as shown in Figure 5 and Figure 6.Component 1# and
Residual stress peak value elimination factor is larger after component 2# vibration stress relief treatment, and place lesser for residual stress is (far from sharp
Optical surface Burn Center region), the elimination factor of residual stress is relatively small after vibration stress relief treatment.As seen from Figure 5, component 1# is passed through
Cross after vibration stress relief treatment its x-axis to y-axis to the burst size of peak value residual stress be respectively 106.8MPa and 131.8MPa；By
Fig. 6 is as it can be seen that component 2# its x-axis after vibration stress relief treatment to y-axis is respectively to the burst size of peak value residual stress
183.6MPa and 226.1MPa shows according to residual stress distribution state and the strain vibration shape regularity of distribution, reasonable selection vibration
The excited frequency of timeliness can obtain ideal Relieving Residual Stress in Libration effect.By Fig. 5 and Fig. 6 as it can be seen that component
There is 2# residual stress after preferably vibration to be homogenized distributed effect, show according to residual stress distribution state and strain vibration shape distribution rule
Rule, the reasonable excited frequency for selecting oscillating aging can obtain ideal oscillating aging homogenizing residual stress effect.It is comprehensive
It is upper described, when carrying out vibration stress relief treatment for component, need the residual stress distribution state and strain vibration according to component surface
The type regularity of distribution, the excited frequency of Rational choice oscillating aging can obtain ideal vibration aging on eliminating and be homogenized residual
Residue stress effect.
Further, the finite element software is ANSYS finite element software.
Content described in this specification embodiment is only enumerating to the way of realization of inventive concept, protection of the invention
Range should not be construed as being limited to the specific forms stated in the embodiments, and protection scope of the present invention is also and in art technology
Personnel conceive according to the present invention it is conceivable that equivalent technologies mean.
Claims (3)
1. the determination method of the oscillating aging excited frequency based on the strain vibration shape, comprising the following steps:
(1), the surface residual stress distribution of component is obtained；
(2), it derives the strain vibration shape and is displaced the relationship of the vibration shape, establish the governing equation of the strain vibration shape；
(3), the finite element model that component is established using finite element software carries out numerical value model analysis to component；
(4), each rank strain vibration shape and intrinsic frequency of component are obtained；
(5), on the basis of step (1), the region where the larger residual stress of component is determined；On the basis of step (4), really
Surely the region where vibration shape large strain is strained；Region and the larger residual stress institute of component where straining vibration shape large strain
Region it is consistent when, record intrinsic frequency corresponding to the strain vibration shape, and be denoted as f0；
(6), with (f0-f0× 20%) Hz is initial excited frequency, frequency sweep vibration is carried out to component, when obtaining vibration amplitude maximum
Excited frequency f of the frequency as oscillating aging1。
2. as described in claim 1 based on the determination method of the oscillating aging excited frequency of the strain vibration shape, it is characterised in that: step
Suddenly in (6), with (f0-f0× 20%) Hz is initial excited frequency, frequency sweep vibration is carried out to component, when obtaining vibration amplitude maximum
Excited frequency f of the frequency as oscillating aging1The following steps are included:
(6.1), the initial excited frequency f of frequency sweep vibration2It is set as (f0-f0× 20%) then Hz gradually increases using 10Hz as step-length
The frequency for adding frequency sweep to vibrate；The vibration amplitude on component is acted on when recording each excited frequency respectively, obtains vibration amplitude
Frequency when maximum, and it is denoted as f21；
(6.2), the initial excited frequency of frequency sweep vibration is set as (f21- 10) Hz is stepped up frequency sweep vibration by step-length of 1Hz
Frequency；Then the process for repeating step (6.1), obtains frequency when vibration amplitude maximum, and be denoted as f22, as oscillating aging
Excited frequency f1。
3. as described in claim 1 based on the determination method of the oscillating aging excited frequency of the strain vibration shape, it is characterised in that: institute
The finite element software stated is ANSYS finite element software.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810282496.4A CN108570553B (en) | 2018-04-02 | 2018-04-02 | A kind of determination method of the oscillating aging excited frequency based on the strain vibration shape |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810282496.4A CN108570553B (en) | 2018-04-02 | 2018-04-02 | A kind of determination method of the oscillating aging excited frequency based on the strain vibration shape |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108570553A CN108570553A (en) | 2018-09-25 |
| CN108570553B true CN108570553B (en) | 2019-08-06 |
Family
ID=63573930
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810282496.4A Expired - Fee Related CN108570553B (en) | 2018-04-02 | 2018-04-02 | A kind of determination method of the oscillating aging excited frequency based on the strain vibration shape |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108570553B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109554532B (en) * | 2018-11-30 | 2020-10-09 | 湖北武铁山桥轨道装备有限公司 | Vibration aging method for reducing and homogenizing residual stress of quenched turnout switch rail |
| CN111325718B (en) * | 2020-01-23 | 2023-07-25 | 深圳大学 | Strain modal analysis method and related device |
| CN113190786B (en) * | 2021-05-13 | 2024-03-15 | 岳聪 | Vibration prediction method for large-scale rotating equipment by utilizing multidimensional assembly parameters |
| CN113836650B (en) * | 2021-08-30 | 2023-04-28 | 北京工业大学 | Modal analysis method based on mechanical structure of laser tracking measurement system |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0889140A1 (en) * | 1997-07-24 | 1999-01-07 | VSR Martin Engineering GmbH | Method of operating a machine for the stress relieving of workpieces |
| CN103488850A (en) * | 2013-10-15 | 2014-01-01 | 北京航空航天大学 | Method for positioning and homogenizing residual stress of thermal shock compounding |
| CN104268342A (en) * | 2014-09-28 | 2015-01-07 | 中航工业哈尔滨轴承有限公司 | Method for analyzing vibration characteristics of bearing retainer based on finite elements |
| CN105779755A (en) * | 2016-04-28 | 2016-07-20 | 上海海事大学 | Quantitative evaluation method for vibration aging effect |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106446443B (en) * | 2016-10-12 | 2020-03-17 | 西南交通大学 | Method and device for identifying resonant frequency of track fastener system |
| CN107038270B (en) * | 2016-10-27 | 2020-01-10 | 华中科技大学 | Method for calculating machining deformation caused by surface machining residual stress field |
-
2018
- 2018-04-02 CN CN201810282496.4A patent/CN108570553B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0889140A1 (en) * | 1997-07-24 | 1999-01-07 | VSR Martin Engineering GmbH | Method of operating a machine for the stress relieving of workpieces |
| CN103488850A (en) * | 2013-10-15 | 2014-01-01 | 北京航空航天大学 | Method for positioning and homogenizing residual stress of thermal shock compounding |
| CN104268342A (en) * | 2014-09-28 | 2015-01-07 | 中航工业哈尔滨轴承有限公司 | Method for analyzing vibration characteristics of bearing retainer based on finite elements |
| CN105779755A (en) * | 2016-04-28 | 2016-07-20 | 上海海事大学 | Quantitative evaluation method for vibration aging effect |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108570553A (en) | 2018-09-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108570553B (en) | A kind of determination method of the oscillating aging excited frequency based on the strain vibration shape | |
| Zhou et al. | Numerical simulation of crack curving and branching in brittle materials under dynamic loads using the extended non-ordinary state-based peridynamics | |
| CN104536941B (en) | A kind of frequency domain load recognition method based on Tikhonov regularizations | |
| CN201392292Y (en) | Tester simulating dynamic load and static load | |
| CN101532931A (en) | Experimental method of simulating dynamic and static load and device thereof | |
| CN105825029B (en) | A method of for optimization design high-frequency vibratory energy amplifying device | |
| Qu et al. | Vertical vibration of piles with square cross‐section | |
| KR101468279B1 (en) | Polycrystalline silicon portion and method for breaking a silicon body | |
| Manzanares-Martínez et al. | Flexural vibrations of a rectangular plate for the lower normal modes | |
| Mead et al. | Measurement of the dynamic moduli and Poisson's ratios of a transversely isotropic fibre-reinforced plastic | |
| CN108195532B (en) | Method for measuring equivalent rigidity of beam structure crack | |
| Wilson et al. | Dynamic behaviour of steel foundations for turbo-alternators | |
| CN107229776A (en) | A kind of diesel engine vibrating isolation system modal analysis method for considering smoke exhaust pipe thermal expansion | |
| Kurahashi et al. | Numerical experiments for cavity shape identification analysis using hammering test data based on level-set type topology optimization | |
| Li et al. | The resonant ultrasound spectroscopy method for determining the Poisson׳ s ratio of spheres over the full range | |
| Sun et al. | Numerical study of monitoring of early-age concrete strength development using PZT induced stress wave | |
| Jingze et al. | Design and mode analysis of random vibration test fixture of some airborne heat exchanger | |
| Zhou et al. | Experimental study on the effect of reinforcement ratio on the ultrasonic velocity and modal frequency of concrete members | |
| Bu et al. | Reduction of dynamic response of a bar using constrained layer damping composite structure | |
| Kannamwar et al. | Free vibration of marine rudder: theoretical and numerical analysis with experimental verification | |
| Wang et al. | Development and Preliminary Application of Resonant Column-Bending Elements Combined Testing Apparatus | |
| Birk et al. | Dynamic response of foundations on three-dimensional layered soil using the scaled boundary finite element method | |
| Song et al. | Modal Experiment and Analysis of HSK Tool-Holder and Thin-Disk Cutter | |
| Zhou et al. | Combined application of three nondestructive testing methods in simulation test of a large concrete-filled steel tube | |
| Deendayal | Finite element analysis of a single pile under cyclic loading |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190806 Termination date: 20210402 |