CN104166747A - Bolting joint part dynamic characteristic analysis method taking surface machining quality into consideration - Google Patents
Bolting joint part dynamic characteristic analysis method taking surface machining quality into consideration Download PDFInfo
- Publication number
- CN104166747A CN104166747A CN201410178375.7A CN201410178375A CN104166747A CN 104166747 A CN104166747 A CN 104166747A CN 201410178375 A CN201410178375 A CN 201410178375A CN 104166747 A CN104166747 A CN 104166747A
- Authority
- CN
- China
- Prior art keywords
- contact
- analysis
- model
- damping
- modeling
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 52
- 238000003754 machining Methods 0.000 title abstract description 7
- 238000013016 damping Methods 0.000 claims abstract description 34
- 238000000034 method Methods 0.000 claims abstract description 20
- 230000003068 static effect Effects 0.000 claims abstract description 14
- 239000000463 material Substances 0.000 claims abstract description 13
- 238000005316 response function Methods 0.000 claims abstract description 5
- 238000012545 processing Methods 0.000 claims description 15
- 239000011159 matrix material Substances 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000013461 design Methods 0.000 claims description 3
- 238000012217 deletion Methods 0.000 claims description 2
- 230000037430 deletion Effects 0.000 claims description 2
- 239000007787 solid Substances 0.000 claims description 2
- 238000012916 structural analysis Methods 0.000 claims description 2
- 238000010408 sweeping Methods 0.000 claims description 2
- 238000011161 development Methods 0.000 abstract description 4
- 230000001050 lubricating effect Effects 0.000 description 6
- 238000004088 simulation Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000012360 testing method Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 2
- 238000000227 grinding Methods 0.000 description 2
- 238000003801 milling Methods 0.000 description 2
- 208000032443 Masked facies Diseases 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a bolting joint part dynamic characteristic analysis method taking the surface machining quality into consideration and belongs to the field of mechanical property analysis. The method comprises the steps of parametric modeling, common joint surface modeling, boundary condition applying, static force analyzing, contact stiffness damping calculating, joint surface modeling reconstruction, bolting joint surface stiffness damping modeling and dynamic characteristic analyzing; automatic modeling is achieved according to model parameters, after static force analysis is selected, by means of the selected material parameters, the contact face machining method, the roughness analysis result and the static force analysis result, the joint face contact stiffness and the damping value are calculated, a contact unit is reconstructed, and inherent frequency, main vibration shapes and frequency-response functions of a joint part are obtained by dynamically analyzing the bolting joint part. According to the method, based on an ANSYS secondary development platform, a parameterized model of the bolting joint part is constructed, and the method is used for analyzing the influence on the dynamic characteristics of the bolting joint part from material attributes, the contact surface machining method, the machining roughness and external loads.
Description
Technical field
The invention belongs to software application technology field, particularly relate to a kind of bolted joint Dynamic Characteristics Analysis Method of considering surface of contact job operation and processing roughness.
Background technology
In lathe, have a large amount of bolted joint, in the lathe course of work, the dynamic perfromance of bolted joint affects the working (machining) efficiency of lathe, machining precision, the wearing and tearing between acceleration components, and then affect serviceable life and the lathe reliability of lathe.Therefore the research of bolted joint is very meaningful to improving machine mass.
The dynamic perfromance of bolted joint is subject to the impact of material properties, surface quality, external load very large.Can analysis of material attribute and the impact of external load on bolted joint dynamic perfromance by existing finite element analysis software, but the impact of surface of contact surface quality on Structure dynamic characteristics but can not be analyzed.
Chinese scholars has been set up contact stiffness and the contact damping model of considering surface of contact surface quality by Fractal Contact theory at present.The contact stiffness and the contact damping model that the present invention is based on Fractal Contact theory, utilize finite element software ANSYS secondary developing platform, sets up the bolted joint parameterized model of considering surface of contact job operation and processing roughness, for analyzing its dynamic perfromance.
Summary of the invention
The object of the invention is to, overcome the defect of existing faying face emulation mode, a kind of contact stiffness and contact damping model based on Fractal Contact theory is provided, utilize finite element software ANSYS secondary developing platform, set up the bolted joint parameterized model of considering surface of contact job operation and processing roughness, for analyzing its dynamic perfromance
The present invention adopts following technological means to realize, and the system that realizes the method comprises
A bolt faying face dynamic emulation method of considering surface of contact job operation and roughness, comprising: parametric modeling, the common modeling of faying face, apply boundary condition, static analysis, contact stiffness and contact damping calculating, faying face Model Reconstruction, the modeling of bolt faying face stiffness and damping, dynamic analysis step;
Parametric modeling step, utilizes body order to generate cube and column structure; Adopt Boolean subtraction calculation and bonding operation to generate bolted joint model; Generate after solid model, adopt gauze lattice control command and the order of grid sweeping, generate finite element grid model; In modeling process, utilize APDL language definition model parameterization variable, generate parameterized model;
The common modeling procedure of faying face, selects faying face by Parametric Design Language; Target face unit is set as the unit into TARGE170, and surface of contact unit is set to the unit into CONTA174; Set soft-soft contact, initial contact is binding;
Apply boundary condition step, comprise bolted joint constrained, load applying position, magnitude of load, loading direction; One section, bolt hole place is set to fixed constraint, and another section applies pressure load; Load applying region utilizes coordinate method to pick up node;
Static analysis step, first analysis type is set is static structural analysis, arranges and solves content, solves; By selecting osculating element, all nodes of selected cell, by loop command, extract node serial number and corresponding pressure values separately with array form according to order from small to large.
Contact stiffness and contact damping calculation procedure, according to the faying face node pressure values of extracting, by Fractal Contact rigidity and contact damping model, used loop command and formula operation, calculates contact stiffness and the contact damping value of each Nodes, with array form, preserves.
Faying face Model Reconstruction, by element deletion order, delete the contact model of setting up in the common modeling procedure of faying face, by loop command, create self-defined matrix unit, setup unit numbering, contact stiffness and contact damping value assignment are to unit real constant, generate contact stiffness and contact damping unit, according to node serial number, sequentially set up from small to large joint surface contact stiffness and and damping unit.
Dynamic analysis step, arranges analysis type, comprises model analysis, humorous response analysis; In model analysis aftertreatment, export natural frequency and the main shape of shaking; In humorous response analysis aftertreatment, export key point frequency response function figure.
According to the above results, user can analyze surface of contact machine-tooled method and the impact of processing roughness on bolted joint dynamic perfromance.
The aforesaid APDL of utilization language definition structural parameters variable comprises: bolted joint structural parameters, material number, bolt pretightening, faying face job operation, processing roughness; Wherein bolted joint structural parameters comprise total length, land length width, overall height, boss height, diameter of bolt hole; Material parameter comprises QT-235, QT250; Faying face job operation comprises milling, grinding; Processing roughness parameter comprises five roughness values (0.4,0.6,0.8,1.2,1.6), and unit is micron.
Aforesaid self-defined matrix unit is MATRIX27 matrix unit.
A kind of method of considering the analysis bolted joint dynamic perfromance of surface of contact job operation and processing roughness of the present invention, compared with prior art, has following obvious advantage and beneficial effect:
The present invention is based on contact stiffness and the contact damping model of Fractal Contact theory, utilize finite element software ANSYS secondary developing platform, set up to consider the bolted joint parameterized model of surface of contact job operation and processing roughness, for analysis of material attribute, surface of contact job operation, surfaceness, external load, impact in conjunction with shape on bolted joint dynamic perfromance.
Accompanying drawing explanation
Fig. 1 is bolted joint structural representation;
Fig. 2 is test specimen physical dimension schematic diagram;
Fig. 3 is bolted joint dynamic analysis process flow diagram;
Fig. 4 is MATRIX27 cell geometry sketch;
Fig. 5 is MATRIX27 cell matrix schematic diagram;
Fig. 6 is MATRIX27 matrix unit coefficient schematic diagram.
In figure, 1, upper test specimen, 2, joint portion osculating element, 3, lower test specimen.
Embodiment
Below in conjunction with Figure of description, specific embodiment of the invention is illustrated.
Referring to shown in Fig. 1, is bolted joint structural representation.Wherein, 1 is upper test specimen, and 2 is joint portion osculating element, and 3 is lower test specimen.
Below in conjunction with Fig. 2, Fig. 3, Fig. 4 further illustrates the present invention.
Wherein, L1 is total length, and L2 is land length, and D is width, and H1 is overall height, and H2 is boss height, and d is diameter of bolt hole.
(1) parametric modeling module, based on ANSYS secondary development, utilizes body order to generate cube and column structure; Adopt Boolean subtraction calculation and bonding operation to generate bolted joint model; Input parameter comprises bolted joint structural parameters (total length, land length width, overall height, boss height, diameter of bolt hole), load parameter (imposed load area, magnitude of load); Select parameter to comprise material (QT-235, QT250), job operation (milling, grinding), processing roughness value (0.4,0.6,0.8,1.2,1.6); After setting completes, automatically generate finite element grid model.
(2) the common MBM of faying face, based on ANSYS secondary development, setting up target face unit is TARGE170 unit, and surface of contact unit is CONTA174 unit, sets soft-soft contact, and initial contact is binding.
(3) apply boundary condition module, based on ANSYS secondary development, automatically according to input load parameter, adopt coordinate position to pick up order and determine that boundary condition applies region, load applying region and load value.
(4) static analysis module, choice structure static analysis, program is carried out static analysis calculating automatically, utilizes loop command to extract faying face node and corresponding pressure, and shows with text and picture form.
(5) contact stiffness Damping calculating module, program calculates in conjunction with material, job operation, the roughness selected the ct clamping value that each node is corresponding by combination node and the corresponding pressure extracting automatically automatically, with array form, preserves.
(6) faying face Model Reconstruction module (shown in Fig. 1), utilize APDL language to delete original faying face osculating element, ct clamping value according to calculating, creates by MATRIX27 unit faying face joint surface contact stiffness and the contact damping unit that each node is corresponding.
(7) dynamic analysis module, arranges analysis type, comprises the harmonious response analysis of model analysis, output modalities frequency and Mode Shape figure in model analysis aftertreatment; In humorous response analysis, also need further to define method for solving and load step option, in aftertreatment, export key point frequency response function figure.
User can select material, surface of contact to connect job operation and roughness by varying sized size, magnitude of load, analysis of material attribute, surface of contact crudy, the impact of external load on bolted joint dynamic perfromance.
Described parametric modeling, system, according to the structural parameters of user's input, is revised existing Parametric Design Language command stream, generates bolted joint model, then automatically carries out grid division.
Described lubricating pad faying face modeling, system, according to the lubricating pad parameter of input, adopt Do-Loop circulation to generate each lubricating pad self geometric center node and lubricating pad plane, and circulation generates supporting lubricating pad faying face and pretension lubricating pad faying face.
Described boundary condition, system, according to the boundary condition parameter of input, is revised APDL command stream, automatically boundary loading condition.
Described ct clamping calculates, and the node pressure that system obtains according to the material of selecting, job operation, roughness and static analysis calculates contact stiffness and the contact damping that each node is corresponding automatically.
Described dynamic analysis, comprises the harmonious response analysis of model analysis, obtains bolted joint natural frequency, principal mode and key point frequency response function.
In system faying face Model Reconstruction module, use self-defined matrix unit matrix27 to build faying face rigidity and damping unit.In ANSYS program, matrix27 unit can represent a kind of any unit.This unit is defined by two nodes and matrix of coefficients.Fig. 4 is MATRIX27 geometric representation.Wherein each node has 6 degree of freedom, is respectively along the translation of nodal coordinate system x, y, z direction with around the rotation of nodal coordinate system x, y, z direction.The geometrical property of this unit is without definition, but its elastokinematics response can specify with rigidity, damping or quality coefficient, and the constant in rigidity, damping and mass matrix can be used as the real constant input of unit.Fig. 6 is MATRIX27 matrix unit coefficient schematic diagram.When analog equivalent rigidity unit, C
1, C
58, C
7simulation directions X rigidity, C
18, C
64, C
19the rigidity of simulation Y-direction, C
24, C
69, C
60the rigidity of simulation Z direction.When analog equivalent damping unit, C
1, C
58, C
7the damping of simulation directions X, C
13, C
64, C
19the damping of simulation Y-direction, C
24, C
69, C
60the damping of simulation Z direction.
Claims (3)
1. a bolted joint Dynamic Characteristics Analysis Method of considering suface processing quality, comprising: parametric modeling, the common modeling of faying face, apply boundary condition, static analysis, contact stiffness Damping calculating, faying face Model Reconstruction, the modeling of faying face stiffness and damping, dynamic analysis step; It is characterized in that:
Described parametric modeling step, utilizes body order to generate cube and column structure; Adopt Boolean subtraction calculation and bonding operation to generate bolted joint model; Generate after solid model, adopt gauze lattice control command and the order of grid sweeping, generate finite element grid model; In modeling process, utilize APDL language definition model parameterization variable, generate parameterized model;
The described common modeling procedure of the face that closes, sets up the initial Binding Model of bolt faying face by Parametric Design Language.Target face unit is TARGE170 unit, soft-soft contact that surface of contact unit is that CONTA174 unit carries out;
The described boundary condition step that applies, comprises bolted joint constrained, load applying position, loading, magnitude of load, loading direction; Wherein load applying position and magnitude of load are defined as input variable;
Described static analysis step, first analysis type is set is static structural analysis, arranges and solves content, solves; By selecting osculating element, all nodes of selected cell, by loop command, extract node serial number and corresponding pressure values separately with array form according to order from small to large;
Described contact stiffness Damping calculating step, according to the faying face node pressure values of extracting, asks for faying face node contact stiffness and contact damping value by Fractal Contact rigidity and contact damping model, with array form, preserves;
Described faying face Model Reconstruction, by element deletion order, delete the contact model of setting up in the common modeling procedure of faying face, by loop command, create self-defined matrix unit, setup unit numbering, the node rigidity that contact stiffness Damping calculating step is obtained and damping value assignment are to unit real constant, generate contact stiffness and contact damping unit, according to node serial number, sequentially set up from small to large joint surface contact stiffness and damping unit;
Described dynamic analysis step, arranges analysis type, comprises model analysis, humorous response analysis; In model analysis aftertreatment, export natural frequency and the main shape of shaking; In humorous response analysis aftertreatment, export key point frequency response function figure.According to the above results, user can analyze surface of contact machine-tooled method and the impact of processing roughness on bolted joint dynamic perfromance.
2. a kind of bolted joint Dynamic Characteristics Analysis Method of considering suface processing quality according to claim 1, is characterized in that: by input load value, selection material, selection surface of contact job operation and processing roughness, can carry out bolted joint dynamic analysis.
3. a kind of bolted joint Dynamic Characteristics Analysis Method of considering suface processing quality according to claim 1, it is characterized in that: carry out static analysis and extract each node pressure, by pressure, calculate each node ct clamping, by MATRIX27, set up faying face osculating element for bolted joint dynamic analysis.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410178375.7A CN104166747B (en) | 2014-04-29 | 2014-04-29 | A kind of bolted joint Dynamic Characteristics Analysis Method for considering suface processing quality |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410178375.7A CN104166747B (en) | 2014-04-29 | 2014-04-29 | A kind of bolted joint Dynamic Characteristics Analysis Method for considering suface processing quality |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104166747A true CN104166747A (en) | 2014-11-26 |
| CN104166747B CN104166747B (en) | 2017-06-16 |
Family
ID=51910559
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410178375.7A Expired - Fee Related CN104166747B (en) | 2014-04-29 | 2014-04-29 | A kind of bolted joint Dynamic Characteristics Analysis Method for considering suface processing quality |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104166747B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106529035A (en) * | 2016-11-09 | 2017-03-22 | 北京工业大学 | Method for analyzing dynamic characteristics of bolted structure based on three-dimensional fractal theory |
| CN106815407A (en) * | 2016-12-22 | 2017-06-09 | 四川大学 | A kind of whole machine dynamic property optimization method of Digit Control Machine Tool |
| CN107679301A (en) * | 2017-09-20 | 2018-02-09 | 西安建筑科技大学 | A kind of segmented heavy duty crossbeam scale model design method |
| CN108563850A (en) * | 2018-03-29 | 2018-09-21 | 内蒙古久和能源装备有限公司 | A kind of Random Response Analysis method and system considering material damping otherness |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108628250B (en) * | 2018-06-01 | 2020-03-31 | 西安交通大学 | Optimization method for matching structure size of machine tool spindle and tool shank |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5377116A (en) * | 1991-07-01 | 1994-12-27 | Valenite Inc. | Method and system for designing a cutting tool |
| CN102063552A (en) * | 2011-01-07 | 2011-05-18 | 西安交通大学 | Method for determining rigidity of bolt connecting piece |
| CN102592017A (en) * | 2011-12-31 | 2012-07-18 | 北京工业大学 | Two-sided locking knife handle/main shaft coupling performance simulating and optimizing method |
| CN103186698A (en) * | 2013-04-12 | 2013-07-03 | 北京工业大学 | Simulation optimization method of dynamic and static performances of static pressure turntable of heavy machine tool |
| CN103366050A (en) * | 2013-06-25 | 2013-10-23 | 华中科技大学 | High-speed electric spindle power-heat coupling modeling method |
-
2014
- 2014-04-29 CN CN201410178375.7A patent/CN104166747B/en not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5377116A (en) * | 1991-07-01 | 1994-12-27 | Valenite Inc. | Method and system for designing a cutting tool |
| CN102063552A (en) * | 2011-01-07 | 2011-05-18 | 西安交通大学 | Method for determining rigidity of bolt connecting piece |
| CN102592017A (en) * | 2011-12-31 | 2012-07-18 | 北京工业大学 | Two-sided locking knife handle/main shaft coupling performance simulating and optimizing method |
| CN103186698A (en) * | 2013-04-12 | 2013-07-03 | 北京工业大学 | Simulation optimization method of dynamic and static performances of static pressure turntable of heavy machine tool |
| CN103366050A (en) * | 2013-06-25 | 2013-10-23 | 华中科技大学 | High-speed electric spindle power-heat coupling modeling method |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106529035A (en) * | 2016-11-09 | 2017-03-22 | 北京工业大学 | Method for analyzing dynamic characteristics of bolted structure based on three-dimensional fractal theory |
| CN106815407A (en) * | 2016-12-22 | 2017-06-09 | 四川大学 | A kind of whole machine dynamic property optimization method of Digit Control Machine Tool |
| CN107679301A (en) * | 2017-09-20 | 2018-02-09 | 西安建筑科技大学 | A kind of segmented heavy duty crossbeam scale model design method |
| CN107679301B (en) * | 2017-09-20 | 2020-11-24 | 西安建筑科技大学 | Sectional heavy-load cross beam scaling model design method |
| CN108563850A (en) * | 2018-03-29 | 2018-09-21 | 内蒙古久和能源装备有限公司 | A kind of Random Response Analysis method and system considering material damping otherness |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104166747B (en) | 2017-06-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| WO2021212680A1 (en) | Dynamic grid division-based method and system for analyzing load of non-standard part | |
| CN104166747A (en) | Bolting joint part dynamic characteristic analysis method taking surface machining quality into consideration | |
| JP2020071887A (en) | Designing mechanical part with topology optimization | |
| Chan | The overgrid interface for computational simulations on overset grids | |
| CN102592017B (en) | Two-sided locking knife handle/main shaft coupling performance simulating and optimizing method | |
| CN104077428A (en) | Remote finite element analysis method serving for industry alliance | |
| CN103488812A (en) | Satellite in-orbit micro-vibration simulation modeling correction method | |
| Palacios et al. | Adjoint-based method for supersonic aircraft design using equivalent area distribution | |
| CN101604353A (en) | The method for designing of airfoil assembly | |
| CN103324800B (en) | Boom-type roadheader cutting head computer-implemented method and CAD system | |
| CN110889166A (en) | Aviation bearing support lightweight design method based on thermal coupling constraint | |
| CN112380650A (en) | Method for designing structural member of working device | |
| CN110008523B (en) | Rock-soil excavation numerical simulation method for CEL method | |
| CN107563085A (en) | The finite element modeling method and strength assessment method of a kind of high-strength bolt attachment structure | |
| CN107784181A (en) | A kind of fluid structurecoupling towards carrier rocket emulates geometric model simplification method | |
| CN109726435B (en) | Method for calculating load of joint nail group | |
| CN103870634A (en) | Quick design optimization method for large structure | |
| CN109614748B (en) | Structure optimization method for improving dynamic characteristics of machine tool by combining test and simulation technology | |
| CN103294866A (en) | Tool-finite-element-analysis-oriented screw nut database creating method | |
| Anand et al. | Process automation of simulation using Toolkit/Tool Command Language (TK/TCL) scripting | |
| CN113609742A (en) | Wind generating set main shaft optimization method for overall structure | |
| CN103177163A (en) | Follow-up forecast method for milling deformation of complex integral structural components | |
| Ju et al. | Fracture analysis for damaged aircraft fuselage subjected to blast | |
| Arun | BUCKLING ANALYSIS AND MATERIAL SELECTION OF FUSELAGE FRAME USING CAD/CAM SOFTWARE | |
| CN118484981A (en) | Fiber reinforced composite material impact-resistant full-flow automatic calculation method and device |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170616 |