Wheel disc finite element calculation method and device based on finite element calculation model

Abstract

The invention relates to the technical field of aero-engines, and provides a wheel disc finite element calculation method and device based on a finite element calculation model.The method comprises the steps that the bolt number and the blade number of a calculation model are determined according to the bolt number and the blade number of an original geometric model of a blade disc, and a cyclic symmetry model is obtained according to the blade number of the calculation model; the cyclic symmetry model is cut from the original geometric model of the blade disc without the original bolt holes, simulated bolt holes are formed in the cyclic symmetry model, the bolts are arranged in the simulated bolt holes, the cyclic symmetry model is imported into limited element analysis software to adjust parameters of the bolts and the simulated bolt holes, and the parameters of the bolts and the simulated bolt holes are obtained; and performing finite element calculation according to the adjusted bolt, the parameters of the simulated bolt hole and preset parameters. Compared with a traditional cyclic symmetry model, the cyclic symmetry model is simple in structure, so that the calculation amount is smaller, and the finite element calculation efficiency is higher. Compared with the technology of directly deleting bolts and bolt holes, the finite element calculation result is higher in precision.

Description

Technical Field

[0001] The present invention belongs to the technical field of aeroengines, and particularly relates to a finite element calculation method, device and storage medium for a disk based on a finite element calculation model. Background Art

[0002] Disks of aeroengines are often connected by bolts. For the finite element displacement and stress calculation of aeroengine disks, if a complete model of multiple disks is calculated using a multiple disk bolt connection method, not only is the calculation volume large, but also, due to the large number of contacts between disks, the calculation is difficult to converge. If a cyclic symmetric model is used, the number of blades and bolts need to have a large greatest common divisor. However, in actual situations, the number of blades and bolts often do not have a large greatest common divisor, or the greatest common divisor is 2 or 3. When the greatest common divisor is 2 or 3, even if a cyclic symmetric model is used, only a 1 / 2 or 1 / 3 cyclic symmetric model can be adopted, and there are still problems of large calculation volume, many contacts, and difficult convergence.

[0003] In response to the above problems, although the prior art has proposed corresponding solutions, however, there are large errors in the calculation results of displacement and stress for the disks of various aeroengines. Summary of the Invention

[0004] In response to the above problems, the present invention proposes a finite element calculation method for a disk based on a finite element calculation model, and the method includes:

[0005] Determine the number of bolts and blades of the calculation model according to the number of bolts and blades of the original geometric model of the blade disk;

[0006] Obtain a cyclic symmetric model according to the number of blades of the calculation model, and the cyclic symmetric model is cut from the original geometric model of the blade disk after deleting the original bolt holes;

[0007] Open simulated bolt holes on the cyclic symmetric model, and place the bolts in the simulated bolt holes;

[0008] Import the cyclic symmetric model into finite element analysis software to adjust the parameters of the bolts and the simulated bolt holes;

[0009] Perform finite element calculation according to the adjusted parameters of the bolts and the simulated bolt holes and preset parameters.

[0010] Optionally, the determining the number of bolts and blades of the calculation model according to the number of bolts and blades of the original geometric model of the blade disk includes:

[0011] When the ratio of the number of bolts to the number of blades of the original geometric model of the blade disk is within a first set value range, the calculation model contains 1 bolt and 1 blade; or

[0012] If the ratio of the number of bolts to the number of blades of the original geometric model of the blisk is less than the minimum value of the first set value range, the calculation model contains 1 bolt and 1 blade or multiple blades, and the number of blades of the calculation model does not exceed the greatest common divisor of the preset value, where the preset value is the ratio of the number of blades to the number of bolts of the original geometric model of the blisk; or

[0013] If the ratio of the number of bolts to the number of blades of the original geometric model of the blisk is greater than the maximum value of the first set value range, the calculation model contains multiple bolts and 1 blade, and the number of bolts of the calculation model is the ratio of the number of bolts to the number of blades of the original geometric model of the blisk. If this ratio is not an integer, round it to the nearest integer.

[0014] Optionally, the cyclic symmetric model is cut from the original geometric model of the blisk after deleting the original bolt holes, and includes:

[0015] Form a first surface according to the axis of the original geometric model of the blisk and a line preset on the rim of the original geometric model of the blisk;

[0016] The first surface rotates around the axis of the original geometric model of the blisk by a preset angle to form a second surface;

[0017] Cut the original geometric model of the blisk, and take the sector model between the first surface and the second surface as the cyclic symmetric model.

[0018] Optionally, the first surface rotates around the axis of the original geometric model of the blisk by a preset angle to form a second surface, including:

[0019] The preset angle is 360°÷y×z;

[0020] y is the number of blades of the original geometric model of the blisk, and z is the number of blades of the calculation model.

[0021] Optionally, simulating bolt holes are opened on the cyclic symmetric model, including:

[0022] Adding simulating bolt holes at the circumferential middle position of the mounting edge of the cyclic symmetric model, where the position of the simulating bolt holes on the cyclic symmetric model is determined according to the radius of the center of the simulating bolt holes being equal to the radius of the center of the original bolt holes;

[0023] The volume of the simulating bolt holes × the number of simulating bolt holes × y÷z = the volume of the corresponding original bolt holes of the original geometric model of the blisk × the number of original bolt holes, where y is the number of blades of the original geometric model of the blisk and z is the number of blades of the calculation model.

[0024] Optionally, importing the cyclic symmetric model into finite element analysis software to adjust the parameters of the bolts and the simulating bolt holes, including:

[0025] Modify the density of the bolt material to obtain the target density of the bolt material;

[0026] The target density × the number of simulated bolt holes × y ÷ z = the true material density of the bolt × j;

[0027] Set the calculated interference amount at the contact position between the bolt and the simulated bolt hole, and the geometric interference amount of the calculation model + the calculated interference amount = the geometric interference amount of the true model;

[0028] Among them, the geometric interference amount of the calculation model = the radius of the bolt in the calculation model - the radius of the bolt hole in the calculation model, the geometric interference amount of the true model = the actual bolt radius - the actual bolt hole radius, j is the number of bolt holes in the original geometric model of the blade disk, y is the number of blades in the original geometric model of the blade disk, and z is the number of blades in the calculation model.

[0029] Optionally, after opening the simulated bolt holes on the cyclic symmetric model, it includes:

[0030] Select the original relief groove on the original geometric model of the blade disk, and generate a sheet body at the position of the original relief groove;

[0031] Copy the sheet body into the calculation model and form an offset sheet body by offsetting a preset distance;

[0032] Rotate and copy the offset sheet body around the axis of the original geometric model of the blade disk and cut the calculation model to obtain the target number of simulated relief grooves.

[0033] Optionally, the rotating and copying the offset sheet body around the axis of the original geometric model of the blade disk and cutting the calculation model to obtain the target number of simulated relief grooves includes:

[0034] The rotation angle between two adjacent sheet bodies = 360° ÷ y × z ÷ p, where p is the target number, y is the number of blades in the original geometric model of the blade disk, and z is the number of blades in the calculation model;

[0035] The volume of each relief groove in the calculation model × p × y ÷ z = the volume of the original relief groove × the number of original relief grooves, where the target number is equal to the number of bolts in the calculation model.

[0036] The present invention also provides a disk finite element calculation device based on a finite element calculation model, including:

[0037] A determination module, configured to determine the number of bolts and the number of blades in the calculation model according to the number of bolts and the number of blades in the original geometric model of the blade disk;

[0038] An acquisition module, configured to acquire a cyclic symmetry model according to the number of blades of the calculation model, where the cyclic symmetry model is cut from an original geometric model of a bladed disk with the original bolt holes removed;

[0039] An opening and assembly module, configured to open simulated bolt holes on the cyclic symmetry model and set the bolts in the simulated bolt holes;

[0040] A pre-processing module for calculation, configured to import the cyclic symmetry model into finite element analysis software to adjust the parameters of the bolts and the simulated bolt holes;

[0041] A calculation module, configured to perform finite element calculation according to the adjusted parameters of the bolts and the simulated bolt holes and preset parameters.

[0042] The present invention also provides a storage medium, including: when the program is processed and executed, it implements the disk finite element calculation method based on a finite element calculation model as described above.

[0043] The disk finite element calculation method based on a finite element calculation model of the present invention has the following advantages compared with the prior art:

[0044] By determining the number of bolts and the number of blades of the calculation model, it is possible to determine the number of blades on the acquired cyclic symmetry model. Moreover, compared with the traditional complete model or 1 / 2 or 1 / 3 cyclic symmetry models, the structure of the cyclic symmetry model is simpler, so the amount of calculation is smaller, thereby improving the efficiency of finite element calculation. In addition, the convergence of finite element calculation gradually deteriorates as the number of contact pairs increases. Since the structure of the cyclic symmetry model is simpler, and the parameters of the bolts and the simulated bolt holes can also be adjusted, the number of contact pairs can be correspondingly reduced, thereby improving the convergence of the calculation.

[0045] Other features and advantages of the present invention will be described in the following description of the specification, and in part, will be obvious from the description of the specification, or will be understood by implementing the present invention. The objectives and other advantages of the present invention can be achieved and obtained through the structures pointed out in the specification and the drawings. Description of the Drawings

[0046] In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings required for use in the description of the embodiments or the prior art. Obviously, the following drawings are some embodiments of the present invention. For those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative efforts.

[0047] Figure 1 Shows a schematic flow chart of the disk finite element calculation method based on a finite element calculation model in an embodiment of the present invention;

[0048] Figure 2 Shows the original schematic diagram of the finite element calculation model device of the wheel in the embodiment of the present invention;

[0049] Figure 3 Shows the schematic diagram of the finite element calculation model device of the wheel in the embodiment of the present invention without holes;

[0050] Figure 4 Shows the schematic diagram of the finite element calculation model device of the wheel in the embodiment of the present invention with holes;

[0051] Figure 5 Shows the schematic diagram of the finite element calculation model device of the wheel in the embodiment of the present invention with grooves;

[0052] Figure 6 Shows the schematic diagram of the finite element calculation model device of the wheel in the embodiment of the present invention with relief grooves;

[0053] Figure 7 Shows the schematic diagram of the assembly of bolts and screw holes of the finite element calculation model device of the wheel in the embodiment of the present invention.

[0054] The following are the reference numerals of the present invention:

[0055] 1. Original geometric model of the blade disc; 2. Number of blades in the calculation model; 21. Simulated bolt holes; 30. Original relief groove; 31. Sheet body, 32. Offset sheet body; 3. Simulated relief groove; 4. Bolt. Detailed implementation manners

[0056] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are some, but not all, of the embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0057] As Figure 1 shown, the present invention provides a finite element calculation method for a wheel based on a finite element calculation model, and the method includes:

[0058] Step S10: Determine the number of bolts 41 and the number of blades in the calculation model according to the number of bolts 4 and the number of blades of the original geometric model of the blisk. By calculating the number of blades 2 in the model, the number of blades on the obtained cyclic symmetry model can be determined. Compared with the traditional complete model or 1 / 2 or 1 / 3 cyclic symmetry model, the structure of the cyclic symmetry model is simpler, reducing the amount of calculation, thereby affecting the calculation efficiency and achieving an improvement in the efficiency of finite element calculation. Additionally, simulate the number of bolts 4 corresponding to the number of bolt holes 21 to determine the number of contact pairs formed between the two. Since the structure of the cyclic symmetry model is simpler, and by adjusting the parameters of the bolts 4 and the simulated bolt holes 21, the number of contact pairs can be reduced. The convergence of finite element calculation gradually deteriorates as the number of contact pairs increases. Therefore, by calculating to determine a reduction in the number of contact pairs, the convergence of the calculation can be improved. After verification, finite element calculation is performed according to the adjusted parameters of the bolts 4, the simulated bolt holes 21, and the preset parameters. The stress values and displacement values of the blisk obtained from the finite element calculation have higher accuracy. The verification can refer to the content in Table 3 below. Optionally, the blade can also be a disk with a tenon groove or a disk in other structural forms. The bolt hole can also be a vent hole.

[0059] In one embodiment, determining the number of bolts 4 and the number of blades 2 in the calculation model according to the number of bolts 4 and the number of blades of the original geometric model of the blisk includes:

[0060] If the ratio of the number of bolts 4 to the number of blades of the original geometric model 1 of the blisk is within the first set value range, the calculation model contains only 1 bolt 4 and 1 blade; or

[0061] If the ratio of the number of bolts 4 to the number of blades of the original geometric model 1 of the blisk is less than the minimum value of the first set value range, the calculation model contains 1 bolt 4 and 1 blade (in some cases, the number of blades may also be one) or multiple blades. The number of blades in the calculation model is the largest divisor not exceeding the preset value, and the preset value is the ratio of the number of blades to the number of bolts 4 of the original geometric model 1 of the blisk. It should be noted that the preset value is the upper limit of the number of blades. Therefore, it is also possible that the number of blades may be one; or

[0062] If the ratio of the number of bolts 4 to the number of blades of the original geometric model 1 of the blisk is greater than the maximum value of the first set value range, the calculation model contains multiple bolts 4 and 1 blade. The number of bolts 4 in the calculation model is the ratio of the number of bolts 4 to the number of blades of the original geometric model 1 of the blisk. If this ratio is not an integer, round it to the nearest integer.

[0063] Specifically, the calculation model is usually a cyclic symmetry model with only 1 bolt 4 and 1 blade, or 1 bolt 4 and multiple blades, or multiple bolts 4 and 1 blade. If the number of bolts 4 in the original geometric model 1 of the blade disc is close to the number of blades, that is, the ratio of the number of bolts 4 to the number of blades in the original geometric model 1 of the blade disc is within the first set value range (the first set value is between 0.66 and 1.5), the calculation model contains only 1 bolt 4 and 1 blade; if the ratio is much smaller than the number of blades, the calculation model contains 1 bolt 4 and 1 or more blades; if the ratio is much larger than the number of blades, the calculation model contains multiple bolts 4 and 1 blade. The number of bolts 4 and blades in the calculation model is shown in Table 1 below.

[0064] Table 1

[0065]

[0066] x is the number of bolts 4 in the original geometric model 1 of the blade disc, and y is the number of blades in the original geometric model 1 of the blade disc.

[0067] Determining the number of bolts 4 and blades in the calculation model in this way can not only ensure that the calculation model is simple enough, but also make the simplification of the calculation model small enough (generally speaking, the smaller the simplification, the higher the calculation accuracy). By adjusting the number of bolt holes according to the number of bolts 4, the number of blades 2 and the number of bolt holes in the calculation model have a simple proportional relationship, such as 1:1, 1:2, 2:1, 1:3, 3:1. In this way, the calculation can be carried out with a simple cyclic symmetry model. Specifically, for example, the blade disc contains 39 blades and 36 bolt holes, and the number of blades is close to the number of bolt holes at 1:1. Therefore, the calculation model is simplified to a cyclic symmetry model with only 1 bolt 4 and 1 blade.

[0068] Step S20: Obtain a cyclic symmetry model according to the number of blades in the calculation model. The cyclic symmetry model is cut from the original geometric model 1 of the blade disc after deleting the original bolt holes. It should be noted that the original geometric model 1 of the blade disc after deleting the original bolt holes is to fill the original bolt holes. The purpose of deleting the original bolt holes is to re-open holes, and the number and diameter size of the holes need to be designed according to actual requirements, and are specifically determined according to the ratio of the number of bolts 4 to the number of blades in the original geometric model 1 of the blade disc. For the cyclic symmetry model formed by cutting the original geometric model 1 of the blade disc after deleting the original bolt holes, subsequent re-opening operations need to be carried out on it. In addition, when obtaining the cyclic symmetry model according to the number of blades in the calculation model, the obtained cyclic symmetry model should have the corresponding number of blades.

[0069] In one embodiment, cutting the cyclic symmetry model from the original geometric model 1 of the blade disc after deleting the original bolt holes includes:

[0070] A first surface is formed based on the axis of the original blade disk geometric model 1 and a line preset on the rim of the original blade disk geometric model 1. The line preset on the rim of the original blade disk geometric model 1 can be designed according to actual needs. For example, it can be a curved line or a straight line. (The process of forming the first surface can be achieved with the help of 3D software tools. Of course, it can also be operated physically. Either way can be chosen, and it is more convenient and faster to form the first surface with the help of 3D software tools).

[0071] The first surface rotates around the axis of the original blade disk geometric model 1 by a preset angle to form a second surface;

[0072] The original blade disk geometric model 1 is cut, and the sector model between the first surface and the second surface is taken as the cyclic symmetric model. It should be noted that the cyclic symmetric model selected by this method is the required cyclic symmetric model, and this selected cyclic symmetric model is more convenient. Optionally, the first surface rotates around the axis of the original blade disk geometric model 1 by a preset angle to form a second surface, where the preset angle is 360°÷y×z, y is the number of blades of the original blade disk geometric model 1, and z is the number of blades of the calculation model.

[0073] Step S30, simulate bolt holes 21 are opened on the cyclic symmetric model, and bolts 4 are arranged in the simulated bolt holes 21. It should be noted that the simulated bolt holes 21 are newly opened holes, and the aperture sizes and the numbers between the original bolt holes and the simulated bolt holes 21 are different. As Figure 7 shown, the bolts 4 (including bolts 4, nuts, washers (if any), etc.) of the original blade disk geometric model 1 are assembled into the bolt holes of the calculation model (cyclic symmetric model), ensuring that the bolts 4 and the simulated bolt holes 21 are coaxial, and all contact surfaces are in contact.

[0074] In one embodiment, opening the simulated bolt holes 21 on the cyclic symmetric model includes:

[0075] Simulated bolt holes 21 are added at the circumferential middle position of the mounting edge of the cyclic symmetric model. Among them, it is necessary to measure the radius of the center of the original bolt hole and the volume of the original bolt hole on the original blade disk geometric model 1, and determine the position of the simulated bolt hole on the cyclic symmetric model according to the radius of the center of the simulated bolt hole 21 being equal to the radius of the center of the original bolt hole. It should be noted that the radius of the bolt hole center is the distance from the bolt hole center to the axis of the blade disk, which is convenient for determining the position of the simulated bolt hole 21 on the cyclic symmetric model. Thus, it can be seen that the radius of the bolt hole center is not equal to the radius of the bolt hole.

[0076] The volume of the simulated bolt hole 21 × the number of simulated bolt holes 21 × y ÷ z = the volume of the original bolt holes corresponding to the original blade disk model × the number of original bolt holes. y is the number of blades of the original geometric model 1 of the blade disk, and z is the number of blades of the calculation model. Specifically, as Figure 2 shown, measure the radius of the center of the mounting edge bolt hole (the distance from the axis of the mounting edge bolt hole to the axis of the blade disk) and the volume of each mounting edge bolt hole (the volume of the material required to fill each mounting edge bolt hole) in the original geometric model 1 of the blade disk. As Figure 3 shown, add simulated bolt holes 21 at the circumferential middle position of the mounting edge of the calculation model to make the bolt holes corresponding to each bolt 4 coaxial. The radius of the center of the simulated bolt hole 21 is equal to the radius of the center of the bolt hole in the original model. The volume of each bolt hole in the calculation model × the number of bolts 4 in the calculation model × y ÷ z = the volume of the original bolt holes corresponding to the original blade disk model × the number of original bolt holes. y is the number of blades of the original geometric model 1 of the blade disk, and z is the number of blades of the calculation model. It should be noted that since the number of bolt holes is adjusted, in order to ensure mass equivalence and stiffness equivalence, the size of the bolt holes needs to be adjusted to ensure that the total volume of the bolt holes is equal to the original to ensure the calculation accuracy. The adjusted bolt holes are as Figure 4 shown.

[0077] Step S40, import the cyclic symmetric model into the finite element analysis software to adjust the parameters of the bolt 4 and the simulated bolt hole 21. Specifically, since the number of bolts 4 is adjusted, in order to ensure mass equivalence, the density of the bolt 4 needs to be adjusted to ensure that the mass of the bolt 4 is equal to ensure the calculation accuracy. Import the cyclic symmetric model into the general finite element analysis software, such as ANSYS software.

[0078] In one embodiment, importing the cyclic symmetric model into the finite element analysis model to adjust the parameters of the bolt 4 and the simulated bolt hole 21 includes:

[0079] Modify the material density of the bolt 4 to obtain the target density of the bolt 4 material;

[0080] The target density × the number of simulated bolt holes 21 × y ÷ z = the true material density of the bolt 4 × j;

[0081] Set the calculated interference amount at the contact position between the bolt 4 and the simulated bolt hole 21, and the geometric interference amount of the calculation model + the calculated interference amount = the geometric interference amount of the true model;

[0082] Among them, the geometric interference of the calculation model = the radius of bolt 4 of the calculation model - the radius of the bolt hole of the calculation model, the geometric interference of the real model = the actual radius of bolt 4 - the actual radius of the bolt hole, j is the number of bolt holes of the original geometric model 1 of the blade disk, y is the number of blades of the original geometric model 1 of the blade disk, and z is the number of blades of the calculation model. It should be noted that since the bolt hole diameter is adjusted, the calculated interference needs to be adjusted so that the total interference (calculated interference + geometric interference) is equal to the original total interference to ensure the calculation accuracy.

[0083] Step S50, perform finite element calculation according to the parameters of the adjusted bolt 4 and the simulated bolt hole 21 and the preset parameters. By adjusting the relevant parameters, the result accuracy of the finite element calculation is affected. Specifically, the interference between the calculated bolt 4 and the simulated bolt hole 21 and the density of the bolt 4 are adjusted, thereby affecting the accuracy of the calculation result. The preset parameters include load parameters such as rotational speed and temperature.

[0084] In one embodiment, after the simulated bolt hole 21 is opened on the cyclic symmetric model, it includes:

[0085] Select the original relief groove 30 on the original geometric model 1 of the blade disk, and generate a sheet body 31 at the position of the original relief groove 30;

[0086] Copy the sheet body 31 to the calculation model and form an offset sheet body 32 by offsetting a preset distance;

[0087] Rotate and copy the offset sheet body 32 around the axis of the original geometric model 1 of the blade disk and cut the calculation model to obtain the target number of simulated relief grooves 3. Specifically, in this step, select a relief groove on the original geometric model 1 of the blade disk (the number of relief grooves on each mounting edge of the original geometric model (if any), the number of relief grooves generally equals the number of bolt holes corresponding to the mounting on the original geometric model 1 of the blade disk, that is, the number of bolts 4 on the mounting edge), and generate a sheet body 31 at the position of this relief groove. Copy this sheet body 31 to Figure 6 the shown calculation model, and generate an offset sheet body 32 by offsetting a predetermined distance. Perform rotational copy (associated copy) of the offset sheet body 32 around the axis (rotor axis) of the original geometric model 1 of the blade disk to obtain the target number of simulated relief grooves 3, and the number of copied copies equals the number of bolts 4 in the calculation model.

[0088] In one embodiment, rotating and copying the offset sheet body 32 around the axis of the original geometric model 1 of the blade disk and cutting the calculation model to obtain the target number of simulated relief grooves 3 includes:

[0089] The rotation angle between two adjacent sheet bodies 31 = 360°÷y×z÷p, where p is the target number, y is the number of blades of the original geometric model 1 of the blade disk, and z is the number of blades of the calculation model;

[0090] The volume of each unloading groove of the calculation model × p × y ÷ z = the volume of the original unloading groove 30 × the number of the original unloading grooves 30. It should be noted that, as Figure 5 shown, the installation edge is trimmed by using the unloading grooves generated after offset and the unloading grooves generated by associated copying to generate a predetermined number of unloading grooves of the calculation model (the two incomplete unloading grooves on both sides form 1 complete unloading groove, and the number of unloading grooves of the calculation model is equal to the number of bolts 4 of the calculation model), where the target number is equal to the number of bolts of the calculation model and also equal to the number of bolt holes of the calculation model.

[0091] The present invention was verified by using a general finite element analysis software. Taking the blade disc shown in Figure 1 as an example, the finite element analysis considered the following 5 calculation models, as shown in Table 2.

[0092] Table 2

[0093]

[0094]

[0095] The corresponding displacement and stress calculation results obtained from the finite elements calculated by each model are shown in Table 3 below. It can be seen that: ① The disadvantage of Model II is that it cannot calculate the displacement and stress of the blade; ② The disadvantages of Model III and Model IV are that there are large errors in the displacement and stress calculation results; ③ The displacement and stress calculation errors of Model V are both small, and the absolute value of the relative error is within 0.5%.

[0096] Table 3

[0097]

[0098] In summary, by implementing this method, the structure can be simplified, and the calculation model can be made into a cyclic symmetric model containing only 1 bolt 4 and 1 blade, or 1 bolt 4 and multiple blades, or multiple bolts 4 and 1 blade. By adjusting the size of the bolt holes of the calculation model, the size of the unloading grooves, the density of the bolts 4, and the calculated interference amount between the bolts 4 and the bolt holes, the mass of each part of the calculation model = the mass of each part of the actual model; the calculation model can consider the clearance between the bolts 4 and the bolt holes.

[0099] The present invention also provides a wheel disc finite element calculation device based on a finite element calculation model, and a wheel disc finite element calculation method based on the finite element calculation model, including:

[0100] A determination module, configured to determine the number of bolts 4 and the number of blades of the calculation model according to the number of bolts 4 and the number of blades of the original geometric model of the blade disc;

[0101] An acquisition module, configured to acquire a cyclic symmetric model according to the number of blades of a calculation model, where the cyclic symmetric model is cut from the original geometric model 1 of a blisk after deleting the original bolt holes;

[0102] An opening and assembling module, configured to open simulated bolt holes 21 on the cyclic symmetric model and dispose bolts 4 in the simulated bolt holes 21;

[0103] A pre-calculation processing module, configured to import the cyclic symmetric model into finite element analysis software to adjust the parameters of the bolts 4 and the simulated bolt holes 21;

[0104] A calculation module, configured to perform finite element calculation according to the adjusted parameters of the bolts 4 and the simulated bolt holes 21 and preset parameters.

[0105] This device can determine the bolts 4 and the number of blades of the calculation model, so as to determine the number of blades on the cyclic symmetric model. Compared with the traditional complete model or 1 / 2 or 1 / 3 cyclic symmetric model, the structure of the cyclic symmetric model is simpler, reducing the amount of calculation, and thus improving the efficiency of finite element calculation. In addition, the convergence of finite element calculation gradually deteriorates as the number of contact pairs increases. In the present invention, by adjusting the parameters of the bolts 4 and the simulated bolt holes 21, the number of simulated bolt holes 21 and the number of contact pairs formed by the corresponding bolts 4 can also be adjusted. Since the structure of the cyclic symmetric model is simpler, the number of contact pairs is correspondingly reduced, thereby improving the calculation convergence. It has been verified that the calculation results of the stress value and displacement value of the blisk obtained by performing finite element calculation according to the adjusted parameters of the bolts 4 and the simulated bolt holes 21 and the preset parameters are more accurate.

[0106] The present invention also provides a storage medium, including: a disk finite element calculation method based on a finite element calculation model implemented when a program is processed and executed. By this storage medium, the disk finite element calculation method based on the finite element calculation model is run, and thus the calculation accuracy is improved through the disk finite element calculation method based on the finite element calculation model.

[0107] Although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that: they can still modify the technical solutions recorded in the foregoing embodiments, or perform equivalent substitution on some of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A finite element calculation method for a roulette based on a finite element calculation model, characterized in that, The method includes: Determining the number of bolts and the number of blades of the calculation model according to the number of bolts and the number of blades of the original geometric model of the blisk; Obtaining a cyclic symmetric model according to the number of blades of the calculation model, the cyclic symmetric model being cut from the original geometric model of the blisk after deleting the original bolt holes; Opening simulated bolt holes on the cyclic symmetric model and arranging the bolts in the simulated bolt holes; Importing the cyclic symmetric model into finite element analysis software to adjust the parameters of the bolts and the simulated bolt holes; Performing finite element calculation according to the adjusted parameters of the bolts and the simulated bolt holes and preset parameters.

2. The finite element calculation method of the roulette based on the finite element calculation model according to claim 1, characterized in that, The determining the number of bolts and the number of blades of the calculation model according to the number of bolts and the number of blades of the original geometric model of the blisk includes: When the ratio of the number of bolts to the number of blades of the original geometric model of the blisk is within the first set value range, the calculation model contains 1 bolt and 1 blade; or When the ratio of the number of bolts to the number of blades of the original geometric model of the blisk is less than the minimum value of the first set value range, the calculation model contains 1 bolt and 1 blade or multiple blades, and the number of blades of the calculation model is the largest divisor not exceeding the preset value, the preset value being the ratio of the number of blades to the number of bolts of the original geometric model of the blisk; or When the ratio of the number of bolts to the number of blades of the original geometric model of the blisk is greater than the maximum value of the first set value range, the calculation model contains multiple bolts and 1 blade, and the number of bolts of the calculation model is the ratio of the number of bolts to the number of blades of the original geometric model of the blisk, and if this ratio is not an integer, it is rounded to an integer.

3. The finite element calculation method of the roulette based on the finite element calculation model according to claim 1, characterized in that The cyclic symmetric model being cut from the original geometric model of the blisk after deleting the original bolt holes includes: Forming a first surface according to the axis of the original geometric model of the blisk and a line preset on the rim of the original geometric model of the blisk; Rotating the first surface around the axis of the original geometric model of the blisk by a preset angle to form a second surface; Cutting the original geometric model of the blisk and taking the sector model between the first surface and the second surface as the cyclic symmetric model.

4. The finite element calculation method of the roulette based on the finite element calculation model according to claim 3, characterized in that, The rotating the first surface around the axis of the original geometric model of the blisk by a preset angle to form a second surface includes: The preset angle is 360°÷y×z; y is the number of blades of the original geometric model of the blisk, and z is the number of blades of the calculation model.

5. The finite element calculation method of the roulette based on the finite element calculation model according to claim 1, characterized in that Opening simulated bolt holes on the cyclic symmetric model includes: Adding simulated bolt holes at the circumferential middle position of the mounting edge of the cyclic symmetric model, wherein the position of the simulated bolt holes on the cyclic symmetric model is determined according to the radius of the center of the simulated bolt holes being equal to the radius of the center of the original bolt holes; The volume of the simulated bolt holes × the number of simulated bolt holes × y÷z = the volume of the corresponding original bolt holes of the original geometric model of the blisk × the number of original bolt holes, y is the number of blades of the original geometric model of the blisk, and z is the number of blades of the calculation model.

6. The finite element calculation method of the roulette based on the finite element calculation model according to claim 1, characterized in that The importing the cyclic symmetric model into finite element analysis software to adjust the parameters of the bolts and the simulated bolt holes includes: Modifying the density of the bolt material to obtain the target density of the bolt material; The target density × the number of simulated bolt holes × y÷z = the true material density of the bolt × j; Set the calculated interference amount at the contact position between the bolt and the simulated bolt hole, and the geometric interference amount of the calculation model + the calculated interference amount = the geometric interference amount of the real model; Among them, the geometric interference amount of the calculation model = the radius of the bolt in the calculation model - the radius of the bolt hole in the calculation model, the geometric interference amount of the real model = the actual bolt radius - the actual bolt hole radius, j is the number of bolt holes in the original geometric model of the blade disk, y is the number of blades in the original geometric model of the blade disk, and z is the number of blades in the calculation model.

7. The finite element calculation method of the roulette based on the finite element calculation model according to claim 1, wherein, After opening the simulated bolt hole on the cyclic symmetric model, it includes: Select the original relief groove on the original geometric model of the blade disk, and generate a sheet body at the position of the original relief groove; Copy the sheet body into the calculation model and form a biased sheet body by offsetting a preset distance; Rotate and copy the biased sheet body around the axis of the original geometric model of the blade disk and cut the calculation model to obtain the target number of simulated relief grooves.

8. The finite element calculation method of the roulette based on the finite element calculation model according to claim 7, characterized in that The step of rotating and copying the biased sheet body around the axis of the original geometric model of the blade disk and cutting the calculation model to obtain the target number of simulated relief grooves includes: The rotation angle between two adjacent sheet bodies = 360°÷y×z÷p, where p is the target number, y is the number of blades in the original geometric model of the blade disk, and z is the number of blades in the calculation model; The volume of each relief groove in the calculation model × p × y÷z = the volume of the original relief groove × the number of original relief grooves, where the target number is equal to the number of bolts in the calculation model.

9. A finite element calculation device for a roulette wheel based on a finite element calculation model, based on the finite element calculation method for a roulette wheel based on a finite element calculation model according to any one of claims 1-8, characterized in that, It includes: A determination module for determining the number of bolts and the number of blades in the calculation model according to the number of bolts and the number of blades in the original geometric model of the blade disk; An acquisition module for acquiring the cyclic symmetric model according to the number of blades in the calculation model, and the cyclic symmetric model is cut from the original geometric model of the blade disk after deleting the original bolt holes; An opening and assembling module for opening a simulated bolt hole on the cyclic symmetric model and arranging the bolt in the simulated bolt hole; A pre-processing module for importing the cyclic symmetric model into finite element analysis software to adjust the parameters of the bolt and the simulated bolt hole; A calculation module for performing finite element calculation according to the adjusted parameters of the bolt and the simulated bolt hole and preset parameters.

10. A storage medium, characterized in that, It includes: When the program is processed and executed, it implements the disk finite element calculation method based on the finite element calculation model as described in any one of claims 1-8.

Images