Finite element simulation method for calculating riveting process parameters and riveting performance test parameters

Abstract

The present invention relates to a method for calculating riveting process parameters and riveting performance test parameters by finite element simulation. The method comprises: using UG software to establish a working process model of rivets and flaring ejector pins; defining analysis methods and workflows in ANSYS software, and defining materials of flaring rivets and flaring ejector pins; performing meshing and contact setting on flaring rivets and flaring ejector pins; defining constraints and loads on flaring rivets and flaring ejector pin models, wherein the constraints are fixed to the bottom surface of the rivets, and fixed loads are applied to the ejector pins respectively; performing operation setting and calculation; comparing the calculation results with the assembly drawing parameters, and gradually adjusting the applied load until it matches the assembly parameters, and finally determining the applied fixed load. The method ensures the reliability and safety of assembly of key components and the scientific nature of the test during rivet manufacturing, and also greatly reduces the cost compared with the physical test method.

Description

Technical Field

[0001] The present invention belongs to the technical field of computer-aided design, and is used for calculating riveting force process parameters and rivet manufacturing test parameters of a turboshaft engine turbine disk and turbine blade, and specifically relates to a method for calculating riveting process parameters and riveting performance test parameters through finite element simulation. Background Art

[0002] Inconel 600 semi-hollow rivets are mainly used for riveting turbine disks and turbine blades of turboshaft engines (see Figure 1 ), which is a must-replace part during aircraft engine overhaul and consumes a lot of energy. Rivets bear axial loads and smaller radial loads during operation, and their failure mode is mainly micro-motion wear. Their failure may cause blades to fall off, rotor blades to be damaged, and even cause serious consequences such as the scrapping of the entire engine. The riveting process parameters during the assembly of turbine disks and turbine blades of turboshaft engines and the riveting performance test parameters during manufacturing are crucial. Using traditional physical tests to determine the parameters is not only time-consuming and may damage the turbine disk, but also unable to accurately determine the applied fixed load. Finite element simulation calculations can avoid these shortcomings. Summary of the invention

[0003] The present invention provides a method for finite element simulation calculation of riveting process parameters during assembly and riveting performance test parameters during manufacturing of Inconel 600 semi-hollow rivets, which solves the problem that traditional physical test to determine parameters is not only time-consuming and may damage the turbine disk, but also unable to accurately determine the applied fixed load.

[0004] The present invention is achieved by the following technical solutions:

[0005] A method for calculating riveting process parameters and riveting performance test parameters by finite element simulation, including Inconel600 semi-hollow rivets, and the specific steps are as follows:

[0006] (1) Use UG software to establish a model of the working process of the flaring rivet, including the rivet and the flaring ejector rod;

[0007] (2) According to the work requirements, define the analysis method and workflow in ANSYS software and adopt the transient dynamics workflow;

[0008] (3) Define the materials of the flaring rivet and flaring ejector pin in ANSYS software;

[0009] (4) Pre-processing of flaring rivets and flaring ejector pins, including meshing and contact setting;

[0010] (5) Define constraints and loads for the flared rivet and flared ejector models. Constraints are applied to the bottom surface of the rivet and different fixed loads are applied to the ejector.

[0011] (6) Operation settings, including step size, number of steps, and solution time;

[0012] (7) Calculate the results for different loads respectively;

[0013] (8) Compare the calculation results with the assembly drawing parameters, and gradually adjust the applied load until it matches the assembly parameters, and finally determine the applied fixed load.

[0014] In step (4), the mesh is divided into hexahedral elements with a total of 2100-2150 elements and 6600-6800 nodes, and the contact setting adopts bound contact.

[0015] In step (5), the different fixed loads are 100kg, 200kg, 400kg, and 600kg respectively.

[0016] In step (6), set the step size to 10, the number of steps to 10, and the solution time to 0.1S.

[0017] The method of the invention can effectively ensure the reliability and safety of assembly of key parts and components and the scientific nature of the test during rivet manufacturing, and greatly reduces the cost compared with the physical test method. BRIEF DESCRIPTION OF THE DRAWINGS

[0018] Figure 1 It is a schematic diagram of the assembly of rivets, blades and discs of the present invention.

[0019] Specific embodiment

[0020] like Figure 1 As shown, the method of calculating the riveting process parameters of Inconel 600 semi-hollow rivets during assembly and the riveting performance test parameters during manufacturing by finite element simulation of a turboshaft engine is used as an example to illustrate.

[0021] Part name: semi-hollow rivet, material: Inconel 600, material standard: AMS 5687, the rivet is solution heat treated and is used for riveting between turbine disks and blades of turboshaft engines.

[0022] According to the finite element simulation calculation steps:

[0023] (1) Using UG software, establish the rivet and flaring ejector model according to the semi-hollow rivet parts drawing and the riveting tooling drawing configured by the riveting equipment.

[0024] (2) According to work requirements, the analysis method and workflow are defined in ANSYS software, and the transient dynamics workflow is adopted.

[0025] (3) Define the materials of the flared rivet and flared ejector in ANSYS software. Rivet material: Inconel600, density: 8430kg / m 3 , elastic modulus E = 205GPa, tensile strength σ b =585MPa, yield strength: σ 0.2 =240MPa, Poisson's ratio 0.34.

[0026] Flaring mandrel 40Cr: Density: 7820kg / m 3 , elastic modulus E = 211 GPa, tensile strength σ b =1085MPa, yield strength: σ 0.2 =978MPa, Poisson's ratio 0.277.

[0027] (4) Pre-processing of the flared rivets and flared ejector pins, including meshing of hexahedral elements with a total of 2128 elements and 6739 nodes; the contact setting uses bound contact.

[0028] (5) Constraints and loads are defined for the expansion rivet and expansion ejector models. The constraints are fixed on the bottom surface of the rivet, and fixed loads are applied to the ejector (the loads are 100 kg, 200 kg, 400 kg, and 600 kg, respectively).

[0029] (6) Operation settings, including step size 10, number of steps 10, solution time 0.1 s, etc.

[0030] (7) Calculate the results for different loads respectively.

[0031] (8) Compare the calculation results with the assembly drawing parameters, and gradually adjust the applied load until it matches the assembly parameters. Finally, the riveting force is determined to be 397kg to 419kg.

Claims

1. A method for calculating riveting process parameters and riveting performance test parameters by finite element simulation, including Inconel 600 semi-hollow rivets, characterized in that: Specific steps as follows: (1) Use UG software to establish a model of the working process of the flaring rivet, including the rivet and the flaring ejector rod; (2) According to the work requirements, define the analysis method and workflow in ANSYS software and adopt the transient dynamics workflow; (3) Define the materials of the flaring rivet and flaring ejector pin in ANSYS software; (4) Pre-processing of flaring rivets and flaring ejector pins, including meshing and contact setting; (5) Define constraints and loads for the flared rivet and flared ejector models. Constraints are applied to the bottom surface of the rivet and different fixed loads are applied to the ejector. (6) Operation settings, including step size, number of steps, and solution time; (7) Calculate the results for different loads respectively; (8) Compare the calculation results with the assembly drawing parameters, and gradually adjust the applied load until it matches the assembly parameters, and finally determine the applied fixed load.

2. The method for calculating riveting process parameters and riveting performance test parameters by finite element simulation according to claim 1 is characterized by: In step (4), the mesh is divided into hexahedral elements with a total of 2100-2150 elements and 6600-6800 nodes, and the contact setting adopts bound contact.

3. The method for calculating riveting process parameters and riveting performance test parameters by finite element simulation according to claim 1, characterized in that: In step (5), the different fixed loads are 100kg, 200kg, 400kg, and 600kg respectively.

4. The method for calculating riveting process parameters and riveting performance test parameters by finite element simulation according to claim 1 is characterized by: In step (6), set the step size to 10, the number of steps to 10, and the solution time to 0.1S.

Images