CN114139278A

CN114139278A - Intelligent optimization design method for steel wheel spoke process curve

Info

Publication number: CN114139278A
Application number: CN202111322169.5A
Authority: CN
Inventors: 袁飞; 黄信; 王贤付; 龚明学
Original assignee: Dongfeng Automotive Wheel Suizhou Co ltd
Current assignee: Dongfeng Automotive Wheel Suizhou Co ltd
Priority date: 2021-11-09
Filing date: 2021-11-09
Publication date: 2022-03-04

Abstract

The invention relates to an intelligent optimization design method for a steel wheel spoke process curve, which comprises the following main steps: s1, carrying out three-dimensional modeling on a drawing die and a reverse drawing die of the spoke, parameterizing key structures of the drawing die and the reverse drawing die, compiling a macro program, and driving the macro program through a specified command to realize automatic modification and updating of the drawing die and the reverse drawing die; s2, carrying out punching process simulation analysis on the drawing die and the reverse drawing die by using finite element analysis software; s3, carrying out secondary development on the relevant programs of the deep drawing and reverse deep drawing stamping procedures, and replacing the program of the original finite element analysis software with the program of the secondary development; s4, performing high-precision approximate model fitting on the simulation model; s5, building an approximate model optimization design flow frame in ISIGHT software, and performing optimization design on the created approximate model, so that the limitation of the existing process curve optimization design is broken through, and the optimal deep drawing and reverse deep drawing process curves can be quickly found.

Description

Intelligent optimization design method for steel wheel spoke process curve

Technical Field

The invention relates to the technical field of automobile wheel spoke design, in particular to an intelligent optimization design method for a steel wheel spoke process curve.

Background

The wheel is an important part of an automobile, the performance of the wheel is mainly influenced by the product structure and the stamping process, when the steel wheel is developed, the product structure needs to be optimally designed, the thickness of a plate is reduced as much as possible on the premise of meeting the product performance and cost requirements, and the lightweight of the wheel is realized to the greatest extent. On the other hand, the production quality of the steel wheel is influenced by the stamping process, the phenomena of wrinkling, cracking, bulging, thinning and the like of a product can be caused due to unreasonable die process design in the production process, and the produced sample piece has larger deviation with the size of a product drawing, so that the product quality and the safety performance are influenced. The steel wheel spoke is characterized in that 9 processes are shared in the production process, wherein 2 processes of deep drawing and reverse deep drawing are important, and the process curves of the 2 processes directly determine the thinning rate, wrinkling rate, cracking rate and the like of key parts in the forming process of the spoke.

When designing the die, most domestic steel wheel enterprises mainly determine the process curves of the 2 procedures of deep drawing and reverse deep drawing through engineering development experience, the reduction rate of key parts of spokes is controlled to be less than or equal to 10% in the production process, but for some products with thinner thickness (less than or equal to 3.2mm) and larger size (more than or equal to 17 inches), cracking and wrinkling are easy to occur in the production process, and the cracking and wrinkling are difficult to solve by repeatedly repairing and testing the die by the engineering experience at this time, although part of the wheel enterprises have the capability of stamping process simulation, the process curves can be guided and modified according to stamping simulation results, so that the problems of repeatedly repairing and testing the die are avoided, the optimal process curves are required to be repeatedly modified and repeatedly CAE modeling and analysis are required to be found, the efficiency is low, and the product development cycle is too long to be rapidly put on the market, the development progress of the project is influenced, and the development cost is increased.

Although few enterprises realize automatic optimization of process curves of software, most of the optimization design systems of the enterprises acquire product information from commercial CAD/CAE systems through interfaces, the efficiency is greatly improved, but the optimization systems and the design systems are mutually independent, and the performance indexes cannot be comprehensively considered to carry out multidisciplinary optimization design, so that the optimization effect is limited.

How to break through the limitation of the prior process curve optimization design is to research an intelligent optimization design method for the process curve of the steel wheel spoke, and the key of the invention is to quickly find the optimal deep drawing and reverse deep drawing process curves.

Disclosure of Invention

Based on the above expression, the invention provides an intelligent optimization design method for a steel wheel spoke process curve, which aims to solve the technical problems that the steel wheel spoke optimization design in the prior art is large in limitation and poor in comprehensive consideration effect of performance indexes.

The technical scheme for solving the technical problems is as follows:

an intelligent optimization design method for a steel wheel spoke process curve comprises the following main steps:

s1, carrying out three-dimensional modeling on a drawing die and a reverse drawing die of the spoke, parameterizing key structures of the drawing die and the reverse drawing die, compiling a macro program, and driving the macro program through a specified command to realize automatic modification and updating of the drawing die and the reverse drawing die;

s2, carrying out punching process simulation analysis on the drawing die and the reverse drawing die by using finite element analysis software;

s3, carrying out secondary development on the relevant programs of the deep drawing and reverse deep drawing stamping procedures, and replacing the program of the original finite element analysis software with the program of the secondary development;

s4, performing high-precision approximate model fitting on the simulation model;

s5, building an approximate model optimization design flow framework in the ISIGHT software, and carrying out optimization design on the created approximate model.

Compared with the prior art, the technical scheme of the application has the following beneficial technical effects:

the intelligent optimization design method integrates the 3D design of the die, the simulation of the stamping process and the optimization of the process curve on the basis of the secondary development of a computer program, realizes the integration automation of three modules of the die design, the simulation of the stamping process and the process optimization, integrates the scientific technologies of various subjects, comprehensively and systematically analyzes from the overall situation, quickly finds the optimal stretching and reverse stretching process curve of the spoke, and prevents the spoke from wrinkling, bulging, serious thinning and other process problems in the forming process.

On the basis of the technical scheme, the invention can be further improved as follows.

Further, in the main step S1, the CATIA software is used to perform three-dimensional modeling, and a CATIA macro program is written, where the designated command is a bat command.

Further, the main step S2 includes the following sub-steps:

s21, mould introduction and plate creation: importing the parameterized numerical model and the formed numerical model of the drawing die and the reverse drawing die into finite element analysis software; converting the introduced solid part into a box body and creating a plate with the same thickness as the box body;

s22, creating a new material type, defining performance parameters of the material, and giving material attributes to the plate with the defined thickness;

s23, assembling a die: assembling the die according to the die closing condition;

s24, a step of analyzing the wound parts: creating a plurality of analysis steps and defining the types of the analysis steps;

s25, defining interactions: creating interaction attributes, and defining an interaction type, a friction type and a friction coefficient;

s26, defining boundary conditions and loads: creating an amplitude curve in finite element analysis software and defining a curve type;

s27, grid division: performing grid division on the plate and the die;

s28, submitting and analyzing: carrying out finite element analysis on the plate and the die for dividing the grids and deriving a minimum thickness value T of the formed spoke_minAnd maximum stress value S_max。

Further, the finite element analysis software is ABAQUS; step S21, the introduction type of the parameterized digital model and the forming digital model of the drawing die and the reverse drawing die is 3D. In the substep S21, converting the imported entity component into a box body specifically includes converting all the imported entity components into a box body through a shell module in the Abaqus software, respectively naming each of the die components as die-draw, punch-draw, binder-draw, di-form, and punch-form, establishing a reference point for each of the die components and respectively naming the reference point as RP-die-draw, RP-punch-draw, RP-binder-draw, RP-di-draw, RP-die-form, and RP-punch-form, where the type of the created sheet is 3d.

Further, defining the performance parameters of the material in substep S22 includes defining at least the density, elastic modulus, poisson' S ratio, and plasticity of the material.

Further, the step S23 of completing the assembly of the mold according to the mold closing condition specifically includes: uniformly distributing the plate and the die at a certain distance along the X axis according to the sequence of the plate, the drawing die, the reverse drawing die and the forming die.

Further, the step S24 of creating a plurality of analysis steps and defining the types of the analysis steps specifically includes: creating 12 analysis steps, wherein the type of the analysis steps is dynamic. explore, the analysis steps respectively correspond to the opening of a drawing die, the sheet entering the drawing die, the edge pressing force application, the drawing, the opening of the drawing die and a reverse drawing die, the sheet entering the reverse drawing die, the edge pressing force application, the reverse drawing, the opening of the reverse drawing die and a forming die, the sheet entering the forming die, the forming and the forming die opening, and the 12 analysis steps are respectively named as a segment 1, a blank position1, a force1, a draw, a segment 2, a blank position2, a force2, a reverse draw, a segment 3, a blank position3, a form and a segment 4, and adding an STH Output in a Field Output Manager to obtain the thickness variation condition of the sheet after the analysis.

Further, defining the interaction type as Contact and the friction type as Penalty in the substep S25, defining the friction coefficient, defining the interaction between the drawing die and the sheet in the draw analysis step, defining the interaction between the reverse drawing die and the sheet in the reverse draw analysis step, and defining the interaction between the forming die and the sheet in the form analysis step; the type of the amplitude curve in the substep S26 is Smooth step.

Further, in the main step S3, a secondary development is performed on the program related to the drawing and reverse drawing stamping process by using Python, and the program of the original finite element analysis software is replaced by the secondary developed Python program, which includes the following sub-steps:

s31, determining the points O (O1, O2, O3), B (B1, B2, B3), then determining the portion of the annular zone where the slab is pressed by function 1, wherein,

the function 1 is:

getByBoundingCylinder(center＝(O1,O2,O3),center(B1,B2,B3),radius＝C)；

s32, determining a fixed point coordinate (x, y, z) in the middle area of the horizontal plane section of the female die, and then finding a plane passing through the point through a FindAt () function;

s33, after spoke forming, extracting three key part elements of a spoke bottom fillet, a spoke top and a radiating hole respectively through a function 2, thereby establishing set1, set2 and set3 respectively, and outputting minimum thicknesses T1 of the three key parts respectively_min、T2_min、T3_minAnd maximum stress S1_max、S2_max、S3_maxWherein, in the step (A),

the function 2 is:

getByBoundingBox(x_min＝a1,x_max＝a2,y_min＝b1,y_max＝b2,z_min＝c1,z_max＝c2)；

s34, replacing the corresponding program in the ABAQUS analysis by all the secondarily developed Python programs, debugging the Python programs, and running the replaced Python programs through run script.

Further, the main step S4 of performing high-precision approximate model fitting on the simulation model includes the following sub-steps:

s41, sampling the N key structure parameters in a design space region by an optimized Latin method, setting the number of the key structure parameters to be N, setting the number of sample data to be M, wherein M, N is positive integer, wherein M, N meets the condition that M is more than or equal to (N +1) × (N +2)/2, operating DOE after the DOE process is built, realizing automatic calculation of the process, and obtaining the minimum thickness T1 by taking the N key structure parameters as independent variables_min、T2_min、T3_minAnd maximum stress S1_max、S2_max、S3_maxThe 6 output results are M groups of space sample data of the dependent variable;

s42, evaluating the obtained sample data, eliminating points with abnormal results or carrying out calculation evaluation again to ensure that all sample points are normal; when the approximate model is established, the approximate model is fitted in a second-order response surface mode, and the expression of the second-order response surface model is as follows:

s43 by R²And average error pair T1_min、T2_min、T3_min、S1_max、S2_max、S3_maxAnd performing error analysis on the six indexes until the precision reaches a preset precision, and performing approximate model fitting again by adding certain sample point data each time under the condition that the precision of the approximate model is smaller than the preset precision.

Further, the step of building an approximate model optimization design flow box in the ISIGHT software in the main step S5 includes setting the N key structure parameters in an Approximation module, discretizing the N key structure parameters within their respective variation ranges, and iterating for 1 time at predetermined intervals;

the main step S5 of optimizing and designing the created approximate model includes selecting a global ASA algorithm and a local gradient Optimization NLPQLP algorithm in an Optimization module of ISIGHT software to perform combined calculation, so as to ensure that the result is an optimal solution; at T1_min、T2_min、T3_minMaximization, S1_max、S2_max、S3_maxAnd the minimization is to calculate an optimization target and find an optimal drawing and reverse drawing curve structure.

Drawings

Fig. 1 is a schematic step diagram of an intelligent optimized design method for a process curve of a steel wheel spoke according to an embodiment of the invention;

FIG. 2 is a schematic view of the substep of S2 in FIG. 1;

FIG. 3 is a schematic diagram of key structural parameters of a drawing die and a reverse drawing die in an embodiment of the invention;

fig. 4 is a flowchart of step S4.

Detailed Description

To facilitate an understanding of the present application, the present application will now be described more fully with reference to the accompanying drawings. Embodiments of the present application are set forth in the accompanying drawings. This application may, however, be embodied in many different forms and should not be construed as limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete.

Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs. The terminology used herein in the description of the present application is for the purpose of describing particular embodiments only and is not intended to be limiting of the application.

It will be understood that spatial relationship terms, such as "under", "below", "beneath", "below", "over", "above", and the like, may be used herein to describe one element or feature's relationship to another element or feature as illustrated in the figures. It will be understood that the spatially relative terms are intended to encompass different orientations of the device in use or operation in addition to the orientation depicted in the figures. For example, if the device in the figures is turned over, elements or features described as "below" or "beneath" other elements or features would then be oriented "above" the other elements or features. Thus, the exemplary terms "under" and "under" can encompass both an orientation of above and below. In addition, the device may also include additional orientations (e.g., rotated 90 degrees or other orientations) and the spatial descriptors used herein interpreted accordingly.

It will be understood that when an element is referred to as being "connected" to another element, it can be directly connected to the other element or be connected to the other element through intervening elements. The "connection" in the following embodiments is understood as "electrical connection", "communication connection", or the like if the connected circuits, modules, units, or the like have electrical signals or data transmission therebetween.

As used herein, the singular forms "a", "an" and "the" may include the plural forms as well, unless the context clearly indicates otherwise. It will be further understood that the terms "comprises/comprising," "includes" or "including," etc., specify the presence of stated features, integers, steps, operations, components, parts, or combinations thereof, but do not preclude the presence or addition of one or more other features, integers, steps, operations, components, parts, or combinations thereof.

An intelligent optimization design method for a steel wheel spoke process curve comprises the main steps from S1 to S5d, and specifically comprises the following steps:

the method comprises the following steps that a drawing die and a reverse drawing die of a spoke both comprise a female die and a male die, common spoke drawing and reverse drawing processes are taken as examples, the key structure parameters of the drawing die and the reverse drawing die refer to quantitative values capable of preparing and expressing key structures of a die structure, generally include but are not limited to drawing size, fillet radius size, die top width, die bottom width and the like, and the corresponding key structure sizes can be determined by technicians in the field during die manufacturing, for example, in one embodiment of the application, as shown in fig. 3, the key structure parameters of the drawing die and the reverse drawing die are 11 in total, and specifically include H1 and H4 related to the drawing sizes; r1, R2, R3, R4, R5, R6, R7 relating to the fillet radius size; h2 for die top width and H3 for die bottom width.

It can be understood that the three-dimensional software used for the three-dimensional modeling may be reasonably selected according to actual conditions, in this embodiment, the CATIA software is used for the three-dimensional modeling, a CATIA macro program is written, and the macro program is driven by a bat command to realize automatic modification and update of the stretching and reverse stretching modes.

Then, the step S2 is executed, and the finite element analysis software is used for carrying out the simulation analysis of the stamping process on the drawing die and the reverse drawing die;

similar to the above, the choice of finite element analysis software can be determined according to practical conditions and the choice of designer, and the present invention provides only an optimized design method, and does not limit the physical or computer tools used in the application, and preferably, the finite element analysis software is selected as ABAQUS.

The main step S2 can be divided into sub-steps S21 to S28 according to the procedure, and specifically, the flow sequence sequentially includes:

the method comprises the steps of establishing a reference point for each die part, and respectively naming RP-di-draw, RP-punch-draw, RP-binder-draw, RP-di-re-draw, RP-punch-form, RP-punch-re-draw, RP-binder-re-draw, RP-di-punch-draw, RP-punch-form and P-punch-form. Wherein the type of the created plate is 3D.

newly building a material type, defining the density, the elastic modulus, the Poisson ratio and the plasticity of the material, defining the thickness of a plate material, and endowing the plate material with the properties of the newly-built material.

specifically, the plate and the die are uniformly distributed at a certain distance along the X axis according to the sequence of the plate, the drawing die, the reverse drawing die and the forming die.

in this embodiment, in order to fully analyze the interaction between the sheet and the die in various states, 12 analysis steps are created, the analysis steps are of the type dynamic, wherein the 12 analysis steps correspond to the opening of a drawing die, the sheet entering the drawing die, the application of a blank pressing force, drawing, the opening of the drawing die and a reverse drawing die, the sheet entering the reverse drawing die, the application of a blank pressing force, reverse drawing, the opening of the reverse drawing die and a forming die, the sheet entering the forming die, forming and the opening of the forming die, and the 12 analysis steps are named as a blank 1, a blank position1, a force1, a draw, a blank 2, a blank position2, a force2, a reverse draw, a blank 3, a blank position3, a form and a blank 4 respectively.

And adding STH Output in the Field Output Manager in order to obtain the thickness change condition of the plate after analysis.

defining the interaction type as Contact, the friction type as Penalty, defining the friction coefficient, defining the interaction between the drawing die and the sheet in a draw analysis step, defining the interaction between the reverse drawing die and the sheet in a reverse draw analysis step, and defining the interaction between the forming die and the sheet in a form analysis step on the basis of the analysis steps.

When the male die surface and the female die surface are selected, more selection is required, otherwise, the strain of the plate material in the analysis process is overlarge, so that the analysis is not converged, less selection is not required, the plate material can penetrate through the die in the analysis process, and the result is inconsistent with the actual result. A mass point and moment of inertia are defined for each mold part, the mass point selecting the reference point established in substep S21.

S26, defining boundary conditions and loads: and creating an amplitude curve and defining a curve type in the finite element analysis software, wherein the amplitude curve type is Smooth step in the step S26.

Taking drawing processes corresponding to the former six analysis steps as an example, in the separate1 analysis step, a Z-direction displacement is applied to the punch-draw, the remaining degrees of freedom are limited, and the mold is opened; applying X-direction displacement to the plate in a blank position1 analysis step, limiting the other degrees of freedom, and putting the plate into a drawing die; releasing displacement limitation in the binderZ direction in the force analysis step and applying a blank holder force in the Z direction, wherein the load point is RP-binder-draw; applying a displacement to the punch in the draw analysis step opposite to that in the step of seperate1 analysis to complete drawing; in addition, in order to prevent rigid displacement of the plate during stamping, the displacement of the center of the plate in the X and Y directions is limited; in the step of seperate2 analysis, displacement is applied to two sets of dies, and a drawing die and a reverse drawing die are opened simultaneously; in the blank position2 analysis step, the sheet material is displaced so that it is taken out of the drawing die and placed in the reverse drawing die.

According to the set logic, the subsequent reverse drawing, die opening and sheet material placing into a forming die are similar to the drawing setting; note that no blank holder is arranged during molding, and blank holder force is not required to be applied, so long as displacement is applied to the mold to realize molding.

S27, grid division: performing grid division on the plate and the die;

it should be noted that the sheet belongs to a deformable body, and the stress-strain condition in the sheet forming process needs to be closely concerned, so the grid needs to be divided densely, and the grid type adopts S4R shell units.

The die belongs to a discrete rigid body and does not deform in the punch forming process, so that a grid with a thick point can be divided to accelerate the calculation speed, and the grid type of the die is an R3D4 rigid body shell unit.

In the analysis and optimization process, as the 3D structure of the stretching and reverse stretching die changes along with the change of key structure parameters, the positions and the numbers of finite element elements of the stretching and reverse stretching male and female dies change along with the introduction of a 3D model every time when CAE analysis is carried out, if a py program automatically generated by ABAQUS software in the whole CAE analysis process is directly applied, the optimization process is wrong, and therefore, the CAE program related to the stretching and reverse stretching stamping processes is secondarily developed by adopting a programming means.

Then, the main step S3 is performed to perform secondary development on the program related to the drawing and reverse drawing stamping process, and the program of the original finite element analysis software is replaced by the program of the secondary development, in this embodiment, the program related to the drawing and reverse drawing stamping process is secondarily developed by Python, and the program of the original finite element analysis software is replaced by the program of the secondary development Python, specifically, the main step S3 includes three substeps according to the step flow, namely, S31 to S33, wherein,

because before punching, the sheet material can be pushed down earlier with the blank holder to tensile, the anti-drawing die horizontal segment, and the die horizontal segment establishes the contact with the sheet material upper surface and is right the relation, prevents to appear wrinkling phenomenon at the in-process of punching press, in order to guarantee how the die structure changes, the contact is right the relation and is unchangeable all the time, carries out substep S31 and S32, promptly: s31, determining the points O (O1, O2, O3), B (B1, B2, B3), then determining the portion of the annular zone where the slab is pressed by function 1, wherein,

the function 1 is:

getByBoundingCylinder(center＝(O1,O2,O3),center(B1,B2,B3),radius＝C)；

s32, by determining a fixed point coordinate (x, y, z) in the middle area part of the horizontal plane segment of the die, and then finding a plane passing through the point by the FindAt () function,

this ensures that the constraint pair relationship does not change all the time.

After the spoke is formed, it is usually necessary to pay attention to the reduced thickness and maximum stress of the spoke base fillet, the spoke top and the heat dissipation hole region, so S33 is performed, namely:

s33, extracting three key part elements of a spoke bottom fillet, a spoke top and a radiating hole respectively through a function 2, establishing set1, set2 and set3 respectively, and outputting minimum thicknesses T1 of the three key parts respectively_min、T2_min、T3_minAnd maximum stress S1_max、S2_max、S3_maxWherein

The function 2 is:

care is taken to ensure that each output is the same part of thickness and stress.

After the program is developed for 2 times, the program developed for the second time needs to be replaced with the original corresponding program, that is, S34, all the programs developed for the second time are replaced with the corresponding programs in the ABAQUS analysis, the Python programs are debugged, and the replaced Python programs are run through run script.

Because the steel wheel stamping forming simulation adopts display dynamics simulation, the calculation speed is very slow, the calculation of a general simulation model needs about 3 hours, although the 3D design and stamping simulation integration automation of the die is realized by integrating CAD and CAE software by using the light software through the python secondary development technology, the complex work of manual repeated modification, repeated modeling and the like is avoided, the whole optimization period cannot be accepted due to the slow calculation speed of the simulation model, and in order to solve the problem, the invention introduces a main step S4 and performs high-precision approximate model fitting on the simulation model; and a high-precision approximate model is used for replacing a simulation model for calculation, so that the optimization calculation time is greatly shortened.

Before the approximate model is created, the DOE design needs to be performed on simulation analysis of three procedures of wheel stretching, counter-stretching and shaping, and the method specifically comprises the following steps:

s41, the optimized latin square samples the N key structural parameters in the design space region, the number of the key structural parameters is set to N (11 in this embodiment), the number of the sample data is M, M, N is positive integers, M, N satisfies that M ≧ N +1 (N +2)/2 (78), the DOE is operated after the DOE process is completed, the process automation calculation is implemented, and the minimum thickness T1 with the 11 key structural parameters as arguments is obtained_min、T2_min、T3_minAnd maximum stress S1_max、S2_max、S3_maxAnd the 6 output results are M groups of space sample data of the dependent variable.

Then, the step S42 is executed, the obtained sample data is evaluated, points with abnormal results are removed or calculation and evaluation are carried out again, and all sample points are ensured to be normal; when the approximate model is created, the approximate model is fitted in a second-order response surface mode, because the second-order RSM (response surface model) has better effect when approximating a nonlinear and arbitrary design space, wherein the expression of the second-order response surface model is as follows:

then, the approximate model obtained by fitting is subjected to error analysis, and step S43 is executed to pass through R²And average error pair T1_min、T2_min、T3_min、S1_max、S2_max、S3_maxAnd performing error analysis on the six indexes until the precision reaches a preset precision, and performing approximate model fitting again by adding certain sample point data each time under the condition that the precision of the approximate model is smaller than the preset precision.

In this embodiment, the predetermined accuracy is not less than 95%, and when the accuracy of the approximate model is less than 95%, 1/5 sample point data is added each time to perform approximate model fitting again; when carrying out approximate model error analysis, the method passes through R²And analyzing six indexes of T1min, T2min, T3min, S1max, S2max and S3max according to the average error, wherein the precision of each index is more than 95%.

Step S4 can be summarized as according to the flow diagram, data is collected first, an approximate model type is selected, then an approximate model is obtained according to the collected data fitting, error analysis is performed on the approximate model, if the model meets a predetermined precision requirement, the model can be accepted, that is, the approximate model is used, if the model does not meet the predetermined precision requirement, a certain sample point data is added, and the above steps are repeated until an acceptable approximate model is obtained.

And after the fitting approximate model meeting the requirements is obtained, executing a main step S5, building an approximate model optimization design flow framework in the ISIGHT software, and performing optimization design on the created approximate model.

Specifically, the step of building an approximate model optimization design flow box in the ISIGHT software in the main step S5 includes setting the 11 key structure parameters in an Approximation module, discretizing the N key structure parameters within their respective variation ranges, and iterating for 1 time at predetermined intervals; in this embodiment, the predetermined interval is preferably 0.1.

The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like that fall within the spirit and principle of the present invention are intended to be included therein.

Claims

1. An intelligent optimization design method for a steel wheel spoke process curve is characterized by comprising the following main steps:

2. The intelligent optimized design method for the technological curve of the steel wheel spoke as claimed in claim 1, wherein in the main step S1, three-dimensional modeling is performed by adopting CATIA software, and a CATIA macro program is written, and the specified command is a bat command.

3. The intelligent optimization design method for the steel wheel spoke process curve according to claim 1, wherein the main step S2 comprises the following sub-steps:

s27, grid division: performing grid division on the plate and the die;

4. The intelligent optimization design method for the steel wheel spoke process curve according to claim 3, wherein the finite element analysis software is ABAQUS; step S21, the introduction type of the parameterized digital model and the forming digital model of the drawing die and the reverse drawing die is 3D. In the substep S21, converting the imported entity component into a box body specifically includes converting all the imported entity components into a box body through a shell module in the Abaqus software, respectively naming each of the die components as die-draw, punch-draw, binder-draw, di-form, and punch-form, establishing a reference point for each of the die components and respectively naming the reference point as RP-die-draw, RP-punch-draw, RP-binder-draw, RP-di-draw, RP-die-form, and RP-punch-form, where the type of the created sheet is 3d.

5. The intelligent optimized design method for the steel wheel spoke process curve according to claim 4, wherein the defining of the performance parameters of the material in the substep S22 at least comprises defining the density, the elastic modulus, the Poisson' S ratio and the plasticity of the material.

6. The intelligent optimization design method for the steel wheel spoke process curve according to claim 4, wherein the step S23 of completing the assembly of the mold according to the mold closing condition specifically comprises the following steps: uniformly distributing the plate and the die at a certain distance along the X axis according to the sequence of the plate, the drawing die, the reverse drawing die and the forming die.

7. The intelligent optimization design method for the steel wheel spoke process curve according to claim 4, wherein the step S24 of creating a plurality of analysis steps and defining the types of the analysis steps specifically comprises: creating 12 analysis steps, wherein the type of the analysis steps is dynamic. explore, the analysis steps respectively correspond to the opening of a drawing die, the sheet entering the drawing die, the edge pressing force application, the drawing, the opening of the drawing die and a reverse drawing die, the sheet entering the reverse drawing die, the edge pressing force application, the reverse drawing, the opening of the reverse drawing die and a forming die, the sheet entering the forming die, the forming and the forming die opening, and the 12 analysis steps are respectively named as a segment 1, a blank position1, a force1, a draw, a segment 2, a blank position2, a force2, a reverse draw, a segment 3, a blank position3, a form and a segment 4, and adding an STH Output in a Field Output Manager to obtain the thickness variation condition of the sheet after the analysis.

8. The intelligent optimization design method for the steel wheel spoke process curve according to claim 6, wherein the interaction type is defined as Contact in the substep S25, the friction type is Penalty, the coefficient of friction is defined, the interaction between the drawing die and the sheet is defined in the draw analysis step, the interaction between the reverse drawing die and the sheet is defined in the reverse draw analysis step, and the interaction between the forming die and the sheet is defined in the form analysis step; the type of the amplitude curve in the substep S26 is Smooth step.

9. The intelligent optimization design method for the process curve of the steel wheel spoke according to the claim 3 is characterized in that in the main step S3, Python is adopted to carry out secondary development on related programs of drawing and reverse drawing stamping procedures, and the secondarily developed Python program is used for replacing the program of the original finite element analysis software, and the method comprises the following steps:

determining the points O (O1, O2, O3) and B (B1, B2 and B3), and determining the pressed annular area part of the plate material through a function 1, wherein

The function 1 is:

getByBoundingCylinder(center＝(O1,O2,O3),center(B1,B2,B3),radius＝C)；

the function 2 is:

s34, replacing corresponding programs in the ABAQUS analysis by all the secondarily developed Python programs, debugging the Python programs, and running the replaced Python programs through run script;

the main step S4 of performing high-precision approximate model fitting on the simulation model includes the following sub-steps:

10. The intelligent optimized design method for the process curve of the steel wheel spoke according to claim 9, wherein the step of constructing an approximate model optimized design flow box in the ISIGHT software in the main step S5 includes setting the N key structure parameters in an Approximation module, discretizing the N key structure parameters in respective variation ranges, and iterating for 1 time at predetermined intervals;