CN112036071A - Shot peening strengthening numerical simulation method considering surface roughness and coupling finite element and discrete element - Google Patents
Shot peening strengthening numerical simulation method considering surface roughness and coupling finite element and discrete element Download PDFInfo
- Publication number
- CN112036071A CN112036071A CN202011027056.8A CN202011027056A CN112036071A CN 112036071 A CN112036071 A CN 112036071A CN 202011027056 A CN202011027056 A CN 202011027056A CN 112036071 A CN112036071 A CN 112036071A
- Authority
- CN
- China
- Prior art keywords
- surface roughness
- shot peening
- shot
- numerical simulation
- abaqus
- 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
Images
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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention relates to a shot peening strengthening numerical simulation method for coupling a finite element and a discrete element considering surface roughness, which comprises the following steps: (1) establishing an autocorrelation function representing the real surface roughness of the part based on Gaussian distribution, and extracting data information of the surface roughness by using MATLAB mathematical software; (2) importing the data information of the surface roughness into PROE three-dimensional modeling software, and performing surface smoothing treatment through fitting to establish a three-dimensional shell model representing the surface roughness of the part; (3) importing the three-dimensional shell model into ABAQUS finite element software, and establishing a part finite element model considering surface roughness on an ABAQUS/CAE platform; (4) the method comprises the steps of creating a particle generator by modifying an inp file of ABAQUS, establishing a discrete element model of a shot flow in a shot peening strengthening process, and controlling related parameters of the shot flow according to actual shot peening working conditions; (5) carrying out shot peening numerical simulation calculation on the shot peening process of the part by using an ABAQUS/Explicit dynamics algorithm, wherein the shot peening numerical simulation calculation is carried out by coupling a finite element and a discrete element considering the surface roughness.
Description
Technical Field
The invention relates to a shot peening strengthening numerical simulation method in which a finite element and a discrete element are coupled and the surface roughness is considered, which considers the influence of the surface roughness of a target body on the residual stress of shot peening and also considers the influence of random discrete shots on shot peening and belongs to the field of aerospace surface engineering.
Background
The aerospace parts have complex structures and working conditions, and particularly, the blades, transmission beams, compressor disks, blades and the like of an aero-engine are easy to fatigue fracture and wear out to lose efficacy. Shot peening is a surface strengthening technique widely used in the domestic and foreign aviation industry, and can be used to improve the fatigue performance and the service life of materials. After shot blasting, the surface of the material can generate plastic deformation, high residual stress is formed on the surface layer, and the residual stress can partially offset the stress applied by the working load, so that the expansion of fatigue cracks is delayed, and the fatigue strength of the part is improved. At present, the influence of the surface roughness of the material on the residual stress of shot blasting is usually ignored in shot blasting strengthening numerical simulation, and the shots are formed by manual modeling, so that the number of the shots is small, and the collision positions are relatively fixed. All of the above factors can bias the prediction result from reality. Therefore, it is necessary to establish a shot peening numerical simulation method in which a finite element and a discrete element are coupled in consideration of surface roughness, and the influence of each parameter on shot peening numerical simulation is comprehensively considered.
The prior document "Ahmad A S, Wu Y, Gong H.coupled fine and discrete element mapping based on Johnson-Cook material model [ J ]. Proceedings of the institute of Mechanical Engineers Part L Journal of Materials Design and Applications,2020(5): 146442072092121." establishing coupled shot peening of finite elements and discrete elements based on ABAQUS software, but without establishing target body surface roughness and therefore without considering its effect on residual stress; "Lechunxie, Juong Zhang, CenboXiong, Lihong Wu, Chuanhai Jiang, Weijie Lu. investment on experiments and numerical modeling of the residual stress distribution in the formed surface layer of Ti-6Al-4V after shot peening [ J ]. Materials and Design,2012, 41" establishes shot finite element peening based on ABAQUS software by directly establishing the positions of the shot and the target body, and thus does not consider the influence of shot randomness on residual stress.
Disclosure of Invention
The technical scheme of the invention is as follows: the method is characterized in that the surface roughness and the finite element are coupled with the discrete element to introduce shot peening strengthening, so that the residual stress under the real working condition can be accurately predicted, and the accuracy of parameters required during the shot peening processing of the aerospace surface engineering can be further improved.
The technical scheme of the invention is as follows: a shot peening strengthening numerical simulation method considering coupling of finite elements and discrete elements of surface roughness is characterized in that the real surface roughness simulation of parts is represented on the basis of an autocorrelation function of Gaussian distribution, and the influence of shot size, shot speed, nozzle flow and the number of shots is considered on the basis of random simulation of discrete element shots, so that the influence of target body roughness and the discrete shots on the residual stress of the surface of a material is considered at the same time. The method comprises the following implementation steps:
the method comprises the steps of firstly, representing the real surface roughness of a part based on an autocorrelation function of Gaussian distribution, and extracting data information of the surface roughness by using MATLAB mathematical software; and generating Gaussian random surface roughness data information by using a probability distribution function and an autocorrelation function, extracting the data information of the surface roughness, and sequentially writing the data information into a text file (xxx. txt) according to the sequence of surface nodes.
And secondly, importing the data information of the surface roughness into PROE three-dimensional modeling software, and establishing a three-dimensional shell model representing the surface roughness of the part through proper smoothing treatment.
And thirdly, introducing the three-dimensional shell model into ABAQUS/CAE, establishing a part considering surface roughness and establishing a shot blasting finite element model of the nozzle under different parameters on an ABAQUS/CAE platform, wherein the parameters of the nozzle comprise the size, the angle and the distance between the nozzle and the surface, and finally assembling and meshing the target body and the nozzle.
Fourthly, carrying out control variables on key parameters of a Particle Generator in ABAQUS, giving the number of the shots, the size of the shots, the speed of the shots and the flow rate of a nozzle; defining the contact of the projectile with the same quality characteristic, the contact of the projectile with the projectile established by different quality characteristics, and the contact of the projectile with the target body respectively by different quality characteristics.
The fifth step: carrying out shot peening numerical simulation calculation considering coupling of a finite element and a discrete element of surface roughness on the shot peening process of the part by utilizing an ABAQUS/Explicit dynamic algorithm, and obtaining the influence of the shot peening numerical simulation on residual stress under different parameters. And (4) predicting the shot blasting coverage rate according to an Avrami equation, calculating the number of the shots required when the shot blasting coverage rate reaches 98%, and setting the time step.
Compared with the prior art, the invention has the advantages that:
(1) the rough surface simulation method with the autocorrelation function based on the Gaussian distribution condition can effectively simulate the roughness of the surface of the target body, and the conventional shot peening numerical simulation is difficult to simulate random surface roughness, so that the change of the surface roughness cannot be taken into account.
(2) The method is based on a finite element and discrete element coupling simulation method, can simulate the flow of an actual nozzle and the randomness and the probability of the shot, is simple and convenient to operate, and better meets the actual situation of shot peening. The existing shot peening simulation method mainly gives shot generation positions and shot directions, the number of shots needs to be manually added, the random probability of shot peening is not considered, the operation of establishing a shot model is complex, and the numerical simulation accuracy of the residual stress on the surface of a target body cannot meet the requirement.
Drawings
FIG. 1 is a flow chart of a shot peening numerical simulation method of the present invention in which finite elements and discrete elements are coupled to take into account surface roughness;
FIG. 2 is a graphical representation of the extraction of different rough surfaces by autocorrelation function simulations under Gaussian distribution conditions; (a) a roughness map (b) with a 19-size surface root mean square of 0.05 was generated, a roughness map (c) with a 50-size surface root mean square of 0.03 was generated, and a roughness map with a 50-size surface root mean square of 0.01 was generated;
FIG. 3 is a processed rough surface data file;
FIG. 4 is a three-dimensional target body model constructed from three surface data; (a) generating a three-dimensional target body model with the size of 19 surfaces and the root mean square of 0.05, (b) generating a three-dimensional target body model with the size of 50 surfaces and the root mean square of 0.03, (c) generating a three-dimensional target body model with the size of 50 surfaces and the root mean square of 0.01;
FIG. 5 is a three-dimensional target volume mesh generation with a 19 surface root mean square of 0.05;
FIG. 6 is a key parameter of Particle Generator in ABAQUS;
FIG. 7 is a shot simulation of the same size shot;
FIG. 8 is a shot simulation of different size shots and an inclined 45 degree nozzle;
FIG. 9 shows a dimple obtained by numerical simulation of shot blasting in the center area (1/36, full size) of the target body;
FIG. 10 shows the residual stress after numerical simulation of shot peening in the center area (1/36, full size) of the target body;
Detailed Description
The shot peening numerical simulation method of the present invention, in which finite elements and discrete elements are coupled to each other with consideration of surface roughness, will be further described by way of example with reference to the accompanying drawings.
As shown in fig. 1, the implementation process of the present invention is as follows:
firstly, based on rough surface simulation with an autocorrelation function under a Gaussian distribution condition, obtaining and extracting surface data of a target body; generating a Gaussian random rough surface value by using a probability distribution function and an autocorrelation function; the graph showing the extraction of different rough surfaces simulated by the autocorrelation function under the Gaussian distribution condition is shown in FIG. 2. The functional model for simulating a gaussian surface is as follows:
in the formula taux＝1,···，N；τy＝1,···，N。η(x+τx,y+τy) A random sequence that is gaussian; h (τ)x,τy) Is the impulse response function of the filter.
The engineering surface profile is given an exponential autocorrelation function expressed as:
R(τx，τy)＝σ2exp(-2.3((τx，βx)2+(τy，βy)2)1/2)
where σ is the root mean square roughness of the surface, βx、βyRespectively, the relevant length in the x direction and the y direction; when beta isx＝βyThe time surface is anisotropic; when beta isx≠βyThe surface is anisotropic.
Secondly, a three-dimensional shell model for representing the surface roughness of the part is established by using PROE three-dimensional modeling software, the data information (xxx. txt file) of the surface roughness in the first step is imported into the PROE, the rough surface data needs to be preprocessed before the import so as to be identified by the PROE software, and the processed rough surface data file is shown in figure 3. Using boundary blending commands in the pro software, the surface node coordinate information is fitted to a spatial three-dimensional shell model, and a three-dimensional target body model with surface roughness is created as shown in fig. 4.
And thirdly, importing a three-dimensional shell model into ABAQUS/CAE, establishing a part considering surface roughness and establishing a shot blasting finite element model of the nozzle under different parameters on an ABAQUS/CAE platform, wherein the parameters of the nozzle comprise the size, the angle and the distance between the nozzle and the surface, controlling the parameters to simulate the shot blasting strengthening result under various working conditions in reality, and finally assembling the target body and the nozzle and meshing the target body as shown in figure 5.
Fourthly, creating a Particle Generator by modifying an inp file of the ABAQUS, establishing a discrete element model of the shot flow in the shot peening process, and carrying out control variables on key parameters of a Particle Generator in the ABAQUS, wherein the key parameters of the Particle Generator are shown in FIG. 6, the number of shots is given, the size, the speed and the nozzle flow rate of the shots are given, the contact between the shots with the same quality characteristics are defined, the contact between the shots with different quality characteristics are established, the contact between the shots with different quality characteristics and the target body is simulated, the shot with the same size is simulated, the shot with different size is shown in FIG. 7, and the angle of a nozzle is simulated, and the angle of the shot with different size is shown in FIG. 8.
The fifth step: carrying out shot peening numerical simulation calculation considering coupling of a finite element and a discrete element of surface roughness on the shot peening process of the part by utilizing an ABAQUS/Explicit dynamic algorithm to obtain the residual stress of shot peening under different variables. Firstly, numerically simulating the shot to impact a target body, calculating the area of a single pit left by the shot, predicting the shot blasting coverage rate according to an Avrami equation, calculating the number of the shots required when the shot blasting coverage rate reaches 98%, setting the time step, and finally carrying out simulation analysis, wherein the shot blasting coverage rate is as shown in figure 9. The influence of different variables on the residual stress is obtained according to the model analysis result, and the residual stress is shown in fig. 10.
Claims (4)
1. A shot peening strengthening numerical simulation method for coupling a finite element and a discrete element considering surface roughness is characterized in that: the method comprises the following steps:
step (1): obtaining and extracting target body surface data based on rough surface simulation with an autocorrelation function under a Gaussian distribution condition; the data information of the surface roughness of the part is obtained by MATLAB software calculation; the extraction of the surface roughness data information is to write a text file (xxx. txt) in sequence according to the sequence of surface nodes;
step (2): establishing a three-dimensional shell model representing the surface roughness of the part by using PROE three-dimensional modeling software, importing the data information (xxx. txt file) of the surface roughness in the step (1) into PROE, and establishing a spatial three-dimensional shell model according to the surface node coordinate information; the surface node coordinate information needs to be preprocessed; the surface node coordinate information can represent the surface roughness of the part;
and (3): introducing the three-dimensional shell model into ABAQUS/CAE, establishing a part considering surface roughness and establishing a shot blasting finite element model of the nozzle under different parameters on an ABAQUS/CAE platform; the nozzle parameters comprise nozzle size, nozzle angle and nozzle-surface distance;
and (4): the method comprises the steps of creating a particle generator by modifying an inp file of ABAQUS, establishing a discrete element model of a shot flow in a shot peening strengthening process, and controlling related parameters of the shot flow according to actual shot peening working conditions; the parameters refer to the number of the shots, the size of the shots, the speed of the shots and the flow rate of the nozzle;
and (5): carrying out shot peening numerical simulation calculation on the shot peening process of the part by using an ABAQUS/Explicit dynamics algorithm, wherein the shot peening numerical simulation calculation is carried out by coupling a finite element and a discrete element considering the surface roughness.
2. A shot peening numerical simulation method of coupling a finite element and a discrete element considering surface roughness according to claim 1, wherein: in the step (1), the surface data of the target body is obtained and extracted as follows:
and generating a Gaussian random rough surface numerical value by using a probability distribution function and an autocorrelation function, and extracting data information of the surface roughness by using MATLAB mathematical software.
3. A shot peening numerical simulation method of coupling a finite element and a discrete element considering surface roughness according to claim 1, wherein: in the step (2), data information (xxx. txt file) of the surface roughness is imported, a rough surface is generated by using boundary mixing, and a three-dimensional shell model representing the surface roughness of the part is established.
4. A shot peening numerical simulation method of coupling a finite element and a discrete element considering surface roughness according to claim 1, wherein: in the step (4), the parameter control is to write control and parameter modification on the keywords in the ABAQUS/inp file.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011027056.8A CN112036071A (en) | 2020-09-25 | 2020-09-25 | Shot peening strengthening numerical simulation method considering surface roughness and coupling finite element and discrete element |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011027056.8A CN112036071A (en) | 2020-09-25 | 2020-09-25 | Shot peening strengthening numerical simulation method considering surface roughness and coupling finite element and discrete element |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112036071A true CN112036071A (en) | 2020-12-04 |
Family
ID=73574447
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011027056.8A Pending CN112036071A (en) | 2020-09-25 | 2020-09-25 | Shot peening strengthening numerical simulation method considering surface roughness and coupling finite element and discrete element |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112036071A (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112668087A (en) * | 2021-01-11 | 2021-04-16 | 武汉理工大学 | Ballasted track modeling analysis method and system based on finite element and discrete element coupling |
| CN115157128A (en) * | 2022-06-15 | 2022-10-11 | 西北工业大学 | Method and device for reconstructing surface appearance of shot blasting part |
| CN115455756A (en) * | 2022-08-23 | 2022-12-09 | 成都飞机工业（集团）有限责任公司 | Method, system and application for acquiring shot blasting speed by using roughness |
| CN115828472A (en) * | 2023-02-15 | 2023-03-21 | 太原理工大学 | Method for simulating residual stress of surface of barreled workpiece |
| WO2023134303A1 (en) * | 2022-01-11 | 2023-07-20 | 上海飞机制造有限公司 | Method for obtaining shot distribution in shot peening process |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111651921A (en) * | 2020-06-02 | 2020-09-11 | 重庆大学 | Shot peening strengthening surface integrity parameter prediction method based on material real state |
-
2020
- 2020-09-25 CN CN202011027056.8A patent/CN112036071A/en active Pending
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111651921A (en) * | 2020-06-02 | 2020-09-11 | 重庆大学 | Shot peening strengthening surface integrity parameter prediction method based on material real state |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112668087A (en) * | 2021-01-11 | 2021-04-16 | 武汉理工大学 | Ballasted track modeling analysis method and system based on finite element and discrete element coupling |
| WO2023134303A1 (en) * | 2022-01-11 | 2023-07-20 | 上海飞机制造有限公司 | Method for obtaining shot distribution in shot peening process |
| CN115157128A (en) * | 2022-06-15 | 2022-10-11 | 西北工业大学 | Method and device for reconstructing surface appearance of shot blasting part |
| CN115157128B (en) * | 2022-06-15 | 2024-01-30 | 西北工业大学 | Method and device for reconstructing surface morphology of shot blasting part |
| CN115455756A (en) * | 2022-08-23 | 2022-12-09 | 成都飞机工业（集团）有限责任公司 | Method, system and application for acquiring shot blasting speed by using roughness |
| CN115828472A (en) * | 2023-02-15 | 2023-03-21 | 太原理工大学 | Method for simulating residual stress of surface of barreled workpiece |
| CN115828472B (en) * | 2023-02-15 | 2023-05-26 | 太原理工大学 | Method for simulating residual stress on surface of barreling finishing workpiece |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112036071A (en) | Shot peening strengthening numerical simulation method considering surface roughness and coupling finite element and discrete element | |
| CN111832102B (en) | Novel composite material structure optimization design method under high-dimensional random field condition | |
| Zhang et al. | A study of airfoil parameterization, modeling, and optimization based on the computational fluid dynamics method | |
| CN110059330B (en) | Method and system for authoring simulated scenes | |
| CN105893694A (en) | Complex system designing method based on resampling particle swarm optimization algorithm | |
| Manap et al. | Computer simulation of cold sprayed deposition using smoothed particle hydrodynamics | |
| CN111797554B (en) | Turbine tongue-and-groove shot blasting discrete element-finite element coupling multi-scale simulation method | |
| WO2007087490A2 (en) | Sector meshing and neighbor searching for object interaction simulation | |
| CN113673027A (en) | Hypersonic aircraft aerodynamic load optimization design method based on proxy model | |
| CN108428256B (en) | Soft tissue deformation simulation method based on self-adaptive grid refinement of softness | |
| Menzel et al. | Representing the change-free form deformation for evolutionary design optimization | |
| CN114036631A (en) | Spacecraft autonomous rendezvous and docking guidance strategy generation method based on reinforcement learning | |
| CN116956684A (en) | Ultrasonic shot peening strengthening residual stress and deformation analysis method | |
| CN115017647A (en) | External foreign body collision and wear assessment method based on DEM solver | |
| Henry et al. | Reduced order geometric mistuning models using principal component analysis approximations | |
| CN107577899B (en) | Three-dimensional discrete element characterization method for rock mass random structural plane | |
| CN105354387B (en) | A kind of Three-Dimensional Atmospheric Turbulence generation method based on improvement white Gaussian noise | |
| Wang et al. | Hypersonic Vehicle Aerodynamic Optimization Using Field Metamodel‐Enhanced Sequential Approximate Optimization | |
| CN111125971A (en) | Method for determining thrust uncertainty of air-breathing hypersonic aircraft | |
| Oliveira Neto et al. | Aerodynamic study of sounding rocket flows using Chimera and patched multiblock meshes | |
| Tai et al. | A numerical study in the interaction of blast wave with a wheeled armoured vehicle | |
| Sung et al. | Parallel Structure Design Method for Mass Spring Simulation | |
| Hattangady et al. | Towards realistic automated 3D modelling of metal forming problems | |
| CN115048730B (en) | Axial symmetry supersonic velocity spray pipe optimization design method and device based on grid displacement | |
| Romlay et al. | Failure prediction for automotive suspension springs using Gaussian and Monte Carlo 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 | ||
| CB03 | Change of inventor or designer information |
Inventor after: Wang Cheng Inventor after: Zhou Bin Inventor after: Li Kaifa Inventor after: Zhou Yijun Inventor after: Yao Shangcheng Inventor before: Li Kaifa Inventor before: Wang Cheng Inventor before: Zhou Yijun Inventor before: Yao Shangcheng |
|
| CB03 | Change of inventor or designer information |