CN106777691B - Rubber o-ring finite element modeling method for Structural Dynamics emulation - Google Patents
Rubber o-ring finite element modeling method for Structural Dynamics emulation Download PDFInfo
- Publication number
- CN106777691B CN106777691B CN201611163679.1A CN201611163679A CN106777691B CN 106777691 B CN106777691 B CN 106777691B CN 201611163679 A CN201611163679 A CN 201611163679A CN 106777691 B CN106777691 B CN 106777691B
- Authority
- CN
- China
- Prior art keywords
- ring
- finite element
- equivalent
- modal
- pore structure
- 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
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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Abstract
The invention discloses a kind of rubber o-ring finite element modeling methods for Structural Dynamics emulation, for solving the technical problem of existing rubber o-ring finite element modeling method complexity.Technical solution is to carry out structure to the system of axle construction, pore structure and O-ring composition to simplify, and establishes its finite element model；Then the modal test under free boundary condition is carried out to system, obtains the modal parameters such as intrinsic frequency, the vibration shape；With the minimum objective function of the error of Modal Analysis calculated result and two modal testing results, the elastic modulus E of O-ring is determined using optimization algorithm, to obtain the O-ring finite element model with degree of precision.The Dynamics Simulation Model of the available shaft hole matching structure with higher simulation accuracy of this method, and then high-precision dynamical simulation results are obtained, dynamical simulation results and test result error are maintained within 10%, and method is simple.
Description
Technical field
The present invention relates to a kind of rubber o-ring finite element modeling methods, in particular to a kind of to emulate for Structural Dynamics
Rubber o-ring finite element modeling method.
Background technique
O-ring is one of most common sealing means of shaft hole matching structure, since its quality and rigidity are much smaller than axle construction
And pore structure, therefore when carrying out shaft hole matching structure integral power simulation analysis, O-ring structure is typically ignored.But
When the rigidity of coupling stiffness and O-ring between axle construction and pore structure is in the same order of magnitude or axle construction and pore structure
Between coupling stiffness when being provided completely by O-ring, ignore O-ring structure, will lead to overall structure kinetic model occur compared with
Big error.
The large-scale generals Nonlinear Finite metaprogram such as ANSYS generallys use Mooney-Rivlin model analysis and calculates rubber
The mechanical property of glue
Document " hydraulic servo actuator O-ring seals experimental study and finite element analysis [D], Harbin Institute of Technology,
2011, P31-47 " disclose a kind of rubber o-ring finite element modeling method, and this method is establishing O shape rubber seal finite element
In model process, the strain energy density function of several parameters expressions is selected to describe hyperelastic body incompressible material characteristic.And
In order to determine Mooney-Rivlin model parameter, document carries out equal twin shafts by the glue dumbbell shaped test specimen to butyronitrile rubber material and draws
Experiment, uniaxial tensile test and shearing experiment are stretched, according to the stress-strain data that experiment measures, is linearly returned using MATLAB software
Return the calculating with ANSYS software, found out the Mooney-Rivlin constant of test specimen, to be built for O shape rubber seal finite element
Mould provides material characteristic parameter.Due to the non-linear of rubber material, Incoercibility and large deformation, test method is used in document
It determines that the process of rubber material mechanical constant is many and diverse, and tests the glue dumb-bell model test specimen used and O-ring structure not
Unanimously, mechanical property parameters also have a deviation, the preparation of testpieces also time and effort consuming.Moreover, this kind of method is tested with statics
Determining material characteristic parameter applies to the static seals such as research material failure, the deformation depressed and reliability and fine motion sealing more
Specificity analysis, for being modeled by the dynamics simulation of seal with O ring connecting shaft and the unit construction system in hole, this kind of method is true
Fixed mechanical property parameters may not be applicable in.
Summary of the invention
In order to overcome the shortcomings of that existing rubber o-ring finite element modeling method is complicated.The present invention provides a kind of for structure
The rubber o-ring finite element modeling method of dynamics simulation.The system that this method forms axle construction, pore structure and O-ring into
Row structure simplifies, and establishes its finite element model；Then the modal test under free boundary condition is carried out to system, obtains intrinsic frequency
The modal parameters such as rate, the vibration shape；With the minimum objective function of the error of Modal Analysis calculated result and two modal testing results, use is excellent
Change the elastic modulus E that algorithm determines O-ring, to obtain the O-ring finite element model with degree of precision.This method can obtain
The Dynamics Simulation Model of the shaft hole matching structure with higher simulation accuracy is taken, and then obtains high-precision dynamics simulation knot
Fruit, method are simple.
The technical solution adopted by the present invention to solve the technical problems is: a kind of RUBBER O for Structural Dynamics emulation
Shape circle finite element modeling method, its main feature is that the following steps are included:
Step 1: 2 structure of O-ring is equivalent.
Perfect condition condition for the ease of establishing the finite element model of O-ring 2, after the setting assembly of O-ring 2 are as follows:
A) after O-ring 2 is packed into the O-ring slot that pore structure 3 and axle construction 1 are constituted, O-ring slot is filled up；
B) 2 outer surface of O-ring and pore structure 3 are without relative motion；
C) 2 inner surface of O-ring and axle construction 1 are without relative motion.
Step 2: establishing the finite element model of system.
Axle construction 1 and pore structure 3 are modeled by practical structures size, elastic modulus E, Poisson's ratio ν and density p parameter
By the physical parameter setting of structural material；O-ring 2 is modeled according to perfect condition condition a), and Poisson's ratio ν and density p press O
The physical parameter setting of 2 material of shape circle, elastic modulus E are variable undetermined；Company between O-ring 2 and axle construction 1 and pore structure 3
Relationship is connect to be handled according to perfect condition condition b) and c).
Step 3: determining modal parameter.
The system formed using axle construction 1, pore structure 3 and O-ring 2 is object, under free boundary condition, carries out to system
Modal test obtains the intrinsic frequency of system, the vibration shape these modal parameters under free boundary condition.
Step 4: determining equivalent 4 equivalent elastic modulus E of O-ringx。
The equivalent elastic modulus E of equivalent O-ring 4 is solved using optimization algorithmx.Define the Equivalent Elasticity mould of equivalent O-ring 4
Measure ExAs design variable, ExValue range be 0.01E~10E, E be material elasticity modulus；It defines Modal Analysis and calculates knot
The minimum objective function f (x) of the error of fruit and two modal testing results, i.e.,
In formula, i indicates the order of intrinsic frequency and Mode Shape；Indicate the i-th rank intrinsic frequency simulation calculation value
Absolute value of the difference is missed with test value；MACiIndicate the correlation of the i-th rank Mode Shape simulation calculation value and test value；αiAnd βiPoint
The error of the i-th rank intrinsic frequency and the weighting coefficient of MAC value residual error are not indicated,And MACiIt is that equivalent O-ring 4 is equivalent
Elastic modulus ExFunction.
The αiAnd βiValue is 1.
The beneficial effects of the present invention are: this method carries out structure letter to the system that axle construction, pore structure and O-ring form
Change, establishes its finite element model；Then the modal test under free boundary condition is carried out to system, obtains intrinsic frequency, the vibration shape
Equal modal parameters；With the minimum objective function of the error of Modal Analysis calculated result and two modal testing results, using optimization algorithm
The elastic modulus E of O-ring is determined, to obtain the O-ring finite element model with degree of precision.This method is available to be had
The Dynamics Simulation Model of the shaft hole matching structure of higher simulation accuracy, and then high-precision dynamical simulation results are obtained, it moves
Mechanics Simulation result and test result error are maintained within 10%, and method is simple.
It elaborates with reference to the accompanying drawings and detailed description to the present invention.
Detailed description of the invention
Fig. 1 is the O-ring structure chart of shaft hole matching structure involved in the method for the present invention.
Fig. 2 be it is equivalent after O-ring structure chart.
In figure, 1- axle construction, 2-O shape circle, 3- pore structure, the equivalent O-ring of 4-.
Specific embodiment
Referring to Fig.1-2.The present invention is using the O-ring of shaft hole matching structure as research object, for Structural Dynamics emulation
Specific step is as follows for rubber o-ring finite element modeling method:
Step 1: 2 structure of O-ring is equivalent.
Perfect condition condition for the ease of establishing the finite element model of O-ring 2, after the setting assembly of O-ring 2 are as follows:
A) after O-ring 2 is packed into the O-ring slot that pore structure 3 and axle construction 1 are constituted, O-ring slot is filled up；
B) 2 outer surface of O-ring and pore structure 3 are without relative motion；
C) 2 inner surface of O-ring and axle construction 1 are without relative motion.
Step 2: establishing the finite element model of system.
Axle construction 1 and pore structure 3 are modeled by practical structures size, the ginseng such as elastic modulus E, Poisson's ratio ν and density p
Number presses the physical parameter setting of structural material；O-ring 2 is modeled according to perfect condition condition a), Poisson's ratio ν and density p
By the physical parameter setting of 2 material of O-ring, elastic modulus E is variable undetermined；Between O-ring 2 and axle construction 1 and pore structure 3
Connection relationship is handled according to perfect condition condition b) and c).
Step 3: determining modal parameter.
The system formed using axle construction 1, pore structure 3 and O-ring 2 is object, under free boundary condition, carries out to system
Modal test obtains the modal parameters such as intrinsic frequency, vibration shape of system under free boundary condition.
Step 4: determining equivalent 4 equivalent elastic modulus E of O-ringx。
The equivalent elastic modulus E of equivalent O-ring 4 is solved using optimization algorithmx.Define the Equivalent Elasticity mould of equivalent O-ring 4
Measure ExAs design variable, value range is defined as 0.01E~10E (elasticity modulus that E is material)；Define Modal Analysis meter
The minimum objective function f (x) of error of result and two modal testing results is calculated, i.e.,
In formula, i indicates the order of intrinsic frequency and Mode Shape；Indicate the i-th rank intrinsic frequency simulation calculation value
Absolute value of the difference is missed with test value；MACiIndicate that the i-th rank Mode Shape simulation calculation value and test value obtain correlation；αiAnd βiPoint
The error of the i-th rank intrinsic frequency and the weighting coefficient of MAC value residual error are not indicated, take 1 under normal conditions.And MACi
For equivalent 4 equivalent elastic modulus E of O-ringxFunction.
Claims (2)
1. a kind of rubber o-ring finite element modeling method for Structural Dynamics emulation, it is characterised in that the following steps are included:
Step 1: O-ring (2) structure is equivalent；
Perfect condition condition for the ease of establishing the finite element model of O-ring (2), after setting O-ring (2) assembly are as follows:
A) after O-ring (2) is packed into the O-ring slot that pore structure (3) and axle construction (1) are constituted, O-ring slot is filled up；
B) O-ring (2) outer surface and pore structure (3) are without relative motion；
C) O-ring (2) inner surface and axle construction (1) are without relative motion；
Step 2: establishing the finite element model of system；
Axle construction (1) and pore structure (3) are modeled by practical structures size, elastic modulus E, Poisson's ratio ν and density p parameter
By the physical parameter setting of structural material；O-ring (2) is modeled according to perfect condition condition a), Poisson's ratio ν and density p
By the physical parameter setting of O-ring (2) material, elastic modulus E is variable undetermined；O-ring (2) and axle construction (1) and pore structure
(3) connection relationship between is according to perfect condition condition b) and c) it is handled；
Step 3: determining modal parameter；
The system formed using axle construction (1), pore structure (3) and O-ring (2) is object, under free boundary condition, to system into
Row modal test obtains the intrinsic frequency of system, the vibration shape these modal parameters under free boundary condition；
Step 4: determining equivalent O-ring (4) equivalent elastic modulus Ex；
The equivalent elastic modulus E of equivalent O-ring (4) is solved using optimization algorithmx；Define the Equivalent Elasticity mould of equivalent O-ring (4)
Measure ExAs design variable, ExValue range be 0.01E~10E, E be material elasticity modulus；It defines Modal Analysis and calculates knot
The minimum objective function f (x) of the error of fruit and two modal testing results, i.e.,
In formula, i indicates the order of intrinsic frequency and Mode Shape；Indicate the i-th rank intrinsic frequency simulation calculation value and test
Value misses absolute value of the difference；MACiIndicate the correlation of the i-th rank Mode Shape simulation calculation value and test value；αiAnd βiIt respectively indicates
The error of i-th rank intrinsic frequency and the weighting coefficient of MAC value residual error,And MACiIt is equivalent O-ring (4) Equivalent Elasticity
Modulus ExFunction.
2. the rubber o-ring finite element modeling method according to claim 1 for Structural Dynamics emulation, feature exist
In: the αiAnd βiValue is 1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611163679.1A CN106777691B (en) | 2016-12-16 | 2016-12-16 | Rubber o-ring finite element modeling method for Structural Dynamics emulation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611163679.1A CN106777691B (en) | 2016-12-16 | 2016-12-16 | Rubber o-ring finite element modeling method for Structural Dynamics emulation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106777691A CN106777691A (en) | 2017-05-31 |
| CN106777691B true CN106777691B (en) | 2019-12-03 |
Family
ID=58891614
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611163679.1A Active CN106777691B (en) | 2016-12-16 | 2016-12-16 | Rubber o-ring finite element modeling method for Structural Dynamics emulation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106777691B (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107292056A (en) * | 2017-07-21 | 2017-10-24 | 山东省科学院海洋仪器仪表研究所 | A kind of O-ring seals fast design method based on Finite Element Simulation Analysis |
| CN109086498B (en) * | 2018-07-17 | 2020-05-08 | 东风商用车有限公司 | Analysis method for fender bracket with rubber pad |
| CN110020485B (en) * | 2019-04-10 | 2022-12-13 | 东北大学 | Method for analyzing inherent characteristics of suspended thin-wall column shell based on bolt connection |
| CN111159934A (en) * | 2019-12-13 | 2020-05-15 | 北京电子工程总体研究所 | Aluminum matrix composite material dynamics simulation method |
| CN113779832B (en) * | 2021-09-06 | 2023-08-29 | 北京强度环境研究所 | High-precision finite element simulation model correction method |
| CN116011124B (en) * | 2023-03-20 | 2023-07-14 | 西安航天动力研究所 | Simulation modeling method of vibration test system |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102063549A (en) * | 2011-01-07 | 2011-05-18 | 西安交通大学 | Method for stimulating fluid-solid-heat coupling value of machine tool spindle |
| CN102661001A (en) * | 2012-05-28 | 2012-09-12 | 中建三局第三建设工程有限责任公司 | Connecting device for membranes penetrating through cylinder and construction method for connecting device |
| CN103632000A (en) * | 2013-11-25 | 2014-03-12 | 中国船舶重工集团公司第七〇五研究所 | Finite element modeling method of underwater vehicle wedge ring structure |
| CN105069241A (en) * | 2015-08-19 | 2015-11-18 | 山东大学 | Step-by-step analysis and prediction method for dynamic performances of rubber material structure |
| CN105740541A (en) * | 2016-01-29 | 2016-07-06 | 厦门大学 | Structural dynamical model modification-based prestress recognition method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8532941B2 (en) * | 2010-08-26 | 2013-09-10 | Siemens Aktiengesellschaft | Fatigue life estimation method and system |
-
2016
- 2016-12-16 CN CN201611163679.1A patent/CN106777691B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102063549A (en) * | 2011-01-07 | 2011-05-18 | 西安交通大学 | Method for stimulating fluid-solid-heat coupling value of machine tool spindle |
| CN102661001A (en) * | 2012-05-28 | 2012-09-12 | 中建三局第三建设工程有限责任公司 | Connecting device for membranes penetrating through cylinder and construction method for connecting device |
| CN103632000A (en) * | 2013-11-25 | 2014-03-12 | 中国船舶重工集团公司第七〇五研究所 | Finite element modeling method of underwater vehicle wedge ring structure |
| CN105069241A (en) * | 2015-08-19 | 2015-11-18 | 山东大学 | Step-by-step analysis and prediction method for dynamic performances of rubber material structure |
| CN105740541A (en) * | 2016-01-29 | 2016-07-06 | 厦门大学 | Structural dynamical model modification-based prestress recognition method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106777691A (en) | 2017-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106777691B (en) | Rubber o-ring finite element modeling method for Structural Dynamics emulation | |
| Silva et al. | Modeling bamboo as a functionally graded material: lessons for the analysis of affordable materials | |
| CN109030202B (en) | Method for rapidly determining discrete element model parameters of rock brittle material | |
| CN108959686A (en) | A kind of correction method for finite element model based on sensitivity analysis | |
| Han et al. | Ultimate bearing capacity of the welded hollow spherical joints in spatial reticulated structures | |
| CN106644711B (en) | A kind of ductile material single shaft constitutive relation test method | |
| CN103838913B (en) | The Finite Element of the curved bridge of curved box girder | |
| CN109115879B (en) | Structural damage identification method based on modal shape and convolutional neural network | |
| Şerban et al. | Behaviour of semi-crystalline thermoplastic polymers: Experimental studies and simulations | |
| CN106326501B (en) | The calculation method of the fabric structure kinematic analysis natural frequency of vibration and the vibration shape | |
| CN106840721B (en) | A kind of Flight Vehicle Structure posting characteristic model test design method | |
| CN110220781A (en) | A kind of plate anisotropy constitutive parameter scaling method and system | |
| CN106815442A (en) | A kind of incompressible hyperelastic body constitutive model of isotropism and its application | |
| CN106055733B (en) | The kinetic parameter of multifunction structure determines method | |
| CN107314938B (en) | The implementation method of nugget region material plastic inverting identification | |
| CN106813816A (en) | Counterweight balance is measured | |
| CN102297745A (en) | GM (0,2) model-based static decoupling method of multidimensional force sensor | |
| CN104008234B (en) | Method for correcting closely spaced mode model with damping structure | |
| Aboudi et al. | High-fidelity micromechanical modeling of continuously reinforced elastic multiphase materials undergoing finite deformations | |
| CN103542894B (en) | Reinforcement stresses, strain measurement method under high temperature change loading speed | |
| CN104765928B (en) | A kind of Plastic Forming frictional behavior measuring method | |
| CN105138746B (en) | Pass through the anti-method for pushing away its Young's modulus of the stiffness factor of loudspeaker vibration component | |
| El Gendy et al. | Finite element formulation of Timoshenko tapered beam-column element for large displacement analysis based on the exact shape functions | |
| CN104077444A (en) | Analysis method of indentation data | |
| CN103412968A (en) | Method for establishing constitutive relation model of hardening delay material |
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 |