CN105779755B - A kind of quantitative evaluation methods of effect of vibration stress relief - Google Patents

A kind of quantitative evaluation methods of effect of vibration stress relief Download PDF

Info

Publication number
CN105779755B
CN105779755B CN201610272774.9A CN201610272774A CN105779755B CN 105779755 B CN105779755 B CN 105779755B CN 201610272774 A CN201610272774 A CN 201610272774A CN 105779755 B CN105779755 B CN 105779755B
Authority
CN
China
Prior art keywords
mtd
msub
mtr
mover
effect
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
Application number
CN201610272774.9A
Other languages
Chinese (zh)
Other versions
CN105779755A (en
Inventor
顾邦平
严小兰
胡雄
黄婉娟
赖金涛
张帆
周冶
季霞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Maritime University
Original Assignee
Shanghai Maritime University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Maritime University filed Critical Shanghai Maritime University
Priority to CN201610272774.9A priority Critical patent/CN105779755B/en
Publication of CN105779755A publication Critical patent/CN105779755A/en
Application granted granted Critical
Publication of CN105779755B publication Critical patent/CN105779755B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • CCHEMISTRY; METALLURGY
    • C21METALLURGY OF IRON
    • C21DMODIFYING 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/00Modifying the physical properties by methods other than heat treatment or deformation
    • CCHEMISTRY; METALLURGY
    • C22METALLURGY; FERROUS OR NON-FERROUS ALLOYS; TREATMENT OF ALLOYS OR NON-FERROUS METALS
    • C22FCHANGING THE PHYSICAL STRUCTURE OF NON-FERROUS METALS AND NON-FERROUS ALLOYS
    • C22F3/00Changing the physical structure of non-ferrous metals or alloys by special physical methods, e.g. treatment with neutrons
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/15Correlation function computation including computation of convolution operations

Landscapes

  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Materials Engineering (AREA)
  • Metallurgy (AREA)
  • Organic Chemistry (AREA)
  • Computational Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Mathematical Physics (AREA)
  • Mechanical Engineering (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Algebra (AREA)
  • Thermal Sciences (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Investigating Or Analysing Biological Materials (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

Quantitative evaluation methods for effect of vibration stress relief, set up first by multivariate function approximation theory and consider excited frequency f simultaneously, the quantitative evaluation function model of exciting amplitude A and exciting time t, secondly the technological parameter f during oscillating aging is gathered in real time, A, t, and dimensionless processing is carried out to these three technological parameters, then the dimensionless number of these three technological parameters is substituted into quantitative evaluation function model, can online quantitative evaluation be carried out to effect of vibration stress relief, finally based on quantitative evaluation function model, the partial derivative of single technological parameter and the partial derivative of hybrid technique parameter are asked it respectively, disclose the affecting laws of independent action and reciprocation to Relieving Residual Stress in Libration effect of technological parameter.The present invention has the advantages that the online quantitative evaluation of effect of vibration stress relief can be realized with disclosing the affecting laws of independent action and reciprocation to effect of vibration stress relief of technological parameter.

Description

A kind of quantitative evaluation methods of effect of vibration stress relief
Technical field
The present invention relates to vibration aging technology field, a kind of quantitative evaluation methods of effect of vibration stress relief are refered in particular to.
Technical background
The evaluation of effect of vibration stress relief is one of key issue of vibration aging technology area research, so-called oscillating aging effect Fruit refers to the effect of Relieving Residual Stress in Libration.The evaluation method of current effect of vibration stress relief mainly includes:Parameter curve Method, residual stress measurement method and precision stability detection method.(1) parameter curve method:Parameter curve method is according to oscillating aging mistake The regular variation phenomenon of parameter generation sums up what is come in journey, and the effect of vibration aging technology can comment in linear location Valency, interpretational criteria used mainly includes:The resonance peak of frequency sweep curve before vibration stress relief treatment than rising after vibration stress relief treatment It is high;The resonant frequency of frequency sweep curve before vibration stress relief treatment than reducing after vibration stress relief treatment;Frequency sweep is bent after vibration stress relief treatment Narrow before the bandwidth ratio vibration stress relief treatment of line.(2) residual stress measurement method:Eliminated to accurately evaluate vibration aging technology The effect of residual stress, it is necessary to assess the residual stress before and after component oscillating aging, but remnants of component inside exactly The distribution of stress distribution situation, particularly some complex component internal residual stress, is difficult to adopt the method for theoretical calculation Analysis and solution is carried out, so the residual stress for assessing component inside using the method for experiment measurement has very important reality meaning Justice.The experimental measurement method of current residual stress mainly includes lossless physical measurement method and the mechanical measurement method damaged, can The offline residual stress assessed before and after component vibration stress relief treatment.(3) precision stability detection method:Precision stability detection method It is the effect that vibration aging technology is evaluated by the precision after detection means vibration stress relief treatment, can be to vibration aging technology Effect qualitatively evaluate offline, mainly including long-term placement precision method, add precision method etc. after dynamic loading.It is long-term to place essence Degree method is by the long-term placement of component after vibration stress relief treatment and the dimensional stability of periodic detection component, concrete operations Flow is detected for the first time when placing 15 days, and, total standing time should be more than half a year later detection every 30 days once; Plus precision method is the change of component dimensional accuracy dynamic load effect after of the detection after vibration stress relief treatment after dynamic loading Amount.
Although parameter curve method is a kind of online, directly perceived, quickly and easily method, but this method is when can only be to vibration The effect of effect is qualitatively described, and is unable to the evaluation result of quantitative;Residual stress measurement is owned by France to be commented in offline Valency method, can only be evaluated timeliness effect after oscillating aging process terminates;Precision stability detection method fall within offline, Qualitatively evaluation method, and the assessment cycle that this method needs is very long.Although we are it is seen that these methods can be right The effect of oscillating aging carries out qualitative or quantitative evaluation, but all can not carry out online quantification to the effect of oscillating aging and comment Valency;Simultaneously technological parameter independent action and reciprocation to the affecting laws of effect of vibration stress relief in existing evaluation method body It can not all be disclosed in system, this selection for also resulting in Vibration Aging Process parameter lacks the guidance of the criterion of unified standard. The present invention proposes a kind of quantitative evaluation methods of effect of vibration stress relief, i.e., based on the theoretical oscillating aging effect of multivariate function approximation The quantitative evaluation methods of fruit, this method can carry out online quantitative evaluation to effect of vibration stress relief, can disclose technique ginseng The affecting laws of several independent action and reciprocation to effect of vibration stress relief.
The content of the invention
In order to which the evaluation method for overcoming existing effect of vibration stress relief can not realize the online quantitative evaluation of effect of vibration stress relief With the independent action and reciprocation of technological parameter to the affecting laws of effect of vibration stress relief in existing evaluation method system The deficiency disclosed can not be all able to, the present invention proposes a kind of quantitative evaluation methods of effect of vibration stress relief, i.e., based on polynary letter The quantitative evaluation methods of the effect of vibration stress relief of number approximation theory.
For the quantitative evaluation methods of effect of vibration stress relief, comprise the following steps:
(1) technological parameter of influence Relieving Residual Stress in Libration effect, is determined, mainly includes excited frequency f, exciting Amplitude A and exciting time t, these three technological parameters are the inputs of quantitative evaluation function model;Vibration aging on eliminating is residual The effect F of residue stress is the output end of quantitative evaluation function model;Setting up quantitative evaluation function model F, (f, A t), that is, are built Vertical technological parameter input and the quantitative evaluation function model of effect of vibration stress relief output end.
(2) it is theoretical using multivariate function approximation, by the multinomial function of f, A and t these three technological parameters composition Approach set up out quantitative evaluation function model F (f, A, t).
(3), excited frequency f unit be kHz, exciting amplitude A unit for μm, exciting time t unit be min, be The influences of the different dimensions of elimination, defines following three dimensionless variable:Dimensionless excited frequencyDimensionless swashs AmplitudeThe dimensionless exciting timeWherein f0For 1kHz, A0For 1 μm, t0For 1min.
(4), (f, A can t) be approached quantitative evaluation function model F with following multinomial:In formula:a1、a2、…、a11For real coefficient, Wherein a11For the constant term of equation.
(5), when sample is without vibration stress relief treatment, i.e., when the value of three technological parameters is 0, the remnants of sample should Power eradicating efficacy should be 0, so constant term a11=0.Therefore, quantitative evaluation function model can be further represented as:
(6), the effect of vibration stress relief F by 10 groups of different dimensionless combination of process parameters and corresponding to them substitutes into fixed Quantitatively evaluating function model can be obtained on undetermined coefficient a1、a2、…、a1010 yuan of 1 equation of n th order n group, by solving this equation Group can solve the numerical value of this 10 undetermined coefficients.This 10 yuan of 1 equation of n th order n group can be expressed as using the form of matrix:
It can be further represented as:Ga=F.
In formula:For dimensionless combination of process parameters matrix, the i-th row in the matrix The combination multinomial of dimensionless technological parameter during expression ith vibration stress relief treatment;F=[Fi] be oscillating aging effect square Battle array;A=[ai] it is undetermined coefficient matrix.
(7), in order to solve 10 undetermined coefficients, 10 groups of different experimental datas are at least needed, in order that obtaining step (6) In matrix equation have unique solution, then dimensionless combination of process parameters matrix G order should meet r (G)=10, i.e. dimensionless work Skill parameter combination matrix should be non-singular matrix.Therefore, undetermined coefficient can be solved by following matrix equation and obtained:
By solve 10 The quantitative evaluation function model that undetermined coefficient is substituted into step (5), you can with set up out quantitative evaluation function model F (f, A,t)。
(8), technological parameter f, A, t in real time during collection oscillating aging, and these three technological parameters are carried out immeasurable Guiding principle processing, then substitutes into the quantitative evaluation function model that step (7) is set up by the dimensionless number of these three technological parameters, Can online quantitative evaluation be carried out to effect of vibration stress relief;The quantitative evaluation Function Modules set up simultaneously with step (7) Based on type, partial derivative of single technological parameter and the partial derivative of hybrid technique parameter are asked it respectively, disclose technological parameter The affecting laws of independent action and reciprocation to Relieving Residual Stress in Libration effect.
Described multinomial function includes:Single factor test first order f, A and t, represent that each technological parameter is independent respectively Influence to effect of vibration stress relief;Single factor test quadratic term f2、A2And t2, quadratic term and the vibration of each technological parameter are represented respectively Interactively between timeliness effect;Confounding factor quadratic term fA, ft and At, represent any two in three technological parameters respectively Influence of the reciprocation to effect of vibration stress relief between individual parameter;Confounding factor cubic term fAt, represents all three technique ginsengs Influence of the reciprocation to effect of vibration stress relief between number.
Described single technological parameter or the partial derivative of hybrid technique parameter be on the occasion of, then it represents that vibration aging on eliminating is remaining The effect of stress is improved with the increase of single technological parameter or hybrid technique parameter.
Described single technological parameter or the partial derivative of hybrid technique parameter are negative value, then it represents that vibration aging on eliminating is remaining The effect of stress is reduced with the increase of single technological parameter or hybrid technique parameter.
Described single technological parameter or the absolute value of the partial derivative of hybrid technique parameter are bigger, then it represents that single technique ginseng Number or hybrid technique parameter are bigger to the influential effect of Relieving Residual Stress in Libration.
The present invention technical concept be:Set up first by multivariate function approximation theory and consider excited frequency f, exciting simultaneously The quantitative evaluation function model of amplitude A and exciting time t, secondly gathers the technological parameter f during oscillating aging in real time, A, t, and dimensionless processing is carried out to these three technological parameters, then the dimensionless number of these three technological parameters is substituted into quantitative Change evaluation function model, you can to carry out online quantitative evaluation to effect of vibration stress relief, finally with quantitative evaluation Function Modules Based on type, partial derivative of single technological parameter and the partial derivative of hybrid technique parameter are asked it respectively, disclose technological parameter The affecting laws of independent action and reciprocation to Relieving Residual Stress in Libration effect.
The beneficial effects of the invention are as follows:
1st, it can be set up by multivariate function approximation theory while considering excited frequency f, exciting amplitude A and exciting time T quantitative evaluation function model, can realize and carry out online quantitative evaluation to effect of vibration stress relief.
2nd, based on quantitative evaluation function model, the partial derivative and hybrid technique of single technological parameter is sought it respectively The partial derivative of parameter, can disclose the independent action and reciprocation of technological parameter to Relieving Residual Stress in Libration effect Affecting laws.
Brief description of the drawings
The schematic diagram of Fig. 1 quantitative evaluation function models
The scale diagrams of Fig. 2 Cr12MoV steel curved beams
Embodiment
Referring to the drawings, the present invention is further illustrated:
For the quantitative evaluation methods of effect of vibration stress relief, comprise the following steps:
(1) technological parameter of influence Relieving Residual Stress in Libration effect, is determined, mainly includes excited frequency f, exciting Amplitude A and exciting time t, these three technological parameters are the inputs of quantitative evaluation function model;Vibration aging on eliminating is residual The effect F of residue stress is the output end of quantitative evaluation function model;Setting up quantitative evaluation function model F, (f, A t), that is, are built Vertical technological parameter input and the quantitative evaluation function model of effect of vibration stress relief output end.
Specifically implementation detail is:Quantitative evaluation function model schematic diagram such as Fig. 1 of Relieving Residual Stress in Libration It is shown.As can be known from Fig. 1, the technological parameter of influence Relieving Residual Stress in Libration effect mainly includes excited frequency f, exciting Amplitude A and exciting time t, these three technological parameters are the inputs of quantitative evaluation function model, and vibration aging on eliminating is residual The effect F of residue stress is the output end of quantitative evaluation function model.It is contemplated that setting up quantitative evaluation function model F (f, A, t), to evaluate it is determined that the effect F of the Relieving Residual Stress in Libration under the conditions of combination of process parameters, that is, set up Technological parameter input and the quantitative evaluation function model of effect of vibration stress relief output end.
Assuming that (f, A, domain of definition t) is bounded closed region D to function F, and is real-valued continuous function in the D of the region.P1 (f,A,t)、P2(f,A,t)、…、Pk(f, A are t) that the real-valued continuous function of k linear independence being defined on D (is usually taken to be Multinomial function), wherein by P1(f,A,t)、P2(f,A,t)、…、Pk(f, A t) combine the linear space P made and claimed For linear interpolation space, then arbitrary function can be expressed as in P:P(f,A,t)α1P1(f,A,t)+α2P2(f,A,t)+…+αkPk(f, A, t), in formula:α1、α2、…、αkFor real coefficient.It can be seen from Weierstrass approximation theorems, if function P (f, A, t) Can for approximating function F (f, A, t), then for any given ε, have | F (f, A, t)-P (f, A, t) | < ε, wherein f, A, t ∈D.Therefore, it is contemplated that by the multinomial function of f, A and t these three technological parameters composition, using the function of many variables Approximation theory approach set up out quantitative evaluation function model F (f, A, t).
(2) it is theoretical using multivariate function approximation, by the multinomial function of f, A and t these three technological parameters composition Approach set up out quantitative evaluation function model F (f, A, t).
(3), excited frequency f unit be kHz, exciting amplitude A unit for μm, exciting time t unit be min, be The influences of the different dimensions of elimination, defines following three dimensionless variable:Dimensionless excited frequencyDimensionless swashs AmplitudeThe dimensionless exciting timeWherein f0For 1kHz, A0For 1 μm, t0For 1min.
(4), (f, A can t) be approached quantitative evaluation function model F with following multinomial:In formula:a1、a2、…、a11For real coefficient, Wherein a11For the constant term of equation.
(5), when sample is without vibration stress relief treatment, i.e., when the value of three technological parameters is 0, the remnants of sample should Power eradicating efficacy should be 0, so constant term a11=0.Therefore, quantitative evaluation function model can be further represented as:
(6), the effect of vibration stress relief F by 10 groups of different dimensionless combination of process parameters and corresponding to them substitutes into fixed Quantitatively evaluating function model can be obtained on undetermined coefficient a1、a2、…、a1010 yuan of 1 equation of n th order n group, by solving this equation Group can solve the numerical value of this 10 undetermined coefficients.This 10 yuan of 1 equation of n th order n group can be expressed as using the form of matrix:
It can be further represented as:Ga=F.
In formula:For dimensionless combination of process parameters matrix, the i-th row in the matrix The combination multinomial of dimensionless technological parameter during expression ith vibration stress relief treatment;F=[Fi] be oscillating aging effect square Battle array;A=[ai] it is undetermined coefficient matrix.
(7), in order to solve 10 undetermined coefficients, 10 groups of different experimental datas are at least needed, in order that obtaining step (6) In matrix equation have unique solution, then dimensionless combination of process parameters matrix G order should meet r (G)=10, i.e. dimensionless work Skill parameter combination matrix should be non-singular matrix.Therefore, undetermined coefficient can be solved by following matrix equation and obtained:
Solve 10 are treated Determine the quantitative evaluation function model that coefficient is substituted into step (5), you can with set up out quantitative evaluation function model F (f, A, t)。
Specifically implementation detail is:In order to obtain the experimental data set up needed for quantitative evaluation function model, to Fig. 2 institutes The sample shown carries out vibration stress relief treatment.10 in quantitative evaluation function model in order to solve step (5) foundation Individual undetermined coefficient, it is necessary to select 10 groups of different experimental datas, and dimensionless combination of process parameters matrix is necessary for non-singular matrix, I.e. dimensionless combination of process parameters rank of matrix is 10.In view of there is certain error in experimental data test process, in order to Improving the quantitative evaluation function model set up out has higher reliability and generality, and choosing 11 groups of experimental datas is used to build Halt quantitatively evaluating function model.Choosing 10 groups therein from 11 groups of experimental datas successively is used to construct dimensionless technological parameter Combinatorial matrix G, and provide rj(G) to remove after the jth group data in 11 groups of experimental datas, the dimensionless technique ginseng constructed Number combinatorial matrix G order, j=1~11.Calculate constructed dimensionless work successively using the rank functions in Matlab softwares Skill parameter combination matrix G order, if solving result meets r1(G)=r2(G)=...=r11(G)=10, the then reality of description selection Test data to can be used for solving 11 groups of undetermined coefficients, if solving result is unsatisfactory for r1(G)=r2(G)=...=r11(G)=10 During this condition, then need to reselect experimental data.
Solve 10 undetermined coefficients during, first every time chosen from 11 groups of experimental datas 10 set constructor therein without Dimension combination of process parameters matrix G;Secondly dimensionless combination of process parameters square is solved using the inv functions in Matlab softwares Battle array G inverse matrix G-1, and by matrix G in Matlab softwares-1Being multiplied with the effect matrix F of oscillating aging can solve Undetermined coefficient matrix, can so obtain 11 groups of different undetermined coefficients;Then it is taken to be averaged respectively this 11 groups of undetermined coefficients It is worth the undetermined coefficient as quantitative evaluation function model;The average of these undetermined coefficients is finally substituted into what step (5) was set up Quantitative evaluation function model is that can obtain quantitative evaluation function model.
Table 1 is used for the dimensionless technological parameter and timeliness effect for calculating undetermined coefficient matrix
The each group dimensionless technological parameter of table 2 and its corresponding multinomial numerical value
11 groups of undetermined coefficients are solved using the method described in step (7), 3 are the results are shown in Table.Last of table 3 is classified as 11 groups and treated Determine the average of coefficient, the average of these undetermined coefficients is substituted into the quantitative evaluation function model that step (5) is set up, then finally built Vertical quantitative evaluation function model can be expressed as:
The numerical value for each group undetermined coefficient that table 3 is solved
(8), technological parameter f, A, t in real time during collection oscillating aging, and these three technological parameters are carried out immeasurable Guiding principle processing, then substitutes into the quantitative evaluation function model that step (7) is set up by the dimensionless number of these three technological parameters, Can online quantitative evaluation be carried out to effect of vibration stress relief;The quantitative evaluation Function Modules set up simultaneously with step (7) Based on type, partial derivative of single technological parameter and the partial derivative of hybrid technique parameter are asked it respectively, disclose technological parameter The affecting laws of independent action and reciprocation to Relieving Residual Stress in Libration effect.
Described multinomial function includes:Single factor test first order f, A and t, represent that each technological parameter is independent respectively Influence to effect of vibration stress relief;Single factor test quadratic term f2、A2And t2, quadratic term and the vibration of each technological parameter are represented respectively Interactively between timeliness effect;Confounding factor quadratic term fA, ft and At, represent any two in three technological parameters respectively Influence of the reciprocation to effect of vibration stress relief between individual parameter;Confounding factor cubic term fAt, represents all three technique ginsengs Influence of the reciprocation to effect of vibration stress relief between number.
Described single technological parameter or the partial derivative of hybrid technique parameter be on the occasion of, then it represents that vibration aging on eliminating is remaining The effect of stress is improved with the increase of single technological parameter or hybrid technique parameter.
Described single technological parameter or the partial derivative of hybrid technique parameter are negative value, then it represents that vibration aging on eliminating is remaining The effect of stress is reduced with the increase of single technological parameter or hybrid technique parameter.
Described single technological parameter or the absolute value of the partial derivative of hybrid technique parameter are bigger, then it represents that single technique ginseng Number or hybrid technique parameter are bigger to the influential effect of Relieving Residual Stress in Libration.
Specifically implementation detail is:Technological parameter f, A, t during collection oscillating aging, and to these three techniques in real time Parameter carries out dimensionless processing, and the dimensionless number of these three technological parameters then is substituted into the quantification that step (7) is set up Evaluation function model, you can to carry out online quantitative evaluation to effect of vibration stress relief.
By quantitative evaluation function model respectively to dimensionless excited frequencyDimensionless exciting amplitudeSwash with dimensionless Shake the timePartial derivative is sought, and the average undetermined coefficient in table 3 is substituted into the expression formula of partial derivative, quantification can be solved Partial derivative expression formula of the evaluation function model to independent process parameter
By the 11 groups of dimensionless technological parameters and its polynomial value in table 2 substitute into respectively partial derivative expression formula (1), (2), (3) calculated in, the results are shown in Table last column in 4, table is the equal of the concrete numerical value of the partial derivative of each independent process parameter Value.From table 4 it was found from the average of the partial derivative of last column,WithIt is just Value, shows that, with excited frequency f, exciting amplitude A and exciting time t increase, the effect of Relieving Residual Stress in Libration is got over Significantly.It is also found that the average of the partial derivative of these three independent process parameters meets following relation from table 4
Show the effect of oscillating aging is influenceed maximum by excited frequency f, next to that exciting amplitude A, by exciting time t's Influence is minimum.
Concrete numerical value and average of the quantitative evaluation function model of table 4 to the partial derivative of each independent process parameter
Quantitative evaluation function model is sought into three technological parameters second-order mixed partial derivative respectively, and will be flat in table 3 Equal undetermined coefficient is substituted into the expression formula of partial derivative, can solve quantitative evaluation function model to the inclined of hybrid technique parameter Derivative expressions
By the 11 groups of dimensionless technological parameters and its polynomial value in table 2 substitute into respectively partial derivative expression formula (5), (6), (7) calculated in, the results are shown in Table last column in 5, table is the equal of the concrete numerical value of the mixed partial derivative of three technological parameters Value.
Concrete numerical value and average of the quantitative evaluation function model of table 5 to the mixed partial derivative of technological parameter
From table 5 it was found from the average of the partial derivative of last column,WithBe on the occasion of, andFor negative value, show excited frequency f and reciprocation, excited frequency f and the exciting time t of exciting amplitude A friendship Interaction is all conducive to improving the effect of Relieving Residual Stress in Libration, and exciting amplitude A and exciting time t reciprocation It is unfavorable for improving the effect of Relieving Residual Stress in Libration.It is also found that the second order mixing of three technological parameters from table 5 The absolute value of partial derivative meets following relation
Show the effect of oscillating aging is influenceed maximum by excited frequency f and the interactive of exciting amplitude A, next to that swashing Vibration frequency f and exciting time t reciprocation, minimum is influenceed by exciting amplitude A and the interactive of exciting time t.
Content described in this specification embodiment is only enumerating to the way of realization of inventive concept, protection of the invention Scope is not construed as being only limitted to the concrete form that embodiment is stated, protection scope of the present invention is also and in art technology Personnel according to present inventive concept it is conceivable that equivalent technologies mean.

Claims (5)

1. for the quantitative evaluation methods of effect of vibration stress relief, comprise the following steps:
(1) technological parameter of influence Relieving Residual Stress in Libration effect, is determined, mainly includes excited frequency f, exciting amplitude A and exciting time t, these three technological parameters are the inputs of quantitative evaluation function model;Vibration aging on eliminating remnants should The effect F of power is the output end of quantitative evaluation function model;Setting up quantitative evaluation function model F, (f, A t), that is, set up work The quantitative evaluation function model of skill parameter input end and effect of vibration stress relief output end;
(2), by the multinomial function of f, A and t these three technological parameters composition, approached using multivariate function approximation theory Set up out quantitative evaluation function model F (f, A, t);
(3), excited frequency f unit be kHz, exciting amplitude A unit for μm, exciting time t unit be min, in order to disappear Except the influence of different dimensions, following three dimensionless variable is defined:Dimensionless excited frequencyDimensionless exciting shakes WidthThe dimensionless exciting timeWherein f0For 1kHz, A0For 1 μm, t0For 1min;
(4), (f, A t) are approached quantitative evaluation function model F with following multinomial:In formula:a1、a2、…、a11For real coefficient, Wherein a11For the constant term of equation;
(5), when sample is without vibration stress relief treatment, i.e., when the value of three technological parameters is 0, the residual stress of sample disappears Except effect should be 0, so constant term a11=0, therefore, quantitative evaluation function model is further represented as:
(6), the effect of vibration stress relief F by 10 groups of different dimensionless combination of process parameters and corresponding to them substitutes into quantification Evaluation function model can be obtained on undetermined coefficient a1、a2、…、a1010 yuan of 1 equation of n th order n group, by solving this equation group i.e. The numerical value of this 10 undetermined coefficients can be solved, this 10 yuan of 1 equation of n th order n group is expressed as in the form of matrix:
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mrow> <mo>(</mo> <mover> <mi>f</mi> <mo>&amp;OverBar;</mo> </mover> <mo>)</mo> </mrow> <mn>1</mn> </msub> </mtd> <mtd> <msub> <mrow> <mo>(</mo> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> <mo>)</mo> </mrow> <mn>1</mn> </msub> </mtd> <mtd> <msub> <mrow> <mo>(</mo> <mover> <mi>t</mi> <mo>&amp;OverBar;</mo> </mover> <mo>)</mo> </mrow> <mn>1</mn> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mrow> <mo>(</mo> <mover> <mi>f</mi> <mo>&amp;OverBar;</mo> </mover> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> <mover> <mi>t</mi> <mo>&amp;OverBar;</mo> </mover> <mo>)</mo> </mrow> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mrow> <mo>(</mo> <mover> <mi>f</mi> <mo>&amp;OverBar;</mo> </mover> <mo>)</mo> </mrow> <mn>2</mn> </msub> </mtd> <mtd> <msub> <mrow> <mo>(</mo> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> <mo>)</mo> </mrow> <mn>2</mn> </msub> </mtd> <mtd> <msub> <mrow> <mo>(</mo> <mover> <mi>t</mi> <mo>&amp;OverBar;</mo> </mover> <mo>)</mo> </mrow> <mn>2</mn> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mrow> <mo>(</mo> <mover> <mi>f</mi> <mo>&amp;OverBar;</mo> </mover> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> <mover> <mi>t</mi> <mo>&amp;OverBar;</mo> </mover> <mo>)</mo> </mrow> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mrow></mrow> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mrow></mrow> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mrow></mrow> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mrow> <mo>(</mo> <mover> <mi>f</mi> <mo>&amp;OverBar;</mo> </mover> <mo>)</mo> </mrow> <mn>10</mn> </msub> </mtd> <mtd> <msub> <mrow> <mo>(</mo> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> <mo>)</mo> </mrow> <mn>10</mn> </msub> </mtd> <mtd> <msub> <mrow> <mo>(</mo> <mover> <mi>t</mi> <mo>&amp;OverBar;</mo> </mover> <mo>)</mo> </mrow> <mn>10</mn> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mrow> <mo>(</mo> <mover> <mi>f</mi> <mo>&amp;OverBar;</mo> </mover> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> <mover> <mi>t</mi> <mo>&amp;OverBar;</mo> </mover> <mo>)</mo> </mrow> <mn>10</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>a</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>a</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>a</mi> <mn>10</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>F</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>F</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>F</mi> <mn>10</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>.</mo> </mrow>
It is further represented as:Ga=F,
In formula:For dimensionless combination of process parameters matrix, the i-th row in the matrix is represented The combination multinomial of dimensionless technological parameter during ith vibration stress relief treatment;F=[Fi] be oscillating aging effect matrix;a =[ai] it is undetermined coefficient matrix;
(7), in order to solve 10 undetermined coefficients, 10 groups of different experimental datas are at least needed, in order that obtaining in step (6) Matrix equation has unique solution, then dimensionless combination of process parameters matrix G order should meet r (G)=10, i.e. dimensionless technique and join Number combinatorial matrix should be non-singular matrix, therefore, and undetermined coefficient is solved by following matrix equation and obtained:
By solve 10 systems undetermined The quantitative evaluation function models that number is substituted into steps (5), you can with set up out quantitative evaluation function model F (f, A, t);
(8), technological parameter f, A, t in real time during collection oscillating aging, and these three technological parameters are carried out at dimensionless Reason, then substitutes into the quantitative evaluation function model that step (7) is set up by the dimensionless number of these three technological parameters, you can To carry out online quantitative evaluation to effect of vibration stress relief;The quantitative evaluation function model set up simultaneously using step (7) as Basis, asks it partial derivative of single technological parameter and the partial derivative of hybrid technique parameter respectively, discloses the independence of technological parameter The affecting laws of effect and reciprocation to Relieving Residual Stress in Libration effect.
2. it is used for the quantitative evaluation methods of effect of vibration stress relief as claimed in claim 1, it is characterised in that:Described is polynary Polynomial function includes:Single factor test first order f, A and t, represent shadow of each technological parameter individually to effect of vibration stress relief respectively Ring;Single factor test quadratic term f2、A2And t2, the effect between the quadratic term of each technological parameter and effect of vibration stress relief is represented respectively Relation;Confounding factor quadratic term fA, ft and At, represent the interaction work between any two parameter in three technological parameters respectively With the influence to effect of vibration stress relief;Confounding factor cubic term fAt, represents the reciprocation pair between all three technological parameters The influence of effect of vibration stress relief.
3. it is used for the quantitative evaluation methods of effect of vibration stress relief as claimed in claim 1, it is characterised in that:Described is single Technological parameter or the partial derivative of hybrid technique parameter be on the occasion of, then it represents that the effect of Relieving Residual Stress in Libration is with single The increase of technological parameter or hybrid technique parameter and improve.
4. it is used for the quantitative evaluation methods of effect of vibration stress relief as claimed in claim 1, it is characterised in that:Described is single Technological parameter or the partial derivative of hybrid technique parameter are negative value, then it represents that the effect of Relieving Residual Stress in Libration is with single The increase of technological parameter or hybrid technique parameter and reduce.
5. it is used for the quantitative evaluation methods of effect of vibration stress relief as claimed in claim 1, it is characterised in that:Described is single The absolute value of technological parameter or the partial derivative of hybrid technique parameter is bigger, then it represents that single technological parameter or hybrid technique parameter pair The influential effect of Relieving Residual Stress in Libration is bigger.
CN201610272774.9A 2016-04-28 2016-04-28 A kind of quantitative evaluation methods of effect of vibration stress relief Expired - Fee Related CN105779755B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610272774.9A CN105779755B (en) 2016-04-28 2016-04-28 A kind of quantitative evaluation methods of effect of vibration stress relief

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610272774.9A CN105779755B (en) 2016-04-28 2016-04-28 A kind of quantitative evaluation methods of effect of vibration stress relief

Publications (2)

Publication Number Publication Date
CN105779755A CN105779755A (en) 2016-07-20
CN105779755B true CN105779755B (en) 2017-10-31

Family

ID=56399802

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610272774.9A Expired - Fee Related CN105779755B (en) 2016-04-28 2016-04-28 A kind of quantitative evaluation methods of effect of vibration stress relief

Country Status (1)

Country Link
CN (1) CN105779755B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106555046B (en) * 2016-11-30 2018-02-23 北京石油化工学院 A kind of method for eliminating welding residual stress
CN108570553B (en) * 2018-04-02 2019-08-06 上海海事大学 A kind of determination method of the oscillating aging excited frequency based on the strain vibration shape
CN109633284B (en) * 2019-01-18 2023-02-10 上海海事大学 Electric pulse aging electric pulse waveform determining system based on data mining technology
CN111730229B (en) * 2020-06-17 2021-11-30 上海海事大学 Method for determining vibration welding technological parameters

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE59700079D1 (en) * 1997-07-24 1999-03-25 Vsr Martin Eng Gmbh Method for operating a machine for the relaxation of workpieces
US7703325B2 (en) * 2007-08-24 2010-04-27 Weite Wu Method for relieving residual stress in an object
CN101691632B (en) * 2009-06-26 2012-11-14 哈尔滨工业大学深圳研究生院 Method for rapidly judging peak
CN101994001A (en) * 2010-09-14 2011-03-30 上海海事大学 Support vector machine algorithm based method for predicting vibration aging effect
CN101967553A (en) * 2010-10-12 2011-02-09 北京翔博科技有限责任公司 Three-dimensional vibration stress relief multitask acquisition and control system and method thereof
CN103589855B (en) * 2013-12-03 2015-01-07 北京航空航天大学 Low temperature treatment-vibration aging combined residual stress homogenization method
CN103757197B (en) * 2014-01-25 2015-09-30 浙江大学 For eliminating high-frequency vibration aging system and the method for small size component unrelieved stress

Also Published As

Publication number Publication date
CN105779755A (en) 2016-07-20

Similar Documents

Publication Publication Date Title
CN105779755B (en) A kind of quantitative evaluation methods of effect of vibration stress relief
Breusch et al. On the fixed-effects vector decomposition
Anastasiou et al. Unitarity cuts and reduction to master integrals in d dimensions for one-loop amplitudes
Tatari et al. The use of the Adomian decomposition method for solving multipoint boundary value problems
CN104730025B (en) Mixture quantitative analysis method based on terahertz spectroscopy
CN112131685B (en) Uncertainty evaluation method for overall thermodynamic performance test of combined cycle unit
CN103500281B (en) Nonlinear system modeling method of process of sugar boiling crystallization
CN105891321A (en) Calibration method for micromagnetic detection of ferromagnetic material structural mechanics performance
CN114065919A (en) Deficiency value completion method and medium based on generation countermeasure network
CN106874251A (en) A kind of measuring method of rock mass RBI and RQD indexs
CN113872239A (en) Knowledge and data combined drive-based grid-connected inverter impedance obtaining method
EP3692362A1 (en) Method for determining the geometry of a defect and for determining a load limit
CN117874955A (en) Gear life fatigue prediction method based on finite element analysis
CN103197692B (en) A kind of flow control methods, Apparatus and system
Moon et al. On the Breitung test for panel unit roots and local asymptotic power
Rogov et al. Fourth-order accurate bicompact schemes for hyperbolic equations
CN109388780A (en) A kind of Uncertainty Analysis Method of falling film evaporator LCA model
Belov Numerical diagnostics of solution blowup in differential equations
WO2019001902A1 (en) Method for creating an evaluation table for an ultrasonic inspection and method for ultrasonic inspection
Sahoo et al. Voltage stability monitoring based on feed forward and layer recurrent neural networks
CN205662571U (en) Online quantification evaluation system of vibration prescription effect
Filaretov et al. Generalized parameter extraction method for analog circuit fault diagnosis
CN106644977A (en) Spectral variable selection method based on bat algorithm
CN106950933B (en) Quality conformance control method and device, computer storage medium
CN110455902A (en) A kind of method that the more standard specimens of environment measuring are quickly calibrated

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20171031

Termination date: 20200428

CF01 Termination of patent right due to non-payment of annual fee