CN105260574B - A kind of all Multiaxial Fatigue Life Prediction methods of height based on critical surface method Fatigue criteria - Google Patents

A kind of all Multiaxial Fatigue Life Prediction methods of height based on critical surface method Fatigue criteria Download PDF

Info

Publication number
CN105260574B
CN105260574B CN201510780883.7A CN201510780883A CN105260574B CN 105260574 B CN105260574 B CN 105260574B CN 201510780883 A CN201510780883 A CN 201510780883A CN 105260574 B CN105260574 B CN 105260574B
Authority
CN
China
Prior art keywords
stress
fatigue
critical surface
width
loading
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510780883.7A
Other languages
Chinese (zh)
Other versions
CN105260574A (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.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN201510780883.7A priority Critical patent/CN105260574B/en
Publication of CN105260574A publication Critical patent/CN105260574A/en
Application granted granted Critical
Publication of CN105260574B publication Critical patent/CN105260574B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)

Abstract

The all Multiaxial Fatigue Life Prediction methods of height that the present invention provides a kind of based on critical surface method Fatigue criteria, are related to fatigue strength field, steps of the method are:(1) the high all non-proportional loading criterion of construction.(2) situation is loaded by substituting into two kinds of uniaxial loadings and pulsating cyclic, derives that parameter in formula (3) reads high all multiaxis Identical loading test datas, determines critical surface position.(4) by the shear stress width on calculated critical surface, direct stress width and mean normal stress, equivalent stress width is calculated using formula.(5) stress- life is utilized, the high-Cycle Fatigue Life Prediction under multiaxis constant amplitude loading is calculated.The method of proposition carries out the high-Cycle Fatigue Life Prediction estimation under multiaxis constant amplitude loading and obtains preferable prediction effect.

Description

A kind of all Multiaxial Fatigue Life Predictions of height based on critical surface method Fatigue criteria Method
Technical field
The present invention relates to metal non-proportional loading strength fields, refer in particular to the high-cycle fatigue life under multiaxis constant amplitude loading Method.
Background technology
The Fatigue Failures of machinery are one of the main reason for causing in commission machine components and component failure.Example Such as in aircraft, ship, engineering machinery field, machine components and component often bear the effect of complicated alternating load for a long time, most Eventually since fatigue failure causes accident to account for 80 or more the percent of total failare factor.Have at present in single shaft fatigue field and compares Perfect theory simultaneously has accumulated mass data and experience, but practical in commission each major mechanical parts be typically subjected to it is complicated The effect of Multiaxial Proportional and nonproportional cycle loading.Due to multi-axis stress state, the variation of direction of crack propagation, load history Complexity, so traditional single shaft fatigue strength theory has been not suitable with present design requirement, instead non-proportional loading Strength theory, especially high all Multiaxial Fatigue Life Prediction methods, which become, solves component of machine Problem of Failure in practical be on active service Important method.The all Multiaxial Fatigue Life Predictions of traditional height often ignore the influence of mean stress, therefore consider mean stress shadow The all Multiaxial Fatigue Life Prediction methods of loud height are worth research.
Invention content
Present invention aims at the requirements for non-proportional loading Intensity Design, propose a kind of based on critical surface method fatigue rupture The all Multiaxial Fatigue Life Prediction methods of height of criterion.
A kind of all Multiaxial Fatigue Life Prediction methods of height based on critical surface method Fatigue criteria provided by the present invention, Its step is:
Step 1):Thin-wall tube test specimen is critical with face where maximum shear stress by drawing-torsion disproportional Sine-wave loading Face utilizes the fatigue damage parameter on the face:Shear stress width Ca, direct stress width Na, mean normal stress NmOne non-linear group of construction Conjunction form isThe all non-proportional loading criterion of height.
Step 2):α in solution formula, β, λ.By substituting into two kinds of uniaxial loading situations, it is respectively uniaxially stretched tired pole Limit f-1Load and pure torsion fatigue limit t-1Load is releasedλ=t-1.Situation is loaded by pulsating cyclic again, It corrects to obtain using Goodman equationsσ in formulauTo be uniaxially stretched strength degree.To sum up release one Kind is based on the high all non-proportional loading failure criteria of critical surface method:
Step 3):Read high all multiaxis Identical loading test datas.Include load direct stress width σ in every group of test dataa, Shear stress width τa, mean normal stress σm, average shearing stress τm, phase difference δ.
Step 4):Determine critical surface position.In order to consider to load influence of the average shearing stress to fatigue damage, it is believed that crackle Plane germinating and extension where maximum shear stress, i.e., using face where maximum shear stress as critical surface.
Step 5):Pass through the shear stress width C on calculated critical surfacea, direct stress width NaWith mean normal stress Nm, utilize Formula calculates equivalent stress width:
Step 6):Utilize equivalent stress width τeqThe S-N curve for being -1 with stress ratio, calculates multiaxis constant amplitude loading Under high-Cycle Fatigue Life Prediction.
The step 1) and step 2) are constructed containing there are three Damage Parameter (shear stress width Ca, direct stress width Na, averagely just Stress Nm) nonlinear combination form beThe all non-proportional loading criterion of height, this method is logical It crosses this Damage Parameter by common maximum (normal) stress and is decomposed into direct stress width and mean normal stress, consider respectively and fatigue is damaged The Different Effects of injured labour life.
Compared with prior art, the present invention has the advantages that:
The present invention proposes a kind of all Multiaxial Fatigue Life Prediction methods of the height based on critical surface method Fatigue criteria.It should Method is decomposed into direct stress width and mean normal stress by this Damage Parameter by common maximum (normal) stress, considers him respectively For fatigue damage generate Different Effects, construct the Fatigue criteria of a tri-consult volume, and the determination of critical surface Consider the influence of average shearing stress.It is applied widely, and it is based on critical surface law theory, explicit physical meaning, the material being related to Expect that constant is less, is convenient for engineer application.By verification experimental verification, the high cycle fatigue longevity under multiaxis constant amplitude loading is carried out using this method It orders estimation and obtains preferable prediction effect.
Description of the drawings
Fig. 1 is all Multiaxial Fatigue Life Prediction sides of the height based on critical surface method Fatigue criteria that the method for the present invention provides The flow chart of method.
Specific implementation mode
Illustrate the specific implementation mode of the present invention in conjunction with attached drawing and verification experimental verification.
High-cycle fatigue life under the multiaxis constant amplitude loading of the invention to be carried out is utilized a kind of based on critical surface The Fatigue criteria and prediction model of method, and pass through its reliability of verification experimental verification and accuracy.
Step 1):Thin-wall tube test specimen is utilized the oblique section formula in the mechanics of materials by drawing-torsion disproportional Sine-wave loading It can calculate under given loading environment, each stress value on determined critical surface, such as shear stress width Ca, direct stress width Na, averagely just Stress Nm
Step 3):α in solution formula, β, λ.By substituting into two kinds of uniaxial loading situations, it is respectively uniaxially stretched tired pole Limit f-1Load and pure torsion fatigue limit t-1Load is releasedλ=t-1.Situation is loaded by pulsating cyclic again, It corrects to obtain using Goodman equationsσ in formulauTo be uniaxially stretched strength degree.To sum up release one The all non-proportional loading failure criteria of height of the kind based on critical surface method:
Step 4):Read the fatigue data of three groups of materials:7075-T651 aluminium alloys, 30CrMnSiA steel, LY12CZ Aluminium alloy.Using Papadopoulos formula, maximum shear stress C is determinedmaxPlace plane, i.e. critical surface position, and calculate this and face Shear stress width C on interfacea, direct stress width NaWith mean normal stress Nm
Step 5):Pass through the shear stress width C on calculated critical surfacea, direct stress width NaWith mean normal stress Nm, utilize Formula calculates equivalent stress width:
Step 6):Utilize equivalent stress width τeqThe S-N curve for being -1 with stress ratio, calculates multiaxis constant amplitude loading Under high-Cycle Fatigue Life Prediction.Such as the expression formula of the S-N curves of 7075-T651 aluminium alloys is τa=719.08 (Nf)-0.1194, formula Middle NfFor fatigue life.Convert multi-axis stress state to uniaxial stress state, τeq=719.08 (Nf)-0.1194, you can it calculates Fatigue life.
It is pre- in order to verify a kind of height based on critical surface method Fatigue criteria proposed by the present invention all non-proportional loading service life Prediction result obtained by this method is compared by the effect of survey method with the test life obtained by multiaxial experiment, the results showed that, Critical surface method Fatigue criteria based on proposition, Multiaxial Proportional that computational methods through the invention obtain, disproportional fatigue Life prediction value is compared with testing actual life, and the error factor overwhelming majority is within 3 times.This method is based on critical surface method, The Different Effects of direct stress width and mean normal stress to high all Multiaxial Fatigue Damages are considered respectively.The method of proposition carries out multiaxis High-Cycle Fatigue Life Prediction estimation under constant amplitude loading obtains preferable prediction effect.

Claims (2)

1. a kind of all Multiaxial Fatigue Life Prediction methods of height based on critical surface method Fatigue criteria, it is characterised in that:Step It is as follows,
Step 1):Thin-wall tube test specimen is by drawing-torsion disproportional Sine-wave loading, using face where maximum shear stress as critical surface, profit With the fatigue damage parameter on the face:Shear stress width Ca, direct stress width Na, mean normal stress NmConstruct a nonlinear combination shape Formula isThe all non-proportional loading criterion of height;
Step 2):α in solution formula, β, λ;By substituting into two kinds of uniaxial loading situations, it is respectively uniaxially stretched fatigue limit f-1 Load and pure torsion fatigue limit t-1Load is releasedSituation is loaded by pulsating cyclic again, utilizes Goodman Equation is corrected to obtainσ in formulauTo be uniaxially stretched strength degree;It to sum up releases a kind of based on critical surface The high all non-proportional loading failure criteria of method:
Step 3):Read high all multiaxis Identical loading test datas;Include load direct stress width σ in every group of test dataa, cut and answer Power width τa, mean normal stress σm, average shearing stress τm, phase difference δ;
Step 4):Determine critical surface position;Crackle plane germinating and extension where maximum shear stress, i.e., with maximum shear stress institute It is critical surface in face, it is contemplated that influence of the average shearing stress to fatigue damage;
Step 5):Pass through the shear stress width C on calculated critical surfacea, direct stress width NaWith mean normal stress Nm, utilize formula Calculate equivalent stress width:
Step 6):Utilize equivalent stress width τeqThe S-N curve for being -1 with stress ratio, calculates under multiaxis constant amplitude loading High-Cycle Fatigue Life Prediction.
2. all Multiaxial Fatigue Life Prediction sides of a kind of height based on critical surface method Fatigue criteria according to claim 1 Method, it is characterised in that:The step 1) and step 2) are constructed is containing the nonlinear combination form there are three Damage ParameterThe all non-proportional loading criterion of height, three Damage Parameters are respectively shear stress width Ca, just answer Power width Na, mean normal stress Nm;This method is decomposed into direct stress width peace by this Damage Parameter by common maximum (normal) stress Equal direct stress considers the Different Effects generated to fatigue damage respectively.
CN201510780883.7A 2015-11-15 2015-11-15 A kind of all Multiaxial Fatigue Life Prediction methods of height based on critical surface method Fatigue criteria Active CN105260574B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510780883.7A CN105260574B (en) 2015-11-15 2015-11-15 A kind of all Multiaxial Fatigue Life Prediction methods of height based on critical surface method Fatigue criteria

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510780883.7A CN105260574B (en) 2015-11-15 2015-11-15 A kind of all Multiaxial Fatigue Life Prediction methods of height based on critical surface method Fatigue criteria

Publications (2)

Publication Number Publication Date
CN105260574A CN105260574A (en) 2016-01-20
CN105260574B true CN105260574B (en) 2018-08-07

Family

ID=55100263

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510780883.7A Active CN105260574B (en) 2015-11-15 2015-11-15 A kind of all Multiaxial Fatigue Life Prediction methods of height based on critical surface method Fatigue criteria

Country Status (1)

Country Link
CN (1) CN105260574B (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106441765A (en) * 2016-11-30 2017-02-22 上海航天精密机械研究所 Setting method for triaxial vibration test conditions according to maximum stress equivalence
CN107066677A (en) * 2017-01-19 2017-08-18 北京工业大学 A kind of determination method in random combined tension and torsion lower critical face
CN107506535B (en) * 2017-08-07 2020-09-01 电子科技大学 Multi-axial fatigue life prediction method based on critical strain damage parameters
CN107748817B (en) * 2017-10-19 2021-04-09 北京工业大学 High-temperature multi-axis constitutive relation determination method considering non-proportional additional reinforcement
CN109033709B (en) * 2018-08-30 2020-03-31 电子科技大学 Component fatigue life evaluation method based on nonlinear fatigue damage accumulation theory
CN109115479A (en) * 2018-08-30 2019-01-01 电子科技大学 A kind of turbine wheel shaft Multiaxial Fatigue Life Prediction method based on critical surface
CN109459955B (en) * 2018-11-21 2021-07-20 华南理工大学 FPGA-based multi-axis servo motor position reading and decoding system and method
CN110018049A (en) * 2019-04-24 2019-07-16 长沙理工大学 A kind of asphalt Fatigue Life Prediction method under Simple stress condition
CN110334405B (en) * 2019-06-11 2020-12-25 南京航空航天大学 High-temperature multi-axis low-cycle fatigue life prediction method based on Chaboche structure and Lemailre damage model
CN111624116B (en) * 2020-05-09 2021-09-21 清华大学 Fatigue life prediction method and device based on weighted average maximum shear stress plane
CN114297893B (en) * 2021-12-28 2023-03-28 北京工业大学 Multi-axial fatigue failure life assessment method for welding structure

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103604688A (en) * 2013-12-01 2014-02-26 北京航空航天大学 Prediction method for multi-axial high-cycle fatigue life of plastic metal material based on critical plane approach
CN104392130A (en) * 2014-11-21 2015-03-04 南京衍达软件科技有限公司 Method for determining multi-axis fatigue most damage load direction and application thereof

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5721227B2 (en) * 2011-08-25 2015-05-20 株式会社Ihi Multiaxial fatigue life evaluation method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103604688A (en) * 2013-12-01 2014-02-26 北京航空航天大学 Prediction method for multi-axial high-cycle fatigue life of plastic metal material based on critical plane approach
CN104392130A (en) * 2014-11-21 2015-03-04 南京衍达软件科技有限公司 Method for determining multi-axis fatigue most damage load direction and application thereof

Also Published As

Publication number Publication date
CN105260574A (en) 2016-01-20

Similar Documents

Publication Publication Date Title
CN105260574B (en) A kind of all Multiaxial Fatigue Life Prediction methods of height based on critical surface method Fatigue criteria
CN104699976B (en) A kind of metal material multiaxis high cycle fatigue failure prediction method influenceed comprising mean stress
Lee et al. Fatigue failure of a composite wind turbine blade at its root end
Noroozi et al. Prediction of fatigue crack growth under constant amplitude loading and a single overload based on elasto-plastic crack tip stresses and strains
Dekker et al. A cohesive XFEM model for simulating fatigue crack growth under mixed‐mode loading and overloading
Bang et al. A modification of UniGrow 2‐parameter driving force model for short fatigue crack growth
CN110274826A (en) A kind of hard metal material multiaxis high cycle fatigue failure prediction method based on single shaft fatigue S-N curve
Demir et al. Evaluation of mixed mode-I/II criteria for fatigue crack propagation using experiments and modeling
CN107977516B (en) It is a kind of to consider that the Notched specimen A LOCAL STRESS-STRAIN of multiaxial loading disproportional degree determines method
WO2021227923A1 (en) Notched part fatigue life prediction method and prediction device
CN109635385A (en) A kind of part life prediction technique comprehensively considering Fatigue Strength Effect factor
CN111090957B (en) Dangerous point stress-strain calculation method for high-temperature structure
Kumar et al. Design of an aircraft wing structure for static analysis and fatigue life prediction
Sander et al. Fatigue crack growth under variable amplitude loading Part I: experimental investigations
Montesano et al. Critical stiffness damage envelopes for multidirectional laminated structures under multiaxial loading conditions
Bruyneel et al. Damage modeling of laminated composites: validation of the intra-laminar damage law of SAMCEF at the coupon level for UD plies
Shahani et al. Damage tolerance approach for analyzing a helicopter main rotor blade
Newman Jr et al. Validation of the two-parameter fracture criterion using finite-element analyses with the critical CTOA fracture criterion
Boljanović et al. Computational fatigue analysis of the pin-loaded lug with quarter-elliptical corner crack
CN110020485B (en) Method for analyzing inherent characteristics of suspended thin-wall column shell based on bolt connection
Pawliczek et al. Influence of the irregularity coefficient of loading on calculated fatigue life
CN107066677A (en) A kind of determination method in random combined tension and torsion lower critical face
Rizov Elastic-plastic fracture of functionally graded circular shafts in torsion
Shlyannikov et al. Nonlinear fracture resistance parameters for elements of aviation structures under biaxial loading
Nozhnitsky et al. Numerical simulation of spin testing for turbo machine disks using energy-based fracture criteria

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