CN109635348A - A kind of spatial complex curved surface grinding temperature Numerical Simulation Analysis method - Google Patents

A kind of spatial complex curved surface grinding temperature Numerical Simulation Analysis method Download PDF

Info

Publication number
CN109635348A
CN109635348A CN201811362258.0A CN201811362258A CN109635348A CN 109635348 A CN109635348 A CN 109635348A CN 201811362258 A CN201811362258 A CN 201811362258A CN 109635348 A CN109635348 A CN 109635348A
Authority
CN
China
Prior art keywords
grinding
curved surface
numerical simulation
simulation analysis
temperature
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
Application number
CN201811362258.0A
Other languages
Chinese (zh)
Other versions
CN109635348B (en
Inventor
杨钱
汪东政
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hubei Sanjiang Aerospace Hongfeng Control Co Ltd
Original Assignee
Hubei Sanjiang Aerospace Hongfeng Control Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hubei Sanjiang Aerospace Hongfeng Control Co Ltd filed Critical Hubei Sanjiang Aerospace Hongfeng Control Co Ltd
Priority to CN201811362258.0A priority Critical patent/CN109635348B/en
Publication of CN109635348A publication Critical patent/CN109635348A/en
Application granted granted Critical
Publication of CN109635348B publication Critical patent/CN109635348B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Grinding Of Cylindrical And Plane Surfaces (AREA)
  • Constituent Portions Of Griding Lathes, Driving, Sensing And Control (AREA)

Abstract

The invention discloses a kind of spatial complex curved surface grinding temperature Numerical Simulation Analysis methods, comprising the following steps: S1: calculating the curved surface according to the per unit system of the curved surface of workpiece to be processed arrow and normal curvature and takes up an official post two principal curvatures and its corresponding two principal direction of a fixed point；S2: long axis, short axle and the area of the instantaneous contact ellipse of workpiece to be processed and grinding wheel are calculated according to principal curvatures and its corresponding principal direction, the first grinding process parameters；S3: determining grinding method, calculates grinding heat flow density according to the characterization coefficient of the grinding method, workpiece material inherent characteristic parameter, the area of instantaneous contact ellipse and the second grinding parameter；S4: grinding temperature is calculated according to grinding heat flow density；The present invention establishes the Three-dimensional simulation analysis model of Grinding Process of Curved Surface parameter and temperature based on model of ellipse, realizes the Three-dimensional simulation analysis to grinding temperature in the components grinding with spatial complex curved surface feature.

Description

A kind of spatial complex curved surface grinding temperature Numerical Simulation Analysis method
Technical field
The invention belongs to field of engineering technology, more particularly, to a kind of spatial complex curved surface grinding temperature numerical simulation Analysis method.
Background technique
Grinding has the characteristics that high specific cutting force, fast grinding linear velocity, the energy in processing almost turn Turn to thermal energy.And the distribution of grinding temperature has highly significant shadow to the thermal damage of workpiece and surface integrity in grinding It rings.Therefore, particularly important in grinding by being distributed in for Numerical Simulation Analysis grinding temperature.
Grinding temperature Numerical Simulation Analysis model achieves initial success earlier than the 1950s, classical grinding temperature Numerical Simulation Analysis is to establish based on rectangle heat source is theoretical or triangle heat source is theoretical in conjunction with finite element method Grinding temperature Numerical Simulation Analysis model, can accurately analytical calculation is one-dimensional, in two-dimentional grinding workpiece mill Cut temperature.
With the rapid development of major machinery industry such as aerospace, the service performance of product is stepped up, spatial complex The product of curved surface is more and more, and the requirement to its quality and service performance is higher and higher.But the grinding temperature in grinding Numerical Simulation Analysis method fail to break through two-dimension analysis always and calculate, the Numerical Simulation Analysis research of grinding temperature stagnation is not Before, often requirement is not achieved because of thermal burn in the product grinding with spatial complex curved surface feature.
Summary of the invention
For at least one defect or Improvement requirement of the prior art, the present invention provides a kind of grindings of spatial complex curved surface Temperature value simulating analysis establishes the spatial complex curved surface grinding work based on model of ellipse in conjunction with space contact theory The three-dimensional numerical simulation analysis model of skill parameter and temperature misses to solve one-dimensional, two-dimentional grinding temperature Numerical Simulation Analysis The big problem of difference realizes imitative to the three-dimensional numerical value of grinding temperature in the components grinding with spatial complex curved surface feature True analysis, and introduce grinding temperature numerical simulation and calculate correction factor, improve grinding temperature Numerical Simulation Analysis precision.
To achieve the above object, according to one aspect of the present invention, a kind of spatial complex curved surface grinding temperature number is provided It is worth simulating analysis, comprising the following steps:
S1: the curved surface is calculated according to the per unit system of the curved surface of workpiece to be processed arrow and normal curvature and is taken up an official post two masters of a fixed point Curvature and its corresponding two principal direction；
S2: ellipsometric analysis mould is established according to two principal curvatures and its corresponding principal direction and the first grinding process parameters Type calculates long axis, short axle and the area of the instantaneous contact ellipse of workpiece to be processed and grinding wheel；The first grinding process parameters packet Include the amount of collapsing and grinding depth；
S3: determining grinding method, according to the characterization coefficient of the grinding method, workpiece material inherent characteristic parameter, instantaneous The area of Contact Ellipse and the second grinding parameter calculate grinding heat flow density；Second grinding process parameters include that grinding connects Touching effective radius, transforms into speed and tangential grinding force at grinding wheel speed；
S4: grinding temperature is calculated according to the grinding heat flow density.
Preferably, above-mentioned spatial complex curved surface grinding temperature Numerical Simulation Analysis method, instantaneous contact ellipse are grinding The Contact Ellipse of process medium plain emery wheel and workpiece to be processed projects to processing curve again after plane projection and is formed by ellipse Curved surface load scope.
Preferably, above-mentioned spatial complex curved surface grinding temperature Numerical Simulation Analysis method, further include after step S4 with Lower step:
S5: calculating the error coefficient of model of ellipse according to cutting angle between processing curve and projection plane, with the mistake Poor coefficient carries out numerical simulation amendment to the grinding temperature as temperature correction coefficient, obtains modified grinding temperature.
Preferably, above-mentioned spatial complex curved surface grinding temperature Numerical Simulation Analysis method includes following son in step S3 Step:
S31: determining grinding method, according to the characterization coefficient of the grinding method and workpiece material inherent characteristic parameter, Grinding Contact effective radius and grinding wheel speed calculate the heat distribution coefficient in grinding process；
S32: according to the heat distribution coefficient, grinding wheel speed, speed, tangential grinding force and instantaneous contact ellipse are transformed into Areal calculation is ground heat flow density.
Preferably, above-mentioned spatial complex curved surface grinding temperature Numerical Simulation Analysis method, in step S1, workpiece to be processed Curved surface is taken up an official post two principal curvatures k of a fixed pointn1, kn2And its corresponding two principal direction λ1, λ2It calculates as follows:
(E·M-F·L)·λ2+ (EN-GL) λ+(FN-GM)=0
Preferably, above-mentioned spatial complex curved surface grinding temperature Numerical Simulation Analysis method, in step S1, workpiece to be processed The per unit system of curved surface r=r (μ, υ) swears n and normal curvature knIt calculates as follows:
N=rμ×rυ/|rμ×rυ|
Wherein, rμIt is surface equation r to the first-order partial derivative of independent variable μ；rυIt is inclined for single order of the surface equation r to independent variable υ Derivative；
The normal curvature k of any fixed point on curved surface r=r (μ, υ) along direction λ=d μ/d υnAre as follows:
Wherein, E, F, G are the first fundamental quantity of curved surface r；L, M, N are the second fundamental quantity of curved surface r.
Preferably, above-mentioned spatial complex curved surface grinding temperature Numerical Simulation Analysis method, in step S2, the instant contact Elliptical long axis a and short axle b calculates as follows:
Wherein, δ is the amount of collapsing；apFor grinding depth；
kn1、kn2Respectively two principal curvatures of space curved surface and grinding wheel contact point, αsFor two principal direction λ1、λ2Between folder Angle.
Preferably, above-mentioned spatial complex curved surface grinding temperature Numerical Simulation Analysis method, in step S3, the heat point Distribution coefficient RwIt calculates as follows:
Wherein, ξ is the characterization coefficient of grinding method；λgFor thermal coefficient；βwFor the intrinsic thermodynamic behaviour ginseng of workpiece material Number；reFor Grinding Contact effective radius；vsFor grinding wheel speed.
The grinding heat flow density q calculates as follows:
Wherein, vwTo transform into speed；S is instantaneous contact area area；" ± " number characterization is ground along inverse mode, reverse to be ground For "+".
Preferably, above-mentioned spatial complex curved surface grinding temperature Numerical Simulation Analysis method, in step S4, the grinding temperature Degree calculates as follows:
Wherein, α is principal direction λ1With the angle of reference axis X；L=(vw·a)/(2α)；Z=(vwZ)/(2 α), z are indicated Surface distance of the temp probe apart from workpiece.
Preferably, above-mentioned spatial complex curved surface grinding temperature Numerical Simulation Analysis method, in step S5, the temperature is repaired Positive coefficient calculates as follows:
Wherein, ψ cutting angle between processing curve and curved projection surfaces；
Modified grinding temperature TEIt calculates as follows:
In general, through the invention it is contemplated above technical scheme is compared with the prior art, can obtain down and show Beneficial effect:
(1) it is multiple to calculate space for a kind of spatial complex curved surface grinding temperature Numerical Simulation Analysis method provided by the invention The instantaneous contact ellipse of miscellaneous curved surface establishes spatial complex curved surface product grinding process parameters based on model of ellipse and temperature Three-dimensional numerical simulation analysis model can carry out the grinding temperature in the components grinding process with spatial complex curved surface feature Numerical Simulation calculates analysis, point of residual stress after product can be assisted to carry out intensity and thermal-structure coupled design, grinding Cloth analysis and grinding process parameters optimization design, improve the production efficiency of spatial complex curved surface product, shorten production life cycle；It starts Grinding temperature Numerical Simulation Analysis by two-dimensional surface Numerical Simulation Analysis enters the new office of three-dimensional space curved surface Numerical Simulation Analysis Face；
(2) a kind of spatial complex curved surface grinding temperature Numerical Simulation Analysis method provided by the invention introduces grinding temperature Degree value simulation calculation correction factor, had both had modified the theoretical error of one-dimensional, two-dimentional heat source model, also had modified three-dimensional elliptical song The analytical error of surface model, to improve grinding temperature Numerical Simulation Analysis precision, Numerical Simulation Analysis resultant error is small, point Analysis method high reliablity about improves 15% to the numerical simulation computational accuracy of grinding temperature.
Detailed description of the invention
Fig. 1 is the flow chart of spatial complex curved surface grinding temperature Numerical Simulation Analysis method provided in an embodiment of the present invention；
Fig. 2 is the curved surface per unit system arrow figure of workpiece to be processed provided in an embodiment of the present invention；
Fig. 3 is that instantaneous contact ellipse provided in an embodiment of the present invention calculates schematic diagram；
Fig. 4 is model of ellipse error analysis figure provided in an embodiment of the present invention.
Specific embodiment
In order to make the objectives, technical solutions, and advantages of the present invention clearer, with reference to the accompanying drawings and embodiments, right The present invention is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, and It is not used in the restriction present invention.As long as in addition, technical characteristic involved in the various embodiments of the present invention described below Not constituting a conflict with each other can be combined with each other.
A kind of spatial complex curved surface grinding temperature Numerical Simulation Analysis method provided by the present invention, for complicated to having Grinding temperature numerical value of the product of curved surface in grinding process carries out simulation calculation；It is space provided in this embodiment shown in Fig. 1 The flow chart of complex-curved grinding temperature Numerical Simulation Analysis method, as shown in Figure 1, method includes the following steps:
S1: according to the basic principle of curved surface wheel, the per unit system arrow n and method of the curved surface r=r (μ, υ) of workpiece to be processed are calculated Curvature kn；
Per unit system swears that the calculation formula of n is as follows:
N=rμ×rυ/|rμ×rυ| (1)
Wherein, rμIt is surface equation r to the first-order partial derivative of independent variable μ；rυIt is inclined for single order of the surface equation r to independent variable υ Derivative.
It is the curved surface per unit system arrow figure of workpiece to be processed provided in this embodiment shown in Fig. 2, as shown in Fig. 2, curved surface r=r The normal curvature k of any fixed point P (μ, υ) on (μ, υ) along direction λ=d μ/d υnIt may be expressed as:
Wherein, E, F, G are the first fundamental quantity of curved surface r；L, M, N are the second fundamental quantity of curved surface r.
For the fixed point P (μ, υ) on curved surface, the first fundamental quantity of curved surface E, F, G and second fundamental quantity L, M, N are normal Number；At this point, λ is variable, therefore knIt is considered as the function of λ.N is sweared according to the per unit system that the solution of formula (1) and formula (2) obtains And normal curvature kn, two principal curvatures of the fixed point P (μ, υ) on the spatial complex curved surface r=r (μ, υ) of workpiece to be processed can be calculated (kn1, kn2) and its corresponding two principal direction (λ1, λ2), it calculates as follows:
(E·M-F·L)·λ2λ+(FN-GM)=0 (4)+(EN-GL)
S2: according to two principal curvatures (k on spatial complex curved surface r=r (μ, υ)n1, kn2) and its corresponding two main sides To (λ1, λ2) and the first grinding process parameters, ellipsometric analysis model is established, the long axis a and short axle of instantaneous contact ellipse are calculated b；First grinding process parameters include the amount of collapsing and grinding depth；
Spatial complex curved surface r=r (μ, υ) includes two symmetroids, two principal curvatures kn1, kn2It is located at two In symmetroid, and two principal direction bending directions are on the contrary, two principal curve value near symmetricals；It, can be any in calculating process Specifying a certain principal direction is the first principal direction, and corresponding principal curvatures is the first principal curvatures；Another principal direction is the second principal direction, Principal curvatures is the second principal curvatures.
It is that instantaneous contact ellipse provided in this embodiment calculates schematic diagram shown in Fig. 3, is provided in this embodiment shown in Fig. 4 Model of ellipse error analysis figure；As shown in figure 4, the ellipsometric analysis model includes contact of the grinding medium plain emery wheel with curve surface work pieces Oval O, while also including projection ellipse O' and the projection ellipse O' quilt that Contact Ellipse O is obtained after plane projection The oval calotte load scope (i.e. projection ellipse O ") of processing curve formation is projected to again；Calculate oval calotte load work With the long axis a and short axle b of domain (i.e. projection ellipse O ").
Calculating process is as follows:
Wherein, δ is the amount of collapsing；apFor grinding depth；
A, B are principal curvatures and the constant that principal direction angle calcu-lation obtains；Specifically,
kn1、kn2The respectively principal curvatures of space curved surface and grinding wheel contact point, αsFor two principal direction λ1、λ2Between angle.
After the long axis a and short axle b that obtain instantaneous contact ellipse, the face of instantaneous contact ellipse is calculated according to long axis a and short axle b Product S, i.e. the instantaneous contact area area S of workpiece to be processed and grinding wheel.
S3: heat distribution coefficients R is calculated according to the second grinding process parameters in grinding processw, calculation formula is as follows:
Wherein, ξ is the characterization coefficient of grinding method；λgFor thermal coefficient；βwFor the intrinsic thermodynamic behaviour ginseng of workpiece material Number；reFor Grinding Contact effective radius；vsFor grinding wheel speed；
The characterization coefficient ξ of grinding method is related with the grinding method of grinding wheel, and the different grinding methods ξ value is different, such as: Characterization coefficient ξ=+ 1 when reverse dry grinding, characterization coefficient ξ=- 1 when forward dry grinding；Thermal coefficient λgIt is intrinsic with workpiece material Thermodynamic behaviour parameter betawDepending on the material of workpiece to be processed, different materials has different thermal coefficient λgAnd thermodynamics Characterisitic parameter βw.Due to being ground heat distribution coefficient RwInfluence factor in calculation formula is more, and the index variation range is usual Less, therefore empirical value, such as water base coolant liquid can be taken to be ground heat distribution coefficient R according to coolant liquid typewUsually take 0.65 ~0.70.
S4: according to heat distribution coefficients Rw, grinding wheel speed vs, transform into speed vw, tangential grinding force Ft and instantaneous contact area Area S calculates the grinding heat flow density q under model of ellipse, and calculation formula is as follows:
Wherein, vwTo transform into speed；S is instantaneous contact area area；" ± " number characterization is ground along inverse mode, reverse to be ground For "+".
Tangential grinding force FtIt can be solved according to the calculation formula under different situations, it is however generally that, tangential grinding force FtTool There are following three kinds of calculations: first is that the estimation formula established based on dimension analytic method；Second is that the warp established based on experimental data Test formula；Third is that combining the calculation formula established based on dimension analytic method and experimental data；The present embodiment preferably uses the third Calculation, specific calculation formula are as follows:Wherein, FpFor unit grinding force, x, y, z is that experience is public Formula index.
S5: grinding temperature T is calculated according to heat flow density q is ground in step S4, calculation formula is as follows:
Wherein, α is the angle of principal direction λ and reference axis X；L=(vwA)/(2 α) is dimensionless coordinate amount；For the unidirectional dimensionless coordinate amount of z, z indicates surface distance of the temp probe apart from workpiece.
S6: according to the analysis of calculation models error coefficient that cuts angle between processing curve and projection plane, with the error system Number carries out numerical simulation amendment to the grinding temperature T that step S5 is obtained as temperature correction coefficient κ, obtains modified grinding temperature TE；
Temperature correction coefficient κ calculates as follows:
Wherein, ψ cutting angle between processing curve and curved projection surfaces；
Modified grinding temperature TEIt calculates as follows:
One-dimensional heat source model is linear heat source, and two-dimentional heat source model is rectangular model, and heat source model is calculating instant contact Projection pattern is not used to be calculated when region area；The grinding temperature correction factor κ is from heat source model error and ellipsometric analysis Model error two parts are modified, and have both had modified the theoretical error of one-dimensional, two-dimentional heat source model, also have modified three-dimensional elliptical song The analytical error of surface model, to improve the reliability and accuracy of Numerical Simulation Analysis.
So far, by Contact Ellipse Numerical Simulation Analysis model, the grinding in spatial complex curved surface grinding has been acquired Temperature.
Compared to the existing grinding temperature Numerical Simulation Analysis method based on heat source model, space provided in this embodiment Complex-curved grinding temperature Numerical Simulation Analysis method, the instant contact of spatial complex curved surface is calculated in conjunction with space contact theory Ellipse establishes the three-dimensional numerical simulation analysis of the spatial complex curved surface product grinding process parameters based on model of ellipse and temperature Model can carry out Numerical Simulation calculating to the grinding temperature in the components grinding process with spatial complex curved surface feature Analysis, the distributional analysis of residual stress and grinding process after product can be assisted to carry out intensity and thermal-structure coupled design, grinding Parameters Optimal Design improves the production efficiency of spatial complex curved surface product, shortens production life cycle；And introduce grinding temperature numerical value Simulation calculation correction factor improves grinding temperature Numerical Simulation Analysis precision, and Numerical Simulation Analysis resultant error is small, analysis side Method high reliablity about improves 15% to the numerical simulation computational accuracy of grinding temperature；Meanwhile it is imitative to have started grinding temperature numerical value True analysis enters the three-dimensional space curved surface Numerical Simulation Analysis new situation by two-dimensional surface Numerical Simulation Analysis.
As it will be easily appreciated by one skilled in the art that the foregoing is merely illustrative of the preferred embodiments of the present invention, not to The limitation present invention, any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should all include Within protection scope of the present invention.

Claims (10)

1. a kind of spatial complex curved surface grinding temperature Numerical Simulation Analysis method, which comprises the following steps:
S1: the curved surface is calculated according to the per unit system of the curved surface of workpiece to be processed arrow and normal curvature and is taken up an official post two principal curvatures of a fixed point And its corresponding two principal direction；
S2: establishing ellipsometric analysis model according to two principal curvatures and its corresponding principal direction and the first grinding process parameters, meter Calculate long axis, short axle and the area of the instantaneous contact ellipse of workpiece to be processed and grinding wheel；First grinding process parameters include pressure Fall into amount and grinding depth；
S3: determining grinding method, according to the characterization coefficient of the grinding method, workpiece material inherent characteristic parameter, instant contact Elliptical area and the second grinding parameter calculate grinding heat flow density；Second grinding process parameters include that Grinding Contact has Effect radius, transforms into speed and tangential grinding force at grinding wheel speed；
S4: grinding temperature is calculated according to the grinding heat flow density.
2. spatial complex curved surface grinding temperature Numerical Simulation Analysis method as described in claim 1, which is characterized in that the wink When Contact Ellipse be that the Contact Ellipse of grinding process medium plain emery wheel and workpiece to be processed projects to processing again after plane projection Curved surface is formed by oval calotte load scope.
3. spatial complex curved surface grinding temperature Numerical Simulation Analysis method as claimed in claim 2, which is characterized in that step S4 It is further comprising the steps of later:
S5: calculating the error coefficient of model of ellipse according to cutting angle between processing curve and projection plane, with the error system Number carries out numerical simulation amendment to the grinding temperature as temperature correction coefficient, obtains modified grinding temperature.
4. spatial complex curved surface grinding temperature Numerical Simulation Analysis method as described in claim 1, which is characterized in that step S3 In include following sub-step:
S31: determining grinding method, according to the characterization coefficient of the grinding method and workpiece material inherent characteristic parameter, grinding Contact the heat distribution coefficient in effective radius and grinding wheel speed calculating grinding process；
S32: according to the heat distribution coefficient, grinding wheel speed, the area for transforming into speed, tangential grinding force and instantaneous contact ellipse Calculate grinding heat flow density.
5. spatial complex curved surface grinding temperature Numerical Simulation Analysis method as described in claim 1, which is characterized in that step S1 In, the curved surface of workpiece to be processed is taken up an official post two principal curvatures k of a fixed pointn1, kn2And its corresponding two principal direction λ1, λ2It calculates such as Under:
(E·M-F·L)·λ2+ (EN-GL) λ+(FN-GM)=0.
6. spatial complex curved surface grinding temperature Numerical Simulation Analysis method as claimed in claim 5, which is characterized in that step S1 In, the per unit system of the curved surface r=r (μ, υ) of workpiece to be processed swears n and normal curvature knIt calculates as follows:
N=rμ×rυ/|rμ×rυ|
Wherein, rμIt is surface equation r to the first-order partial derivative of independent variable μ；rυIt is surface equation r to the single order local derviation of independent variable υ Number；
The normal curvature k of any fixed point on curved surface r=r (μ, υ) along direction λ=d μ/d υnAre as follows:
Wherein, E, F, G are the first fundamental quantity of curved surface r；L, M, N are the second fundamental quantity of curved surface r.
7. spatial complex curved surface grinding temperature Numerical Simulation Analysis method as claimed in claim 5, which is characterized in that step S2 In, the long axis a and short axle b of the instantaneous contact ellipse calculate as follows:
Wherein, δ is the amount of collapsing；Ap is grinding depth；
kn1、kn2Respectively two principal curvatures of space curved surface and grinding wheel contact point, αsFor two principal direction λ1、λ2Between angle.
8. the spatial complex curved surface grinding temperature Numerical Simulation Analysis method as described in claim 4 or 7, which is characterized in that step In rapid S3, the heat distribution coefficients RwIt calculates as follows:
Wherein, ξ is the characterization coefficient of grinding method；λgFor thermal coefficient；βwFor workpiece material intrinsic heat mechanics parameters；re For Grinding Contact effective radius；Vs is grinding wheel speed.
The grinding heat flow density q calculates as follows:
Wherein, vwTo transform into speed；S is the area of instantaneous contact ellipse；" ± " number characterization is ground along inverse mode, and reverse grinding is “+”。
9. spatial complex curved surface grinding temperature Numerical Simulation Analysis method as claimed in claim 8, which is characterized in that step S4 Described in grinding temperature calculate it is as follows:
Wherein, α is the angle of principal direction λ and reference axis X；L=(vw·a)/(2α)；Z=(vwZ)/(2 α), z indicate that temperature is visited Surface distance of the head apart from workpiece.
10. spatial complex curved surface grinding temperature Numerical Simulation Analysis method as claimed in claim 9, which is characterized in that step Temperature correction coefficient described in S5 calculates as follows:
Wherein, ψ cutting angle between processing curve and curved projection surfaces；
Modified grinding temperature TEIt calculates as follows:
CN201811362258.0A 2018-11-15 2018-11-15 Numerical simulation analysis method for grinding temperature of spatial complex curved surface Active CN109635348B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811362258.0A CN109635348B (en) 2018-11-15 2018-11-15 Numerical simulation analysis method for grinding temperature of spatial complex curved surface

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811362258.0A CN109635348B (en) 2018-11-15 2018-11-15 Numerical simulation analysis method for grinding temperature of spatial complex curved surface

Publications (2)

Publication Number Publication Date
CN109635348A true CN109635348A (en) 2019-04-16
CN109635348B CN109635348B (en) 2023-04-07

Family

ID=66068098

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811362258.0A Active CN109635348B (en) 2018-11-15 2018-11-15 Numerical simulation analysis method for grinding temperature of spatial complex curved surface

Country Status (1)

Country Link
CN (1) CN109635348B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113486474A (en) * 2021-07-28 2021-10-08 大连理工大学 Method for designing shape of grinding wheel for grinding complex curved surface part
CN114036637A (en) * 2021-11-05 2022-02-11 沈阳飞机设计研究所扬州协同创新研究院有限公司 Rapid design method for active cooling channel on complex curved surface

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11272321A (en) * 1998-03-23 1999-10-08 Shigiya Machinery Works Ltd Grinding method of non-completely round work piece
CN107116404A (en) * 2017-04-09 2017-09-01 北京工业大学 A kind of constant speed bent axle follow-up grinding machining experiment system and process
CN108151885A (en) * 2017-12-15 2018-06-12 湖南科技大学 Non-circular profile Work piece high-speed grinding temperature Forecasting Methodology based on variable heat source model

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11272321A (en) * 1998-03-23 1999-10-08 Shigiya Machinery Works Ltd Grinding method of non-completely round work piece
CN107116404A (en) * 2017-04-09 2017-09-01 北京工业大学 A kind of constant speed bent axle follow-up grinding machining experiment system and process
CN108151885A (en) * 2017-12-15 2018-06-12 湖南科技大学 Non-circular profile Work piece high-speed grinding temperature Forecasting Methodology based on variable heat source model

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113486474A (en) * 2021-07-28 2021-10-08 大连理工大学 Method for designing shape of grinding wheel for grinding complex curved surface part
CN113486474B (en) * 2021-07-28 2024-06-11 大连理工大学 Shape design method of grinding wheel for grinding complex curved surface part
CN114036637A (en) * 2021-11-05 2022-02-11 沈阳飞机设计研究所扬州协同创新研究院有限公司 Rapid design method for active cooling channel on complex curved surface

Also Published As

Publication number Publication date
CN109635348B (en) 2023-04-07

Similar Documents

Publication Publication Date Title
CN104298818B (en) A kind of end mill processing surface error prediction and emulation mode
CN105117547B (en) The emulation mode of micro- milling nickel base superalloy prediction residue stress
Kaszynski et al. Uncertainties of an automated optical 3d geometry measurement, modeling, and analysis process for mistuned integrally bladed rotor reverse engineering
Huang et al. A systematic approach for online minimizing volume difference of multiple chambers in machining processes based on high-definition metrology
CN109635348A (en) A kind of spatial complex curved surface grinding temperature Numerical Simulation Analysis method
Li et al. Quadratic curve heat flux distribution model in the grinding zone
CN109783968A (en) The threedimensional FEM method of metal cutting process based on Simulation Based On Multi-step
CN108304687B (en) Method for predicting turning deformation of thin-wall complex curved surface rotating member
Deja et al. Application of reverse engineering technology in part design for shipbuilding industry
CN103366069A (en) Hierarchical algorithm of selective laser sintering
Wang et al. The temperature field study on the annular heat source model in large surface grinding by cup wheel
Guo et al. An adaptive sampling methodology for measuring blades with CMM based on dominant feature points
Liu et al. Numerical and experimental investigation on temperature field during belt grinding considering elastic contact
Xiao et al. A measurement method of the belt grinding allowance of hollow blades based on blue light scanning
Zhang et al. Understanding kinematics of the orthogonal cutting using digital image correlation—measurement and analysis
CN105956326B (en) Response surface model-based technological parameter optimization method for roller cooling system
CN117473802B (en) Method for rapidly solving polishing residence time of large-caliber optical element
Zhang et al. Adaptive NC path generation from massive point data with bounded error
JP3935278B2 (en) Method for manufacturing a blade of a fluid machine
CN110879926B (en) Method for cutting heat simulation calculation in brittle material cutting process
Wang et al. Optimization of parameters for bonnet polishing based on the minimum residual error method
Yu et al. Surface modeling method for aircraft engine blades by using speckle patterns based on the virtual stereo vision system
CN115964594B (en) Method and system for rapidly acquiring vertex curvature of grid model
Wan et al. Constructing process models of engine blade surfaces for their adaptive machining: an optimal approach
Zhou-jie et al. Recent progress on 3D shape and deformation measurement based on fringe projection

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant