CN106525418B - A kind of engaging portion stiffness analysis method for sided locking knife handle-main shaft under high speed conditions - Google Patents

A kind of engaging portion stiffness analysis method for sided locking knife handle-main shaft under high speed conditions Download PDF

Info

Publication number
CN106525418B
CN106525418B CN201610985606.4A CN201610985606A CN106525418B CN 106525418 B CN106525418 B CN 106525418B CN 201610985606 A CN201610985606 A CN 201610985606A CN 106525418 B CN106525418 B CN 106525418B
Authority
CN
China
Prior art keywords
rigidity
engaging portion
face
stiffness
faying face
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.)
Active
Application number
CN201610985606.4A
Other languages
Chinese (zh)
Other versions
CN106525418A (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.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
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 Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN201610985606.4A priority Critical patent/CN106525418B/en
Publication of CN106525418A publication Critical patent/CN106525418A/en
Application granted granted Critical
Publication of CN106525418B publication Critical patent/CN106525418B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • G01M13/02Gearings; Transmission mechanisms
    • G01M13/025Test-benches with rotational drive means and loading means; Load or drive simulation

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Prostheses (AREA)

Abstract

The engaging portion stiffness analysis method that the invention discloses a kind of for sided locking knife handle-main shaft under high speed conditions, more specifically a kind of bound fraction shape rigidity modeling method for considering double combination surface pressure uneven distribution, and system engaging portion rigidity under high speed conditions is analyzed based on this method.This method initially sets up three-dimensional fractal normal direction and shear stiffness model, and each node pressure values on the double combination face obtained based on static analysis, the corresponding equivalent stiffness value of calculate node, and obtains engaging portion torsion and radial rigidity value by transformation model.The final effect tendency for disclosing broach force under different rotating speeds and high revolving speed, disc spring rigidity with disc spring pretightning force to engaging portion rigidity.The characteristics of this method, considers elastic-plastic deformation in fractal theory and the influence of the factor is expanded in domain, it is considered in finite element analysis simultaneously and is unevenly distributed situation in conjunction with surface pressure, to carry out Accurate Model and analysis for BT40 knife handle-main shaft engaging portion rigidity under high speed conditions.

Description

A kind of engaging portion rigidity for sided locking knife handle-main shaft under high speed conditions Analysis method
Technical field
The invention belongs to knife handle-axis system dynamic analysis fields, are related to a kind of for sided locking knife handle-main shaft Engaging portion stiffness analysis method under high speed conditions, more specifically a kind of double combination surface pressure uneven distribution that considers Bound fraction shape rigidity modeling method, and system engaging portion rigidity under high speed conditions is analyzed based on this method.
Background technique
Knife handle is the important component of connecting machine tool main shaft and cutter, and main shaft-cutter handle combining part is the weak ring of the system One of section, engaging portion stiffness characteristics directly affect system machining accuracy.BT40 type knife handle-axis system is a kind of traditional list Face contact knife handle can generate local detachment in high speed conditions inferior pyramidal faying face, seriously affect its overall performance, sided locking Knife handle-axis system effectively raises the engaging portion rigidity under high speed conditions, improves system machining accuracy.Sided locking Knife handle is made of knife handle ontology, disc spring, tapered sleeve three parts, by applying broach force together with main shaft close-coupled, in height Disc spring can push tapered sleeve to adjust between its male cone (strobilus masculinus) and main shaft inner conical surface and knife handle end face and main shaft end face automatically under speed conditions Between night distribution, fitting closely for double combination face is remained, so as to improve overall dynamics performance.But it is being higher than certain turn Big end local detachment still can occur for taper faying face when the fast upper limit, to reduce engaging portion rigidity, influence dynamic performance. It is that the interference effect of the external environments such as noise in Experimental Identification method is affected to rigidity estimation results due to traditional, and the party Method is only limitted to analysis without knife handle under speed conditions-axis system engaging portion rigidity Research.In order to knife handle main shaft engaging portion rigidity It is analyzed, fractal theory provides a kind of effective way with the mode that emulation combines for the engaging portion rigidity modeling.Divide shape Theory is a kind of effective modeling method for analyzing faying face, in this method faying face characterization parameter have Scaling, not by The advantages such as instrumental resolution and sample length influence, are widely used in and estimate the exposure parameters such as joint surface contact stiffness damping. The three-dimensional fractal model that considers in the present invention while considering elastic-plastic deformation and the influence of factor pair point shape rigidity is expanded in domain, and Three-Dimensional contact normal direction and shear stiffness model are established based on this.By carrying out three-dimensional modeling and finite element to sided locking knife handle Static analysis obtains the non-uniform node pressure distribution situation in double combination face, and three-dimensional fractal model is combined with finite element analysis Mode can effectively improve the estimate accuracy of engaging portion rigidity.It is multiple that this method can be used for solving current sided locking knife handle structure It is miscellaneous, it is difficult to which that the bottleneck for analyzing engaging portion stiffness analysis realizes sided locking knife handle-axis system engaging portion under high speed conditions Rigidity modeling analysis.
Summary of the invention
The present invention is intended to provide a kind of engaging portion stiffness analysis for sided locking knife handle-main shaft under high speed conditions Method.This method is mainly characterized by using in conjunction with three-dimensional fractal microscopic contact rigidity modeling theory and Finite Element Static Analysis Method models the sided locking knife handle under high speed conditions-axis system engaging portion torsion with radial rigidity, and is divided Analysis, becomes to disclose the influence of different technical parameters (revolving speed, broach force, disc spring rigidity and disc spring pretightning force) to engaging portion rigidity Gesture.
The present invention is realized using following technological means:
S1, consider that elastic-plastic deformation and domain expand the influence of the factor simultaneously, obtained accurate three-dimensional fractal theory, and base Contact normal direction and shear stiffness model are established in the theory.
S2, to sided locking knife handle-axis system 3-D geometric model, and add osculating element, carry out finite element static(al) Each node pressure values in double combination face are extracted in analysis, pass through the faying face of regional area where obtaining each node multiplied by grid area Pressure value.
S3, it is based on three-dimensional fractal theory, it is corresponding to calculate each node using node regional area faying face pressure value as input Faying face contact normal direction and shear stiffness value.
S4, the modulus of conversion that each node corresponds to normal stiffness, shear stiffness and double combination face torsion stiffness, radial rigidity is established Type calculates engaging portion torsion and radial rigidity value.
It is S5, corresponding by seeking under different rotating speeds and high speed conditions different broach forces, disc spring pretightning force, disc spring rigidity Engaging portion rigidity value discloses revolving speed, broach force, disc spring rigidity and disc spring pretightning force to the effect tendency of engaging portion rigidity.
Method of the features of the present invention in conjunction with three-dimensional fractal microscopic contact rigidity modeling theory and Finite Element Static Analysis, Sided locking knife handle under high speed conditions-axis system engaging portion torsion is modeled with radial rigidity, to disclose different works Skill parameter (revolving speed, broach force, disc spring rigidity and disc spring pretightning force) is to the effect tendency of engaging portion rigidity.
Detailed description of the invention
The torsion of the engaging portion Fig. 1 and radial rigidity schematic diagram;
Effect tendency of Fig. 2 different rotating speeds to engaging portion rigidity;
Effect tendency of the broach force to engaging portion rigidity under Fig. 3 25000r/min high speed conditions;
Effect tendency of the disc spring rigidity to engaging portion rigidity under Fig. 4 25000r/min high speed conditions;
Effect tendency of the disc spring pretightning force to engaging portion rigidity under Fig. 5 25000r/min high speed conditions;
Specific embodiment
Below in conjunction with attached drawing 1-5, invention is further described in detail.
Step (1) establishes three-dimensional normal direction and shear stiffness model;
Based on M-B fractal theory, in conjunction with Hertz theory, while considering that elastic-plastic deformation and domain expand factor ψ, by place Single micro-bulge normal load in different distortion region, which carries out integral, can be obtained proof resilience normal load, elastoplasticity normal load And total plastic property normal load difference is as follows:
Wherein HG1,HG2It is coefficient relevant to material properties and faying face fractal parameter,
H be compared with The hardness of soft material, H=2.8Y (Y is yield strength value);K is parameter relevant to Poisson's ratio, k=0.454+0.41 ν;a′1c, a′2cCritical cross-section product respectively between elasticity, elastoplasticity and plastic deformation.
Then the total normal load of faying face is represented by F=Fe+Fep+Fp
In normal stiffness modeling, it is respectively as follows: in flexible deformation and elastic-plastic deformation region normal contact stiffness
In shear stiffness modeling, it is respectively as follows: in flexible deformation and the tangential contact stiffness in elastic-plastic deformation region
G' is engaging portion equivalent shear modulus, 1/G'=(2- ν in formula1)/G1+(2-ν2)/G2;H1,H2Pass through formulaIt obtains.
To sum up, faying face normal direction and tangential total contact stiffness are respectively as follows: Kn=Kne+Knep, Kt=Kte+Ktep
Step (2) Finite Element Static Analysis;
3-D geometric model is established to sided locking BTF40 knife handle-axis system, by using TARGE object element with CONTAC osculating element establishes double combination face contact pair, using mapping mode to two-sided knife handle grid division, to main shaft into The constraint of row axial restraint applies broach force in knife handle small end, integrally applies revolving speed constraint to system, carry out static analysis.Dividing The node pressure values on conical surface faying face and end face faying face are extracted in analysis result respectively, and are calculated in grid institute where each node In regional area faying face pressure value.
Step (3) calculates the corresponding equivalent normal direction of each node and shear stiffness on double combination face;
By partial pressure value where the conical surface and end face faying face node substitute into faying face fractal method to shear stiffness model In, calculate equivalent normal direction and shear stiffness value corresponding to each node on double combination face.
Step (4) calculations incorporated portion radial direction and torsion stiffness;
Since taper faying face and end face faying face are rotationally symmetrical surf, the radial direction established respectively in the y-direction is rigid It spends and as follows around the torsion stiffness model of y-axis:
N in formulaT, nERespectively taper-face contact face and end face contact surface interstitial content;Respectively conical surface faying face Node normal stiffness and shear stiffness;Respectively end face faying face normal stiffness and shear stiffness;θ is contact surface section Point position angle;Respectively conical surface faying face and end face faying face node between x-axis at a distance from.
Step (5) discloses the influence of revolving speed and broach force to engaging portion rigidity;
Engaging portion is reversed to become with radial rigidity value, influence of the announcement revolving speed to engaging portion rigidity under the conditions of calculating different rotating speeds Gesture;Different broach forces under 25000r/min high speed conditions, different disc spring rigidity and different disc spring pretightning forces are calculated to engaging portion The influence of torsion and radial rigidity, and make tendency chart.

Claims (1)

1. a kind of engaging portion stiffness analysis method for sided locking knife handle-main shaft under high speed conditions, this method is used In conjunction with the method for three-dimensional fractal microscopic contact rigidity modeling theory and Finite Element Static Analysis, to the double-sided lock under high speed conditions Tight knife handle-axis system engaging portion torsion is modeled with radial rigidity, and is analyzed, to disclose different rotating speeds, broach force, dish The effect tendency of spring rigidity and disc spring pretightning force technological parameter to engaging portion rigidity;
It is characterized by:
S1, consider that elastic-plastic deformation and domain expand the influence of the factor simultaneously, it is theoretical to have obtained accurate three-dimensional fractal, and being based on should Theory establishes contact normal direction and shear stiffness model;
S2,3-D geometric model is established to sided locking knife handle-axis system, and adds osculating element, carry out finite element static(al) Each node pressure values in double combination face are extracted in analysis, pass through the faying face of regional area where obtaining each node multiplied by grid area Pressure value;
S3, it is based on three-dimensional fractal microscopic contact rigidity modeling theory, using node regional area faying face pressure value as in terms of inputting Calculate the corresponding faying face contact normal direction of each node and shear stiffness value;
S4, the transformation model that each node corresponds to normal stiffness, shear stiffness and double combination face torsion stiffness, radial rigidity is established, Calculate engaging portion torsion and radial rigidity value;
S5, by seeking under different rotating speeds and high speed conditions different broach forces, disc spring pretightning force, the corresponding combination of disc spring rigidity Portion's rigidity value discloses revolving speed, broach force, disc spring rigidity and disc spring pretightning force to the effect tendency of engaging portion rigidity;
Step (1) establishes three-dimensional normal direction and shear stiffness model;
Based on M-B fractal theory, in conjunction with Hertz theory, while considering that elastic-plastic deformation and domain expand factor ψ, by not Single micro-bulge normal load with deformed region carries out integral proof resilience normal load, elastoplasticity normal load and total can be obtained Plasticity normal load difference is as follows:
Wherein HG1,HG2It is coefficient relevant to material properties and faying face fractal parameter,
H is The hardness of softer material, H=2.8Y, Y are yield strength value;K is parameter relevant to Poisson's ratio, k=0.454+0.41 ν;Then The total normal load of faying face is expressed as F=Fe+Fep+Fp
In normal stiffness modeling, it is respectively as follows: in flexible deformation and elastic-plastic deformation region normal contact stiffness
In shear stiffness modeling, it is respectively as follows: in flexible deformation and the tangential contact stiffness in elastic-plastic deformation region
G' is engaging portion equivalent shear modulus, 1/G'=(2- ν in formula1)/G1+(2-ν2)/G2;H1,H2Pass through formulaIt obtains;
To sum up, faying face normal direction and tangential total contact stiffness are respectively as follows: Kn=Kne+Knep, Kt=Kte+Ktep;Step (2) finite element Static analysis;
3-D geometric model is established to sided locking knife handle-axis system, is contacted by using TARGE object element with CONTAC Unit establishes double combination face contact pair, using mapping mode to two-sided knife handle grid division, carries out axial restraint constraint to main shaft, Apply broach force in knife handle small end, integrally applies revolving speed constraint to system, carry out static analysis;It is extracted respectively in analysis result Node pressure values on conical surface faying face and end face faying face, and calculate in regional area faying face where grid where each node Pressure value;
Step (3) calculates the corresponding equivalent normal direction of each node and shear stiffness on double combination face;
Partial pressure value where the conical surface and end face faying face node is substituted into faying face fractal method to shear stiffness model, is counted Calculate equivalent normal direction and shear stiffness value corresponding to each node on double combination face;
Step (4) calculations incorporated portion radial direction and torsion stiffness;
Since taper faying face and end face faying face are rotationally symmetrical surf, establish respectively radial rigidity in the y-direction and Torsion stiffness model around y-axis is as follows:
N in formulaT, nERespectively taper-face contact face and end face contact surface interstitial content;Respectively conical surface faying face nodal method To rigidity and shear stiffness;Respectively end face faying face normal stiffness and shear stiffness;θ is contact surface node location Angle;Respectively conical surface faying face and end face faying face node between x-axis at a distance from;
Step (5) discloses the influence of revolving speed and broach force to engaging portion rigidity;
Engaging portion is reversed and radial rigidity value, effect tendency of the announcement revolving speed to engaging portion rigidity under the conditions of calculating different rotating speeds.
CN201610985606.4A 2016-11-09 2016-11-09 A kind of engaging portion stiffness analysis method for sided locking knife handle-main shaft under high speed conditions Active CN106525418B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610985606.4A CN106525418B (en) 2016-11-09 2016-11-09 A kind of engaging portion stiffness analysis method for sided locking knife handle-main shaft under high speed conditions

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610985606.4A CN106525418B (en) 2016-11-09 2016-11-09 A kind of engaging portion stiffness analysis method for sided locking knife handle-main shaft under high speed conditions

Publications (2)

Publication Number Publication Date
CN106525418A CN106525418A (en) 2017-03-22
CN106525418B true CN106525418B (en) 2019-09-27

Family

ID=58350542

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610985606.4A Active CN106525418B (en) 2016-11-09 2016-11-09 A kind of engaging portion stiffness analysis method for sided locking knife handle-main shaft under high speed conditions

Country Status (1)

Country Link
CN (1) CN106525418B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108303251B (en) * 2018-04-19 2020-03-06 清华大学 Rigidity modeling and indirect detection method in electric spindle rotation state
CN108628250B (en) * 2018-06-01 2020-03-31 西安交通大学 Optimization method for matching structure size of machine tool spindle and tool shank
CN108760275B (en) * 2018-06-29 2024-04-05 上海第二工业大学 Device and method for analyzing static rigidity of combination part of cutter, cutter handle and main shaft system

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102592017A (en) * 2011-12-31 2012-07-18 北京工业大学 Two-sided locking knife handle/main shaft coupling performance simulating and optimizing method
CN104294745A (en) * 2014-10-16 2015-01-21 武汉科技大学 Detection method of skid resistance of asphalt pavement based on three-dimensional fractal dimension

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5466451B2 (en) * 2009-07-31 2014-04-09 株式会社豊田中央研究所 Numerical analysis device, numerical analysis method, program, and test method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102592017A (en) * 2011-12-31 2012-07-18 北京工业大学 Two-sided locking knife handle/main shaft coupling performance simulating and optimizing method
CN104294745A (en) * 2014-10-16 2015-01-21 武汉科技大学 Detection method of skid resistance of asphalt pavement based on three-dimensional fractal dimension

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
依据各向异性分形几何理论的固定结合部法向接触力学模型;田红亮等;《机械工程学报》;20131130;第49卷(第21期);全文 *
高转速下主轴_刀柄结合部接触刚度分析;赵永胜等;《华中科技大学学报(自然科学版)》;20160731;第44卷(第7期);论文第92-93页第1-2节 *

Also Published As

Publication number Publication date
CN106525418A (en) 2017-03-22

Similar Documents

Publication Publication Date Title
CN106525418B (en) A kind of engaging portion stiffness analysis method for sided locking knife handle-main shaft under high speed conditions
Piomelli et al. Turbulent structures in accelerating boundary layers
Wang et al. A study of effects of roller path profiles on tool forces and part wall thickness variation in conventional metal spinning
Jeswiet et al. Single point and asymmetric incremental forming
Houck III et al. A tuned holder for increased boring bar dynamic stiffness
US10140416B2 (en) Process-orientated design method for machine tool structures
Hussain et al. Guidelines for tool-size selection for single-point incremental forming of an aerospace alloy
CN106802630B (en) The Parallel frequency domain method of milling parameter stability prediction based on thin-wall part
Kim et al. Finite element inverse analysis for the design of intermediate dies in multi-stage deep-drawing processes with large aspect ratio
Hirt et al. Rolling of functional metallic surface structures
WO2018187291A3 (en) Topology optimization on the cloud for computer aided design
Lang et al. Study on hydromechanical deep drawing with uniform pressure onto the blank
ATE551174T1 (en) METHOD FOR PRODUCING A COMPONENT
CN106777452A (en) A kind of sided locking knife handle axis system joint portion stiffness characteristics optimization method based on particle cluster algorithm
Ucun et al. 3D numerical modelling of micro-milling process of Ti6Al4V alloy and experimental validation
Suresh et al. Spatially variant microstructured adhesive with one-way friction
Huang et al. Occupant comfort evaluation and wind-induced serviceability design optimization of tall buildings
CN106777463A (en) A kind of single-contact handle of a knife main shaft joint portion stiffness analysis method based on three-dimensional modification fractal theory
CN104036079A (en) Quantitative type static pressure turntable optimization design method based on sensitivity analysis
Bajpai et al. Orthogonal micro-grooving of anisotropic pyrolytic carbon
MX2017009261A (en) Method for estimating surface shape defect generating regions, method for estimating surface shape defect source regions, surface shape defect generating region estimating device, shape defect source region estimating device, program and recording media.
CN108052756A (en) A kind of method that faying face exposure parameter is determined based on FFT
Luo et al. Micro-hydromechanical deep drawing of metal cups with hydraulic pressure effects
Guo et al. A frequency-domain grinding force model-based approach to evaluate the dynamic performance of high-speed grinding machine tools
CN114722677A (en) Method for calculating natural frequency of bolted structure under normal and tangential load coupling effect

Legal Events

Date Code Title Description
C06 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