CN113987868A - Simulation analysis method for multilayer thin-wall rubber-metal composite revolving body structure - Google Patents
Simulation analysis method for multilayer thin-wall rubber-metal composite revolving body structure Download PDFInfo
- Publication number
- CN113987868A CN113987868A CN202111207879.3A CN202111207879A CN113987868A CN 113987868 A CN113987868 A CN 113987868A CN 202111207879 A CN202111207879 A CN 202111207879A CN 113987868 A CN113987868 A CN 113987868A
- Authority
- CN
- China
- Prior art keywords
- simulation analysis
- dimensional
- axial
- load
- torsional
- 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.)
- Pending
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 169
- 238000004088 simulation Methods 0.000 title claims abstract description 154
- 239000002905 metal composite material Substances 0.000 title claims abstract description 60
- 229920001971 elastomer Polymers 0.000 claims abstract description 46
- 238000013507 mapping Methods 0.000 claims abstract description 4
- 238000004364 calculation method Methods 0.000 claims description 31
- 238000000034 method Methods 0.000 claims description 23
- 230000036316 preload Effects 0.000 claims description 17
- 239000000463 material Substances 0.000 claims description 12
- 238000007781 pre-processing Methods 0.000 claims description 3
- 238000012916 structural analysis Methods 0.000 claims description 2
- 230000015572 biosynthetic process Effects 0.000 claims 1
- 230000001360 synchronised effect Effects 0.000 claims 1
- 239000002131 composite material Substances 0.000 abstract description 2
- 230000035882 stress Effects 0.000 description 21
- 238000005516 engineering process Methods 0.000 description 8
- 230000007547 defect Effects 0.000 description 6
- 239000002184 metal Substances 0.000 description 6
- 238000013461 design Methods 0.000 description 5
- 239000000806 elastomer Substances 0.000 description 5
- 238000012545 processing Methods 0.000 description 4
- 239000007787 solid Substances 0.000 description 3
- 230000006872 improvement Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012827 research and development Methods 0.000 description 2
- 238000004904 shortening Methods 0.000 description 2
- 125000006850 spacer group Chemical group 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 230000032683 aging Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013329 compounding Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005538 encapsulation Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000003348 petrochemical agent Substances 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
A simulation analysis method of a multilayer thin-wall rubber metal composite revolving body structure comprises the following steps of firstly carrying out two-dimensional mesh subdivision operation on the multilayer thin-wall rubber metal composite revolving body; after the mesh subdivision of the two-dimensional axisymmetric cross section is completed, carrying out axial or torsional load simulation analysis; after the simulation analysis of the axial load or the torsional load is finished, converting the two-dimensional axisymmetric grid into a three-dimensional grid, and synchronously mapping a stress-strain result of the two-dimensional axisymmetric model on the newly generated three-dimensional grid to be used as initial stress strain of the simulation analysis of the deflection or radial load working condition; and then, the three-dimensional model simulation analysis under the working condition of deflection or radial load after axial or torsional preloading is carried out continuously, and the simulation analysis under the working condition of multiple loads of axial, radial, torsional and deflection is completed. The invention intercepts a part of section of the composite revolving body, and carries out simulation analysis of a three-dimensional model after mesh subdivision and load analysis, and realizes simulation analysis of axial and torsional load working conditions by using a simple two-dimensional axisymmetric model.
Description
Technical Field
The invention relates to a simulation analysis method of a stressed part, in particular to a simulation analysis method of a multilayer thin-wall rubber-metal composite revolving body structure, which greatly improves the simulation analysis efficiency of the multilayer thin-wall rubber-metal composite revolving body structure under the multi-load working conditions of radial, axial, deflection, torsion and the like by performing simulation analysis in a 180-degree or 360-degree rotation mode.
Background
On one hand, compared with the rubber metal composite revolving body with the traditional structure, the multilayer thin-wall rubber metal composite revolving body has higher bearing capacity, aging resistance and stability; tests also show that in the elastic metal structure of the multilayer thin-wall rubber-metal composite revolving body, the unit consumption of the elastomer can be reduced by one or more than ten orders of magnitude, and the use coefficient of the elastomer in the multilayer thin-wall rubber-metal composite revolving body structure is increased by more than nine times, so that the multilayer thin-wall rubber-metal composite revolving body is more and more widely applied at present. However, the multilayer thin-wall rubber-metal composite revolving body is different from the traditional rubber-metal composite revolving body in deformation characteristics, structural characteristics, and processing and using characteristics, so that the multilayer thin-wall rubber-metal composite revolving body is much more difficult to perform stress analysis.
On the other hand, with the increasing market competition and the shortening of the product update period, the enterprise demands new technology more urgently, the finite element simulation analysis technology is an effective means for improving the product quality, shortening the design period and improving the product competitiveness, therefore, with the development of computer technology and calculation method, finite element simulation analysis gets more and more extensive attention and application in the engineering design and scientific research field, has become an effective way to solve the problem of complex engineering analysis and calculation, almost all designs and manufactures from automobiles to space shuttles can not leave the finite element simulation analysis and calculation, its widespread use in various fields such as machine building, material processing, aerospace, automobiles, civil construction, electronic and electrical appliances, defense and military industry, ships, railways, petrochemicals, energy and scientific research has made a qualitative leap in the level of design. Therefore, a simulation analysis method is generally adopted when the stress analysis of the multilayer thin-wall rubber-metal composite revolving body is carried out, and the simulation analysis can help designers to optimize the optimal structure in research and development by means of a virtual prototype through computer and finite element simulation analysis, so that some schemes without advantages are eliminated, the test times are reduced, and the research and development cost and the sample piece cost are saved.
However, the existing general simulation analysis is to integrally grid the parts, and according to the conventional grid division method, each layer of rubber and each layer of spacer are divided into two-dimensional and three-dimensional grids separately, and the method has low grid division efficiency for the rod end bearing with the multilayer structure. Especially, the simulation analysis efficiency under the multi-load working conditions of radial, axial, deflection, torsion and the like is low, the three-dimensional mesh generation technology is difficult, and the defects of calculation non-convergence are easily caused, so that the improvement is needed.
Patent documents in which the same technology as that of the present invention is not found through patent search are reported, and the following patents which have a certain relationship with the present invention are mainly included:
1. the application number is CN202110222321.6, the name is 'a finite element parametric modeling method for a bolt and a nut divided by hexahedron meshes', the application is Chinese invention patent of Tianjin university, the patent discloses a finite element parametric modeling method for a bolt and a nut divided by hexahedron meshes, which comprises the following steps: s1: determining the geometrical characteristics of the bolt and the nut; s2: determining the grid characteristics of the bolt and the nut; s3: dividing the bolt into a nut area, a screw area, a transition area and a thread area along the axial direction, and constructing a bolt and nut single-section node coordinate perpendicular to the axial direction of the bolt by using a segmented expression method; s4: copying, translating, numbering and rotating the obtained single-section nodes of the bolts and the nuts along the axis direction of the bolts so as to construct a whole bolt and nut model node coordinate matrix; s5: the node coordinates of the bolts and the nuts obtained in the step S4 are regulated according to the connection sequence of the nodes of the hexahedron unit with eight nodes in the finite element software; s6: the node coordinates and cell number matrix obtained in step S3 and step S5 are derived.
2. The application number is CN201210267815.7, which is named as a modeling method in three-dimensional simulation of a complex motion system, and the applicant is a Chinese invention patent of the 92232 army of the people's liberation army of China, and the patent discloses a modeling method in three-dimensional simulation of a complex motion system, and the method comprises the following steps: (1) the system composition is split in a layered mode, and each level of subsystem is split step by step to a component unit in the complex motion system; (2) respectively determining position coordinates, rotation angles, relative position relations and relative angle rotation relations among the sub-systems and the component units after being disassembled; (3) recording the independent and relative motion relationship of each level of subsystem and component unit; (4) and (4) constructing a whole system and a motion model of each level of subsystem according to the data relation of the steps (1), (2) and (3).
3. The application number is CN201110302715.9, which is named as 'a bolt finite element parametric modeling method capable of realizing hexahedral mesh division', the application is a Chinese patent invention of the Western university of transportation, and the patent discloses a bolt finite element parametric modeling method capable of realizing hexahedral mesh division, wherein the key geometric dimension of a bolt model is parameterized in finite element software, a thread single-section model vertical to the axial direction of a bolt is constructed by a segmented expression method, each thread section model on a single pitch is constructed by a key point translation and rotation method, thread segment models with the number of threads are generated along the axial direction of the bolt by a body generation, copy and translation method, screw rod models without threads are divided according to the division limit of the thread body models along the axial line, a screw rod and thread joint part are subjected to individual body generation processing, the bolt head body models are divided according to the division limit of the thread body models along the axial line, and dividing the whole bolt body model by using the Mapped map to finally generate a bolt hexahedron grid unit model.
Through careful analysis of the above patents, although some simulation analysis methods and some improvement solutions have been proposed, through careful analysis, these patents still adopt the conventional method of integral modeling analysis, and though different methods are adopted in grid division, the problem that the modeling workload is large according to the conventional simulation analysis for complex parts, especially for multilayer thin-wall rubber-metal composite revolving bodies which are subjected to various loads and are provided by the present invention, still exists, and therefore, the problems still need to be further researched and improved.
Disclosure of Invention
The invention aims to overcome the defects that the simulation analysis efficiency of the conventional simulation analysis method of a multilayer thin-wall rubber-metal composite revolving body structure under the multi-load working conditions of radial, axial, deflection, torsion and the like is low, the difficulty of a three-dimensional mesh generation technology is high, and the non-convergence of calculation is easily caused, and provides a method for performing simulation analysis on meshes and stress strain of a two-dimensional axisymmetric section in a 180-degree or 360-degree rotation mode.
In order to achieve the purpose, the invention provides a simulation analysis method of a multilayer thin-wall rubber-metal composite revolving body structure, which comprises the following steps of firstly carrying out two-dimensional mesh subdivision operation on the multilayer thin-wall rubber-metal composite revolving body; after the mesh subdivision of the two-dimensional axisymmetric cross section is completed, carrying out axial or torsional load simulation analysis calculation; after the simulation analysis of the axial load or the torsional load is finished, converting the two-dimensional axisymmetric grid into a three-dimensional grid, and simultaneously synchronously mapping the stress-strain result of the two-dimensional axisymmetric model on the newly generated three-dimensional grid as the initial stress-strain of the simulation analysis of the deflection or radial load working condition; and then, the three-dimensional model simulation analysis under the working condition of deflection or radial load after axial or torsional preloading is carried out continuously, and finally, the simulation analysis calculation under the working condition of multiple loads of axial, radial, torsional and deflection is completed.
Furthermore, the two-dimensional mesh generation operation is to extract a two-dimensional axisymmetric section from the multilayer thin-wall rubber-metal composite revolving body structure and then perform the two-dimensional mesh generation operation by using the section; the two-dimensional mesh generation operation is to extract a two-dimensional axisymmetric section and then generate mesh generation required by simulation analysis based on large deformation of rubber on the section.
Furthermore, the mesh division adopts a finite element mesh division method, when the finite element method is adopted for structural analysis, firstly, the structure must be dispersed to form a finite element mesh, and various information corresponding to the mesh is given.
Further, the various information includes unit information, node coordinates, material information, constraint information, and load information; the load information includes deflection load, radial, axial load, or torsional load information.
Furthermore, the discretization of the structure is to discretize a boundary curve, namely, to arrange points on the boundary according to the requirement of a density control function; the finite element mesh forming is to generate units in a target area, and the units and nodes are generated according to specific conditions and divided into boundary unit division and internal unit division; and optimizing the finally generated grid unit by using various information.
Further, the axial or torsional load simulation analysis calculation is to continue to set up material properties, boundary conditions and other conventional simulation analysis preprocessing work after the mesh subdivision of the two-dimensional axisymmetric cross section is completed, apply an axial or torsional load suitable for the two-dimensional axisymmetric cross section model, and perform simulation analysis based on the axial or torsional load on the two-dimensional axisymmetric cross section model.
Further, when the two-dimensional axisymmetric cross-section model is subjected to simulation analysis based on axial or torsional loads, whether radial or deflection loads have axial preload or not is judged, and different three-dimensional model simulation analysis is performed according to whether the axial preload exists or not.
Further, the three-dimensional model simulation analysis under the working condition of deflection or radial load after axial or torsional preloading is that axial preload exists after the axial or torsional load simulation analysis is completed; and rotating the two-dimensional axisymmetric grid by 180 degrees or 360 degrees to generate a three-dimensional grid, and simultaneously synchronously rotating the stress-strain result of the two-dimensional axisymmetric model by 180 degrees or 360 degrees to map the stress-strain result on the newly generated 180-degree or 360-degree three-dimensional grid as initial stress strain of deflection or radial load working condition simulation analysis. And (4) continuing to perform simulation analysis on the 180-degree or 360-degree three-dimensional model under the deflection or radial load working condition after axial or torsional preloading, and finally completing simulation analysis calculation under the axial, radial, torsional and deflection multi-load working condition.
Further, the three-dimensional model simulation analysis under the working condition of deflection or radial load after axial or torsional preloading is that no axial preload exists after the axial or torsional load simulation analysis is completed; and (3) rotating the two-dimensional axisymmetric section model by 180 degrees or 360 degrees to generate a three-dimensional grid, and then performing simulation analysis based on radial or deflection load working conditions.
Further, the simulation analysis method of the multilayer thin-wall rubber-metal composite revolving body structure comprises the following steps:
step 1: carrying out two-dimensional axisymmetric analysis on the to-be-split grid of the multilayer thin-wall rubber-metal composite revolving body, and selecting a two-dimensional axisymmetric section;
step 2: extracting a two-dimensional axisymmetric section from a multilayer thin-wall rubber-metal composite revolving body structure for performing two-dimensional mesh subdivision operation on the section; and carrying out mesh subdivision required by simulation analysis based on large deformation of rubber on the section; generating a two-dimensional grid beneficial to large deformation of rubber on a two-dimensional axisymmetric section;
and step 3: after the mesh subdivision of the two-dimensional axisymmetric cross section is completed, conventional simulation analysis pretreatment work such as material properties, boundary conditions and the like is continuously set;
and 4, step 4: carrying out simulation analysis based on axial or torsional load on the two-dimensional axisymmetric section model; applying an axial or torsional load suitable for the two-dimensional axisymmetric section model, and then carrying out simulation analysis calculation under the axial or torsional load;
and 5: after the simulation analysis under the axial load or the torsional load is finished, rotating the two-dimensional axisymmetric grid by 180 degrees or 360 degrees to generate a three-dimensional grid, and simultaneously synchronously rotating the stress-strain result of the two-dimensional axisymmetric model by 180 degrees or 360 degrees to map the newly generated 180-degree or 360-degree three-dimensional grid as the initial stress-strain of the deflection load working condition simulation analysis;
step 6: continuing to perform simulation analysis on the newly generated 180-degree or 360-degree three-dimensional model under the deflection load working condition of the 180-degree three-dimensional grid;
and 7: finally, axial or torsional and deflection multi-load working condition simulation analysis and calculation are completed.
The invention has the advantages that:
the simulation analysis method for the multilayer thin-wall rubber-metal composite revolving body structure has the advantages that the simulation analysis of the axial and torsional load working conditions can be realized by using a simple two-dimensional axisymmetric model, meanwhile, the simulation analysis of the deflection and radial load working conditions is realized by adopting a 180-degree or 360-degree rotation mode, and the analysis efficiency of the multilayer thin-wall rubber-metal composite revolving body structure is greatly improved. The method overcomes the defects that the simulation analysis efficiency of the existing simulation analysis method of the multilayer thin-wall rubber-metal composite revolving body structure is low under the multi-load working conditions of radial, axial, deflection, torsion and the like, the difficulty of the three-dimensional mesh generation technology is high, and the calculation is not converged easily.
Drawings
FIG. 1 is a flow chart of one embodiment of a simulation analysis method suitable for a multilayer thin-walled rubber-metal composite rotary body structure of the present invention;
FIG. 2 is a schematic representation of a three-dimensional geometric model (full 360 degree model) of a multilayer thin-walled rubber-metal composite solid of revolution structure according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of a three-dimensional geometric model (270 degree model) of a multi-layer thin-walled rubber-metal composite solid of revolution structure according to an embodiment of the present invention;
FIG. 4 is a two-dimensional axisymmetric cross-section of a multi-layer thin-walled rubber-metal composite solid of revolution structure according to an embodiment of the present invention;
FIG. 5 is a two-dimensional axisymmetric cross-sectional mesh generation model of an embodiment of the present invention;
FIG. 6 is a three-dimensional mesh model generated by 180 degree rotation of a mesh of an axisymmetric cross-section according to an embodiment of the present invention;
FIG. 7 is a stress cloud of a two-dimensional axisymmetric cross-section after axial preload of an embodiment of the present invention;
FIG. 8 is a three-dimensional initial stress cloud after 180 degrees rotation of an axisymmetric cross-section of an embodiment of the present invention;
FIG. 9 is a strain cloud of a two-dimensional axisymmetric cross-section after axial preload of an embodiment of the present invention;
FIG. 10 is a three-dimensional initial strain cloud after 180 degrees rotation of an axisymmetric cross-section of an embodiment of the present invention;
FIG. 11 is a three-dimensional stress cloud representation of an axisymmetric cross-section of an embodiment of the present invention after 180 degrees of rotation and a simulation analysis based on deflection load conditions;
FIG. 12 is a three-dimensional strain cloud representation of an axisymmetric cross-section of an embodiment of the present invention after 180 degrees rotation and simulation analysis based on deflection loading conditions.
Reference numerals: 1. a top plate; 2. an elastomer; 3. a metal separator; 4. a base plate; 5. rubber; 6. two-dimensional axisymmetric cross section, 7, multilayer thin-wall rubber metal composite rubber pile.
Detailed Description
The invention is further illustrated with reference to the following figures and specific examples.
Example one
As can be seen from the attached drawings 1 and 2, the invention relates to a simulation analysis method of a multilayer thin-wall rubber metal composite rubber pile, wherein the multilayer thin-wall rubber metal composite rubber/7 comprises a top plate 1, a bottom plate 4 and an elastic body 2; the elastomer 2, the top plate 1 and the bottom plate 4 are integrally vulcanized together, and the elastomer 2 is a composite part (shown in figure 2) formed by compounding a plurality of layers of rubber 5 and a metal separator 3 and is symmetrically distributed around a central axis; firstly, a two-dimensional axisymmetric section 6 (shown in figures 3 and 4) is extracted from a multilayer thin-wall rubber-metal composite rubber pile by cutting through the center, and two-dimensional mesh generation operation is carried out on the two-dimensional axisymmetric section 6 (shown in figure 4). After extracting a two-dimensional axisymmetric section, the mesh subdivision required by simulation analysis based on large deformation of rubber is carried out on the section (as shown in figure 5). After the mesh subdivision of the two-dimensional axisymmetric cross section is completed, conventional simulation analysis pretreatment work such as material properties, boundary conditions and the like is continuously set, and axial or torsional load suitable for the two-dimensional axisymmetric cross section model is applied to perform simulation analysis calculation. After the axial or torsional load simulation analysis is completed, the two-dimensional axisymmetric grid is rotated by 180 degrees to generate a three-dimensional grid (as shown in fig. 6), and simultaneously, the stress-strain result of the two-dimensional axisymmetric model is synchronously rotated by 180 degrees and mapped on the newly generated 180-degree three-dimensional grid to serve as initial stress strain of deflection or radial load working condition simulation analysis. And (4) continuing to perform simulation analysis on the 180-degree three-dimensional model under the deflection or radial load working condition after axial or torsional preloading, and finally completing simulation analysis and calculation under the axial, radial, torsional and deflection multi-load working condition.
The specific simulation analysis steps are as follows:
step 1: extracting a two-dimensional axisymmetric cross section from the multilayer thin-wall rubber-metal composite rubber pile structure for performing two-dimensional mesh generation operation on the cross section;
step 2: mesh subdivision required by simulation analysis based on large deformation of rubber is carried out on the section;
and step 3: after the mesh subdivision of the two-dimensional axisymmetric cross section is completed, conventional simulation analysis pretreatment work such as material properties, boundary conditions and the like is continuously set; because the multilayer thin-wall rubber-metal composite rubber pile has no torsional load, only an axial load suitable for a two-dimensional axisymmetric section model is applied, and then simulation analysis calculation under the axial load is carried out;
and 4, step 4: after the axial load simulation analysis is completed, rotating the two-dimensional axisymmetric grid by 180 degrees to generate a three-dimensional grid, and simultaneously synchronously rotating the stress-strain result of the two-dimensional axisymmetric model by 180 degrees to map the stress-strain result on the newly generated 180-degree three-dimensional grid as the initial stress strain of the deflection load working condition simulation analysis;
and 5: continuing to perform simulation analysis on the 180-degree three-dimensional model under the deflection load working condition after the axial preloading;
step 6: finally, axial and deflection multi-load working condition simulation analysis calculation is completed.
Wherein:
in the step 1, the structural characteristics of the multilayer thin-wall rubber metal composite rubber pile are firstly analyzed, the multilayer thin-wall rubber metal composite rubber pile is subjected to plane segmentation in the radius direction along the central axis, two-dimensional sections are extracted from the rubber and the spacer, a dwg file is generated, and Cad software is adopted for processing. And simplifying the two-dimensional section to different degrees according to the calculation requirement. When the rigidity performance is calculated, the molded surface encapsulation with the thickness less than 1mm can be removed; when calculating the fatigue performance, the original design structure is maintained as much as possible, and the processed two-dimensional cross section is shown in fig. 3. Stress cloud analysis of the two-dimensional axisymmetric cross-section after post-axial preload, the analysis results are shown in FIG. 6.
In the step 2, a dxf format file is derived from the processed two-dimensional cross section of the axisymmetric part, the two-dimensional cross section is introduced into finite element simulation software (such as ABAQUS) to be subjected to meshing, the two-dimensional cross section is divided into parts for meshing when the meshing is performed, the whole cross section is divided into a plurality of different meshes according to stress conditions, and the result is shown in fig. 4. The exported dxf format can also be imported into hypermesh software for meshing.
In the step 3, after material properties and boundary conditions are set according to the stress condition of the multilayer thin-wall rubber-metal composite rubber pile, the two-dimensional axisymmetric grid model is subjected to simulation analysis, so as to obtain a stress-strain result of the two-dimensional axisymmetric grid under the condition of axial preload, as shown in fig. 7 and 9.
In the step 4, the two-dimensional axisymmetric grid can be rotated by 180 degrees by a revalve command of the ABAQUS software to generate a three-dimensional grid as an initial grid model for deflection load condition simulation analysis, and simultaneously, a deformed stress-strain result page of the two-dimensional axisymmetric model under the load condition in the step 3 can be synchronously rotated and mapped onto the newly generated 180-degree three-dimensional grid as an initial stress-strain for deflection load condition simulation analysis by a symmetry stresses transfer command, as shown in fig. 8 and 10.
And the simulation analysis calculation under the axial load and the simulation analysis of the 180-degree three-dimensional model are carried out according to a conventional simulation analysis calculation method. The simulation analysis comprises stress cloud picture analysis of a two-dimensional axisymmetric section after axial preloading, three-dimensional initial stress cloud picture analysis of the two-dimensional axisymmetric section after rotating 180 degrees, strain cloud picture analysis of the two-dimensional axisymmetric section after axial preloading, three-dimensional stress cloud picture analysis of the three-dimensional initial strain cloud picture after the axisymmetric section rotates 180 degrees and then simulation analysis based on the deflection load working condition and three-dimensional strain cloud picture analysis of the three-dimensional initial strain cloud picture after the axisymmetric section rotates 180 degrees and then simulation analysis based on the deflection load working condition.
The stress cloud graph of the two-dimensional axisymmetric section after the axial preload is shown in the attached figure 7;
the three-dimensional initial stress cloud picture of the two-dimensional axisymmetric section after rotating 180 degrees is shown in the attached figure 8;
the strain cloud graph of the two-dimensional axisymmetric section after the axial preload is shown in the attached figure 9;
the three-dimensional initial strain cloud picture of the axisymmetric cross section after rotating 180 degrees is shown in the attached figure 10;
the three-dimensional stress cloud picture after the axial symmetry section is rotated by 180 degrees and subjected to simulation analysis based on deflection load working conditions is shown in the attached figure 11;
the three-dimensional strain cloud picture after the axial symmetry section is rotated by 180 degrees and subjected to simulation analysis based on deflection load working conditions is shown in the attached figure 12;
example two
The principle of the second embodiment is the same as that of the first embodiment, only the part structure is slightly different, and the second embodiment relates to a simulation analysis method of a multilayer thin-wall rubber metal spherical hinge, wherein firstly, a two-dimensional axisymmetric section is extracted by cutting the multilayer thin-wall rubber metal spherical hinge along the central axis, and the two-dimensional axisymmetric section is used for carrying out two-dimensional mesh subdivision operation; after a two-dimensional axisymmetric section is extracted, mesh subdivision required by simulation analysis based on large deformation of rubber is carried out on the section. After the mesh subdivision of the two-dimensional axisymmetric cross section is completed, conventional simulation analysis pretreatment work such as material properties, boundary conditions and the like is continuously set, and axial or torsional load suitable for the two-dimensional axisymmetric cross section model is applied to perform simulation analysis calculation. After the axial or torsional load simulation analysis is completed, the two-dimensional axisymmetric grid is rotated by 360 degrees to generate a three-dimensional grid, and simultaneously, the stress-strain result of the two-dimensional axisymmetric model is synchronously rotated by 360 degrees and mapped on the newly generated 360-degree three-dimensional grid to serve as the initial stress strain of the deflection or radial load working condition simulation analysis. And (4) continuing to perform simulation analysis on the 360-degree three-dimensional model under the deflection or radial load working condition after axial or torsional preloading, and finally completing simulation analysis and calculation under the axial, radial, torsional and deflection multi-load working condition.
The simulation analysis of the multilayer thin-wall rubber metal spherical hinge specifically comprises the following steps:
step 1: extracting a two-dimensional axisymmetric section from the multilayer thin-wall rubber-metal composite revolving body structure for performing two-dimensional mesh generation operation on the section;
step 2: mesh subdivision required by simulation analysis based on large deformation of rubber is carried out on the section;
and step 3: after the mesh subdivision of the two-dimensional axisymmetric cross section is completed, the conventional simulation analysis pretreatment work such as material properties, boundary conditions and the like is continuously set, an axial load suitable for the two-dimensional axisymmetric cross section model is applied, and then simulation analysis calculation under the axial load is carried out;
and 4, step 4: after the simulation analysis under the axial or torsional load is finished, rotating the two-dimensional axisymmetric grid for 360 degrees to generate a three-dimensional grid, and simultaneously synchronously rotating the stress-strain result of the two-dimensional axisymmetric model for 360 degrees and mapping the stress-strain result on the newly generated 360-degree three-dimensional grid as the initial stress strain of the deflection load working condition simulation analysis;
and 5: continuing to perform simulation analysis on the newly generated 360-degree three-dimensional model under the deflection load working condition of the 360-degree three-dimensional grid;
step 6: finally, axial or torsional and deflection multi-load working condition simulation analysis and calculation are completed.
The above listed embodiments are only for clear and complete description of the technical solution of the present invention with reference to the accompanying drawings; it should be understood that the embodiments described are only a part of the embodiments of the present invention, and not all embodiments, and the terms such as "upper", "lower", "front", "back", "middle", etc. used in this specification are for clarity of description only, and are not intended to limit the scope of the invention, which can be implemented, and the changes or modifications of the relative relationship thereof are also regarded as the scope of the invention without substantial technical changes. Meanwhile, the structures, the proportions, the sizes, and the like shown in the drawings are only used for matching with the contents disclosed in the specification, so as to be understood and read by those skilled in the art, and are not used for limiting the conditions under which the present invention can be implemented, so that the present invention has no technical essence, and any structural modification, changes in proportion relation, or adjustments of the sizes, can still fall within the range covered by the technical contents disclosed in the present invention without affecting the effects and the achievable purposes of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention has the advantages that:
the simulation analysis method for the multilayer thin-wall rubber-metal composite revolving body structure has the advantages that the simulation analysis of the axial and torsional load working conditions can be realized by using a simple two-dimensional axisymmetric model, meanwhile, the simulation analysis of the deflection and radial load working conditions is realized by adopting a 180-degree or 360-degree rotation mode, the simulation analysis efficiency of the multilayer thin-wall rubber-metal composite revolving body structure under the multi-load working conditions of radial, axial, deflection and torsion and the like can be greatly improved, and the defect that the efficiency is low because the axial and torsion are calculated by adopting the axisymmetric model in the conventional simulation analysis method, and the deflection and radial are calculated by adopting the 180-degree or 360-degree model or all the multi-load working conditions of axial, radial, deflection and torsion are calculated by adopting the 360-degree model is overcome. Has the following advantages:
1. the method overcomes the defects that the simulation analysis efficiency of the existing simulation analysis method of the multilayer thin-wall rubber-metal composite revolving body structure is low under the multi-load working conditions of radial, axial, deflection, torsion and the like, the difficulty of the three-dimensional mesh generation technology is high, and the calculation is not converged easily;
2. the method can effectively overcome the defect that the axial direction and the torsion in the conventional simulation analysis method are calculated by adopting an axisymmetric model;
3. according to the method, the 180-degree model or 360-degree model is adopted for calculation in the deflection and radial directions, so that the workload of grid division can be greatly reduced, and the grid division time can be saved by more than 50% for complex symmetrical products.
The noun explains:
1) multilayer thin wall rubber metal composite structure: a rubber metal composite product comprising a plurality of layers of rubber and a plurality of layers of separators;
2) a revolving body structure: the axially symmetric model in engineering application refers to a product three-dimensional model which can be formed by rotating a certain section by 360 degrees around a central shaft and covers products such as rubber spherical hinges, rubber bushings, rubber nodes, conical springs and the like.
Claims (10)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111207879.3A CN113987868A (en) | 2021-10-18 | 2021-10-18 | Simulation analysis method for multilayer thin-wall rubber-metal composite revolving body structure |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111207879.3A CN113987868A (en) | 2021-10-18 | 2021-10-18 | Simulation analysis method for multilayer thin-wall rubber-metal composite revolving body structure |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113987868A true CN113987868A (en) | 2022-01-28 |
Family
ID=79739028
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111207879.3A Pending CN113987868A (en) | 2021-10-18 | 2021-10-18 | Simulation analysis method for multilayer thin-wall rubber-metal composite revolving body structure |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113987868A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116663190A (en) * | 2023-06-06 | 2023-08-29 | 嘉丰盛精密电子科技(孝感)有限公司 | Method for identifying splicing strength of stamping parts in shielding cover |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106599503A (en) * | 2016-12-23 | 2017-04-26 | 大工（青岛）新能源材料技术研究院有限公司 | Finite element analysis-based winter-heating parallel jointed steel-structured greenhouse simulation method |
-
2021
- 2021-10-18 CN CN202111207879.3A patent/CN113987868A/en active Pending
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106599503A (en) * | 2016-12-23 | 2017-04-26 | 大工（青岛）新能源材料技术研究院有限公司 | Finite element analysis-based winter-heating parallel jointed steel-structured greenhouse simulation method |
Non-Patent Citations (2)
| Title |
|---|
| 孙存政: "波纹换热管的研究与开发", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》, vol. 2010, no. 7, 15 July 2010 (2010-07-15), pages 3 * |
| 李林: "鲹科鱼类游动力学特性的理论与仿真研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》, vol. 2011, no. 2, 15 December 2011 (2011-12-15), pages 4 - 5 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116663190A (en) * | 2023-06-06 | 2023-08-29 | 嘉丰盛精密电子科技(孝感)有限公司 | Method for identifying splicing strength of stamping parts in shielding cover |
| CN116663190B (en) * | 2023-06-06 | 2023-11-07 | 嘉丰盛精密电子科技(孝感)有限公司 | Method for identifying splicing strength of stamping parts in shielding cover |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Shi et al. | Vibration and buckling analysis of curvilinearly stiffened plates using finite element method | |
| CN106021644B (en) | A Method for Determining Coefficients of Interfacial Constraint Equations in Mixed-Dimensional Models | |
| CN108984829B (en) | Calculation method and system for accumulation process of rockfill concrete rockfill | |
| CN112836316B (en) | A Finite Element Parametric Modeling Method for Bolts and Nuts Based on Hexahedral Meshing | |
| Yuan et al. | An equivalent modeling method for honeycomb sandwich structure based on orthogonal anisotropic solid element | |
| Zhang et al. | A construction method of digital twin model for contact characteristics of assembly interface | |
| Dewangan et al. | Numerical deflection and stress prediction of cutout borne damaged composite flat/curved panel structure | |
| CN113505435A (en) | Carbon fiber shell analysis method based on response surface model | |
| CN117725859A (en) | Fatigue crack growth prediction method and system for aircraft fatigue damage key part | |
| Wang et al. | Voronoi polygonal hybrid finite elements with boundary integrals for plane isotropic elastic problems | |
| CN112528415B (en) | Method for analyzing macro-micro failure mode of composite material shaft structure | |
| CN113987868A (en) | Simulation analysis method for multilayer thin-wall rubber-metal composite revolving body structure | |
| CN107368660B (en) | An Analytical Method for Predicting the Yield Surface of Periodic Lattice Materials | |
| Ho-Huu et al. | An efficient procedure for lightweight optimal design of composite laminated beams | |
| Larwood et al. | Domain decomposition approach for parallel unstructured mesh generation | |
| Park et al. | Multidisciplinary design optimization of a structurally nonlinear aircraft wing via parametric modeling | |
| Ansari et al. | Compression buckling of elastically supported cylindrical shells with Bend/Twist coupling | |
| Al Mahmoud et al. | Implementation of different types of meshfree technique in computational solid mechanics: a comprehensive review across nano, micro, and macro scales | |
| Venkataraman et al. | Structural optimization: what has Moore's law done for us? | |
| CN113191016A (en) | Body expression model-based multi-material product modeling and analyzing integrated method | |
| CN112507445A (en) | Geometrical structure reconstruction method of topological optimization result | |
| CN116305589A (en) | Method, system, equipment and medium for reduced order analysis of helicopter blade structure | |
| CN115910239A (en) | Laminated plate modeling method and device based on composite material layering design | |
| Chang | Simulation Analysis of Finite Element Preprocessing based on HyperMesh | |
| CN110781621B (en) | A Fast Reanalysis Method for Bearing Response of Thin-Walled Structures with Geometric Defects |
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 |