CN1869640A - Method for investigating fatigue crack expansion - Google Patents

Method for investigating fatigue crack expansion Download PDF

Info

Publication number
CN1869640A
CN1869640A CN 200610051602 CN200610051602A CN1869640A CN 1869640 A CN1869640 A CN 1869640A CN 200610051602 CN200610051602 CN 200610051602 CN 200610051602 A CN200610051602 A CN 200610051602A CN 1869640 A CN1869640 A CN 1869640A
Authority
CN
China
Prior art keywords
curve
crack
test
angular frequency
fatigue
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN 200610051602
Other languages
Chinese (zh)
Other versions
CN100439896C (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CNB2006100516025A priority Critical patent/CN100439896C/en
Publication of CN1869640A publication Critical patent/CN1869640A/en
Application granted granted Critical
Publication of CN100439896C publication Critical patent/CN100439896C/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

The invention discloses a fatigue crack extension testing method, making a sweeping test in the fatigue test course and obtaining the variation condition of vibration response of the test system varies with test piece crack extension on the forced vibration condition; using a numeric method to make sweeping data fitting and obtaining variation curve of resonant frequency with number of circulating times; after the test piece is broken, according to actually observed crack size and shape, prefabricating crack of a CAD model and making finite model analysis and obtaining a curve of resonant angular frequency of the system against test piece size, using the curve and able to obtain dynamic extension rule of test piece crack and then obtaining its curve of variation rate; and finally combining with calculation of stress strength factor range Delta K, and able to obtain a curve of crack extension rate.

Description

A kind of method of investigating fatigue crack expansion
Technical field
The present invention relates to a kind of crack growth test method that is applied to the component of machine torture test.
Background technology
The fatigure failure of mechanical component is its a kind of principal mode that inefficacy takes place.Conclude from various statistical data, its ratio is about about 80%~95%.Therefore problem of fatigue damage is major fields of Mechanical Reliability research.Fatigure failure is a very complicated problems, and the factor that influences fatigue strength is numerous, and most influence factor also can't provide strict Physical Mechanism at present and explain, also is difficult to provide more careful quantitative description from mathematics.
People from different viewing angles, have formed the engineering approximation describing method of some fatigure failures gradually in a large amount of production practices.Wherein, the rate curve of crack Propagation has been obtained good effect in the fatigue reliability estimation of a lot of actual engineering problems.Fatigue crack growth rate (da/dN-Δ K) is not only an important indicator of exosyndrome material fatigue crack-resistant toughness, also is to the traditional fatigue strength design and the important supplement of content of the test.Along with the progressively development of Theory of Fracture Mechanics and perfect, this parameter has obtained application more and more widely in Fatigue Design.
One of main task of fatigue crack propagation test is exactly to measure crack growth rate.The mensuration of these data then depends on the kinetic measurement of crack size a, and promptly crack size a loads the dynamic changing curve of times N with test.This measurement accurately whether test findings is had significant impact, be the key issue of fatigue crack propagation test.
The dynamic measurement method of existing crack size is a lot, and commonly used in the laboratory have portable measuring microscope direct-reading methods, potential method, flexibility method, eddy-current method, acoustic-emission and based on the mechanics method of experimental modal analysis etc.; And on-the-spot parts crackle is measured commonly used have eddy-current method, fluorescence method and magnaflux etc.For different subjects, the process of the test in the varying environment, various test methods in various degree application all arranged.But these methods all have a same limitation, promptly often can only be applied to material test or simple test specimen.For numerous actual parts, as engine crankshaft etc., because the complicacy of its structure and crack position and form is special, these methods all are difficult to implement.
Summary of the invention
The objective of the invention is for overcoming the limitation of original fatigue crack test method, a kind of method of novel investigating fatigue crack expansion is provided, makes that some have labyrinth and special process or crack position, form is special and parts that be not easy to observe can directly carry out the test of crack growth rate.
The method of investigating fatigue crack expansion of the present invention may further comprise the steps:
1) test specimen is installed on the load maintainer of forming by master, passive swing arm and eccentric wheel vibrator;
2) torture test is carried out in loading, and every is carried out frequency sweep test at regular intervals one time, and the resulting exciting angular frequency of frequency sweep-vibratory response acceleration information is carried out match by the prototype function that formula (1) is established by least square method:
α = A + B ω 2 1 - C ω 2 , ω ∈ ( 0.3 p , 0.85 p ) - - - ( 1 )
In the formula: ω is the exciting angular frequency; α is the vibratory response acceleration; P is the resonant angular frequency of vibrational system; A, B and C are undetermined coefficient, obtain the dynamic change N-p curve of resonant angular frequency p with test cycle index N after the match;
3) finish in torture test, after test specimen ruptures, obtain the morphological feature and the size a of actual crack by actual observation fracture apperance fatigue striation;
4) it is prefabricated to utilize the actual observation data to carry out cad model, and the resonant angular frequency p of each crack size a drag is carried out finite element analysis, obtains the a-p relation curve;
5) utilize the corresponding relation of a-p curve and N-p curve, draw the a-N curve; According to the a-N curve, increase progressively quadratic polynomial or secant method calculates the da/dN-N curve with 7, the stress intensity factor range Δ K in conjunction with Finite Element Method calculates obtains the crack growth rate curve, i.e. da/dN-Δ K curve.
Beneficial effect of the present invention is:
Flow process by this test method, it is after test specimen ruptures, and under the situation that the particular location, form and the size characteristic that make crackle can be surveyed fully, the monitor message in the control test process utilizes the finite element simulation technology that the spreading rate of crackle is calculated accurately again.Therefore, the inventive method has been given full play to the technological merit of resonance loaded type fatigue tester and ingenious utilization the " posteriority " data, can tolerate farthest that the complicacy and the polytrope of position, crackle form and crack propagation mode appear in crackle in the torture test process.In fact under the enough meticulous situation of finite element model, its measuring accuracy is subjected to the influence of these factors hardly.The inventive method on its technology implementation route, embodied some with test method in the past different peculiar thinking.Make that some have labyrinth and special process or crack position, form is special and parts that be not easy to observe can directly carry out the test of crack growth rate.
The inventive method is with respect to traditional straight survey method, electrical measuring method and based on the crack growth test methods such as mechanics method of experimental modal analysis.It is not subjected to the restriction of sample dimensions, structure, technology characteristics and crack position, form, so long as utilize the torture test of resonance loading principle, can carry out crack growth test to any parts.
Description of drawings
Fig. 1 is the synoptic diagram of resonant mode torture test load maintainer;
Fig. 2 is that frequency sweep method crankshaft cracks spreading rate is measured flow process;
Fig. 3 is the frequency sweep illustrated example;
Fig. 4 is the prefabricated synoptic diagram of finite element analysis model crackle, a) the resonance mechanism three-dimensional perspective behind the precrack wherein, b) be a) in the enlarged drawing of semielliptic surface crack.
Embodiment
The method of investigating fatigue crack expansion of the present invention may further comprise the steps:
1) with test specimen, bent axle 1 for example shown in Figure 1 is installed on the resonant mode load maintainer of being made up of master arm 2, passive arm 3 and eccentric wheel vibrator 4.
2) load and to carry out torture test, and every at regular intervals, as 50,000 times or 100,000 times, carry out frequency sweep test one time.Frequency sweep test is in order to obtain the resonant frequency of vibrational system under the forced vibration condition, and its method is by applying the excitation of different frequency, and writes down under each frequency excitation the response acceleration value of system.Promptly obtain data sequence { ω i, α i| I=1,2,3..., the sequence of actual measurement is drawn on the figure, can obtain frequency sweep figure, Fig. 3 has provided some examples of frequency sweep figure constantly.Along with the increase of test cycle index, the frequency sweep curve moves to left and shows the decline of its resonant frequency.
The resulting exciting angular frequency of frequency sweep-vibratory response acceleration information is carried out match by the prototype function that formula (1) is established by least square method:
α = A + B ω 2 1 - C ω 2 , ω ∈ ( 0.3 p , 0.85 p ) - - - ( 1 )
In the formula: ω is the exciting angular frequency; α is the vibratory response acceleration; P is the resonant angular frequency of vibrational system; A, B and C are undetermined coefficient, use least square fitting and obtain after coefficient A, B and the C, and the calculating formula of resonant angular frequency p is:
p = 1 C
Obtain the dynamic change situation of resonant angular frequency p after the match, i.e. the N-p curve with test cycle index N;
3) finish in torture test, after test specimen ruptures, obtain the morphological feature and the size a of actual crack by actual observation fracture apperance fatigue striation.
4) it is prefabricated to utilize the morphological feature of crackle of actual observation and size a that cad model is carried out crackle.Concrete grammar is: on the CAD of former resonator system solid model, survey the crackle form factually in the actual crack position and size a draws an ellipse, with this ellipse is laminar entity of profile stretching, then round-corner transition being done at its edge handles so that its three-dimensional configuration with actual crack is similar to, at last, with the flawless cad model for being operated entity, with this thin slice it is done entity boolean reducing and finish the prefabricated of crackle, Fig. 4 is the partial enlarged drawing of prefabricated three-dimensional perspective of cad model crackle and crackle.Resonant angular frequency p to each crack size a drag carries out finite element analysis, obtains the a-p relation curve.
5) utilize the corresponding relation of a-p curve and N-p curve, draw the a-N curve; According to the a-N curve, increase progressively quadratic polynomial or secant method calculates change rate curve with 7, i.e. the da/dN-N curve.Stress intensity factor can be the 4th) calculate with Finite Element Method on the prefabricated model of CAD crackle finished of step, computing application 1/4 displacement method, the unit adopts 1/4 unusual yuan of 15 node.Stress intensity factor range Δ K in conjunction with calculating obtains the crack growth rate curve, i.e. da/dN-Δ K curve.

Claims (1)

1. the method for an investigating fatigue crack expansion may further comprise the steps:
1) test specimen is installed on the load maintainer of forming by master, passive swing arm and eccentric wheel vibrator;
2) torture test is carried out in loading, and every is carried out frequency sweep test at regular intervals one time, and the resulting exciting angular frequency of frequency sweep-vibratory response acceleration information is carried out match by the prototype function that formula (1) is established by least square method:
α = A + Bω 2 1 - Cω 2 , ω∈(0.3p,0.85p) (1)
In the formula: ω is the exciting angular frequency; α is the vibratory response acceleration; P is the resonant angular frequency of vibrational system; A, B and C are undetermined coefficient, obtain the dynamic change N-p curve of resonant angular frequency p with test cycle index N after the match;
3) finish in torture test, after test specimen ruptures, obtain the morphological feature and the size a of actual crack by actual observation fracture apperance fatigue striation;
4) it is prefabricated to utilize the actual observation data to carry out cad model, and the resonant angular frequency p of each crack size a drag is carried out finite element analysis, obtains the a-p relation curve;
5) utilize the corresponding relation of a-p curve and N-p curve, draw the a-N curve; According to the a-N curve, increase progressively quadratic polynomial or secant method calculates the da/dN-N curve with 7, the stress intensity factor range Δ K in conjunction with Finite Element Method calculates obtains the crack growth rate curve, i.e. da/dN-Δ K curve.
CNB2006100516025A 2006-05-24 2006-05-24 Method for investigating fatigue crack expansion Expired - Fee Related CN100439896C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2006100516025A CN100439896C (en) 2006-05-24 2006-05-24 Method for investigating fatigue crack expansion

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2006100516025A CN100439896C (en) 2006-05-24 2006-05-24 Method for investigating fatigue crack expansion

Publications (2)

Publication Number Publication Date
CN1869640A true CN1869640A (en) 2006-11-29
CN100439896C CN100439896C (en) 2008-12-03

Family

ID=37443362

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2006100516025A Expired - Fee Related CN100439896C (en) 2006-05-24 2006-05-24 Method for investigating fatigue crack expansion

Country Status (1)

Country Link
CN (1) CN100439896C (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101975704A (en) * 2010-10-08 2011-02-16 中国船舶重工集团公司第七二五研究所 Method for predicting vibration fatigue performance of viscoelastic material
CN102253087A (en) * 2011-06-22 2011-11-23 南京航空航天大学 Device and method for automatically measuring fatigue crack propagation velocity
CN102426137A (en) * 2011-08-31 2012-04-25 中国船舶重工集团公司第七二五研究所 Testing and measuring method capable of axially loading low cycle fatigue crack initiation
CN103868786A (en) * 2014-03-18 2014-06-18 东北大学 Method for predicting fatigue crack propagation rule
CN106338445A (en) * 2016-08-25 2017-01-18 浙江工业大学 Fatigue crack propagation test image online acquisition method based on stroboscopic illumination
CN107884290A (en) * 2017-10-16 2018-04-06 北京理工大学 A kind of fretting fatigue cracks extension life-span prediction method under consideration effect of attrition
CN108827804A (en) * 2018-07-12 2018-11-16 浙江工业大学 A kind of resonant mode fatigue tester dynamic load error online compensation method
CN109541037A (en) * 2018-11-21 2019-03-29 扬州大学 A kind of loading machine structure trunk crackle screening technique based on compromise with mechanism of competition
CN110376062A (en) * 2019-07-22 2019-10-25 中国航发沈阳发动机研究所 A kind of resonant fatigue test crackle pre-setting method
CN112665810A (en) * 2020-12-28 2021-04-16 湖北亿咖通科技有限公司 Method and system for determining chip vibration falling, storage medium and electronic equipment
CN117894407A (en) * 2023-12-25 2024-04-16 西南交通大学 Material fatigue damage tolerance determination method based on near field dynamics theory

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102645385A (en) * 2012-04-12 2012-08-22 西安石油大学 Method for detecting influence of different parameters on material fatigue crack expansion performance

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5071492A (en) * 1989-12-19 1991-12-10 Parker Research Inc. Method for reducing the fatigue crack growth rate of cracks in the aluminum alloy fuselage skin of an aircraft structure
JPH08278241A (en) * 1995-04-07 1996-10-22 Toshiba Corp Crack-development-characteristic evaluation apparatus
JPH095220A (en) * 1995-06-21 1997-01-10 Ishikawajima Harima Heavy Ind Co Ltd Fatigue crack introducing method for disc-like testpiece
CN1598528A (en) * 2004-09-06 2005-03-23 浙江大学 Crankshaft bend fatigue test system and method
CN2729691Y (en) * 2004-09-06 2005-09-28 浙江大学 Crankshaft bending fatigue experiment system

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101975704B (en) * 2010-10-08 2013-03-06 中国船舶重工集团公司第七二五研究所 Method for predicting vibration fatigue performance of viscoelastic material
CN101975704A (en) * 2010-10-08 2011-02-16 中国船舶重工集团公司第七二五研究所 Method for predicting vibration fatigue performance of viscoelastic material
CN102253087A (en) * 2011-06-22 2011-11-23 南京航空航天大学 Device and method for automatically measuring fatigue crack propagation velocity
CN102426137A (en) * 2011-08-31 2012-04-25 中国船舶重工集团公司第七二五研究所 Testing and measuring method capable of axially loading low cycle fatigue crack initiation
CN102426137B (en) * 2011-08-31 2013-06-05 中国船舶重工集团公司第七二五研究所 Testing and measuring method capable of axially loading low cycle fatigue crack initiation
CN103868786A (en) * 2014-03-18 2014-06-18 东北大学 Method for predicting fatigue crack propagation rule
CN106338445B (en) * 2016-08-25 2019-07-05 浙江工业大学 Fatigue crack propagation test image online acquisition method based on stroboscopic illumination
CN106338445A (en) * 2016-08-25 2017-01-18 浙江工业大学 Fatigue crack propagation test image online acquisition method based on stroboscopic illumination
CN107884290B (en) * 2017-10-16 2019-11-19 北京理工大学 It is a kind of to consider that the fretting fatigue cracks under the influence of abrasion extend life-span prediction method
CN107884290A (en) * 2017-10-16 2018-04-06 北京理工大学 A kind of fretting fatigue cracks extension life-span prediction method under consideration effect of attrition
CN108827804A (en) * 2018-07-12 2018-11-16 浙江工业大学 A kind of resonant mode fatigue tester dynamic load error online compensation method
CN108827804B (en) * 2018-07-12 2021-04-06 浙江工业大学 Dynamic load error online compensation method for resonant fatigue testing machine
CN109541037A (en) * 2018-11-21 2019-03-29 扬州大学 A kind of loading machine structure trunk crackle screening technique based on compromise with mechanism of competition
CN109541037B (en) * 2018-11-21 2021-01-01 扬州大学 Loader structure main crack screening method based on compromise and competition mechanism
CN110376062A (en) * 2019-07-22 2019-10-25 中国航发沈阳发动机研究所 A kind of resonant fatigue test crackle pre-setting method
CN112665810A (en) * 2020-12-28 2021-04-16 湖北亿咖通科技有限公司 Method and system for determining chip vibration falling, storage medium and electronic equipment
CN112665810B (en) * 2020-12-28 2023-05-30 亿咖通(湖北)技术有限公司 Method and system for determining vibration shedding of chip, storage medium and electronic equipment
CN117894407A (en) * 2023-12-25 2024-04-16 西南交通大学 Material fatigue damage tolerance determination method based on near field dynamics theory

Also Published As

Publication number Publication date
CN100439896C (en) 2008-12-03

Similar Documents

Publication Publication Date Title
CN100439896C (en) Method for investigating fatigue crack expansion
Yang et al. Experimental study on peak shear strength criterion for rock joints
Andreaus et al. Cracked beam identification by numerically analysing the nonlinear behaviour of the harmonically forced response
Xiang et al. Identification of crack in a rotor system based on wavelet finite element method
Jassim et al. A review on the vibration analysis for a damage occurrence of a cantilever beam
Ostachowicz Damage detection of structures using spectral finite element method
Wu et al. Experimental and numerical studies on the evolution of shear behaviour and damage of natural discontinuities at the interface between different rock types
Yang et al. An automatic finite element modelling for deformation analysis of composite structures
Alfano et al. Determining the elastic constants of isotropic materials by modal vibration testing of rectangular thin plates
CN112326787A (en) Beam bridge identification method based on multipoint rapid static acquisition of exclusive test car
De Falco et al. Sensitivity analysis of frequency-based tie-rod axial load evaluation methods
Żółtowski et al. The use of modal analysis in the evaluation of welded steel structures.
Dos Santos et al. An overview of experimental strain-based modal analysis methods
Dilena et al. Crack identification in rods and beams under uncertain boundary conditions
Lepik Exploring vibrations of cracked beams by the Haar wavelet method.
Zhang et al. Natural frequency response evaluation for RC beams affected by steel corrosion using acceleration sensors
Wan et al. Three-dimensional micromechanical finite element analysis on gauge length dependency of the dynamic modulus of asphalt mixtures
Reddy et al. Digital image correlation for structural health monitoring–a review
Huang et al. Numerical simulation of internal stress in pavement concrete under rolling fatigue load
Shokouhifard et al. Inverse dynamic analysis of an inclined FGM beam due to moving load for estimating the mass of moving load based on a CGM
El-Kafrawy Crack detection by modal analysis in 3D beams based on FEM
Ali et al. Constrained piezoelectric thin film for sensing of subsurface cracks
Juhászová et al. Evaluation of fatigue crack growth rates in an IPE beam made of AISI 304 under various stress ratios
Yu et al. Measurement of the M-integral for a hole in an aluminum plate or strip
An et al. Identification of beam crack using the dynamic response of a moving spring-mass unit

Legal Events

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

Granted publication date: 20081203

Termination date: 20110524