CN105825025A - Plate-shaped part processing deformation simulation prediction method - Google Patents
Plate-shaped part processing deformation simulation prediction method Download PDFInfo
- Publication number
- CN105825025A CN105825025A CN201610191266.8A CN201610191266A CN105825025A CN 105825025 A CN105825025 A CN 105825025A CN 201610191266 A CN201610191266 A CN 201610191266A CN 105825025 A CN105825025 A CN 105825025A
- Authority
- CN
- China
- Prior art keywords
- plate part
- cutting
- deformation
- model
- depth
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23Q—DETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
- B23Q15/00—Automatic control or regulation of feed movement, cutting velocity or position of tool or work
- B23Q15/20—Automatic control or regulation of feed movement, cutting velocity or position of tool or work before or after the tool acts upon the workpiece
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mechanical Engineering (AREA)
- Numerical Control (AREA)
Abstract
The invention discloses a plate-shaped part processing deformation simulation prediction method, which comprises the following steps of calculating the surface density of equivalent acting forces required for causing deformation; determining a plurality of groups of cutting parameters including a plurality of practical cutting depth values; performing cutting residue stress simulation analysis according to the plurality of groups of cutting parameters; performing fitting on the residue stress values of a plurality of points so as to obtain a function relationship; performing grid division on a model of the plate-shaped part so as to obtain a plurality of units; determining a residue stress state of each unit; loading the cutting residue stress field and distribution force load; submitting a limit element analysis; obtaining a deformation result; if the deformation result does not meet the requirements, changing the cutting parameters and performing the former steps until the deformation result meets the requirements. The plate-shaped part processing deformation simulation prediction method according to the embodiment of the invention has the advantage that the plate-shaped part can effectively and precisely realize a deformation prediction.
Description
Technical field
The present invention relates to mechanical field, in particular to plate part machining deformation simulated prediction method.
Background technology
Along with the development of aerospace industry, in order to alleviate weight, promote mobility etc., thin-walled parts has been obtained for increasingly being widely applied.But the weak rigidity characteris of thin-walled parts also makes it be easily deformed in the course of processing, affect follow-up assembling and use link.Elongated sheet part is the most representational thin-walled parts, with it as basic element, may be constructed the integral structure component of complexity, complicated blade construction etc..
In the course of processing of elongated sheet part, in order to realize the indexs such as preferable flatness, generally require two sides reprocessabilty, rely on very much the experience of technical staff.Existing school shape means are also easy to adversely affect the intensity etc. of elongated sheet part.Therefore, there is the problem that milling deformation is relatively big and is difficult to predict and control in elongated sheet part.
Summary of the invention
The application is to make following facts and the discovery of problem and understanding based on inventor: the release of plate part flatness error before processing, the deformation that clamping causes and cutting residual stress all can affect final machining deformation.If the machining deformation of plate part reliably can be predicted, just further it can be carried out actively control.It is to say, initial deformation and the impact of clamping effect may be compensated for by controlling cutting residual stress.
It is contemplated that one of technical problem solved the most to a certain extent in correlation technique.To this end, the present invention proposes a kind of plate part machining deformation simulated prediction method.
Plate part machining deformation simulated prediction method according to embodiments of the present invention comprises the following steps: A) beat the upper surface profile of plate part described in table record；B) plate part described in clamping, then beats the upper surface profile of plate part described in table record, calculates the deformation caused because of described clamping, and then calculates the surface density causing the equivalent action power needed for described deformation；C) the actual milling depth of described each position of plate part is determined according to the upper surface profile of the described plate part after clamping and name milling depth, then according to the span of described actual milling depth, determining multiple real cutting depth numerical value, multiple described real cutting depth numerical value are constituted with the speed of mainshaft, feed engagement and axial milling depth organizes cutting parameter more；D) simulation analysis utilizing finite element analysis software to carry out cutting residual stress according to described many group cutting parameters, simulation analysis extracts the residual stress numerical value of multiple point after completing, be fitted the residual stress numerical value of the plurality of point to obtain the functional relationship of multiple components of stress and the distance of the finished surface apart from described plate part；E) in finite element analysis software, set up the model of described plate part, to described model partition grid to obtain multiple unit, the actual cut degree of depth of the different length position according to described model, select corresponding described functional relationship, determine its corresponding residual stress state according to the distance of the upper surface of model described in each described cell distance；F) load cutting residual stress field, and apply the distributed force load equal, in opposite direction with the surface density numerical value of described equivalent action power, submission finite element analysis at the lower surface of described model, and obtain deformation result；And G) if described deformation result meets requirement, then select described step C) in cutting parameter for processing cutting parameter, if described deformation result is unsatisfactory for requirement, then change described step C) in cutting parameter and repeat described D) to described step F), until described deformation result meets requirement.
Plate part machining deformation simulated prediction method according to embodiments of the present invention is by implementing based on working angles finite element simulation and cutting residual stress field and the loading of distributed force load, such that it is able to effectively, accurately to plate part carry out Deformation Prediction.
It addition, plate part machining deformation simulated prediction method according to the above embodiment of the present invention can also have a following additional technical characteristic:
According to one embodiment of present invention, the length of described plate part is more than or equal to 5:1 more than or equal to 3:1, the width of described plate part with the ratio of the thickness of described plate part with the ratio of the width of described plate part.
According to one embodiment of present invention, described plate part is ferromagnetic part, described plate part magnetic suction disc clamping, the processing mode of described plate part is face milling, described face milling is the milling cutter length direction feeding along described plate part, uses the mode of symmetrical milling, and wherein the diameter of milling cutter is more than the width of described plate part, advantageously, the diameter of described milling cutter and the ratio of the width of described plate part are more than or equal to 3:1.
According to one embodiment of present invention, in described step A) in, described plate part is placed on magnetic suction disc, beat the upper surface profile of plate part described in table record, described step B) comprise the following steps: B-1) open described magnetic suction disc, so that plate part described in clamping, after clamping is stable, beat the upper surface profile of plate part described in table record；B-2) according to the upper surface profile of plate part described before and after clamping, the deformation caused because of described clamping is calculated；And B-3) described plate part is reduced to simply supported beam model, a length of l of described plate part, width are b, thickness is h, elasticity modulus of materials is E, then the line density of equivalent action powerThe surface density of equivalent action power
According to one embodiment of present invention, described step C) comprise the following steps: C-1) determine the speed of mainshaft, feed engagement and axial milling depth；C-2) the actual milling depth of described plate part each position along its length is determined according to the upper surface profile of the described plate part after clamping and name milling depth；C-3) according to the span of described actual milling depth, determine multiple real cutting depth numerical value, multiple described real cutting depth numerical value are constituted with the fixing speed of mainshaft, feed engagement and axial milling depth organizes cutting parameter more, preferably, determine 5-10 real cutting depth numerical value, and then constitute 5-10 group cutting parameter.
According to one embodiment of present invention, described step D) comprise the following steps: D-1) utilize threedimensional model and the threedimensional model of described plate part of the blade of CAD software drafting Integral milling cutter or indexable milling cutter, intercept part that the actual processing of described blade uses and the part that the actual processing of described plate part is used, import in finite element analysis software；D-2) material properties of cutter, contact attribute and cutting parameter and the material properties of described plate part are set, carry out three-dimensional RAP material；D-3) after described three-dimensional RAP material terminates, load removal, extract finished surface and the residual stress numerical value of following multiple points of described plate part, the spacing of adjacent two described points is 15 microns-25 microns, preferably, the spacing of adjacent two described points is 20 microns, is fitted the residual stress numerical value of the plurality of point to obtain the functional relationship of multiple components of stress and the distance of the finished surface apart from described plate part；And D-4) change the cutting parameter in phantom and repeat described three-dimensional RAP material, in order to complete the described simulation analysis organizing cutting parameter more.
According to one embodiment of present invention, utilizing 6 rank multinomials to be fitted the plurality of point, obtain the functional relationship of multiple components of stress and the distance of the finished surface apart from described plate part, described functional relationship is
σij=k0+k1z+k2z2+k3z3+k4z4+k5z5+k6z6
Wherein, σijFor the components of stress, z is the distance of the finished surface apart from described plate part, k0-k6For undetermined constant.
According to one embodiment of present invention, described step E) comprise the following steps: in finite element analysis software, E-1) set up the new model of described plate part, described model be shaped as regular hexahedron, the length and width of described model is respectively equal to the length and width of described plate part, and the thickness of described model deducts described name milling depth equal to the original depth of described plate part；E-2) to described model partition grid to obtain multiple hexahedral element, it is preferable that the height of the grid of the upper surface of neighbouring described model is less than the height of the grid of the upper surface away from described model；E-3) from finite element analysis software, described hexahedral element and nodal information are extracted；E-4) coordinate of each node is extracted, according to described node coordinate on the length direction of described plate part and described actual milling depth, find immediate described cutting parameter with it from organize described cutting parameter more, and then determine corresponding described functional relationship, the numerical value of the corresponding described components of stress is then determined according to described hexahedral element coordinate on the thickness direction of described plate part；And E-5) each components of stress of whole nodes of each described hexahedral element are averaged as the residual stress state of described hexahedral element.
According to one embodiment of present invention, described step F) comprise the following steps: F-1) load the stress state of described hexahedral element；F-2) lower surface at described model applies the distributed force load equal, in opposite direction with the surface density numerical value of described equivalent action power；And F-3) submit finite element analysis to, and obtain deformation result.
Accompanying drawing explanation
Fig. 1 is the flow chart of plate part machining deformation simulated prediction method according to embodiments of the present invention；
Fig. 2 is the operating process of plate part machining deformation simulated prediction method according to embodiments of the present invention；
Fig. 3 is the schematic diagram before and after clamping of the plate part according to embodiments of the present invention；
Fig. 4 is the schematic diagram of the actual milling depth of plate part according to embodiments of the present invention；
Fig. 5 is cutting residual stress field and the loading schematic diagram of distributed force load of plate part machining deformation simulated prediction method according to embodiments of the present invention；
Fig. 6 is the cutting simulation schematic diagram of plate part machining deformation simulated prediction method according to embodiments of the present invention；
Fig. 7 is the machining deformation simulation result schematic diagram of plate part machining deformation simulated prediction method according to embodiments of the present invention.
Detailed description of the invention
Embodiments of the invention are described below in detail, and the example of described embodiment is shown in the drawings.The embodiment described below with reference to accompanying drawing is exemplary, it is intended to is used for explaining the present invention, and is not considered as limiting the invention.
Below with reference to the accompanying drawings plate part machining deformation simulated prediction method according to embodiments of the present invention is described.As shown in Fig. 1-Fig. 7, plate part machining deformation simulated prediction method according to embodiments of the present invention comprises the following steps:
A) the upper surface profile of table record plate part 100 is beaten；
B) clamping plate part 100, then beats the upper surface profile of table record plate part 100, calculates the deformation caused because of this clamping, and then calculates the surface density causing the equivalent action power needed for this deformation；
C) the actual milling depth of each position of plate part 100 is determined according to the upper surface profile of the plate part 100 after clamping and name milling depth, then according to the span of this actual milling depth, determining multiple real cutting depth numerical value, these real cutting depth numerical value multiple are constituted with the speed of mainshaft, feed engagement and axial milling depth organizes cutting parameter more；
D) finite element analysis software is utilized to carry out cutting the simulation analysis of residual stress according to this many groups cutting parameter, simulation analysis extracts the residual stress numerical value of multiple point after completing, be fitted the residual stress numerical value of the plurality of point to obtain the functional relationship of multiple components of stress and the distance of the finished surface of distance plate part 100；
E) in finite element analysis software, set up the model of plate part 100, to this model partition grid to obtain multiple unit, the actual cut degree of depth of the different length position according to this model, select this corresponding functional relationship, determine its corresponding residual stress state according to the distance of the upper surface of this this model of cell distance each；
F) load cutting residual stress field, and apply the distributed force load equal, in opposite direction with the surface density numerical value of this equivalence active force, submission finite element analysis at the lower surface of this model, and obtain deformation result；With
G) if this deformation result meets requirement, then select this step C) in cutting parameter for processing cutting parameter, if this deformation result is unsatisfactory for requirement, then change this step C) in cutting parameter and repeat this D) to this step F), until this deformation result meets requirement.
Plate part machining deformation simulated prediction method according to embodiments of the present invention is by implementing based on working angles finite element simulation and cutting residual stress field and the loading of distributed force load, such that it is able to effectively, accurately to plate part 100 carry out Deformation Prediction, thus can realize the active to plate part 100 Milling Process rear surface profile to control, it is also possible to develop and be applied to prediction and the control of complex thin-wall component machining deformation further.
That is, plate part machining deformation simulated prediction method according to embodiments of the present invention is by being predicted cutting residual stress and machining deformation, and control cutting residual stress by regulation cutting parameter, to cut residual stress payment part initial deformation and the impact of clamping effect.
As shown in Fig. 1-Fig. 7, comprise the following steps according to the plate part machining deformation simulated prediction method of some embodiments of the present invention:
Plate part 100 is placed on magnetic suction disc, beats the upper surface profile of table record plate part 100.Wherein, plate part 100 is ferromagnetic part, and the length of plate part 100 is more than or equal to 5:1 more than or equal to 3:1, the width of plate part 100 with the ratio of the thickness of plate part 100 with the ratio of the width of plate part 100.It is to say, the length of plate part 100 is at least three times of the width of plate part 100, the width of plate part 100 is at least five times of the thickness of plate part 100.
Open this magnetic suction disc, in order to clamping plate part 100, after clamping is stable, beat the upper surface profile of table record plate part 100.Upper surface profile according to plate part before and after clamping 100, calculates the deformation caused because of this clamping.Wherein, the deformation caused because of this clamping can be amount of deflection δ in Fig. 3.
Plate part 100 is reduced to simply supported beam model, and a length of l of plate part 100, width are b, thickness is h, elasticity modulus of materials is E, then the line density of equivalent action powerThe surface density of equivalent action power
Then, it is determined that the speed of mainshaft, feed engagement and axial milling depth, and determine the actual milling depth of plate part 100 each position along its length according to the upper surface profile of the plate part 100 after clamping and name milling depth.Wherein, this actual milling depth is as shown in the dash area in Fig. 4.
According to the span of this actual milling depth, determine multiple real cutting depth numerical value.For example, it is assumed that the span of this actual milling depth is 0.9 millimeter-1.1 millimeters, then can with 0.02 millimeter as step-length, determine 10 real cutting depth numerical value (0.91,0.93,0.95 ..., 1.07,1.09).
It is to say, the actual milling depth of 0.9 millimeter-1.1 millimeters is reduced to 10 real cutting depth numerical value.Such as, when actual milling depth is 0.90 millimeter-0.92 millimeter, actual milling depth can be reduced to 0.91 millimeter, when actual milling depth is 0.92 millimeter-0.94 millimeter, actual milling depth can be reduced to 0.93 millimeter, when actual milling depth is 0.94 millimeter-0.96 millimeter, actual milling depth can be reduced to 0.95 millimeter, the rest may be inferred, finally when actual milling depth is 1.08 millimeters-1.10 millimeters, actual milling depth can be reduced to 1.09 millimeters, to improve the speed of simulation calculation.
These real cutting depth numerical value multiple are constituted with the fixing speed of mainshaft, feed engagement and axial milling depth organizes cutting parameter more.Preferably, determine 5-10 real cutting depth numerical value, and then constitute 5-10 group cutting parameter.In other words, the speed of mainshaft of different groups, feed engagement and axial milling depth are to maintain constant.
Utilize threedimensional model and the threedimensional model of plate part 100 of the blade of CAD software drafting Integral milling cutter or indexable milling cutter, intercept part (part i.e. contacted) and the actual part (i.e. the part being milled of plate part 100) used of processing of plate part 100 that the actual processing of this blade is used with plate part 100, in importing finite element analysis software.
Arrange cutting parameter, cutter and plate part 100 contacts attribute, the material properties of cutter and geometric parameter and the material properties of plate part 100 and constitutive equation, utilizes this finite element analysis software to carry out three-dimensional RAP material.Wherein, this finite element analysis software can be the general finite element analysis softwares such as ANSYS, Abaqus, it is also possible to is the specialty finite element analysis softwares such as AdvantEdge, Deform.Such as, simulation example based on AdvantEdge is as shown in Figure 6.
Advantageously, the processing mode of plate part 100 is face milling, and this face milling is the milling cutter length direction feeding along plate part 100, uses the mode of symmetrical milling.Wherein, the diameter of milling cutter is more than the width of plate part 100.Preferably, the diameter of this milling cutter and the ratio of the width of plate part 100 are more than or equal to 3:1.In other words, the diameter of this milling cutter is at least three times of width of plate part 100.
After this three-dimensional RAP material terminates, load removal, the finished surface (i.e. the upper surface of plate part 100) of extraction plate part 100 and the residual stress numerical value of following multiple points are to be analyzed, and the spacing of adjacent two these points is 15 microns-25 microns.Advantageously, the spacing of adjacent two these points is 20 microns.
The plurality of point is fitted (being i.e. fitted multiple data) to obtain the functional relationship of multiple components of stress and the distance of the finished surface of distance plate part 100.In other words, this functional relationship is the functional relationship of multiple components of stress and the point on the thickness direction of plate part 100.
In one embodiment of the invention, utilizing 6 rank multinomials to be fitted the plurality of point, obtain the functional relationship of multiple components of stress and the distance of the finished surface of distance plate part 100, this functional relationship is
σij=k0+k1z+k2z2+k3z3+k4z4+k5z5+k6z6
Wherein, σijFor the components of stress, z is the distance of the finished surface of distance plate part 100, k0-k6For undetermined constant.
Change the cutting parameter in phantom and repeat this three-dimensional RAP material, in order to completing the simulation analysis of this many groups cutting parameter.It is to say, change this real cutting depth numerical value, repeat this three-dimensional RAP material to complete the simulation analysis of this many groups cutting parameter.
Then, the new model of plate part 100 is set up in finite element analysis software, this model be shaped as regular hexahedron, the length of this model is equal to the length of plate part 100, the width of this model deducts this name milling depth equal to the width of plate part 100, the thickness of this model equal to the original depth of plate part 100
To this model partition grid to obtain multiple hexahedral element (such as C3D8T).In an example of the present invention, the height of the grid of the upper surface of this model neighbouring is less than the height of the grid of the upper surface away from this model.Each grid is exactly a hexahedral element.
It is to say, the grid near the upper surface of this model needs to divide more delicate, in order to the loading of cutting residual stress.Such as, the height of grid is 10 microns-20 microns.Grid away from the upper surface of this model can divide more sparse, to improve the speed of simulation calculation.
This hexahedral element and nodal information is extracted from finite element analysis software.Extract the coordinate of each node, according to this node coordinate on the length direction of plate part 100 and this actual milling depth, from organize, this cutting parameter finds this cutting parameter immediate with it more, and then determine this functional relationship corresponding.
Specifically, this node coordinate on the length direction of this model is exactly this node coordinate on the length direction of plate part 100.By this node coordinate on the length direction of plate part 100, this corresponding actual milling depth can be obtained.That is, it is determined that a point on the length direction of plate part 100, it is possible to determine this actual milling depth at this point.This real cutting depth numerical value (as mentioned above) can be obtained by this actual milling depth, can from organize this cutting parameter find with it immediate (corresponding) this cutting parameter by this real cutting depth numerical value more, and then determine this functional relationship corresponding.
Then the numerical value of these components of stress corresponding is determined according to this hexahedral element coordinate on the thickness direction of plate part 100.Each components of stress of whole nodes of this hexahedral element each are averaged as the residual stress state of this hexahedral element.Specifically, the residual stress state of this hexahedral element is write the file needed for the primary stress input of corresponding finite element analysis software.
Loading the stress state of this hexahedral element, the stress state of these hexahedral elements multiple constitutes this cutting residual stress field (as shown in the dash area in Fig. 5).Specifically, the file of the residual stress state comprising this hexahedral element exported by previous step is loaded in this model (the new model of plate part 100 the most described above), as Abaqus software can define the initial stress state of this model by editor's keyword (keywords), the hexahedral element containing this initial stress state is distributed mainly on the upper surface layer of this model.
Lower surface at this model applies the distributed force load equal, in opposite direction with the surface density numerical value of this equivalence active force, then submits finite element analysis to, and obtains deformation result.
If this deformation result meets requirement, then this step C) in cutting parameter for processing cutting parameter, if this deformation result is unsatisfactory for requirement, then change this step C) in cutting parameter and repeat this D) to this step F), until this deformation result meets requirement.Wherein, processing cutting parameter is the speed of mainshaft, feed engagement and axial milling depth.
It is an object of the invention to problem that is relatively big for elongated sheet part milling deformation and that be difficult to predict and control, the present invention solves elongated sheet part Deformation in Milling Process forecasting problem, and can realize the control of the active to elongated sheet part Milling Process rear surface profile further in smaller range.Analyze method accordingly, it is also possible to develop and be applied to prediction and the control of complex thin-wall component machining deformation further.
In describing the invention, it will be appreciated that, term " " center ", " longitudinally ", " laterally ", " length ", " width ", " thickness ", on " ", D score, " front ", " afterwards ", " left ", " right ", " vertically ", " level ", " push up ", " end " " interior ", " outward ", " clockwise ", " counterclockwise ", " axially ", " radially ", orientation or the position relationship of the instruction such as " circumferential " are based on orientation shown in the drawings or position relationship, it is for only for ease of the description present invention and simplifies description, rather than indicate or imply that the device of indication or element must have specific orientation, with specific azimuth configuration and operation, therefore it is not considered as limiting the invention.
Additionally, term " first ", " second " are only used for describing purpose, and it is not intended that instruction or hint relative importance or the implicit quantity indicating indicated technical characteristic.Thus, define " first ", the feature of " second " can express or implicitly include at least one this feature.In describing the invention, " multiple " are meant that at least two, such as two, three etc., unless otherwise expressly limited specifically.
In the present invention, unless otherwise clearly defined and limited, term " is installed ", " being connected ", " connection ", the term such as " fixing " should be interpreted broadly, and connects for example, it may be fixing, it is also possible to be to removably connect, or integral；Can be mechanically connected, it is also possible to be electrical connection or each other can communication；Can be to be joined directly together, it is also possible to be indirectly connected to by intermediary, can be connection or the interaction relationship of two elements of two element internals, unless otherwise clear and definite restriction.For the ordinary skill in the art, above-mentioned term concrete meaning in the present invention can be understood as the case may be.
In the present invention, unless otherwise clearly defined and limited, fisrt feature second feature " on " or D score can be that the first and second features directly contact, or the first and second features are by intermediary mediate contact.And, fisrt feature second feature " on ", " top " and " above " but fisrt feature directly over second feature or oblique upper, or be merely representative of fisrt feature level height higher than second feature.Fisrt feature second feature " under ", " lower section " and " below " can be fisrt feature immediately below second feature or obliquely downward, or be merely representative of fisrt feature level height less than second feature.
In the description of this specification, the description of reference term " embodiment ", " some embodiments ", " example ", " concrete example " or " some examples " etc. means that the specific features, structure, material or the feature that combine this embodiment or example description are contained at least one embodiment or the example of the present invention.In this manual, the schematic representation of above-mentioned term is necessarily directed to identical embodiment or example.And, the specific features of description, structure, material or feature can be to combine in one or more embodiments in office or example in an appropriate manner.Additionally, in the case of the most conflicting, the feature of the different embodiments described in this specification or example and different embodiment or example can be combined and combine by those skilled in the art.
Although above it has been shown and described that embodiments of the invention, it is understandable that, above-described embodiment is exemplary, it is impossible to be interpreted as limitation of the present invention, and above-described embodiment can be changed, revises, replace and modification by those of ordinary skill in the art within the scope of the invention.
Claims (9)
1. a plate part machining deformation simulated prediction method, it is characterised in that comprise the following steps:
A) the upper surface profile of plate part described in table record is beaten；
B) plate part described in clamping, then beats the upper surface profile of plate part described in table record, calculates the deformation caused because of described clamping, and then calculates the surface density causing the equivalent action power needed for described deformation；
C) the actual milling depth of described each position of plate part is determined according to the upper surface profile of the described plate part after clamping and name milling depth, then according to the span of described actual milling depth, determining multiple real cutting depth numerical value, multiple described real cutting depth numerical value are constituted with the speed of mainshaft, feed engagement and axial milling depth organizes cutting parameter more；
D) simulation analysis utilizing finite element analysis software to carry out cutting residual stress according to described many group cutting parameters, simulation analysis extracts the residual stress numerical value of multiple point after completing, be fitted the residual stress numerical value of the plurality of point to obtain the functional relationship of multiple components of stress and the distance of the finished surface apart from described plate part；
E) in finite element analysis software, set up the model of described plate part, to described model partition grid to obtain multiple unit, the actual cut degree of depth of the different length position according to described model, select corresponding described functional relationship, determine its corresponding residual stress state according to the distance of the upper surface of model described in each described cell distance；
F) load cutting residual stress field, and apply the distributed force load equal, in opposite direction with the surface density numerical value of described equivalent action power, submission finite element analysis at the lower surface of described model, and obtain deformation result；With
G) if described deformation result meets requirement, then select described step C) in cutting parameter for processing cutting parameter, if described deformation result is unsatisfactory for requirement, then change described step C) in cutting parameter and repeat described D) to described step F), until described deformation result meets requirement.
Plate part machining deformation simulated prediction method the most according to claim 1, it is characterized in that, the length of described plate part is more than or equal to 5:1 more than or equal to 3:1, the width of described plate part with the ratio of the thickness of described plate part with the ratio of the width of described plate part.
Plate part machining deformation simulated prediction method the most according to claim 1, it is characterized in that, described plate part is ferromagnetic part, described plate part magnetic suction disc clamping, the processing mode of described plate part is face milling, described face milling is the milling cutter length direction feeding along described plate part, use the mode of symmetrical milling, wherein the diameter of milling cutter is more than the width of described plate part, advantageously, the diameter of described milling cutter and the ratio of the width of described plate part are more than or equal to 3:1.
Plate part machining deformation simulated prediction method the most according to claim 1, it is characterized in that, in described step A) in, described plate part is placed on magnetic suction disc, beat the upper surface profile of plate part, described step B described in table record) comprise the following steps:
B-1) open described magnetic suction disc, in order to plate part described in clamping, after clamping is stable, beat the upper surface profile of plate part described in table record；
B-2) according to the upper surface profile of plate part described before and after clamping, the deformation caused because of described clamping is calculated；With
B-3) described plate part being reduced to simply supported beam model, a length of l of described plate part, width are b, thickness is h, elasticity modulus of materials is E, then the line density of equivalent action powerThe surface density of equivalent action power
Plate part machining deformation simulated prediction method the most according to claim 4, it is characterised in that described step C) comprise the following steps:
C-1) speed of mainshaft, feed engagement and axial milling depth are determined；
C-2) the actual milling depth of described plate part each position along its length is determined according to the upper surface profile of the described plate part after clamping and name milling depth；
C-3) according to the span of described actual milling depth, determine multiple real cutting depth numerical value, multiple described real cutting depth numerical value are constituted with the fixing speed of mainshaft, feed engagement and axial milling depth organizes cutting parameter more, preferably, determine 5-10 real cutting depth numerical value, and then constitute 5-10 group cutting parameter.
Plate part machining deformation simulated prediction method the most according to claim 5, it is characterised in that described step D) comprise the following steps:
D-1) threedimensional model and the threedimensional model of described plate part of the blade of CAD software drafting Integral milling cutter or indexable milling cutter are utilized, intercept part that the actual processing of described blade uses and the part that the actual processing of described plate part is used, import in finite element analysis software；
D-2) material properties of cutter, contact attribute and cutting parameter and the material properties of described plate part are set, carry out three-dimensional RAP material；
D-3) after described three-dimensional RAP material terminates, load removal, extract finished surface and the residual stress numerical value of following multiple points of described plate part, the spacing of adjacent two described points is 15 microns-25 microns, preferably, the spacing of adjacent two described points is 20 microns, is fitted the residual stress numerical value of the plurality of point to obtain the functional relationship of multiple components of stress and the distance of the finished surface apart from described plate part；With
D-4) change the cutting parameter in phantom and repeat described three-dimensional RAP material, in order to completing the described simulation analysis organizing cutting parameter more.
Plate part machining deformation simulated prediction method the most according to claim 6, it is characterized in that, utilizing 6 rank multinomials to be fitted the plurality of point, obtain the functional relationship of multiple components of stress and the distance of the finished surface apart from described plate part, described functional relationship is
σij=k0+k1z+k2z2+k3z3+k4z4+k5z5+k6z6
Wherein, σijFor the components of stress, z is the distance of the finished surface apart from described plate part, k0-k6For undetermined constant.
Plate part machining deformation simulated prediction method the most according to claim 6, it is characterised in that described step E) comprise the following steps:
E-1) in finite element analysis software, set up the new model of described plate part, described model be shaped as regular hexahedron, the length and width of described model is respectively equal to the length and width of described plate part, and the thickness of described model deducts described name milling depth equal to the original depth of described plate part；
E-2) to described model partition grid to obtain multiple hexahedral element, it is preferable that the height of the grid of the upper surface of neighbouring described model is less than the height of the grid of the upper surface away from described model；
E-3) from finite element analysis software, described hexahedral element and nodal information are extracted；
E-4) coordinate of each node is extracted, according to described node coordinate on the length direction of described plate part and described actual milling depth, find immediate described cutting parameter with it from organize described cutting parameter more, and then determine corresponding described functional relationship, the numerical value of the corresponding described components of stress is then determined according to described hexahedral element coordinate on the thickness direction of described plate part；With
E-5) each components of stress of whole nodes of each described hexahedral element are averaged as the residual stress state of described hexahedral element.
Plate part machining deformation simulated prediction method the most according to claim 8, it is characterised in that described step F) comprise the following steps:
F-1) stress state of described hexahedral element is loaded；
F-2) lower surface at described model applies the distributed force load equal, in opposite direction with the surface density numerical value of described equivalent action power；With
F-3) submit finite element analysis to, and obtain deformation result.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610191266.8A CN105825025B (en) | 2016-03-30 | 2016-03-30 | Plate part machining deformation simulated prediction method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610191266.8A CN105825025B (en) | 2016-03-30 | 2016-03-30 | Plate part machining deformation simulated prediction method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105825025A true CN105825025A (en) | 2016-08-03 |
| CN105825025B CN105825025B (en) | 2019-01-15 |
Family
ID=56525336
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610191266.8A Active CN105825025B (en) | 2016-03-30 | 2016-03-30 | Plate part machining deformation simulated prediction method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105825025B (en) |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106529086A (en) * | 2016-12-07 | 2017-03-22 | 同济大学 | Static load intensity finite element simulation precision determination method used for mechanical structural part |
| CN106599496A (en) * | 2016-12-20 | 2017-04-26 | 上海电气电站设备有限公司 | Method for measuring residual stress based on ring core method of numerical calculation |
| CN106709148A (en) * | 2016-11-25 | 2017-05-24 | 中国商用飞机有限责任公司 | Finite element simulation method for roll bending-milling process of large-scale stepped thin-walled workpiece |
| CN106777630A (en) * | 2016-12-07 | 2017-05-31 | 同济大学 | For the determination method of the mode finite element simulation precision of machine components |
| CN107145662A (en) * | 2017-05-04 | 2017-09-08 | 北京理工大学 | A kind of meso-scale turning Deformation Prediction method |
| CN109213082A (en) * | 2017-11-24 | 2019-01-15 | 中国航空制造技术研究院 | Complex curved surface parts machining Deformation Prediction method |
| CN109581971A (en) * | 2017-09-28 | 2019-04-05 | 中国商用飞机有限责任公司 | The deformation control method of aircraft entirety escape door numerical control processing |
| CN110532708A (en) * | 2019-09-02 | 2019-12-03 | 哈尔滨工业大学（深圳） | A kind of residual stress prediction technique, device, server and storage medium |
| CN112036057A (en) * | 2020-07-28 | 2020-12-04 | 成都飞机工业（集团）有限责任公司 | Hole position control method for hole making of aircraft structural part |
| CN112307575A (en) * | 2020-09-24 | 2021-02-02 | 成都飞机工业（集团）有限责任公司 | Method for checking qualification of composite material wallboard part for judging assembly result |
| CN113084233A (en) * | 2021-04-06 | 2021-07-09 | 湖南大学 | High-precision cutting machining method for section bar |
| CN116992778A (en) * | 2023-08-30 | 2023-11-03 | 西北工业大学 | Milling residual stress deformation prediction method, storage medium and electronic equipment |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120051853A1 (en) * | 2009-04-30 | 2012-03-01 | Seco Tools Ab | cutting tool and an insert holder for a cutting tool |
| CN103500252A (en) * | 2013-09-26 | 2014-01-08 | 西安工业大学 | Thin-walled workpiece cutting condition optimization machining method based on finite element analysis |
| CN104077442A (en) * | 2014-06-24 | 2014-10-01 | 南京晨光集团有限责任公司 | Method for controlling machining accuracy of large integrated thin-walled parts based on finite element analysis |
| CN104899355A (en) * | 2015-05-11 | 2015-09-09 | 沈阳航空航天大学 | Finite-element prediction method for machining distortion of large complex structural member |
-
2016
- 2016-03-30 CN CN201610191266.8A patent/CN105825025B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120051853A1 (en) * | 2009-04-30 | 2012-03-01 | Seco Tools Ab | cutting tool and an insert holder for a cutting tool |
| CN103500252A (en) * | 2013-09-26 | 2014-01-08 | 西安工业大学 | Thin-walled workpiece cutting condition optimization machining method based on finite element analysis |
| CN104077442A (en) * | 2014-06-24 | 2014-10-01 | 南京晨光集团有限责任公司 | Method for controlling machining accuracy of large integrated thin-walled parts based on finite element analysis |
| CN104899355A (en) * | 2015-05-11 | 2015-09-09 | 沈阳航空航天大学 | Finite-element prediction method for machining distortion of large complex structural member |
Cited By (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106709148B (en) * | 2016-11-25 | 2020-02-14 | 中国商用飞机有限责任公司 | Finite element simulation method for roll bending-milling process of large-size thin-walled part with step |
| CN106709148A (en) * | 2016-11-25 | 2017-05-24 | 中国商用飞机有限责任公司 | Finite element simulation method for roll bending-milling process of large-scale stepped thin-walled workpiece |
| CN106777630A (en) * | 2016-12-07 | 2017-05-31 | 同济大学 | For the determination method of the mode finite element simulation precision of machine components |
| CN106529086A (en) * | 2016-12-07 | 2017-03-22 | 同济大学 | Static load intensity finite element simulation precision determination method used for mechanical structural part |
| CN106599496A (en) * | 2016-12-20 | 2017-04-26 | 上海电气电站设备有限公司 | Method for measuring residual stress based on ring core method of numerical calculation |
| CN106599496B (en) * | 2016-12-20 | 2019-12-17 | 上海电气电站设备有限公司 | Method for measuring residual stress by ring core method based on numerical calculation |
| CN107145662A (en) * | 2017-05-04 | 2017-09-08 | 北京理工大学 | A kind of meso-scale turning Deformation Prediction method |
| CN107145662B (en) * | 2017-05-04 | 2019-09-17 | 北京理工大学 | A kind of meso-scale turning Deformation Prediction method |
| CN109581971A (en) * | 2017-09-28 | 2019-04-05 | 中国商用飞机有限责任公司 | The deformation control method of aircraft entirety escape door numerical control processing |
| CN109213082A (en) * | 2017-11-24 | 2019-01-15 | 中国航空制造技术研究院 | Complex curved surface parts machining Deformation Prediction method |
| CN110532708A (en) * | 2019-09-02 | 2019-12-03 | 哈尔滨工业大学（深圳） | A kind of residual stress prediction technique, device, server and storage medium |
| CN110532708B (en) * | 2019-09-02 | 2023-01-06 | 哈尔滨工业大学（深圳） | Residual stress prediction method, residual stress prediction device, server and storage medium |
| CN112036057A (en) * | 2020-07-28 | 2020-12-04 | 成都飞机工业（集团）有限责任公司 | Hole position control method for hole making of aircraft structural part |
| CN112307575A (en) * | 2020-09-24 | 2021-02-02 | 成都飞机工业（集团）有限责任公司 | Method for checking qualification of composite material wallboard part for judging assembly result |
| CN112307575B (en) * | 2020-09-24 | 2022-05-06 | 成都飞机工业（集团）有限责任公司 | Method for checking qualification of composite material wallboard part for judging assembly result |
| CN113084233A (en) * | 2021-04-06 | 2021-07-09 | 湖南大学 | High-precision cutting machining method for section bar |
| CN113084233B (en) * | 2021-04-06 | 2022-04-29 | 湖南大学 | High-precision cutting machining method for section bar |
| CN116992778A (en) * | 2023-08-30 | 2023-11-03 | 西北工业大学 | Milling residual stress deformation prediction method, storage medium and electronic equipment |
| CN116992778B (en) * | 2023-08-30 | 2024-04-26 | 西北工业大学 | Milling residual stress deformation prediction method, storage medium and electronic equipment |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105825025B (en) | 2019-01-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105825025A (en) | Plate-shaped part processing deformation simulation prediction method | |
| CN108182325B (en) | Prediction analysis method for machining deformation of thin-wall structural part | |
| CN105760624B (en) | It is a kind of support extensive three dimensional integrated circuits heat emulation and Thermal design | |
| Villalpando et al. | An optimization approach for components built by fused deposition modeling with parametric internal structures | |
| CN108062432A (en) | A kind of method for numerical simulation of selective laser fusion process | |
| Górski et al. | Computation of mechanical properties of parts manufactured by fused deposition modeling using finite element method | |
| Yang et al. | Towards a numerical approach of finding candidates for additive manufacturing-enabled part consolidation | |
| CN101641699A (en) | Analytical equipment, analytical approach and routine analyzer | |
| Shimada | Current issues and trends in meshing and geometric processing for computational engineering analyses | |
| Matras et al. | Analysis of machining accuracy during free form surface milling simulation for different milling strategies | |
| Kahhal et al. | Multi-objective optimization of sheet metal forming die using FEA coupled with RSM | |
| Damanpack et al. | Boundary Element Method applied to the bending analysis of thin functionally graded plates | |
| Xi et al. | A prediction model of the cutting force–induced deformation while considering the removed material impact | |
| KR20050047237A (en) | Method and apparatus for generating 3-dimensional finite element mesh | |
| Scurtu et al. | Design optimization method used in mechanical engineering | |
| Pierce et al. | A method for integrating form errors into geometric tolerance analysis | |
| Kamali Nejad et al. | 3D simulation of manufacturing defects for tolerance analysis | |
| Cao et al. | The strategy of datum reference frame selection based on statistical learning | |
| CN113094960B (en) | Method for quickly generating electromagnetic model of multi-core aviation connector based on moment method | |
| CN102819237A (en) | Method for generating simulation target blank in solid milling simulation process | |
| Gallage et al. | Formation of the three-dimensional geometry of the red blood cell membrane | |
| CN114222659A (en) | Method of generating a structural grid, method of using a structural grid, computer program and computer readable medium | |
| CN107273558A (en) | The method and apparatus that mechanical performance for threedimensional model is verified | |
| Haas | Prediction of Structural Reliability Through an Alternative Variability-Based Methodology | |
| JP6969757B2 (en) | Analysis result data reduction device, analysis result data reduction method and analysis result data reduction program |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |