CN111460714B - Auxiliary stiffness approximation calculation method and system suitable for multi-body static contact problem - Google Patents
Auxiliary stiffness approximation calculation method and system suitable for multi-body static contact problem Download PDFInfo
- Publication number
- CN111460714B CN111460714B CN202010244971.6A CN202010244971A CN111460714B CN 111460714 B CN111460714 B CN 111460714B CN 202010244971 A CN202010244971 A CN 202010244971A CN 111460714 B CN111460714 B CN 111460714B
- Authority
- CN
- China
- Prior art keywords
- displacement
- stress
- stiffness
- rigidity
- auxiliary
- 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
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides an auxiliary stiffness approximation calculation method and system suitable for a multi-body static contact problem, which comprises the following steps: step 1: establishing a finite element model of each part; step 2: fixedly connecting each contact component with the ground space according to the auxiliary rigidity unit; and step 3: endowing the auxiliary stiffness unit with stiffness, and calculating displacement and stress; and 4, step 4: obtaining a function of displacement and rigidity and a function of stress and rigidity according to the displacement and the stress; and 5: and obtaining the displacement and the stress under the condition of no auxiliary rigidity according to the function of the displacement and the rigidity and the function of the stress and the rigidity. The invention can reflect the contact state between the parts faithfully and improve the confidence coefficient of the calculation result; the problem that parts are difficult to converge in a free contact state is effectively avoided, the convergence speed of contact calculation is improved, and the overall difficulty of the contact calculation is reduced; the method is suitable for multi-task synchronous calculation and improves simulation efficiency.
Description
Technical Field
The invention relates to the technical field of simulation calculation, in particular to an auxiliary stiffness approximation calculation method and system suitable for a multi-body static force contact problem.
Background
The most of parts in engineering are structurally formed by mechanical connection, and the force between different parts is transmitted by means of contact, extrusion and even impact among parts, which is called as a contact problem in mechanics. The contact problem is widely existed in various discipline fields such as civil engineering, machinery, aerospace and the like and practical engineering problems as a boundary nonlinear problem. In most cases, in order to design and assemble a machine, it is necessary to know the stress on the contact portions of the components, the size of the contact area, and other mechanical indexes. Therefore, the research on the contact problem has very important engineering practical significance. However, in the actual calculation process, due to the limitation of the calculation method, the multi-static contact calculation is often difficult to converge, and has a large calculation difficulty.
Therefore, a calculation method capable of faithfully reflecting the actual contact state of the component, increasing the convergence rate of the multi-body static force contact calculation and reducing the calculation difficulty is needed in the industry.
Patent document CN107330158a (application number: 201710426821.5) discloses a bearing radial stiffness identification method based on a local fine contact model, which includes the following steps: 1) Establishing a local fine contact model of the bearing; 2) Establishing an integral equivalent model of the bearing; 3) Respectively carrying out static calculation on the two models, and observing the radial displacement of the central points of the two models; 4) Taking the radial displacement of the central point of the local fine contact model as a standard, adjusting the rigidity value of the spring unit of the integral equivalent model to ensure that the radial displacement of the central point of the spring unit is the same as the radial displacement of the central point of the local fine contact model, and considering that the total rigidity value of the spring unit at the moment is the radial rigidity value of the whole bearing under the load Fr; 5) Changing the external load Fr to obtain recognition rigidity values under multiple groups of loads, fitting a curve of the radial rigidity of the bearing along with the change of the loads, and finally recognizing to obtain a bearing equivalent radial rigidity change diagram; 6) And a reasonable bearing modal test scheme is designed, and the accuracy of the identification rigidity is verified.
Disclosure of Invention
Aiming at the defects in the prior art, the invention aims to provide an auxiliary stiffness approximation calculation method and system suitable for a multi-body static contact problem.
The auxiliary stiffness approximation calculation method suitable for the multi-body static force contact problem provided by the invention comprises the following steps:
step 1: establishing a finite element model of each part;
step 2: fixedly connecting each contact component with the ground space according to the auxiliary rigidity unit;
and step 3: endowing the auxiliary stiffness unit with stiffness, and calculating displacement and stress;
and 4, step 4: obtaining a function of displacement and rigidity and a function of stress and rigidity according to the displacement and the stress;
and 5: and obtaining the displacement and the stress under the condition of no auxiliary rigidity according to the function of the displacement and the rigidity and the function of the stress and the rigidity.
Preferably, the step 2 includes: a virtual auxiliary stiffness unit is established between the contact member and the ground space.
Preferably, the step 3 comprises: giving the auxiliary stiffness unit stiffness K1, completing the 1 st round of calculation, and obtaining displacement W1 and stress S1;
imparting auxiliary stiffness Unit stiffness K 2 And completing the 2 nd round calculation to obtain the displacement W 2 And stress S 2 ；
Imparting stiffness K to an auxiliary stiffness unit n And completing the nth calculation to obtain the displacement W n And stress S n ；
Wherein, K 2 <<K 1 ，K n <<K n-1 。
Preferably, the step 4 comprises: the obtained data were fitted to obtain a displacement and stiffness function W = f (K) and a stress and stiffness function S = g (K).
Preferably, the step 5 comprises: bringing K =0 into W = f (K) and S = g (K), resulting in the unassisted stiffness conditionDisplacement W of lower part 0 = f (0) and stress S 0 ＝g(0)。
The invention provides an auxiliary stiffness approximation computing system suitable for a multi-body static force contact problem, which comprises:
a module M1: establishing a finite element model of each part;
a module M2: fixedly connecting each contact component with the ground space according to the auxiliary rigidity unit;
a module M3: endowing the auxiliary stiffness unit with stiffness, and calculating displacement and stress;
a module M4: obtaining a function of displacement and rigidity and a function of stress and rigidity according to the displacement and the stress;
a module M5: and obtaining the displacement and the stress under the condition of no auxiliary rigidity according to the function of the displacement and the rigidity and the function of the stress and the rigidity.
Preferably, the module M2 comprises: a virtual auxiliary stiffness unit is established between the contact member and the ground space.
Preferably, the module M3 comprises: giving the auxiliary stiffness unit stiffness K1, completing the 1 st round of calculation, and obtaining displacement W1 and stress S1;
imparting auxiliary stiffness Unit stiffness K 2 And completing the 2 nd round calculation to obtain the displacement W 2 And stress S 2 ；
Imparting stiffness K to an auxiliary stiffness unit n And completing the nth calculation to obtain the displacement W n And stress S n ；
Wherein, K 2 <<K 1 ，K n <<K n-1 。
Preferably, the module M4 comprises: the obtained data were fitted to obtain a displacement versus stiffness function W = f (K) and a stress versus stiffness function S = g (K).
Preferably, the module M5 comprises: bringing K =0 into W = f (K) and S = g (K) yields the displacement W without auxiliary stiffness 0 = f (0) and stress S 0 ＝g(0)。
Compared with the prior art, the invention has the following beneficial effects:
1. the invention can reflect the contact state between the parts faithfully and improve the confidence coefficient of the calculation result;
2. the method effectively avoids the problem that the parts are difficult to converge in a free contact state, improves the convergence speed of contact calculation, and reduces the overall difficulty of the contact calculation;
3. the method is suitable for multi-task synchronous calculation and improves the simulation efficiency.
Drawings
Other features, objects and advantages of the invention will become more apparent upon reading of the detailed description of non-limiting embodiments with reference to the following drawings:
FIG. 1 is a flow chart of an auxiliary stiffness approximation calculation method applicable to the multi-body static contact problem.
Detailed Description
The present invention will be described in detail with reference to specific examples. The following examples will assist those skilled in the art in further understanding the invention, but are not intended to limit the invention in any way. It should be noted that it would be obvious to those skilled in the art that various changes and modifications can be made without departing from the spirit of the invention. All falling within the scope of the present invention.
The invention provides an auxiliary stiffness approximation calculation method suitable for a multi-body static force contact problem, which comprises the following steps:
step one, completing the establishment of finite element models of all parts;
step two, fixedly connecting each contact part with the ground space through a plurality of auxiliary rigidity units;
step three, endowing rigidity K for the auxiliary rigidity unit 1 Completing the 1 st calculation to obtain the displacement W 1 And stress S 1 ；
Step four, endowing rigidity K for the auxiliary rigidity unit 2 (K 2 <<K 1 ) And completing the 2 nd round calculation to obtain the displacement W 2 And stress S 2 ；
Step five, endowing rigidity K for the auxiliary rigidity unit n (K n <<K n-1 )，Completing the nth calculation to obtain the displacement W n And stress S n ；
Step six, fitting the obtained data to obtain a function W = f (K) of displacement and rigidity and a function S = g (K) of stress and rigidity;
step seven, substituting K =0 into W = f (K) and S = g (K), and obtaining the displacement W under the condition of no auxiliary rigidity 0 = f (0) and stress S 0 ＝g(0)。
And step two, establishing a virtual auxiliary rigidity unit between the contact component and the ground space.
And in the second step to the fifth step, different rigidity characteristics are sequentially given to the auxiliary rigidity unit and calculation is carried out.
The data are fitted in the sixth step and the seventh step, and finally the displacement W without auxiliary rigidity is obtained 0 = f (0) and stress S 0 ＝g(0)。
The invention provides an auxiliary stiffness approximation computing system suitable for a multi-body static force contact problem, which comprises:
a module M1: establishing a finite element model of each part;
a module M2: fixedly connecting each contact component with the ground space according to the auxiliary rigidity unit;
a module M3: endowing the auxiliary stiffness unit with stiffness, and calculating displacement and stress;
a module M4: obtaining a function of displacement and rigidity and a function of stress and rigidity according to the displacement and the stress;
a module M5: and obtaining the displacement and the stress under the condition of no auxiliary rigidity according to the function of the displacement and the rigidity and the function of the stress and the rigidity.
Those skilled in the art will appreciate that, in addition to implementing the systems, apparatus, and various modules thereof provided by the present invention in purely computer readable program code, the same procedures can be implemented entirely by logically programming method steps such that the systems, apparatus, and various modules thereof are provided in the form of logic gates, switches, application specific integrated circuits, programmable logic controllers, embedded microcontrollers and the like. Therefore, the system, the device and the modules thereof provided by the present invention can be considered as a hardware component, and the modules included in the system, the device and the modules thereof for implementing various programs can also be considered as structures in the hardware component; modules for performing various functions may also be considered to be both software programs for performing the methods and structures within hardware components.
The foregoing description of specific embodiments of the present invention has been presented. It is to be understood that the present invention is not limited to the specific embodiments described above, and that various changes or modifications may be made by one skilled in the art within the scope of the appended claims without departing from the spirit of the invention. The embodiments and features of the embodiments of the present application may be combined with each other arbitrarily without conflict.
Claims (4)
1. An auxiliary stiffness approximation calculation method suitable for a multi-body static force contact problem is characterized by comprising the following steps:
step 1: establishing a finite element model of each part;
and 2, step: fixedly connecting each contact component with the ground space according to the auxiliary rigidity unit;
and step 3: endowing the auxiliary stiffness unit with stiffness, and calculating displacement and stress;
and 4, step 4: obtaining a function of displacement and rigidity and a function of stress and rigidity according to the displacement and the stress;
and 5: obtaining the displacement and the stress under the condition without auxiliary rigidity according to the function of the displacement and the rigidity and the function of the stress and the rigidity;
the step 3 comprises the following steps: giving the auxiliary stiffness unit stiffness K1, completing the 1 st round of calculation, and obtaining displacement W1 and stress S1;
imparting auxiliary stiffness Unit stiffness K 2 And completing the 2 nd round calculation to obtain the displacement W 2 And stress S 2 ；
Imparting stiffness K to an auxiliary stiffness unit n And completing the nth calculation to obtain the displacement W n And stress S n ；
Wherein, K 2 <<K 1 ，K n <<K n-1 ；
The step 4 comprises the following steps: fitting the obtained data to obtain a function W = f (K) of displacement and rigidity and a function S = g (K) of stress and rigidity;
the step 5 comprises the following steps: bringing K =0 into W = f (K) and S = g (K) yields the displacement W without auxiliary stiffness 0 = f (0) and stress S 0 ＝g(0)。
2. The method for calculating the aided stiffness approximation for the multi-body static contact problem according to claim 1, wherein the step 2 comprises: a virtual auxiliary stiffness unit is established between the contact member and the ground space.
3. An aided stiffness approximation computing system for use in a multi-volume static contact problem, comprising:
a module M1: establishing a finite element model of each part;
a module M2: fixedly connecting each contact component with the ground space according to the auxiliary rigidity unit;
a module M3: endowing the auxiliary stiffness unit with stiffness, and calculating displacement and stress;
a module M4: obtaining a function of displacement and rigidity and a function of stress and rigidity according to the displacement and the stress;
a module M5: obtaining the displacement and the stress under the condition without auxiliary rigidity according to the function of the displacement and the rigidity and the function of the stress and the rigidity;
the module M3 comprises: giving the auxiliary stiffness unit stiffness K1, completing the 1 st round of calculation, and obtaining displacement W1 and stress S1;
imparting auxiliary stiffness Unit stiffness K 2 And completing the 2 nd round calculation to obtain the displacement W 2 And stress S 2 ；
Imparting stiffness K to an auxiliary stiffness unit n And completing the nth calculation to obtain the displacement W n And stress S n ；
Wherein, K 2 <<K 1 ，K n <<K n-1 ；
The module M4 comprises: fitting the obtained data to obtain a displacement and rigidity function W = f (K) and a stress and rigidity function S = g (K);
the module M5 comprises: bringing K =0 into W = f (K) and S = g (K) yields the displacement W without auxiliary stiffness 0 = f (0) and stress S 0 ＝g(0)。
4. The system for calculating the aided stiffness approximation suitable for the multi-body static contact problem according to claim 3, wherein the module M2 comprises: a virtual auxiliary stiffness unit is established between the contact member and the ground space.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010244971.6A CN111460714B (en) | 2020-03-31 | 2020-03-31 | Auxiliary stiffness approximation calculation method and system suitable for multi-body static contact problem |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010244971.6A CN111460714B (en) | 2020-03-31 | 2020-03-31 | Auxiliary stiffness approximation calculation method and system suitable for multi-body static contact problem |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111460714A CN111460714A (en) | 2020-07-28 |
| CN111460714B true CN111460714B (en) | 2022-11-29 |
Family
ID=71679443
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010244971.6A Active CN111460714B (en) | 2020-03-31 | 2020-03-31 | Auxiliary stiffness approximation calculation method and system suitable for multi-body static contact problem |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111460714B (en) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106202802A (en) * | 2016-07-22 | 2016-12-07 | 中国航空工业集团公司西安飞机设计研究所 | A kind of undercarriage course stiffness simulation method |
| CN106650073A (en) * | 2016-12-14 | 2017-05-10 | 燕山大学 | Balanced force system boundary condition submodel analysis method |
| CN107330158A (en) * | 2017-06-08 | 2017-11-07 | 南京航空航天大学 | The equivalent radial rigidity recognition methods of bearing based on local fine contact model |
| CN109726436A (en) * | 2018-12-04 | 2019-05-07 | 中国航空工业集团公司西安飞机设计研究所 | A kind of wing-box rigidity determines method |
| CN110147628A (en) * | 2019-05-29 | 2019-08-20 | 台州学院 | Consider the tire cornering stiffness zoning method for calculating of multifactor variation |
| CN110162829A (en) * | 2019-07-15 | 2019-08-23 | 华南理工大学 | A kind of three ball pins-Rzeppa constant velocity drive shaft axially derives from the calculation method of power |
| CN110287593A (en) * | 2019-06-25 | 2019-09-27 | 东北大学 | One kind being bolted model interface parameter identification method |
| CN110399628A (en) * | 2019-05-17 | 2019-11-01 | 中国人民解放军军事科学院评估论证研究中心 | Consider that the level pressure that vibration influences pre-tightens angular contact bearing rigidity data calibration method |
| CN110516384A (en) * | 2019-08-30 | 2019-11-29 | 西北工业大学 | A kind of multiple-input and multiple-output gear train assembly dynamic modeling method |
| CN110555229A (en) * | 2019-07-12 | 2019-12-10 | 北京航空航天大学 | Meshless solid mechanics simulation method, electronic equipment and storage medium |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2029996A2 (en) * | 2006-06-02 | 2009-03-04 | Mts Systems Corporation | Measurement of properties of thin specimens based on experimentally acquired force-displacement data |
-
2020
- 2020-03-31 CN CN202010244971.6A patent/CN111460714B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106202802A (en) * | 2016-07-22 | 2016-12-07 | 中国航空工业集团公司西安飞机设计研究所 | A kind of undercarriage course stiffness simulation method |
| CN106650073A (en) * | 2016-12-14 | 2017-05-10 | 燕山大学 | Balanced force system boundary condition submodel analysis method |
| CN107330158A (en) * | 2017-06-08 | 2017-11-07 | 南京航空航天大学 | The equivalent radial rigidity recognition methods of bearing based on local fine contact model |
| CN109726436A (en) * | 2018-12-04 | 2019-05-07 | 中国航空工业集团公司西安飞机设计研究所 | A kind of wing-box rigidity determines method |
| CN110399628A (en) * | 2019-05-17 | 2019-11-01 | 中国人民解放军军事科学院评估论证研究中心 | Consider that the level pressure that vibration influences pre-tightens angular contact bearing rigidity data calibration method |
| CN110147628A (en) * | 2019-05-29 | 2019-08-20 | 台州学院 | Consider the tire cornering stiffness zoning method for calculating of multifactor variation |
| CN110287593A (en) * | 2019-06-25 | 2019-09-27 | 东北大学 | One kind being bolted model interface parameter identification method |
| CN110555229A (en) * | 2019-07-12 | 2019-12-10 | 北京航空航天大学 | Meshless solid mechanics simulation method, electronic equipment and storage medium |
| CN110162829A (en) * | 2019-07-15 | 2019-08-23 | 华南理工大学 | A kind of three ball pins-Rzeppa constant velocity drive shaft axially derives from the calculation method of power |
| CN110516384A (en) * | 2019-08-30 | 2019-11-29 | 西北工业大学 | A kind of multiple-input and multiple-output gear train assembly dynamic modeling method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111460714A (en) | 2020-07-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Menghal et al. | Real time simulation: Recent progress & challenges | |
| CN111324991B (en) | Reconstruction method of aerodynamic model in ground flutter test | |
| CN109376378B (en) | Front-back processing method for rapidly extracting rigidity result of white car body mounting point | |
| CN109902415B (en) | Method for calculating local stress strain of notch piece under high-temperature multiaxial load | |
| CN104809273A (en) | Creep deformation describing method | |
| CN106934127B (en) | Method for acquiring stress state of metal framework of power assembly suspension system | |
| CN103838913B (en) | The Finite Element of the curved bridge of curved box girder | |
| CN112414668B (en) | Wind tunnel test data static bomb correction method, device, equipment and medium | |
| Cella et al. | Geometric parameterization strategies for shape optimization using RBF mesh morphing | |
| CN112818462A (en) | Method and device for generating wheel parameter model, storage medium and computer equipment | |
| He et al. | A Coupled Newton--Krylov Time Spectral Solver for Flutter Prediction | |
| Peng et al. | Numerical investigation of the effects of structural geometric and material nonlinearities on limit-cycle oscillation of a cropped delta wing | |
| CN115422654A (en) | CFD/CSD technology-based efficient high-precision flutter time domain analysis method for cross/supersonic aircraft | |
| CN111460714B (en) | Auxiliary stiffness approximation calculation method and system suitable for multi-body static contact problem | |
| CN108595769B (en) | Shock absorber rigidity simulation method based on optimization algorithm | |
| CN100570619C (en) | A kind of towards products grade simulated model of parts order reducing method | |
| CN110989497B (en) | Multi-axis multi-excitation vibration control combination selection method and system based on iterative optimization | |
| CN110794813B (en) | In-loop real-time linkage testing method and system for electromechanical parts | |
| CN116151156A (en) | Aerodynamic force order reduction modeling method applicable to any aerodynamic parameters and modal shape | |
| CN116141314A (en) | Method and system for identifying dynamic parameters of robot based on projective geometry algebra | |
| CN116561885A (en) | Data-driven supercritical airfoil and mapping method and system of flow field of data-driven supercritical airfoil | |
| CN113065269A (en) | Composite material multi-scale concurrent simulation method and system | |
| Li et al. | Structural dynamic reanalysis method for transonic aeroelastic analysis with global structural modifications | |
| Yang et al. | Improved Fluid-Structure Interface for Aeroelastic Computations with Non-Matching Outer Mold Lines | |
| Mihaylova et al. | On the improvement of concept modeling of joints within simplified finite element models with application to structural dynamics |
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 |