CN104166747B - A kind of bolted joint Dynamic Characteristics Analysis Method for considering suface processing quality - Google Patents
A kind of bolted joint Dynamic Characteristics Analysis Method for considering suface processing quality Download PDFInfo
- Publication number
- CN104166747B CN104166747B CN201410178375.7A CN201410178375A CN104166747B CN 104166747 B CN104166747 B CN 104166747B CN 201410178375 A CN201410178375 A CN 201410178375A CN 104166747 B CN104166747 B CN 104166747B
- Authority
- CN
- China
- Prior art keywords
- analysis
- contact
- model
- faying face
- bolted joint
- 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.)
- Expired - Fee Related
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 53
- 238000012545 processing Methods 0.000 title claims abstract description 30
- 238000013016 damping Methods 0.000 claims abstract description 32
- 238000000034 method Methods 0.000 claims abstract description 27
- 230000003068 static effect Effects 0.000 claims abstract description 14
- 239000000463 material Substances 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims abstract description 5
- 238000005316 response function Methods 0.000 claims abstract description 5
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000013461 design Methods 0.000 claims description 3
- 238000012805 post-processing Methods 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 3
- 238000012916 structural analysis Methods 0.000 claims description 3
- 238000012217 deletion Methods 0.000 claims description 2
- 230000037430 deletion Effects 0.000 claims description 2
- 238000010408 sweeping Methods 0.000 claims description 2
- 230000001050 lubricating effect Effects 0.000 description 6
- 238000012360 testing method Methods 0.000 description 6
- 238000011161 development Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000003672 processing method Methods 0.000 description 4
- 230000004087 circulation Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000003801 milling Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 208000032443 Masked facies Diseases 0.000 description 1
- 238000005299 abrasion Methods 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
- 230000000694 effects Effects 0.000 description 1
- 238000000227 grinding Methods 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003746 surface roughness Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
A kind of bolted joint Dynamic Characteristics Analysis Method for considering suface processing quality, category Mechanical Property Analysis field.It includes：Parametric modeling, faying face are commonly modeled, apply boundary condition, static analysis, contact stiffness Damping calculating, faying face Model Reconstruction, the modeling of bolt faying face stiffness and damping, dynamic analysis step；Model is set up according to model parameter automatically, after selection static analysis, by the material parameter, contact method for processing surface, roughness and the static analysis result that select, calculations incorporated face ct clamping value, osculating element is re-created, dynamically analysis obtains joint portion intrinsic frequency, main shake shape and frequency response function to carry out bolted joint.The present invention is based on ANSYS secondary developing platforms, sets up bolted joint parameterized model, the influence for analysis of material attribute, contact method for processing surface, processing roughness, external load to bolted joint dynamic characteristic.
Description
Technical field
It is more particularly to a kind of to consider that contact method for processing surface and processing are coarse the invention belongs to software application technology field
The bolted joint Dynamic Characteristics Analysis Method of degree.
Background technology
There is substantial amounts of bolted joint in lathe, during lathe work, the dynamic characteristic influence of bolted joint
The processing efficiency of lathe, machining accuracy, the abrasion between acceleration components, and then influence the service life and lathe reliability of lathe
Property.Therefore the research of bolted joint is very meaningful to improving machine mass.
The dynamic characteristic of bolted joint is influenceed very big by material properties, surface quality, external load.By existing
Finite element analysis software can the influence of analysis of material attribute and external load to bolted joint dynamic characteristic, but but can not
Influence of the analysis contact surface surface quality to Structure dynamic characteristics.
Current domestic and foreign scholars are established by Fractal Contact theory and consider the contact stiffness of contact surface surface quality and connect
Touch damper model.Contact stiffness and contact damping model of the present invention based on Fractal Contact theory, using FEM-software ANSYS
Secondary developing platform, foundation considers the bolted joint parameterized model of contact method for processing surface and processing roughness, for dividing
Analyse its dynamic characteristic.
The content of the invention
It is an object of the present invention to overcome the defect of existing faying face emulation mode, there is provided a kind of based on Fractal Contact reason
The contact stiffness and contact damping model of opinion, using FEM-software ANSYS secondary developing platform, set up and consider contact surface processing
The bolted joint parameterized model of method and processing roughness, for analyzing its dynamic characteristic
The present invention realizes that realizing the system of the method includes using following technological means
A kind of bolt faying face dynamic emulation method for considering contact method for processing surface and roughness, including：Parametrization is built
Mould, faying face commonly modeled, apply boundary condition, static analysis, contact stiffness and contact damping are calculated, with reference to surface model weight
Build, bolt the modeling of faying face stiffness and damping, dynamic analysis step；
Parametric modeling step, cube and cylindrical structure are generated using body order；Grasped using Boolean subtraction calculation and bonding
Make generation bolted joint model；After generation physical model, using wire grid control command and grid sweeping order, generate limited
First grid model；In modeling process, using APDL language definition model parameterization variables, generation parameter model；
The common modeling procedure of faying face, faying face is selected by Parametric Design Language；Target face unit be set as
TARGE170 units, Contact surface element is set to be CONTA174 units；Soft-soft face contact is set, it is binding to initially contact；
Apply boundary condition step, including bolted joint constrained, load apply position, magnitude of load, load side
To；One section is set to fixed constraint, another section of applying pressure load at bolt hole；Load applying zone picks up section using coordinate method
Point；
Static analysis step, it is static structural analysis that analysis type is set first, sets and solves content, is solved；It is logical
Selection osculating element is crossed, all nodes of select unit, by loop command, extract each with array formats according to order from small to large
From node serial number and corresponding pressure values.
Contact stiffness and contact damping calculation procedure, it is firm by Fractal Contact according to the faying face node pressure values extracted
Degree and contact damping model, using loop command and formula operation, calculate the contact stiffness and contact damping value at each node, with
Array form is preserved.
Faying face Model Reconstruction, the contact mould set up in the common modeling procedure of faying face is deleted by element deletion order
Type, self-defined matrix unit is created by loop command, and setup unit is numbered, and contact stiffness and contact damping value are assigned to unit
Real constant, generates contact stiffness and contact damping unit, and joint surface contact stiffness is sequentially set up from small to large according to node serial number
With and damping unit.
Dynamic analysis step, sets analysis type, including model analysis, harmonic responding analysis；Model analysis is post-processed
Middle output intrinsic frequency and main shape of shaking；Key point frequency response function figure is exported in harmonic responding analysis post processing.
According to the above results, user can analyze contact surface machine-tooled method and processing roughness to bolted joint dynamic
The influence of characteristic.
Foregoing utilization APDL language definition structural parameters variables include：Bolted joint structural parameters, material number, spiral shell
Bolt pretightning force, with reference to method for processing surface, processing roughness；Wherein bolted joint structural parameters include total length, land length
Width, total height, boss height, diameter of bolt hole；Material parameter includes QT-235, QT250；Include milling with reference to method for processing surface
Cut, be ground；Processing roughness parameter includes five roughness values (0.4,0.6,0.8,1.2,1.6), and unit is micron.
Foregoing self-defined matrix unit is MATRIX27 matrix units.
A kind of side of the analysis bolted joint dynamic characteristic for considering contact method for processing surface and processing roughness of the present invention
Method, compared with prior art, with following obvious advantage and beneficial effect：
Contact stiffness and contact damping model of the present invention based on Fractal Contact theory, using FEM-software ANSYS two
Secondary development platform, sets up the bolted joint parameterized model for considering contact method for processing surface and processing roughness, for analyzing
Material properties, contact method for processing surface, surface roughness, external load, with reference to shape to the shadow of bolted joint dynamic characteristic
Ring.
Brief description of the drawings
Fig. 1 is bolted joint structural representation；
Fig. 2 is test specimen physical dimension schematic diagram one；
Fig. 3 is test specimen physical dimension schematic diagram two；
Fig. 4 is bolted joint dynamic analysis flow chart；
Fig. 5 is MATRIX27 geometric representations；
Fig. 6 is MATRIX27 matrix unit coefficient schematic diagrames.
In figure, 1, upper test specimen, 2, joint portion osculating element, 3, lower test specimen.
Specific embodiment
Below in conjunction with Figure of description, specific implementation of the invention is illustrated.
Refer to shown in Fig. 1, be bolted joint structural representation.Wherein, 1 is upper test specimen, and 2 is that joint portion contacts list
Unit, 3 is lower test specimen.
With reference to Fig. 2, Fig. 3, the present invention is furture elucidated for Fig. 4.
Wherein, L1 is total length, and L2 is land length, and D is width, and H1 is total height, and H2 is boss height, and d is bolt
Bore dia.
(1) parametric modeling module, based on ANSYS secondary development, cube and cylindrical structure is generated using body order；Adopt
Bolted joint model is generated with Boolean subtraction calculation and bond operation；|input paramete includes bolted joint structural parameters (overall length
Degree, land length width, total height, boss height, diameter of bolt hole), load parameter (imposed load area, magnitude of load)；
Selection parameter include material (QT-235, QT250), processing method (milling, grinding), process roughness value (0.4,0.6,0.8,
1.2、1.6)；After being provided with, FEM mesh is automatically generated.
(2) the common MBM of faying face, based on ANSYS secondary development, sets up target face unit for TARGE170 units,
Contact surface element is CONTA174 units, sets soft-soft face contact, and it is binding to initially contact.
(3) boundary condition module is applied, based on ANSYS secondary development, automatically according to input load parameter, using coordinate bit
Put pickup order and determine boundary condition applying zone, load applying zone and load value.
(4) static analysis module, selects static structural analysis, and program carries out static analysis calculating automatically, is ordered using circulation
Faying face node and corresponding pressure are extracted in order, and are shown with text and graphic form.
(5) contact stiffness Damping calculating module, program combines selection automatically by the combination node and correspondence pressure for extracting
Material, processing method, roughness calculate the corresponding ct clamping value of each node automatically, with array form preserve.
(6) faying face Model Reconstruction module (referring to shown in Fig. 1), deletes original faying face contact single using APDL language
Unit, according to the ct clamping value for calculating, creates the corresponding faying face faying face of each node and connects by MATRIX27 units
Touch rigidity and contact damping unit.
(7) dynamic analysis module, sets analysis type, including model analysis and harmonic responding analysis, after model analysis
Output modalities frequency and Mode Shape figure in treatment；Also need to further define method for solving and load step choosing in harmonic responding analysis
, key point frequency response function figure is exported in post processing.
User can be by changing size, magnitude of load, and selection material, contact surface connect processing method and roughness,
The influence of analysis of material attribute, contact surface crudy, external load to bolted joint dynamic characteristic.
Described parametric modeling, system changes existing Parametric Design Language according to the structural parameters of user input
Command stream, generates bolted joint model, and mesh generation is then carried out automatically.
Described lubricating pad faying face modeling, system generates each lubricating pad according to the lubricating pad parameter of input using Do-Loop circulations
Itself geometric center node and lubricating pad plane, and circulate generation supporting lubricating pad faying face and pretension lubricating pad faying face.
Described boundary condition, system changes APDL command streams according to the boundary condition parameter of input, and border is loaded automatically
Condition.
Described ct clamping is calculated, and system is according to material, processing method, roughness and the static analysis for selecting
The node pressure of acquisition calculates the corresponding contact stiffness of each node and contact damping automatically.
Described dynamic analysis, including model analysis and harmonic responding analysis, obtain bolted joint intrinsic frequency, master
The vibration shape and key point frequency response function.
In system faying face Model Reconstruction module faying face rigidity and resistance are built using self-defined matrix unit matrix27
Buddhist nun's unit.Matrix27 units can represent a kind of any unit in ANSYS programs.The unit is by two nodes and coefficient square
Battle array definition.Fig. 5 is MATRIX27 geometric representations.Wherein each node has 6 frees degree, be respectively along nodal coordinate system x, y,
The translation in z directions and the rotation around nodal coordinate system x, y, z direction.The geometrical property of the unit is without definition, but its elastic movement
Learning response can be specified with rigidity, damping or quality coefficient, i.e., the constant in rigidity, damping and mass matrix can be as list
The real constant input of unit.Fig. 6 is MATRIX27 matrix unit coefficient schematic diagrames.In analog equivalent rigidity unit, C1、C58、C7
Simulate the rigidity of X-direction, C13、C64、C19Simulate the rigidity of Y-direction, C24、C69、C30Simulate the rigidity of Z-direction.In simulation etc.
During effect damping unit, C1、C58、C7Simulate the damping of X-direction, C13、C64、C19Simulate the damping of Y-direction, C24、C69、C30Simulation Z
The damping in direction.
Claims (2)
1. it is a kind of consider suface processing quality bolted joint Dynamic Characteristics Analysis Method, including：Parametric modeling, faying face
Common modeling, applying boundary condition, static analysis, contact stiffness Damping calculating, faying face Model Reconstruction, faying face stiffness and damping
Modeling, dynamic analysis step；It is characterized in that：
Described parametric modeling step, cube and cylindrical structure are generated using body order；Using Boolean subtraction calculation and bonding
Operation generation bolted joint model；After generation physical model, using wire grid control command and grid sweeping order, generation has
The first grid model of limit；In modeling process, using APDL language definition model parameterization variables, generation parameter model；
The common modeling procedure of faying face, bolt faying face initial bind model is set up by Parametric Design Language；Target face list
Unit is TARGE170 units, and Contact surface element carries out soft-soft face contact for CONTA174 units；
Described applying boundary condition step, including bolted joint constrained, load apply position, loading, load
Size, loading direction；Wherein load applies position and magnitude of load is defined as input variable；
Described static analysis step, it is static structural analysis that analysis type is set first, sets and solves content, is solved；
, by loop command, extracted according to order from small to large with array formats by selecting osculating element, all nodes of select unit
Respective node serial number and corresponding pressure values；
Described contact stiffness Damping calculating step, according to the faying face node pressure values extracted, by Fractal Contact rigidity and
Contact damping model asks for faying face node contact rigidity and contact damping value, is preserved with array form；
Faying face Model Reconstruction, the contact model set up in the common modeling procedure of faying face is deleted by element deletion order, is led to
Cross loop command and create self-defined matrix unit, setup unit numbering, the node that contact stiffness Damping calculating step is obtained is firm
Degree and damping value are assigned to unit real constant, generate contact stiffness and contact damping unit, suitable from small to large according to node serial number
Sequence sets up joint surface contact stiffness and damping unit；
The dynamic analysis step, sets analysis type, including model analysis, harmonic responding analysis；Model analysis is post-processed
Middle output intrinsic frequency and main shape of shaking；Key point frequency response function figure is exported in harmonic responding analysis post processing；According to key point frequency response
Functional arrangement, user can analyze the influence of contact method for processing surface and processing roughness to bolted joint dynamic characteristic；
Bolted joint can be carried out by input load value, selection material, selection contact method for processing surface and processing roughness to move
Step response is analyzed.
2. it is according to claim 1 it is a kind of consider suface processing quality bolted joint Dynamic Characteristics Analysis Method, its
It is characterised by：Carry out static analysis and extract each node pressure, each node contact rigidity and damping are calculated by pressure, pass through
MATRIX27 sets 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 CN104166747A (en) | 2014-11-26 |
| CN104166747B true 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 (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108628250A (en) * | 2018-06-01 | 2018-10-09 | 西安交通大学 | A kind of optimization method of machine tool chief axis and handle of a knife fit structure size |
Families Citing this family (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 |
| 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 |
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 (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108628250A (en) * | 2018-06-01 | 2018-10-09 | 西安交通大学 | A kind of optimization method of machine tool chief axis and handle of a knife fit structure size |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104166747A (en) | 2014-11-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104166747B (en) | A kind of bolted joint Dynamic Characteristics Analysis Method for considering suface processing quality | |
| JP2020071887A (en) | Designing mechanical part with topology optimization | |
| CN109359353B (en) | Machine tool upright post structure optimization design method based on topology optimization | |
| Löhner | Recent advances in parallel advancing front grid generation | |
| CN104239624A (en) | Optimal design method for internal structure of machine tool body | |
| CN106650045A (en) | Optimized design and strength analysis method of secondary box section structure | |
| CN106777820A (en) | A kind of big part Dynamic performance Optimization method for designing of machine tool basic | |
| CN104462725B (en) | A kind of light weight optimal design method of veneer punching press lower control arm | |
| CN108984829A (en) | The calculation method and system of rock-fill concrete rockfill banking process | |
| CN112380650A (en) | Method for designing structural member of working device | |
| CN104484511A (en) | Simulation analysis based dynamic characteristic design method for robot structures | |
| CN110008523B (en) | Rock-soil excavation numerical simulation method for CEL method | |
| CN105468824B (en) | Rise the parameter simulation method of device assembly based on the moon | |
| Annarumma et al. | Methodology development of human task simulation as PLM solution related to OCRA ergonomic analysis | |
| CN103294866A (en) | Tool-finite-element-analysis-oriented screw nut database creating method | |
| CN101393582A (en) | Non-support grinding and high stiffness principal axis system simulation and analysis method | |
| CN110334459B (en) | Rapid and refined modeling system and method for power transmission tower line system | |
| Anand et al. | Process automation of simulation using Toolkit/Tool Command Language (TK/TCL) scripting | |
| Mueller et al. | Automated coupling of MDI/ADAMS and MSC. CONSTRUCT for the topology and shape optimization of flexible mechanical systems | |
| Baco et al. | Design of experiments used in computer trials: a supportive method for product development | |
| CN113609742A (en) | Wind generating set main shaft optimization method for overall structure | |
| US20230401351A1 (en) | Method and Apparatus for Obtaining a Composite Laminate | |
| Zhang et al. | A practical design process to optimize fatigue performance for chassis components | |
| Ding et al. | Cutting tool CAD/CAE integrating system modeling | |
| Li et al. | Assistant Teaching Platform for Creative Mechanical Design Course |
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 |
Granted publication date: 20170616 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |