CN117709171A - High cycle fatigue failure numerical simulation method and system - Google Patents
High cycle fatigue failure numerical simulation method and system Download PDFInfo
- Publication number
- CN117709171A CN117709171A CN202410167867.XA CN202410167867A CN117709171A CN 117709171 A CN117709171 A CN 117709171A CN 202410167867 A CN202410167867 A CN 202410167867A CN 117709171 A CN117709171 A CN 117709171A
- Authority
- CN
- China
- Prior art keywords
- stage
- model
- numerical simulation
- target component
- point
- 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
- 238000000034 method Methods 0.000 title claims abstract description 73
- 238000004088 simulation Methods 0.000 title claims abstract description 55
- 239000000463 material Substances 0.000 claims abstract description 70
- 230000006378 damage Effects 0.000 claims abstract description 32
- 238000004364 calculation method Methods 0.000 claims abstract description 30
- 125000004122 cyclic group Chemical group 0.000 claims abstract description 24
- 230000008569 process Effects 0.000 claims description 29
- 238000010276 construction Methods 0.000 claims description 27
- 238000010899 nucleation Methods 0.000 claims description 24
- 239000007787 solid Substances 0.000 claims description 18
- 239000000126 substance Substances 0.000 claims description 17
- 239000011343 solid material Substances 0.000 claims description 12
- 238000006073 displacement reaction Methods 0.000 claims description 9
- 239000002245 particle Substances 0.000 claims description 7
- 208000027418 Wounds and injury Diseases 0.000 claims description 6
- 208000014674 injury Diseases 0.000 claims description 6
- 239000002243 precursor Substances 0.000 claims description 6
- 230000000977 initiatory effect Effects 0.000 claims description 5
- 239000013013 elastic material Substances 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 4
- 239000002184 metal Substances 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 10
- 230000009471 action Effects 0.000 description 4
- 238000009661 fatigue test Methods 0.000 description 4
- 238000009434 installation Methods 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- 239000007769 metal material Substances 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 230000002269 spontaneous effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a high-cycle fatigue failure numerical simulation method, which comprises the steps of obtaining parameter information of a target component; establishing a three-dimensional geometric model of the target component and discretizing; constructing a fracture criterion and corresponding material parameters of a target component under high-cycle fatigue failure; and (3) carrying out cyclic load calculation under high-cycle fatigue fracture on the target component and completing numerical simulation of the target component under high-cycle fatigue fracture. The invention also discloses a system for realizing the high cycle fatigue failure numerical simulation method. The invention establishes a complete high-cycle fatigue damage numerical simulation method under near-field dynamics, can simulate the fatigue damage of most metal components under high-cycle load, and has high reliability, good accuracy and good effect.
Description
Technical Field
The invention belongs to the field of data signal simulation, and particularly relates to a high-cycle fatigue failure numerical simulation method and system.
Background
Near field dynamics (PD) is a mechanical theory based on non-local action ideas set forth in the professor sizing in 2000; the theory describes discrete object particle motion by adopting a space integral equation so as to solve the limitation of classical media mechanics in researching the problems of crack initiation, crack propagation and crack bifurcation of solid material prediction due to continuity assumption and local contact principle. In near field dynamics theory, interactions are generated between substance points into which an object is discrete through bonds, and the substance points are only connected with bonds (functions of interactions) between other substance points in the near field range with the near field radius d (as d tends to be 0, the interactions tend to be localized, so classical elastic medium mechanics can be regarded as the limit case when near field dynamics theory approaches 0; when the object is deformed, the key is deformed.
Near field dynamics can be divided into: bond-based near field dynamics (the force state between two material points is the same as the direction of the bond, and the magnitude and the direction are equal to each other), normal state-based near field dynamics (the force state between two material points is the same as the direction of the bond, and the magnitude and the direction are different from each other), and unconventional state-based near field dynamics (the force state between two material points is any direction and the magnitude and the direction are different from each other). And simulating failure of the microscopic lower key according to the unified key breaking criteria (such as key breaking criteria based on critical elongation and key breaking criteria based on critical energy density) by definition, so as to further obtain spontaneous initiation and expansion processes of the macroscopic lower crack.
Although near field dynamics solves the problems in many fracture mechanics fields, numerical simulation of three-dimensional cracks of a component under the action of high Zhou Hezai in engineering is still remained in theoretical perfection and experimental demonstration stages at present, and a numerical simulation method for fatigue life of the component under the action of high Zhou Hezai is still lacking.
Disclosure of Invention
The invention aims to provide a high-cycle fatigue failure numerical simulation method with high reliability, good accuracy and good effect.
The second object of the invention is to provide a system for realizing the high cycle fatigue failure numerical simulation method.
The high cycle fatigue failure numerical simulation method provided by the invention comprises the following steps:
s1, acquiring parameter information of a target component;
s2, establishing a three-dimensional geometric model of the target component based on the geometric parameters and the spatial topological relation of the target component according to the parameter information acquired in the step S1;
s3, discretizing the three-dimensional geometric model constructed in the step S2 based on the voxelization idea, so as to convert the three-dimensional geometric model into equidistant substance points in a three-dimensional space;
s4, determining a fracture criterion and corresponding material parameters of a target component under high-cycle fatigue failure based on the linear near-field dynamics solid material model;
and S5, calculating the cyclic load of the target component under the high-cycle fatigue failure based on the material parameters determined in the step S4 and the discrete model obtained in the step S3, so as to complete the numerical simulation of the target component under the high-cycle fatigue failure.
The step S2 specifically comprises the following steps:
according to the parameter information obtained in the step S1, based on the geometric parameters and the space topological relation of the target component, a three-dimensional geometric model of the target component is established by modeling software; the modeling software comprises SolidWorks and AutoCAD; the model file format is the file format of the importing component supported by the finite element analysis software.
The step S3 specifically comprises the following steps:
importing a three-dimensional geometric model into finite element analysis software, and creating a three-dimensional cube entity capable of completely wrapping the imported part; and carrying out grid division on the created three-dimensional cube entity, adjusting and deleting the grid division result, and finally dispersing the three-dimensional geometric model into equidistant substance points in the three-dimensional space.
The step S3 specifically comprises the following steps:
importing the three-dimensional geometric model of the target component obtained in the step S2 into finite element analysis software as a component;
in the preprocessing part of finite element analysis software, a three-dimensional cube entity capable of completely wrapping the leading-in part is newly built; the three-dimensional cube entity has the length ofWidth of->High->Wherein->Is a positive integer which is used for the preparation of the high-voltage power supply,the distance between adjacent material points is set;
to be used forFor the seed distribution of the global size, adopting a linear hexahedral unit to divide the three-dimensional solid rectangular component into grids by a structural network dividing method;
assembling a three-dimensional geometric model of the target component and the three-dimensional solid rectangular component subjected to grid division: the space position of the three-dimensional solid rectangular component is adjusted, so that the three-dimensional solid rectangular component safely wraps the three-dimensional geometric model constructed by the target; then traversing all hexahedral units in the three-dimensional solid rectangular component, and selecting a hexahedral unit if all nodes of the hexahedral unit are positioned in the three-dimensional geometric model of the target member;
creating a new network component, hiding the selected hexahedral unit in the view, deleting all units in the view, and extracting node space information in the grid model to obtain a voxelized target component particle model;
and finally, deleting the three-dimensional geometric model and the three-dimensional solid rectangular part of the target component, creating a node set for applying constraint and load on the specific interval node according to the required component boundary condition, and writing an input file containing the discretization information of the target construction.
The step S4 specifically comprises the following steps:
determining that the calculation model is a linear near-field dynamics solid material model;
defining the extreme difference of the bond stretching amount under the single cyclic load based on the maximum value and the minimum value of the bond stretching amount under the single cyclic load;
determining a key breaking criterion of a stage I and a stage II in the key breaking process of the fatigue crack;
and determining material parameters of the stage I and the stage II in the bond breaking process of the fatigue crack.
The step S4 specifically comprises the following steps:
determining that the calculation model is a linear near-field dynamics solid material model, wherein corresponding component material parameters comprise density, poisson ratio and bulk modulus;
based on the maximum value and the minimum value of the bond stretching amount under the single cyclic load, the extreme difference of the bond stretching amount under the single cyclic load is defined：
In->Is the maximum value of the bond stretching amount under single cycle load, and +.>，/>Is a deformation precursor particle->And substance point in near field range->Is (are) displacement vector>Is a deformation precursor particle->And substance point in near field range->Is a relative vector of (2); />Is the minimum value of the bond stretching amount under single cycle load, and +.>；RIs stress ratio, and->；
Determining a key breaking criterion of a stage I and a stage II in the key breaking process of the fatigue crack:in->The remaining life of the key is that,Nis the number of load cycles, and->；
The bond breaking process of the fatigue crack comprises a stage I-nucleation stage and a stage II-expansion stage; the phase I-nucleation phase is used for determining the position of fatigue crack initiation, and the following formula is adopted as a corresponding life calculation formula:
in->A first material parameter for a stage I-nucleation stage; />A second material parameter for the stage I-nucleation stage; />The fatigue limit elongation for the stage I-nucleation stage;
the following equation is used as the life calculation equation for the phase II-expansion phase:
in->A first material parameter for a phase II-expansion phase; />A second material parameter that is a phase II-expansion phase;
the following formula was used as an expansion rate formula of fatigue crack:
in the middle oflIs crack length; />Is a geometric parameter; />Is the stretching amount of the core bond;
the following formula is used as a damage function of the object point:
in->Is the point of substance->Injury at time t; />Is the point of substance->Is +.o.f>Is a near field range of (2);is a binary function and->，/>Is the point of substance->Is defined by the volume of (2);
the sign that the object point bond enters the stage II-expansion stage from the stage I-nucleation stage is as follows: object pointIs>The damage of any material point is larger than the set value; mass point->Is>When the damage of any material point is larger than the set value,breaking the corresponding core bond, and reinitializing the residual life of other bonds to 1;
determining material parameters of a stage I-nucleation stage in the bond breaking process of the fatigue crack:
；
strain of material-lifetime N double index table with shape +.>The fitting function of (2) is:
in the middle ofbSlope for the low lifetime portion of the material in the fit function;aconstant coefficients for the fitting function; />Is a strain amplitude; />A first constant that is a fitting function;Ca second constant that is a fitting function;
determining material parameters of a stage II-expansion stage in the breaking process of the fatigue crack:
the fracture mechanics Paris formula of fatigue crack propagation of the linear elastic material is as follows:
in the middle ofcIs the first constant in the Paris formula;Mis the second constant in the Paris formula; />Is a stress intensity factor;
first material parameters of stage II-expansion stageThe calculation formula of (2) is as follows:
in->Is a set arbitrary positive number; />Actual crack growth rate calculated for using the Paris formula; />The expansion rate of the single-side notch sample uniaxial tensile failure numerical simulation is carried out by using the set positive number; />The stress intensity factor is calculated according to the material type and the national standard GB/T6398-2017 fatigue crack propagation method for fatigue test of metallic materials.
The step S5 specifically comprises the following steps:
A. importing the discrete model obtained in the step S3, and inputting model parameters; the model parameters comprise the number, coordinates, material parameters, boundary conditions and near-field radius of each unit node;
B. for the phase I-nucleation phase, the cyclic load calculation was started, and the simulation was performed using a linear time map:
if the current time step is the last step of the load loading, calculating the extremely poor key stretching amountAnd residual lifetime of bond->And breaking the corresponding bond according to the bond breaking criterion;
if the current time step is not the sameThe last step of secondary load loading is to solve the material point in the current time stepNear field range->Another object point->Maximum value of bond stretching amount of inter bond i-j +.>Minimum value of bond stretch->And near field force f, adopting a display center difference method; /> In->Is the point of substance->Density at; />Is the point of substance->Displacement at time step n; />Is the time step;totalis the point of substance->Is>The number of other material points in the container；Is the point of substance->Total near field force experienced, +.>Is the point of substance->Position in space, +.>Is the point of substance->Volume of->Is the point of substance->External force applied at n time steps;
repeating the calculation process of the step B until traversingAll other matter particles in the near field range;
C. calculating the point of matter at the end of the current time stepIs a function of the total near field force, displacement and injury;
repeating the calculation of the step B-step C until all the object points in the component are traversed;
after traversing all the object points, if the damage value of the object point is larger than a set value, entering a stage II-expansion stage; resetting the residual life of the unbroken bond to 1, and recording the data of the current time step;
D. b, entering the next time step, and repeating the steps B-D until the set total time step is reached; finally, the numerical simulation of the target member under high cycle fatigue failure is completed.
The invention also provides a system for realizing the high cycle fatigue failure numerical simulation method, which comprises an information acquisition module, a model construction module, a model discrete module, a condition construction module and a numerical simulation module; the information acquisition module, the model construction module, the model discrete module, the condition construction module and the numerical simulation module are sequentially connected in series; the information acquisition module is used for acquiring parameter information of the target component and uploading the data to the model construction module; the model construction module is used for establishing a three-dimensional geometric model of the target component based on the geometric parameters and the space topological relation of the target component according to the received data, and uploading the data to the model discrete module; the model discrete module is used for dividing the constructed three-dimensional geometric model into grids based on the voxelized idea according to the received data, further dispersing the three-dimensional geometric model into equidistant substance points in the three-dimensional space, and uploading the data to the condition construction module; the condition construction module is used for constructing a fracture criterion and corresponding material parameters of a target component under high-cycle fatigue damage based on the linear near-field dynamic solid material model according to the received data, and uploading the data to the numerical simulation module; the numerical simulation module is used for importing the constructed material parameters into a discrete model according to the received data, and performing cyclic load calculation on the target component under high-cycle fatigue damage, so that the numerical simulation of the target component under the high-cycle fatigue damage is completed.
The high cycle fatigue failure numerical simulation method and the system provided by the invention establish a complete near-field dynamics high cycle fatigue failure numerical simulation method, can simulate fatigue failure of most metal components under high cycle load, and are high in reliability, good in accuracy and good in effect.
Drawings
FIG. 1 is a schematic flow chart of the method of the present invention.
Fig. 2 is a cloud of spring bar stress in a normal installed state of an embodiment of the method of the present invention.
FIG. 3 is a schematic diagram of a broken physical diagram of a spring strip fatigue test according to an embodiment of the method of the invention.
FIG. 4 is a schematic diagram of a simulated spring strip fatigue failure process according to an embodiment of the present invention, wherein (a) is a schematic diagram of a first stage of the fatigue failure process, (b) is a schematic diagram of a second stage of the fatigue failure process, (c) is a schematic diagram of a third stage of the fatigue failure process, and (d) is a schematic diagram of a fourth stage of the fatigue failure process.
FIG. 5 is a schematic diagram of functional modules of the system of the present invention.
Detailed Description
A schematic process flow diagram of the method of the present invention is shown in fig. 1: the high cycle fatigue failure numerical simulation method provided by the invention comprises the following steps:
s1, acquiring parameter information of a target component;
s2, establishing a three-dimensional geometric model of the target component based on the geometric parameters and the spatial topological relation of the target component according to the parameter information acquired in the step S1; the method specifically comprises the following steps:
according to the parameter information obtained in the step S1, based on the geometric parameters and the space topological relation of the target component, a three-dimensional geometric model of the target component is established by modeling software; the modeling software comprises SolidWorks and AutoCAD; the model file format is the file format of the import component (such as sat, igs, iges, stp, step, etc.) supported by finite element analysis software (preferably ABAQUS);
s3, discretizing the three-dimensional geometric model constructed in the step S2 based on the voxelization idea, so as to convert the three-dimensional geometric model into equidistant substance points in a three-dimensional space; the method specifically comprises the following steps:
importing a three-dimensional geometric model into finite element analysis software, and creating a three-dimensional cube entity capable of completely wrapping the imported part; performing grid division on the created three-dimensional cube entity, adjusting and deleting the grid division result, and finally dispersing the three-dimensional geometric model into equidistant substance points in the three-dimensional space;
the specific implementation method comprises the following steps:
importing the three-dimensional geometric model of the target component obtained in the step S2 into finite element analysis software as a component;
in the preprocessing part of finite element analysis software, a three-dimensional cube entity capable of completely wrapping the leading-in part is newly built; the three-dimensional cube entity has the length ofWidth of->High->Wherein->Is a positive integer which is used for the preparation of the high-voltage power supply,the distance between adjacent material points is set;
to be used forFor the seed distribution of the global size, adopting a linear hexahedral unit to divide the three-dimensional solid rectangular component into grids by a structural network dividing method;
assembling a three-dimensional geometric model of the target component and the three-dimensional solid rectangular component subjected to grid division: the space position of the three-dimensional solid rectangular component is adjusted, so that the three-dimensional solid rectangular component safely wraps the three-dimensional geometric model constructed by the target; then traversing all hexahedral units in the three-dimensional solid rectangular component, and selecting a hexahedral unit if all nodes of the hexahedral unit are positioned in the three-dimensional geometric model of the target member;
creating a new network component, hiding the selected hexahedral unit in the view, deleting all units in the view, and extracting node space information in the grid model to obtain a voxelized target component particle model;
finally, deleting the three-dimensional geometric model and the three-dimensional solid rectangular part of the target component, creating a node set for applying constraint and load on a specific interval node according to the required component boundary condition, and writing an input file containing target construction discretization information;
s4, determining a fracture criterion and corresponding material parameters of a target component under high-cycle fatigue failure based on the linear near-field dynamics solid material model; the method specifically comprises the following steps:
determining that the calculation model is a linear near-field dynamics solid material model;
defining the extreme difference of the bond stretching amount under the single cyclic load based on the maximum value and the minimum value of the bond stretching amount under the single cyclic load;
determining a key breaking criterion of a stage I and a stage II in the key breaking process of the fatigue crack;
determining material parameters of a stage I and a stage II in the key breaking process of the fatigue crack;
the specific implementation method comprises the following steps:
because the stress of the material at the high Zhou Hezai does not exceed the yield strength, and the fatigue fracture of the component does not have obvious plastic deformation, namely the brittle fracture, the calculation model is determined to be a linear near-field dynamics solid-state material model, and the corresponding component material parameters comprise density, poisson ratio and bulk modulus;
based on the maximum value and the minimum value of the bond stretching amount under the single cyclic load, the extreme difference of the bond stretching amount under the single cyclic load is defined：
In->Is the maximum value of the bond stretching amount under single cycle load, and +.>，/>Is a deformation precursor particle->And substance point in near field range->Is (are) displacement vector>Is a deformation precursor particle->And substance point in near field range->Is a relative vector of (2); />Is the minimum value of the bond stretching amount under single cycle load, and +.>；RIs stress ratio, and->The method comprises the steps of carrying out a first treatment on the surface of the For the purpose of calculating the efficiency of the process +.>And->，/>Is the minimum of cyclic load, +.>Is the maximum value of the cyclic load;
determining a key breaking criterion of a stage I and a stage II in the key breaking process of the fatigue crack:
in->The remaining life of the key is that,Nis the number of load cycles, and->；
The bond breaking process of the fatigue crack comprises a stage I-nucleation stage and a stage II-expansion stage; the phase I-nucleation phase is used for determining the position of fatigue crack initiation, and the following formula is adopted as a corresponding life calculation formula:
in->A first material parameter for a stage I-nucleation stage; />A second material parameter for the stage I-nucleation stage; />The fatigue limit elongation for the stage I-nucleation stage;
the following equation is used as the life calculation equation for the phase II-expansion phase:
in->A first material parameter for a phase II-expansion phase; />A second material parameter that is a phase II-expansion phase;
the following formula was used as an expansion rate formula of fatigue crack:
in the middle oflIs crack length; />Is a geometric parameter; />Is the stretching amount of the core bond;
the following formula is used as a damage function of the object point:
in->Is the point of substance->Injury at time t; />Is the point of substance->Is +.o.f>Is a near field range of (2);is a binary function and->，/>Is the point of substance->Is defined by the volume of (2);
the sign that the object point bond enters the stage II-expansion stage from the stage I-nucleation stage is as follows: object pointIs>The damage of any material point is larger than the set value (preferably 0.5); mass point->Is>When the damage of any material point is larger than a set value, breaking and destroying the corresponding core bond, and reinitializing the residual life of other bonds to 1;
fitting in a double pair meter of comparative member materialsFunction->And a life calculation formula for determining material parameters of a stage I-nucleation stage in the bond breaking process of the fatigue crack:
；
strain of material-lifetime N double index table with shape +.>The fitting function of (2) is:
in the middle ofbSlope for the low lifetime portion of the material in the fit function;aconstant coefficients for the fitting function; />Is a strain amplitude; />A first constant that is a fitting function;Ca second constant that is a fitting function;
determining material parameters of a stage II-expansion stage in the breaking process of the fatigue crack:
the fracture Paris formula of fatigue crack propagation of the linear elastic material is as follows:
in the middle ofcIs the first constant in the Paris formula;Mis the second constant in the Paris formula; />Is a stress intensity factor;
due toComparison formula->And->Can get->Thus phase II-first material parameter of the expansion phase +.>The calculation formula of (2) is as follows:
in->Is a set arbitrary positive number; />Actual crack growth rate calculated for using the Paris formula; />The expansion rate of the single-side notch sample uniaxial tensile failure numerical simulation is carried out by using the set positive number; />As stress intensity factors, according to the material types, according to the national standard GB/T6398-2017 fatigue crack propagation method for fatigue test of metallic materials;
s5, calculating the cyclic load of the target component under the high-cycle fatigue failure based on the material parameters determined in the step S4 and the discrete model obtained in the step S3, so as to complete the numerical simulation of the target component under the high-cycle fatigue failure; the method specifically comprises the following steps:
A. importing the discrete model obtained in the step S3, and inputting model parameters; the model parameters comprise the number, coordinates, material parameters, boundary conditions and near-field radius of each unit node;
B. for the phase I-nucleation phase, the cyclic load calculation was started, and the simulation was performed using a linear time map:
if the current time step is the last step of the load loading, calculating the extremely poor key stretching amountAnd residual lifetime of bond->And breaking the corresponding bond according to the bond breaking criterion;
if the current time step is not the last step of the load loading, solving the material point in the current time stepNear field range->Another object point->Maximum value of bond stretching amount of inter bond i-j +.>Minimum value of bond stretch->And near field force f, adopting a display center difference method; /> In->Is the point of substance->Density at; />Is the point of substance->Displacement at time step n; />Is the time step;totalis the point of substance->Is>The number of all other material points in the container;is the point of substance->Total near field force experienced, +.>Is the point of substance->Position in space，/>Is the point of substance->Volume of->Is the point of substance->External force applied at n time steps;
repeating the calculation process of the step B until traversingAll other matter particles in the near field range;
C. calculating the point of matter at the end of the current time stepIs a function of the total near field force, displacement and injury;
repeating the calculation of the step B-step C until all the object points in the component are traversed;
after traversing all the object points, if the damage value of the object point is larger than a set value, entering a stage II-expansion stage; resetting the residual life of the unbroken bond to 1, and recording the data of the current time step;
D. b, entering the next time step, and repeating the steps B-D until the set total time step is reached; finally, the numerical simulation of the target member under high cycle fatigue failure is completed.
The effect of the method of the invention is further illustrated in the following in connection with one example:
taking the fatigue damage of the simulated railway WJ-8 fastener elastic strip under the action of high-frequency wheel load as an example in the normal installation state;
firstly, adopting finite element software ABAQUS to simulate stress distribution of the WJ-8 spring strip in a normal installation state. The contact part of the heel end of the spring strip and the backing plate is simulated by adopting a longitudinal grounding spring and a transverse grounding spring with the rigidity of 30kN/mm, and the contact part of the toe end of the spring strip and the insulating plate is vertically fixed and restrained. The pre-tightening force applied by the bolt on the top of the spring strip is 10kN, and the pressure of the top contact area is used for replacing the pre-tightening force. The calculation results are shown in FIG. 2, where the maximum stress occurs at the heel end of the spring strip, about 732MPa, where the stress concentration area of the spring strip.
The fatigue damage of the WJ-8 elastic strip under the vertical load of high Zhou Lieche in the normal installation state is simulated by adopting the near-field dynamics high-cycle fatigue damage system. Vertical train load adopts and uses wheel load amplitudeAnd frequency->The form of the sine function described:
，for static wheel load 63kN +.>Taking the maximum value of 90kN of the vertical load under the running speed of 300km/h of the train. The final failure condition is shown in fig. 3 and 4;
comparing the elastic strip fatigue test result (corresponding to fig. 3) with the near field dynamics simulation result (corresponding to fig. 4), the elastic strip fracture positions are consistent, the elastic strip fracture positions are positioned at the bending position of the outer arm near the heel end and are positioned in the stress concentration area, and the bending-torsion fatigue load fracture characteristics are met.
According to the embodiment, the method has good reliability and accuracy, and can be used for carrying out high-cycle fatigue failure numerical simulation of the metal component.
FIG. 5 is a schematic diagram of functional modules of the system of the present invention: the system for realizing the high cycle fatigue failure numerical simulation method comprises an information acquisition module, a model construction module, a grid division module, a condition construction module and a numerical simulation module; the information acquisition module, the model construction module, the grid division module, the condition construction module and the numerical simulation module are sequentially connected in series; the information acquisition module is used for acquiring parameter information of the target component and uploading the data to the model construction module; the model construction module is used for establishing a three-dimensional geometric model of the target component based on the geometric parameters and the space topological relation of the target component according to the received data, and uploading the data to the grid division module; the grid division module is used for carrying out grid division on the constructed three-dimensional geometric model based on the voxelized idea according to the received data, so that the three-dimensional geometric model is scattered into material points with equal intervals in the three-dimensional space, and the data is uploaded to the condition construction module; the condition construction module is used for constructing a fracture criterion and corresponding material parameters of a target component under high-cycle fatigue damage based on the linear near-field dynamic solid material model according to the received data, and uploading the data to the numerical simulation module; the numerical simulation module is used for importing the constructed material parameters into a discrete model according to the received data, and performing cyclic load calculation on the target component under high-cycle fatigue damage, so that the numerical simulation of the target component under the high-cycle fatigue damage is completed.
Claims (8)
1. The high cycle fatigue failure numerical simulation method is characterized by comprising the following steps of:
s1, acquiring parameter information of a target component;
s2, establishing a three-dimensional geometric model of the target component based on the geometric parameters and the spatial topological relation of the target component according to the parameter information acquired in the step S1;
s3, discretizing the three-dimensional geometric model constructed in the step S2 based on the voxelization idea, so as to convert the three-dimensional geometric model into equidistant substance points in a three-dimensional space;
s4, determining a fracture criterion and corresponding material parameters of a target component under high-cycle fatigue failure based on the linear near-field dynamics solid material model;
and S5, calculating the cyclic load of the target component under the high-cycle fatigue failure based on the material parameters determined in the step S4 and the discrete model obtained in the step S3, so as to complete the numerical simulation of the target component under the high-cycle fatigue failure.
2. The method for simulating high cycle fatigue failure numerical simulation according to claim 1, wherein the step S2 comprises the following steps:
according to the parameter information obtained in the step S1, based on the geometric parameters and the space topological relation of the target component, a three-dimensional geometric model of the target component is established by modeling software; the modeling software comprises SolidWorks and AutoCAD; the model file format is the file format of the importing component supported by the finite element analysis software.
3. The method for simulating high cycle fatigue failure numerical simulation according to claim 2, wherein the step S3 comprises the following steps:
importing a three-dimensional geometric model into finite element analysis software, and creating a three-dimensional cube entity capable of completely wrapping the imported part; and carrying out grid division on the created three-dimensional cube entity, adjusting and deleting the grid division result, and finally dispersing the three-dimensional geometric model into equidistant substance points in the three-dimensional space.
4. The method for simulating high cycle fatigue failure numerical simulation according to claim 3, wherein the step S3 comprises the following steps:
importing the three-dimensional geometric model of the target component obtained in the step S2 into finite element analysis software as a component;
in the preprocessing part of finite element analysis software, a three-dimensional cube entity capable of completely wrapping the leading-in part is newly built; the three-dimensional cube entity has the length ofWidth of->High->Wherein->Is a positive integer>The distance between adjacent material points is set;
to be used forFor the seed distribution of the global size, adopting a linear hexahedral unit to divide the three-dimensional solid rectangular component into grids by a structural network dividing method;
assembling a three-dimensional geometric model of the target component and the three-dimensional solid rectangular component subjected to grid division: the space position of the three-dimensional solid rectangular component is adjusted, so that the three-dimensional solid rectangular component safely wraps the three-dimensional geometric model constructed by the target; then traversing all hexahedral units in the three-dimensional solid rectangular component, and selecting a hexahedral unit if all nodes of the hexahedral unit are positioned in the three-dimensional geometric model of the target member;
creating a new network component, hiding the selected hexahedral unit in the view, deleting all units in the view, and extracting node space information in the grid model to obtain a voxelized target component particle model;
and finally, deleting the three-dimensional geometric model and the three-dimensional solid rectangular part of the target component, creating a node set for applying constraint and load on the specific interval node according to the required component boundary condition, and writing an input file containing the discretization information of the target construction.
5. The method for simulating high cycle fatigue failure numerical simulation according to claim 4, wherein the step S4 comprises the following steps:
determining that the calculation model is a linear near-field dynamics solid material model;
defining the extreme difference of the bond stretching amount under the single cyclic load based on the maximum value and the minimum value of the bond stretching amount under the single cyclic load;
determining a key breaking criterion of a stage I and a stage II in the key breaking process of the fatigue crack;
and determining material parameters of the stage I and the stage II in the bond breaking process of the fatigue crack.
6. The method for simulating high cycle fatigue failure numerical simulation according to claim 5, wherein the step S4 comprises the following steps:
determining that the calculation model is a linear near-field dynamics solid material model, wherein corresponding component material parameters comprise density, poisson ratio and bulk modulus;
based on the maximum value and the minimum value of the bond stretching amount under the single cyclic load, the extreme difference of the bond stretching amount under the single cyclic load is defined：/>In->Is the maximum value of the bond stretching amount under single cycle load, and +.>，/>Is a deformation precursor particle->And substance point in near field range->Is (are) displacement vector>Is a deformation precursor particle->And substance point in near field range->Is a relative vector of (2); />Is the minimum value of the bond stretching amount under single cycle load, and +.>；RIs stress ratio, and->；
Determining a key breaking criterion of a stage I and a stage II in the key breaking process of the fatigue crack:in->The remaining life of the key is that,Nis the number of load cycles, and->；
The bond breaking process of the fatigue crack comprises a stage I-nucleation stage and a stage II-expansion stage; the phase I-nucleation phase is used for determining the position of fatigue crack initiation, and the following formula is adopted as a corresponding life calculation formula:
in->A first material parameter for a stage I-nucleation stage; />A second material parameter for the stage I-nucleation stage; />The fatigue limit elongation for the stage I-nucleation stage;
the following equation is used as the life calculation equation for the phase II-expansion phase:
in->A first material parameter for a phase II-expansion phase; />A second material parameter that is a phase II-expansion phase;
the following formula was used as an expansion rate formula of fatigue crack:
in the middle oflIs crack length; />Is a geometric parameter; />Is the stretching amount of the core bond;
the following formula is used as a damage function of the object point:
in->Is the point of substance->Injury at time t; />Is the point of substance->Is +.o.f>Is a near field range of (2);is a binary function and->，/>Is the point of substance->Is defined by the volume of (2);
the sign that the object point bond enters the stage II-expansion stage from the stage I-nucleation stage is as follows: object pointIs in the near field range of (2)The damage of any material point is larger than the set value; mass point->Is>When the damage of any material point is larger than the set value, breaking the corresponding core bond, and re-initializing the residual life of other bonds1；
Determining material parameters of a stage I-nucleation stage in the bond breaking process of the fatigue crack:
；
strain of material-lifetime N double index table with shape +.>The fitting function of (2) is:
in the middle ofbSlope for the low lifetime portion of the material in the fit function;aconstant coefficients for the fitting function; />Is a strain amplitude; />A first constant that is a fitting function;Ca second constant that is a fitting function;
determining material parameters of a stage II-expansion stage in the breaking process of the fatigue crack:
the fracture mechanics Paris formula of fatigue crack propagation of the linear elastic material is as follows:
in the middle ofcIs the first constant in the Paris formula;Mis the second constant in the Paris formula; />Is a stress intensity factor;
first material parameters of stage II-expansion stageThe calculation formula of (2) is as follows:
in->Is a set arbitrary positive number; />Actual crack growth rate calculated for using the Paris formula; />The expansion rate of the uniaxial tensile failure numerical simulation was performed for the single-side notched specimen using the set positive number.
7. The method for simulating high cycle fatigue failure numerical simulation according to claim 6, wherein the step S5 comprises the following steps:
A. importing the discrete model obtained in the step S3, and inputting model parameters; the model parameters comprise the number, coordinates, material parameters, boundary conditions and near-field radius of each unit node;
B. for the phase I-nucleation phase, the cyclic load calculation was started, and the simulation was performed using a linear time map:
if the current time step is the last step of the load loading, calculating the extremely poor key stretching amountAnd the remaining life of the keyAnd breaking the corresponding part according to the key breaking criterionIs a bond to (a);
if the current time step is not the last step of the load loading, solving the material point in the current time stepNear field range->Another object point->Maximum value of bond stretching amount of inter bond i-j +.>Minimum value of bond stretch->And near field force f, adopting a display center difference method; /> In->Is the point of substance->Density at; />Is the point of substance->Displacement at time step n; />Is the time step;totalis the point of substance->Is>The number of all other material points in the container;is the point of substance->Total near field force experienced, +.>Is the point of substance->Position in space, +.>Is the point of substance->Volume of->Is the point of substance->External force applied at n time steps;
repeating the calculation process of the step B until traversingAll other matter particles in the near field range;
C. calculating the point of matter at the end of the current time stepIs a function of the total near field force, displacement and injury;
repeating the calculation of the step B-step C until all the object points in the component are traversed;
after traversing all the object points, if the damage value of the object point is larger than a set value, entering a stage II-expansion stage; resetting the residual life of the unbroken bond to 1, and recording the data of the current time step;
D. b, entering the next time step, and repeating the steps B-D until the set total time step is reached; finally, the numerical simulation of the target member under high cycle fatigue failure is completed.
8. A system for implementing the high cycle fatigue failure numerical simulation method according to any one of claims 1 to 7, characterized by comprising an information acquisition module, a model construction module, a model discrete module, a condition construction module and a numerical simulation module; the information acquisition module, the model construction module, the model discrete module, the condition construction module and the numerical simulation module are sequentially connected in series; the information acquisition module is used for acquiring parameter information of the target component and uploading the data to the model construction module; the model construction module is used for establishing a three-dimensional geometric model of the target component based on the geometric parameters and the space topological relation of the target component according to the received data, and uploading the data to the model discrete module; the model discrete module is used for dividing the constructed three-dimensional geometric model into grids based on the voxelized idea according to the received data, further dispersing the three-dimensional geometric model into equidistant substance points in the three-dimensional space, and uploading the data to the condition construction module; the condition construction module is used for constructing a fracture criterion and corresponding material parameters of a target component under high-cycle fatigue damage based on the linear near-field dynamic solid material model according to the received data, and uploading the data to the numerical simulation module; the numerical simulation module is used for importing the constructed material parameters into a discrete model according to the received data, and performing cyclic load calculation on the target component under high-cycle fatigue damage, so that the numerical simulation of the target component under the high-cycle fatigue damage is completed.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410167867.XA CN117709171A (en) | 2024-02-06 | 2024-02-06 | High cycle fatigue failure numerical simulation method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410167867.XA CN117709171A (en) | 2024-02-06 | 2024-02-06 | High cycle fatigue failure numerical simulation method and system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117709171A true CN117709171A (en) | 2024-03-15 |
Family
ID=90144715
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410167867.XA Pending CN117709171A (en) | 2024-02-06 | 2024-02-06 | High cycle fatigue failure numerical simulation method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117709171A (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106096134A (en) * | 2016-06-13 | 2016-11-09 | 北京航空航天大学 | A kind of Structural Metallic Fatigue fail-safe analysis based on damage mechanics and Optimization Design |
| WO2022121203A1 (en) * | 2020-12-08 | 2022-06-16 | 江苏科技大学 | Method for calculating spherical shell surface three-dimensional crack propagation fatigue life |
-
2024
- 2024-02-06 CN CN202410167867.XA patent/CN117709171A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106096134A (en) * | 2016-06-13 | 2016-11-09 | 北京航空航天大学 | A kind of Structural Metallic Fatigue fail-safe analysis based on damage mechanics and Optimization Design |
| WO2022121203A1 (en) * | 2020-12-08 | 2022-06-16 | 江苏科技大学 | Method for calculating spherical shell surface three-dimensional crack propagation fatigue life |
Non-Patent Citations (4)
| Title |
|---|
| S. A. SILLING等: "Peridynamic model for fatigue cracking", SANDIA NATIONAL LABORATORIES, 30 September 2014 (2014-09-30), pages 1 - 42 * |
| 张恒等: "近场动力学在断裂力学领域的研究进展", 力学进展, 31 December 2022 (2022-12-31), pages 852 - 873 * |
| 潘雨崧等: "金属材料腐蚀疲劳的近场动力学模拟", 固体力学学报, 23 January 2024 (2024-01-23), pages 1 - 17 * |
| 石宏顺等: "基于近场动力学理论的疲劳断裂分析", 机械强度, 31 July 2018 (2018-07-31), pages 954 - 960 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhang et al. | A review of numerical techniques approaching microstructures of crystalline rocks | |
| Kaveh et al. | Damage detection based on MCSS and PSO using modal data | |
| Bolander Jr et al. | Simulation of shrinkage induced cracking in cement composite overlays | |
| CN106446402B (en) | A kind of soil body dehydration cracking multi- scenarios method discrete element Fast simulation modeling method | |
| Dubois et al. | The non smooth contact dynamic method: recent LMGC90 software developments and application | |
| CN106570268A (en) | Temperature-deformation coupling analysis method and system for concrete beam structure | |
| Lone et al. | A coupled finite element-element free Galerkin approach for modeling frictional contact in engineering components | |
| CN112949065B (en) | Double-scale method, device, storage medium and equipment for simulating mechanical behavior of layered rock mass | |
| CN111259593B (en) | Method for realizing visualization of rock breaking process based on EDEM | |
| CN103454206B (en) | Simulation method of metastable state pitting of stainless steel under loading condition | |
| Zeng et al. | An adaptive peridynamics material point method for dynamic fracture problem | |
| CN117709171A (en) | High cycle fatigue failure numerical simulation method and system | |
| Xie et al. | A life prediction method of mechanical structures based on the phase field method and neural network | |
| Lukashevich et al. | Modelling and numerical solution of problems of structural mechanics with unilateral constraints and friction | |
| CN117150804A (en) | Bridge crane tamping operation simulation method, device, equipment and medium | |
| Wang et al. | Multi-physics coupling analysis on the time-dependent localized corrosion behavior of carbon steel in CO2-H2O environment | |
| Zeng et al. | An explicit phase field material point method for modeling dynamic fracture problems | |
| Sukumar et al. | Voronoi-based interpolants for fracture modelling | |
| EP0840243A2 (en) | Method capable of accurately simulating ion implantation at a high speed | |
| Martinez et al. | Distributed force simulation for arbitrarily shaped IPMC actuators | |
| CN106934729A (en) | Building Testing and appraisal method and device | |
| CN109408977B (en) | Distance potential function based deformable three-dimensional convex polyhedron block discrete unit method | |
| Hu et al. | Numerical simulation of elastic deformation based on peridynamic differential operator | |
| Ji et al. | A micromechanical finite element model for predicting the fatigue life of heterogenous adhesives | |
| CN103235865B (en) | Based on the moving contact closing speed qualification rate Forecasting Methodology of the contactor of Monte Carlo simulation |
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 |