CN107391799A - A kind of low-cycle fatigue data analysing method based on matlab programmings - Google Patents
A kind of low-cycle fatigue data analysing method based on matlab programmings Download PDFInfo
- Publication number
- CN107391799A CN107391799A CN201710481640.2A CN201710481640A CN107391799A CN 107391799 A CN107391799 A CN 107391799A CN 201710481640 A CN201710481640 A CN 201710481640A CN 107391799 A CN107391799 A CN 107391799A
- Authority
- CN
- China
- Prior art keywords
- column matrix
- fatigue
- matrix
- stress
- strain
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
A kind of low-cycle fatigue data analysing method based on matlab programmings, the present invention relates to the low-cycle fatigue data analysing method programmed based on matlab.The invention aims to solve existing method using when need to be manually entered experiment value, easily there is input error, influence analysis of experiments result, data-handling efficiency is low；And the problem of obtained image can not carry out secondary editor, and the variate-value of processing procedure is not directly perceived, and matched curve numerical value can not be exported directly.Process is：Generate total data matrix D ata；Calculate the bound numerical value of Dfe, Ef, K corresponding to 95% confidential interval；Result is stored in result and write in the entitled result of file form document；Font, font size in setting image；Draw three curves and scatterplot on fatigue life strain pattern；Plot functions are called to draw plastic strain steady cycle stress curve map.The present invention is used for low-cycle fatigue data analysis field.
Description
Technical field
The present invention relates to the low-cycle fatigue data analysing method programmed based on matlab.
Background technology
Under alternate load effect, fatigue rupture will occur under the stress level far below yield strength for material.
The failure mode of more than 80% component of machine is all fatigue rupture.So the material to needing to bear alternate load, such as turns
Sub- steel, blade, spring supporting etc. must carry out fatigue test.Material low cycle fatigue property is evaluated and component's life prediction is dependent on
In material or the low cycle fatigue test of component, therefore material low cycle fatigue test tested number is very big, and test result analysis treat
Need the test data amount that inputs very big in journey.
Traditional low-cycle fatigue data analysing method need to be manually entered typically using basic language, C language etc., during use
Experiment value, and microstrain needs after input decimal point four, and life value may need to input 5 bit digitals, therefore, easily occur
Input error, influences analysis of experiments result, and data-handling efficiency is low；And the image obtained by conventional programming techniques can not carry out two
Secondary editor, the variate-value of processing procedure is not directly perceived, and matched curve numerical value can not be exported directly.
The content of the invention
The invention aims to solve existing method using when need to be manually entered experiment value, it is wrong easily to there is input
By mistake, analysis of experiments result is influenceed, data-handling efficiency is low；And obtained image can not carry out secondary editor, the change of processing procedure
The problem of value is not directly perceived, and matched curve numerical value can not be exported directly, and propose a kind of low-cycle fatigue based on matlab programmings
Data analysing method.
It is a kind of based on matlab programming low-cycle fatigue data analysing method detailed process be：
Step 1: prediction testing of materials log file is obtained according to MTS testing machines during low cycle fatigue test, to obtaining
Prediction testing of materials log file in the dependent variable, fatigue life, maximum stress, minimum stress, the modulus of elasticity data that record
Information is copied directly to generate total data matrix D ata in matlab program human windows；
Overall strain amount is total data matrix D ata first row Et；Fatigue life is total data matrix D ata secondary series
Nf；The 3rd row max_stress that maximum stress is total data matrix D ata；Minimum stress is the 4 of total data matrix D ata
Arrange min_stress；The 5th row E that modulus of elasticity is total data matrix D ata；
Step 2: distinguished with elastic strain column matrix Ee, plastic strain column matrix Ep and stress average column matrix stress
To store elastic strain array, plastic strain array and stress average array,
Wherein stress=(max_stress-min_stress)/2, Ee=stress/E, Ep=Et-Ee；2 times of fatigues
Life-span column matrix Nf2=2*Nf；
Elastic strain column matrix Ee, plastic strain column matrix Ep, 2 times of fatigue life column matrix Nf2 and stress average are arranged
Each element seeks the result of calculation after denary logarithm in matrix stress, is separately stored in column matrix Ee1, column matrix
In Ep1, column matrix Nf1 and column matrix stress1；Plastic strain is exactly that overall strain subtracts elastic strain；
Step 3: calling the built-in polyfit functions in matlab, as abscissa and indulged using column matrix Nf1 and Ee1 respectively
Coordinate, column matrix Nf1 and Ep1 are abscissa and ordinate, Ep1 and stress1 is that abscissa and ordinate do 3 groups of Linear Quasis
Close, numerical value fatigue strength coefficient Dfe and tired will be converted into using the parameter that column matrix Nf1 and Ee1 is abscissa and ordinate output
Labor intensity index b；Fatigue ductile coefficient Ef will be converted into by the parameter that abscissa and ordinate export of column vector Nf1 and Ep1
With fatigue ductility index c；Intensity of circulation COEFFICIENT K will be converted into by the parameter that abscissa and ordinate export of Ep1 and stress1
With cyclic strain hardening exponent n；
Step 4: the fatigue strength coefficient Dfe obtained according to step 3, fatigue ductile coefficient Ef and intensity of circulation COEFFICIENT K,
Calculate fatigue strength coefficient Dfe corresponding to 95% confidential interval, fatigue ductile coefficient Ef, the bound number of intensity of circulation COEFFICIENT K
Value；
Fatigue strength coefficient Dfe bound numerical value is up_Dfe and down_Dfe；
Fatigue ductile coefficient Ef bound numerical value is up_Ef and down_Ef；
The bound numerical value of intensity of circulation COEFFICIENT K is up_K and down_K；
Step 5: setting empty matrix of consequence result, the numerical value obtained in step 3 and step 4 is stored in what is set
In empty matrix of consequence result and write in the entitled result of file form document；
Step 6: the newly-built drawing window of figure orders in matlab is called, font, font size in setting image；
Step 7: draw three curves on fatigue life-strain pattern；
Step 8: draw the scatterplot on fatigue life-strain pattern；
Step 9: a figure orders newly-built drawing window again in matlab is called, with column matrix Ep and column matrix
Stress is respectively the abscissa value and ordinate value put on curve, and calling plot functions to draw plastic strain-stable circulation should
Force curve figure.
Beneficial effects of the present invention are：
The present invention proposes a kind of low-cycle fatigue data analysing method based on matlab programmings, to realize to large-scale data
Reproducible input, improve data-handling efficiency, can obtain directly carrying out the image model of secondary editor in result,
The variate-value of processing procedure is directly perceived, can export matched curve numerical value, and low-cycle fatigue analysis result and program drawing editor can be with
Disposably completed using matlab, ensure that the accuracy of data input, data-handling efficiency improves 2 times.Solve existing
Method using when need to be manually entered experiment value, easily there is input error, influence analysis of experiments result；And obtained image can not
Secondary editor is carried out, the problem of variate-value of processing procedure is not directly perceived, and matched curve numerical value can not be exported directly, improves low week
The efficiency and reliability of fatigue data analysis.
Using matlab Programming with Pascal Language efficiency highs, C language etc. is avoided in numerical fitting and calculates confidential interval bound
When using recursive language complexity, program debugging is very convenient；The data in the initial report of duplication are deposited using matlab
Store up as original matrix, greatly reduce the workload of artificial input data；Data are saved as each column vector, data processed result
It is more directly perceived, convenient export；Matlab programs calling figure from after drawing as tool box, utilizes the powerful picture editting's work(in tool box
It can handle and directly handle curve, avoid the trouble that export data are edited again to drawing softwares such as origin/excel.
Brief description of the drawings
Fig. 1 is flow chart of the present invention；
Fig. 2 be step 7 step 8 of the present invention draw fatigue life-strain pattern on three curves and fatigue life-should
The scatter diagram become on image, strain are strain, and unit mm/mm, cycle are the cycle, and ε t are overall strain amount, and ε e should for elasticity
Become column matrix, ε p are plastic strain column matrix；
Fig. 3 is that step 9 of the present invention calls plot functions to draw plastic strain-steady cycle stress curve map；
Fig. 4 is the design sketch after the adjustment of matlab Curve Tools case in the embodiment of the present invention.
Embodiment
Embodiment one：Illustrate present embodiment with reference to Fig. 1, present embodiment it is a kind of based on matlab programmings
Low-cycle fatigue data analysing method detailed process is：
Step 1: prediction testing of materials log file is obtained according to MTS testing machines during low cycle fatigue test, to obtaining
Prediction testing of materials log file in the dependent variable, fatigue life, maximum stress, minimum stress, the modulus of elasticity data that record
Information is copied directly to generate total data matrix D ata in matlab program human windows when this program performs；
Overall strain amount is total data matrix D ata first row Et；Fatigue life is total data matrix D ata secondary series
Nf；The 3rd row max_stress that maximum stress is total data matrix D ata；Minimum stress is the 4 of total data matrix D ata
Arrange min_stress；The 5th row E that modulus of elasticity is total data matrix D ata；
Step 2: distinguished with elastic strain column matrix Ee, plastic strain column matrix Ep and stress average column matrix stress
To store elastic strain array, plastic strain array and stress average array,
Wherein stress=(max_stress-min_stress)/2, Ee=stress/E, Ep=Et-Ee；2 times of fatigues
Life-span column matrix Nf2=2*Nf；
Elastic strain column matrix Ee, plastic strain column matrix Ep, 2 times of fatigue life column matrix Nf2 and stress average are arranged
Each element seeks the result of calculation after denary logarithm in matrix stress, is separately stored in column matrix Ee1, column matrix
In Ep1, column matrix Nf1 and column matrix stress1；Plastic strain is exactly that overall strain subtracts elastic strain；
Step 3: calling the built-in polyfit functions in matlab, as abscissa and indulged using column matrix Nf1 and Ee1 respectively
Coordinate, column matrix Nf1 and Ep1 are abscissa and ordinate, Ep1 and stress1 is that abscissa and ordinate do 3 groups of Linear Quasis
Close, numerical value fatigue strength coefficient Dfe and tired will be converted into using the parameter that column matrix Nf1 and Ee1 is abscissa and ordinate output
Labor intensity index b；Fatigue ductile coefficient Ef will be converted into by the parameter that abscissa and ordinate export of column vector Nf1 and Ep1
With fatigue ductility index c；Intensity of circulation COEFFICIENT K will be converted into by the parameter that abscissa and ordinate export of Ep1 and stress1
With cyclic strain hardening exponent n；
Step 4: the fatigue strength coefficient Dfe obtained according to step 3, fatigue ductile coefficient Ef and intensity of circulation COEFFICIENT K,
Calculate fatigue strength coefficient Dfe corresponding to 95% confidential interval, fatigue ductile coefficient Ef, the bound number of intensity of circulation COEFFICIENT K
Value；
Fatigue strength coefficient Dfe bound numerical value is up_Dfe and down_Dfe；
Fatigue ductile coefficient Ef bound numerical value is up_Ef and down_Ef；
The bound numerical value of intensity of circulation COEFFICIENT K is up_K and down_K；
Step 5: setting empty matrix of consequence result, the numerical value obtained in step 3 and step 4 is stored in what is set
In empty matrix of consequence result and write in the entitled result of file form document；
Step 6: the newly-built drawing window of figure orders in matlab is called, font, font size in setting image；
Step 7: draw three curves on fatigue life-strain pattern；
Step 8: draw the scatterplot on fatigue life-strain pattern；
Step 9: a figure orders newly-built drawing window again in matlab is called, with column matrix Ep and column matrix
Stress is respectively the abscissa value and ordinate value put on curve, and calling plot functions to draw plastic strain-stable circulation should
Force curve figure.
With matlab Curve Tool casees, the selection of line color line style, label editing, transverse and longitudinal coordinate setting, font are completed
The image processing requirements such as editor, derived curve data.It is corresponding under each total variable except the numerical value that command windows directly export
Elastic strain matrix Ee, plastic strain matrix Ep and figure in each curve on point transverse and longitudinal coordinate value can be by matlab
Directly export carries out analysis contrast in each variable of workplace windows.Each total variable refers to predicting testing of materials record data
The dependent variable of middle record.
Embodiment two：Present embodiment is unlike embodiment one：Column matrix in the step 2
Ee1 is to seek each element in elastic strain column matrix Ee the column matrix after denary logarithm；Column matrix Ep1 is to plasticity
Each element seeks the column matrix after denary logarithm in strain column matrix Ep；Column matrix Nf1 is to 2 times of fatigue life row squares
Each element seeks the column matrix after denary logarithm in battle array Nf2；Column matrix stress1 is to stress average column matrix
Each element seeks the column matrix after denary logarithm in stress.
Other steps and parameter are identical with embodiment one.
Embodiment three：Present embodiment is unlike embodiment one or two：Painted in the step 7
Three curves on fatigue life processed-strain pattern；Detailed process is：
Each point on three curves is determined in image according to 2 times of fatigue life column matrix Nf2 maximum and minimum value respectively
Abscissa value, wherein fatigue life-elastic strain curve, fatigue life-plastic strain curve, fatigue life-overall strain
The ordinate value of each point is respectively Dfe* (Nf2) on curveb, Ef* (Nf2)c, Dfe* (Nf2)b+Ef*(Nf2)c；Call plot letters
Number is drawn, and sets the color of curve, linear, line width；
Dfe is fatigue strength coefficient, and Ef is fatigue ductile coefficient, and Nf2 is 2 times of fatigue life column matrix, and b is fatigue strength
Index, c are fatigue ductility index, and * is multiplication sign.
Other steps and parameter are identical with embodiment one or two.
Embodiment four：Unlike one of present embodiment and embodiment one to three：The step 8
The middle scatterplot drawn on fatigue life-strain pattern；Detailed process is：
Keep drawing interface motionless, draw three groups of scatterplots, each point is in length and breadth wherein in fatigue life-elastic strain scatterplot group
Coordinate value is respectively derived from column matrix Nf2 and column matrix Ee, the coordinate value in length and breadth of each point in fatigue life-plastic strain scatterplot group
It is respectively derived from column matrix Nf2 and column matrix Ep, the difference of the coordinate value in length and breadth source of each point in fatigue life-overall strain scatterplot group
In column matrix Nf2 and column matrix Et, the shape of scatterplot is set, and each group scatterplot is set by the graph color set in step 7
Color.
Other steps and parameter are identical with one of embodiment one to three.
Beneficial effects of the present invention are verified using following examples：
Embodiment one：
A kind of low-cycle fatigue data analysing method based on matlab programmings of the present embodiment is specifically to enter according to following steps
Capable：
One group of low cycle fatigue test data instance only with certain steel under a certain experiment condition：
1st, prediction testing of materials log file is obtained according to MTS testing machines during low cycle fatigue test, to the pre- of acquisition
Measure and monitor the growth of standing timber the overall strain amount that records in material test data sheet file, fatigue life, maximum stress, minimum stress, modulus of elasticity data letter
Breath is copied directly to generate total data matrix D ata in matlab program human windows when this program performs：Overall strain amount is
Total data matrix D ata first row Et；Fatigue life is total data matrix D ata secondary series Nf；Maximum stress is total data
Matrix D ata the 3rd row max_stress；The 4th row min_stress that minimum stress is total data matrix D ata；Springform
Measure and arrange E for the 5th of total data matrix D ata；
Table 1MTS testing machine record datas
This experimental data obtains after processing
2nd, stored respectively with elastic strain column matrix Ee, plastic strain column matrix Ep and stress average column matrix stress
Elastic strain array, plastic strain array and stress average array, stress=(max_stress-min_stress)/2, Ee
=stress/E, Ep=Et-Ee,；Column matrix Nf2=2*Nf；To elastic strain column matrix Ee, plastic strain row Ep, 2 times of life-spans
Each element seeks the result calculated after denary logarithm in column matrix Nf2 and stress average column matrix stress, stores up respectively
Exist in column matrix Ee1, column matrix Ep1, column matrix Nf1 and column matrix stress1；Plastic strain is exactly that overall strain subtracts elasticity
Strain；This experimental data obtains Ee1=[- 2.64029 after processing；-2.64354；-2.65647；-2.65025；-
2.66763；-2.66238；-2.66186；-2.67038；-2.68187；-2.6794；-2.68121；-2.68411；-
2.70329；-2.71003；-2.74538；- 2.76909] Ep1=[- 2.49341；-2.4911；-2.55372；-2.55869；-
2.62887；-2.63373；-2.73955；-2.72957；-2.8478；-2.85146；-3.03787；-3.03133；-
3.28416；-3.25941；-3.69317；- 3.52548] Nf1=[3.164353；3.203033；3.360215；3.317227；
3.486005；3.476252；3.610234；3.625724；3.798789；3.784902；3.909556；3.919392；
4.30311；4.341909]
Stress1=[2.566001；2.563297；2.548144；2.554138；2.547799；2.542041；
2.539741；2.534851；2.524344；2.526972；2.521177；2.516734；2.492677；2.491523]
3rd, the built-in polyfit functions in matlab are called, respectively using column matrix Nf1 and Ee1 as abscissa and ordinate,
Column matrix Nf1 and Ep1 are abscissa and ordinate, Ep1 and stress1 is that abscissa and ordinate do 3 groups of linear fits, with row
Vectorial Nf1 and Ee1 are that the parameter that abscissa and ordinate export can be converted into numerical value fatigue strength coefficient Dfe and fatigue strength refers to
Number b；Fatigue ductile coefficient Ef can be converted into as the parameter that abscissa and ordinate export using row vector Nf1 and Ep1 and fatigue is prolonged
Sex index c；Intensity of circulation COEFFICIENT K can be converted into as the parameter that abscissa and ordinate export using Ep1 and stress1 and circulation should
Become hardenability value n；This experimental data obtains Dfe=0.0036865 after processing；B=-0.0651045；Ef=0.565657；c
=-0.69478；K=607.2929；N=0.08961；
4th, according to step 3, fatigue strength coefficient, fatigue ductile coefficient, intensity of circulation corresponding to 95% confidential interval are calculated
The bound numerical value of coefficient；Bound the numerical value up_Dfe and down_Dfe of fatigue strength coefficient, fatigue ductile coefficient bound
Numerical value up_Ef and down_Ef, intensity of circulation coefficient bound numerical value up_K and down_K；This experimental data obtains after processing
To up_Dfe=0.003823, up_Ef=0.7398, up_K=627.90, down_Dfe=0.003554；Down_Ef=
0.4325, down_K=587.36；
The 5th, empty matrix of consequence is set, the numerical value obtained in step 3 and step 4 is stored in the empty matrix of consequence set
In result and write in the entitled result of file form document, the matrix r esult=that this Data Processing in Experiment obtains
[0.0036865-0.06510450.565657-0.69478607.29290.08961
0.003823-0.06510450.7398-0.69478627.900.08961
0.003554-0.06510450.4325-0.69478587.360.08961]；
6th, Step 6: calling the newly-built drawing window of figure orders in matlab, the preliminary font set in image,
Font size；Reference axis font size is arranged to 14 in this Data Processing in Experiment；
7th, three curves on fatigue life-strain pattern are drawn.Distinguished according to column matrix Nf2 maximum and minimum value
To determine the abscissa value of each point on three curves in first image, wherein fatigue life-elastic strain curve, tired longevity
The ordinate value of each point is respectively Dfe* (Nf2) ^b, Ef* (Nf2) ^ in life-plastic strain curve, labor life-span-overall strain curve
C, Dfe* (Nf2) ^b+Ef* (Nf2) ^c；Plot function plottings are called, the color of curve, linear, line width are set；Taken in this program
The each point of abscissa fatigue life point is at intervals of 1；
8th, the scatterplot on fatigue life-strain pattern is drawn.Keep drawing interface motionless, three groups of scatterplots are drawn, wherein tired
Each point is vertical in labor life-span-elastic strain scatterplot group, fatigue life-plastic strain scatterplot group, fatigue life-overall strain scatterplot group
Abscissa value is respectively derived from column matrix Nf2 and column matrix Ee, column matrix Nf2 and column matrix Ep, column matrix Nf2 and column matrix
Et, the shape of scatterplot is set, and the color of each group scatterplot is set by the color in step 7；Obtained after this Data Processing in Experiment
Image results such as Fig. 2；
9th, a figure orders newly-built drawing window again in matlab is called, is divided with column matrix Ep and column matrix stress
The abscissa value and ordinate value that Wei do not put on curve, call plot functions to draw plastic strain-steady cycle stress curve map；
The so treated obtained result such as Fig. 3 of this experimental result；
10th, with matlab Curve Tool casees, the selection of line color line style, label editing, transverse and longitudinal coordinate setting, word are completed
The image processing requirements such as body editor, derived curve data.It is right under each total variable except the numerical value that command windows directly export
The point transverse and longitudinal coordinate value on each curve in elastic strain matrix Ee, plastic strain matrix Ep and the figure answered can be by workplace
The directly export of window variable carries out analysis contrast.This experimental data obtains Ee=[0.002289 after processing；0.002272；
0.002205；0.002237435；0.002150；0.002175；0.002178；0.002136；0.002080；0.002092；
0.002083；0.002069；0.00198；0.001950；0.001797；0.001702] Ep=[0.003211；0.003228；
0.002794；0.002763；0.002350；0.002324；0.001822；0.001864；0.001420；0.001408；
0.000917；0.000930；0.000520；0.000550；0.000203；0.000298], derived fatigue life-overall strain,
Fatigue life-elastic strain, the horizontal value of fatigue life-plastic strain curve are x=[1460；1461；1462；1463；1464；
1465；1466……17838；17839；17840；17841；17842；17843], ordinate value is respectively y33=
[0.005875；0.005873；0.005871；0.005870；0.005868；0.005866……0.002014862；
0.002014858；0.002014854；0.002014850；0.002014847；0.002014843], y11=
[0.002294062；0.002293960；0.002293858；0.002293756；0.002293654；0.002293552……
0.001776323；0.001776321；0.00177632；0.001776318；0.001776317；0.001776315], y22=
[0.003580981；0.003579277；0.003577576；0.003575877；0.003574180；0.003572485……
0.000233632；0.000233630；0.000233628；0.000233626；0.000233623；0.000233621],
The effect such as Fig. 4 of this figure after the adjustment of matlab Curve Tools case：
Attached configuration processor is as follows：
The present invention can also have other various embodiments, in the case of without departing substantially from spirit of the invention and its essence, this area
Technical staff works as can make various corresponding changes and deformation according to the present invention, but these corresponding changes and deformation should all belong to
The protection domain of appended claims of the invention.
Claims (4)
- A kind of 1. low-cycle fatigue data analysing method based on matlab programmings, it is characterised in that：Methods described detailed process is：Step 1: prediction testing of materials log file is obtained according to MTS testing machines during low cycle fatigue test, to the pre- of acquisition Measure and monitor the growth of standing timber material test data sheet file in record dependent variable, fatigue life, maximum stress, minimum stress, modulus of elasticity data message It is copied directly to generate total data matrix D ata in matlab program human windows；Overall strain amount is total data matrix D ata first row Et；Fatigue life is total data matrix D ata secondary series Nf；Most The 3rd row max_stress that big stress is total data matrix D ata；The 4th row min_ that minimum stress is total data matrix D ata stress；The 5th row E that modulus of elasticity is total data matrix D ata；Step 2: stored up respectively with elastic strain column matrix Ee, plastic strain column matrix Ep and stress average column matrix stress Elastic strain array, plastic strain array and stress average array are deposited,Wherein stress=(max_stress-min_stress)/2, Ee=stress/E, Ep=Et-Ee；2 times of fatigue lives Column matrix Nf2=2*Nf；To elastic strain column matrix Ee, plastic strain column matrix Ep, 2 times of fatigue life column matrix Nf2 and stress average column matrix Each element seeks the result of calculation after denary logarithm in stress, is separately stored in column matrix Ee1, column matrix Ep1, row In matrix N f1 and column matrix stress1；Step 3: the built-in polyfit functions in matlab are called, respectively using column matrix Nf1 and Ee1 as abscissa and vertical seat Mark, column matrix Nf1 and Ep1 are abscissa and ordinate, Ep1 and stress1 is that abscissa and ordinate do 3 groups of linear fits, Numerical value fatigue strength coefficient Dfe will be converted into as the parameter that abscissa and ordinate export using column matrix Nf1 and Ee1 and fatigue is strong Spend index b；Fatigue ductile coefficient Ef and tired will be converted into by the parameter that abscissa and ordinate export of column vector Nf1 and Ep1 Labor ductility index c；Intensity of circulation COEFFICIENT K will be converted into as the parameter that abscissa and ordinate export using Ep1 and stress1 and followed Ring strain hardening exponent n；Step 4: the fatigue strength coefficient Dfe obtained according to step 3, fatigue ductile coefficient Ef and intensity of circulation COEFFICIENT K, are calculated Fatigue strength coefficient Dfe, fatigue ductile coefficient Ef, the bound numerical value of intensity of circulation COEFFICIENT K corresponding to 95% confidential interval；Fatigue strength coefficient Dfe bound numerical value is up_Dfe and down_Dfe；Fatigue ductile coefficient Ef bound numerical value is up_Ef and down_Ef；The bound numerical value of intensity of circulation COEFFICIENT K is up_K and down_K；Step 5: setting empty matrix of consequence result, the numerical value obtained in step 3 and step 4 is stored in the empty of setting In matrix of consequence result and write in the entitled result of file form document；Step 6: the newly-built drawing window of figure orders in matlab is called, font, font size in setting image；Step 7: draw three curves on fatigue life-strain pattern；Step 8: draw the scatterplot on fatigue life-strain pattern；Step 9: a figure orders newly-built drawing window again in matlab is called, with column matrix Ep and column matrix stress The abscissa value and ordinate value put respectively on curve, plot functions are called to draw plastic strain-steady cycle stress curve Figure.
- A kind of 2. low-cycle fatigue data analysing method based on matlab programmings according to claim 1, it is characterised in that：Institute It is to seek each element in elastic strain column matrix Ee the column matrix after denary logarithm to state column matrix Ee1 in step 2； Column matrix Ep1 is to seek each element in plastic strain column matrix Ep the column matrix after denary logarithm；Column matrix Nf1 is Column matrix after denary logarithm is sought each element in 2 times of fatigue life column matrix Nf2；Column matrix stress1 is pair Each element seeks the column matrix after denary logarithm in stress average column matrix stress.
- A kind of 3. low-cycle fatigue data analysing method based on matlab programmings according to claim 2, it is characterised in that：Institute State three curves drawn in step 7 on fatigue life-strain pattern；Detailed process is：Determine the horizontal stroke of each point on three curves in image respectively according to 2 times of fatigue life column matrix Nf2 maximum and minimum value Coordinate value, wherein fatigue life-elastic strain curve, fatigue life-plastic strain curve, fatigue life-overall strain curve The ordinate value of upper each point is respectively Dfe* (Nf2)b, Ef* (Nf2)c, Dfe* (Nf2)b+Ef*(Nf2)c；Plot functions are called to paint Figure, the color of curve, linear, line width are set；Dfe is fatigue strength coefficient, and Ef is fatigue ductile coefficient, and Nf2 is 2 times of fatigue life column matrix, and b refers to for fatigue strength Number, c is fatigue ductility index, and * is multiplication sign.
- A kind of 4. low-cycle fatigue data analysing method based on matlab programmings according to claim 3, it is characterised in that：Institute State the scatterplot drawn in step 8 on fatigue life-strain pattern；Detailed process is：Keep drawing interface motionless, draw three groups of scatterplots, the coordinate in length and breadth of each point wherein in fatigue life-elastic strain scatterplot group Value is respectively derived from column matrix Nf2 and column matrix Ee, the difference of coordinate value in length and breadth of each point in fatigue life-plastic strain scatterplot group From column matrix Nf2 and column matrix Ep, the coordinate value in length and breadth of each point is respectively derived from row in fatigue life-overall strain scatterplot group Matrix N f2 and column matrix Et, the shape of scatterplot is set, and the face of each group scatterplot is set by the graph color set in step 7 Color.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710481640.2A CN107391799A (en) | 2017-06-22 | 2017-06-22 | A kind of low-cycle fatigue data analysing method based on matlab programmings |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710481640.2A CN107391799A (en) | 2017-06-22 | 2017-06-22 | A kind of low-cycle fatigue data analysing method based on matlab programmings |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107391799A true CN107391799A (en) | 2017-11-24 |
Family
ID=60332554
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710481640.2A Pending CN107391799A (en) | 2017-06-22 | 2017-06-22 | A kind of low-cycle fatigue data analysing method based on matlab programmings |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107391799A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108073397A (en) * | 2017-12-21 | 2018-05-25 | 沈阳鼓风机集团自动控制***工程有限公司 | Automatic programming method and electronic equipment |
| CN110909412A (en) * | 2019-12-09 | 2020-03-24 | 江南大学 | MATLAB-based force displacement curve main characteristic point batch processing method |
| CN111721647A (en) * | 2020-06-24 | 2020-09-29 | 四川大学 | Low-cycle fatigue test data processing and internal stress evaluation method |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102967512A (en) * | 2012-10-25 | 2013-03-13 | 中南大学 | Low cycle fatigue life prediction method based on asymmetrical cyclic stress control load |
| CN103344515A (en) * | 2013-07-05 | 2013-10-09 | 北京航空航天大学 | Damage calculation method for low-cycle fatigue and high-strength impact coupling based on local stress strain method |
-
2017
- 2017-06-22 CN CN201710481640.2A patent/CN107391799A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102967512A (en) * | 2012-10-25 | 2013-03-13 | 中南大学 | Low cycle fatigue life prediction method based on asymmetrical cyclic stress control load |
| CN103344515A (en) * | 2013-07-05 | 2013-10-09 | 北京航空航天大学 | Damage calculation method for low-cycle fatigue and high-strength impact coupling based on local stress strain method |
Non-Patent Citations (3)
| Title |
|---|
| 刘长虹 等: ""对低周疲劳寿命有重要影响试验参数的确定方法"", 《机械工程材料》 * |
| 李武镇: ""基于MATLAB疲劳工具箱的开发"", 《中国优秀硕士学位论文全文数据库-基础科学辑》 * |
| 黄利军 等: ""Ti-1023钛合金的低周疲劳行为"", 《稀有金属材料与工程》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108073397A (en) * | 2017-12-21 | 2018-05-25 | 沈阳鼓风机集团自动控制***工程有限公司 | Automatic programming method and electronic equipment |
| CN110909412A (en) * | 2019-12-09 | 2020-03-24 | 江南大学 | MATLAB-based force displacement curve main characteristic point batch processing method |
| CN110909412B (en) * | 2019-12-09 | 2021-06-01 | 江南大学 | MATLAB-based force displacement curve main characteristic point batch processing method |
| CN111721647A (en) * | 2020-06-24 | 2020-09-29 | 四川大学 | Low-cycle fatigue test data processing and internal stress evaluation method |
| CN111721647B (en) * | 2020-06-24 | 2021-12-28 | 四川大学 | Low-cycle fatigue test data processing and internal stress evaluation method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107391799A (en) | A kind of low-cycle fatigue data analysing method based on matlab programmings | |
| CN104236893B (en) | Performance parameter test system and performance parameter test method of hydraulic damper | |
| CN105718641A (en) | Method for establishing CAE two-dimensional grids in batches | |
| CN109163997A (en) | A kind of rock surface strength detection method based on sonograph deep learning | |
| CN109815855B (en) | Electronic equipment automatic test method and system based on machine learning | |
| CN109711707B (en) | Comprehensive state evaluation method for ship power device | |
| JPH01134624A (en) | Data name standardizing system in program | |
| CN110347570B (en) | Analysis method for code automatic generation tool in IDE environment | |
| CN115393645A (en) | Automatic soil classification and naming method and system, storage medium and intelligent terminal | |
| CN106991050B (en) | False positive identification method for reference defect of static test null pointer | |
| CN117266936A (en) | Mining pressure display characteristic monitoring method and device based on support resistance | |
| CN108703762A (en) | A kind of Mental Workload based on mental arithmetic and working memory etc. induces technology | |
| MX2022013313A (en) | Selecting chromatography parameters for manufacturing therapeutic proteins. | |
| CN111145056A (en) | Character cultivation method and device based on data mining and machine learning | |
| CN116012456A (en) | Visual positioning method for electric power inspection robot | |
| JP7184093B2 (en) | Skilled Index Providing Device, Skilled Index Providing Method, and Program | |
| CN106685727A (en) | Screening design scheme of performance test data | |
| CN109165108B (en) | Failure data reduction method and test method for software reliability accelerated test | |
| CN112732549A (en) | Test program classification method based on cluster analysis | |
| Reimers | Simulation of mixed mode fatigue crack growth | |
| Strubergs et al. | Use of StanForD2010 data for determination of effect of harvester operator periodic training on productivity and fuel economy | |
| CN105510963B (en) | A kind of method for analyzing seismic prospecting auxiliary data file | |
| CN108596167A (en) | A kind of elongated optical character word system for rapidly identifying and method, model | |
| CN114526127B (en) | Turbine valve flow characteristic test system | |
| CN115170001B (en) | Bridge construction risk assessment method based on improved Kent index method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20171124 |
|
| WD01 | Invention patent application deemed withdrawn after publication |